三阶系统综合分析与设计

课程设计任务书

学生姓名： 专业班级： 指导教师： 工作单位： 自动化学院 题 目: 三阶系统综合分析与设计

初始条件：某单位反馈系统结构图如图1所示：

R(s)-Ks(s?2)(s?4)C(s)R(s)-Me025KC(s)ss(s82))(s?4)(s?? 图1 图2

要求完成的主要任务: （包括课程设计工作量及其技术要求，以及说明书撰写等具体

要求）

1、 试绘制随根轨迹

2、 当-6为闭环系统的一个极点时，K=?

3、 求取主导极点阻尼比为0.7时的K值（以下取这个值）

4、 分别求取位置误差系数、速度误差系数、加速度误差系数及输入信号为

r(t)?1(t)?2t?t单位阶跃信号、斜坡信号及单位加速度信号时的稳态误差

25、 用Matlab绘制单位阶跃相应曲线

6、 绘制Bode图和Nyquist曲线，求取幅值裕度和相角裕度

7、 如在比较点与开环传递函数之间加1个非线性环节，如图2所示，其中

e0?0.7,M?1.7，试求取非线性环节的描述函数，并根据负倒描述函数和

Nyquist图判断系统的稳定性 8、 认真撰写课程设计报告。

时间安排：

任务 审题、查阅相关资料 分析、计算 编写程序 撰写报告 论文答辩 时间（天） 1 1.5 1 1 0.5 指导教师签名： 年 月 日 系主任（或责任教师）签名： 年 月 日

武汉理工大学《自动控制原理》课程设计说明书

目录

摘要 ...................................................................... 1 1 设计意义及要求 ......................................................... 2 1.1 设计意义 ............................................................ 2 1.2 设计要求 ............................................................ 2 2 设计过程 ............................................................... 3 2.1 绘制根轨迹 ......................................................... 3 2.1.1 理论计算 ........................................................ 3 2.1.2 MATLAB绘制根轨迹 ................................................ 4 2.2 极点-6时的K值的求取 ................................................ 4 2.3主导极点阻尼比为0.7时的K值求取 .................................... 5 2.4稳态误差 ............................................................ 6 2.4.1 系统的误差系数分析 .............................................. 6 2.4.2 系统的稳态误差分析 .............................................. 6 2.5 用MATLAB绘制单位阶跃响应曲线 ........................................ 6 2.6绘制BODE图和NYQUIST曲线，求取幅值裕度和相角裕度 ..................... 7 2.6.1 绘制Bode图 ..................................................... 7 2.6.2绘制Nyquist曲线 ................................................. 8 2.6.3 幅值裕度和相角裕度 .............................................. 9 2.7 系统加入非线性环节的稳定性分析 ...................................... 9 2.7.1 非线性环节的描述函数的求取 ..................................... 10 2.7.2 负倒描述函数的求取 ............................................. 11 2.7.3 系统稳定性的判据及原理 ......................................... 12 2.7.4 系统稳定性判断 ................................................. 13 结束语 ................................................................... 14 参考文献 ................................................................. 15

摘要

《自动控制原理》是为了培养学生统筹运用自动控制原理课程中所学的理论知识，掌握反馈控制系统的基本理论和基本方法，对工程实际系统进行完整的全面的分析和综合而开设的重要教学环节，此次课程设计可以锻炼学生的动手能力和解决问题的原因，把课本知识运用的实际中，同时也以更为自主创新的形式检验了学生对所学知识的掌握程度。

三阶系统是以三级微分方程为运动方程的控制系统。在控制工程中，三阶系统非常普遍，其动态性能指标的确定是比较复杂。在工程上常用闭环主导极点的概念对三级系统进行分析，或直接用MATLAB软件进行高级系统分析。在课程设计中，我们不仅要掌握用MATLAB绘制闭环系统根轨迹和和系统响应曲线，用系统的闭环主导极点来估算三阶系统的动态性能，还要掌握BODE图和Nyquist曲线的绘制。以及在比较点与开环传递函数之间加一个非线性环节后用负倒描述函数和Nyquist曲线判断系统的稳定性。

关键字： 三阶系统 闭环主导极点 MATLAB

1 设计意义及要求 1.1 设计意义

本次设计主要是让学生将自动控制原理中所学的理论知识与实践结合起来，对工程实际系统进行完整全面分析和综合，掌握利用MATLAB对控制理论进行分析，研究和仿真技能，提高分析问题和解决问题的能力。本次的课程设计是对我们平时学习的理论知识的一个检验，也是让我们更加熟练的运用MATLAB软件，更好的解决自动控制方面的一些问题。

1.2 设计要求

初始条件：某单位反馈系统结构图如下图所示：

R(s)-Ks(s?2)(s?4)C(s)R(s)-Me025KC(s)ss(s82))(s?4)(s?? 图1 图2

要求完成的主要任务: （包括课程设计工作量及其技术要求，以及说明书撰写等具体

要求）

1、 试绘制随根轨迹

2、 当-6为闭环系统的一个极点时，K=?

3、 求取主导极点阻尼比为0.7时的K值（以下取这个值）

4、 分别求取位置误差系数、速度误差系数、加速度误差系数及输入信号为

r(t)?1(t)?2t?t2单位阶跃信号、斜坡信号及单位加速度信号时的稳态

误差

5、 用Matlab绘制单位阶跃相应曲线

6、 绘制Bode图和Nyquist曲线，求取幅值裕度和相角裕度

7、 如在比较点与开环传递函数之间加1个非线性环节，如图2所示，其中

e0?0.7,M?1.7，试求取非线性环节的描述函数，并根据负倒描述函数

和Nyquist图判断系统的稳定性

8、 认真撰写课程设计报告。

2 设计过程

2.1 绘制根轨迹

某单位反馈系统结构图如下图所示：

R(s)-Ks(s?2)(s?4)C(s)

2.1.1 理论计算

（1）根轨迹的起点和终点。根轨迹起于开环极点(包括无限极点)，终于开环零点（包括无限零点）。根据系统开环传递函数可得：系统的开环极点分别为0、-2、-4，开环零点为无限远。

（2） 根轨迹的分支数。n=3，m=0，所以分支数为3 。且它们是连续的并且对称于实轴。

（3） 根轨迹在实轴上的分布。实轴上的某一个区域，若其右边开环零、极点的个数之和为奇数，则该区域必是根轨迹。因此实轴上(-?，-4]、 [-2,0]必为根轨迹。

nm?（4）根轨迹的渐进线。渐近线与实轴的交点?a?i?1pi??zj?1jn?m = = -2，即

3?0?2?4与实轴的交点为（-2,j0）,与实轴的交角为?a?（k=0，1，2）。

（5）确定根轨迹的分离点。 分离点的方程为:离点d=-0.85,d=-3.15（不合题意，舍去）

1d?1d?2?2k?1??n?m?2k?1??3??3，?,5?3

?1d?4?0 ，因此可以求得分

（6）根轨迹与虚轴的交点。 由开环传递函数写出系统的闭环特征方程式为 s3?6s2?8s?K=0

将s?j?代入上式，可得实部方程为?6?2?k?0，虚部方程为?j?3?8j??0 解得 ???22，k?48 。 所以与虚轴的交点为?0,j22?,?0,?j22?。

2.1.2 MATLAB绘制根轨迹

MATLAB为绘制根轨迹编程如下： num=[48]; den=[1 6 8 0]; rlocus(num,den);

绘制出的根轨迹如图1所示：

Root Locus1086420-2-4-6-8-10-14Imaginary Axis (seconds-1)-12-10-8-6-4-2024Real Axis (seconds-1)

图1 闭环根轨迹

根据理论计算和MATLAB绘制的根轨迹示意图，可以知道理论计算和MATLAB绘制的根轨迹完全相符。从而可以知道所绘制的根轨迹是正确的。

2.2 极点-6时的K值的求取

闭环系统的特征方程为：

s?6s?8s?K?0

32将闭环极点s=-6代入方程式中，从而可以得到k?48。所以当-6为闭环系统的一个极点时，k等于48。

2.3主导极点阻尼比为0.7时的K值求取

当主导极点阻尼比为0.7时，先做出??0.7的等阻尼比线，使这条直线与负实轴方向的夹角为??cos?1??cos?10.7=45.6?,此直线的斜率为k=-tan45.6?=-1.02,在MATLAB中画出此直线，并能找到与根轨迹的交点s1即是满足??0.7的闭环主导极点之一。编写的此程序为：

k=-1.02; x=-25:5; y=k*x; plot(x,y); hold on num=[48];

den=[1 6 8 0]; rlocus(num,den);

Root Locus3020Imaginary Axis (seconds-1)100-10-20-30-25-20-15-10Real Axis (seconds-1)-505

图2

由图2我们可以得出s1=-0.753+j0.768，由根轨迹的对称性，可求得另一个极点为

s2=-0.753-j0.768，由幅值条件可知，闭环极点s1对应的根轨迹的增益为

K= |s1|| s1+2|| s1+4|=|-0.753+j0.768||1.247+j0.768||3.247+j0.768|=5.27

经验证，s1和s2满足主导极点的条件，另一极点实部的模比主导极点实部的模大三倍以上，不是主导极点，所以该系统可近似成一个由主导极点构成的二阶系统传递函数为：

?（s）??n222s?2??ns??ns1s21.16??2（s?s1）（s?s2）s?1.51s?1.16对应的系统的开环增益为Kr= K/8=0.66。

2.4稳态误差

2.4.1 系统的误差系数分析

位置误差系数

Kp?limG(s)H(s)?lim??；s(s?2)(s?4）

5.27 速度误差系数

s?0s?0Kv?limsG(s)H(s)?lims 加速度误差系数

s?0s?0?0.66；s(s?2)(s?4）

5.27

Ka?limsG(s)H(s)?lims22

s?0s?0?0；s(s?2)(s?4）

5.272.4.2 系统的稳态误差分析

当输入为r1（t）=1（t）时，稳态误差为ess1= 当输入为r2（t）=2t时，稳态误差为ess2= 当输入为r3（t）=t2时，稳态误差为ess3= 2Kv211?K?2 ?p11???0;Ka?3.03;0.662???;0

则当输入为r(t)= r1（t）+ r2（t）+ r3（t）时，总的稳态误差为 ess= ess1+ ess2+ ess3=?.

2.5 用Matlab绘制单位阶跃响应曲线

系统的闭环传递函数为

?(s)?5.27s?6s?8s?5.2732 其分子系数为5.27，分母系数分别为1，6，8，5.27

利用MATLAB程序可编以下程序： num=[5.27]; den=[1 6 8 5.27]; step(num,den)

在MATLAB程序中输入此程序，运行后得到单位阶跃响应曲线如图3所示。

Step Response1.41.21Amplitude0.80.60.40.200123456789Time (seconds)

图3 单位阶跃响应曲线

2.6绘制Bode图和Nyquist曲线，求取幅值裕度和相角裕度

2.6.1 绘制Bode图

系统的开环传递函数为

G(s)?5.27s?6S?8S32 其分子系数为5.27，分母系数分别为1，6，8，0。 利用MATLAB程序可编以下程序： num=[5.27]; den=[1 6 8 0]; bode(num,den)

在MATLAB程序中输入此程序，运行后得到Bode图如图4所示。

Bode Diagram50Magnitude (dB)Phase (deg)0-50-100-150-90-135-180-225-27010-210-1100101102Frequency (rad/s)

图4 Bode图

2.6.2绘制Nyquist曲线

系统的开环传递函数为

G(s)?5.27s?6S?8S32 可得其分子系数其分子系数为5.27，分母系数分别为1，6，8，0。利用MATLAB程序可编以下程序： num= [5.27]; den=[1 6 8 0]; nyquist(num,den)

在MATLAB程序中输入此程序，运行后得到Nyquist如图5所示。

其中Nyquist图与实轴的交点：将s?j?代入开环传递函数中，令虚部等于0，得到

??2.83，此时曲线与负实轴的交点为（-0.11，j0）。

结束语

通过这次自动控制原理课程设计，我掌握并巩固了自动控制原理知识 ，这次我做的是一个三阶系统的综合分析和设计，主要完成的任务是绘制根轨迹﹑绘制单位阶跃响应﹑求稳态误差﹑绘制Bode图和Nyquist曲线﹑求幅值裕度和相角裕度以及加入非线性死区利用负倒描述函数和Nyquist曲线来判断系统的稳定性。这些问题在课堂上老师大部分都讲过，但是做起课程设计时还是遇到了一点困难，对于加入非线性死区利用负倒描述函数和Nyquist曲线来判断系统的稳定性这方面我比较陌生。但是通过翻阅资料，认真分析，我最终还是解决了问题。

在这次设计中，我不仅学会了用理论来分析检验实践，学到了关于MATLAB方面的很多知识，MATLAB它有着强大的数据处理能力，处理速度快，精度高，它不仅可以用来绘制曲线，而且可以用来帮助解方程，以及做仿真处理，帮助验证理论分析的真确性。

这次课程设计是对平时学习的一个检验，也学会了用实践来验证理论，总的来说，通过这次课程设计我受益匪浅。

参考文献

[1]胡寿松.《自动控制原理（第五版）》.科学出版社 [2]王万良.《自动控制原理》.高等教育出版社

[3]胡寿松.《自动控制原理习题解析（第五版）》.科学出版社

[4]薛定宇.《反馈控制系统设计与分析—MATLAB语言应用》.清华大学出版社 [5]张静.《MATLAB在控制系统中的应用》，电子工业出版社，2007

本文来源：https://www.bwwdw.com/article/kjvg.html