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

1、学号：32课程设计题目用MATLAB进行控制系统的超前校正设计学院自动化学院专业自动化班级自动化11班姓名指导教师谭思云2013年 12 月 25 日课程设计任务书学生姓名：刘嘉雯专业班级：自动化 1102班指导教师：谭思云工作单位：自动化学院题目：用MATLAB进行控制系统的超前校正设计。初始条件：已知一单位反馈系统的幵环传递函数是KG（s）s（10.1s）（10.3s）要求系统的静态速度误差系数 Kv 6,45 o要求完成的主要任务:（包括课程设计工作量与其技术要求，以与说明书撰写 等具体要求）1、用MATLAB作出满足初始条件的K值的系统伯德图，计算系统的幅值裕 度和相位裕度。2、在系统

2、前向通路中插入一相位超前校正，确定校正网络的传递函数，并用MATLAB进行验证。3、用MATLAB画出未校正和已校正系统的根轨迹。4、 课程设计说明书中要求写清楚计算分析的过程，列出MATLAB程序和 MATLAB输出。说明书的格式按照教务处标准书写。时间安排：1、 课程设计任务书的布置，讲解（半天）2、根据任务书的要求进行设计构思。（半天）3、熟悉MATLAB中的相关工具（一天）4、系统设计与仿真分析。（三天）5、撰写说明书。 （二天）6、课程设计答辩（半天）指导教师签名： 系主任（或责任教师）签名：目录摘要 01 课程设计目的 12 设计条件与任务要求 12.1 设计条件 12.2 设计任

3、务要求 13 设计基本原理 23.1 超前校正 23.2 根轨迹法 44 设计过程 54.1 基本思路与步骤 54.2 校正前系统分析 64.2.1 开环增益 64.2.2 相位裕度和幅值裕度 64.2.3 伯德图 74.2.4 根轨迹 84.3 超前校正系统设计 94.3.1 理论分析 94.3.2 参数计算 1.0.4.3.3 编程设计 1.1.4.4 校正后系统分析 1.2.4.4.1 伯德图 1.2.4.4.2 根轨迹 1.3.5 结果对比与分析 1.4.5.1 校正前后阶跃响应曲线 1.4.5.2 结果分析 1.5.6 总结 1.7.参考文献 1.8.摘要在自动控制理论中，超前校正是

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

5、提高实 际动手能力。使用 MATLAB软件，学习利用根轨迹与伯德图等对系统进行分析。 旨在培养学生的系统分析能力、设计能力与编程能力，巩固理论所学。2设计条件与任务要求2.1设计条件初始条件如下：已知一单位反馈系统的幵环传递函数是KG（s）s（10.1s）（10.3s）要求系统的静态速度误差系数 Kv 645。2.2设计任务要求（1 ）用MATLAB作出满足初始条件的 K值的系统伯德图，计算系统的幅值裕 度和相位裕度。（2）在系统前向通路中插入一相位超前校正，确定校正网络的传递函数，并 用MATLAB进行验证。（3 ）用MATLAB画出未校正和已校正系统的根轨迹。（4）课程设计说明书中要求写清

6、楚计算分析的过程，列出MATLAB程序和MATLAB输出。说明书的格式按照教务处标准书写。3设计基本原理3.1超前校正超前校正的基本原理就是利用超前相角补偿系统的滞后相角，改善系统的动态 性能，增加相位裕度，提高系统的稳定性。超前校正利用超前网络或 PD控制器进行串联校正，在前向通道中串联传递函 数为：Gc s =，>1其中参数超前校正装置的伯德图如图3-2所示lOlgay*aOdB/d兮/ !/T !诞初吗1叫叫曲i110°图3- 2超前校正网络的伯德图超前校正的目的是改善系统的动态性能，实现在系统静态性能不受损的前提下,提高系统的动态性能。通过加入超前校正环节，利用其相位超

7、前特性来增大系统的 相位裕度，改变系统的幵环频率特性，提高其稳定性。一般使校正环节的最大相位超前角出现在系统新的穿越频率点。如图3-3所示，典型的无源超前装置由阻容元件组成。装置的传递函数为：-I®1（1）C式中°°图3-3典型无源超前装通常a为分度系数，T叫时间常数，由式（1 ）可知，采用无源超前网络进行串 联校正时，整个系统的幵环增益要下降a倍，因此需要提高放大器增益交易补偿。根据式（1），可以得无源超前网络aQ（s）的对数频率特性，超前网络对频率在 1/aT至1/T之间的输入信号有明显的微分作用，在该频率范围内，输出信号相角比输入信号相角超前，超前网络的名称

