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

1、武汉理工大学自动控制原理课程设计说明书课程设计任务书学生姓名： 专业班级：指导教师：程 平工作单位：自动化学院题目：用MATLABS行控制系统的滞后超前校正设计初始条件：已知一单位反馈系统的开环传递函数是KG (s)s(0.1s + 1)( 0.02 s + 1)要求系统的静态速度误差系数K/50S-1，丫色40 w弋1°rad /s要求完成的主要任务：(包括课程设计工作量及其技术要求，以及说明书撰写等具体要求)1、MATLAB出满足初始条件的最小 K值的系统伯德图，计算系统的幅值 裕量和相位裕量。2、前向通路中插入一相位滞后一超前校正，确定校正网络的传递函数。3、用MATLA画出未

2、校正和已校正系统的根轨迹。4、用Matlab对校正前后的系统进行仿真分析，画出阶跃响应曲线5、 课程设计说明书中要求写清楚计算分析的过程，列出MATLAB?序和 MATLABt出。说明书的格式按照教务处标准书写。时间安排:任务时间(天)审题、查阅相关资料 2分析、计算3编写程序2撰写报告2论文答辩1指导教师签名：系主任(或责任教师)签名摘要串联滞后-超前校正兼有滞后校正和超前校正的优点，即已校正系统的响应 速度较快，超调量较小，抑制高频噪声的性能也较好。当校正系统不稳定，且要 求校正后系统的响应速度，相角裕度和稳态精度较高时，以采用串联滞后-超前校正为宜。其基本原理是利用滞后-超前网络的超前部

3、分来增大系统的相角裕度， 同时利用滞后部分来改善系统的稳态性能。此次课程设计就是利用 MATLAB对一单位反馈系统进行滞后-超前校正。 通过运用MATLAB的相关功能，绘制系统校正前后的伯德图、根轨迹和阶跃响 应曲线，并计算校正后系统的时域性能指标。关键字：超前-滞后校正MATLAB伯德图 时域性能指标武汉理工大学自动控制原理课程设计说明书目录1滞后-超前校正设计目的和原理1.1.1 滞后-超前校正设计目的 1.1.2滞后-超前校正设计原理 1.2滞后-超前校正的设计过程 3.2.1校正前系统的参数3.2.1.1用MATLAB绘制校正前系统的伯德图42.1.2用MATLAB求校正前系统的幅值裕

4、量和相位裕量 42.1.3用MATLAB绘制校正前系统的根轨迹52.1.4对校正前系统进行仿真分析 62.2滞后-超前校正设计参数计算 7.2.2.1选择校正后的截止频率-'c 82.2.2确定校正参数8.2.3滞后-超前校正后的验证 9.2.3.1用MATLAB求校正后系统的幅值裕量和相位裕量 92.3.2用MATLAB绘制校正后系统的伯德图102.3.3用MATLAB绘制校正后系统的根轨迹112.3.4用MATLAB对校正前后的系统进行仿真分析 123心得体会14参考文献16武汉理工大学自动控制原理课程设计说明书用MATLA进行控制系统的滞后超前校正设计1滞后-超前校正设计目的和原

5、理1.1滞后-超前校正设计目的所谓校正就是在系统不可变部分的基础上， 加入适当的校正元部件，使系统 满足给定的性能指标。校正方案主要有串联校正、并联校正、反馈校正和前馈校 正。确定校正装置的结构和参数的方法主要有两类：分析法和综合法。分析法 是针对被校正系统的性能和给定的性能指标，首先选择合适的校正环节的结 构，然后用校正方法确定校正环节的参数。在用分析法进行串联校正时，校正环节的结构通常采用超前校正、滞后校正和滞后-超前校正这三种类型。超前校正通常可以改善控制系统的快速性和超调量， 但增加了带宽，而滞后 校正可以改善超调量及相对稳定度，但往往会因带宽减小而使快速性下降。滞后 -超前校正兼用两

