用MATLAB进行控制系统的超前校正设计

1、武汉理工大学自动控制原理课程设计说明书目录摘要11课程设计目的22设计条件及任务要求22.1设计条件22.2设计任务要求23设计基本原理33.1超前校正33.2根轨迹法54设计过程64.1基本思路及步骤64.2校正前系统分析64.2.1开环增益64.2.2相位裕度和幅值裕度74.2.3伯德图84.2.4根轨迹94.3超前校正系统设计104.3.1 理论分析104.3.2参数计算114.3.3编程设计124.4校正后系统分析134.4.1伯德图134.4.2根轨迹145结果对比与分析145.1校正前后阶跃响应曲线145.2结果分析166总结16参考文献17附录18摘要在自动控制理论中，超前校正是

2、相当重要的一环，对于系统的优化有很重要的意义。用MATLAB进行控制系统的超前校正设计是对所学的自动控制原理的初步运用。本课程设计先针对校正前系统的稳定性能，用MATLAB画出其根轨迹及伯德图进行分析，是否达到系统的要求，然后对校正装置进行参数的计算和选择，串联适当的超前校正装置。最后通过阶跃响应曲线对比前后系统，论证校正后的系统是否达到要求。关键字：超前校正 MATLAB 根轨迹 伯德图用MATLAB进行控制系统的超前校正设计1课程设计目的本次课程设计通过让学生自主设计控制系统，进一步加深对自动控制原理基本知识的掌握和认识，更深入地研究相关理论。并通过对系统的分析和设计，提高实际动手能力。使

3、用MATLAB软件，学习利用根轨迹及伯德图等对系统进行分析。旨在培养学生的系统分析能力、设计能力及编程能力，巩固理论所学。2设计条件及任务要求2.1设计条件初始条件如下：已知一单位反馈系统的开环传递函数是要求系统的静态速度误差系数，幅值裕度。2.2设计任务要求（1）用MATLAB作出满足初始条件的K值的系统伯德图，计算系统的幅值裕度和相位裕度。（2）在系统前向通路中插入一相位超前校正，确定校正网络的传递函数，并用MATLAB进行验证。（3）用MATLAB画出未校正和已校正系统的根轨迹。（4）课程设计说明书中要求写清楚计算分析的过程，列出MATLAB程序和MATLAB输出。说明书的格式按照教务处

4、标准书写。3设计基本原理3.1超前校正超前校正的基本原理就是利用超前相角补偿系统的滞后相角，改善系统的动态性能，增加相位裕度，提高系统的稳定性。超前校正利用超前网络或PD控制器进行串联校正，在前向通道中串联传递函数为：=，>1其中参数、为可调，如图3-1所示。图3-1超前校正装置结构图超前校正装置的伯德图如图3-2所示。图3- 2超前校正伯德图超前校正的目的是改善系统的动态性能，实现在系统静态性能不受损的前提下，提高系统的动态性能。通过加入超前校正环节，利用其相位超前特性来增大系统的相位裕度，改变系统的开环频率特性，提高其稳定性。一般使校正环节的最大相位超前角出现在系统新的穿越频率点。如

5、图3-3所示，典型的无源超前装置由阻容元件组成。装置的传递函数为 （1） 式中 图3-3典型无源超前装置 通常a为分度系数，T叫时间常数，由式（1）可知，采用无源超前网络进行串联校正时，整个系统的开环增益要下降a倍，因此需要提高放大器增益交易补偿。根据式（1），可以得无源超前网络的对数频率特性，超前网络对频率在1/aT至1/T之间的输入信号有明显的微分作用，在该频率范围内，输出信号相角比输入信号相角超前，超前网络的名称由此而得。在最大超前角频率处，具有最大超前角。超前网络（2-1）的相角为： （2）将上式对求导并且令其为零，得到最大超前角频率： 将上式代入（2），得最大超前角： 同时还容易得到

