matlab2014里simulink画bode图

Version MA TLAB 2014a

1 首先在simulink里建好model，输入输出要用input port和output port。这样系统会识别它们之间的输入输出关系。如图1所示：

图 1 model

2 依次选择工具栏上的Analysis、control design、linear analysis。如图2所示：

图2 打开Linear Analysis

3 选择分析端口为Root level inports and outports。如图3所示：

图3 选择端口

4 选择画出图的类型，New Step为在一个新的图中绘出阶跃响应，New Bode为在一个新的图中绘出bode图。EXISTING PLOTS表示在一个已经存在的图中画出bode或其它图，可以将多个bode放在同一个图中用来对比波形。这里能画出很多曲线，可以根据需要选择。

5 显示设置，根据需要选择。添加网格，图例。

6 也可以发送到figure对画出的图进行更加详细的编辑。

7 查看状态空间方程，零极点增益和传递函数。

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图

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)+O2) 0]; my_sys=tf( nu m,de n); 打开Matlab7.0软件，并新建一个空文档，将程序复制到文档内，如图1所示： 图1 2.运行程序并保存运行结果。如图2所示: 图2

图4 3. 打开 Start-Toolboxes — Control System — SISO Design Tool 启动SISO Design,如图3所示 4. 将my_sys 程序导入到SISO Design Tool 中， 如图4所示 啪号TW Hi^vicn Toni P il ■ ErLE ? Vi, mr r za-ipMi-i k t Dqri ￡x> al > 1* E □'l'l. ?冷日■丁 11*』]1卜 l>] X o 4 T 11; ?M s K 1 m i 4||'=4Dp 山 watL a ft — 15-S-IQ Jt 午肌 "■存 -i-s — 15-5-in 上午恥 ny_2F2 c !x| F L 1? Edu I 上LIM tmp viiK?l&rl 姿kl 尸L C. Tadul llaJp II H n y. L R ann 曲 闻 出田画田田刚 Dur r ■n.l 卫a r nry CtiaiTijinfl Hi ?f -15-5-JQ 上牛奔 ■■■jip.syg G ii — 15-5-10 上午和 i SB (TT = EJS > 尺jgtrt.Hci nr ||T ^ nkm 1rn mfi B * njrtlini ￡>■ |si90 ?ui. fi 1 v f Ja S D .EM 1 色 ■'i -Fi'Mii ^mpp+if-slnr Ci|r TBH.*: n i. B-Bf I arv ^(Bft-olKl fir tiv? for “偷史 俺±&n )埔1口 曲 riw^uefCf rijSr.1 4 fftcrt

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图了。是不是很简单？！

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) 这时，我们将会看到以下的伯德图：

波特图的画法

