《单位负反馈系统校正自动控制原理课程设计》由会员分享，可在线阅读，更多相关《单位负反馈系统校正自动控制原理课程设计（12页珍藏版）》请在金锄头文库上搜索。
1、1. 设计题目已知单位负反馈系统被控制对象的开环传递函数G (s) = K0s (0.5s +1)用相应的频率域校正方法对系统进行校正设计，使系统满足如下动态及静态性能 指标：(1) 选取相应的频率域校正方法(2) 在斜坡信号r(t) = 2t作用下，系统的稳态误差e 5Oo。c(4) 当时，系统开环对数频率特性，不应有斜率超过-40dB /十c低频的线段。要求：(1) 分析设计要求，说明校正的设计思路(滞后校正，超前校正或滞后-超前校正)；(2) 详细设计(包括的图形有：校正结构图，校正前系统的Bode图，校正 装置的 Bode 图，校正后系统的 Bode 图)；(3) MATLAB编程代码
2、及运行结果(包括图形、运算结果)；(4) 校正实现的电路图及结果(校正前后系统的阶跃响应图)；2. 设计报告正文2.1摘要利用超前网络或PD控制器进行串联校正的基本原理，是利用超前网络和PD 控制器的相角超前特性。只要正确的将超前网络的交接频率1/aT和1/T选在待 校正系统截止频率的两旁，并适当选取a和T,就可以是已校正系统的截止频率 和相角裕度满足性能指标的要求，从而改善闭环系统的动态性能。闭环系统的稳 态性能要求，可通过选择已校正系统的开环增益来保证。关键词：稳态误差e，相位裕量厂，超前校正ss2.2设计思路1）根据稳态误差e要求，确定开环增益K。ss2）利用已确定的开环增益K，计算待校
3、正系统的相位裕度r。3）根据截止频率wc”的要求，计算超前网络参数a和T。在本步骤中，关键 是选择最大超前角频率等于要求的系统截止频率，即wm=wc，以保证系统的响 应速度，并充分利用网络的相角超前特性。显然，wm=wc成立的条件是-L (wm ”) = Lc (wm) =10lga根据上式不难求出a值，然后由wm. a确定T。4）验证已校正系统的相位裕度丫。由于超前网络的参数书根据满足系统截 止频率要求选择的，因此相角裕度是否满足要求，必须验证。验算时，由已知a值由. a 10 二arcsina +1求的a值，再由已知的wc”算出代校正系统在wc时的相角裕度丫（wc”）。2.3报告正文第一步
4、：由题意可知，开环传递函数在斜坡信号r（t）=21作用下，要满足系统的稳态误差e 0.02。根据稳态误差公式sse =lim sE(s) =sslimSR(S)1 + G (s) H (s)可得到系统的稳态误差e S) = 100。题目给出的开环传递函数的二阶函数，故不能用劳斯判定法来进 一步缩小K的取值范围。取K=100 (rad) -1将其带入原来的开环传递函数得待 校正系统开环传递函数为G(s)=100s(0.5s +1)第二步：上式代表最小系统，可通过MATLAB绘制出待校正系统的bode图和 待校正系统的单位阶跃响应分别如下图1和图2。00Bode DiaaramGm = Inf d
5、B (at Inf radXsec) , Pm = 8.09 deg (at 1 4.1 rad/sec)41:1ooO I2 4 9 一 _10 101Frequency (radXsec)102图1待校正系统的bode图itep Response234Tim已(sec)5e64.21-ou o.ILIp-fcJ-dE=l6 4 o.o.o2 o_0.wma1 + aTs1 + Ts9.8281 Gc1 + 0.1255s1 + 0.0127s图2待校正系统的单位阶跃响应由图1和图2可知待校正系统的频域性能指标如下：幅值稳定裕度Lh=+ g dB；-兀穿越频率wx=+ g rad/s;相角稳
6、定裕度r=8.09度；截止频率wc=14.1rad/s。第三步:假设wm=wc =25rad/s，由待校正系统的bode图查的L(wm)=-9.93dB, 再根据公式-L (wm)=10lga,计算出 a= 9.8281, T=一二=0.0127。因此，再根据超前网络传递函数公式a G (s) =c得为了补偿无源超前网络产生的增益衰减，放大器的增益需提高9.8281倍，否则 不能保证稳态误差要求。超前网络参数确定后，已校正系统的开环传递函数为G (s)= G (s) G (s) = I。0。+ 0.1255s)_ 0cs (0.5s +1)(1 + 0.0127 s)算出待校正系统的r (wc
7、)=3.8由最大超前角的计算公式0 二arcsin可求得0 =55.4，故已知系统的相角裕度r=0 + r (wc”)=55.4 + 3.8 =59.250。再用MATLAB绘制出校正装置的bode图(图3)ESBede DiagramGm = Int , Frn = -180 deg (at Int rad/sec)o100 121 0 1 0 1 0103Frequency (rad/sec)3CEiQpjQJHlnJZd10*图3校正装置的bode图第四步：验证校正系统的相角裕度r”是否满足要求。Bode DiagramGm = Inf dB (at Inf rad/sec) , Pm
8、= 59.2 deg (at 25 rad/sec)(EP)电品匚尉壬n- o OO108 1 2 10 10 10101QFrequency frad/aec)图4校正后系统的bode图44O1 1图5校正后系统的阶跃响应图由MATLAB绘制的校正后系统的bode图（图4）可知，相角裕度 r”=59.2。50。，此时相位裕量性能指标满足。通过对前后系统的bode图（图6）的比较可知，当时，校正后的开 c环对数频率特性没有斜率超过-40dB /十低频的线段。因此，全部性能指标都满足。正前后系统的bode图如下：woBode Diagram101io2Frequency (rad/sec)3i1
9、0 10505L(EiQp.jILIwasJZdo o o A-1 -1 10 10图6校正前后系统bode图（虚线校正前；实线校正后）3. 设计总结在本次课程设计中，我不仅加深了对自动控制原理课程的理解，将理论应用 到实际中去，而且我还学会了如何去培养我们的创新精神，从而不断地战胜自己、 超越自己。创新，是我们学会将理论很好地联系实际，并不断地去开动自己的大 脑，做自己所能及的，但别人没有想到的事。设计过程中，也好比是我们人类成 长的过程，常有不如意的事，也许这就是在对我们提出挑战，勇敢地面对他，胜 利的钟声也一定会为我们儿敲响。通过本次设计，我了解了自动控制原理中校正系统的基本概念及其对系
10、统设 计的相关应用，通过对初步知识的了解，对校正系统各种方案的比较，进一步了 解了校正系统的合理性和实用性。什么样的课程设计都离不开理论与实际相结合 的真理，设计过程中的方案选择和参数设定使我进一步深刻认识到自控原理中校 正环节对整个系统的重要作用。一个细小的参数设定出现偏差，可能导致最后的 性能指标不和标准。所以选择一个优良的方案结于实验至关重要。我认为，在设计时应该怎样少走一些弯路，怎样能够非常透彻的理解系统并 用简单方法设计校正装置，我想这是这次课程设计最锻炼人的地方。然而这也要 求我们有相当厚实的理论基础，并能很好地运用到实际中去。这是我们学习和掌 握好自控原理最重要的。我们运用Mat
11、lab软件进行系统仿真验证，这不仅对我 们设计带来了方便，也能很准确地为我们改动参数提供依据，同时也让我们对 Ma tlab软件进行了又一步的学习，也为我们再次熟练运用Mat lab打下了基础。 这次课程设计也锻炼了我们的细心和耐心，这尤其体现在我们课程设计报告的格 式上，很多人格式上不过关，关键是他们太浮躁，不能细心地对待。格式很重要， 不然，美观切不说，报告看上去很混乱，不易于阅读。严格的格式要求也是在培 养我们的做事态度。我也很感谢有课程设计这次机会，让我明白了这么多从前不知道的东西。这 些对我日后的工作都是十分宝贵的。同时也让我看清了自己，明白了自己哪里欠 缺。理论知识的不足在这次课程
12、设计中给我带来了很多麻烦，这也算是提醒。今 后，在学习中，我要端正自己的态度努力学习，只有这样我们才能真正的掌握好 知识。我用的是串联超前校正装置系统进行校正，从中我认识到了超前校正的重要 作用，通过它可以提高系统的截止频率和相角裕量，从而减小了阶跃响应的超调 量和调节时间，更进一步改善系统的动态性能指标。通过此次的课程设计，让我 进一步巩固了超前校正、滞后校正和超前滞后校正三种校正的具体步骤。 参考文献1 胡寿松自动控制原理（第五版） 北京：科学出版社.20072 张静.MATLAB在控制系统中的应用.北京：电子工业出版社.20075 黄坚主. 自动控制原理及其应用. 北京：高等教育出版社
13、2004附录：MATLAB编程代码G0=tf(100,0.5, 1 0); figure(1);margin(G0);hold on figure(2);sys=feedback(G0,1);step(sys) wm=25;L=bode(G0,wm);Lwc=20*log10(L) a=10(-0.1 * Lwc)T=1/(wm*sqrt(a);phi=asin(a-1)/(a+1)Gc=(1/a)*tf(a*T 1,T 1);figure(3);margin(Gc);Gc=a*Gc;figure(4);margin(Gc);G=Gc*G0;bode(G,r,G0,b-);grid figur
14、e(5);sys=feedback(G,1);step(sys) figure(6);margin(G);T=1/(wm*sqrt(a)Gc=(1/a)*tf(a*T 1,T 1)MATLAB运行结果G0=tf(100,0.5 ,1 0);figure(1);margin(G0);hold on figure(2);sys=feedback(G0,1);step(sys) wm=25;L=bode(G0,wm);Lwc=20*log10(L) a=10(0.1 * Lwc)T=1/(wm*sqrt(a); phi=asin(a-1)/(a+1) Gc=(1/a)*tf(a*T 1,T 1); figure(3);margin(Gc);Gc=a*Gc; figure(4);margin(Gc);G=Gc*G0; bode(G,r,G0,b-);grid figure(5);sys=feedback(G,1);step(sys) figure(6);margin(G);T=1/(wm*sqrt(a)Gc=(1/a)*tf(a*T 1,T 1)Lwc =-9.9247a =
