完整word版MATLAB计算结果可视化Word文档格式.docx

1、t=（0:pi/100:pi）; %长度为101的时间采样列向量 y1=sin（t）*1,-1; %包络线函数值，是（101x2）的矩阵 y2=sin（t）.*sin（9*t）; %长度为101的调制波列向量 t3=pi*（0:9）/9; y3=sin（t3）.*sin（9*t3）;plot（t,y1,r:,t,y2,b,t3,y3,bo） axis（0,pi,-1,1） %控制轴的范围 图6-3【例6-4】fplot 与一般绘图指令的绘图效果比较。x,y=fplot（cos（tan（pi*x）,-0.4,1.4,0.2e-3）;n=length（x）;subplot（1,2,1）,plot（

2、x,y）title（fontsize20fontname隶书泛函绘图指令效果t=（-0.4:1.8/n:1.4）subplot（1,2,2）,plot（t,cos（tan（pi*t）fontsize20fontname隶书等分采样绘图 图6-4 右图曲线中底部有明显绘图缺陷【例6-5】用图形演示平面上一个方块四个顶点在仿射投影（Affine Projection）下的位置、形状变化。%平面上的四个点和它们构成的方块p1=-0.5,0,1p2=-0.5,1,1p3=0.5,1,1p4=0.5,0,1Sq=p1,p2,p3,p4,p1;%平移投影：沿x轴移动0.5 , 沿y轴移动1 。dx=0.5

3、;dy=1;T=1,0,dx;0,1,dy;0,0,1;%旋转投影：逆时针旋转30度。th=pi/6;R=cos（th）,-sin（th）,0;sin（th）,cos（th）,0;%刻度投影：x方向放大到2倍，y方向放大到3倍。alpha=2;beta=3;S=alpha,0,0;0,beta,0;E=eye（3,3）;%为编程方便，设计一个单位阵。在它投影下，任何形状都不变。TRS=E,T,R,S; %用元胞数组存放四个变换矩阵 ss=rrdrprh; %用元胞数组存放数据点形设置值 tt=Original SquareTranslationRotationScaling %用元胞数组存放四

4、张子图的图名 for i=1:4 W=TRSi*Sq; %进行仿射变换 subplot（2,2,i） for k=1: plot（W（1,k）,W（2,k）,ssk）; %用不同点形标志四个顶点 axis（-3,3,-1,5）,axis equal hold on %使以后图形画在当前子图上 end plot（W（1,:）,W（2,:） %连接四个顶点使图形封闭 grid on %画过坐标格线 title（tti） %给子图题写图名 hold off %使以后图形不再画在当前子图上end 图 6-5 仿射投影演示【例6-6】观察各种轴控制指令的影响。演示采用长轴为3.25，短轴为1.15的椭圆

5、。注意：采用多子图表现时，图形形状不仅受“控制指令”影响，而且受整个图面“宽高比”及“子图数目”的影响。t=0:2*pi/99:2*pi;x=1.15*cos（t）;y=3.25*sin（t）; %y为长轴，x为短轴subplot（2,3,1）,plot（x,y）,axis normal,grid on,Normal and Grid onsubplot（2,3,2）,plot（x,y）,axis equal,grid on,title（Equalsubplot（2,3,3）,plot（x,y）,axis square,grid on,title（Squaresubplot（2,3,4）,pl

6、ot（x,y）,axis image,box off,title（Image and Box offsubplot（2,3,5）,plot（x,y）,axis image fill,box offImage and Fillsubplot（2,3,6）,plot（x,y）,axis tight,box off,title（Tight图6-6 各种轴控制指令的不同影响【例6-7】绘制二阶系统阶跃响应。clf;t=6*pi*（0:100）/100;y=1-exp（-0.3*t）.*cos（0.7*t）;tt=t（find（abs（y-1）0.05）;ts=max（tt）; subplot（1,2,

7、1）,plot（t,y,r-LineWidth,3）,grid on axis（0,6*pi,0.6,max（y） y=1exp（-alpha*t）*cos（omega*t） text（11,1.25,alpha=0.3）;text（11,1.15,omega=0.7hold on;plot（ts,0.95,10）;hold off text（ts+1.5,0.95,ts= num2str（ts）xlabel（t -）,ylabel（y -） subplot（1,2,2）,plot（t,y,3） axis（-inf,6*pi,0.6,inf） set（gca,Xtick,2*pi,4*pi,6

8、*pi,Ytick,0.95,1,1.05,max（y） grid on it y = 1 - e -alphatcosomegat） text（13.5,1.2,fontsize12alpha=0.3） text（13.5,1.1,fontsize12omega=0.7cell_string1=fontsize12uparrowcell_string2=fontsize16 fontname隶书镇定时间 cell_string3=fontsize6 cell_string4=fontsize14rmt_s = num2str（ts）;text（ts,0.85,cell_string） fon

9、tsize14 bft rightarrowylabel（fontsize14 bfy rightarrow图 6-7 二阶阶跃响应图用MATLAB4.x和5.x版标识时的差别【例6-8】利用hold绘制离散信号通过零阶保持器后产生的波形。t=2*pi*（0:20）/20;y=cos（t）.*exp（-0.4*t）;stem（t,y,gstairs（t,y,rhold off 图6-8 离散信号的重构【例6-9】受热压力容器的期望温度是120度，期望压力是0.25Mpa。在同一张图上画出它们的阶跃响应曲线。S1=tf（1 1,1 3 2 1）; %温度的传递函数对象模型。S2=tf（1,1 1

