三阶系统综合分析与设计

1、武汉理工大学自动控制理论课程设计说明书目录摘要11 设计内容21.1 设计题目21.2 设计任务2 绘制三阶系统的根轨迹32.1 常规方法绘制根轨迹32.2用MATLAB绘制根轨迹43 不同条件下K的取值53.1 当-8为闭环系统的一个极点时，K的取值53.2 主导极点阻尼比为0.7时的k值54 求系统的稳态误差64.1 位置误差系数74.2 速度误差系数74.3 加速度误差系数84.4 输入信号为时的稳态误差85 绘制单位阶跃响应曲线86 频域特性分析96.1绘制Bode图和Nyquist曲线96.2相角裕度和幅值裕度127 加入非线性环节判断稳定性137.1 求死区特性环节的描述函数137

2、.2 根据负倒描述函数和Nyquist图判断系统的稳定性148 设计体会15参考文献.17 摘要三阶系统是以三级微分方程为运动方程的控制系统。在控制工程中，三阶系统非常普遍，其动态性能指标的确定是比较复杂。在工程上常采用闭环主导极点的概念对三阶系统进行近似分析，或直接用MATLAB软件进行高阶系统分析。在课程设计中，要掌握用MATLAB绘制闭环系统根轨迹和系统响应曲线，用系统的闭环主导极点来估算三系统的动态性能，以及在比较点与开环传递函数之间加一个非线性环节判断其稳定性。1 设计内容1.1 设计题目三阶系统的综合分析和设计初始条件：某单位反馈系统结构图如图1-1所示： 图1-1 图1-21.2

3、 设计任务要求完成的主要任务: （包括课程设计工作量及其技术要求，以及说明书撰写等具体要求）1、试绘制随根轨迹2、当-8为闭环系统的一个极点时，K=?3、求取主导极点阻尼比为0.7时的K值（以下K取这个值）4、分别求取位置误差系数、速度误差系数、加速度误差系数及输入信号为单位阶跃信号、斜坡信号及单位加速度信号时的稳态误差5、用Matlab绘制单位阶跃相应曲线6、绘制Bode图和Nyquist曲线，求取幅值裕度和相角裕度7、如在比较点与开环传递函数之间加1个死区非线性环节，如图1-2所示，其中，试求取非线性环节的描述函数，并根据负倒描述函数和Nyquist图判断系统的稳定性8、认真撰写课程设计报

4、告。 绘制三阶系统的根轨迹图2 系统结构图由图1可得，三阶系统的开环传递函数为：G（s）= 2.1 常规方法绘制根轨迹根据绘制根轨迹的规则，可知该系统的根轨迹绘制步骤如下：（1）根轨迹的起点和终点。根轨迹起于开环极点，终于开环零点。开环极点分别为0、-3、-6，无开环零点。（2）根轨迹的分支数。n=3，m=0，所以分支数为3 。且它们是连续的并且对称于实轴。（3）根轨迹的渐进线。本系统根轨迹的渐近线有三条，据其与实轴的夹角公式: 把n=3,m=0代入求得：渐近线与实轴的交点为： (4)根轨迹在实轴上的分布。实轴上的某一个区域，若其右边开环零、极点的个数之和为奇数，则该区域必是根轨迹。因此实轴上

5、-，-6 -3,0必为根轨迹。(5)确定根轨迹的分离点。该系统中没有有限零点，由法则五得：于是分离点方程为：因此可以求得分离点d=-1.3,d=-4.7（不合题意，舍去）(6)根轨迹与虚轴的交点。闭环特征方程式为 对上式应用劳斯判据，有： 1 18 9 k k令劳斯表中行的首项为零，得k= 162,根据行的系数，得辅助方程 9+k=0代k=162并令s=jw,解得交点坐标2.2用MATLAB绘制根轨迹MATLAB为绘制根轨迹编程如下：num=1; den=1 9 18 0; syms=tf(num,den); rlocus(syms)MATLAB产生的根轨迹如图2所示： 图3 闭环系统的根轨迹