6、者优点，并在结构设计时设法限制它们的缺点。1.2滞后-超前校正设计原理滞后-超前校正RC网络电路图如图1所示。C1+R1C2T+e(t)m(t)R2图1滞后-超前校正RC网络F面推导它的传递函数:艮1sC2R丄sc1r 1R 丄SC2sC1_ 21 亠 iRG R2C2 RiC2 s R1C1R2C2s1 R1C1s 1 R2C2s令= R1C1 兀二 R2C2 ,AT2 二 R1C1 R2C2R1C2 J 1，则八（l+TsXlfs）s）=1+s （1 十陆s）其中Ti为超前部分的参数，T2为滞后部分。对控制系统进行串联滞后-超前校正的基本原理是利用滞后-超前校正装置 的滞后部分改善控制系统

7、的稳态性能，同时利用其超前部分改善控制系统的动态 性能。在确定校正装置的参数时，两部分基本上可以独立进行。设计串联滞后 - 超前校正装置的基本步骤大致如下：1. 根据控制系统稳态性能指标的要求，确定校正后系统的开环放大数K。2. 基于步骤（1）中的开环放大倍数 K,绘制待校正系统的开环对数幅频特 性L（），并求出此时系统的截止角频率；和相角稳定裕度'。3. 在待校正系统的对数幅频特性上，取斜率从 -20dB/dec变为-40dB/dec的转折角频率-b作为校正装置超前部分的第一转折角频率1/G 2T2），这实际上是利用校正装置中超前部分的比例微分项和未校正系统中的惯性环节相对消的原理。

8、 这样做即可以降低校正后系统的阶次，又可以保证中频区域的斜率为-20dB/dec , 并占有较宽的频带。4. 根据响应的速度的要求，确定校正后系统的截止角频率-'c，并由该角频率处待校正系统的对数幅值L'C c）计算校正装置滞后部分的分度系数:i。要保证 校正后系统的截止角频率为以确定的-c，当- 'b时，下列等式成立20lg ：J L ( c) 20lg( c / b) = 05. 根据滞后部分在校正后系统的截止角频率.c处产生的相角滞后越小越好的原则，确定滞后部分的第二转折角频率 1/G 1T1)，一般可以取1/(：祐)=(1/51/10) c进而计算出T1的值，则

