典型3阶系统的二阶参考

1、编 号：_审定成绩：_ 课程设计报告课程设计题目:_ 典型3阶系统的二阶参考_ _ _ 模型设计及仿真研究_ 单 位系别：_ 自动化_ 学 生 姓 名：_ 李春斌_ 专 业：_电气工程及其自动化_ _ 班 级：_1班_ 学 号：_0511110110_ 指 导 教 师：_汪纪锋_ 填表时间： 2021 年5 月重庆邮电大学移通学院教务处制重庆邮电大学移通学院毕业设计(论文)任务书设计(论文) 题目28 典型3阶系统的二阶参考模型设计及仿真研究 学生姓名 李春斌 系别 自动化系 专业 电气工程及其自动化 班级 05111101指导教师 汪纪锋 职称 教授 联系 42871150 教师单位 自动化

2、系 主 要 研 究 内 容 、 方 法 和 要 求研究典型型3阶系统动静态性能特性并完成综合设计，以到达使该系统满足工程实际性能指标的要求。根本要求：1.运用经典控制理论中的频域理论方法，分析给定的典型系统根本特性；2.根据提出的系统实际工艺指标，运用二阶参考模型法综合设计一套满足系统性能指标闭环系统方案；3.应用MATLAB/SIMULINK(或物理模拟)对设计系统进行仿真验证，给出仿真结果图。 进 度 计 划第4周第6周：完成系统建模，性能分析并提出校正设计方案；第4周第9周：按性能指标进行系统的综合设计；第10周第14周：完成原系统及校正后系统仿真研究，并作比拟研究；第15周第17周：撰

3、写论文、修改论文，完成辩论。 主 要 参 考 文 献1.汪纪锋、党晓圆，?现代控制理论?，人民邮电出版社，2021.11；2.郑大钟，?线性系统理论?第2版，清华大学出版社，2021.04；3.Chi-Tsong Chen,?Linear System Theory and Design?,HOLT.RINEHAET AND WINSTON,2001.03；4.熊晓君?自动控制原理实验教程?硬件模拟与MATLAB仿真，机械工业出版社，2021.01；5.方水良?现代控制理论及其MATLAB实践?，浙江大学出版社，2021.6；6.相关学术期刊等。指导教师签字： 汪纪锋 年 月 日教研室主任签字

4、： 年 月 日备注：此任务书由指导教师填写，并于毕业设计(论文)开始前下达给学生。重庆邮电大学移通学院毕业设计任务书简明技术资料一、设计题目：题目28 典型3阶系统的二阶参考模型设计及仿真研究二、系统说明：设三阶系统开环结构如下G0sysus=1ss+12s+1三、系统参量：系统输入信号：ut；系统输出信号：yt；四、设计指标：1.设定：在输入为r(t) = u(t) = a + bt,(其中：a=4 b=1/sec)2.在保证稳态误差essv0.2的前提下，要求动态期望指标：%5%；ts5s五、设计要求：1.分析原系统性能；2.试用频率特性法按二阶参考模型设计满足系统性能指标的闭环系统；3.

5、绘制系统校正前后的物理模拟仿真图。重庆邮电大学移通学院 自动化系指导教师：汪纪锋2021.12摘 要目录前言第一章 控制系统的简介第一节 自动控制系统的性能指标性能指标，是在分析一个控制系统的时候，评价系统性能好坏的标准。对于一个控制系统首要的要求是系统的绝对稳定性。在实际控制系统中，往往由于具体的对象不同或控制任务的指标不同，而对控制系统性能指标的要求有所不同。话虽如此，但仍可以对各种控制系统的性能指标概括为答题三个方面，即，快快速性、稳稳定性、准准确性。具体而言，对于定值控制系统，要能够迅速克服外界干扰的影响，使被控对象能够准确地恢复至给定值。系统性能的描述，可以分为动态性能和稳态性能。被

6、控量不随时间变化的平衡状态称为“稳态；而被控量随时间变化的不平衡状态称为“动态。在系统的动态过程中，其被控量是不断变化的，这一随时间变化的过程称为动态过程，也叫做控制过程。在此次设计中，主要研究三阶定常线性系统。其性能指标包括动态和静态指标两局部。静态指标即误差系数，动态性能指标包括上升时间tr、峰值时间tp、超调量Mp、调节时间ts等。第二节 系统概述对于原系统物理仿真结构首先第一个环节将是一个比拟环节，将输入信号与由输出反应回的信号进行比拟，以改善放大器的静态和动态性能；第二个是比例环节，对输出信号进行一定量的放大；第三个是积分环节，可以是系统的跟踪能力增强，积分环节是当输入信号为零时，输

