大數(shù)據(jù)可視化部分問題詳解

1、數(shù)據(jù)可視化（即畫圖）部分作業(yè)x a cost1.已知橢圓的長、短軸a 4,b2，用“小紅點線”畫橢圓。(參見圖P5-1)y bs int程序：t=2*pi*(0:)/;a=4;b=2;x=a*cos(t);y=b*si n(t);plot(x,y,'r.','MarkerSize',15)axis equal-3-2-101234圖 p5-12.根據(jù)表達式 1 cos繪制如圖p5-2的心臟線。（提示：采用極坐標繪線指令 polar） clft = 0:.01:2*pi; P=1-cos(t);pli ne=polar(t,P,'r'), set(

2、pli ne,'Li neWidth',5)title('P=1-costheta')圖 p5-23. A,B,C三個城市上半年每個月的國民生產(chǎn)總值見表p5.1。試畫出如圖p5-3所示的三城市上半年每月生產(chǎn)總值的累計直方圖。（提示：使用指令bar）表p5.1各城市生產(chǎn)總值數(shù)據(jù)（單位：億元）城市1月2月3月4月5月6月A170120180200190220B120100110180170180C70508010095120clearx=1:6;Y=170 120180200190220;120100110180170180;70508010095120bar(x,

3、Y','grouped');bar(x，丫 ','stacked');colormap(cool);% lege nd('A','B',C,'Locatio n'，Northwest') lege nd('A','B',C,2)圖 p5-34.二階線性系統(tǒng)的歸一化（即令n 1 ）沖激響應(yīng)可表示為：=e t sin( t)01y(t)te t1 ，其中,|12 | ,為阻尼1e (" e()t1系數(shù)。試在同一圖上，畫出0.2 : 0.2 : 1.4不同

4、取值時的各條曲線，時間區(qū)間為t 0,18。此外，要求（1）1的各條曲線用藍色,1的用紅色，1的用黑要求圖形參見圖p5-4)% exmp504.m clc,clf,clear;t=(0:0.05:18)：N=le ngth(t); zeta=0.2:0.2:1.4;L=le ngth(zeta);y=zeros(N,L);hold onfor k=1:Lzk=zeta(k);beta=sqrt(abs(1-zkA2);if zk-1<-2*epsy=1/beta*exp(-zk*t).*si n(beta*t); plot(t,y,'b') if zk<0.4text

5、(2.2,0.63,'zeta = 0.2')endelseif (zk-1)<2*eps%y=t.*exp(-t);plot(t,y,'k','Li neWidth',2)elsey=(exp(-(zk-beta)*t)-exp(-(zk+beta)*t)/(2*beta); plot(t,y,'r') if zk>1.2text(0.3,0.14,'zeta = 1.4')endendendtext(10,0.7,'Deltazeta=0.2')axis(0,18,-0.4,0.8)

6、hold offbox ongrid on圖 p5-45.用綠實線繪制x sin(t) , y cos(t), z t的三維曲線，曲線如圖p5-5所示。(提示:使用plot3指令)elft=4*pi*(O:1OO)/1OO;x=si n(t);y=cos(t);z=t;plot3(x,y,z,'g','Li neWidth',3),box on圖 p5-56.在區(qū)域x, y 3,3，繪制z24xe x的如圖p5-6的三維(透視)網(wǎng)格曲面。(不得使用 ezmesh)elfx=-3:0.1:3;y=x;X,Y=meshgrid(x,y);Z=4.*X.*exp(-X

7、.A2-Y.A2); mesh(X,Y,Z) hidde n off% colormap(cool),% shad ing in terp,syms x y z% z=4.*x.*exp(-x.A2-y.A2); z=4*x*exp(-xA2-yA2); ezmesh(z,-3,3)hidde n of圖 p5-67在x, y 4 ,4 區(qū)間里，根據(jù)表達式z泄，繪制如圖p5-7所示的曲面。x yclear allx=4*pi*(_50：50)/50；y=x; X,Y=meshgrid(x,y);Z=si n(X+Y)./(X+Y+(X+Y=O)*eps);surf(X,Y,Z)view(21,

8、32)shad ing in terp圖 p5-7&試用圖解法回答:(1)方程組廠x2不0.1有多少個實數(shù)解?sin(x cos(y) 0(2)求出離x 0, y0最近、且滿足該方程組的一個近似解。(1)ezplot('y/(1+xA2+yA2)-0.1',-2*pi,2*pi,-pi/3,3.5*pi)hold onezplot('si n(x+cos(y)',-2*pi,2*pi,-pi/3,3.5*pi)圖 p5-8可看到6個交點，即方程組有6個實數(shù)解(2)edit sushufunction F=myfun 8(x,y)F=y/(1+xA2+yA2)-0.1;si n(x+cos(y)end zoom onxy= gin put(1) f=fsolve(myfu n8,xy) X,Y = gin put(1) f=fs