1、武汉理工大学自动控制原理课程设计说明书课程设计任务书学生姓名: 李萌 专业班级: 电气 1101 指导教师: 肖纯 工作单位: 自动化学院 题 目: 基于频率法的超前校正设计 初始条件:已知系统的传递函数模型为:)2(104)(0ssG要求完成的主要任务: (包括课程设计工作量及其技术要求,以及说明书撰写等具体要求)(1) 掌握采用频率法设计超前校正装置的具体步骤;(2) 设计超前校正环节,使其校正后系统的静态速度误差系数 ,相角裕度6vK大于 ;48(3) 采用 Matlab 工具进行分析设计,并绘制校正前后系统的单位阶跃响应曲线,开环 Bode 图和 Nyquist 图;(4) 分析比较采
2、用校正前后的 Bode 图和 Nyquist 图,说明其对系统的各项性能指标的影响。总结频率法校的优缺点及其适应条件; (5) 对上述任务写出完整的课程设计说明书,说明书中必须写清楚分析计算的过程,并包含 Matlab 源程序或 Simulink 仿真模型,说明书的格式按照教务处标准书写。时间安排:任务 时间(天)指导老师下达任务书,审题、查阅相关资料 1.5分析、计算 2.5编写程序 2撰写报告 1.5论文答辩 0.5指导教师签名: 年 月 日武汉理工大学自动控制原理课程设计说明书系主任(或责任教师)签名: 年 月 日武汉理工大学自动控制原理课程设计说明书目录摘要 .11 串联超前(微分)校
3、正 21.1RC 超前网络 .21.2 串联超前校正 32 控制系统校正前分析 42.1 用 MATLAB 做出校正前系统的阶跃响应曲线、伯德图和奈奎斯特图 .42.1.1 系统的开环传递函数 42.1.2 校正前系统的单位阶跃响应曲线 .42.1.4 校正前系统的奈奎斯特图 .62.2 稳态误差分析 63 控制系统超前校正分析设计 73.1 串联超前校正原理分析 73.2 理论计算 .83.3 采用 MATLAB 工具进行串联超前校正设计 93.3.1 利用 MATLAB 进行超前校正设计的程序 93.3.2 利用 MATLAB 工具设计超前校正结果 113.3.2.1 校正环节传递函数及校
4、正后系统传递函数 114 控制系统校正前后的对比 .154.1 控制系统校正前后的阶跃响应曲线、伯德图和奈奎斯特图对比 154.1.1 系统校正前后单位阶跃曲线对比及分析 .154.1.2 系统校正前后伯德图与奈奎斯特图对比 .165 频率法校正优缺点及适用条件 185.1 频率法超前校正的优缺点及适用条件 .185.1.1 频率法超前校正的优缺点: 185.1.2 频率法超前校正的适用条件: 185.2 频率法校正的其他情况 185.3 频率法校正的优缺点及适用条件 196 课程设计总结 207 参考文献 21武汉理工大学自动控制原理课程设计说明书0摘要自动控制理论是研究自动控制共同规律的技
5、术科学。它的发展初期,是以反馈理论为基础的自动调节原理,并主要用于工业控制。控制系统通常是由被控对象、控制器和检测环节组成,一般情况下,仅仅依靠被控对象本身的特性不可能同时满足对系统所提出的各项性能指标的要求。因此需引入校正装置改善系统的性能。本文是利用自动控制原理中所学的知识,结合 MATLAB 软件,对系统进行频率法超前校正以使其性能指标满足要求,并对超前校正进行一定分析。关键字:传递函数 相位裕度 频率法 超前校正 武汉理工大学自动控制原理课程设计说明书11 串联超前(微分)校正1.1RC 超前网络RC 超前网络的传递函数为(1)1)(TsaUrc式中 (2);221RCR为了讨论问题的
6、方便,在网络前(或后)附加一个放大器,使其放大系数等于 。这样的超前网络作为校正装置,其传递函数可看成为(3)1)(TsaGc超前网络对频率在 1/aT1/T 之间的输入信号有明显的微分作用,在该频率范围内,输出信号的相角超前于输入信号的相角。超前网络的名称也由此而来。当频率等于最大超前角频率 时,相角超前量最大,以 表示。而 又恰好是频率aT 和 1/T 的几何mmm中点,即(4)aTaTm1lg)l1(lg2l 因此 (5)m最大超前角为(6)Tarctgarctgmmm根据三角函数两角求和公式,可解得或 (7)artm211rsinam上式表明, 仅与 值有关, 值选得越大,则超前网络的
7、微分效应越强。实际选用的 值必须考虑到网络物理结构的限制及附加放大器的放大系数等原因,一般取值不大于 20。此外, 处的对数幅值为m ()ajGLmcmlg0|)(|lg20武汉理工大学自动控制原理课程设计说明书21.2 串联超前校正利用超前网络进行串联校正的基本原理,乃是利用超前网络相角超前特性。只要正确地将超前网络的交接频率 1/aT 和 1/T 设置在待校正系统截止频率 的两边,就可以使c已校正系统的截止频率 和相裕量满足性能指标要求,从而达到改善系统动态性能的目c的。串联超前校正设计的一般步骤为(1)根据稳态误差要求,确定开环增益 K。(2)在已确定 K 值条件下,计算未校正系统的相裕
8、量。(3)根据指标要求,确定在系统中需要增加的相角超前量。(4)由式(7)确定 值及 值,在未校正系统的对数幅频特性曲线上找到幅值等mL于 点所对应的频率,即为 ,这一频率为所选网络的 ,并且在此频率上将产生mLcc最大超前相角值 。m()确定超前网络的交接频率 1/aT , 1/T 。12(6)验算。武汉理工大学自动控制原理课程设计说明书32 控制系统校正前分析2.1 用 MATLAB 做出校正前系统的阶跃响应曲线、伯德图和奈奎斯特图2.1.1 系统的开环传递函数系统的开环传递函数如下: )2(104)(0ssG化成时间常数形式如下:)5.)(.)(0ss2.1.2 校正前系统的单位阶跃响应
9、曲线用 MATLAB 绘制校正前系统单位阶跃响应曲线的程序如下:num=2;den=conv(1 0,conv(0.1 1,0.5 1); %分母多项式展开W=tf(num,den); %开环传递函数Wy_c=feedback(W,1) %校正前开环系统传递函数step(Wy_c,r,5); %校正前单位阶跃响应曲线grid on;系统校正前单位阶跃响应曲线。图 1 校正前系统的阶跃响应曲线Step ResponseTime (seconds)Amplitude0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 500.20.40.60.811.21.4System: Wy_cRise
10、 time (seconds): 0.766System: Wy_cPeak amplitude: 1.25Overshoot (%): 25.4At time (seconds): 1.85System: Wy_cSettling time (seconds): 4.4System: Wy_cFinal value: 1武汉理工大学自动控制原理课程设计说明书4由图可知,上升时间 ,超调时间 ,调节时间 s,超调量str76.0st85.1p4.st。%4.252.1.3 校正前系统的波德图用 MATLAB 绘制校正前系统波德的程序如下:num=2;den=conv(1 0,conv(0.1
11、1,0.5 1); %分母多项式展开W=tf(num,den); %开环传递函数Wy_c=feedback(W,1) %校正前开环系统传递函数bode(W,r); %校正前开环系统伯德图grid on;系统校正前的波德图。图 2 校正前系统伯德图系统校正前的波德图如图 2 所示。由图可知,此时截止频率 Wc=1.56rad/s,相位裕度Pm=43.2,相位穿越频率 Wg=4.47rad/s,幅值裕度 GM=15.6dB。Bode DiagramFrequency (rad/s)-150-100-50050 System: WGain Margin (dB): 15.6At frequency
12、(rad/s): 4.47Closed loop stable? YesMagnitude (dB)10-1 100 101 102 103-270-225-180-135-90System: WPhase Margin (deg): 43.2Delay Margin (sec): 0.484At frequency (rad/s): 1.56Closed loop stable? YesPhase (deg)武汉理工大学自动控制原理课程设计说明书52.1.4 校正前系统的奈奎斯特图用 MATLAB 绘制校正前系统波德的程序如下:num=2;den=conv(1 0,conv(0.1 1,0.
13、5 1); %分母多项式展开W=tf(num,den); %开环传递函数Wy_c=feedback(W,1) %校正前开环系统传递函数figure(3);nyquist(W,r); %校正前开环系统奈奎斯特图grid on;校正前系统的奈奎斯特图。图 3 校正前系统的奈氏曲线由图可知原系统稳定。2.2 稳态误差分析已知系统为 1 型系统,当为斜坡信号输入时,则有: 2)15.0)(.lim)15.0)(.2lim)(li 000 sssHsGKssv,满足题意,且 K 不变,则 不变。6v v-1.4 -1.2 -1 -0.8 -0.6 -0.4 -0.2 0-20-15-10-5051015
14、20 0 dB-6 dB4-2 dB6 dB4 dB2 dBNyquist DiagramReal AxisImaginary Axis武汉理工大学自动控制原理课程设计说明书63 控制系统超前校正分析设计3.1 串联超前校正原理分析利用超前校正网络或 PD 控制器进行串联校正的基本原理,是利用超前校正网络或 PD控制器的相角超前特性,只要正确地将校正网络的交接频率 1/ T 和 1/T 选在待校正系统截止频率的两旁,并适当选择参数 和 T,就可以使已校正系统的截止频率和相角裕度满足性能指标的要求,从而改善闭环系统的动态性能。闭环系统的稳态性能要求,可以通过选择已校正系统的开环增益来保证。用频域
15、法设计无源超前网络的步骤如下:(1)根据稳态误差要求,确定开环增益 K;(2)利用已确定的开环增益,计算未校正系统的相角裕度;(3)由校正后的相角裕度确定 a 值;(4)由 a 值计算出 的值,然后计算出 T 以及 aT 的值;c(5)得出校正函数和校正后的传递函数。在本步骤中,关键是选择最大超前相角频率等于要求的系统截止频率,即 ,cm以保证系统的响应速度,并充分利用网络的相角超前特性。显然 成立的条件是cm(9)lg10)()( mcL根据上式求出 值,然后由mT1(10)确定 T 值。如果对校正后系统的截止频率 未提出要求,可以从给出的相角裕度 出发,通过c下式求得:(11)mm式中,
16、为利用超前校正网络产生的最大超前相角; 为系统所要求的相角裕度;m 为未校正系统在 时的相角裕度;选择 为 是由于加入超前校正装置后,cc125对数幅频特性向右移动,为补偿由此而造成的相角滞后,需要在相角裕度上加以修正。武汉理工大学自动控制原理课程设计说明书7求出校正装置的最大超前相角 后,根据下式:m(12)masin1求得 值。在未校正系统的对数幅频特性上计算其幅值等于所对应的频率就是校正后系统的截止 ,且 。ccm由以上各式计算出校正装置的参数 ,T 和 T,并以此写出校正装置应具有的传递c函数。校正系统的传递函数为:(13)asGc1)(校正后系统的传递函数为:(14))()(Tssc
17、当完成校正装置设计后,需要进行实际系统调试工作,或者进行计算机仿真以检查实际系统的响应特性。这时,需将系统建模时省略的部分尽可能加入系统,以保证仿真结果的逼真度。如果由于系统各种固有非线性因素影响,或由于系统噪声和负载效应等因素的影响,使已校正系统不能满足全部性能指标要求,则需要适当调整校正装置的形式或参数,直到满足全部性能指标要求为止。3.2 理论计算(1)给定 K 值K=2(2)未校正系统的相位裕度 1)()(cccjHjGA180(求得: =1.56rad/sc=43.2(3)根据未校正前系统的相角裕度和题目初始条件规定的相角裕度值,可以得到最大超前角,由 计算,取规定相角裕度为 48.
18、2,取补偿角为 10。mm可得最大超前角为:武汉理工大学自动控制原理课程设计说明书8=48.2-43.2+10=15mm(4)由 求得asin11.70(5)根据公式 求得lg10)()(mccL=1.88rad/sc根 据公式求得mT1T=0.408s所以系统的超前校正装置的传递函数为: sTsGsc408.169校正后系统的开环传递函数为 ).)(5.)(.020 ssssc(6)校验:系统校正后的相位裕度: )(c180-90-arctan0.1*1.88-arctan0.3*1.88-arctan0.408*1.88+arctan0.694*1.8852.533.3 采用 MATLAB
19、 工具进行串联超前校正设计3.3.1 利用 MATLAB 进行超前校正设计的程序num=2;den=conv(1 0,conv(0.1 1,0.5 1); %分母多项式展开W=tf(num,den); %开环传递函数kc=1; %稳态误差系数不变yPm=48.2+10; %增加量取 10deg武汉理工大学自动控制原理课程设计说明书9W=tf(W); %超前校正环节mag,pha,w=bode(W*kc); %开环频率特性的幅值和相位值Mag=20*log10(mag); %幅值的对数值的 20 倍Gm,Pm,Wcg,Wcp=margin(W*kc); %幅值稳定裕度 Gm,相位稳定裕度 Pm
20、和相应的交接频率 Wcg 和 Wcpphi=(yPm-Pm)*pi/180; %确定 m 值alpha=(1+sin(phi)/(1-sin(phi); %确定 a 的值Mn=-10*log10(alpha); %a 的对数值的负 10 倍Wcgn=spline(Mag,w,Mn); %确定最大相角位移频率T=1/Wcgn/sqrt(alpha); %求 T 值Tz=alpha*T;Wc=tf(Tz 1,T 1) %超前校正环节的传递函数Wy_c=feedback(W*kc,1) %校正前开环系统传递函数Wx_c=feedback(W*kc*Wc,1) %校正后开环系统传递函数figure(1
21、);step(Wy_c,r,5); %校正前开环单位阶跃响应曲线hold on;step(Wx_c,b,5); %校正后开环单位阶跃响应曲线grid on;figure(2);bode(W*kc,r); %校正前开环系统伯德图hold on;bode(W*kc*Wc,b); %校正后开环系统伯德图grid on;figure(3);nyquist(W*kc,r); %校正前开环系统奈奎斯特图hold on;nyquist(W*kc*Wc,b); %校正后开环系统奈奎斯特图grid on;武汉理工大学自动控制原理课程设计说明书103.3.2 利用 MATLAB 工具设计超前校正结果3.3.2.1
22、 校正环节传递函数及校正后系统传递函数校正环节传递函数:0.6965 s + 1-0.4102 s + 1校正前系统闭环传递函数:2-0.05 s3 + 0.6 s2 + s + 2校正后系统闭环传递函数:1.393 s + 2-0.02051 s4 + 0.2961 s3 + 1.01 s2 + 2.393 s + 23.2.3.2 校正后系统的阶跃响应曲线、伯德图和奈奎斯特图用 MATLAB 绘制校正后系统阶跃响应曲线的程序如下 :num=2;den=conv(1 0,conv(0.1 1,0.5 1); %分母多项式展开W=tf(num,den); %开环传递函数kc=1; %稳态误差系
23、数不变yPm=48.2+10; %增加量取10degW=tf(W); %超前校正环节mag,pha,w=bode(W*kc); %扩大系数后的开环频率特性的幅值和相位值Mag=20*log10(mag); %幅值的对数值的20倍Gm,Pm,Wcg,Wcp=margin(W*kc); %幅值稳定裕度Gm,相位稳定裕度Pm和相应的交接频率Wcg和Wcpphi=(yPm-Pm)*pi/180; %确定m值alpha=(1+sin(phi)/(1-sin(phi); %确定a的值武汉理工大学自动控制原理课程设计说明书11Mn=-10*log10(alpha); %a的对数值的负10倍Wcgn=spli
24、ne(Mag,w,Mn); %确定最大相角位移频率T=1/Wcgn/sqrt(alpha); %求T值Tz=alpha*T;Wc=tf(Tz 1,T 1) %超前校正环节的传递函数Wx_c=feedback(W*kc*Wc,1) %校正后开环系统传递函数step(Wx_c,b,5); %校正后开环单位阶跃响应曲线grid on;校正后系统的阶跃响应曲线。图 4 校正后系统的阶跃响应曲线由图可知,上升时间 ,超调时间 ,调节时间 3.29s,超调量str652.0stp5.1st。%2.16用 MATLAB 绘制校正后系统阶跃响应曲线的程序如下 :num=2;den=conv(1 0,conv(
25、0.1 1,0.5 1); %分母多项式展开W=tf(num,den); %开环传递函数kc=1; %稳态误差系数不变yPm=48.2+10; %增加量取10degW=tf(W); %超前校正环节Step ResponseTime (seconds)Amplitude0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 500.20.40.60.811.21.4System: Wx_cSettling time (seconds): 3.29System: Wx_cPeak amplitude: 1.16Overshoot (%): 16.2At time (seconds): 1.5S
26、ystem: Wx_cRise time (seconds): 0.652System: Wx_cFinal value: 1武汉理工大学自动控制原理课程设计说明书12mag,pha,w=bode(W*kc); %扩大系数后的开环频率特性的幅值和相位值Mag=20*log10(mag); %幅值的对数值的20倍Gm,Pm,Wcg,Wcp=margin(W*kc); %幅值稳定裕度Gm,相位稳定裕度Pm和相应的交接频率Wcg和Wcpphi=(yPm-Pm)*pi/180; %确定m值alpha=(1+sin(phi)/(1-sin(phi); %确定a的值Mn=-10*log10(alpha);
27、 %a的对数值的负10倍Wcgn=spline(Mag,w,Mn); %确定最大相角位移频率T=1/Wcgn/sqrt(alpha); %求T值Tz=alpha*T;Wc=tf(Tz 1,T 1) %超前校正环节的传递函数Wx_c=feedback(W*kc*Wc,1) %校正后开环系统传递函数bode(W*kc*Wc,b); %校正后开环系统伯德图grid on;校正后的bode 图。Bode DiagramFrequency (rad/s)10-1 100 101 102 103-270-225-180-135-90System: untitled1Phase Margin (deg):
28、51.3Delay Margin (sec): 0.479At frequency (rad/s): 1.87Closed loop stable? YesPhase (deg)-150-100-50050System: untitled1Gain Margin (dB): 15.3At frequency (rad/s): 5.55Closed loop stable? YesMagnitude (dB)图 5 校正后系统的 bode 图由图可知,此时截止频率 Wc=1.87rad/s,相位裕度 Pm=51.3,相位穿越频率Wg=5.55rad/s,幅值裕度 GM=15.3dB。用 MATL
29、AB 绘制校正后系统奈氏曲线的程序如下 :num=2;武汉理工大学自动控制原理课程设计说明书13den=conv(1 0,conv(0.1 1,0.5 1); %分母多项式展开W=tf(num,den); %开环传递函数kc=1; %稳态误差系数不变yPm=48.2+10; %增加量取10degW=tf(W); %超前校正环节mag,pha,w=bode(W*kc); %扩大系数后的开环频率特性的幅值和相位值Mag=20*log10(mag); %幅值的对数值的20倍Gm,Pm,Wcg,Wcp=margin(W*kc); %幅值稳定裕度Gm,相位稳定裕度Pm和相应的交接频率Wcg和Wcpphi
30、=(yPm-Pm)*pi/180; %确定m值alpha=(1+sin(phi)/(1-sin(phi); %确定a的值Mn=-10*log10(alpha); %a的对数值的负10倍Wcgn=spline(Mag,w,Mn); %确定最大相角位移频率T=1/Wcgn/sqrt(alpha); %求T值Tz=alpha*T;Wc=tf(Tz 1,T 1) %超前校正环节的传递函数Wx_c=feedback(W*kc*Wc,1) %校正后开环系统传递函数nyquist(W*kc*Wc,b); %校正后开环系统奈奎斯特图grid on;校正后系统的奈奎斯特曲线。图 6 校正后系统的奈奎斯特曲线-1
31、 -0.9 -0.8 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 0-15-10-5051015 0 dB-10 dB-6-4 d-2 dB10 dB6 dB4 2 dBNyquist DiagramReal AxisImaginaryAxis武汉理工大学自动控制原理课程设计说明书144 控制系统校正前后的对比4.1 控制系统校正前后的阶跃响应曲线、伯德图和奈奎斯特图对比采用扩大系数之后的校正前系统曲线与校正后进行对比。这样在控制 K 值相等的前提下,可以更好地看到超前校正对系统各项性能指标的影响。4.1.1 系统校正前后单位阶跃曲线对比及分析4.1.1.1 系统校正
32、前后单位阶跃响应曲线系统校正前后单位阶跃响应曲线。Step ResponseTime (seconds)Amplitude0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 500.20.40.60.811.21.4System: 前前前Rise time (seconds): 0.652System: 前前前Rise time (seconds): 0.766System: 前前前Peak amplitude: 1.16Overshoot (%): 16.2At time (seconds): 1.5System: 前前前Peak amplitud: 1.25Overshoot (%
33、): 25.4At time (seconds): 1.85System: 前前前Settling time (seconds): 3.29 System: 前前前Settling time (seconds): 4.4前前前前前前图 7 系统校正前后单位阶跃响应曲线虚线为校正前图像,实线为校正后图像。由图可知,校正前:上升时间 ,超调时间: ,调节时间str6.0stp85.14.4s,超调量 %=25.4% 。st校正后:上升时间 ,超调时间: ,调节时间 3.29s,超调量str652. stp5.st%=16.2%。武汉理工大学自动控制原理课程设计说明书154.1.2 系统校正前后伯德
34、图与奈奎斯特图对比4.1.2.1 系统校正前后伯德图系统校正前后伯德图。Bode DiagramFrequency (rad/s)-150-100-50050 System: 前前前Gain Margin (dB): 15.3At frequency (rad/s): 5.55Closed loop stable? YesSystem: 前前前Gain Margin (dB): 15.6At frequency (rad/s): 4.47Closed loop stable? YesMagnitude (dB)10-1 100 101 102 103-270-225-180-135-90Sys
35、tem: 前前前Phase Margin (deg): 43.2Delay Margin (sec): 0.484At frequency (rad/s): 1.56Closed loop stable? YesSystem: 前前前Phase Margin (deg): 51.3Delay Margin (sec): 0.479At frequency (rad/s): 1.87Closed loop stable? YesPhase (deg)前前前前前前图 8 系统校正前后伯德图虚线为校正前图像,实线为校正后图像。由图可知,校正前:截止频率 Wc=1.56rad/s,相位裕度 Pm=43
36、.2,相位穿越频率Wg=4.47rad/s,幅值裕度 GM=15.6dB。校正后:截止频率 Wc=1.87rad/s,相位裕度 Pm=51.3,相位穿越频率 Wg=5.55rad/s,幅值裕度 GM=15.3dB。武汉理工大学自动控制原理课程设计说明书163.1.2.2 系统校正前后奈奎斯特图系统校正前后奈奎斯特图。-1.4 -1.2 -1 -0.8 -0.6 -0.4 -0.2 0-20-15-10-505101520 0 dB-6 dB4 -2 dB6 dB4 dB2 dBNyquist DiagramReal AxisImaginary Axis前前前前前前图 9 系统校正前后奈奎斯特图
37、虚线为校正前图像,实线为校正后图像。4.1.2.3 频域角度分析超前校正对系统性能指标的影响从伯德图与奈奎斯特图中可以看出系统校正后,(1)截止频率增大,系统频带宽增加,增强了系统的快速性,但同时也造成了对高频干扰较敏感。(2)截止频率附近的对数幅频特性斜率减小,相位裕度增大,系统的相对稳定性增强。4.1.2.4 时域角度分析超前校正对系统性能指标的影响从系统校正前后的单位阶跃响应可以看出,(1)校正后系统的超调量减小,稳定裕度增加。(2)校正后系统的上升时间减少,调节时间减少,说明系统的快速性变好。武汉理工大学自动控制原理课程设计说明书175 频率法校正优缺点及适用条件5.1 频率法超前校正
38、的优缺点及适用条件5.1.1 频率法超前校正的优缺点:本方案采用超前校正来改善系统性能,从前面的分析设计中可以总结出其作用。优点:(1)在 附近使对数幅频特性斜率减小,增大系统相位裕度和增益裕度;c(2)截止频率增大,频带宽增加;(3)由于稳定裕度增加,单位阶跃响应的超调量减小;缺点:(1)由于频带加宽,对高频干扰较敏感;5.1.2 频率法超前校正的适用条件:适用范围:(1)在靠近 处,随 变化,相位滞后增加较缓慢的情况;c(2)要求有大的带宽和较快的动态响应的情况;(3)高频干扰不是主要问题的情况;不适用情况:(1)在 附近,相位滞后随 迅速增大的情况;c(2)相位超前要求过大;(3)因高频
39、干扰指标所限,不能增大高频增益的情况。5.2 频率法校正的其他情况频率法校正通常可分为三种基本的校正方式:串联校正、反馈校正和前馈校正。串联校正是最常用的校正方式,按串联装置的特点来分,又分为串联超前(微分)校正、串联滞后(积分)校正和串联滞后-超前(积分-微分)校正。本方案采用的超前校正,是用来提高系统的动态性能,而又不影响系统的稳态精度的一种校正方法。它是在系统中加入一个相位超前的校正装置,使之在穿越频率处相位超前,以增加相位裕度,这样既能使开环增益足够大,武汉理工大学自动控制原理课程设计说明书18又能提高系统的稳定性。滞后校正是在系统动态品质满意的情况下,为了改善系统稳态性能的一种校正方
40、法。从这种方法的频率特性上看,是在低频段提高其增益,而在穿越频率附近,保持其相位移的大小几乎不变。同时,如果稳态性能满足指标,而动态性能未满足要求,并希望降低频带宽度时,也可利用滞后校正来降低穿越频率。超前校正会使带宽增加,加快系统的动态响应速度,滞后校正可改善系统的稳态特性,减少稳态误差。如果需要同时改善系统的动态品质和稳态精度,则可采用串联滞后-超前校正。总之,每种方法的应用要根据系统的具体情况而定。5.3 频率法校正的优缺点及适用条件控制系统校正的方法有三种:时域法、根轨迹法和频率法。时域法的基本思想是在系统的传递函数基础上增加一定的校正环节,使校正后的传递函数的系统在性能上满足要求。根
41、轨迹法在系统中加入校正装置,就是加入了新的开环零,极点,这些新的零,极点将使校正后的闭环根轨迹,也就是闭环极点,向有利于改善系统性能的方向改变,这样可以做到使闭环零,极点重新布置,从而满足闭环系统的性能要求。频率法主要是应用开环 Bode 图,它的基本做法是利用恰当的校正装置,配合开环增益的调整,来修改原有的开环系统的 Bode 图,使得开环系统经校正于增益调整后的 Bode 图符合性能指标的要求。频率法校正的应用较为普遍,由于应用伯德图,使它具有清晰、直观、简便的优点。频率法通常以频率指标来衡量和调整控制系统的动态响应性能,而不是按时域指标,因此是一种间接的方法。在控制系统设计中,采用的设计
42、方法一般依据性能指标的形式而定.如果性能指标以单位阶跃响应的峰值时间,调节时间,超调量,阻尼比,稳态误差等时域特征量给出时,可采用时域法校正或根轨迹法校正;如果性能指标以系统的相角裕度,幅值裕度,谐振裕度,闭环带宽,稳态误差系数等频域特征量给出时,一般采用频率法校正.时域法与频域法是两种常用的方法,其性能指标可以互换。武汉理工大学自动控制原理课程设计说明书196 课程设计总结 这次课程设计,我选择的是超前校正网络的设计,通过这次课程设计我对课本知识又有了更深的理解,对校正过程中的静态速度误差系数,相角裕度,截止频率,超前(滞后)角频率,分度系数,时间常数等参数有了更深理解并应用到了设计当中。这
43、次设计中比较麻烦的就是超前网络参数的确定就是这个占用了大部分时间。在这参数确定以后,此次设计也就基本完事了。 在这次设计的过程中也有不少的失误出现,这对我的耐心也提出了考验。在已开始的时候,就没认真审题,把闭环 Nyquist 看成开环,在截取了图以后才发现。还有其他的一些错误。在不断的修改中,成果终于出现了。我们需要更多的思考,仔细想想所学知识的应用,各因素的相互制约关系,在工程计算中各性能指标的平衡等;我们更需要实践,实践不仅能够检验理论,更可以帮助我们发现更多的问题,我们必须在不断解决问题的过程中积累经验,在专业方面有所成长;当然还有不可或缺的耐心和坚持,是任何事成功的基础。在此我也感谢
44、在这次设计过程中帮助过我的老师和同学。武汉理工大学自动控制原理课程设计说明书207 参考文献 1 孟华.自动控制原理 .北京:机械工业出版社,2007 2 胡寿松.自动控制原理(第五版).北京:科学出版社,2007 3 薛定宇.控制系统仿真与计算机辅助设计.北京: 机械工业出版社,2006 4 刘卫国.MATLAB 程序设计与应用.北京:高等教育出版社,20075 朱衡君 .MATLAB 语言及实践教程(第二版).北京:清华大学出版社,20086 张静 .MATLAB 在控制系统中的应用 .北京:电子工业出版社,2007武汉理工大学自动控制原理课程设计说明书21附录用 MATLAB 绘制校正前
45、系统各曲线的程序num=2;den=conv(1 0,conv(0.1 1,0.5 1); %分母多项式展开W=tf(num,den); %开环传递函数Wy_c=feedback(W,1) %校正前开环系统传递函数figure(1);step(Wy_c,r,5); %校正前单位阶跃响应曲线grid on;figure(2);bode(W,r); %校正前开环系统伯德图grid on;figure(3);nyquist(W,r); %校正前开环系统奈奎斯特图grid on;武汉理工大学自动控制原理课程设计说明书222 利用 MATLAB 进行超前校正设计的程序num=2;den=conv(1 0
46、,conv(0.1 1,0.5 1); %分母多项式展开W=tf(num,den); %开环传递函数kc=1; %稳态误差系数不变yPm=48.2+10; %增加量取 10degW=tf(W); %超前校正环节mag,pha,w=bode(W*kc); %开环频率特性的幅值和相位值Mag=20*log10(mag); %幅值的对数值的 20 倍Gm,Pm,Wcg,Wcp=margin(W*kc); %幅值稳定裕度 Gm,相位稳定裕度 Pm 和相应的交接频率 Wcg 和 Wcpphi=(yPm-Pm)*pi/180; %确定 m 值alpha=(1+sin(phi)/(1-sin(phi); %
47、确定 a 的值Mn=-10*log10(alpha); %a 的对数值的负 10 倍Wcgn=spline(Mag,w,Mn); %确定最大相角位移频率T=1/Wcgn/sqrt(alpha); %求 T 值Tz=alpha*T;Wc=tf(Tz 1,T 1) %超前校正环节的传递函数Wy_c=feedback(W*kc,1) %校正前开环系统传递函数Wx_c=feedback(W*kc*Wc,1) %校正后开环系统传递函数figure(1);step(Wy_c,r,5); %校正前开环单位阶跃响应曲线hold on;step(Wx_c,b,5); %校正后开环单位阶跃响应曲线grid on;figure(2);bode(W*kc,r); %校正前开环系统伯德图hold on;bode(W*kc*Wc,b); %校正后开环系统伯德图grid on;武汉理工大学自动控制原理课程设计说明书23figure(3);nyquist(W*kc,r); %校正前开环系统奈奎斯特图hold on;nyquist(W*kc*Wc,b); %校正后开环系统奈奎斯特图grid on;