7、出信号才能保持不变，而且能保持在任何位置上。在控制系统中，引用积分环节可以消除被控量的偏差。第四个是惯性环节，由于惯性环节系统的阻力，一开始输出并不与输入同步按比例变化，直到过渡过程结束,输出才能与输出保持比例，从而保证了控制过程作无差控制。第五个是惯性比例环节。第六个是反应环节，根据输入与输出在广义上是否相等来调节系统使之误差减小。第三节 设计根本要求基于频率特性法，按二阶参考模型法即=0.707设计满足系统性能指标的闭环系统。第四节 系统根本指标1.设定：在输入为rt=ut=a+bt,(其中：a=4 b=1/sec)2.在保证稳态误差essv0.2的前提下，要求动态期望指标：%5%；ts5

8、s第二章 系统建模第一节 各环节模型建立1.比例环节系统结构图如下列图1，结构框图如图2传递函数:G1=-RfR0=-K1 图1 图22.积分环节系统结构图如下列图3，结构框图如图4传递函数: G2s=-1R0C0S 图3 图43.惯性环节系统结构图如下列图5，结构框图如图6传递函数： G3s=-R1R01+R1C1S 图5 图64.反相器系统结构图如下列图7，结构框图如图8传递函数： G4s=-R2R0=-1 图7 图8第二节 系统数学模型三阶系统结构如下：G0sysus=1ss+12s+1开环传递函数为：G0s=1ss+12s+1闭环传递函数为：Gs=G0s1+G0s=1ss+12s+1+

9、1特征方程为：Ds= ss+12s+1+1=0第三章 系统分析第一节 二阶系统的数学模型二阶系统的开环传递函数为：G0s=2ss+2n闭环传递函数为：Gcs=n2s2+2ns+n2闭环传递函数的分母多项式等于零的代数方程式称为二阶系统的闭环特征方程，即：s2+2ns+n2=0闭环特征方程的两个根称为二阶系统的特征根，即s1,2=-n±n2-1上述二阶系统的数学模型中有两个特征参数 和 n，其中 称为二阶系统的阻尼比， n称为二阶系统的无阻尼振荡频率rad/s。二阶系统的系统分析和性能描述，根本上是以这两个特征参数来表示的。上述二阶系统的特征根表达式中，随着阻尼比 的不同取值，特征根s

10、i有不同类型的值，或者说特征根si在s平面上位于不同的位置，共有以下五种情况。>1时，特征根为一对不相等的负实根，位于s平面的负实轴上，使得系统的响应表现为过阻尼的。=1时，特征根为一对相等的负实根，也是位于s平面的负实轴上，系统的响应表现为临界阻尼的。0<<1时，特征根为一对带有负实部的共轭复数根，位于s平面的左半平面上，使得系统的响应表现为欠阻尼的。=0时，特征根为一对纯虚根，位于s平面的虚轴上，系统的响应表现为无阻尼的。<0时，特征根位于s平面的右半平面上，系统的响应是发散的。阻尼比取不同值时其特征根在s平面上的不同位置如下图图3.1.1 特征根在s平面上的不同位

11、置图阻尼比的大小决定了闭环极点在s平面的位置，反映了根的性质，极点的是不大小，决定了指数衰减的快慢，极点的虚部大小决定了系统响应震荡的快慢。第二节 稳定性分析频率稳定性判据有代数稳定性判据和频域稳定性判据又称Nyquist稳定性判据，简称奈氏判据。两种稳定性判别方法不同的是，代数稳定性判据是基于控制系统的闭环特征方程的判别方法。根本上提供的是控制系统绝对稳定性的信息，而对于系统的相对稳定性信息提供较少。而频域稳定性判据所依据的是控制系统的开环频率特性，也就是仅仅利用系统的开环信息，不仅可以确定系统的绝对稳定性，而且还可以提供相对稳定性的信息。也就是说，系统如果是稳定的，那么动态性能是否好；或者

12、如果系统是不稳定的，那么与稳定情况还差多少等。所以频域稳定性判据不仅用于系统稳定性分析，而且还可以方便的用于控制系统的设计与综合。1Nyquist判据负反应系统稳定的充分必要条件是：系统开环传递函数在G(s)H(s)平面上，Nyquist围线的象曲线逆时针绕(-1,j0)点的圈数R与G(s)H(s)在右半平面极点的个数P相同。即：系统在右半s闭环极点个数 Z = P R = 0由于G(s)H(s)曲线的对称性，因此可以用系统的开环频率特性曲线G(jw)H(jw)对(-1,j0)的包围情况来判断。设特性曲线G(jw)H(jw)对(-1,j0)的逆时针包围次数为N那么R=2N注意补充积分环节Nyq