6、。最大超前角仅仅与衰减因子a有关，a值越大，超前网络的微分效果越强。但是a的最大值还受到超前网络物理结构的制约，通常情况下，a取为20左右，这也就意味着超前网络可以产生的最大相位超前约为65°，如果所需要的大于65°的相位超前角，那么就可以采用两个超前校正网络串联实现，并且在串联的两个网络之间加入隔离放大器，借以消除它们之间的负载效应。所以通过以上的分析发现，利用超前网络进行串联校正的基本原理，是利用超前网络的相角超前特性。只要正确的将超前网络的交接频率1/aT或1/T选在待校正系统截止频率的两旁，并适当的选择参数a和T，就可以使已校正系统的截止频率和相角裕度满足性能指标的

7、要求，从而改善系统的动态性能，使校正后的系统具有以下特点：（1）低频段的增益满足稳态精度的要求；（2）中频段对数幅频特性的斜率为-20db/dec，并且具有较宽频带，使系统具有满足的动态性能；（3）高频段要求幅值迅速衰减，以减少噪声的影响。3.2根轨迹法根轨迹法是一种求特征根的简单方法，在控制系统的分析与设计中得到广泛的应用。这一方法不直接求解特征方程，用作图的方法表示特征方程的根与系统某一参数的全部数值关系，当这一参数取特定值时，对应的特征根可在上述关系图中找到。这种方法叫根轨迹法。根轨迹法具有直观的特点，利用系统的根轨迹可以分析结构和参数已知的闭环系统的稳定性和瞬态响应特性，还可分析参数变

8、化对系统性能的影响。在设计线性控制系统时，可以根据对系统性能指标的要求确定可调整参数以及系统开环零极点的位置，即根轨迹法可以用于系统的分析与综合。在控制系统的分析和综合中，往往只需要知道根轨迹的粗略形状。由相角条件和幅值条件所导出的 8条规则，为粗略地绘制出根轨迹图提供方便的途径。根轨迹绘制方法如下：（1）根轨迹的分支数等于开环传递函数极点的个数。（2）根轨迹的始点（相应于K=0）为开环传递函数的极点，根轨迹的终点（相应于K=）为开环传递函数的有穷零点或无穷远零点。（3）根轨迹形状对称于坐标系的横轴（实轴）。（4）实轴上的根轨迹按下述方法确定：将开环传递函数的位于实轴上的极点和零点由右至左顺序

9、编号，由奇数点至偶数点间的线段为根轨迹。（5）实轴上两个开环极点或两个开环零点间的根轨迹段上，至少存在一个分离点或会合点，根轨迹将在这些点产生分岔。（6）在无穷远处根轨迹的走向可通过画出其渐近线来决定。渐近线的条数等于开环传递函数的极点数与零点数之差。（7）根轨迹沿始点的走向由出射角决定，根轨迹到达终点的走向由入射角决定。（8）根轨迹与虚轴（纵轴）的交点对分析系统的稳定性很重要，其位置和相应的K值可利用代数稳定判据来决定。4设计过程4.1基本思路及步骤超前校正的方法有很多种，本次课程设计中采用频域法步骤如下：（1）根据稳态误差要求，确定开环增益K。（2）利用已确定的开环增益，计算待校正系统的相

10、角裕度和幅值裕度。（3）根据截止频率的要求，计算a和T。令，以保证系统的响应速度，并充分利用网络的相角超前特性。显然成立的条件是 根据上式不难求出a值，然后确定T。（4）验算已校正系统的相角裕度。验算时，求得，再由已知的算出待校正系统在时的相角裕度。最后，按下式算出 如果验算结果不满足指标要求，要重选，一般使增大，然后重复以上步骤。4.2校正前系统分析4.2.1开环增益校正前系统开环传递函数为：题目要求系统的静态速度误差系数，幅值裕度。取Kv=20s-1，根据静态速度误差系数Kv的计算公式：Kv=lims0sGsHs得： Kv=lims0sKss+4=K4=20可知校正前的开环增益为：K=80

