三阶系统综合分析与设计

1、引言在控制工程中，三阶系统非常普遍，但是三阶系统属于高阶系统,其动态性能指标的确定是比较复杂，不能像二阶系统那样可以用特定的公式计算。因此，我们可以借助于MATLAB软件对高阶系统进行分析。在课程设计中，我们不仅要掌握用MATLAB绘制闭环系统根轨迹和和系统响应曲线，还要掌握BODE图和Nyquist曲线的绘制。以及在比较点与开环传递函数之间加一个非线性环节后用负倒描述函数和Nyquist曲线判断系统的稳定性。1 设计内容1.1 设计题目: 三阶系统综合分析与设计初始条件：某单位反馈系统结构图如图1所示： EMBED Visio.Drawing.6 图1-1 图1-21.2 设计任务要求完成的

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

3、us函数：1、首先根轨迹绘制需要明确的是根轨迹起于开环极点（包括无限极点），终于开环零点（包括无限零点）。根据系统开环传递函数 可知：系统根轨迹有3条分支，开环零点无限远，有三个开环极点，分别是0，-3，-6。2、实轴上的根轨迹是那些在其右侧开环实数零点和开环实数极点总数和为奇数的区间，所以该系统根轨迹在实轴上的区间为(-,-6、-3,0。3、该系统根轨迹有3条渐近线，求渐近线与实轴交点,得-3.0j,而渐近线与实轴正方向的夹角为，既有。则可绘制系统根轨迹渐近线。4、根轨迹与虚轴的交点；闭环系统的特征方程： 即 令代入上式，可得： 解得 5、根轨迹的分离点： 解得 MATLAB为绘制根轨迹编程

4、如下：num=1; den=1 9 18 0; syms=tf(num,den); rlocus(syms)绘制出的根轨迹如图2-1所示：图2-1 闭环根轨迹2.2求取当闭环系统的一个极点为-8时的k值闭环系统的特征方程为：令=0，将闭环极点代入方程式中，从而可以得到=80.所以当-8为闭环系统的一个极点时，等于80。2.3求主导极点阻尼比为0.7时的K值 系统的动态性能基本上由接近虚轴的闭环极点确定，这样的极点被称为主导极点。故主导极点定义为对整个时间相应过程中起到主要作用的闭环极点，只有既接近虚轴，又不十分接近闭环零点的闭环极点才能成为主导极点。 在全部闭环极点中，选取最靠近虚轴又不十分靠

5、近闭环零极点的一个或几个闭环极点作为主导极点，略去不十分接近原点的偶极子，以及比主导极点距虚轴远6倍以上的闭环零、极点。 该系统中只有无限远零点，极点分别是0、-3、-6，其中0及-3比-6更加靠近虚轴，可以作为主导极点，-6则可忽略，对系统的影响较小。 故系统的闭环传递函数可写为： 又典型二阶系统的闭环传递函数为：比较上式，可得 ， ，又系统主导极点阻尼比 ，可以求得 ， 2.4稳态误差及稳态误差系数的计算控制系统的稳态误差数值，与开环传递函数G(s)H(s)的结构和输入信号R(s)的形式密切相关，对于一个给定的稳定系统，当输入信号形式一定时，系统是否存在稳态误差就取决于开环传递函数描述的系

6、统结构。一般，分子阶次为m，分母阶次为n的开环传递函数课表示为：则可根据已知信号的输入形式，判断系统是否存在原理性稳态误差及稳态误差的大小，根据已知条件可知系统为型系统。 误差系数的求取1、若,R为输入阶跃信号的幅值，则R(s)=R/s,此时有稳态误差 静态位置误差系数=2、若,R为输入速度信号的幅值，则R(s)=R/ ,此时有稳态误差 静态速度误差系数 3、若,R为输入加速度信号的幅值。则，此时有稳态误差 静态加速度误差系数 4、则当输入信号为时，稳态误差的求法如下：首先将输入信号分解为三个信号的叠加：, , 。输入信号R(t)的稳态误差就是信号、的稳态误差的和。输入信号为时，稳态误差输入信

7、号为，时，稳态误差输入信号为时，稳态误差所以输入信号的稳态误差为。2.5 绘制单位阶跃响应曲线绘制单位阶跃响应利用的函数是step函数。MATLAB绘制单位阶跃响应曲线编程如下：num=4.59; den=1 9 18 4.59;syms=tf(num,den);step(tf(num,den)绘制出的单位阶跃响应曲线如图2-2所示：图2-2 单位阶跃响应曲线由图中可以看到单位阶跃响应最终趋向于1，所以稳态误差为0，所以该系统是稳定的。2.6 绘制BODE图和Nyquist曲线及求取幅值裕度和相角裕度2.6.1绘制BODE图Bode图又称对数频率特性曲线图，由对数幅频曲线和对数相频曲线组成，是

8、工程中广泛使用的一组曲线。绘制轨迹利用的函数是bode函数：1、将传递函数写成伯德标准型，确定开环传递系数和各转折频率。 系统的开环频率特性为 由伯德标准型容易看书，开环传递系数为，转折频率为，。2、 确定低频段：由传递函数可知该系统为1型系统，即有微分环节，所以绘制低频段，可过， 作一条斜率为-20dB/dec的斜线。3、 绘制开环对幅频特性的渐近线：将低频段延伸到第一个转折点频率处。因为第一个转折频率是惯性环节的转折频率，所以，开环对数频率特性的渐近线下降20dB/dec,再延伸到第二个转折频率处，因为也是惯性环节，所以再下降20dB/dec。4、 绘制相频特性：绘制各个环节的对数相频特性

9、曲线，然后逐点叠加。一般在一些特征点上进行叠加，如各个转折频率处。5、 修正对数幅频特性。MATLAB绘制BODE图编程如下：num=4.59; den=1 9 18 0;bode(num,den);grid on绘制出的BODE图如图2-3所示：图2-3 BODE图2.6.2绘制Nyquist曲线 Nyquist曲线图又称开环幅相频率特性曲线图，是利用Nyquist判别系统稳定性的依据。绘制Nyquist曲线利用的函数是nyquistMATLAB绘制Nyquist曲线编程如下：num=4.59; den=1 9 18 0; syms=tf(num,den);nyquist(syms)绘制出的

10、Nyquist曲线如图2-4所示： 图2-4 Nyquist曲线2.6.3 求取幅值裕度h和相角裕度 设为系统的截止频率则 定义相角裕度为 相角裕度的含义是：对于闭环稳定系统，如果系统开环相频特性再滞后度，则系统处于临界稳定状态。设为系统的穿越频率则 ；定义幅值裕度为 幅值裕度的含义是: 对于闭环稳定系统，如果系统开环幅频特性再增大h倍，则系统将处于临界稳定状态。MATLAB求取幅值裕度h及相角裕度编程如下： num4.59; den1 9 18 0; sysm=tf(num,den); margin(sysm);grid on得到幅值裕度h及相位裕度如下图25 图2-5 幅值裕度及相角裕度图

11、 h=-20lg|=-29.8 ; 2.7 非线性控制系统的分析 描述函数的求取 图2-6 加入非线性环节的反馈结构图根据系统结构图可知该非线性类型为死区特性，又其中代入描述函数公式 即得到描述函数为 系统稳定性的判定由描述函数可知负倒数函数为 对于线性环节， 可解得穿越频率负倒函数 则 即曲线此时在实轴上由-指向-0.5的直线，而此时奈氏曲线过负实轴交于-0.02833，没有包围曲线，根据非线性系统稳定判据，该系统稳定。 图2-6 奈氏曲线放大图 3 体会和总结通过这次自动控制原理课程设计，我学到了关于MATLAB方面的很多知识，MATLAB它有着强大的数据处理能力，处理速度快，精度高，它不

12、仅可以用来绘制曲线，而且可以用来帮助解方程，以及做仿真处理，帮助验证理论分析的真确性。这次我做的是一个三阶系统的综合分析和设计，主要完成的任务是绘制根轨迹绘制单位阶跃响应求稳态误差绘制BODE图和Nyquist曲线求幅值裕度和相角裕度以及加入非线性死区利用负倒描述函数和Nyquist曲线来判断系统的稳定性。这些问题在课堂上老师都讲过，但是做起课程设计时还是遇到了一点困难，对于加入非线性死区利用负倒描述函数和Nyquist曲线来判断系统的稳定性这方面我比较陌生。但是通过翻阅资料，认真分析，我最终还是解决了问题。这次课程设计是对平时学习的一个检验，我不仅学会了用理论来分析检验实践，也学会了用实践来验证理论，总的来说，通过这次课程设计我受益匪