13、uist围线上小1/4圆的象也可用G(jw)H(jw)曲线对(-, -1)实轴段的穿越计算NN+ 正穿越(由上到下) N- 负穿越(由下到上)闭合曲线F包围原点圈数的计算根据GH包围-1,j0的圈数，计算2劳斯Routh判据线性定常系统的特征方程为Ds=ansn+an-1sn-1+a1s+a0=0首先，作劳斯表，将方程的个系数间隔填入前两行，表如下所示“依照下式计算其余各项：， ， 将计算各项依照上述法那么全部计算完毕，填入劳斯表。计算完毕的劳斯表呈上三角形，系统稳定的充分必要条件为：劳斯表中，如果第一列元素全部大于零，系统就是稳定的；否那么系统是不稳定的。系统的特征方程为：Ds= ss+12

14、s+1+1=0作劳斯表：s3 2 1s2 3 1s1 13s0 1根据劳斯表可知第一列中全部大于零，所以原系统是稳定的。第三节 稳态误差分析不同系统在不同输入信号下，所产生的稳态误差也是不同形式的。第一种情况：输入信号为单位阶跃信号。10型系统 前向通路积分环节的个数为零，即=0，单位阶跃信号输入时，由于0型系统的静态位置误差系数Kp等于开环增益的大小K0，所以0型系统在阶跃信号输入作用时的稳态误差也为常数。2型系统 前向通路积分环节的个数为1，即=1，单位阶跃信号输入时，由于型系统的静态位置误差系数Kp等于无穷大，所以型系统得稳态误差为零，型系统是一阶无差系统。第二种情况：输入信号为单位斜坡

15、信号。10型系统 前向通路积分环节的个数为零，即=0，0型系统施加斜坡信号，当时间趋于无穷大时，其稳态误差的值也是无穷大。也就是说，0型系统不能跟踪等速率信号。2型系统 前向通路积分环节的个数为1，即=1，型系统施加斜坡信号，当时间趋于无穷大时，其稳态误差趋于常数值，且大小等于系统的开环增益K0的倒数。也就是说，型系统有跟踪速率信号的能力，但是在跟踪过程中，只能实现有差跟踪。可以加大开环增益K0来减小稳态误差，可是不能消除它。3型系统 前向通路积分环节的个数为1，即=2，如果系统的前向通路中有两个积分环节，那么在跟踪等速率信号时，由于稳态误差为零，所以可以实现无差跟踪。也就是说，只要系统是稳定

16、系统，那么系统的响应在过了暂态时间之后，就与等速率信号相同了。所以=2得系统又称为二阶无差系统。第三种情况：输入信号为加速度信号。当输入信号为加速度信号时，0型系统、型系统的稳态误差都是无穷大，只有型系统的稳态误差等于常数，且与系统开环增益的大小成反比，有限的实现对于加速度信号的跟踪。上面的分析汇成总表如下：系统静态误差系数稳态误差lss类型KpKvKa1ttt220型K00011+Kp型K0001Kv型K0001Ka表3.3.1 稳态误差分析简表由系统的开环传递函数得知，此系统是型3阶系统。用MATLAB进行系统稳定性判定系统的闭环传递函数为：Gs=G0s1+G0s=1ss+12s+1+1&

17、gt;> num=1;>> den=2 3 1 1;>> G=tf(num,den) Transfer function: 1-2 s3 + 3 s2 + s + 1 >> p=eig(G)p = -1.3982 -0.0509 + 0.5958i -0.0509 - 0.5958i>> p1=pole(G)p1 = -1.3982 -0.0509 + 0.5958i -0.0509 - 0.5958i>> r=roots(den)r = -1.3982 -0.0509 + 0.5958i -0.0509 - 0.5958i分析

18、：系统的特征根全部位于s平面的左半平面，很显然该系统是稳定的。第四节 动态性能分析系统的阶跃响应一般曲线如下图。图3.4.1 二阶欠阻尼系统的单位阶跃响应根据图中所展示的响应特性，可以定义如下的性能指标。1上升时间tr上升时间tr指阶跃响应ct上升至稳态值所需要的时间。考虑到不敏感区或者允许误差，有时取为从稳态值的10%上升至90%时所需要的时间。2峰值时间tp峰值时间tp指阶跃响应ct从运动开始至到达第一峰的时间。3超调量Mp超调量Mp指系统响应的第一峰值超出稳态值的局部，将其取百分比，可以表示为Mp=ctp-cc×100%=ctp-11×100%=-1-2×1