11、。则校正前的开环传递函数为：GS=20ss4+14.2.2相位裕度和幅值裕度使开环频率特性的相角为的频率，称为相位穿越频率，即 使开环频率特性的幅值为1或者的频率，称为增益穿越频率或者截止频率，即或者 相角裕度和幅值裕度或的定义如下： 或者 根据和可列出如下方程20c164+c2=1 tan-1g4=90°求得校正前系统截止频率c=8.51rad/s，相位穿越频率g=。于是得出校正前系统的相位裕度和幅值裕度为Pm=0=180°+GjcHjc=180°-90°-tan-1c4=25.2°Gm=Kg=1|GjgHjg|=+4.2.3伯德图由伯德图可

12、得校正前系统的相位裕度和幅值裕度。由校正前系统开环传递函数GS=20ss4+1用MATLAB绘制伯德图的程序如下： num=20;>> den=1/4,1,0; %开环传递函数的分子、分母多项式>> sys=tf(num,den); %系统建模>> mag,phase,w=bode(sys); >> margin(sys); %画伯德图并求相位穿越频率和截止频率>> grid由MATLAB绘制出的系统校正前的伯德图如图4-1所示。图4-1校正前系统伯德图由图可知，校正前系统：相位裕度Pm=25.2°截止频率c=8.51rad

13、/s 幅值裕度Gm=+穿越频率g=+用MATLAB得出的结果与计算结果一致。4.2.4根轨迹用MATLAB绘制根轨迹的程序如下：num=20;>> den=1/4,1,0;>> rlocus(num,den);>> grid图4-2校正前系统根轨迹右图为校正前系统的根轨迹。4.3超前校正系统设计4.3.1 理论分析先求取超前校正的最大超前相角及取得最大超前相角的频率。m=-0+超前校正时系统的截止频率c增大，未校正系统的相角一般是较大的负相角。这里加上一个附加角，用于补偿。超前校正的频率特性 相频特性 令，求得 于是 即当时，超前相角最大为。可以看出，只与有

14、关。由上面的结果可见，只与有关。越大，越大，对系统的相角补偿越大，超前网络的微分效应越强，高频干扰越严重，为了保持较高的系统信噪比，实际选用的值一般不超过20。进一步分析还可以得到，当(即)，的增加就不显著了。当时，超前校正最大补偿相角为，其实是不可能实现的。同时越大，需要增加的放大倍数越大，校正环节的物理实现就越困难。所以，综合考虑上述因素，一般取，即用超前校正补偿的相角一般不超过。4.3.2参数计算由校正前系统伯德图可得系统的开环截止频率为c=8.51rad/s，相位裕度为g=+，幅值裕度Kg=+，相位裕度0=25.2°因此，取附加角=10°。=-0+=50°

15、-25.2°+10°=34.8°令m=34.8°，则=1+sinm1-sinm=1+sin1-sin=1+sin34.8°1-sin34.8°=3.66令20lg|20c'c'4+1|=-10lg3.66，则求得校正后开环截止频率为c'=12.05。也可在校正前系统伯德图上作-10lg3.66db直线与未校正系统对数幅频特性曲线相交， 即可求出校正后开环截止频率为c'=m=12.05rad/s，因此T=1m=0.0434超前校正装置的传递函数为Gcs=1+0.1588s1+0.0434s4.3.3编程设

