数据可视化部分答案.docVIP

数据可视化(即画图)部分作业 已知椭圆的长、短轴，用“小红点线”画椭圆。（参见图p5-1） 程序： 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 图 p5-1 根据表达式绘制如图p5-2的心脏线。（提示：采用极坐标绘线指令polar） clf t = 0:.01:2*pi; P=1-cos(t); pline=polar(t,P,r), set(pline,LineWidth,5) title(P=1-cos\theta) 图 p5-2 A,B,C三个城市上半年每个月的国民生产总值见表p5.1。试画出如图p5-3所示的三城市上半年每月生产总值的累计直方图。（提示：使用指令bar） 表p5.1 各城市生产总值数据（单位：亿元） 城市 1月 2月 3月 4月 5月 6月 A 170 120 180 200 190 220 B 120 100 110 180 170 180 C 70 50 80 100 95 120 clear x=1:6; Y=[170 120 180 200 190 220; 120 100 110 180 170 180; 70 50 80 100 95 120] bar(x,Y,grouped); bar(x,Y,stacked); colormap(cool); % legend(A,B,C,Location,NorthWest) legend(A,B,C,2) 图 p5-3 二阶线性系统的归一化（即令）冲激响应可表示为：，其中，为阻尼系数。试在同一张图上，画出不同取值时的各条曲线，时间区间为。此外，要求（1）的各条曲线用蓝色，的用红色，的用黑色，且线较粗；（2）和的两条曲线要有醒目的阻尼系数标志。（本题要求图形参见图p5-4） % exmp504.m clc,clf,clear; 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-zk^2)); if zk-1-2*eps y=1/beta*exp(-zk*t).*sin(beta*t); plot(t,y,b) if zk0.4 text(2.2,0.63,\zeta = 0.2) end elseif (zk-1)2*eps% y=t.*exp(-t); plot(t,y,k,LineWidth,2) else y=(exp(-(zk-beta)*t)-exp(-(zk+beta)*t))/(2*beta); plot(t,y,r) if zk1.2 text(0.3,0.14,\zeta = 1.4) end end end text(10,0.7,\Delta\zeta=0.2) axis([0,18,-0.4,0.8]) hold off box on grid on 图 p5-4 用绿实线绘制，，的三维曲线，曲线如图p5-5所示。（提示：使用plot3指令） clf t=4*pi*(0:100)/100; x=sin(t);y=cos(t);z=t; plot3(x,y,z,g,LineWidth,3),box on 图 p5-5 在区域，绘制的如图p5-6的三维（透视）网格曲面。（不得使用ezmesh） clf x=-3:0.1:3;y=x;[X,Y]=meshgrid(x,y); Z=4.*X.*exp(-X.^2-Y.^2); mesh(X,Y,Z) hidden off % colormap(cool), % shading interp, syms x y z % z=4.*x.*exp(-x.^2-y.^2); z=4*x*exp(-x^2-y^2); ezmesh(z,[-3,3]) hidden off 图 p5-6 在区间里，根据表达式，绘制如图p5-7所示的曲面。 clear all x=4*pi*(-50:50)/50;y=x;[X,Y]=meshgrid(x,y); Z=sin(X