19、00%4调节时间ts调节时间ts指阶跃响应ct到达稳定值的时间。评价系统整个过程的响应速度，是速度和阻尼程度的综合指标。5稳态误差lss稳态误差是当时间t趋于无穷时，系统希望的输出与实际的输出之差，表示为lss。从上述系统阶跃响应性能指标可以看出，各个时间指标反映了系统的快速性。其中，上升时间tr、峰值时间tp反映了系统的初始快速性；而调节时间ts反映了系统的总体快速性；超调量Mp描述了系统的平稳性；稳态误差lss描述了系统的准确性。3.4.1动态平稳分析这里指超调量Mp是否满足要求。特征方程为：Ds= ss+12s+1+1=0那么有：Ds=2s3+3s2+s+1=Ds为闭环特征方程为阻尼比n

20、为系统频率3.4.2动态快速性分析3.4.3原系统的图像分析在频率分析法中，以控制系统的频率特性作为数学模型，以波德图或者其他图表作为分析工具，来研究。分析控制系统的动态性能与温蔼性能。频率特性曲线有两种种表示方法：极坐标图、对数坐标图。1极坐标图极坐标图又称为幅相图、奈奎斯特图Nyquist.Nyquist曲线。频率特性Gj是输入信号频率w的复变函数，用实坐标与虚坐标可表示为Gj=ReGj+ImGj，当频率从-变到+时， Gj在由实轴与虚轴构成的复平面上走过的根轨迹就称为Gj的极坐标图。MATLAB在绘制奈奎斯特曲线时，频率是从连续变化的，而在自动控制原理的教材中，一般只绘制频率从局部的曲线

21、。可以分析得出，曲线在范围与内，是以横轴为镜像的。极坐标图常用于频域稳定性分析的作图中。2对数坐标图对数坐标图又称为波德图。频率特性的矢量表达式为Gj=GjargGj=A其中A称为频率特性的幅频特性，是的函数。当由0+时，它展示了Gj幅值的变化规律。称为频率特性的相频特性，也是的函数。当由0+时，它展示了Gj相位角的变化规律。对数幅频特性L和对视相频特性两条曲线统称对数频率特性，又称波德图。原系统的波德图如下：系统开环结构如下G0sysus=1ss+12s+1使用MATLAB编写程序观察原系统的频率特性及阶跃响应。程序如下：>> s=tf('s');>>

22、 G=1/(s*(s+1)*(2*s+1) Transfer function: 1-2 s3 + 3 s2 + s >> bode(G)>> nyquist(G)程序运行结果得到系统Bode图和Nyquist曲线，分别如下列图 3.4.2和3.4.3所示。图 3.4.2 校正前的系统Bode图图3.4.3 校正前的系统Nyquist曲线由原系统波德图分析可知：原系统的相频特性曲线一直位于-270之上，所以稳定裕度c>0，且幅频特性曲线Lg>0，显然系统是稳定的。由原系统的Nyquist曲线分析可知：原系统不存在右极点，即P=0，且在-1，j0点，正穿越和负

23、穿越个数为0，所以系统是稳定的。3.3.4截止频率的计算c令对数幅频特性L过0dB的频率为c，那么有20lg2lg1-lgc=-60解得c=1.23.3.5相位裕度计算c定义相位裕度c，计算公式如下c=180+cc=180+1.2c=40°第四章 系统的综合设计第一节 超前校正网络具有超前特性，电路图如下图。图4.1.1超前校正网络传递函数为GDs=1+s1+s其中T=R1C，=R2R1+R2<1频率特性为GDj=1+jT1+jT极坐标图的起点为A0=0=0终点为A=1=0可以证明：1：0>0，即相位角总是超前的；2轨迹为上半圆。使用相位超前校正可以改善系统的动态品质，减

24、小超调量，缩短响应时间。超前校正的目的是改善系统的动态性能，实现在系统静态性能不受损的前提下，提高系统的动态性能。通过参加超前校正环节，利用其相位超前特性来增大系统的相位裕度，改变系统的开环频率特性。一般使校正环节的最大相位超前角出现在新的穿越频率点。第二节 滞后校正网络具有相位滞后特性，滞后校正网络电路图如下图。图4.2.1 滞后校正网络传递函数为G1s=1+Ts1+Ts其中T=R2C，=R1+R2R2>1那么频率特性为G1j=1+jT1+jT极坐标起点为A0=10=0终点为A=1=0可以证明：1：0<0，即相位总是滞后的；2轨迹为下半圆。在滞后校正中，我们利用的是滞后网络在高频

25、段的衰减特性，而不是相位滞后特性。滞后校正通过参加滞后校正环节，使系统的开环增益有较大幅度增加，同时又使校正后的系统动态指标保持原系统的良好状态。他利用滞后校正环节的低通滤波特性，在不影响校正后系统低频特性的情况下，使校正后系统中高频段增益降低，从而使其穿越频率前移，到达增加系统相位裕度的目的。第三节 滞后超前校正图4.3.1 超前滞后校正（1） 在0<<0时，校正网络起滞后网络的作用；在0<<时，它起超前网络的作用；（2） 网络的超前局部增加了相位超前角，增加了系统的相位裕量，改善了系统的动特性；（3） 网络的滞后局部，在剪切频率附近和以上，将产生衰减，因而它允许在低

26、频段增加增益，故改善了系统的稳态特性;（4） 因而系统在动态和稳态均不能满足要求时,可以采用滞后-超前校正。（5） 滞后-超前校正的根本原理是利用校正装置的超前局部增大系统的相位裕度，同时利用其滞后局部来改善系统的稳态性能。滞后超前校正适用于对于校正后系统的动态和静态性能有更多更高要求的场合。施加滞后超前校正环节，主要是利用其超前局部增大系统的相位裕度，以改善系统的动态性能；利用其滞后局部改善系统的静态性能。以上3种不同的校正方法的一般性设计步骤如下：1) 根据静态性能指标，计算开环系统的增益。之后求取校正前后系统的频率特性指标。并与设计要求进行比拟。2) 确定校正后期望的穿越频率，具体值的选

