matlaba第五章.docx
- 文档编号:29682954
- 上传时间:2023-07-26
- 格式:DOCX
- 页数:11
- 大小:89.49KB
matlaba第五章.docx
《matlaba第五章.docx》由会员分享,可在线阅读,更多相关《matlaba第五章.docx(11页珍藏版)》请在冰豆网上搜索。
matlaba第五章
matlab2012a第五章
x,acost,a,4,b,21,已知椭圆的长、短轴,用“小红点线”画如图p5-1所示的椭圆。
y,bsint,
,提示:
参量,点的大小,axisequal,t3
2
1
0
y
-1
-2
-3
-4-3-2-101234
x图p5-1clear
t=2*pi*(0:
149)/149;
a=4;b=2;
x=a*cos(t);y=b*sin(t);plot(x,y,'r.','MarkerSize',15)axisequal
3
2
1
0
-1
-2
-3
-3-2-101234
2,根据表达式,,1,cos,绘制如图p5-2的心脏线。
,提示:
polar,注意title中特殊字符,
线宽,axissquare。
可以用plot试试。
,
=1-cos,
902120601.5
115030
0.5
1800
210330
240300
270
图p5-2
clear
t=0:
.01:
2*pi;
P=1-cos(t);
f=polar(t,P,'r');
axissquare;
set(f,'LineWidth',4)
title('P=1-cos\theta')
P=1-cos
90212060
1.5
115030
0.5
1800
210330
240300
270
3,A,B,C三个城市上半年每个月的国民生产总值表p5.1。
试画出如图p5-3所示的三城市上半
年每月生产总值的累计直方图。
,提示:
bar(x,Y,'style');colormap(cool);legend。
,
表p5.1各城市生产总值数据,单位:
亿元,
城市1月2月3月4月5月6月
A170120180200190220
B120100110180170180600C70508010095120
A
B500C
400
300
200
100
0123456
图p5-3clear
X=(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);
legend('A','B','C',2)
600
A
B
500C
400
300
200
100
0123456
4,二阶线性系统的归一化,即令,冲激响应可表示为:
,1n
1,,,tetsin(,),0,,,1,,,t2ytte(),,,1,其中,,|1,,|,为阻尼系数,,1,(,,,)t,(,,,)t,,,1,,ee,,2,,
。
,1,希望在同一张图上,绘制区间内不同取值时的各条t,[0,18],,0.2:
0.2:
1.4
曲线,参见图p5-4,。
在此图上,的各条曲线为细蓝线,为粗黑线,为细,,1,,1,,1
红线,并且对最上方及最下方的两条曲线给出和的醒目标志。
,2,读者运,,0.2,,1.4
行题下程序exmp504.m,可以发现该程序画出的曲线中没有“粗黑线”。
你能讲出原因吗,0.8
如何对exmp504.m作最少的修改,比如只改一条指令,,就可画出所需图形。
,提示:
该题深,,=0.2层次地暴露数值计算可能存在的隐患。
,,=0.20.6
0.4
0.2,=1.4
0
-0.2
-0.4024681012141618
图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);
holdon
fork=1:
L
zk=zeta(k);
beta=sqrt(abs(1-zk^2));
ifzk-1<-3*eps%满足此条件,绘蓝色线修改此处
y=1/beta*exp(-zk*t).*sin(beta*t);
plot(t,y,'b')
ifzk<0.4
text(2.2,0.63,'\zeta=0.2')
end
elseifzk-1 y=t.*exp(-t); plot(t,y,'k','LineWidth',2) else%其余,绘红色线 y=(exp(-(zk-beta)*t)-exp(-(zk+beta)*t))/(2*beta); plot(t,y,'r') ifzk>1.2 text(0.3,0.14,'\zeta=1.4') end end end text(10,0.7,'\Delta\zeta=0.2')axis([0,18,-0.4,0.8])holdoff boxon gridon 0.8 ,=0.2 =0.20.6 0.4 0.2 =1.4 0 -0.2 -0.4024681012141618 5,用绿实线绘制,,的三维曲线,曲线如图p5-5所示。 ,提示: 参x,sin(t)y,cos(t)z,t 15变量,plot3,线色线粗。 , 10 5 01 10.500-0.5-1-1 图p5-5clear t=(0: 0.02: 4)*pi; x=sin(t);y=cos(t);z=t;plot3(x,y,z,'g-','LineWidth',3)boxon 15 10 5 01 0.51 0.500-0.5-0.5-1-1 22,x,y6,采用两种不同方法绘制在的如图p5-6的三维,透视,网格曲面。 z,4xex,y,[,3,3] ,提示: 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]) hiddenoff 方法二: clear symsxyz; z=4*x*exp(-x^2-y^2);ezmesh(z,[-3,3]) colormap(jet) axis([-3,3,-3,3,-2,2])hiddenoff sin(x,y)7,在区间里,根据表达式,绘制如图p5-7所示的曲面。 z,x,y,[,4,,4,]x,y ,提示: 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) shadinginterp title('z=sin(x+y)/(x+y)') 8,试用图解法回答: y,,0.1,22,1,方程组有多少个实数解,,提示: 图解法;ezplot;ginput,(1,x,y), sin(x,cos(y)),0, ,2,求出离最近、且满足该方程组的一个近似解。 x,0,y,0 clear ezplot('y/(1+x^2+y^2)-0.1',[-2*pi,2*pi,-pi/4,4*pi])holdon ezplot('sin(x+cos(y))',[-2*pi,2*pi,-pi/4,4*pi])gridon sin(x+cos(y))=012 10 8 6y 4 2 0 -6-4-20246 x 由图可知,共有6个交点,故方程组共有六组实数解求离最近、且满足该方程组的一个近似解x,0,y,0 [x,y]=ginput (1) x= -0.9797 y= 0.2002 sin(x+cos(y))=00.2006 0.2005 0.2004 0.2003 y0.2002 0.2001 0.2 0.1999 0.1998-0.98-0.9799-0.9798-0.9797-0.9796-0.9795-0.9794 x 9,制作如文件prob509.p,在光盘的mfiles文件夹上,运行时那样的色图变幻,参见图p5-8,。 ,提示: [jet;flipud(jet)];colormap;spinmap, 图p5-8 clear [X,Y,Z]=sphere(30); colormap(jet) surf(X,Y,Z) axisoff,axisequal,shadinginterplight('position',[0-101.5],'style','infinite') lightingphong,materialshinylight; lightingflat,set(gcf,'Color','w')view([-160,30]),shg,C=jet;CC=[C;flipud(C)]; colormap(CC) disp('按任意键,观察图像变换') pause,spinmap(40,8) 观察图像变换按任意键
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- matlaba 第五