单位负反馈系统超前校正

1、指导教师评定成绩：审定成绩：自动化学院自动控制原理课程设计报告设计题目：单位负反馈系统的校正设计单位（二级学院）：学生姓名：专业：班级：学号：指导教师：负责项目：设计时间：2013 年12 月目录一、自动控制原理课程设计题目错误 !未定义书签。二、摘要错误 !未定义书签。三、未校正系统的分析错误 !未定义书签。1.系统分析错误 !未定义书签。2.单位阶跃信号下系统输出响应33.未校正系统伯德图4四、系统校正设计51.校正方法52.设计总体思路53.参数确定54.校正装置65.校正后系统76.验算结果7五、结果81.校正前后阶跃响应对比图82.结果分析8六、总结体会错误 !未定义书签。七、参考文

2、献10八、附录101.Matlab 程序101一、自动控制原理课程设计题目题目：单位负反馈系统的校正设计原理：设单位负反馈系统的开环传递函数为：G(s)K /(s(s 1)(0.01s 1)用相应的频率域校正方法对系统进行校正设计， 使系统满足如下动态和静态性能：1）相角裕度450 ；2）在单位斜坡输入下的稳态误差为ess0.0625 ；3）系统的穿越频率大于2rad/s。要求：1）分析设计要求，说明校正的设计思路（超前校正，滞后校正或滞后 - 超前校正）；2）详细设计（包括的图形有：校正结构图，校正前系统的 Bode 图，校正装置的 Bode 图，校正后系统的 Bode 图）；3）用 Mat

3、lab 编程代码及运行结果（包括图形、运算结果）；4）校正前后系统的单位阶跃响应图。1二、摘要用频率法对系统进行超前校正的实质是将超前网络的最大超前角补在校正后系统开环频率特性的截止频率处，提高校正后系统的相角裕度和截止频率，从而改善系统的动态性能。 为此，要求校正网络的最大相位超前角出现在系统的截止频率处。只要正确地将超前网络的交接频率1/aT 和 1/T 设置在待校正系统截止频率 Wc 的两边，就可以使已校正系统的截止频率Wc 和相裕量满足性能指标要求，从而改善系统的动态性能。 串联超前校正主要是对未校正系统在中频段的频率特性进行校正。确保校正后系统中频段斜率等于20dB/dec，使系统具

4、有45°60°的相角裕量。以加快系统的反应速度，但同时它也削弱了系统抗干扰的能力。在工程实践中一般不希望系数a 值很大，当 a20 时，最大超前角为 60°，如果需要 60°以上的超前相角时，可以考虑采用两个或两个以上的串联超前校正网络由隔离放大器串联在一起使用。在这种情况下， 串联超前校正提供的总超前相角等于各单独超前校正网络提供的超前相角之和。关键词：串联超前校正动态性能相角裕度2三、未校正系统的分析1.系统分析校正前单位负反馈系统的开环传递函数为：G(s)K /(s(s 1)(0.01s 1)校正要求：450ess 0.0625 WX 2rad/

5、s1解：因 ess，得 k16，取 k=16。k则：开环传递函数变为G(s) = 16 /(s(s +1)(0.01 s +1)在 Simulink 中建立模型为图 3.1 未校正系统模型图2.单位阶跃信号下系统输出响应在 matlab 中编写 m 程序（ test01jieyue.m），作出阶跃响应曲线：图 3.2 为校正系统单位阶跃响应图3由上图可以看出，系统在阶跃输入下还是可以稳定输出的， 但是开始时 振荡，比较大，超调量也比较大，系统的动态性能不佳。3.未校正系统伯德图在 matlab 中编写 m 程序（ test01bode.m），作出未校正系统的伯德图：图 3.3 未校正新系统伯德

6、图由上图可知，该系统幅值裕度h=16dB，穿越频率Wx=10 rad/s，相角裕度（wc）=12°，截止频率 Wc=3.94 rad/s。4四、系统校正设计1.校正方法由未校正系统伯德图，幅值裕度h=16dB，穿越频率 Wx=10 rad/s，相角裕度（wc）=12°，截止频率 Wc=3.94 rad/s，说明该系统稳定，穿越频率大于要求值，唯相角裕度较小，故采用串联超前校正方法是合适的。超前校正传递函数：2.设计总体思路(1)根据稳态误差 ess 的值，确定开环增益K 。(2) 利用已确定的开环增益 K，确定校正前系统的传递函数，画出其bode 图，从图上看出校正前系统的

7、相位裕度和截止频率 Wc。(3) 根据相位裕度的要 求，计算 出滞后校正装置的参数a 和T。 即得校正装置的传递函数，然后得到校正后系统的开环传递函数。(4)验证已校正系统的相位裕度和截止频率 Wc。3.参数确定（ 1）根据稳态误差要求确定开环k由 ess1，得 k16，取 k=16k（ 2） 求所需最大超前角m因m（ wc ）则m- （wc）45 -12 6 39式中m为校正装置最大超前角为目标相角裕度（wc） 为系统相角裕度为调整量5（ 3） 计算 a 值因 marcsin a - 1 a 11sin得 a1- sinmm4.3955（ 4） 图解最大超前角频率Wm(m)10 lg a为保

8、证响应速度，根据L从未校正系统伯德图中得Wm=5.74 rad/s（ 5）求参数 TT10.0831Wm a4.校正装置则校正装置传递函数为a G c1aTs4.3955 Gc1 4.3955 0.0831s1Ts1 0.0831s在 matlab 中编写 m 程序（ test02bode.m），作出校正系统的伯德图：图 4.1 校正装置伯德图65.校正后系统经超前校正后，系统开环传递函数为：16(10.3653S)G （ s）GC ( S ) * G (S)S(S1)(0.01S1)(10.0831S)6.验算结果在 matlab 中编写 m 程序（ test03bode.m），作出校正后系

9、统的伯德图：图 4.2 校正后系统伯德图由图中可知：校正后系统幅值裕度h=24.1dB，穿越频率 Wx=31.48rad/s，相角裕度（wc）°，截止频率 Wc=5.75 rad/s。=45.67五、结果1.校正前后阶跃响应对比图在 matlab 中编写 m 程序（ test04jieyue.m），作出校正前后系统的阶跃响应对比图：图 5.1 校正前后单位阶跃响应对比图2.结果分析校正后系统相角裕度（wc）=45.6 °>45 °,穿越频率 Wx=31.48 rad/s > 2rad/sess 0.0625，故校正装置满足要求，能完成校正要求。但是相角

10、裕度（wc）=45.6度，略小。故实际应用时，可以考虑将裕量度值取大一点，使a 值再大点。另，因采用超前校正， 故为了补偿无源超前网络的增益衰减，放大器的增益需提高 a 倍。8六、总结体会通过这次的课程设计，再次熟悉我们所学的理论知识，同时为将理论知识运用于实际搭建了一个很好的平台； 使课堂学习的理论知识联系实际， 进一步掌握自动控制原理的时域分析法，根轨迹法，频域分析法，以及各种补偿（校正）装置的作用及用法， 同时也学会使用 MATLAB 语言及 Simulink 动态仿真工具进行系统仿真与调试。在这次课程设计过程中， 遇到一些麻烦。 如对某些概念印象模糊； 软件的开发运用不熟悉。但通过复习

11、书本和向同学请教来了解模糊了的概念以及其他不懂得问题；借阅基本软件指导书籍，借鉴里面的例子快速的掌握此软件的使用方法。 从课程设计的入手到最后分析，对于我们来说是个不小的挑战，要清楚的注意到每个细节是不可行的， 所以遇到问题时和班上的同学一起讨论， 是解决各种难题的有效方法。这次课程设计使我深深地体会团队精神的重要性。通过这次课程设计使我懂得了理论与实际相结合是很重要的， 仅有理论知识是远远不够的，只有把所学的理论知识与实践相结合起来，从理论中得出结论，才能真正为社会服务， 从而提高自己的实际动手能力和独立思考的能力。 在设计过程中遇到的问题很多， 但这是难免的， 只有发现了自己的不足之处并加

12、以改正就可以达到我们需要的要求，不断进步。9七、参考文献1 胡寿松 .自动控制原理 M. 北京：科学出版社 ,2007.6.2 梅晓蓉 .自动控制元件及线路 M. 哈尔滨：哈尔滨工业大学出版社 ,2001.3 张静 .MATLAB 在控制系统中的应用 M. 电子工业出版社 ,2007.4 李友善 .自动控制原理 360 题M. 哈尔滨：哈尔滨工业大学出版社 ,2003.八、附录1. Matlab 程序 test01jieyue.m%* 阶跃响应输出*%s=tf( 's');%建模G=16/(s*(s+1)*(0.01*s+1);G1=feedback(G,1);t=0:0.01

13、:5;step(G1,t);%作图gridxlabel('t');ylabel('c(t)' );title( '单位阶跃响应'); %* 未校正系统伯德图*%s=tf( 's');%建模G=16/(s*(s+1)*(0.01*s+1)figure(1)margin(G);%作伯德图gridHm,Ym,Wc,Wx=margin(G);%参数 test02bode.m%* 校正装置伯德图*%s=tf( 's');%建模10G=(1+0.3653*s)/(1+0.0831*s)figure(1)margin(G);grid%作伯德图gm,pm,wg,wp=margin(G); test03bode.m%* 校正后伯德图*%s=tf( 's');%建模G=(16*(1+0.3653*s)/(s*(s+1)*(0.01*s+1)*(1+0.0831*s)figure(1)margin(G);grid%作伯德图gm,pm,wg,wp=margin(G);test02jieyue.m%* 校正前后阶跃响应对比图*%s=tf( 's');%建模G1=16/(s*(s+1)*(0.01*s+1);G2=feedback(G1,1);t=