9、滞后-超前校正装置滞后部分的传递函数即可求出为 a1T1 s +1Gc1(s)= _a v1S +1，w S6. 根据系统相角稳定裕度的要求确定超前部分的时间常数T2，这时下列等式成立=180 -( J arctac - arctariT! c arctan： 2T2 c -arctanT2 c式中，(c)是待校正系统在校正后系统截止角频率 'c处的相角。由于上式中只有T2还未知，因此可以解出来，从而得到超前部分的传递函数为,-S2T2S - 1G c2 (s)£ 2，1T2s 1这样，串联滞后-超前校正装置的传递函数可以求出为(qT s +1)(。2泾 +1)2(9毛(沁2

10、(9二丄T 2 q <1口2拥"s -1)(T>s 1)7.校验校正后系统的性能指标。2滞后-超前校正的设计过程2.1校正前系统的参数根据初始条件，调整开环传递函数:Ks 0.1s10.02 s 1当系统的静态速度误差系数Kv = 50S-1时，K=Kv。则K 二 50s满足初始条件的最小K值时的开环传递函数为50s 0.1s10.02 s 12.1.1用MATLAB制校正前系统的伯德图绘制伯德图可用命令bode(num,den)程序：num=50;den=0.002,0.12,1,0;bode( nu m,de n)grid得到的伯德图如图2所示。-100-M-135-

11、1BO-270Fi Muency £iard/MC>e-nor： ningrAims-sr-s曽图2校正前系统的伯德图2.1.2用MATLAB校正前系统的幅值裕量和相位裕量用命令margin(G)可以绘制出G的伯德图，并标出幅值裕量、相位裕量和对 应的频率。用函数kg,r,wg,wc=margin(G) 可以求出G的幅值裕量、相位裕量和 幅值穿越频率。程序：num=50;den=0.002,0.12,1,0;G=tf( nu m,de n);margi n(G)kg,r,wg,wc=margi n(G)得到的幅值裕量和相位裕量如图 3所示Bode DiagramGm - 1SB

12、 dB (at 22.4 rad/&ec). Pm - 3.94 deg (al 20.4启唸m)00-1.£E wpn 毛cnEWif：1010Frequency (raifj'sec)图3校正前系统的幅值裕量和相位裕量又matlab运行界面可以得到运行结果：kg=1.2000wg=22.3607r=3.9431wc=20.3882即幅值裕量：h =20lg1.2 =1.58dB相角裕量：卞-3.9431，穿越频率： =22.3607rad / s ,截止频率：c =20.3882rad /s2.1.3用MATLA绘制校正前系统的根轨迹MATLAB专门提供了绘制根轨

13、迹的有关函数。p,z=pzmap( num,de n)的功能是绘制连续系统的零、极点图。r,k=rlocus(num,den)的功能是绘制k=0；：部分的根轨迹。程序：num=50;den=0.002,0.12,1,0;rlocus( nu m,de n)得到校正前系统的根轨迹如图4所示Root LocusReal Axis图4校正前系统的根轨迹2.1.4对校正前系统进行仿真分析Simulink是可以用于连续、离散以及混合的线性、非线性控制系统建模、 仿真和分析的软件包，并为用户提供了用方框图进行建模的图形接口，很适合 于控制系统的仿真。仿真后得到的结果如图5和图6所示。32* >+50

14、1k时A0 025+1'rSubtractTtarsfer FenTrsnsfei|=cn1Scope图5校正前系统的仿真图图6校正前系统仿真的阶跃响应曲线由Simulink仿真得到的图形可以看出原系统的阶跃响应并不理想。其响应的速度慢而且超调量明显过大，因此要将此系统串联滞后-超前校正来改善其性 能指标。2.2滞后-超前校正设计参数计算221校正装置的传递函数Gc(s)-G 1T1S 1)(-辺T2s 1)(Tis 1)(T2S 1),：-110武汉理工大学自动控制原理课程设计说明书#武汉理工大学自动控制原理课程设计说明书2.2.2确定校正参数先确定领先部分比例微分项的时间常数： 2

15、T2 0为了计算方便和简单起见，又考虑到降低校正后系统的阶次，取2T2二0.1，以对消原系统中的一个惯性环节，这样，图中对数幅频特性斜率从-20dB/dec变为-40dB/dec的转折角频率-_10rad/s 处曲线仍为-20dB/dec,至 =50rad/s变为-40dB/dec。根据设计指标要求，取校正后系统的截止角频率c=15rad/s。在c处下20lg ,l'( .c) - 20lg15 =010成立，可以测算得L'(c) =6.936dB ,于是可以求得1 = 0.3确定校正装置的滞后部分的传递函数。取10.1 c -1.5, :讪则11=0.6675,11 =2.2

16、2s0.667s 12.22s 1则滞后部分的传递函数为这样校正后系统的传递函数为Go(s)=50(0.667s +1)s(0.02s 1)(2.22 s 1)(T2S - 1)根据相角稳定裕度的要求计算出T2。由于= 180 -90 arctan0.66 7 c-arctan0.02 carctan 2.22 carctan T2 c = 45在，c = 15 rad /s时，求得T2 = 0.03s。于是最后得到串联滞后-领先校正装置的传递函数为Gc(s)(0.667s 1)(0.1s 1)(2.22s 1)(0.03s 1)50(0.667 s 1)G0 (s)-校正后的开环传递函数为s

17、(0.02s1)(2.22s1)(0.03s1)2.3滞后-超前校正后的验证由于校正过程中，多处采用的是近似计算，可能会造成滞后 -超前校正后得 到的系统的传递函数不满足题目要求的性能指标。所以需要对滞后 -超前校正后 的系统进行验证。下面用 MATLA求已校正系统的相角裕量和幅值裕量。2.3.1用MATLA求校正后系统的幅值裕量和相位裕量程序：nu m=33.35,50;den=0.001332,0.1116,2.27,1,0;G=tf( nu m,de n);margi n(G)kg,r,wg,wc=margi n(G)得到的校正后系统的幅值裕量和相位裕量如图7所示。运行结果：kg= 5.

18、2549r= 48.3850wg= 39.7404 wc= 13.5142即校正后系统的幅值裕量：h =20lg 5.2549 =14.41dB，相角裕量：=48.385，穿越频率：39.7404 rad / s ,截止频率：飞=13.5142rad /s由以上数据可以看出满足题目要求的性能指标一一系统的静态速度误差系数Kv_5os-1，-40 ,c - 10rad / s。因此设计的滞后-超前校正系统满足要求。Bo-de DiagramGm 皂 14.4 dB (at 3-9.7 rad/sec), Pm » 43.4 deg (at 13.5 rad/呂匕u)10-010 10F

19、requeney (rad/sec)10°22(WE #SEL_d图7校正后系统的幅值裕量和相位裕量假设验证结果不满足指标，重新选择校正后的截止频率，重复上述过程, 直到满足性能指标为止。232用MATLAB制校正后系统的伯德图程序：nu m=33.35,50;den=0.001332,0.1116,2.27,1,0;bode( nu m,de n)grid13武汉理工大学自动控制原理课程设计说明书得到的伯德图如图8所示1Q0Bode Diagram0045&& lW-! F rf-1-I > I 4141I I T MOO1图8校正后系统的伯德图233用MAT

20、LA绘制校正后系统的根轨迹程序：nu m=33.35,50;den=0.001332,0.1116,2.27,1,0;rlocus( nu m,de n)得到的校正后系统的根轨迹如图 9所示。ROOt LOCUBrHa -120 -loo-60Real AXiB-40-ZD20 ODD5 0 5 5I I14武汉理工大学自动控制原理课程设计说明书#武汉理工大学自动控制原理课程设计说明书图9校正后系统的根轨迹#武汉理工大学自动控制原理课程设计说明书234用MATLA对校正前后的系统进行仿真分析用Simulink对校正后的系统仿真。仿真后得到的结果如图10和图11所示。图10校正后系统的仿真图图1

21、1校正后系统仿真的阶跃响应曲线由图像可以看到与图6校正前系统仿真的阶跃响应曲线有明显区别， 其超调 量以及响应速度，调节时间都有明显的改善。还有串联此校正装置得到的性能指 标满足题目的要求。用MATLA编程绘出阶跃响应曲线程序：k=50;num=co nv(0.667,1,0,1);de n=co nv(co nv(co nv(1,0,0.02,1),2.22,1),0.03,1);sys=tf(k* nu m,de n);Lsys=feedback(sys,1,-1);y,t,x=step(Lsys);plot(t,y);grid图12校正后阶跃响应曲线16武汉理工大学自动控制原理课程设计说

22、明书用MATLA编程绘出系统校正前后阶跃响应曲线程序：k1=50;num仁con v(0.667,1,0,1);den 1=co nv(co nv(co nv(1,0,0.02,1),2.22,1),0.03,1);sys仁 tf(k1* nu m1,de n1);Lsys 1=feedback(sys1,1,-1);y1,t1,x1=step(Lsys1);k2=50;num2=co nv(0,1,0,1);den 2=co nv(co nv(1,0,0.1,1),0.02,1);sys2=tf(k2* nu m2,de n2);Lsys2=feedback(sys2,1,-1);y2,t2

23、,x2=step(Lsys2);plot(t1,y1,'-r',t2,y2,'-');grid得到的系统校正前后阶跃响应曲线见图13图13系统校正前后阶跃响应曲线如图13所示，红线表示系统校正后的阶跃响应，蓝线表示未校正前的阶跃响应曲线。很明显串联滞后-超前校正装置以后起动态性能有了明显的改善。从图中可以看出校正后的 超调量c%有了明显的改善，而且上升时间 tr、峰值时间tp、调节时间ts都有了较大的改善。因此此滞后-超前校正装置起到了系统改善的作用。3心得体会此次课程设计主要是运用matlab这个数学工具来解决一些算法和仿真的问 题。能精确地反映系统的阶跃响应， 还能通过图形和数据一起分析，大大减少了 我的工作量。而且也加深了我对 matlab的了解，并从中产生了自己对 matlab学 习的兴趣。此