实验五MATLAB绘制Bode图(精)

实验五

MATLAB绘制Bode图、Nyquist图实验学科：自动控制理论

姓名：XXX

学号：XXX

班级：XXX

实验五 MATLAB绘制Bode图、Nyquist图实验

一、实验目的

1. 掌握计算机绘制系统的Bode图、Nyquist图、根轨迹曲线方法；

2. 观察系统对单位阶跃信号的时域响应曲线；

3. 通过仿真结果和理论计算的对照，加深对系统特性的理解。

二、实验设备

PC、MATLAB

三、实验结果

1. 各传递函数的Bode图

(1 K (K = 1, K = 0.5

GH(s = 1

GH(s = 0.5

(2 (T = 0.2, T = 0.4

(3 (K = 0.5, T = 0.4; K = 1.0; T = 0.2

(4 (K = 1.0, T = 0.4; K = 1.0; T = 0.2

2. 各传递函数的Nyquist图

(1 K (K = 1, K = 0.5

GH(s = 1

GH(s = 0.5

(2 (T = 0.2, T = 0.4

(3 (K = 0.5, T = 0.4; K = 1.0; T = 0.2

(4 (K = 1.0, T = 0.4; K = 1.0; T = 0.2

3. 各传递函数的根轨迹曲线

(1 K (K = 1, K = 0.5

GH(s = 1, GH(s = 0.5

(2 (T = 0.2, T = 0.4

(5 (K = 0.5, T = 0.4; K = 1.0; T = 0.2

(6 (K = 1.0, T = 0.4; K = 1.0; T = 0.2

4. 各传递函数对单位阶跃信号的时域响应曲线

(1 K (K = 1, K = 0.5

GH(s = 1

GH(s = 0.5

(2 (T = 0.2, T = 0.4

(7 (K = 0.5, T = 0.4; K = 1.0; T = 0.2

(8 (K = 1.0, T = 0.4; K = 1.0; T = 0.2

5. 各传递函数的特征量如下表所示

传递函数GH(S

剪切频率

(rad/s 相角余量

(°

特征方程的根

GH(s = 1

GH(s = 0.5

5 90 -5

2.5 90 -2.5

2.89 135 -

∞135 -2.5

2.5 135 -5 5 135 -10

Matlab中Bode图的绘制技巧(精)

Matlab中Bode图的绘制技巧 我们经常会遇到使用Matlab画伯德图的情况，可能我们我们都知道bode这个函数是用来画bode图的，这个函数是Matlab内部提供的一个函数，我们可以很方便的用它来画伯德图，但是对于初学者来说，可能用起来就没有那么方便了。 譬如我们要画出下面这个传递函数的伯德图： 1.576e010 s^2 H(s= ------------------------------------------------------------------------------------------ s^4 + 1.775e005 s^3 + 1.579e010 s^2 + 2.804e012 s + 2.494e014 (这是一个用butter函数产生的2阶的，频率范围为[20 20K]HZ的带通滤波器。 我们可以用下面的语句： num=[1.576e010 0 0]; den=[1 1.775e005 1.579e010 2.804e012 2.494e014]; H=tf(num,den; bode(H 这样，我们就可以得到以下的伯德图： 可能我们会对这个图很不满意，第一，它的横坐标是rad/s，而我们一般希望横坐标是HZ；第二，横坐标的范围让我们看起来很不爽；第三，网格没有打开（这点当然我们可以通过在后面加上grid on解决）。 下面，我们来看看如何定制我们自己的伯德图风格： 在命令窗口中输入：bodeoptions

我们可以看到以下内容：ans = Title: [1x1 struct] XLabel: [1x1 struct] YLabel: [1x1 struct] TickLabel: [1x1 struct] Grid: 'off' XLim: {[1 10]} XLimMode: {'auto'} YLim: {[1 10]} YLimMode: {'auto'} IOGrouping: 'none' InputLabels: [1x1 struct] OutputLabels: [1x1 struct] InputVisible: {'on'} OutputVisible: {'on'} FreqUnits: 'rad/sec' FreqScale: 'log' MagUnits: 'dB' MagScale: 'linear' MagVisible: 'on' MagLowerLimMode: 'auto' MagLowerLim: 0 PhaseUnits: 'deg' PhaseVisible: 'on' PhaseWrapping: 'off'

MATLAB中bode图绘制技巧(精)

Matlab中Bode图的绘制技巧学术收藏2010-06-04 21:21:48 阅读54 评论0 字号：大中小订阅我们经常会遇到使用Matlab画伯德图的情况，可能我们我们都知道bode这个函数是用来画bode图的，这个函数是Matlab内部提供的一个函数，我们可以很方便的用它来画伯德图，但是对于初学者来说，可能用起来就没有那么方便了。譬如我们要画出下面这个传递函数的伯德图： 1.576e010 s^2 H(s= ------------------------------------------------------------------------------------------ s^4 + 1.775e005 s^3 + 1.579e010 s^2 + 2.804e012 s + 2.494e014 (这是一个用butter函数产生的2阶的，频率范围为[20 20K]HZ的带通滤波器。我们可以用下面的语句：num=[1.576e010 0 0]; den=[1 1.775e005 1.579e010 2.804e012 2.494e014]; H=tf(num,den; bode(H 这样，我们就可以得到以下的伯德图： 可能我们会对这个图很不满意，第一，它的横坐标是rad/s，而我们一般希望横坐标是HZ；第二，横坐标的范围让我们看起来很不爽；第三，网格没有打开（这点当然我们可以通过在后面加上grid on解决）。下面，我们来看看如何定制我们自己的伯德图风格：在命令窗口中输入：bodeoptions 我们可以看到以下

内容：ans = Title: [1x1 struct] XLabel: [1x1 struct] YLabel: [1x1 struct]TickLabel: [1x1 struct]Grid: 'off' XLim: {[1 10]}XLimMode: {'auto'}YLim: {[1 10]} YLimMode: {'auto'}IOGrouping: 'none'InputLabels: [1x1 struct]OutputLabels: [1x1 struct]InputVisible: {'on'} OutputVisible: {'on'}FreqUnits: 'rad/sec'FreqScale: 'log' MagUnits: 'dB' MagScale: 'linear'MagVisible: 'on' MagLowerLimMode: 'auto'MagLowerLim: 0PhaseUnits: 'deg'PhaseVisible: 'on'PhaseWrapping: 'off' PhaseMatching: 'off'PhaseMatchingFreq: 0 PhaseMatchingValue: 0我们可以通过修改上面的每一 项修改伯德图的风格，比如我们使用下面的语句画我 们的伯德图：P=bodeoptions;P.Grid='on'; P.XLim={[10 40000]};P.XLimMode={'manual'};P.FreqUnits='HZ'; num=[1.576e010 0 0];den=[1 1.775e005 1.579e010 2.804e012 2.494e014];H=tf(num,den; bode(H,P 这时，我们将会看到以下的伯德图： 上面这张图相对就比较好了，它的横坐标单位 是HZ，范围是[10 40K]HZ，而且打开了网格，便于我 们观察-3DB处的频率值。当然，你也可以改变bodeoptions中的其它参数，做出符合你的风格的伯

典型环节的Bode图资料

典型环节的B o d e图

控制系统的开环频率特性 目的：掌握开环Bode图的绘制 根据Bode图确定最小相位系统的传递函数 重点：开环Bode图的绘制、根据Bode图确定最小相位系统的传递函数 1 开环伯德图手工作图的一般步骤： 1）将开环传递函数表示为时间常数表达形式，计算各个典型环节的交接频率 2）求20lgK的值，并明确积分环节的个数ν3）通过（1,20lgK）绘制斜率为-20vdB/dec 低频段 4）随着频率增加，每遇到一个典型环节的交接频率，就改变一次斜率 最小相位系统定义：递函数的零点、极点全部位于S 左半平面，同时又无纯滞后环节的系统称为最小相位系统。否则就是非最小相位系统。 对数幅频特性与相频特性之间存在确定的对应关系。对于一个最小相位系统，我们若知道了其幅频特性，它的相频特性也就唯一地确定了。也就是说：只要知道其幅频特性，就能写出此最小相位系统所对应的传递函数，而无需再画出相频特性。 非最小相位系统高频时相角迟后大，起动性能差，响应缓慢。对响应要求快的系统，不宜采用非最小相位元件。 Tf函数用来建立实部或复数传递函数模型或将状态方程、或零级增益模型转化成传递函数形式。sys = tf(num,den)命令可以建立一个传递函数，其中分子和分母分别为num和den。输出sys是储存传递函数数据的传递函数目标。单输入单输出情况下，num和den是s的递减幂级数构成的实数或复数行向量。这两个向量并不要求维数相同。如h = tf([1 0],1)就明确定义了纯导数形式h(s)=s。若要构建多输入多输出传递函数，要分别定义每一个单输入单输出系统的端口的分子与分母。 2 典型环节的伯德图 绘制曲线在MATLAB中实现，利用下述的程序段： num=[b2 b1 b0]; den=[1 a2 a1 a0];H=tf(num,den); bode(H) margin(H) hold on 2.1 比例环节 传递函数：() G s K = 频率特性：() G j K ω= 对数幅频特性：()20lg L j K ω= 对数相频特性：()0 ?ω= 程序段： num=[0 10]; den=[0 1]; H=tf(num,den); bode(H) margin(H) hold on 结论：放大环节的对数幅频特性是一条幅值为20lgK分贝，且平行于横轴的直线，相频特性是一条和横轴重合的直线。 K>1时，20lgK>0dB；K<1时，20lgK<0dB。2.2 惯性环节（低通滤波特性） 传递函数：1 () 1 G s sτ = + 频率特性：() ()()j G j A e?ω ωω = 对数幅频特性： 2 () 1() Lω τω = + 对数相频特性：()arctan ?ωτω =- 绘制1 () 10.1 G s s = + 的Bode图 程序段： num=[0 1]; den=[0.1 1];H=tf(num,den); bode(H) margin(H) hold on

使用simulink bode图的绘制

在Matlab中，大多时候，我们都是用M语言，输入系统的传递函数后，用bode函数绘制bode图对系统进行频率分析，这样做，本人觉得效率远不如Simulink建模高。如何在Matlab/Simulink中画bode图，以前也在网上查过些资料，没看到太多有用的参考。今天做助教课的仿真，又要画电机控制中电流环的bode图，模型已经建好，step response也很容易看出来，可这bode图怎么也出不来，又不愿意用m语言写出传递函数再画。baidu和google 了好一阵，几乎没有一个帖子说的清清楚楚的，经过一番摸索，终于掌握了Simulink里画bode图的方法。.其实，Simulink里画bode图，非常的easy，也很方便。写此文的目的是希望对那些常用Simulink进行仿真希望画bode图又不愿用M语言的新手有所帮助。 以下均是以Matlab R2008a为例。 首先，在simulink里建好model。如图1，这里需要注意的是，输入和输出要用input port 和output port，这样以后画bode图的时候，系统就会知道是这两个变量之间的关系。 图1 建好model 其次，选择线性分析。Tools->Control Design ->Linear Analysis。如图2。 图2 选择Linear Ansysis 将出现如图3所示的Control and Estimation Tools Manager窗口。

图3 Control and Estimation Tools Manager窗口 第三步，激动人心的时刻到了，哈哈。如果你是按照前面的步骤来的，那么这时候，你就应该可以直接画出bode图，在窗口的下方，将“Plot linear analysis result in a ”前面的方框打上勾，已打的就不用管了，再在后面的下拉框里选择“bode response plot”，即画output port和input port之间的bode图，再点击“Linearize Model”按钮，就OK了。其实除了bode图，还可以画其他很多响应曲线，比如step response、impulse response和Nyquist图等等，只需选择相应的step response plot，inpulse response plot或者Nyquist plot等等。方法都是相同的。选择选择“bode response plot”,如图4所示。 图4 画出bode图

第5章MATLAB绘图_习题答案

第5章MATLAB绘图 习题5 一、选择题 1．如果x、y均为4×3矩阵，则执行plot(x,y)命令后在图形窗口中绘制（）条曲线。D A．12B．7C．4D．3 2．下列程序的运行结果是（）。A x=0:pi/100:2*pi; forn=1:2:10 plot(n*sin(x),n*cos(x)) holdon end axissquare 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．两行两列的右上图 x的曲线绘制成直线，应采用的绘图函数是（）。C 5．要使函数y=2e A．polarB．semilogxC．semilogyD．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．ezmeshB．ezsurfC．ezplotD．plot3 8．下列程序运行后，看到的图形（）。C t=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);

第5章--MATLAB绘图-习题答案上课讲义

第5章--M A T L A B绘图-习题答案

第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:0.5:8);

BODE图 画图过程

电机定位系统校正(BODE图) MATLAB软件具有强大的计算能力和绘图功能，能够快速、准确地做出频域特性曲线。利用MATLAB绘制系统的Bode图，为控制系统设计和分析提供了极大的方便。 1. 创建M-file文挡，并输入如下程序，运行后生成LTI对象my_sys： J=3.2284e-6; b=3.5077e-6; K=0.0274; R=4; L=2.75e-6; num=[0 0 0 K]; den=[(J*K) (J*R+(L*b)) ((b*R)+K^2) 0]; my_sys=tf(num,den); 打开Matlab7.0软件，并新建一个空文档，将程序复制到文档内，如图1所示： 图1 2.运行程序并保存运行结果。如图2所示： 图2

3.打开Start-Toolboxes—Control System—SISO Design Tool。启动SISO Design，如图3所示 图3 4.将my_sys程序导入到SISO Design Tool中，如图4所示 图4

5.在View菜单中，关闭根轨迹显示，只显示开环的Bode图。如图5所示 图5 6. 加积分环节；加零点（60角频率）将各个参数进行积分：空白处右键—Add Pole/Zero—Integrator。如图6所示： 图6

7.在magnitude曲线加零点，然后Analysis菜单下Response to Step Command 指令。如图7所示： 图7 8.在管理反馈界面中，只显示闭环的r与y的关系—LT1 Viewer For SISO Design Tool界面空白处右键—Systems—Closed Loop ：r to u （green），如图8所示： 图8

第四讲 、Matlab绘图

第四讲 Matlab绘图 4.1 二维图形 4.2 数据分析图 4.3 三维图形 4.1 二维图形 1、基本图形的绘制 plot(x,y) 对向量x绘制向量y。以x为横坐标，y为纵坐标，按照坐标(xi ,yi)的有序排列绘制曲线。 plot(...,str) 使用字符串str指定的颜色和线型进行绘图。 例1：>> x=-pi:0.02*pi:pi; >> y=sin(x).*x.^2; >> plot(x,y) ezplot(f,xmin,xmax) 绘制函数f在区间[xmin,xmax]上的图形。如果省略xmin和xmax参数，区间将大概取在-2pi——2pi之间。由于ezplot命令使用算法来判断该函数变化显著的区间，因此区间的选取是不固定的。 例2：>> ezplot('sin(x).*x.^2') 2、图形控制 figure(gcf) 显示当前图形窗口。只键入figure命令则创建新的图形窗口； shg 显示当前图形窗口，等价于figure(gcf)。 hold on 保持当前图形。允许在当前图形状态下，使用同样的缩放比例加入另一个图形。 hold off 释放图形窗口，这样下一个图形将称为当前图形。这是缺省状态。 hold 在hold on和hold off之间进行切换。 subplot(m,n,p) 将图形窗口分割成m行n列,并设置p所指定的子窗口为当前窗口。子窗口按行由左至右，由上至下进行编号。这一命令在Matlab的当前版本中也被写作subplot(mnp)。axis…）用行向量中给出的值，设置坐标轴的最大和最小值。对于二维图形，该向量中含有元素： [xmin, xmax, ymin, ymax]。对于三维图形，是[xmin, xmax, ymin, ymax，zmin, zmax]。axis ~~ ~~的不同参数将给出不同的结果： 1.manual 固定坐标轴刻度。如果当前图形窗口为hold on状态，则后面的图形将采用同样的刻度 2.auto 把坐标轴刻度重新设置为缺省状态值。 3.equal 设置x轴和y轴为同样的刻度增量。 4.tight 以数据的大小为坐标轴的范围。 5.ij 翻转y轴，使得正数在下，负数在上。 6.xy 复位y轴，使正数在上。 7.off 坐标轴消隐。 8.on 绘制坐标轴。 title(txt) 在图形窗口顶端的中间位置输出字符串txt作为标题。 xlabel(txt) 在x轴下的中间位置输出字符串txt作为标注。 ylabel(txt) 在y轴边上的中间位置输出字符串txt作为标注。 zlabel(txt) 在z轴边上的中间位置输出字符串txt作为标注。 text(x,y,txt) 在图形窗口的(x,y)处写字符串txt。坐标x和y按照与所绘制图形相同的刻度给出。对于向量x和y，字符串txt写在(xi,yi)的位置上。如果t x t是一个字符串向量，即一个字符矩阵，且与x, y有相同的行数，则第i行的字符串将写在图形窗口的(xi,yi)的位置上。 gtext(txt) 通过使用鼠标或方向键，移动图形窗口中的十字光标，让用户将字串t xt放置在图形窗口中。当十字光标走到所期望的位置时，用户按下任意键或鼠标上的任意按钮，字符串将会写入在窗口中。

第5章--MATLAB绘图-习题答案

》 第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;

第二讲 Matlab编程与作图

第二讲Matlab编程与作图 第一部分Matlab程序设计初步 Matlab除了指令行操作的直接交互外，作为一种高级应用软件还提供了自己的编程语言。通过编写Matlab程序，可以更加方便地调用Matlab提供的各种功能强大的函数库，使得程序能完成复杂的运算处理大量的数值数据。 1、M文件简介 Matlab提供了丰富的编程语言，使得用户可以将一连串的命令写入文件，然后使用简单的函数来执行这些命令。文件被保存为文本文件，后缀为.m，比如说dblquad.m，因此Matlab的程序通常被称为M 文件。 M文件是一个文本文件，可以使用各种文本编辑器对它进行编辑和修改，比如Windows操作系统自带的记事本，也可以用Matlab 内建的M文件编辑器。 M文件分为两类，一类称为脚本（Scripts），类似于批处理文件，相当于将在Matlab命令窗口中执行的一系列指令放在一个文件中，当在命令窗口调用该文件名时，则按顺序执行其中的命令集。 例：编写求10！的程序。 另一类M文件称为函数（Function），它可以接受输入变量，并将运算结果送至输出变量，类似于数学中的函数y=f(x)。 函数M文件的基本结构： function f=fact(n) 函数定义行

%Compute a factorial value. %FACT(N) returns the factorial of N, 帮助文档%usually denoted by N! %Put simply,FACT(N) is PROD(1:N), 注释 f=prod(1:n); 函数体2、运算符 关系运算符：<, <=, >, >=, = =, ~= 逻辑运算符：与（&），或（|），非（~） 例：编写分段函数 21 () 1 -1<1 321 x x f x x x x ?> ? =≤ ? ?+≤- ? %myfun1.m function y=myfun1(x) y=(x.^2).*(x>1)+(x>-1& x<=1)+(3+2*x).*(x<=-1); 注意：1.函数名与变量名的命名法则相同，要求以字母开头，后接字母或下划线；2.函数名与保存的文件名最好一致。 3、控制流 所有的计算机编程语言都提供了控制程序流执行程序的语法，Matlab也不例外。所有的控制流语法都以end 结尾。 ⑴for 循环语句 语法：for 循环变量=数组 指令组；

MATLAB中bode图绘制技巧

Matlab中Bode图的绘制技巧 学术收藏 2010-06-04 21:21:48 阅读54 评论0 字号：大中小订阅 我们经常会遇到使用Matlab画伯德图的情况，可能我们我们都知道bode这个函数是用来画bode图的，这个函数是Matlab内部提供的一个函数，我们可以很方便的用它来画伯德图，但是对于初学者来说，可能用起来就没有那么方便了。 譬如我们要画出下面这个传递函数的伯德图： 1.576e010 s^2 H(s)= ------------------------------------------------------------------------------------------ s^4 + 1.775e005 s^3 + 1.579e010 s^2 + 2.804e012 s + 2.494e014 (这是一个用butter函数产生的2阶的，频率范围为[20 20K]HZ的带通滤波器。) 我们可以用下面的语句： num=[1.576e010 0 0]; den=[1 1.775e005 1.579e010 2.804e012 2.494e014]; H=tf(num,den); bode(H) 这样，我们就可以得到以下的伯德图： 可能我们会对这个图很不满意，第一，它的横坐标是rad/s，而我们一般希望横坐标是HZ；第二，横坐标的范围让我们看起来很不爽；第三，网格没有打开（这点当然我们可以通过在后面加上grid on解决）。

下面，我们来看看如何定制我们自己的伯德图风格： 在命令窗口中输入：bodeoptions 我们可以看到以下内容： ans = Title: [1x1 struct] XLabel: [1x1 struct] YLabel: [1x1 struct] TickLabel: [1x1 struct] Grid: 'off' XLim: {[1 10]} XLimMode: {'auto'} YLim: {[1 10]} YLimMode: {'auto'} IOGrouping: 'none' InputLabels: [1x1 struct] OutputLabels: [1x1 struct] InputVisible: {'on'} OutputVisible: {'on'} FreqUnits: 'rad/sec' FreqScale: 'log' MagUnits: 'dB' MagScale: 'linear' MagVisible: 'on' MagLowerLimMode: 'auto' MagLowerLim: 0 PhaseUnits: 'deg' PhaseVisible: 'on' PhaseWrapping: 'off' PhaseMatching: 'off' PhaseMatchingFreq: 0 PhaseMatchingValue: 0 我们可以通过修改上面的每一项修改伯德图的风格，比如我们使用下面的语句画我们的伯德图： P=bodeoptions; P.Grid='on'; P.XLim={[10 40000]}; P.XLimMode={'manual'}; P.FreqUnits='HZ'; num=[1.576e010 0 0]; den=[1 1.775e005 1.579e010 2.804e012 2.494e014]; H=tf(num,den); bode(H,P) 这时，我们将会看到以下的伯德图：

MatlabSimulink中bode图的画法

Matlab/Simulink中bode图的画法 在Matlab中，大多时候，我们都是用M语言，输入系统的传递函数后，用bode函数绘制bode图对系统进行频率分析，这样做，本人觉得效率远不如Simulink建模高。如何在Matlab/Simulink中画bode图，以前也在网上查过些资料，没看到太多有用的参考。今天做助教课的仿真，又要画电机控制中电流环的bode图，模型已经建好，step response也很容易看出来，可这bode图怎么也出不来，又不愿意用m语言写出传递函数再画。baidu和google了好一阵，几乎没有一个帖子说的清清楚楚的，经过一番摸索，终于掌握了Simulink里画bode图的方法。.其实，Simulink里画bode图，非常的easy，也很方便。写此文的目的是希望对那些常用Simulink进行仿真希望画bode图又不愿用M语言的新手有所帮助。 以下均是以Matlab R2008a为例。 首先，在simulink里建好model。如图1，这里需要注意的是，输入和输出要用input port和output port，这样以后画bode图的时候，系统就会知道是这两个变量之间的关系。 图1 建好model 其次，选择线性分析。Tools->Control Design ->Linear Analysis。如图2。

图2 选择Linear Ansysis 将出现如图3所示的Control and Estimation Tools Manager窗口。

图3 Control and Estimation Tools Manager窗口 第三步，激动人心的时刻到了，哈哈。如果你是按照前面的步骤来的，那么这时候，你就应该可以直接画出bode图，在窗口的下方，将“Plot linear an alysis result in a ”前面的方框打上勾，已打的就不用管了，再在后面的下拉框里选择“bode response plot”，即画output port和input port之间的bode图，再点击“Linearize Model”按钮，就OK了。其实除了bode图，还可以画其他很多响应曲线，比如step response、impulse response和Nyquist图等等，只需选择相应的step response plot，inpulse response plot或者Nyquist plot等等。方法都是相同的。选择选择“bode response plot”,如图4所示。 图4 画出bode图 稍等片刻，便出现了图1中output port和input port的bode图了。是不是很简单？！

第5章 MATLAB绘图_习题答案

第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:0.5:8); z=sin(sqrt(x.^2+y.^2))./sqrt(x.^2+y.^2+eps);

matlab绘制bode图技巧

我们经常会遇到使用Matlab画伯德图的情况，可能我们我们都知道bode这个函数是用来画bode图的，这个函数是Matlab内部提供的一个函数，我们可以很方便的用它来画伯德图，但是对于初学者来说，可能用起来就没有那么方便了。 譬如我们要画出下面这个传递函数的伯德图： 1.576e010 s^2 H(s)= ------------------------------------------------------------------------------------------ s^4 + 1.775e005 s^3 + 1.579e010 s^2 + 2.804e012 s + 2.494e014 (这是一个用butter函数产生的2阶的，频率范围为[20 20K]HZ的带通滤波器。) 我们可以用下面的语句： num=[1.576e010 0 0]; den=[1 1.775e005 1.579e010 2.804e012 2.494e014]; H=tf(num,den); bode(H) 这样，我们就可以得到以下的伯德图： 可能我们会对这个图很不满意，第一，它的横坐标是rad/s，而我们一般希望横坐标是HZ；第二，横坐标的范围让我们看起来很不爽；第三，网格没有打开（这点当然我们可以通过在后面加上grid on解决）。 下面，我们来看看如何定制我们自己的伯德图风格： 在命令窗口中输入：bodeoptions 我们可以看到以下内容： ans = Title: [1x1 struct] XLabel: [1x1 struct]

YLabel: [1x1 struct] TickLabel: [1x1 struct] Grid: 'off' XLim: {[1 10]} XLimMode: {'auto'} YLim: {[1 10]} YLimMode: {'auto'} IOGrouping: 'none' InputLabels: [1x1 struct] OutputLabels: [1x1 struct] InputVisible: {'on'} OutputVisible: {'on'} FreqUnits: 'rad/sec' FreqScale: 'log' MagUnits: 'dB' MagScale: 'linear' MagVisible: 'on' MagLowerLimMode: 'auto' MagLowerLim: 0 PhaseUnits: 'deg' PhaseVisible: 'on' PhaseWrapping: 'off' PhaseMatching: 'off' PhaseMatchingFreq: 0 PhaseMatchingValue: 0 我们可以通过修改上面的每一项修改伯德图的风格，比如我们使用下面的语句画我们的伯德图：P=bodeoptions; P.Grid='on'; P.XLim={[10 40000]}; P.XLimMode={'manual'}; P.FreqUnits='HZ'; num=[1.576e010 0 0]; den=[1 1.775e005 1.579e010 2.804e012 2.494e014]; H=tf(num,den); bode(H,P) 这时，我们将会看到以下的伯德图：

matlab绘制bode图技巧(可编辑修改word版)

我们经常会遇到使用Matlab 画伯德图的情况，可能我们我们都知道bode 这个函数是用来画bode 图的，这个函数是Matlab 内部提供的一个函数，我们可以很方便的用它来画伯德图，但是对于初学者来说，可能用起来就没有那么方便了。 譬如我们要画出下面这个传递函数的伯德图： 1.576e010 s^2 H(s)= s^4 + 1.775e005 s^3 + 1.579e010 s^2 + 2.804e012 s + 2.494e014 (这是一个用butter 函数产生的2 阶的，频率范围为[20 20K]HZ 的带通滤波器。) 我们可以用下面的语句： num=[1.576e010 0 0]; den=[1 1.775e005 1.579e010 2.804e012 2.494e014]; H=tf(num,den); bode(H) 这样，我们就可以得到以下的伯德图： 可能我们会对这个图很不满意，第一，它的横坐标是rad/s，而我们一般希望横坐标是HZ；第二，横坐标的范围让我们看起来很不爽；第三，网格没有打开（这点当然我们可以通过在后面加上grid on 解决）。 下面，我们来看看如何定制我们自己的伯德图风格： 在命令窗口中输入：bodeoptions 我们可以看到以下内容： ans = Title: [1x1 struct] XLabel: [1x1 struct]

YLabel: [1x1 struct] TickLabel: [1x1 struct] Grid: 'off' XLim: {[1 10]} XLimMode: {'auto'} YLim: {[1 10]} YLimMode: {'auto'} IOGrouping: 'none' InputLabels: [1x1 struct] OutputLabels: [1x1 struct] InputVisible: {'on'} OutputVisible: {'on'} FreqUnits: 'rad/sec' FreqScale: 'log' MagUnits: 'dB' MagScale: 'linear' MagVisible: 'on' MagLowerLimMode: 'auto' MagLowerLim: 0 PhaseUnits: 'deg' PhaseVisible: 'on' PhaseWrapping: 'off' PhaseMatching: 'off' PhaseMatchingFreq: 0 PhaseMatchingValue: 0 我们可以通过修改上面的每一项修改伯德图的风格，比如我们使用下面的语句画我们的伯德图：P=bodeoptions; P.Grid='on'; P.XLim={[10 40000]}; P.XLimMode={'manual'}; P.FreqUnits='HZ'; num=[1.576e010 0 0]; den=[1 1.775e005 1.579e010 2.804e012 2.494e014]; H=tf(num,den); bode(H,P) 这时，我们将会看到以下的伯德图：

MATLAB中bode图绘制技巧(精)

M A T L A B中b o d e图绘制技 巧(精) -标准化文件发布号：（9556-EUATWK-MWUB-WUNN-INNUL-DDQTY-KII

Matlab中Bode图的绘制技巧 学术收藏 2010-06-04 21:21:48 阅读54 评论0 字号：大中小订阅 我们经常会遇到使用Matlab画伯德图的情况，可能我们我们都知道bode这个函数是用来画bode图的，这个函数是Matlab内部提供的一个函数，我们可以很方便的用它来画伯德图，但是对于初学者来说，可能用起来就没有那么方便了。 譬如我们要画出下面这个传递函数的伯德图： 1.576e010 s^2 H(s= ------------------------------------------------------------------------------------------ s^4 + 1.775e005 s^3 + 1.579e010 s^2 + 2.804e012 s + 2.494e014 (这是一个用butter函数产生的2阶的，频率范围为[20 20K]HZ的带通滤波器。 我们可以用下面的语句： num=[1.576e010 0 0]; den=[1 1.775e005 1.579e010 2.804e012 2.494e014]; H=tf(num,den; bode(H 这样，我们就可以得到以下的伯德图：

可能我们会对这个图很不满意，第一，它的横坐标是rad/s，而我们一般希望横坐标是HZ；第二，横坐标的范围让我们看起来很不爽；第三，网格没有打开（这点当然我们可以通过在后面加上grid on解决）。 下面，我们来看看如何定制我们自己的伯德图风格： 在命令窗口中输入：bodeoptions 我们可以看到以下内容： ans = Title: [1x1 struct] XLabel: [1x1 struct] YLabel: [1x1 struct] TickLabel: [1x1 struct] Grid: 'off' XLim: {[1 10]} XLimMode: {'auto'} YLim: {[1 10]} YLimMode: {'auto'} IOGrouping: 'none' InputLabels: [1x1 struct] OutputLabels: [1x1 struct] InputVisible: {'on'} OutputVisible: {'on'} FreqUnits: 'rad/sec' FreqScale: 'log' MagUnits: 'dB'