数据可视化部分答案.docxVIP

数据可视化（即画图）部分作业 工一a cos t 己知椭圆的长、短轴a = 4,b = 2，用“小红点线”画椭圆彳~ 。（参见图p5-l） y - bsint 程序： 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 ?3 ?2 0 12 3 4 图 p5-l 根据表达式p = l-cos绘制如图p5-2的心脏线。（提示：采用极朋标绘线指令polar） elf t = 0:.01:2*pi; P=l-cos(t); pline=polar(t,P/r*), set(pline,LineWidth*,5) title(,P=l-cos\thetaf) Q Figure 1 [“回 图 p5-2 A,B,C三个城市上半年每个刀的国民生产总值见表p5.1o试画出如图p5-3所示的三城市 上半年每月生产总值的累计直方图。（提示：使用指令bar） 表p5」各城市生产总值数据（单位：亿元） 城市 1刀 2 n 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=l:6; Y= [170 120 180 200 190 220; 120 100 110 180 170 180; 70 50 80 100 95 120] bar(x#Y1#1 grouped1); bar(x#Y1, 1 stacked1); colormap(cool); % legend (AJBJC, location ■ z 1 Northwest1) legend( -A1 , -B1 z ,2) 图 p5-3 图 p5-3 e~^1 sin(0 r)丘-(0二阶线性系统的归一化（即令? = 1 e~^1 sin(0 r) 丘-(0 0^l ^ = 1 ,其中0 = Jll —厂|, $为阻尼 $1 系数。试在同一张图上，画出^ = 0.2:0.2:14不同取值时的各条曲线，时间区间为 /g[0,18]o此外，要求（1） $vl的各条曲线用蓝色，$1的用红色，$ = 1的用黑 色，且线较粗；（2） $ = 0.2和$ = 1.4的两条III]线要有解目的阻尼系数标志。（本题 耍求图形参见图p5-4) % exmp504.nl clc,clf,clear; 1=(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=l:L zk=zeta(k); beta=sqrt(abs( 1 -zk A 2)); if zk-l-2*eps y=l/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-l)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；r9 ifzk1.2 text(0.3,0.14/\zeta = 1.4*) end end end text( 10,0.7,*\DeIta\zeta=0.2*) axis([09l 8,-0.4,0.8]) hold off box on grid on 图 p5-4 用绿实线绘制x = sin(f) , y = cos(f), z = f的三维|11|线，|11|线如图p5-5所示。(提示: 使用plot3指令) elf t=4*pi*(0:100)/100; x=sin(t);y=cos(t);z=t; plot3(x,y,z/gf/LineWidth\3),box on -1 -1 图p5?5 在区域x,yw[—3,3],绘制z = 4牝一宀声的如图p5?6的三维(透视)网格曲面。(不得 使川 ezmesh) elf x=-3:0.1:3;y=x;[X,Y]=meshgrid(x,y); Z=4.*X.*exp(-X.A2-Y.A2); mesh(X,Y,Z) hidden off % colormap(cool), % shading interp, syms x y z % z=4.*x.*exp(-x.A2-y.A2); z=4*x*exp(-xA2-yA2); ezmesh(zJ-3,3]) hidden off