數(shù)據(jù)可視化部分答案

1、數(shù)據(jù)可視化(即畫(huà)圖)部分作業(yè)1 已知橢圓的長(zhǎng)、短軸，用“小紅點(diǎn)線”畫(huà)橢圓。（參見(jià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)axis equal圖 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-23 A,B,C三個(gè)城市上半年每個(gè)月的國(guó)民生

2、產(chǎn)總值見(jiàn)表p5.1。試畫(huà)出如圖p5-3所示的三城市上半年每月生產(chǎn)總值的累計(jì)直方圖。（提示：使用指令bar）表p5.1 各城市生產(chǎn)總值數(shù)據(jù)（單位：億元）城市1月2月3月4月5月6月A170120180200190220B120100110180170180C70508010095120 clearx=1:6;Y=170120180200190220; 120100110180170180; 70508010095120 bar(x,Y,grouped); bar(x,Y,stacked);colormap(cool); % legend(A,B,C,Location,NorthWest)lege

3、nd(A,B,C,2)圖 p5-34 二階線性系統(tǒng)的歸一化（即令）沖激響應(yīng)可表示為：，其中，為阻尼系數(shù)。試在同一張圖上，畫(huà)出不同取值時(shí)的各條曲線，時(shí)間區(qū)間為。此外，要求（1）的各條曲線用藍(lán)色，的用紅色，的用黑色，且線較粗；（2）和的兩條曲線要有醒目的阻尼系數(shù)標(biāo)志。（本題要求圖形參見(jiàn)圖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);hold onfor k=1:L zk=zeta(k); beta=sqrt(abs(1-zk2);if zk-1

4、-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) endelseif (zk-1)1.2 text(0.3,0.14,zeta = 1.4) end endendtext(10,0.7,Deltazeta=0.2)axis(0,18,-0.4,0.8)hold offbox ongrid on 圖 p5-45 用綠實(shí)線繪制，的三維曲線，曲線如圖p5-5所示。（提示：使用plot3指令）clft=4*pi*(0:100)/100; x=sin(t);y=cos(t);z=

5、t; plot3(x,y,z,g,LineWidth,3),box on圖 p5-56 在區(qū)域，繪制的如圖p5-6的三維（透視）網(wǎng)格曲面。（不得使用ezmesh） clfx=-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(-x2-y2);ezmesh(z,-3,3)hidden off圖 p5-67 在區(qū)間里，根據(jù)表達(dá)式，繪制如圖p5-7所示

6、的曲面。 clear allx=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) shading interp圖 p5-78 試用圖解法回答：（1）方程組有多少個(gè)實(shí)數(shù)解？（2）求出離最近、且滿足該方程組的一個(gè)近似解。 （1） ezplot(y/(1+x2+y2)-0.1,-2*pi,2*pi,-pi/3,3.5*pi) hold on ezplot(sin(x+cos(y),-2*pi,2*pi,-pi/3,3.5*pi)圖 p5-8 可看到6個(gè)交點(diǎn)，即方程組有6個(gè)實(shí)數(shù)解（2） edit sushufunction F=myfun8(x,y)F=y/(1+x2+y2)-0.1;sin(x+cos(y)endzoom onxy= ginput(1)f=fsolve(myfun8,xy)X,Y = ginput(1)f=fsolv