基于MATLAB在自动控制域 中稳定性分析的应用
学院：装备工程学院
学号： 姓名：
沈阳理工大学MATLAB论文
基于MATLAB在自动控制频域中稳定
性分析的应用
摘 要: 自动控制系统主要利用MATLAB高级语言对其进行计算机分析。 MATLAB是一套高性能的数值计算和可视化软件,它集数值分析、矩阵运算、信号处理和图形绘制集于一体,构成了一个方便的、界面友好的用户环境。本文主要介绍了利用MATLAB在自动控制中对频域中系统稳定性的判定，通过MATLAB建立某一系统的模型,并分析该系统的性能。根据响应曲线判断系统的稳定性，当系统的性能不能满足所要求的性能指标时,通过调整系统参数和增添校正装置来改善系统性能并展示方便灵活的动态仿真结果。 关键词:自动控制系统；稳定性；频域分析；频率响应；稳定裕度
引言
频域分析法是应用频域特性研究线性控制系统的一种经典方法，采用这种 方法可以直观的表达出系统的频率特性，利用系统的传递函数绘制系统的bode 图、nyquist曲线和nichols图，然后进行系统稳定性的判定，这样可以很明确 的判定系统频域的稳定性。频域分析法是自动控制领域中应用又一种数学工具 ———频率特性来研究系统控制过程性能，即稳定性、快速性及稳态精度的一种 方法。这种方法不必直接求解系统的微分方程，而是间接的运用系统的开环频率 特性曲线，分析闭环系统的响应，因此它是一种图解的方法。本文介绍了应用 MATLAB在bode图、nyquist曲线和nichols图等控制系统频域分析中的主要方 法，通过具体实例叙述了MATLAB在频域分析中的应用过程。频域分析里主要用 到三种曲线（或叫图）：Bode图、Nyquist曲线图和 Nichols（尼柯尔斯）曲线 图。这三种曲线就是频率分析的三种工具。Bode图可以用于分析相角稳定裕度、 幅（或模）值稳定裕度、剪切频率、??穿越频率、带宽、扰动抑制及其稳定性 等，所以Bode图在频域分析里占有重要的地位。Nyquist与Nichols曲线图在频域分析里也很有用。本文中最重要的函数命令有bode、nyquist、pade、 nichols、margin等。
一.频域分析法的基础
1.有关频率分析的几个概念
2
沈阳理工大学MATLAB论文
（1）频率响应
当正弦函数信号作用于线性系统时，系统稳定后输出的稳态分量仍然是同频率的正弦信号，这种过程叫做系统的频率响应。 （2）频率特性
设有稳定的线性定常系统，在正弦信号作用下，系统输出的稳态分量为同频率的正弦信号，其振幅与输入正弦信号振幅比相对于正弦信号角频率间的关系
A(?)叫做幅频特性；其相位与输入正弦信号的相位之差相对于正弦信号角频率
间的关系?(?)叫做相频特性。系统频率响应与输入正弦信号的复数比叫做系统的频率特性。记作： G(j?)?A(?)ej?(?) (1) 2. Bode图
Bode图即对数频率特性曲线。Bode图有两条曲线，分别是对数幅频特性与对数相频特性。横坐标都是角频率?，是按常用对数（以10为底的对数）刻度的即lg?，其单位是rad/s。对数相频特性的纵坐标是?(?)，等分刻度，其单位为角度或弧度（rad）；而对数幅频特性的纵坐标是 L(?)?20lgA(?)，也是等分刻度的，其单位为dB。 3.幅相选择性
系统的频率特性G(j?)?A(?)ej?(?)里既有振幅信息又有相位信息，所以又叫做系统的幅相特性。幅相特性图形化的形式，即是幅相特性曲线。 4.频率特性的极坐标图
频率特性G(j?)是?的复变函数，在G(j?)复平面上对于某一?可以用以向量或其端点（坐标）来表示。当?从0??时，G(j?)端点的极坐标轨迹即是频率特性的极坐标图，叫做Nyquist曲线图或幅相特性曲线。 5.Nichols曲线图.
对于单位负反馈系统其闭环与开环传递函数G(s)间的关系为?(s)?G(j?) (3)
1?G(j?)3
G(s)(2)
1?G(s)则频率特性之间的关系为 ?(j?)?沈阳理工大学MATLAB论文
若以模幅式G(j?)表示为 G(j?)?Gej?G (4) 则 ?(j?)?M(?)ej?(?) (5) 将闭环幅频M(?)与闭环相频?(?)合成幅相特性为M(?)ej?(?)?由此幅相特性可以求得两个函数表达式，即
M?f1(G,?G) (7)
??f2(G,?G) (8)
这两个函数表达式的图形化曲线就叫做Nichols曲线图。 6.频域性能指标
1.）峰值：它是幅频特性A(?)的最大值。
2．）频带：它是幅频特性A(?)的数值衰减到0.707A(0)时对应的角频率。 3.）相频宽：它是相频特性?(?)等于??/2时对应的角频率。
4.）剪切频率：系统开环对数幅频特性曲线20lgG与横坐标轴(?)交点的角频率，常用?c（或?cp ）来标识。
5.）?? 穿越频率：系统开环对数相频特性曲线?(?)与?? 线交点所对应的角频率（即开环幅相特性曲线G(j?)与负实轴的交点所对应的角频率），常用?g （或?cg ）来标识。 6.）稳度裕度
1.相角稳定裕度：系统开环幅相特性曲线G(j?)上模值等于1的向量与负实轴的夹角，即 ???(?c)?(??) (9)
2．幅值稳定裕度：开环幅相特性曲线G(j?)与负实轴交点(?cg)模值G(?cg)的倒数，即 h?1G(?cg)Gej?G1?Gej?G (6)
(10)
或者是其交点模值倒数的分贝值，即 Lh?20lgh (11)
4
沈阳理工大学MATLAB论文
7.三频段概念
系统开环对数振幅频率特性曲线20lgG按横坐标大致分为3个频段。在开环对数振幅特性曲线20lgG上第一个转折频率以左的区段，是由积分环节与开环增益所确定的，此为低频段；在20lgG上第一个转折频率以右的剪切频率?c附近的区段叫做中频段，这段特性集中反映了闭环系统动态响应的稳定性与快速性。在中频段以右??10?c的区段叫做高频段，此部分是由系统中时间常数很小、频带很高的元件决定的，系统在高频段的分贝值，反映了系统对输入高频干扰信号的抑制能力。 8.Nyquist稳定判据
系统开环幅相特性G(j?)如果满足以下关系式：??[1?G(j?)]?p?
?:0??则系统闭环稳定。式中，P为系统开环特征方程的不稳定根个数。 9.对数频率稳定判据
当P=0时，在开环对数幅相特性曲线20lgG?0的范围内，相频特性曲线
?(?)对??线的正穿（由下向上）次数与负穿（由上向下）次数相等，则系统闭环稳定；当P?0时，在开环对数幅相特性曲线20lgG?0的范围内，若相频特性曲线?(?)对??线正穿次数与负穿次数之差为P/2，则系统闭环稳定。
二．控制系统的频域特性
当正弦函数信号作用于线性系统时，系统稳定后输出的稳态分量仍然是同频率的正弦信号，这种过程叫做系统的频率响应。根据自动控制原理，系统输出的
G(j?)?A(?)ej?(?)稳态分量与输入正弦信号的复数比叫做系统的频率特性记作：
系统的频率特性与系统传递函数之间有着简单而直接的关系：
G(j?)?G(s)s?j? (12)
1.若系统的单位阶跃响应为c(t)?1?1.5e?3t?1.2e?5t （t?0） 。试确定系统的频率响应。
其实现的程序代码如下：
5