6、3 不同条件下K的取值3.1 当-8为闭环系统的一个极点时，K的取值由图1的系统的闭环传递函数为 闭环特征方程式为 把s=-8代入上式中解得 K=803.2 主导极点阻尼比为0.7时的k值 在控制工程实践中，通常要求控制系统即具有较快的响应速度，又具有一定的阻尼程度，此外还要求减少死区间、间隙和库仑摩擦等非线性因素对系统性能的影响，因此高阶系统的增益常常调整到是系统具有一对闭环共轭主导极点。这时，可以用二阶系统的动态性能指标来估算高阶系统的性能。 由于主导极点阻尼比=0.7<1,属于欠阻尼系统。由公式得： = 设系统的自然频率为，阻尼比，由上述用闭环主导极点分析高阶系统的方法可知，距虚轴

7、最近的一对闭环共轭主导极点为： 代入数据： 闭环特征方程式为 代入s的方程化简得: 0.686-12.6+k+(0.686-8.82+12.6)=0 分别令实部和虚部为零得到两个方程： 0.686-12.6+k=00.686-8.82+12.6=0解得 （舍） 所以取。所以，开环传递函数为4 求系统的稳态误差 系统的误差 e(t)一般定义为输出量的希望值与实际值之差，一般情况下采用从系统输入端定义的误差e(t)来进行计算分析。控制系统的稳定系统误差信号的稳态分量称为系统的稳态误差，以表示 4.1 位置误差系数对于单位阶跃输入，R(s)=1/s,求得系统的稳态误差为 令，称为稳态位置误差系数。稳

8、态误差可表示为 对于0型系统，=0，则=；当1时，=。 由图1系统得： 由上式得=1，则=04.2 速度误差系数 对于单位斜坡输入，此时系统的稳态误差为 称为稳态速度误差系数。 于是稳态误差可表示为 由图1系统得：=0.978=1.022 4.3 加速度误差系数对于单位抛物线输入，此时系统的稳态误差为 令 称为稳态加速度误差系数。于是稳态误差可表示为于是稳态误差可表示为 则对于图1系统得：=0=4.4 输入信号为时的稳态误差系统的开环传递函数为 (=17.617)当时，=，=0 当时，=0.978，=2.556当时，=0，=则当输入信号为时 =+=5 绘制单位阶跃响应曲线MATLAB绘制单位阶