8、由此而得。在最大超前角频率m处，具有最大超前角m。超前网络（2-1 ）的相角为：c（ ）=arctgaT arctgT（ 2）将上式对 求导并且令其为零，得到最大超前角频率：m “a将上式代入（2），得最大超前角:丄 a 1. a 1m arctgarcs in2、aa 1同时还容易得到m c。最大超前角m仅仅与衰减因子a有关，a值越大，超前网络的微分效果越强。但是a的最大值还受到超前网络物理结构的制约，通常情况下，a取为20左右，这也就意味着超前网络可以产生的最大相位超前约为 65。，如果所需要的大于65 的相位超前角，那么就可以米用两个超前校正网络串联实现，并且在串联的两个网 络之间加入隔

9、离放大器，借以消除它们之间的负载效应。所以通过以上的分析发现，利用超前网络进行串联校正的基本原理，是利用超 前网络的相角超前特性。只要正确的将超前网络的交接频率1/aT或1/T选在待校正系统截止频率的两旁，并适当的选择参数a和T,就可以使已校正系统的截止频率和相角裕度满足性能指标的要求，从而改善系统的动态性能，使校正后的系统具 有以下特点：(1 )低频段的增益满足稳态精度的要求；(2) 中频段对数幅频特性的斜率为 -20db/dec ，并且具有较宽频带，使系统 具有满足的动态性能；(3) 高频段要求幅值迅速衰减，以减少噪声的影响。3.2根轨迹法根轨迹法是一种求特征根的简单方法，在控制系统的分析

10、与设计中得到广泛的 应用。这一方法不直接求解特征方程，用作图的方法表示特征方程的根与系统某一 参数的全部数值关系，当这一参数取特定值时，对应的特征根可在上述关系图中找 到。这种方法叫根轨迹法。根轨迹法具有直观的特点，利用系统的根轨迹可以分析 结构和参数已知的闭环系统的稳定性和瞬态响应特性，还可分析参数变化对系统性 能的影响。在设计线性控制系统时，可以根据对系统性能指标的要求确定可调整参 数以与系统幵环零极点的位置，即根轨迹法可以用于系统的分析与综合。在控制系统的分析和综合中，往往只需要知道根轨迹的粗略形状。由相角条件 和幅值条件所导出的 8条规则，为粗略地绘制出根轨迹图提供方便的途径。根轨迹绘

11、制方法如下：(1 )根轨迹的分支数等于幵环传递函数极点的个数。(2) 根轨迹的始点(相应于K=0 )为幵环传递函数的极点， 根轨迹的终点(相 应于K= %)为幵环传递函数的有穷零点或无穷远零点。(3) 根轨迹形状对称于坐标系的横轴(实轴)。(4) 实轴上的根轨迹按下述方法确定：将幵环传递函数的位于实轴上的极点 和零点由右至左顺序编号，由奇数点至偶数点间的线段为根轨迹。(5) 实轴上两个幵环极点或两个幵环零点间的根轨迹段上，至少存在一个分 离点或会合点，根轨迹将在这些点产生分岔。（ 6）在无穷远处根轨迹的走向可通过画出其渐近线来决定。渐近线的条数等 于开环传递函数的极点数与零点数之差。（ 7 ）

12、根轨迹沿始点的走向由出射角决定，根轨迹到达终点的走向由入射角决 定。（ 8）根轨迹与虚轴（纵轴）的交点对分析系统的稳定性很重要，其位置和相 应的 K 值可利用代数稳定判据来决定。4 设计过程4.1 基本思路与步骤超前校正的方法有很多种，本次课程设计中采用频域法。步骤如下：（1 ）根据稳态误差要求，确定开环增益K。（2）利用已确定的开环增益，计算待校正系统的相角裕度和幅值裕度。nn（3）根据截止频率c的要求，计算a和T。令m c，以保证系统的响应速II 度，并充分利用网络的相角超前特性。显然 m c 成立的条件是m （ c ）L（ c ） Lc（ m） 10lg a根据上式不难求出 a 值，然后

13、确定 T。（4 ）验算已校正系统的相角裕度'。验算时，求得m，再由已知的c"算出待校 正系统在时的相角裕度（c"）。最后，按下式算出m（ c ）如果验算结果不满足指标要求，要重选m，一般使m增大，然后重复以上步骤。4.2校正前系统分析421幵环增益校正前系统幵环传递函数为:G(s)Ks(10.1s)(1 0.3s)题目要求系统的静态速度误差系数Kv 6取Kv=6，根据静态速度误差系数 Kv的计算公式:Kv sim0sG(s)得：Kv=6可知校正前的幵环增益为：K=6 o则校正前的幵环传递函数为：G(s)6s(1 0.1s)(1 0.3s)相位裕度和幅值裕度使幵环频率

