资源描述:
《MATLAB教程2012a第5章习题解答-张志涌.doc》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、第5章数据和函数的可视化习题5及解答1已知椭圆的长、短轴,用“小红点线”画椭圆。(参见图p5-1)l〖解答〗clfa=4;b=2;t=0:pi/80:2*pi;x=a*cos(t);y=b*sin(t);plot(x,y,'r.','MarkerSize',15)axisequalxlabel('x')ylabel('y')shg图p5-12根据表达式绘制如图p5-2的心脏线。(提示:采用极坐标绘线指令polar)〖解答〗clftheta=0:pi/50:2*pi;rho=1-cos(theta);h=polar(theta,rho,'-r');%极坐标绘线
2、指令。h是所画线的图柄。set(h,'LineWidth',4)%利用set设置h图形对象的“线宽”axissquare%保证坐标的圆整性title('rho=1-costheta')%采用特殊字符映射图p5-21A,B,C三个城市上半年每个月的国民生产总值如见表p5.1。试画出如图p5-3所示的三城市上半年每月生产总值的累计直方图。表p5.1各城市生产总值数据(单位:亿元)城市1月2月3月4月5月6月A170120180200190220B120100110180170180C70508010095120〖目的〗l借助MATLAB的帮助系统,学习直方图
3、指令polar的使用。lbar指令常用格式之一:bar(x,Y,'style')。x是自变量列向量;Y是与x行数相同的矩阵,Y的每一行被作为“一组”数据;style取stacked时,同一组数据中每个元素对应的直方条被相互层叠。l在本例中,Y中的一列代表一个城市。绘图时,各列的颜色,取自“色图”;图形窗的默认色图是jet。〖解答〗x=(1:6)';Y=[170,120,180,200,190,220;120,100,110,180,170,180;70,50,80,100,95,120]';%以下两条指令是为了让读者看清x,Y中对应数据关系而写的。它们不影响
4、画图。disp([blanks(5),'x',blanks(4),'YA',blanks(4),'YB',blanks(4),'YC'])disp([x,Y])%bar(x,Y,'stacked');%层叠直方条colormap(cool);%采用cool色图legend('A','B','C',2);%该指令的第4个输入量,用来控制图例的位置。axis([0,7,0,600])shgxYAYBYC11701207021201005031801108042001801005190170956220180120图p5-31二阶线性系统的归一化(即令)冲激响应可
5、表示为:,其中,为阻尼系数。(1)希望在同一张图上,绘制区间内不同取值时的各条曲线(参见图p5-4)。在此图上,的各条曲线为细蓝线;为粗黑线;为细红线;并且对最上方及最下方的两条曲线给出和的醒目标志。(2)读者运行题下程序exmp504.m,可以发现该程序画出的曲线中没有“粗黑线”。你能讲出原因吗?如何对exmp504.m作最少的修改(比如只改一条指令),就可画出所需图形。(提示:该题深层次地暴露数值计算可能存在的隐患。)〖解答〗clc,clf,clear;t=(0:0.05:18)';N=length(t);zeta=0.2:0.2:1.4;%可能画不出黑
6、线。<3>%zeta=linspace(0.2,1.4,7);L=length(zeta);y=zeros(N,L);holdonfork=1:Lzk=zeta(k);beta=sqrt(abs(1-zk^2));ifzk<1y=1/beta*exp(-zk*t).*sin(beta*t);plot(t,y,'b')ifzk<0.4text(2.2,0.63,'zeta=0.2')endelseifzk==1y=t.*exp(-t);plot(t,y,'k','LineWidth',2)elsey=(exp(-(zk-beta)*t)-exp(-(zk+b
7、eta)*t))/(2*beta);plot(t,y,'r')ifzk>1.2text(0.3,0.14,'zeta=1.4')endendendtext(10,0.7,'Deltazeta=0.2')axis([0,18,-0.4,0.8])holdoffboxongridon图p5-41用绿实线绘制,,的三维曲线,曲线如图p5-5所示。(提示:使用plot3指令)〖解答〗t=(0:0.01:4)*pi;x=sin(t);y=cos(t);z=t;plot3(x,y,z,'-b.','linewidth',3)boxon图p5-5〖解答〗l不管是平面
8、曲线,还是三维立体曲线,它们的描述函数中只有一个自由