27、取与所选择的校正方式相适应。3) 根据待设计的校正环节的形式和转折频率，计算相关参数，进而确定校正环节。4) 得出校正后系统。检验系统是否满足设计要求。如果不满足那么从第二步从新开始。第四节 校正方案确实定由于调节时间，超调量都不满足要求，所以可以使用相位超前校正方案。第五节 设计校正量与校正装置4.5.1计算最大相位超前角mm=c-c0+520=4.5.2计算低频衰减=1-sinm1+sinm=4.5.3确定新的开环截止频率cL0=-1220log101=4.5.4确定截止频率1=1T=c=2=1T=1c=4.5.6补偿增益KcKc=1=4.5.7校正后的开环频率特性4.5.8校正后的截止频

28、率和稳定裕度4.5.9校正后的系统波德图和Nyquist曲线图第五章 物理模拟仿真第一节 原系统物理仿真物理仿真，又称物理效应仿真，是指仿真的过程以物理性质和几何形状相似为根底，其他性质不变的仿真。控制系统是由各种元件或部件按照一定的要求组合在一起的。从形式结构上看，有各种各样的原件，但从动态性能或数字模型上看，却可以分为几个根本环节，这样对于我们研究系统和分析系统带来很大的便利，便于理解各个部件对于各个环节或是系统的性能影响。第一章已经对于各个环节进行了讲解，由于各环节之间没有负载效应，故各环节之间可以直接连接，并且个环节之间的放大电路可以合并为一个放大电路，从而可得出原系统物理模拟结构图如

29、下列图所示；经过校正后的系统物理模拟结构图如下列图所示：第二节 MATLAB仿真5.2.1MATLAB的简介Matlab是起源于美国MathWorks公司发布主要面向数值计算、科学数据可视化以及交互式程序设计的高技术计算语言。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境之中，为科学研究、为适应众多专业科技人员的需要；MathWorks同时提供了数十个应用工具箱为科学和工程领域各类特殊问题及应用定制MATLAB运行环境；并为全面解决复杂数值计算问题以及CAD研究等提供了综合解决方案。MATLAB仿真在科学研究中的地位越来越高，如何

30、利用MATLAB仿真出理想的结果，关键在于如何准确的选择MATLAB的仿真。Matlab的应用范围非常广，包括信号和图像处理、通讯、控制系统设计、调试和测量、财务建模和分析以及计算机生物等众多应用领域。附加的工具箱单独提供的专用Matlab函数集扩展了Matlab环境，以解决这些应用领域内特定类型的问题。Matlab可以用来进行以下各种工作：1. 数值分析2. 数值和符号计算3. 工程与科学绘图4. 控制系统的设计与仿真5. 数字图象处理6. 数字信号处理7. 通讯系统设计与仿真8. 财务与金融工程5.2.2MATLAB/SIMULINKSimulink是面向系统信号流图进行组态的仿真软件，通过Simulink