第5章--MATLAB绘图-习题答案
第5章 MATLAB绘图_习题答案
第5章 MATLAB绘图习题5一、选择题1．如果x、y均为4×3矩阵，则执行plot(x,y)命令后在图形窗口中绘制（）条曲线。
DA．12 B．7 C．4 D．32．下列程序的运行结果是（）。
Ax=0:pi/100:2*pi;for n=1:2:10plot(n*sin(x),n*cos(x))hold onendaxis squareA．5个同心圆B．5根平行线C．一根正弦曲线和一根余弦曲线D．5根正弦曲线和5根余弦曲线3．命令text(1,1,'{\alpha}+{\beta}')执行后，得到的标注效果是（）。
CA．{\alpha}+{\beta} B．{\α}+{\β} C．α+βD．\α+\β4．subplot(2,2,3)是指（）的子图。
AA．两行两列的左下图B．两行两列的右下图C．两行两列的左上图D．两行两列的右上图5．要使函数y=2e x的曲线绘制成直线，应采用的绘图函数是（）。
CA．polar B．semilogx C．semilogy D．loglog6．下列程序的运行结果是（）。
B[x,y]=meshgrid(1:5);surf(x,y,5*ones(size(x)));A．z=x+y平面B．与xy平面平行的平面C．与xy平面垂直的平面D．z=5x平面7．下列函数中不能用于隐函数绘图的是（）。
DA．ezmesh B．ezsurf C．ezplot D．plot38．下列程序运行后，看到的图形（）。
Ct=0:pi/20:2*pi;[x,y]=meshgrid(-8:0.5:8);z=sin(sqrt(x.^2+y.^2))./sqrt(x.^2+y.^2+eps);surf(x,y,z)view(0,90);axis equalA．像墨西哥帽子B．是空心的圆C．边界是正方形D．是实心的圆9．下列程序运行后得到的图形是（）。
A[x,y]=meshgrid(-2:2);z=x+y;i=find(abs(x)<1 & abs(y)<1);z(i)=NaN;surf(x,y,z);shading interpA．在一个正方形的正中心挖掉了一个小的正方形B．在一个正方形的正中心挖掉了一个小的长方形C．在一个正方形的上端挖掉了一个小的正方形D．在一个正方形的下端挖掉了一个小的正方形10．在使用MA TLAB“绘图”选项卡中的命令按钮绘图之前，需要（）。
MATLAB语言基础与应用(第二版)第5章 习题答案
第5章习题与答案5.1用矩阵三角分解方法解方程组123123123214453186920x x x x x x x x x +-=⎧⎪-+=⎨⎪+-=⎩ 解答：>>A=[2 1 -1;4 -1 3;6 9 -1] A =2 1 -1 4 -13 6 9 -1 >>b=[14 18 20]; b =14 18 20 >> [L, U, P]=lu(A) L =1.0000 0 0 0.6667 1.0000 0 0.3333 0.2857 1.0000 U =6.0000 9.0000 -1.0000 0 -7.0000 3.6667 0 0 -1.7143 P =0 0 1 0 1 0 1 0 0 >> y=backsub(L,P*b’) y =20.0000 4.6667 6.0000 >> x=backsub(U,y) x =6.5000 -2.5000 -3.5000 5.2 Cholesky 分解方法解方程组123121332352233127x x x x x x x ++=⎧⎪+=⎨⎪+=⎩ 解答：>> A=[3 2 3;2 2 0;3 0 12] A =3 2 32 2 03 0 12>> b=[5;3;7]b =537>> L=chol(A)L =1.7321 1.1547 1.73210 0.8165 -2.44950 0 1.7321>> y=backsub(L,b)y =-11.6871 15.7986 4.0415>> x=backsub(L',y)x =-6.7475 28.8917 49.93995.3解答:观察数据点图形>> x=0:0.5:2.5x =0 0.5000 1.0000 1.5000 2.0000 2.5000 >> y=[2.0 1.1 0.9 0.6 0.4 0.3]y =2.0000 1.1000 0.9000 0.6000 0.4000 0.3000 >> plot(x,y)图5.1 离散点分布示意图从图5.1观察数据点分布，用二次曲线拟合。
Get清风MATLAB教程a第5章习题解答张志涌
MATLAB教程2012a第5章习题解答-张志涌第5章 数据和函数的可视化习题5及解答1 椭圆的长、短轴2,4==b a ，用“小红点线〞画椭圆⎩⎨⎧==tb y ta x sin cos 。
〔参见图p5-1〕〖解答〗 clf a=4;b=2;t=0:pi/80:2*pi; x=a*cos(t); y=b*sin(t);plot(x,y,'r.','MarkerSize',15) axis equal xlabel('x') ylabel('y')shg-4-3-2-101234-3-2-1123xy2 根据表达式θρcos 1-=绘制如图p5-2的心脏线。
〔提示：采用极坐标绘线指令polar 〕〖解答〗 clftheta=0:pi/50:2*pi;rho=1-cos(theta);h=polar(theta,rho,'-r');%极坐标绘线指令。
h 是所画线的图柄。
set(h,'LineWidth',4) %利用set 设置h 图形对象的“线宽〞axis square %保证坐标的圆整性0.51 1.523021060240902701203001503301800ρ=1-cos θ3 A,B,C 三个城市上半年每个月的国民生产总值如见表p5.1。
试画出如图p5-3所示的三城市上半年每月生产总值的累计直方图。
表p5.1 各城市生产总值数据〔单位：亿元〕城市 1月 2月 3月 4月 5月 6月 A 170 120 180 200 190 220 B 120 100 110 180 170 180 C 70508010095120〖目的〗● 借助MATLAB 的帮助系统，学习直方图指令polar 的使用。
● bar 指令常用格式之一：bar(x,Y,'style') 。
x 是自变量列向量；Y 是与x 行数相同的矩阵，Y 的每一行被作为“一组〞数据；style 取stacked 时，同一组数据中每个元素对应的直方条被相互层叠。
第5章MATLAB绘图_习题答案
第5章MATLAB绘图习题5一、选择题1．如果x、y均为4×3矩阵，则执行plot(x,y)命令后在图形窗口中绘制（）条曲线。
DA．12B．7C．4D．32．下列程序的运行结果是（）。
Ax=0:pi/100:2*pi;forn=1:2:10plot(n*sin(x),n*cos(x))holdonendaxissquareA．5个同心圆B．5根平行线C．一根正弦曲线和一根余弦曲线D．5根正弦曲线和5根余弦曲线3．命令text(1,1,'{\alpha}+{\beta}')执行后，得到的标注效果是（）。
C A．{\alpha}+{\beta}B．αβ}C．α+βD．αβ4．subplot(2,2,3)是指（）的子图。
AA．两行两列的左下图B．两行两列的右下图C．两行两列的左上图D．两行两列的右上图x的曲线绘制成直线，应采用的绘图函数是（）。
C5．要使函数y=2eA．polarB．semilogxC．semilogyD．loglog6．下列程序的运行结果是（）。
B[x,y]=meshgrid(1:5);surf(x,y,5*ones(size(x)));A．z=x+y平面B．与xy平面平行的平面C．与xy平面垂直的平面D．z=5x平面7．下列函数中不能用于隐函数绘图的是（）。
DA．ezmeshB．ezsurfC．ezplotD．plot38．下列程序运行后，看到的图形（）。
Ct=0:pi/20:2*pi;[x,y]=meshgrid(-8:0.5:8);z=sin(sqrt(x.^2+y.^2))./sqrt(x.^2+y.^2+eps);surf(x,y,z)view(0,90);axisequalA．像墨西哥帽子B．是空心的圆C．边界是正方形D．是实心的圆9．下列程序运行后得到的图形是（）。
A[x,y]=meshgrid(-2:2);z=x+y;i=find(abs(x)<1&abs(y)<1);z(i)=NaN;surf(x,y,z);shadinginterpA．在一个正方形的正中心挖掉了一个小的正方形B．在一个正方形的正中心挖掉了一个小的长方形C．在一个正方形的上端挖掉了一个小的正方形D．在一个正方形的下端挖掉了一个小的正方形10．在使用MATLAB“绘图”选项卡中的命令按钮绘图之前，需要（）。
MATLAB基础习题第五章习题答案
case{-1,2}
disp('您赢了');
end
%%
%是否继续判别模块
n=input('是否继续玩该游戏？否（0），是（1）');
while n~=0&n~=1
disp('您输入的不是正确数字，请正确输入');
n=input('是否继续玩该游戏？否（0），是（1）');
end
（3）程序要具有友好性
答：
%%
%该模块实现的功能是：求解六元线性方程组，方程组的系数矩阵由用户通过键盘输入；得到系数矩阵后给出方程的解；
clc;
clear all;
close all;
%%
%方程输入模块
x=inputdlg({'第一个方程系数','第二个','第三个','第四个','第五个','第六个'});
end
>>
请选择，石头（1），剪刀（2），布（3）：1
您的选择是：
石头
电脑的选择是：
、剪
您赢了
是否继续玩该游戏？否（0），是（1）1
请选择，石头（1），剪刀（2），布（3）：2
您的选择是：
、剪
电脑的选择是：
石头
您输了
是否继续玩该游戏？否（0），是（1）0
>>
5.编写一个日程提醒程序实现如下功能：
（1）如果当前时间为7-8点则提醒用户，该吃早饭了；
end
end
%%
%%显示模块
disp('90分以上的人有： Nhomakorabea)disp(x)
matlab课后答案完整版
matlab课后答案完整版ones表⽰1矩阵zeros表⽰0矩阵ones(4)表⽰4x4的1矩阵zeros(4)表⽰4x4的0矩阵zeros(4,5)表⽰4x5的矩阵eye(10,10)表⽰10x10的单位矩阵rand(4,5)表⽰4x5的伴随矩阵det(a)表⽰计算a的⾏列式inv(a)表⽰计算a的逆矩阵Jordan(a)表⽰求a矩阵的约当标准块rank(a)表⽰求矩阵a的秩[v,d]=eig(a)对⾓矩阵b=a’表⽰求a矩阵的转置矩阵sqrt表⽰求平⽅根exp表⽰⾃然指数函数log⾃然对数函数abs绝对值第⼀章⼀、5（1）b=[97 67 34 10;-78 75 65 5;32 5 -23 -59]; >> c=[97 67;-78 75;32 5;0 -12]; >> d=[65 5;-23 -59;54 7];>> e=b*ce =5271 11574-11336 6641978 3112（2）a=50:1:100⼆、1 、x=-74;y=-27;z=(sin(x.^2+y.^2))/(sqrt(tan(abs(x+y)))+pi) z =-0.09012、a=-3.0:0.1:3.0;>> b=exp(-0.3*a).*sin(a+0.3)y =0.7218 1.0474-0.2180 1.15624、a*b表⽰a矩阵和b矩阵相乘a.*b表⽰a矩阵和b矩阵单个元素相乘A(m,n)表⽰取a矩阵第m⾏，第n列A(m,:)表⽰取a矩阵第m⾏的全部元素A(:,n)表⽰取a矩阵的第n列全部元素A./B表⽰a矩阵除以b矩阵的对应元素，B.\A等价于A./BA.^B表⽰两个矩阵对应元素进⾏乘⽅运算A.^2表⽰a中的每个元素的平⽅A^2表⽰A*A例：x=[1,2,3];y=[4,5,6];z=x.^yz=1 32 729指数可以是标量（如y=2）.底数也可以是标量（如x=2）5、a=1+2i;>> b=3+4i;>> c=exp((pi*i)/6)c =0.8660 + 0.5000id=c+a*b/(a+b)d =1.6353 + 1.8462i第⼆章⼆、4、（1）y=0;k=0;>> while y<3k=k+1;>> display([k-1,y-1/(2*k-1)])ans =56.0000 2.9944第三章⼆1（1） x=0:pi/10:2*pi; >> y=x-x.^3/6; >> plot(x,y)1234567-40-35-30-25-20-15-10-505（2）x=0:pi/10:2*pi; y=(exp(-x.^2/2))/2*pi;plot(x,y)012345670.20.40.60.811.21.41.6（3）x=-8:0.01:8; y=sqrt((64-x.^2)/2);plot(x,y)-8-6-4-2024680123456（4）t=0:0.1:8*pi; >> x=t.*sin(t); >> y=t.*cos(t);-25-20-15-10-50510152025-30-20-10102030例3.4x=0:pi/100:2*pi; y1=exp(-0.5*x);y2=exp(-0.5*x).*sin(2*x); plot(x,y1,x,y2)>> title('x from 0 to 2{\pi} '); >> xlabel('variable x'); >> ylabel('variable y'); >> text(1.5,0.5,'曲线y1=e^(-0.5x)'); >> text(3,0.1,'曲线y2=cos(4{\pi}x)e^{-0.5x}'); >> legend('y1','y2')1234567-0.4-0.20.20.40.60.81x from 0 to 2πvariable xv a r i a b l e y曲线y1=e (-0.5x)曲线y2=cos(4πx)e -0.5xy1y22、(1)y1=2*x-0.5;t=linspace(0,pi,100); x=sin(3*t).*cos(t); y=sin(3*t).*sin(t);>> k=find(abs(y-x)<1e-2); >> t1=t(k) t1 =0 0.7933 1.04722.0944>> z=sin(3.*(t1)).*cos(t1) z =0 0.4841 0.0000 0.0000 -0.0000>> plot(t,x,t,y,'k:',t1,z,'bp');0.511.522.533.5-1-0.8-0.6-0.4-0.200.20.40.60.81（2）subplot(1,2,1); >> scatter(x1,y1,10)； >> title('y=2x-0.5')； >> subplot(1,2,2); >> scatter(x,y,10)-1-0.8-0.6-0.4-0.200.20.40.60.81-1-0.8-0.6-0.4-0.200.20.40.63、subplot(1,2,1); x=0:0.01:pi; y=sin(1./x); plot(x,y)subplot(1,2,2);fplot('sin(1./x)',[1,100])1234-1-0.8-0.6-0.4-0.200.20.40.60.81204060801000.10.20.30.44、t=0:pi:2*pi; y=1./(1+exp(-t));subplot(2,2,1);%图形窗⼝的分割bar(t,'group'); %绘制柱形图（分组） subplot(2,2,2);barh(t,'stack');%绘制柱形图（堆积） subplot(2,2,3);loglog(t,y); %函数使⽤全对数坐标，x,y 均采⽤常⽤对数刻度 subplot(2,2,4); semilogy(t,y); %函数使⽤半对数坐标，y 轴为常⽤对数刻度，x 轴仍为线性刻度1230246802468123100.5100.710-0.01810-0.0010246810-0.310-0.210-0.15、（1）theta=linspace(-pi,pi,100); ro=5.*cos(theta)+4; polar(theta,ro); （2）x=linspace(0,2*pi,100);a=1>> r=a.*(1+cos(x)); polar(x,r);3021060240902701203001503301806、（1）t=0:pi/10:2*pi;>> x=exp((-t)/20).*cos(t); >> y=exp((-t)/20).*sin(t); >> z=t; >> plot3(x,y,z);-1-0.50.51-1-0.50.5102468（2）t=0:0.01:1; x=t;>> y=t.^2; >> z=t.^3;>> plot3(x,y,z);0.20.40.60.800.20.40.60.817、x=-30:0.1:0; >> y=0:0.1:30;>> [x,y]=meshgrid(x,y); >>z=10.*sin(sqrt(x.^2+y.^2))./sqrt(1+x.^2+y.^2);>> meshc(x,y,z);绘制曲⾯图和等⾼线-30-20-10102030-4-202468、x=linspace(-3,3,100); >> y=linspace(-3,3,100); >> [x y]=meshgrid(x,y); %可以将向量转化为矩阵 >> fxy=-5./(1+x.^2+y.^2); >> i=find(abs(x)<=0.8 & abs(y)<=0.5); >> fxy(i)=NaN; >>surf(x,y,fxy) %绘制三维曲⾯图-4-224-4-224-4-3-2-19、u=linspace(1,10,100); v=linspace(-pi,pi,100);[u v]=meshgrid(u,v); x=3.*u.*sin(v); y=2.*u.*cos(v); z=4*u.^2; surf(x,y,z); shading interp;-40-20-1010200100200300400第五章⼆1、a=rand(1,30000);mean(a) %求平均数 ans =0.5010 >>b=std(a) %求标准差 b =0.2882 >> c=max(a) c =0.9999 >> d=min(a) d =3.5706e-005size(find(a>0.5))/size(a) %求⼤于0.5的随机数个数占总数的百分⽐ans =0.50322、h=[466,715,950,1422,1635]; >> w=[7.04,4.28,3.40,2.52,2.13]; >> hh=[500,900,1500]; >> ww=interp1(h,w,hh,'spline')ww =6.4903 3.5226 2.3845 3、x=linspace(1,10,50); y=log(x);f=polyfit(x,y,5); %求曲线的拟合 >> yy=polyval(f,x); >> plot(x,y,'r-',x,yy,'g.') 123456789100.511.522.55、（1）、（2） p1=[1,2,0,7]; p2=[1,-2]; p3=[1,0,5,1]; p12=conv(p1,p2); >>p=p12+[zeros(1,size(p12,2)-size(p3,2)),p3]; >> roots(p) ans =-3.4656 0.6128 + 1.6278i 0.6128 - 1.6278i 1.2400-29 291 95 19 -3 697 -13 697 1427 >>y2=polyvalm(p,a)%以矩阵a 为⾃变量 y2 =391 2084 3273 502 2693 4207 720 3775 5892 6、（1）z=fzero('3*x-sin(x)+1',0) %求x=0时附近的根 z =-0.4903 第⼋章⼆、2t=0:pi/20:2*pi; x=sin(t); y=cos(t); x1=sin(7*t); y1=cos(7*t);h=plot(x,y,x1,y1);set(h,'marker','x','linewidth',2); set(gca,'xtick',-1:0.1:1); title('篮筐')-1-0.9-0.8-0.7-0.6-0.5-0.4-0.3-0.2-0.100.10.20.30.40.50.60.70.80.91 -1-0.8-0.6-0.4-0.200.20.40.60.81篮筐3、x=0:pi/10:5*pi;y=exp(-0.2*x).*cos(x)+2; h=plot(x,y);set(gca,'color','red','linestyle','-','linewidth',3);text(5,2.4,'y=exp(-0.2*x).*cos(x)+2');02468101214161.41.61.822.22.42.62.83y=exp(-0.2*x).*cos(x)+24、t=-pi:pi/100:pi; x=cos(t); y=sin(t); z=t;h=plot(t,x,t,y,t,z);set(h,'linestyle','-','linewidth',3);-4-3-2-101234-4-3-2-101234字符串例ch='Welcome to Beijing';subch=ch(12:18) 选12~18个字符串（空格也算）ans =WELCOME TO BEIJING >> length(k)统计⼩写字母的个数ans = 14 例：已知y=1-1/2+1/3-1/4.........-1/100求y 的值y=0; >> n=100; >> for i=1:100; y=y+(-1)^(i-1)/i; end>> disp(y)0.6882绘制⼆维曲线图x=0:pi/100:2*pi; >> y1=0.2*exp(-0.5*x).*cos(4*pi*x); >> y2=1.5*exp(-0.5*x).*cos(pi*x); >> plotyy(x,y1,x,y2); 7-0.20.20123456-202绘制三维图像例：x=sint+tcost y=cost-tsint z=tt=0:pi/10:10*pi; x=sin(t)+t.*cos(t); y=cos(t)-t.*sin(t); z=t; plot3(x,y,z); axis([-30 30 -30 30 0 35]); 坐标轴的最⼤值与最⼩值title('line in 3-D space'); 图形的题⽬ >> xlabel('x');ylabel('y');zlabel('z'); 标注坐标>> grid on; 加⽹格线 -30-20-10102030-20205101520253035xline in 3-D spaceyz三维例]2/,0[],,0[,cos sin 22ππ∈∈+=y x y x z [x,y]=meshgrid(0:pi/100:pi,0:pi/100:pi/2);>> z=sin(x.^2)+cos(y.^2);>> mesh(x,y,z);>> axis([0 4 0 1.8 -1.5 1.5]); 012340.511.5-1.5-1-0.500.511.5例3.16t=0:pi/20:2*pi; subplot(1,2,1);[x,y,z]=cylinder(sin(t),30);surf(x,y,z); 绘制三维曲⾯图subplot(1,2,2);>> [x,y,z]=peaks(100);>> mesh(x,y,z); 绘制三维⽹格图-11-10100.20.40.60.81-55-505-10-5510多项式求导例：f(x)=1/x^2+5 p=[1];>> q=[1,0,5];>> [p,q]=polyder(p,q)注：c=conv(a,b) 表⽰a 多项式与b 多项式乘积[p,r]=deconv(a,b) 表⽰a 多项式与b 多项式相除其中p 为商向量 r 为余数向量p=polyder(p) 表⽰求p 的导数 p=poleder(p,q) 表⽰求p 乘以q 的导数[p,q]=poleder(p,q) 表⽰p 除以q 的导数多项式求根例：f(x)=2x^4-12x^3+3x^2+5 p=[2,-12,3,0,5]; >> x=roots(p); >> p=[2,-12,3,0,5]; x=roots(p) 求⽅程f(x)=0的根 x =5.7246 0.8997 -0.3122 + 0.6229i -0.3122 - 0.6229i>> g=poly(x) 已知多项式的根求多项式 g =1.0000 -6.0000 1.5000 -0.00002.5000符号求导例7.3x=a(t-tsint)y=b(1-cost) 求y 对x 的⼀阶导数 syms x y a b t;>> f21=a*(t-sin(t)); >> f22=b*(1-cos(t));>> diff(f22)/diff(f21) 求y 对x 的⼀阶导数ans =b*sin(t)/a/(1-cos(t))注：diff(f1,x,2) 表⽰f1对x 的⼆阶导数diff （f3,x ）表⽰z 对x 的偏导 diff （f3,y ）表⽰z 对y 的偏导求不定积分int(f) 求f 的不定积分 f1=int(f,a,b) 求f 在a ，b 之间的定积分eval （f1）计算积分值符号求极限例7.2 syms x h>> f=(sin(x+h)-sin(x))/h;>> limit(f,h,0) h 趋向于0ans =cos(x)例2f=sym('(1+t/x)^x');limit(f,inf) f趋向于⽆穷ans =exp(t)例3f=sym('x*(sqrt(x^2+1)-x)');limit(f,sym('x'),inf,'left') x 趋向于正⽆穷ans =1/2⼤⼩写ch='Welcome to Beijing';subch=ch(12:18)subch =Beijing>> k=find(ch>='A'&ch<='Z'); ch(k)=ch(k)-('A'-'a');>> char(ch)ans =welcome to beijing>> length(k)ans =2。
第5章MATLAB绘图99930资料
程序如下：
x=(0:pi/100:2*pi)';
y1=2*exp(-0.5*x)*[1,-1];
y2=2*exp(-0.5*x).*sin(2*pi*x);
plot(x,y1,'b:'); axis([0,2*pi,-2,2]); %设置坐标 hold on; %设置图形保持状态
plot(x,y2,'k');
参数的grid命令在两种状态之间进行切换。 box on/off命令控制是加还是不加边框线，不带
参数的box命令在两种状态之间进行切换。
例5.5 绘制分段函数曲线并添加图形标注。
程序如下：
x=linspace(0,10,100);
y=[];
for x0=x
if x0>=8
y=[y,1];
elseif x0>=6
subplot(4,4,8);
%选择4×4个区中的8号区
plot(x,ct);title('cotangent(x)');axis ([0,2*pi,-40,40]);
对图形窗口灵活分割。
5.1.3 绘制二维图形的其他函数 1. 其他形式的线性直角坐标图
在线性直角坐标系中，其他形式的图形有 条形图、阶梯图、杆图和填充图等，所 采用的函数分别是：
当前图形窗口，图形窗口原有图形将不复 存在。若希望在已存在的图形上再继续添 加新的图形，可使用图形保持命令hold。 hold on/off命令控制是保持原有图形还是刷 新原有图形，不带参数的hold命令在两种状 态之间进行切换。
MATLAB第五章(刘卫国编)
4/72
说明： （1）当x,y是同维矩阵时，则以x,y对应列元素为 横、纵坐标分别绘制曲线，曲线条数等于矩阵 的列数。 （2）当x是向量，y是有一维与x同维的矩阵时， 则绘制出多根不同色彩的曲线。曲线条数等于 y矩阵的另一维数，x被作为这些曲线共同的横 坐标。 （3）plot函数最简单的调用格式是只包含一个输 入参数：plot(x)。
16/72
5.1.3 绘制二维图形的其他函数 1. 其他形式的线性直角坐标图 在线性直角坐标系中，其他形式的图形有 条形图、阶梯图、杆图和填充图等，所 采用的函数分别是： bar(x,y,选项) stairs(x,y,选项) stem(x,y,选项) fill(x1,y1,选项1,x2,y2,选项2,…)
18/72
2．极坐标图 polar函数用来绘制极坐标图，其调用格式为： polar(theta,rho,选项) 其中theta为极坐标极角，rho为极坐标矢径，选 项的内容与plot函数相似。 例5.9 绘制ρ=sin(2θ)cos(2θ)的极坐标图。 程序如下： theta=0:0.01:2*pi; rho=sin(2*theta).*cos(2*theta); polar(theta,rho,'k');
8/72
例5.4 用不同标度在同一坐标内绘制曲线y1=e0.5xsin(2πx)及曲线y2=1.5e-0.1xsin(x)。 程序如下： x1=0:pi/100:2*pi; x2=0:pi/100:3*pi; y1=exp(-0.5*x1).*sin(2*pi*x1); y2=1.5*exp(-0.1*x2).*sin(x2); plotyy(x1,y1,x2,y2);
10/72
2. 坐标控制 函数的调用格式为： axis([xmin xmax ymin ymax zmin zmax]) axis函数功能丰富，常用的用法还有： axis equal 纵、横坐标轴采用等长刻度 axis square 产生正方形坐标系(缺省为矩形) axis auto 使用缺省设置 axis off 取消坐标轴 axis on 显示坐标轴 grid on/off命令控制是画还是不画网格线，不带 参数的grid命令在两种状态之间进行切换。 box on/off命令控制是加还是不加边框线，不带 参数的box命令在两种状态之间进行切换。 11/72
实验06_高层绘图操作(第5章)
x=linspace(-5,5,21);
y=linspace(0,10,31);
[X,Y]=meshgrid(x,y);
Z=cos(X).*cos(Y).*exp(-sqrt(X.^2+Y.^2)/4);
subplot(2,1,1);surf(X,Y,Z);
subplot(2,1,2);contour3(X,Y,Z);
subplot(2,2,4);fill(x,y3,'g');%填充图
图形:
3. 绘制分段函数的曲线
已知
在-5≤x≤5区间绘制函数曲线。（注意：本曲线不连续！）
程序：
fplot('(x+sqrt(pi))/exp(2)',[-5,0]);
holdon;
fplot('log(x+sqrt(1+x^2))/2',[0,5]);
x
\0
∅
\rceil
ù
\surd
√
\midபைடு நூலகம்
附参考答案：
实验06高层绘图操作
（第5章MATLAB绘图）
一、实验目的
1. 掌握绘制二维图形的常用函数。
2. 掌握绘制三维图形的常用函数。
3. 掌握绘制图形的辅助操作。
二、实验内容
1. 绘制函数的曲线
设 ，在x=0~2π区间取101点，绘制函数的曲线。
程序：
x=linspace(0,2*pi,101);
例
%一个单位圆，一个复数参数
t=0:0.01:2*pi;
x=exp(i*t);%cos(t)+i*sin(t)
plot(x); axisequal;
MATLAB-第五章
例如: factor： 因式分解
>>factor(x^3-6*x^2+11*x-6)
第五章 MATLAB的符号计算
五、符号运算
1 初等代数运算 （3）符号表达式化简（page48，表3-4）
例如: simplify: 对表达式化简
>>simplify(x^3-6*x^2+11*x-6) >>simplify(sin(x)^2 + cos(x)^2) >>simplify(exp(c*log(sqrt(a+b)))) >>simplify((x^2+5*x+6)/(x+2)) >>simplify(sqrt(16))
第五章 MATLAB的符号计算
三、符号表达式的定义
建立符号表达式有以下2种方法： (1)用sym函数建立符号表达式。 >> f=sym('a*x^2+b*x+c'); (2) 使用已经定义的符号变量组成符号表达式。 >> syms x y a b c >> f=a*x^2+b*x+c (?)利用单引号来生成符号表达式。 >> f='a*x^2+b*x+c'
第五章 MATLAB的符号计算
五、符号运算
1 初等代数运算 （2）符号表达式的加减乘除幂次方运算
例如： >> f1=sym('1/(a-b)'); >> f2=sym('2*a/(a+b)'); >> f3=sym('(a+1)*(b-1)*(a-b)');
第五章matlab习题
已知 h(t)=e-t,t≥0,取 N=64， 对 t 从 0~5 秒采样， 用 FFT 作快速傅里 叶变换，并绘制相应的振幅-频率图。
4.当矩阵 A 为自变量时，求 P(X)的值。其中 A 的值与题 3 中的相同
将 100 个学生的 5 门功课的成绩存入矩阵 P 中，进行以下处理 1． 2． ． 4． 分别求每门课的最高分，最低分及相应学生序号 分别求每门课的平均分和标准方差。 5 门课总分的最高分，最低分及相应学生序号 将 5 门课总分按从大到小的顺序存入 ZCJ 中，相应学生序号存
入 XSXH
有 3 个多项式 P1(X)=X4+2X3+4X2+5，P2(X)=X+2，P3(X)=X2+2X+3，是进 行以下操作： 1. 求 P(X)=P1(X) +P2(X)P3(X) 2. 求 P(X)的根 3. 当 X 取矩阵 A 的每一行元素时，求 P(X)的值。其中 A= −1 1.2 −1.4 0.75 2 3.5 0 5 2.5
第5章 MATLAB绘图_习题答案
第5章 MATLAB绘图习题5一、选择题1．如果x、y均为4×3矩阵，则执行plot(x,y)命令后在图形窗口中绘制（）条曲线。
DA．12 B．7 C．4 D．32．下列程序的运行结果是（）。
Ax=0:pi/100:2*pi;for n=1:2:10plot(n*sin(x),n*cos(x))hold onendaxis squareA．5个同心圆B．5根平行线C．一根正弦曲线和一根余弦曲线D．5根正弦曲线和5根余弦曲线3．命令text(1,1,'{\alpha}+{\beta}')执行后，得到的标注效果是（）。
CA．{\alpha}+{\beta} B．{\α}+{\β} C．α+βD．\α+\β4．subplot(2,2,3)是指（）的子图。
AA．两行两列的左下图B．两行两列的右下图C．两行两列的左上图D．两行两列的右上图5．要使函数y=2e x的曲线绘制成直线，应采用的绘图函数是（）。
CA．polar B．semilogx C．semilogy D．loglog6．下列程序的运行结果是（）。
B[x,y]=meshgrid(1:5);surf(x,y,5*ones(size(x)));A．z=x+y平面B．与xy平面平行的平面C．与xy平面垂直的平面D．z=5x平面7．下列函数中不能用于隐函数绘图的是（）。
DA．ezmesh B．ezsurf C．ezplot D．plot38．下列程序运行后，看到的图形（）。
Ct=0:pi/20:2*pi;[x,y]=meshgrid(-8:0.5:8);z=sin(sqrt(x.^2+y.^2))./sqrt(x.^2+y.^2+eps);surf(x,y,z)view(0,90);axis equalA．像墨西哥帽子B．是空心的圆C．边界是正方形D．是实心的圆9．下列程序运行后得到的图形是（）。
A[x,y]=meshgrid(-2:2);z=x+y;i=find(abs(x)<1 & abs(y)<1);z(i)=NaN;surf(x,y,z);shading interpA．在一个正方形的正中心挖掉了一个小的正方形B．在一个正方形的正中心挖掉了一个小的长方形C．在一个正方形的上端挖掉了一个小的正方形D．在一个正方形的下端挖掉了一个小的正方形10．在使用MA TLAB“绘图”选项卡中的命令按钮绘图之前，需要（）。
MATLAB语言与控制系统仿真-参考答案-第5章
控制系统的时域响应MATLAB 仿真实训实训目的1. 学会利用MATLAB 绘制系统的单位阶跃响应曲线，掌握读取系统动态性能指标的方法；2. 学会利用MATLAB 绘制系统的单位脉冲响应曲线的方法；3. 掌握利用MATLAB 绘制系统的零输入响应曲线的方法；4. 掌握利用MATLAB 绘制系统的一般输入响应曲线的方法；5. 学会通过仿真曲线读取相关信息，并依据有关信息进行系统的时域分析。
实训内容1.编写程序求取下列各系统的单位阶跃响应，完成表5-5并记录相关曲线。
162.316)(21++=s s s G 164.216)(22++=s s s G 166.116)(23++=s s s G 1616)(24++=s s s G 解：>> n1=16; >> d1=[1,,16]; >> sys1=tf(n1,d1); >> step(sys1)>> n2=16; >> d2=[1,,16]; >> sys2=tf(n2,d2); >> step(sys2)>> n3=16;>> d3=[1,,16]; >> sys3=tf(n3,d3); >> step(sys3)>> n4=16;>> d4=[1,1,16]; >> sys4=tf(n4,d4); >> step(sys4)表5-5序号ξn ωm ax cp ts t (%5=∆)计算值实验计算值实验计算值实验值1 42 43 44 4w=4;cmax1=1+exp(-z1*pi/sqrt(1-z1^2)); tp1=pi/(w*sqrt(1-z1^2)); ts1=(z1*w); [cmax1,tp1,ts1] ans =>> z2=;w=4;cmax2=1+exp(-z2*pi/sqrt(1-z2^2)); tp2=pi/(w*sqrt(1-z2^2)); ts2=(z2*w); [cmax2,tp2,ts2] ans =>> z3=; w=4;cmax3=1+exp(-z3*pi/sqrt(1-z3^2)); tp3=pi/(w*sqrt(1-z3^2)); ts3=(z3*w); [cmax3,tp3,ts3] ans =>> z4=; w=4;cmax4=1+exp(-z4*pi/sqrt(1-z4^2)); tp4=pi/(w*sqrt(1-z4^2)); ts4=(z4*w); [cmax4,tp4,ts4] ans =说明：对于二阶欠阻尼系统（10<<ξ），若系统的闭环传递函数为2222)(nn ns s s Φωξωω++=则系统单位阶跃响应的输出最大值21max 1ξξπ--+=ec峰值时间21ξωπ-=n p t调整时间估算值ns t ξω5.3=（以5%为误差带）ns t ξω4.4=（以2%为误差带）2.已知二阶系统的闭环传递函数如下，编程求取系统的单位阶跃响应并完成表5-6，记录相关曲线。
matlab第5章
③返回响应输出值：[y, x, t]=step(num, den)
[y, x, t]=step(num,den,t)
t：仿真时间，由系统模型的特性自动生成。
状态变量x返回为空矩阵。 可省略不写。
3
s 1 【例】已知传递函数模型 G ( s) 2 s s5 绘制单位阶跃响应曲线。
（1）max(A) ：返回一个行向量，向量的第i个元素是矩
阵A的第i列上的最大值。
（2）[Y,U]=max(A) ：返回行向量Y和U，Y表示A中每 列的最大值，U对应每列最大值的行号。
30
31
超调量
C= dcgain (G)
[Y,k]=max(y)
%求系统的稳态值
%求y的峰值及峰值采样点
percentovershoot=100*(Y-C)/C %计算超调量
减法：math operations →subtract
传递函数： continuous →transfer Fcn 输入、输出信号组合：Signal Routing→Mux 示波器模块：sinks→scope
图5.13 simulink模型
20
图5.14 锯齿波参数设置界面
21
模型连好后仿真，仿真结束后双击示波器，
f (t ) e2t 时，系统的零状态响应曲线。 求当输入信号为
9
5.3.1.2 时域响应应用举例
【例5.1】已知系统的闭环传递函数为
1 G( s) 2 s 0.4s 1
求单位阶跃响应曲线和单位斜坡响应曲线。
10
解：MATLAB程序代码如下：
第5章 MATLAB绘图
例5-7 在0≤x≤2区间内，绘制曲线y1=2e-0.5x和 y2=cos(4πx)，并给图形添加图形标注。
程序如下：
x=0:pi/100:2*pi;
y坐标数据。
例5-1 在0≤x≤2区间内，绘制曲线 y=2e-0.5xcos(4πx)
程序如下： x=0:pi/100:2*pi; y=2*exp(-0.5*x).*cos(4*pi*x); plot(x,y)
例5-2 绘制曲线。 程序如下： t=0:0.1:2*pi; x=t.*sin(3*t); y=t.*sin(t).*sin(t); plot(x,y);
例5-9 用fplot函数绘制f(x)=cos(tan(πx))的曲线。 命令如下： fplot('cos(tan(pi*x))',[ 0,1],1e-4)
5.1.7 图形窗口的分割
subplot函数的调用格式为：
subplot(m,n,p)
该函数将当前图形窗口分成m×n个绘图区， 即每行n个，共m行，区号按行优先编号， 且选定第p个区为当前活动区。在每一个绘 图区允许以不同的坐标系单独绘制图形。
第5章 MATLAB绘图 5.1 二维数据曲线图 5.2 其他二维图形 5.3 隐函数绘图 5.4 三维图形 5.5 图形修饰处理 5.6 图像处理与动画制作
5.1 二维数据曲线图 5.1.1 绘制单根二维曲线 plot函数的基本调用格式为：
plot(x,y) 其中x和y为长度相同的向量，分别用于存储x坐标和
(2) 当x,y是同维矩阵时，则以x,y对应列元素为横、 纵坐标分别绘制曲线，曲线条数等于矩阵的列数。
MATLAB课后习题集附标准答案
第2章MATLAB概论1、与其他计算机语言相比较，MA TLAB 语言突出的特点是什么？答：起点高、人机界面适合科技人员、强大而简易的作图功能、智能化程度高、功能丰富，可扩展性强.2、MA TLAB 系统由那些部分组成？答：开发环境、MATLAB数学函数库、MATLAB语言、图形功能、应用程序接口3、安装MATLAB 时，在选择组件窗口中哪些部分必须勾选，没有勾选的部分以后如何补安装？答：在安装MATLAB时，安装内容由选择组件窗口中各复选框是否被勾选来决定，可以根据自己的需要选择安装内容，但基本平台（即MATLAB选项）必须安装. 第一次安装没有选择的内容在补安装时只需按照安装的过程进行，只是在选择组件时只勾选要补装的组件或工具箱即可. 矚慫润厲钐瘗睞枥庑赖。
4、MATLAB 操作桌面有几个窗口？如何使某个窗口脱离桌面成为独立窗口？又如何将脱离出去的窗口重新放置到桌面上？聞創沟燴鐺險爱氇谴净。
答：在MATLAB 操作桌面上有五个窗口，在每个窗口的右下角有两个小按钮，一个是关闭窗口的Close 按钮，一个是可以使窗口称为独立的Undock 按钮，点击Undock 按钮就可以使该窗口脱离桌面称为独立窗口，在独立窗口的view 菜单中选择Dock，菜单项就可以将独立的窗口重新防止的桌面上.残骛楼諍锩瀨濟溆塹籟。
5、如何启动M 文件编辑/调试器？答：在操作桌面上选择“建立新文件”或“打开文件”操作时，M 文件编辑/调试器将被启动.在命令窗口中键入edit 命令时也可以启动M 文件编辑/调试器.酽锕极額閉镇桧猪訣锥。
6、存储在工作空间中的数组能编辑吗？如何操作？答：存储在工作空间的数组可以通过数组编辑器进行编辑：在工作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输入修改内容即可.彈贸摄尔霁毙攬砖卤庑。
7、命令历史窗口除了可以观察前面键入的命令外，还有什么用途？答：命令历史窗口除了用于查询以前键入的命令外，还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中. 謀荞抟箧飆鐸怼类蒋薔。
第5章_MATLAB绘图
参考教材 : MATLAB程序设计与应用(第二版) 刘卫国 主编 高等教育出版社
计算机学院 郭江鸿 21B434 82519604
5.1 二维图形 5.1.1 绘制二维曲线 1. plot函数的基本调用格式为： plot(x , y) 其中x和y为长度相同的向量，分别用于 存储二维曲线的x坐标和y坐标数据。
90 0.5 120 0.4 60
0.3 150 0.2 30
0.1
180
0
210
330
240 270
300
3．对数坐标图形 MATLAB提供了绘制对数和半对数坐标曲线 的函数，调用格式为： semilogx(x1,y1,选项1,x2,y2,选项2,…) semilogy(x1,y1,选项1,x2,y2,选项2,…) loglog(x1,y1,选项1,x2,y2,选项2,…) 例5.10 绘制y=10x2的对数坐标图并与直角线性 坐标图进行比较。
Line in 3-D Space
5
Z
0
origin
-5 -5 -5 0 5 5 X Y 0
0.8 0.6 0.4
t=0:0.15:2*pi; x=sin(t); y=cos(t); z=cos(2*t); plot3(x,y,z,'b-',x,y,z, 'rd')
axis([0 10 0 2.5]) %设置坐标轴刻度范围
例5.5 绘制分段函数曲线并添加图形标注。
x=linspace(0,10,100); y=[ ]; for x0=x if x0>=8 y = [y,1]; elseif x0>=6 y = [y,5-x0/2]; elseif x0>=4 y = [y , 2]; elseif x0>=0 y = [y,sqrt(x0)]; end end plot(x,y) axis([0 10 0 2.5]) title('分段函数曲线'); xlabel('Variable X'); ylabel('Variable Y');
自动控制原理(第2版)(余成波_张莲_胡晓倩)习题全解及MATLAB实验第5章习题解答(DOC)
频率特性法频域分析法是一种图解分析法，可以根据系统的开环频率特性去判断闭环系统的性能，并能较方便地分 析系统参量对系统性能的影响，从而指出改善系统性能的途径，已经发展成为一种实用的工程方法，其主要 内容是:1）频率特性是线性定常系统在正弦函数作用下，稳态输出与输入的复数之比对频率的函数关系。
频率特性是传递函数的一种特殊形式，也是频域中的数学模型。
频率特性既可以根据系统的工作原理，应用机理 分析法建立起来，也可以由系统的其它数学模型（传递函数、微分方程等）转换得到，或用实验法来确定。
2）在工程分析和设计中，通常把频率特性画成一些曲线。
频率特性图形因其采用的坐标不同而分为幅相特性（Nyquist 图）、对数频率特性（Bode 图）和对数幅相特性（Nichols 图）等形式。
各种形式之间是互通的， 每 种形式有其特定的适用场合。
开环幅相特性在分析闭环系统的稳定性时比较直观，理论分析时经常采用；波 德图可用渐近线近似地绘制，计算简单，绘图容易，在分析典型环节参数变化对系统性能的影响时最方便; 由开环频率特性获取闭环频率指标时，则用对数幅相特性最直接。
3）开环对数频率特性曲线 （波德图）是控制系统分析和设计的主要工具。
开环对数幅频特性的斜率表征了系统的型别 （V ）,其高度则表征了开环传递系数的大小，因而低频段表征系统稳态性能; 中频段的斜率、宽度以及幅值穿越频率，表征着系统的动态性能；高频段则表征了系统抗高频干扰的能力。
对于最小相位系统，幅频特性和相频特性之间存在着唯一的对应关系，根据对数幅频特性，可以唯一地 确定相应的相频特性和传递函数。
平面中的（一I , j0）点来判断闭环系统的稳定性。
利用奈奎斯特稳定判据，可根据系统的开环频率特性来判断 闭环系统的稳定性，并可定量地反映系统的相对稳定性，即稳定裕度。
稳定裕度通常用相角裕量和幅值裕量 来表示。
5）利用开环频率特性或闭环频率特性的某些特征量，均可对系统的时域性能指标作出间接的评估。
第5章 MATLAB绘图_习题答案
第5章MＡTLAB绘图习题5一、选择题1.如果x、ｙ均为４×3矩阵，则执行plot（ｘ,ｙ）命令后在图形窗口中绘制( ）条曲线。
DA。
12 B。
７ C.4 Ｄ.32.下列程序得运行结果就是（）。
Aｘ=0:pi/１00:２＊ｐi;for n=1:2：１０ｐlot（n＊ｓiｎ（x），n＊cos（x)）hold oneｎｄaｘｉｓsｑuarｅA。
5个同心圆Ｂ.5根平行线Ｃ．一根正弦曲线与一根余弦曲线 D.5根正弦曲线与5根余弦曲线3．命令ｔｅxt（1,1，'{\ａｌpｈa｝+｛\ｂeta｝’）执行后,得到得标注效果就是( )。
CＡ.{\ａlｐha}+｛\bｅta｝B。
{\α}＋{\β｝ C.α+βＤ.＼α＋\β4。
subplot(２，2，3）就是指( )得子图。
AA.两行两列得左下图Ｂ.两行两列得右下图Ｃ.两行两列得左上图Ｄ．两行两列得右上图5。
要使函数ｙ=2ｅx得曲线绘制成直线，应采用得绘图函数就是(）.ＣA.pｏｌａｒB。
ｓeｍiloｇｘC。
semiｌogy D。
loｇlｏg６.下列程序得运行结果就是( )。
B[x，ｙ］＝meshｇｒｉd(1：５);surf（x,ｙ,5＊ｏnes(sｉze（x)）)；A.z=x+y平面B.与xy平面平行得平面Ｃ。
与xy平面垂直得平面 D.z=5x平面７.下列函数中不能用于隐函数绘图得就是（)。
ＤA.ezｍeｓhB.eｚsurｆC.ｅzｐloｔD．pｌoｔ38.下列程序运行后，瞧到得图形（）.Ct=0：pi/20：2*pi；［ｘ，ｙ］＝ｍeｓhgrid（-8：0、5:8）;ｚ=siｎ（sqrt（ｘ、＾2+y、＾2）)、／sｑrt（x、^2＋ｙ、^2+eps)；surf（ｘ，y,z)viｅw（0,9０)；axiｓ eｑualA．像墨西哥帽子 B.就是空心得圆C。
边界就是正方形D．就是实心得圆９。
下列程序运行后得到得图形就是( ）.Ａ[x，y］=meshｇrid（－2：2）；z=ｘ＋y;i=ｆind（aｂs（ｘ）<1 ＆ abs（y）<1);z（i）=ＮaＮ;ｓurf（x,ｙ，z)；ｓhａdinｇｉｎterpA。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
》
第5章 MATLAB绘图
习题5
一、选择题
1．如果x、y均为4×3矩阵，则执行plot(x,y)命令后在图形窗口中绘制（）条曲线。
D
A．12 B．7 C．4 D．3
2．下列程序的运行结果是（）。
A
x=0:pi/100:2*pi;
for n=1:2:10
…
plot(n*sin(x),n*cos(x))
hold on
end
axis square
A．5个同心圆 B．5根平行线
C．一根正弦曲线和一根余弦曲线 D．5根正弦曲线和5根余弦曲线3．命令text(1,1,'{\alpha}+{\beta}')执行后，得到的标注效果是（）。
C A．{\alpha}+{\beta} B．{\α}+{\β} C．α+β D．\α+\β
；
4．subplot(2,2,3)是指（）的子图。
A
A．两行两列的左下图 B．两行两列的右下图
C．两行两列的左上图 D．两行两列的右上图
5．要使函数y=2e x的曲线绘制成直线，应采用的绘图函数是（）。
C A．polar B．semilogx C．semilogy D．loglog 6．下列程序的运行结果是（）。
B
[x,y]=meshgrid(1:5);
surf(x,y,5*ones(size(x)));
"
A．z=x+y平面 B．与xy平面平行的平面
C．与xy平面垂直的平面 D．z=5x平面
7．下列函数中不能用于隐函数绘图的是（）。
D
A．ezmesh B．ezsurf C．ezplot D．plot3 8．下列程序运行后，看到的图形（）。
C
t=0:pi/20:2*pi;
[x,y]=meshgrid(-8::8);
z=sin(sqrt(x.^2+y.^2))./sqrt(x.^2+y.^2+eps);
…
surf(x,y,z)
view(0,90);axis equal
A．像墨西哥帽子 B．是空心的圆
C．边界是正方形 D．是实心的圆
9．下列程序运行后得到的图形是（）。
A
[x,y]=meshgrid(-2:2);
z=x+y;
i=find(abs(x)<1 & abs(y)<1);
、
z(i)=NaN;
surf(x,y,z);shading interp
A．在一个正方形的正中心挖掉了一个小的正方形
B．在一个正方形的正中心挖掉了一个小的长方形
C．在一个正方形的上端挖掉了一个小的正方形
D．在一个正方形的下端挖掉了一个小的正方形
10．在使用MATLAB“绘图”选项卡中的命令按钮绘图之前，需要（）。
B
A．在命令行窗口中输入绘图命令 B．在工作区中选择绘图变量
）
C．打开绘图窗口 D．建立M文件
二、填空题
1．执行以下命令：
x=0:pi/20:pi;
y=sin(x);
以x为横坐标、y为纵坐标的曲线图绘制命令为，给该图形加上“正弦波”标题的命令为，给该图形的横坐标标注为“时间”，纵坐标标注为“幅度”的命令分别为和。
plot(x,y)，title(‘正弦波’)，xlabel(‘时间’)，ylabel(‘幅度’)
2．在同一图形窗口中绘制y1和y2两条曲线，并对图形坐标轴进行控制，请补充程序。
x=-3::3;
%
y1=2*x+5;
y2=x.^2-3*x;
plot(x,y1) %绘制曲线y1
① ;
plot(x,y2) %绘制曲线y2
m1=max([y1,y2]);
m2=min([y1,y2]);
② ; %用axis 调制坐标轴，横坐标在[-3,3]之间，纵坐标在[-5,20]之间
#
①hold on ②axis([-3,3,-5,20])
3．下列命令执行后得到的图形是 。
（单位）圆 x=@(t) sin(t);
y=@(t) cos(t);
ezplot(x,y)
4．某工厂2015年度各季度产值分别为、、、，为了表示各季度产值占全年总产值的比例，可以绘制 ，其命令是 。
饼图，pie([,,,]) 或pie3([,,,])
三、应用题
)
1．绘制下列曲线。
（1）2
1100x y += （2）2221x e y -=π （3）122
=+y x （4）⎩⎨⎧==32
5t y t x (1)
x=-1:1/180:1;
y=100./(1+x.^2);
plot(x,y)
(2)
…
x=-1:1/180:1; y=1/(2*pi).*exp(-x.^2./2);
plot(x,y)
(3)
ezplot('x^2+y^2-1=0')
(4)
ezplot('t^2','5*t^3',[,])
2．分别用plot 和fplot 函数绘制函数x y 1sin
=的曲线，分析两曲线的差别。
@
3．绘制下列极坐标图。
（1）4cos 5+=θρ （2）θρ12
=
（3）7cos 5-=
θρ （4）23
θπρ= (1) theta=0::2*pi;
rho=5.*cos(theta)+4;
polar(theta,rho)
(2)
%
theta=::2*pi; rho=12./(sqrt(theta));
polar(theta,rho,'r')
(3)
theta=0:pi/90:2*pi;
t=cos(theta);
a=find(t<;
t(a)=;
】
rho=5./t-7; polar(t,rho)
(4)
theta=0:pi/90:2*pi;
rho=pi/3.*(theta.^2);
polar(theta,rho)
4．在同一坐标轴中绘制下列两条曲线并标注两曲线交叉点。
（1）y=
)
（2）π≤≤⎩⎨⎧==t 0，sin )3sin(cos )3sin(t t y t
t x
t=0::2*pi;
x=sin(3*t).*cos(t);
y=sin(3*t).*sin(t);
plot(x,y);
hold on;
x=-1::1;
y=2*;
—
plot(x,y); hold off
5．绘制下列三维图形。
（1）⎪⎩⎪⎨⎧===t z t y t x sin cos （2）⎪⎩
⎪⎨⎧=+=+=u z v u y v u x sin sin )cos 1(cos )cos 1( （3）5=z （4）半径为10的球面
(1)
t=0:pi/180:2*pi;
x=cos(t);
y=sin(t);
z=t;
plot3(x,y,z)
(2)
u=0:pi/50:2*pi;
v=0:pi/50:2*pi;
ezsurf('(1+cos(u))*cos(v)','(1+cos(u))*sin(v)','sin(u)')
(4) 为什么用不了ezplot3('x^2+y^2+z^2-100')
u=0:pi/50:2*pi;
v=0:pi/50:2*pi;
ezsurf('cos(u)*cos(v)','cos(u)*sin(v)','sin(u)')。