1、1 已知椭圆的长、短轴 ，用“小红点线”画如图 p5-1 所示的椭2,4ba圆 。（提示：参量 ；点的大小；axis equal）tbyxsincot-4 -3 -2 -1 0 1 2 3 4-3-2-10123xy图 p5-1cleart=2*pi*(0:149)/149; a=4;b=2;x=a*cos(t);y=b*sin(t); plot(x,y,r.,MarkerSize,15)axis equal -3 -2 -1 0 1 2 3 4-3-2-101232 根据表达式 绘制如图 p5-2 的心脏线。（提示： polar；注意 title 中特殊字符；线cos宽；axis squar
2、e。可以用 plot 试试。）0.511.52302106024090270120300150330180 0=1-cos图 p5-2cleart = 0:.01:2*pi; P=1-cos(t); f=polar(t,P,r);axis square;set(f,LineWidth,4)title(P=1-costheta) 0.511.52302106024090270120300150330180 0P=1-cos3 A,B,C 三个城市上半年每个月的国民生产总值表 p5.1。试画出如图 p5-3 所示的三城市上半年每月生产总值的累计直方图。（提示：bar(x,Y,style); col
3、ormap(cool); legend。）表 p5.1 各城市生产总值数据（单位：亿元）城市 1 月 2 月 3 月 4 月 5 月 6 月A 170 120 180 200 190 220B 120 100 110 180 170 180C 70 50 80 100 95 1201 2 3 4 5 60100200300400500600ABC图 p5-3clearX=(1:6);Y=170,120,180,200,190,220;120,100,110,180,170,180;70,50,80,100,95,120;bar(X,Y,stacked);colormap(cool);legen
4、d(A,B,C,2) 1 2 3 4 5 60100200300400500600ABC4 二阶线性系统的归一化（即令 ）冲激响应可表示为：1n，其中 ， 为阻尼系1021) si()( )()( tttetty |1|2数。（1）希望在同一张图上，绘制 区间内 不同取值时的各条曲8,0t 4.:02.线（参见图 p5-4）。在此图上， 的各条曲线为细蓝线 ； 为粗黑线； 为细红11线；并且对最上方及最下方的两条曲线给出 和 的醒目标志。（2）读者运行题.下程序 exmp504.m，可以发现该程序画出的曲线中没有“粗黑线”。你能讲出原因吗？如何对exmp504.m 作最少的修改（比如只改一条指
5、令），就可画出所需图形。（提示：该题深层次地暴露数值计算可能存在的隐患。）0 2 4 6 8 10 12 14 16 18-0.4-0.200.20.40.60.8 = 0.2 = 1.4=0.2图 p5-4% exmp504.m 供第 4 道习题使用的程序clc,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:Lzk=zeta(k);beta=sqrt(abs(1-zk2);if zk-11.2text(0.3,0.14,zeta = 1.4)end
6、endendtext(10,0.7,Deltazeta=0.2)axis(0,18,-0.4,0.8)hold offbox ongrid on 0 2 4 6 8 10 12 14 16 18-0.4-0.200.20.40.60.8 = 0.2 = 1.4=0.25 用绿实线绘制 ， ， 的三维曲线，曲线如图 p5-5 所示。（提示：参)sin(tx)cos(tytz变量；plot3；线色线粗。）-1 -0.5 00.5 1-101051015图 p5-5cleart=(0:0.02:4)*pi;x=sin(t);y=cos(t);z=t;plot3(x,y,z,g-,LineWidth,
7、3)box on -1 -0.50 0.51-1-0.500.510510156 采用两种不同方法绘制 在 的如图 p5-6 的三维（透视）网格曲24yxez3,面。（提示：ezmesh; mesh; hidden）图 p5-6方法一：x=-3:0.1:3;y=x;X,Y=meshgrid(x,y);Z=4*X.*exp(-X.2-Y.2);mesh(X,Y,Z)colormap(jet)axis(-3,3,-3,3,-2,2)hidden off 方法二：clearsyms x y z;z=4*x*exp(-x2-y2);ezmesh(z,-3,3)colormap(jet)axis(-3,
8、3,-3,3,-2,2)hidden off 7 在 区间里，根据表达式 ，绘制如图 p5-7 所示的曲面。（提4,yx yxz)sin(示：NaN 的处理）图 p5-7clear x=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(20,18)shading interptitle(z=sin(x+y)/(x+y) 8 试用图解法回答：（1） 方程组 有多少个实数解？（提示：图解法; ezplot; ginput）0)cos(in1.1(2yx（2） 求出离 最近、且满足
9、该方程组的一个近似解。,0clearezplot(y/(1+x2+y2)-0.1,-2*pi,2*pi,-pi/4,4*pi)hold onezplot(sin(x+cos(y),-2*pi,2*pi,-pi/4,4*pi)grid onxysin(x+cos(y) = 0-6 -4 -2 0 2 4 6024681012由图可知，共有 6 个交点,故方程组共有六组实数解求离 最近、且满足该方程组的一个近似解0,yxx,y=ginput(1) x =-0.9797y =0.2002xysin(x+cos(y) = 0-0.98 -0.9799 -0.9798 -0.9797 -0.9796 -
10、0.9795 -0.97940.19980.19990.20.20010.20020.20030.20040.20050.20069 制作如文件 prob509.p（在光盘的 mfiles 文件夹上）运行时那样的色图变幻（参见图 p5-8）。（提示：jet; flipud(jet); colormap; spinmap）图 p5-8clearX,Y,Z=sphere(30);colormap(jet)surf(X,Y,Z)axis off,axis equal,shading interplight(position,0 -10 1.5,style,infinite)lighting phong,material shinylight;lighting flat,set(gcf,Color,w)view(-160,30),shg,C=jet;CC=C;flipud(C);colormap(CC)disp(按任意键，观察图像变换)pause,spinmap(40,8) 按任意键，观察图像变换