14、特性的相角为180的频率，称为相位穿越频率g，即G(j g)H(j g)180使幵环频率特性的幅值为1或者0dB的频率，称为增益穿越频率或者截止频率g，即|G(j c)H(j c)| 1或者20lg |G(j c)H(j c)| 0相角裕度 和幅值裕度Kg或h的定义如下：180 G(j JH(j J1Kg|G(j g)H(j g)|或者h 20lg Kg20lg|G(j g)H(j g)|根据|G(j c)H(j c)| 1和G(j g)H(j g)180可列出如下方程：G(jw) 90arctan0.1w arctan0.3w=-180 °L( c) 20lgA( c) 0求得校正

15、前系统截止频率c 3.74r/s，相位穿越频率g 5.77r/s。于是得出校正前系统的相位裕度和幅值裕度为° 21.2, he 6.94dB 。伯德图由伯德图可得校正前系统的相位裕度和幅值裕度由校正前系统幵环传递函数：G(s)s(1 0.1s)(1 0.3s)用MATLAB绘制伯德图的程序如下：clear all;k0=6;n1=1;d1=co nv(1 0,co nv(0.1 1,0.3 1);mag,phase,w=bode(kO* n1,d1);figure(1);margi n( mag,phase,w);hold onfigure(2);s1=tf(k0* n1,d1);s

16、ys=feedback(s1,1);step(sys)由MATLAB绘制出的系统校正前的伯德图如图4-1所示555羽豈.509035前25-1 畀-1-2Bode DiagramGm - 6.94 dB (at 5.77 rad/aec) , Pm - 21.2 deg (st 3.74 radfeec)Q 0 Q 05 5 Q-1SS355-10u1011O?图4-1校正前系统伯德图由图可知，校正前系统：相位裕度o 21.2 °截止频率c 3.74r/s幅值裕度h。6.94dB穿越频率g 5.77r/s用MATLAB得出的结果与计算结果一致。根轨迹用MATLAB绘制根轨迹的程序如下

17、：num 二1;den 二0.03 0.4 1 0;rlocus( nu m,de n)图4-2为校正前系统的根轨迹图4-2校正前系统根轨迹4.3超前校正系统设计理论分析先求取超前校正的最大超前相角m与取得最大超前相角的频率m。ImImc超前校正时系统的截止频率增大， 未校正系统的相角一般是较大的负相角这里加上一个附加角，用于补偿。超前校正的频率特性相频特性c( ) arctanaT arctanTd c( ) aTTd 1 (aT )21 (T )2令，求得于是丁1丁1m arctanaT arctanT.aT、 aTarctan、. aarcta n即当时，超前相角最大为。可以看出，m只与

18、a有关。由上面的结果可见，m只与a有关。a越大，m越大，对系统的相角补偿越大，超前网络的微分效应越强，高频干扰越严重，为了保持较高的系统信噪比， 实际选用的a值一般不超过20。进一步分析还可以得到，当a 20（即m 65），m 的增加就不显著了。当a时，超前校正最大补偿相角为-，其实是不可能实现的。同时a越大，需要增 加的放大倍数越大，校正环节的物理实现就越困难。所以，综合考虑上述因素,一般取a 520，即用超前校正补偿的相角一般不超过65。28432参数计算取附加角最大超前相角4521.22851.8则求得校正后幵环截止频率为c 6.54rad/s.T=0.0529s也可在校正前系统伯德图上

19、作直线与未校正系统对数幅频特性曲线相交，即可求出校正后幵环截止频率。因此可以得出超前校正的传递函数为：1 si n51.81 -s in 51.88.34编程设计可利用MATLAB编程，直接得到超前校正装置的传递函数程序如下：>>k=6;n1=1;d1=co nv(1 0,co nv(0.1 1,0.3 1);mag,phase,w=bode(k* n1,d1);gama=50;gama1=gama-21.2+30;gam=gama1*pi/180;a=(1-si n( gam)/(1+si n( gam);abd=20*log10(mag);am=20*log10(sqrt(a)

20、;wc=spli ne(abd,w,am);T=1/(wc*sqrt(a);Gc=tf(T 1,a*T 1)MATLAB编程结果如图4-3所示。图4-3 MATLAB显示结果4.4校正后系统分析441伯德图加入超前校正装置后，系统的传递函数为:G(s)6 0.44s 1s 1 0.1s 1 0.3s 0.053s 1用MATLAB绘制伯德图的程序如下:>> n1=6;d1=co nv(1 0,co nv(0.1 1,0.3 1); s1=tf( n1,d1);s2=tf(0.4401 1,0.05278 1);sope二s1*s2;mag,phase,w=bode(sope);ma

21、rgi n( mag,phase,w)由MATLAB绘制出的系统校正后的伯德图如图4-4所示50Bode DiagramOm = 117 dB 14.8 radfeec). Pm = 45.5 deg 6.54 radfcec)O501010Frequericy (rad/sec)5豹£2留匸d10°1010图4-4校正后系统伯德图45.5，截止频率c 6.54rad/s。由图可知，校正后系统：相位裕度442根轨迹用MATLAB绘制根轨迹的程序如下:>> num=6*0.4401 1;de n1=co nv(1 0,0.1 1);de n2=co nv(0.3

22、1,0.05278 1);den二con v(de n1,de n2);G1=tf( nu m,de n);rlocus(G1)图4-5为系统校正后的根轨迹图图4-5校正后系统根轨迹5 结果对比与分析5.1 校正前后阶跃响应曲线用 MATLAB 画出校正前后系统的单位阶跃响应的程序为 num1=6;den1=0.03,0.4,1,0;num3=2.6406,6;den3=0.00158 ， 0.0512,0.453,1,0; t=0:0.02:5numc1,denc1=cloop(num1,den1) y1=step(numc1,denc1,t) numc3,denc3=cloop(num3,

23、den3) y3=step(numc3,denc3,t) plot(t,y1,y3);gridgtext('before') gtext('after') 所得阶跃响应曲线如图 5-1 所示。图5-1校正前后阶跃响应曲线5.2结果分析由校正后系统伯德图和根轨迹可以看出，由于采用了串联超前校正，相角裕度由原来的21.2 °增大到45.5。，满足系统的静态速度误差系数 Kv 6，45的要求由系统单位阶跃响应曲线可以看出：（1）加入校正装置后，校正后系统单位阶跃响应的调节时间大大减小，大大 提升了系统的响应速度。（2）校正后系统的超调量明显减小了，阻尼比增大

24、，动态性能得到改善。（3）校正后系统的上升时间减小很多，从而提升了系统的响应速度。综上，串入超前校正装置后，明显地提升了系统的动态性能指标，增强了系统 的稳定性。6总结通过这次对控制系统的超前校正的设计与分析，我对串联超前校正环节有了更 清晰的认识，加深了对课本知识的进一步理解，也让我更进一步熟悉了相关的MATLAB 软件的基本编程方法和使用方法。 通过复习课本知识以与查阅有关资料确 定了整体思路，通过自己演算确定校正网络参数，然后运用 MATLAB 软件编程验 证，作图。在 word 编辑以与 MATLAB 软件遇到了一系列问题，通过上网查询或 者请教同学都得到了解决， 因而设计的过程中也提

25、高了我的各方面能力， 收获很大。本次课程设计的核心之一是 MATLAB 软件的使用， MATLAB 是一款功能十分 强大的软件，不仅在信号与系统上有充分的应用，在自动控制上也显示出了它强大 的功能。用 MATLAB 可以绘制跟轨迹，伯德图，那奎斯特图，还可以仿真。平时 上课只知道课本上的理论知识，其实 MATLAB 就是一个理论与实际联系的很好的 工具，通过学会用 MATLAB 我可以更好的理解课本上的理论知识。 总之，通过本次课程设计我体会到了学习自动控制原理，不仅要掌握书本 上的内容，还要灵活思考，善于变换，在提出问题、分析问题、解决问题的过程中 提高自己分析和解决实际问题的能力。要把理论

26、知识与实践相结合起来，从而提高 自己的实际动手能力和独立思考的能力。并且以后一定要多读多写 MATLAB 的各 种程序，以便在将来的工作学习中能更好的应用！参考文献1 胡寿松 .自动控制原理（第四版） .北京：科学出版社， 20022 王万良 .自动控制原理 .北京：科学出版社， 20013 宋乐鹏 .自动控制原理 .北京：清华大学出版社， 20124 马莉 .MATLAB 语言实用教程 .北京：清华大学出版社， 20105 黄坚 .自动控制原理与其应用 .北京：高等教育出版社， 2004附录常用 MATLAB 函数简介(1 )step 函数用 MATLAB 作控制系统单位阶跃响应时，如果已知系统的传递函数的系数，则 可以用 step(num,den) 或者 step(num,den,t) ，就得到系统的单位阶跃响应曲线。 step(num,den) 中没有指定时间 t ，系统会自动生成时间向量，响应曲线图的坐标 也是自动标注的。若采用命令 step(num,den,t) 求系统的单位阶跃响应，其中的时 间t由用户指定，MATLAB会根据用户给定的时间t，算出对应的坐标值。执行该 命令不能自动画出系统的单位阶跃响应图，而要令加 plot 绘图命令。除了求单位阶跃响应的 step 函数， MATLAB 中还提供了求系统各种响应的函 数，例如求脉冲响应的 implus