10、 1）; %压力的传递函数对象模型。Y1,T1=step（S1）; %计算阶跃响应Y2,T2=step（S2）;plotyy（T1,120*Y1,T2,0.25*Y2,stairsplot图6-9 双纵坐标图演示【例6-10】演示subplot指令对图形窗的分割。t=（pi*（0:1000）/1000）y1=sin（t）;y2=sin（10*t）;y12=sin（t）.*sin（10*t）;subplot（2,2,1）,plot（t,y1）;axis（0,pi,-1,1）subplot（2,2,2）,plot（t,y2）;subplot（position,0.2,0.05,0.6,0.45）

11、plot（t,y12,b-,t,y1,-y1,axis（0,pi,-1,1） 图6-10 多子图的布置【例6-11】简单例题。0.02:2）*pi;x=sin（t）;y=cos（t）;z=cos（2*t）;plot3（x,y,z,x,y,z,bd）,view（-82,58）,box on,legend（链宝石图 6-11 宝石项链【例6-12】用曲面图表现函数clf,x=-4:4;y=x;X,Y=meshgrid（x,y）; %生成 x-y 坐标“格点”矩阵Z=X.2+Y.2; %计算格点上的函数值surf（X,Y,Z）;hold on,colormap（hot）stem3（X,Y,Z,） %

12、用来表现在格点上计算函数值图 6-12 曲面图和格点, barh, bar3, bar3h【例6-13】二维直方图有两种图型：垂直直方图和水平直方图。而每种图型又有两种表现模式：累计式：分组式。本例选其两种加以表现。x=-2:2; %注意：自变量要单调变化Y=3,5,2,4,1;3,4,5,2,1;5,4,3,2,5; %各因素的相对贡献份额subplot（1,2,1）,bar（x,Ystacked） %“累计式”直方图xSigma y）,colormap（cool）%控制直方图的用色legend（因素A因素B因素Csubplot（1,2,2）,barh（xgrouped） %“分组式”水平直

13、方图y图 6-13 二维直方图, pie3【例6-14】饼图指令pie , pie3 用来表示各元素占总和的百分数。该指令第二输入宗量为与第一宗量同长的0-1向量，1使对应扇块突出。a=1,1.6,1.2,0.8,2.1;subplot（1,2,1）,pie（a,1 0 1 0 0）,legend（12345）subplot（1,2,2）,pie3（a,a=min（a）,colormap（cool） 图 6-14 饼形统计图6.3图形窗的图形编辑功能【例6-15】运行以下指令shg,t=（pi*（0:）,axis（0,pi,-1,1） 图6-15 MATLAB 5.3版图形窗【例6-16】编辑

14、图形图 6-16 经交互编辑后的图形6.4函数绘图的简捷指令【例6-17】绘制和它的积分在间的图形。syms t tao;y=2/3*exp（-t/2）*cos（sqrt（3）/2*t）; %定义符号函数s=subs（int（y,t,0,tao）,tao,t）; %获得积分函数subplot（1,2,1）,ezplot（y,0,4*pi）;grid %一元函数简捷绘图指令subplot（1,2,2）,ezplot（s,0,4*pi）;grids = inty（t）dt） %重写右子图图名图 6-17 ezplot使用示例【例6-18】画出所表示的三维曲面。的取值范围是解：clear; %清内存

15、x=-8:0.5:8; %定义自变量x的一维刻度向量y=x %定义自变量y的一维刻度向量，“”为矩阵转置X=ones（size（y）*x; %计算自变量平面上取值点坐标的二维数组Y=y*ones（size（x）;R=sqrt（X.2+Y.2）+eps; %计算中间变量 Z=sin（R）./R; %计算与自变量二维数组相应的函数值mesh（Z）; %绘制三维网格图colormap（hot） %指定网格图用hot色图绘制 图 6-18 三维网线图 常用绘图函数plot 绘制二维图形viewmtx定义观察点的矩阵Plot3绘制三维图形contour等高线图fplot绘图函数Contour3三维等高线

16、图subplot将当前图分为若干子图meshgrid产生网格semilogx,semilogy对数绘图Surf曲面图view观察点的位置或角度Mesh网格曲面图waterfallFigure设为当前图或创建图surfl,surfc，surfnorm带有特殊光照，等高线和法线的曲面图surface曲面对象fill,fill3填充多边形Line线对象colormap读或设置颜色表Text文本对象newplot设定下一个图的属性whitebg背景颜色slice三个变量的绘图函数Clf清除图形hold保留当前图形Axis坐标轴刻度显示或隐藏网格线comet,comet3动态绘制二维、三维图形p o l

17、 a r在极坐标中绘图semilogx,semil ogy,loglog习题1绘制函数时的曲线。2. 绘制函数y=sin（t） 的曲线。3.编写程序作图：函数y=*sin（3*x）+cos（x），区间：0x4, -5y5。（参考：fplot）4.利用peaks函数产生矩阵x、y、z，在一个图形窗中作出曲面图和等值线图（曲面图位于界面上部，等值线图位于界面下部）。peaks,subplot,surf,contour）仅供个人用于学习、研究；不得用于商业用途。For personal use only in study and research; not for commercial use.Nur fr den persnlichen fr Studien, Forschung, zu kommerziellen Zwecken verwendet werden.Pour l tude et la recherche uniquement des fins personnelles; pas des fins commerciales. , , . 以下无正文