9、跃响应曲线编程如下：num=17.6; den=1 9 18 17.6;syms=tf(num,den);step(tf(num,den) 即系统单位阶跃相应曲线为图3所示： 图4 系统单位阶跃相应曲线6 频域特性分析6.1绘制Bode图和Nyquist曲线6.1.1 绘制Bode图Bode图又称对数频率特性曲线图，由对数幅频曲线和对数相频曲线组成，是工程中广泛使用的一组曲线。手工绘制伯德图的步骤如下：1、将传递函数写成伯德标准型，确定开环传递系数和各转折频率。 系统的开环频率特性为 由伯德标准型容易看书，开环传递系数为，转折频率为，。2、 确定低频段：由传递函数可知该系统为1型系统，即有微分

10、环节，所以绘制低频段，可过， 作一条斜率为-20dB/dec的斜线。3、 绘制开环对幅频特性的渐近线：将低频段延伸到第一个转折点频率处。因为第一个转折频率是惯性环节的转折频率，所以，开环对数频率特性的渐近线下降20dB/dec,再延伸到第二个转折频率处，因为也是惯性环节，所以再下降20dB/dec。4、 绘制相频特性：绘制各个环节的对数相频特性曲线，然后逐点叠加。一般在一些特征点上进行叠加，如各个转折频率处。5、 修正对数幅频特性。Matlab编写程序时可以利用函数Bode和margin,但是margin函数可以直接读出相角裕度和幅值裕度，所以下面用margin函数绘制BODE图MATLAB绘

11、制BODE图编程如下：num=17.6; den=1 9 18 0; syms=tf(num,den);margin(num,den)绘制出的BODE图如图2-3所示： 图5 开环传递函数的Bode图从图中可以看出：幅值裕度h=19.3dB,穿越频率=4.24rad/sec。 相角裕度r=64deg，截止频率=0.924rad/sec。从幅值裕度和相角裕度可以看出这是一个稳定的系统。6.1.2 绘制Nyquist曲线Nyquist曲线图又称开环幅相频率特性曲线图，是利用Nyquist判别系统稳定性的依据。MATLAB绘制Nyquist曲线编程如下：num=17.6; den=1 9 18 0;

12、 syms=tf(num,den);nyquist(syms)绘制出的Nyquist曲线如图图6 开环传递函数Nyquist曲线 6.2相角裕度和幅值裕度相角裕度的含义是，对于闭环稳定系统的，如果系统开环相频特性再后度，则系统将处于临界稳定状态。设为系统的截止频率，显然：由上式可得： =1解得： =0.92 rad/s相角裕度： =64.1幅值裕度的含义是，对于闭环稳定系统，如果系统开环幅频特性再增大倍，则系统将处于临界稳定状态。设为系统的穿越频率，则在处的相角： = ； 由上式得 =180 4.24 rad/s幅值裕度为 =19.3上述为笔算结果，通过MATLAB绘制的Bode图中得到的相角

13、裕度 为，幅值裕度为19.3，两种方法得到结果基本吻合。7 加入非线性环节判断稳定性 图7 加非线性环节系统结构图在比较点与开环传递函数之间加一个死区特性非线性环节，如图6所示。7.1 求死区特性环节的描述函数由正弦输入信号、死区特性可得死区特性环节输出的数学表达式为 由于为奇函数，所以=0,而又为半周期内对称，故 死区特性的描述函数为 代入得： 即得到描述函数为 7.2 根据负倒描述函数和Nyquist图判断系统的稳定性 据图2得，对于线性环节， 解得穿越频率为： 非线性环节为死区特性，负倒描述函数为即的曲线如图7中横轴射线所示。即曲线此时在实轴上由-指向-0.5的直线，而此时奈氏曲线过负实

14、轴交于-0.02833，没有包围曲线，根据非线性系统稳定判据（若曲线不包围负倒数描述函数曲线，则非线性系统稳定;若包围，则系统不稳定），该系统稳定。图8 负倒描述函数曲线与的Nyquist曲线8 设计体会通过这次自动控制原理课程设计，我学到了关于MATLAB方面的很多知识，MATLAB它有着强大的数据处理能力，处理速度快，精度高，它不仅可以用来绘制曲线，而且可以用来帮助解方程，以及做仿真处理，帮助验证理论分析的真确性。这次我做的是一个三阶系统的综合分析和设计，主要完成的任务是绘制根轨迹绘制单位阶跃响应求稳态误差绘制BODE图和Nyquist曲线求幅值裕度和相角裕度以及加入非线性死区利用负倒描述

15、函数和Nyquist曲线来判断系统的稳定性。这些问题在课堂上老师都讲过，但是做起课程设计时还是遇到了一点困难，对于加入非线性死区利用负倒描述函数和Nyquist曲线来判断系统的稳定性这方面我比较陌生。但是通过翻阅资料，认真分析，我最终还是解决了问题。根轨迹法是一种十分便捷得分析和设计线性定常控制系统的图解方法，了解到根轨迹的基本概念、根轨迹与系统性能之间的关系，并从闭环零、极点与开环零、极点之间的关系推导出根轨迹方程，把他转化为常用的相角条件和模值条件形式，最终会出根轨迹。对于不能采用一、二及系统近似的高阶系统来说，其动态性能指标的确定比较复杂，通过用闭环主导极点的概念对高阶系统进行近似分析，从而得到高阶系统动态性能指标的估算公式。由于非线性系统形式多样，一般难以求得非线性微分方程的解析解，只能采用工程上的近似方法，如相平面法、描述函数法、逆系统法等，本次课设中运用 描述函数法分析含有死区特性环节的非线性系统。这次课程设计是对平时学习的一个检验，我不仅学会了用理论来分析检验实践，也学会了用实践来验证理论，总的来说，通过这次课程设计我受益匪