MATLAB_图形绘制.pptVIP

哈尔滨理工大学机械动力工程学院 蔺勇智制作 MATLAB　图形绘制 蔺勇智制作 2008.4.9 本次课程主要介绍内容 基本图形绘制(平面) 三维图像绘制 特殊图形绘制 MATLAB编程绘制曲线 仿真系统分析图形绘制 基本图形绘制 三维图形绘制 准备绘图数据。数据可以是实际工作中采集的数据，也可以由各种命令或函数创建。如峰形函数peaks来生成 实验数据。 创建出图形窗口。可以用figure函数创建，也可以由灰 图函数自动创建窗口。 调用二维绘图命令。例如plot、mesh、surf等函数。当然，也可以自己编写灰图函数命令，根据需要绘制图 形。 修饰图形。例如标注轴、设置视角、设置光源、改变着着色模式等，当然，也可以使用MATLAB提供的默认 值。 特殊图形绘制 MATLAB编程绘制曲线 仿真系统分析图形绘制 MATLAB也可以运用循环语句来绘制曲线，本节主要介绍如果根据已知数学表达式来运用MATLAB循环语句绘制曲线的方法，现在向大家介绍一个用switch-case语句绘制 曲线的例子。 例：绘制 的曲线。 程序如下： t=[0:0.1:18]; for x=0.2:0.2:0.8; 返回 b=sqrt([1-x^2]);z=atan(b/x); y1=-t*x;y2=t*b+z;y=1-exp(y1).*sin(y2)/b; switch round(10*x) %将x圆整成整数 case 2 plot(t,y,‘r-’),hold on %采用switch结构，根据不同x值， 曲线用不同颜色、数据点画出。 case 4 plot(t,y,b*),hold on case 6 plot(t,y,k.),hold on otherwise plot(t,y,g+),hold on end xlabel(t(秒)),ylabel(y),title(二阶系统阶跃响应) legend({\xi}=0.2,{\xi}=0.4,{\xi}=0.6,{\xi}=0.8) 在介绍此类曲线的绘制之前，先向大家 介绍一些这部分的基本知识。 1.传递函数模型： 对于一个连续单输入单输出 (SISO) 的LTI系统， 设输入量为Xi(s)，输出量为Xo(s)则系统的传递 函数G(s)可以表示为 返回 对应的脉冲传递函数为 在MATLAB中，传递函数描述法是通过传递函数分子和分母关于s降序排列的多项式系数来表示的，并用向量num和den表示： num = [b0 b1 … bm－1 bm] den = [a0 a1 … am－1 am] 在本节中，主要向大家介绍一下连续系统 频率特性的MATLAB函数。所以这里有必要 给大家介绍一下如何运用MATLAB对频率响应 函数的计算方法。 2.频率响应函数的计算 (1)polyval()因为G(jw)的分子、分母均为有理 多项式，可用多项式计算指令ployval()计算系 统的频率响应，其调用格式为 Y=polyval(P,X) 【说明】 P是多项式系数向量(降幂排列)，X是自变量(应设为jw)，Y是返回的计算结果 (复数数组)。 对Y利用abs()、angle()即可求出系统的 幅频特性和相频特性。 例如：绘制系统 的幅频 特性曲线和相频特性曲线。 程序如下： num=[11 11];den=[1 15 4 0]; w=0.05:0.01:0.5*pi;%产生频率向量 Gw=polyval(num,j*w)./polyval(den,j*w);%产生频率特性 mag=abs(Gw);%计算幅频特性 theta=angle(Gw);%计算相频特性 subplot(2,1,1),plot(w,mag) grid,title(幅频特性) ylabel(|G|) subplot(2,1,2),plot(w,theta) grid,title(相频特性) xlabel(\omega(rad/s)),ylabel(deg) (2) freqs()若已知系统传递函数，还可以用 freqs求系统的频率相应。其有以下几种调 用格式 h=freqs(b,a,w) 指定实式角频率,返回 响应值。 [h,w]=freqs(b,a)自动确定200个频率 点，返回响应值和