16、计可利用MATLAB编程，直接得到超前校正装置的传递函数。程序如下：G=tf(20,1/4 1 0); phy=50-25.2+10; %计算phy1=phy*pi/180; %换成弧度a=(1+sin(phy1)/(1-sin(phy1); %求M1=1/sqrt(a);m,p,w=bode(G); %参数返回wm=spline(m,w,M1); %计算mT=1/(wm*sqrt(a);Gc=tf(a*T 1,T 1)其中幅值响应为m，相角响应为p，频率向量为wspline为3次曲线插值函数，通过插值求新的截止频率wm结果：Gcs=1+0.1588s1+0.04339s所得结果与计算相符。M

17、ATLAB编程结果如图4-3所示。图4-3MATLAB程序及结果4.4校正后系统分析4.4.1伯德图加入超前校正装置后，系统的传递函数为：Gcs=1+0.1588s1+0.0434s 80ss+4=12.704s+800.0434s3+1.1736s2+4s用MATLAB绘制伯德图的程序如下： num=12.704,80;>> den=0.0434,1.1736 ,4,0; %开环传递函数的分子、分母多项式>> sys=tf(num,den); %系统建模>> mag,phase,w=bode(sys); >> margin(sys); %画伯德

18、图并求相位穿越频率和截止频率>> grid 由MATLAB绘制出的系统校正后的伯德图如图4-4所示。图4-4校正后系统伯德图由图可知，校正后系统：相位裕度Pm=53.2°截止频率c=12.1rad/s 幅值裕度Gm=+穿越频率g=+4.4.2根轨迹用MATLAB绘制根轨迹的程序如下：num=12.704,80;>> den=0.0434,1.1736 ,4,0;>> rlocus(num,den);>> grid 图4-5为校正后系统的根轨迹。图4-5校正后系统根轨迹5结果对比与分析5.1校正前后阶跃响应曲线用MATLAB画出校正前后系

19、统的单位阶跃响应的程序为num1=20;den1=1/4,1,0;num2=12.704,80;den2=0.0434,1.1736 ,4,0;t=0:0.02:5；numc1,denc1=cloop(num1,den1) y1=step(numc1,denc1,t)numc2,denc2=cloop(num2,den2)y2=step(numc2,denc2,t)plot(t,y1,y2)gridgtext('before')gtext('after')所得阶跃响应曲线如图5-1所示图5-1校正前后阶跃响应曲线5.2结果分析由校正后系统伯德图和根轨迹可以看出，

20、由于采用了串联超前校正，相角裕度由原来的0=25.2°增大到=53.2°，幅值裕度为，满足系统的静态速度误差系数，幅值裕度的要求。由系统单位阶跃响应曲线可以看出：（1）加入校正装置后，校正后系统单位阶跃响应的调节时间大大减小，大大提升了系统的响应速度。（2）校正后系统的超调量明显减小了，阻尼比增大，动态性能得到改善。（3）校正后系统的上升时间减小很多，从而提升了系统的响应速度。综上，串入超前校正装置后，明显地提升了系统的动态性能指标，增强了系统的稳定性。6总结通过这次对控制系统的超前校正的设计与分析，让我对串联滞后校正环节有了更清晰的认识，加深了对课本知识的进一步理解，也让

21、我更进一步熟悉了相关的MATLAB软件的基本编程方法和使用方法。通过复习课本知识以及查阅有关资料确定了整体思路，通过自己演算确定校正网络参数，然后运用MATLAB软件编程验证，作图。在word编辑以及MATLAB软件遇到了一系列问题，通过上网查询或者请教同学都得到了解决，因而设计的过程中也提高了我的软件使用能力。本次课程设计的核心之一是MATLAB软件的使用，通过软件编程，我对MATLAB的语言和应用有了更进一步的了解，特别熟悉了一些对自控很使用的命令，如根轨迹绘制函数rlocus()，伯德图绘制函数bode()等等。利用MATLAB对控制系统进行频域分析，大大简化了计算和绘图步骤，是一款很实

22、用的软件，今后利用课余时间也可以拓展一下自己在MATLAB软件编程的能力。 通过本次课程设计，我体会到了学习自动控制原理，不仅要掌握书本上的内容，还要灵活思考，善于变换，在提出问题、分析问题、解决问题的过程中提高自己分析和解决实际问题的能力。要把理论知识与实践相结合起来，从而提高自己的实际动手能力和独立思考的能力。参考文献1 胡寿松 .自动控制原理（第四版）.北京：科学出版社，20022 王万良 .自动控制原理.北京：科学出版社，20013 黄坚 .自动控制原理及其应用.高等教育出版社，20044 邹伯敏 .自动控制理论.北京：机械工业出版社，20065 李素玲. 自动控制原理.西安：西安电子科技大学出版社，2007附录常用MATLAB函数简介：（1）step函数用MATLAB作控制系统单位阶跃响应时，如果已知系统的传递函数的系数，则可以用step(num,den)或者step(num,den,t)，就得到系统的单位阶跃响应曲线。step(num,den)中没有指定时间t，系统会自动生成时间向量，响应曲线图的坐标也是自动