数据可视化部分答案.docx
数据可视化部分答案本页仅作为文档封面，使用时可以删除This document is for reference only-rar21 year.March数据可视化（即画图）部分作业1 .已知椭的长、短轴a=4. b =2用“小红点线”画椭r= AC°SZq（参y - bsint见图P5-I)程序：t=2*pi*(0:199)/199;a=4;b=2;x=a*cos(t);y=b*sin(t);plotfxA/r.'/MarkerSizeAIS)axis equal3-二-立.二手：:Y11II1II '-3-2-10：H1235-12 .根据表达式Q = I-COS&绘制如图p5- 2的心脏线。(提示：采用极坐标绘线指令 polar) elft = 0:.01:2*pi; P=l-cos(t);pline=polar(t/P/T')/ set(pline/LinWidth； 5) title( P=l-costheta )M Figure 1File Edit View Insert Tools Desktop Window Help:J tJ a | Q戈公塔要£ 11 0因 BP=1-COS0图p5-23 - ABC三个城市上半年每个月的国民生产总值见表。试画出如图p53所示的三城市上半年每月生产总值的累计直方图。（提示：使用指令bar）表各城市生产总值数据（单位：亿元）Y却moHLLISO 260190280;城市1月2月3月4JJ5月6月A170120180200190220B120100110180170180C70508010095120clearx=l:6;barfxA，/grouped,);barfxX'stacked1); colormap(cool);% legendCA'/B'/C/Location'/NorthWest'） legend（'ATBTC； 2）图p5-34.二阶线性系统的归一化（即令®二1 ）冲激响应可表示为:0<A<1 y(0 =±严歹二1 ,其中0 = 歹为阻V>1尼系数。试在同一张图上，画出v=0- 2:0- 2:1- 4不同取值时的各条曲线，时间区间为（0,18 o此外，要求（1 ）歹vl的各条曲线用蓝色，：>1的用红色，：=1的用黑色，且线较粗；（2）（本题要求图形参<=0.2和歹=1.4的两条曲线要有醒的阻尼系数标志。见图P5-4 )%cl八df,dear; t=(0:18)' N=length(t); zeta=:; L=length(zeta); y=zeros(N/L);hold onfor k=l:Lzk=zeta(k); beta=sqrt(abs(l-zkA2);if zk-l<-2*eps y=l/beta*exp(-zk*t).*sin(beta*t); plot(t,y/b') if zk<text./Xzeta =')endelseif (zk-l)<2*eps% y=t.*exp(-t);plotltjy/k7LineWidth1A)y=(expelse(-(zk-b(zk+beta)*t)/(2*beta);plotlty/r')if zk>text./Xzeta -)endendendtext(10/,Deltazeta=l) axis(0,18 ) hold offbox ongrid onn au.O.ivA-O.20 X|A.j 1Ln?fl$-1.、0 2OA图 p5- 45用绿实线绘制x = sin(r), y = cos(r), z = t的三维曲线,曲线如图p5-5所示。(提示：使用plot3指令)elft=4*pi*(0:100)/100;x=sin(t);y=cos(t);z=t;plotSlxAz/g'/LineWidthASAbox on图 p5-56在区域x, y e -3,3 绘制z = 4xe 的如图p5-6的三维(透视)网格曲面。(不得使用ezmesh)elfx=-3:3;y=x;X/Y=meshgrid(x/y);Z=4.*X.*exp(-X.A2-Y.A2); mesh(XMZ) hidden off% colormapfcool),% shading interp, syms xy z% z=4.*x.*exp(-x.A2-y.A2); z=4*x*exp(-xA2-yA2); ezmesh(zz-33)hidden off7.在x, ye -4不4刃区间里，根据表达式“曲十刃,绘制如图p5-7所示x + y的曲面。clear 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(2132)shading interp图 p5- 78.试用图解法回答: _01方程组4(1+/ + F)有多少个实数解sin(x + cos(y) = 0(2)求出离x = 0,y = 0最近、且满足该方程组的一个近似解。ezplot(y/(l+xA2+yA2),/-2*pi/2*pi/-pi/3/*pi)hold onezplot('sin(x+cos(y),/-2*pi/2*pi/-pi/3j*pi)图p58可看到6个交点，即方程组有6个实数解edit sushufunction F=myfun8(x/y)F=y/(l+xA2+yA2);sin(x+cos(y) endzoom onxy= ginput(l) f=fsolve() myfun8/xy)X,Y = ginput(l)f=f sol ve( myf u n 8, X, Y)clear allsyms x ys=solg(V/(： L+m2+yA2)Tsin(x+cos(Y)r)离x = 0, y = 0最近的一个近似解为xOO =yOO =