二、 对数频率特性 假设：) ()()(ω?ωωj e j H j H =。对其取对数： [][] [])()()()(ln )(ln )(ln ) (ω?ωω?ωωωω?j G j j H e j H j H j +=+== 其虚部正是系统的相频特性，而实部： [])(ln )(ωωj H G = 称为对数增益，反映了系统幅频特性，单位奈培（Np, Neper ）。 一般情况下不用自然对数，而取常用对数，定义： [])(log 20)(ωωj H G = 单位：分贝(Deci-Bel,dB)。 奈培与分贝的转换关系：1 Np = 8.686 dB 在理论分析中，一般使用Np ；在实际应用中，一般使用dB 用分贝表示增益，解决了信号动态范围与精度之间的矛盾。如果在频率坐标中同样使用对数坐标，则同样可以解决频率的范围与精度之间的矛盾。这样一来就形成了波特图。 ? 波特图的横坐标可以用ωlog ，也可以用f log ; ? 在波特图的横坐标上，一般直接标注频率值； ? 波特图的横坐标上只能表示0>ω或者0>f 频率下的系统特性。 图中的二、三象限并非表示频率小于零的部分，而是表示频率小于1（大于零）部分频率特性。 ? 根据系统频率特性的共扼对称性,不难得到频率小于零部分的 特性。 ? 在波特图的纵坐标上，可以标注系统幅频特性值（如图中红字所 示），也可以标注分贝值。 ? 为了方便参数的判读，实际工程中的波特图中的刻度也不是按照等 间隔设置的，而是按照对数间隔设置。例如下图。

有专用的对数坐标图纸可以用于手工绘制波特图。 波特图的纵坐标上同样也只表示了系统幅频特性中大于零的部分。 图中的三、四象限并非表示系统的幅频特性小于零，而是表示系统的幅频特性小于1（大于零）。 三、 线性系统的波特图 1、一般系统的波特图 ??? ? ??-==∑ ∑==∏∏--=n i i m i i j n i i m i i e p j z j H j H 111 10 )(αβωωω ∑∑∑∑====-+=---+== =n i pi m i zi n i i m i i G G H p j z j H j H G 1 1 01 10) ()(log 20log 20log 20log 20)(log 20)(ωωωωωω 所以，不仅系统的相频特性是各个零点或极点的相频特性的叠加，而且系统的幅频特性是各个零点或极点的相频特性的叠加。所以，可以根据各个零点或极点的波特图的叠加得到系统的波特图。 2、一次因式的波特图 1） 单个零点的波特图： )1(1 )1()(i i i i i zi T j T z j z z j j H ωω ωω+= -+-=-= （1）幅频特性 ()[] 2 1log 10log 201log 201 log 20)(log 20)(i i i i zi zi T T T j T j H j G ωωωω++-=++==

bode图怎么绘制

不过在系统辨识工具箱中有相应的函数能够完成该功能。 应该是invfrqs函数。但这个是由频率响应得到的，所以你的对输入输出信号做傅里叶变换得到频率响应后在用这个函数拟合传递函数。 使用matlab工具箱更为方便和直观： 1. 把u，y信号导入到工作空间里。 2. 用ident命令打开matlab系统辨识工具箱，然后点击import data，从新打开界面里导入工作空间的数据。然后可以通过图形查看该输入输出信号，或者在proprocess进行信号预处理。 3. 根据你的模型在estimate里选择linear parameter models，个人觉得你应该选择ARX 结构，确定阶数，然后进行估计。 4. 在主界面里查看估计模型，并且可以和实际输出比较，看看拟合度。 详细使用方面参考帮助文档System Identification Toolbox User's Guide 也只能是拟合吧，拟合的目标也只能是有限的一些典型传递函数。 做过这个，就是拟合吧 在sisotool工具里放置合适的零极点，尽量使曲线吻合，得到传递函数 我是大概知道实测的曲线的零极点分布的，然后去拟合的，这样好做点 如果已经得到Bode plot，就幅度vs频率，相位vs频率曲线，根据+/-20db/dec,+/-40db/dec 渐近线先预测几个pole,zero,double pole double zero，然后再去近似了。

还有几个比较困难的地方是： 1）实际系统引入的杂散参数和群延迟特性，才是拟合和预测的难点。 2）实际测量时仪器已经采用了离散化的数字处理手段，还有测量系统处理误差和测量时注入的扰动幅度所引起的误差，对低频还好，对高频都是很大的影响。 3）所以目前已有的模型，多是对1/2开关频率之前的预测，V2等新模型有所改进，但是还是有其局限性 以我测试的为例吧， R=1e3; C=1e-7; L=0.1; Rl=1000; magdb_LCR=zeros(1,1e5);% -sL-R-|(1/sC) mag_LCR=zeros(1,1e5); phase_LCR=zeros(1,1e5); freq=(1:1e5); RCplot=zeros(1,1e5); LCRplot=zeros(1,1e5); for f=1:1:1e5 LCR=1/((1i*2*pi*f)^2*L*C+1i*2*pi*f*C*Rl+1); LCRplot(f)=LCR; mag_LCR(f)=abs(LCR); magdb_LCR(f)=20*log10(abs(LCR)); phase_LCR(f)=angle(LCR)*180/pi; end 我用上方的函数写一个已知的传递函数，生成对应频率下的幅值和相位。频率响应曲线如下，需要注意的是，幅值的纵坐标不是dB。

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'

bode图画图过程

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

图3 4. 将my_sys程序导入到SISO Design Tool中，File-Import 如图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 9. 添加零点和极点，如图8所示，并移动极点、零点和线的位置，调整LT1 Viewer For SISO Design Tool窗口中函数图像的变化直到符合Bode图，如图9所示：

新版本matlab中绘制bode图技巧

百度文库里有一篇《matlab绘制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) 这样，我们就可以得到以下的伯德图： 下面对其进行更改，使得： 一、横坐标是HZ； 二、调整横坐标的范围； 三、打开网格（这点当然我们可以通过在后面加上grid on解决）。 下面，我们来看看如何定制我们自己的伯德图风格： 在较新版本的matlab的命令窗口中输入： num=[1.576e010 0 0]; den=[1 1.775e005 1.579e010 2.804e012 2.494e014]; H=tf(num,den); h=bodeplot(H);

p=getoptions(h); p.Grid=’on’; p.XLim={[10 40000]}; p.XLimMode={'manual'}; p.FreqUnits='Hz'; setoptions(h,p); 这时，我们可以得到以下的伯德图：