数据可视化部分答案

1、数据可视化（即画图）部分作业x = a cost e1.已知椭圆的长、短轴a=4,b=2，用“小红点线”画椭圆丿。(参见图P5-1)y = bsi nt程序：t=2*pi*(0:199)/199;a=4;b=2;x=a*cos(t);y=b*sin(t);plot(x,y,r.,MarkerSize,15)axis equal-2-101234图 p5-12.根据表达式=1-C0S绘制如图p5-2的心脏线。（提示：采用极坐标绘线指令polar） clft = 0:.01:2*pi;P=1-cos(t); pline=polar(t,P,r), set(pline,LineWidth,5)tit

2、le(P=1-costheta)图 p5-23. A,B,C三个城市上半年每个月的国民生产总值见表p5.1。试画出如图p5-3所示的三城市上半年每月生产总值的累计直方图。（提示：使用指令bar）表p5.1各城市生产总值数据（单位：亿元）城市1月2月3月4月5月6月A170120180200190220B120100110180170180C70508010095120clearx=1:6;Y=170 120180200190220;120100110180170180;70508010095120bar(x,Y,grouped);bar(x，丫 ,stacked);colormap(cool)

3、;% lege nd(A,B,C,Locatio n，Northwest)lege nd(A,B,C,2)图 p5-34.二阶线性系统的归一化（即令-1 ）冲激响应可表示为y(t)二齐亠 sin( P t)te丄J_edm2 :0 _ ： 1=1 ，其中1 f|1 -，2 |， 为阻尼1系数。试在同一张图上，画出=0.2 : 0.2: 1.4不同取值时的各条曲线，时间区间为r 0,18。此外，要求（1）：1的各条曲线用蓝色，1的用红色，二-1的用黑色，且线较粗；（2）=0.2和 =1.4的两条曲线要有醒目的阻尼系数标志。（本题要求图形参见图p5-4)% exmp504.mclc,clf,cle

4、ar;t=(0:0.05:18);N=length(t);zeta=0.2:0.2:1.4;L=length(zeta); y=zeros(N,L); hold on for k=1:L zk=zeta(k); beta=sqrt(abs(1-zkA2); if zk-1-2*epsy=1/beta*exp(-zk*t).*sin(beta*t); plot(t,y,b) if zkv0.4text(2.2,0.63,zeta = 0.2) endelseif (zk-1)1.2text(0.3,0.14,zeta = 1.4)endendendtext(10,0.7,Deltazeta=0.

5、2) axis(0,18,-0.4,0.8) hold offbox ongrid on图 p5-45.用绿实线绘制x=sin(t) , y = cos(t) , z=t的三维曲线，曲线如图p5-5所示。(提示:使用plot3指令) elf t=4*pi*(O:1OO)/1OO; x=sin(t);y=cos(t);z=t;.x2 _y26.在区域x,y-3,3，绘制z=4xe 的如图p5-6的三维(透视)网格曲面。 (不得使用 ezmesh)elfx=-3:0.1:3;y=x;X,Y=meshgrid(x,y);Z=4.*X.*exp(-X.A2-Y.A2);mesh(X,Y,Z)hidde

6、n off% colormap(cool),% shading interp,syms x y z% z=4.*x.*exp(-x.A2-y.A2);z=4*x*exp(-xA2-yA2);ezmesh(z,-3,3)hidden offx expf-xy2) 4.0图 p5-67.在x,y-4：,4二区间里，根据表达式z=y)，绘制如图p5-7所示的曲面。x + yclear allx=4*pi*(-50:50)/50;y=x;X,Y=meshgrid(x,y);Z=sin(X+Y)./(X+Y+(X+Y=O)*eps);surf(X,Y,Z)view(21,32)shading inter

7、p图 p5-7&试用图解法回答:(1) 方程组丿(d +y2) =0.1有多少个实数解？sin(x cos(y) = 0(2) 求出离x =0,y =0最近、且满足该方程组的一个近似解。(1)ezplot(y心+xA2+yA2)-0.1,-2*pi,2*pi,-pi/3,3.5*pi) hold onezplot(sin(x+cos(y),-2*pi,2*pi,-pi/3,3.5*pi)图 p5-8可看到6个交点，即方程组有6个实数解(2)edit sushufunction F=myfun8(x,y)F=y/(1+xA2+yA2)-0.1;sin(x+cos(y)end zoom onxy= ginput(1) f=fsolve(myfun8,xy) X,Y = ginput(1) f=fsolv