資源描述：
《數(shù)據(jù)可視化部分答案.doc》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在教育資源-天天文庫。
1、數(shù)據(jù)可視化(即畫圖)部分作業(yè)1．已知橢圓的長、短軸，用“小紅點(diǎn)線”畫橢圓。（參見圖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)axisequal圖p5-12．根據(jù)表達(dá)式繪制如圖p5-2的心臟線。（提示：采用極坐標(biāo)繪線指令polar）clft=0:.01:2*pi;P=1-cos(t);pline=polar(t,P,'r'),set(pline,'LineWidth',5)title('P=1-costheta')圖p5-21．A,B,C三個(gè)城市
2、上半年每個(gè)月的國民生產(chǎn)總值見表p5.1。試畫出如圖p5-3所示的三城市上半年每月生產(chǎn)總值的累計(jì)直方圖。（提示：使用指令bar）表p5.1各城市生產(chǎn)總值數(shù)據(jù)（單位：億元）城市1月2月3月4月5月6月A170120180200190220B120100110180170180C70508010095120clearx=1:6;Y=[170120180200190220;120100110180170180;70508010095120]bar(x,Y','grouped');bar(x,Y','stacked');colormap(cool);%legend('A','B','C
3、','Location','NorthWest')legend('A','B','C',2)圖p5-31．二階線性系統(tǒng)的歸一化（即令）沖激響應(yīng)可表示為：，其中，為阻尼系數(shù)。試在同一張圖上，畫出不同取值時(shí)的各條曲線，時(shí)間區(qū)間為。此外，要求（1）的各條曲線用藍(lán)色，的用紅色，的用黑色，且線較粗；（2）和的兩條曲線要有醒目的阻尼系數(shù)標(biāo)志。（本題要求圖形參見圖p5-4）%exmp504.mclc,clf,clear;t=(0:0.05:18)';N=length(t);zeta=0.2:0.2:1.4;L=length(zeta);y=zeros(N,L);holdonfork=1:L
4、zk=zeta(k);beta=sqrt(abs(1-zk^2));ifzk-1<-2*epsy=1/beta*exp(-zk*t).*sin(beta*t);plot(t,y,'b')ifzk<0.4text(2.2,0.63,'zeta=0.2')endelseif(zk-1)<2*eps%y=t.*exp(-t);plot(t,y,'k','LineWidth',2)elsey=(exp(-(zk-beta)*t)-exp(-(zk+beta)*t))/(2*beta);plot(t,y,'r')ifzk>1.2text(0.3,0.14,'zeta=1.4')en
5、dendendtext(10,0.7,'Deltazeta=0.2')axis([0,18,-0.4,0.8])holdoffboxongridon圖p5-41．用綠實(shí)線繪制，，的三維曲線，曲線如圖p5-5所示。（提示：使用plot3指令）clft=4*pi*(0:100)/100;x=sin(t);y=cos(t);z=t;plot3(x,y,z,'g','LineWidth',3),boxon圖p5-52．在區(qū)域，繪制的如圖p5-6的三維（透視）網(wǎng)格曲面。（不得使用ezmesh）clfx=-3:0.1:3;y=x;[X,Y]=meshgrid(x,y);Z=4.*X.
6、*exp(-X.^2-Y.^2);mesh(X,Y,Z)hiddenoff%colormap(cool),%shadinginterp,symsxyz%z=4.*x.*exp(-x.^2-y.^2);z=4*x*exp(-x^2-y^2);ezmesh(z,[-3,3])hiddenoff圖p5-61．在區(qū)間里，根據(jù)表達(dá)式，繪制如圖p5-7所示的曲面。clearallx=4*pi*(-50:50)/50;y=x;[X,Y]=meshgrid(x,y);Z=sin(X+Y)./(X+Y+(X+Y==0)*eps);surf(X,Y,Z)view([21,32])shadingi
7、nterp圖p5-71．試用圖解法回答：（1）方程組有多少個(gè)實(shí)數(shù)解？（2）求出離最近、且滿足該方程組的一個(gè)近似解。（1）ezplot('y/(1+x^2+y^2)-0.1',[-2*pi,2*pi,-pi/3,3.5*pi])holdonezplot('sin(x+cos(y))',[-2*pi,2*pi,-pi/3,3.5*pi])圖p5-8可看到6個(gè)交點(diǎn)，即方程組有6個(gè)實(shí)數(shù)解（2）editsushufunctionF=myfun8(x,y)F=[y/(1+x^2+y^2)-0.1;sin(x+cos
