自控课设串联滞后超前校正要点

TOC\o"1-5"\h\z\o"CurrentDocument"绪论2\o"CurrentDocument"一课程设计的目的及题目3\o"CurrentDocument"1.1课程设计的目的3\o"CurrentDocument"1.2课程设计的题目3\o"CurrentDocument"二课程设计的任务及要求42.1课程设计的任务42.2课程设计的要求4\o"CurrentDocument"三校正函数的设计53.1理论知识5\o"CurrentDocument"3.2设计部分6\o"CurrentDocument"四传递函数特征根的计算104.1校正前系统的传递函数的特征根错误！未定义书签。4.2校正后系统的传递函数的特征根错误！未定义书签。\o"CurrentDocument"五系统动态性能的分析12\o"CurrentDocument"5.1校正前系统的动态性能分析125.2校正后系统的动态性能分析15六系统的根轨迹分析错误!未定义书签。6.1校正前系统的根轨迹分析错误！未定义书签。6.2校正后系统的根轨迹分析20\o"CurrentDocument"七系统的奈奎斯特曲线图21\o"CurrentDocument"7.1校正前系统的奈奎斯特曲线图217.2校正后系统的奈奎斯特曲线图错误！未定义书签。2八系统的对数幅频特性及对数相频特性错误！未定义书签。\o"CurrentDocument"8.1校正前系统的对数幅频特性及对数相频特性238.2校正后系统的对数幅频特性及对数相频特性错误！未定义书签。总结错误!未定义书签。6参考文献错误!未定义书签。绪论在控制工程中用得最广的是电气校正装置，它不但可应用于电的控制系统，而且通过将非电量信号转换成电量信号，还可应用于非电的控制系统。控制系统的设计问题常常可以归结为设计适当类型和适当参数值的校正装置。校正装置可以补偿系统不可变动部分（由控制对象、执行机构和量测部件组成的部分）在特性上的缺陷，使校正后的控制系统能满足事先要求的性能指标。常用的性能指标形式可以是时间域的指标，如上升时间、超调量、过渡过程时间等（见过渡过程），也可以是频率域的指标，如相角裕量、增益裕量（见相对稳定性）、谐振峰值、带宽（见频率响应）等。常用的串联校正装置有超前校正、滞后校正、滞后-超前校正三种类型。在许多情况下，它们都是由电阻、电容按不同方式连接成的一些四端网络。各类校正装置的特性可用它们的传递函数来表示，此外也常采用频率响应的波德图来表示。不同类型的校正装置对信号产生不同的校正作用，以满足不同要求的控制系统在改善特性上的需要。在工业控制系统如温度控制系统、流量控制系统中，串联校正装置采用有源网络的形式，并且制成通用性的调节器，称为PID（比例-积分-微分）调节器，它的校正作用与滞后-超前校正装置类同。-课程设计的目的及题目1.1课程设计的目的1、掌握自动控制原理的时域分析法，根轨迹法，频域分析法，以及各种补偿(校正)装置的作用及用法，能够利用不同的分析法对给定系统进行性能分析，能根据不同的系统性能指标要求进行合理的系统设计，并调试满足系统的指标。2、学会使用MATLAB语言及Simulink动态仿真工具进行系统仿真与调试。3、提高分析问题解决问题的能力。1.2课程设计的题目已知单位负反馈系统的开环传递函数G(S)=——00——，试用频率法设S(S+2)(S+40)计串联滞后一一超前校正装置，使系统的相角裕量y>40,静态速度误差系数K=20s-1。v二课程设计的任务及要求2.1课程设计的任务设计报告中，根据给定的性能指标选择合适的校正方式对原系统进行校正（须写清楚校正过程），使其满足工作要求。然后利用MATLAB对未校正系统和校正后系统的性能进行比较分析，针对每一问题分析时应写出程序，输出结果图和结论。最后还应写出心得体会与参考文献等。2.2课程设计的要求1、首先，根据给定的性能指标选择合适的校正方式对原系统进行校正，使其满足工作要求。要求程序执行的结果中有校正装置传递函数和校正后系统开环传递函数，校正装置的参数T，a等的值。2、利用MATLAB函数求出校正前与校正后系统的特征根，并判断其系统是否稳定，为什么？3、利用MATLAB作出系统校正前与校正后的单位脉冲响应曲线，单位阶跃响应曲线，单位斜坡响应曲线，分析这三种曲线的关系？求出系统校正前与校正后的动态性能指标。％、，、，、，以及稳态误差的值，并分析其有何变化？rps4、绘制系统校正前与校正后的根轨迹图，并求其分离点、汇合点及与虚轴交点的坐标和相应点的增益K*值，得出系统稳定时增益K*的变化范围。绘制系统校正前与校正后的Nyquist图，判断系统的稳定性，并说明理由？5、绘制系统校正前与校正后的Bode图，计算系统的幅值裕量，相位裕量，幅值穿越频率和相位穿越频率。判断系统的稳定性，并说明理由？三校正函数的设计要求：首先，根据给定的性能指标选择合适的校正方式对原系统进行校正，使其满足工作要求。要求程序执行的结果中有校正装置传递函数和校正后系统开环传递函数，校正装置的参数T，a等的值。3.1校正步骤应用频率特性法设计串联滞后-超前校正装置的步骤如下：根据稳态误差要求，确定开环增益K。利用已确定的开环增益，作出未校正系统的对数频率特性曲线，确定未校正系统的剪切频率①°。，相角裕度Y0和幅值裕度h以检验性能指标是否满足要求。若不满足要求，则执行下一步。确定滞后校正器传递函数的参数：1+bTs1—1+Ts1式中b<1,一要距①较远为好。Tc11选择一个新的系统剪切频率o，使在这一点超前校正所提供的相位超前c2量达到系统相位稳定裕量的要求。又要使得在这一点原系统加上滞后校正器综合幅频特性衰减为0dB，即L曲线在o点穿越横坐标。c2确定超前校正器传递函数的参数G(s)=1+aT2S

c21+Ts2式中a>1。由以下表达式:20lga=-L(W2)L(o)为原系统加上滞后校正器后幅频分贝值。还有公式c2T=—求出参数a、Tvaom

画出校正后的系统的Bode图，并验算己校正系统相角裕度和幅值裕度。3.2设计部分已知单位负反馈系统的开环传递函数G(S)=匕，试用频率法设S(S+2)(S+40)计串联滞后一一超前校正装置，使系统的相角裕量y>40,静态速度误差系数K=20s-io⑴由静态速度误差系数K=20s-io可知K=1600oV0所以，该系统的开环传递函数为G&)=1600

s(v+2)G+40)化为：G(v)=v，KA所以，该系统的开环传递函数为G&)=1600

s(v+2)G+40)s(O.5s+lH0.O25s+1)(2)确定未校正系统的剪切频率co,相角裕度y八和幅值裕度//。cQ0MATLAB程序如下：»dl=[l,0];»d2=[l,2];»d3=[l,40];»denl=conv(dl,d2);>>den=conv(deni,d3);>>num=[1600]:»gO=tf(num,den);»[gm,pm,weg,wep]=margin(gO);>>margin(gO);>>figure(1)校正前的Bode图为:BodeDiagramGm=6.44dB(at&4rad/sec).Pm=&.S5d&g(at6.13rati/s&c)EOOO1G-1w1苻Frequencyfract/sec)浏25-1-Z(BlllPJlllwELICLID幅值裕度：h=6.44dB相角裕度：yo=9.35deg穿越频率：①=8.9443rad/sec截止频率：①=6.13rad/sec（3）计算出滞后校正器的传递函数。1+bTs1+Ts因为1=f1〜1\，可以取土=-x6.13=1.226,bT=0.82。bT"510)cbT5由于b<1，取b=0.5，则T=1.64，得到滞后校正补偿器传递函数:气（s）-0.82s+11.64s+1串联滞后校正器的系统传递函数为:G(s)•G(s)=cl200.82s+1s(0.5s+1)6.025s+1).1.64s+1由MATLAB语言计算出串联滞后校正后的Y]=11。G(s)•G(s)=cl>>n1=20;>>d1=conv(conv([10],[0.51]),[0.0251]);>>s1=tf(n1,d1);>>s2=tf([0.821],[1.641]);>>sope=s1*s2;>>margin(sope)(4)用MATLAB语言计算出超前校正器的传递函数。G(s)=1+aT2s，以〉1

c21+T2s题目要求伫40。，取丫=41。，所以最大超前相角：甲=1一匕+(5。〜12。)=41。—11。+11。=41。则校正装置的参数：以=1+sin气=4.8151—sin中m用MATLAB语言计算出o=6.4756，T=0.07:>>n1=20;>>d1=conv(conv([10],[0.51]),[0.0251]);>>g=tf(n1,d1);>>s2=tf([0.821],[1.641]);>>sope=g*s2;>>v=41;>>v1=v-11+11;>>v1=v1*pi/180;>>a=(1+sin(v1))/(1-sin(v1));>>k=10*log10(a);>>[mag,phase,w]=bode(sope);>>kdB=20*log10(mag);>>wm=spline(kdB,w,-k);>>T=1/(wm^sqrt(a));1f1CD=CD=———，T=|c2mjar<'acom解方程组得：a=4.815,T=0.07,得到超前校正器的传递函数:2GGk1+03375

c21+0.07s(5)校验系统校正后系统是否满足题目要求。校正后的传递函数为：()16000.82^+10.3375+1GXsJCj\sJCj\sJ=—fw\cic2++4071.64s+10.07s+1用MATLAB语言校正如下：>>nl=1600;»dl=conv(conv([l0],[12]),[140]):»sl=tf(nl,dl);»s2=tf([0.821],[1.641]);»s3=tfCO.3371],[0.071]);»sope二sl*s2*s3;>>margin(sope)校正后的Bode图为：1Q0BodeDiagramGm=16.7dIB{at22.1rad/sec),Pm=43.7deg(at6.46rad/sec)O5O5Q■-OO-1o1Id1Q1Q1d2w3Frequencyfrad/sec)oO59sc

-(5IDPI)ciibmjzcl幅值裕度：h=16.7dB相角裕度：y0=43.7deg，满足题目要求y>40截止频率：oc=6.46rad/sec穿越频率：o=22.1rad/sec四传递函数特征根的计算要求：利用MATLAB函数求出校正前与校正后系统的特征根，并判断其系统是否稳定，为什么？4.1校正前系统的传递函数的特征根校正前的开环传递函数为:G(G(s)=1600

s(s+2)S+40)MATLAB程序如下:>>num=1600;>>den=conv(conv([10],[12]),[140]):>>g=tf(num,den);>>sys=feedback(g,1);»pzmap(g);»[p,z]=pzmap(g):>>den=sys.den{1};>>r=roots(den);»disp(r)-41.0006-0.4997+6.2269i-0.4997-6.2269i系统没有零极点在右边，所以系统开环稳定。4.2校正后系统的传递函数的特征根校正后的开环传递函数为：()16000.82s+10.337s+1G\s)G\sJCj\sJ=—fw\cic2+27^+4071.64s+10.07s+1MATLAB程序如下：»k=1600;>>num=conv([0.821],[0.3371]);>>den=conv(conv(conv(conv([10],[12]),[140]),[1.641]),[0.071]):>>g=tf(k*num,den);>>sys=feedback(g,1);»pzmap(g);»[p,z]=pzmap(g):>>den=sys.den{1};>>r=roots(den);»disp(r)

-43.0016-4.4129+6.9169i-4.4129-6.9169i-3.8016-1.2665系统没有零极点在右边，所以系统开环稳定。五系统动态性能的分析要求：利用MATLAB作出系统校正前与校正后的单位脉冲响应曲线，单位阶跃响应曲线，单位斜坡响应曲线，分析这三种曲线的关系？求出系统校正前与校正后的动态性能指标a苏，、，、，以及稳态误差的值，并分析其有何变化？rps5.1校正前系统的动态性能分析校正前的开环传递函数为：G(G(s)=1600

s(s+2)S+40)单位脉冲响应MATLAB程序为：>>num=1600;>>den=conv(conv([1,0],[1,2]),[1,40]);>>G1=tf(num,den);>>sys1=feedback(G1,1);>>impulse(sys1);

校正前单位脉冲响应曲线单位阶跃响应MATLAB程序为：>>num=1600;>>den=conv(conv([1,0],[1,2]),[1,40]);>>G1=tf(num,den);>>sys1=feedback(G1,1);>>step(sys1);1.3校正前单位阶跃响应曲线由图可知：tp=0.544(3)单位斜坡响应t=0.178,t=7.67,。％=76.5%,e=1-1=0MATLAB程序为:>>n1=1600;>>d1=conv(conv([10],[12]),[140]);>>sope1=tf(n1,d1);>>sys1=feedback(sope1,1);>>n2=1;>>d2=conv([10],[11]);>>rs=tf(n2,d2);>>sys2二sys1*rs;>>step(sys2)B10121416B101214161820Time(sec)OJpnwdLUV----------Moos642os642o2111115.2校正后系统的动态性能分析校正后的开环传递函数为：八""C16000.82s+10.337s+1GsJGsGs)=v\C1c2s\s+2)s+40)1.64s+10.07s+1(1)单位脉冲响应MATLAB程序为：>>k=1600;>>num=conv([0.821],[0.3371]);>>den=conv(conv(conv(conv([10],[12]),[140]),[1.641]),[0.071]);>>g=tf(k*num,den);>>sys=feedback(g,1);>>impulse(sys);

-2Q.511Time(sec)校正后单位脉冲响应曲线-2Q.511Time(sec)(2)单位阶跃响应MATLAB程序为：>>k=1600;>>num=conv([0.821],[0.3371]);>>den=conv(conv(conv(conv([10],[12]),[140]),[1.641]),[0.071]);>>g=tf(k*num,den);>>sys=feedback(g,1);>>step(sys)

OEIJQ.*.s.64o,o,o.,2oo.OEIJQ.*.s.64o,o,o.,2oo.校正后单位阶跃响应曲线由图可知：成0.438,0.169，八0.816,°%=29.9%,七=1-1=。单位斜坡响应MATLAB程序为:>>k=1600;>>num=conv([0.821],[0.3371]);>>den=conv(conv(conv(conv([10],[12]),[140]),[1.641]),[0.071]);>>g=tf(k*num,den);>>sys1=feedback(g,1);>>n2=1;>>d2=conv([10],[11]);>>rs=tf(n2,d2);>>sys2=sys1*rs;>>step(sys2)

tD*-d与校正后单位斜坡响应曲线tD*-d与三条曲线关系：斜坡响应曲线的导数是阶跃，阶跃响应曲线的导数是脉冲响应曲线。校正后的t、t、t、q%与校正前均减小，e不变，都是0。prsss六系统的根轨迹分析要求：绘制系统校正前与校正后的根轨迹图，并求其分离点、汇合点及与虚轴交点的坐标和相应点的增益K*值，得出系统稳定时增益K*的变化范围。6.1校正前系统的根轨迹分析校正前的开环传递函数为：G(s)G(s)=1600

s(s+2)S+40)num=[1600];den二conv(conv([10],[12]),[140]);rlocus(num,den)[k,poles]=rlocfind(num,den)得到校正前系统的根轨迹：RootLocusSO6040oo-2w_xv胃LnBEE-20oo-soSyst&m:sysGain:1.9SPole:0.000152+3.7iDamping:-1.74e-QQ5OvershootSO6040oo-2w_xv胃LnBEE-20oo-soSyst&m:sysGain:1.9SPole:0.000152+3.7iDamping:-1.74e-QQ5Overshoot(%):1Q0Frequency(rad/sec):S.7-6Q-40-2dRealAxisQ2Q4a0+8.6957ik=1.985poles:-41.9-0.050111+8.7062i-0.050111-8.7062iselected_point=0-8.1988ik=1.7657poles=-41.7060

-0.1470+8.2291i-0.1470-8.2291i分离点：d=—1.52K*=0.0176与虚轴交点：0.000152土8.7jK*=1.986.2校正后系统的根轨迹分析校正后的开环传递函数为：八""C16000.82s+10.337s+1GsJGsGs)=wC1c2s\s+2)s+40)1.64s+10.07s+1MATLAB程序如下：>>n1=conv(conv([01600],[0.821]),[0.3371]);>>den=conv(conv(conv(conv([10],[12]),[140]),[1.641]),[0.071]);>>g=tf(n1,den);>>rlocus(g)>>[k,poles]=rlocfind(g)得到校正后系统的根轨迹：RootLocus604020|»_点^r-rac-araE-System:gGain:InfPole:-2.97Damping:1Overshoot(%.):0Frequency(rad/sec):2.97System:gGain:6.69Pole:-0.00695+21.SDamping:O.OdQSISOvershoot(%):Frequency(rad/sec):21.S-20OO80■-6Q-40-2QRealAxisSystem:gGain:0.00946Pole:-0342Damping:1Overshoot(%):0Frequency(rad/sec):0.342Q2Q604020|»_点^r-rac-araE-System:gGain:InfPole:-2.97Damping:1Overshoot(%.):0Frequency(rad/sec):2.97System:gGain:6.69Pole:-0.00695+21.SDamping:O.OdQSISOvershoot(%):Frequency(rad/sec):21.S-20OO80■-6Q-40-2QRealAxisSystem:gGain:0.00946Pole:-0342Damping:1Overshoot(%):0Frequency(rad/sec):0.342Q2Q40要求：绘制系统校正前与校正后的Nyquist图，判断系统的稳定性，并说明理由。7.1校正前系统的奈奎斯特曲线图校正前的开环传递函数为:G(s)G(s)=1600

s(s+2)S+40)>>num=1600;>>den=conv(conv([10],[12]),[140]);>>g=tf(num,den);>>nyquist(g)

o2050505Q

■-oo1-50-12Q-1-2QNyquistDiagram-6RealAxis-4-20o2050505Q

■-oo1-50-12Q-1-2QNyquistDiagram-6RealAxis-4-20o7.2校正后系统的奈奎斯特曲线图校正后的开环传递函数为:G(sG(sG(s)=

clc216000.82s+10.337s+1s(s+2)S+40)1.64s+10.07s+1MATLABG(sG(sG(s)=

clc2>>k=1600;>>num=conv([0.821],[0.3371]);>>den=conv(conv(conv(conv([10],[12]),[140]),[1.641]),[0.071]);>>g=tf(k*num,den);>>nyquist(g)

o40ooooow_>v^r-racaraE-oo-=-2o-4NyquistCiagramRealAxisQ-2a由图可知N=0o40ooooow_>v^r-racaraE-oo-=-2o-4NyquistCiagramRealAxisQ-2a第八章系统的对数幅频特性及对数相频特性要求：绘制系统校正前与校正后的Bode图，计算系统的幅值裕量，相位裕量，幅值穿越频率和相位穿越频率。判断系统的稳定性，并说明理由。8.1校正前系统的对数幅频特性及对数相频特性校正前的开环传递函数为:G(s)=1600

s(s+2G(s)=1600

s(s+2)S+40)>>d1=[1,0];>>d2=[1,2];>>d3=[1,40];>>den1=conv(d1,d2);>>den=conv(den1,d3);

>>num=[1600];>>g0=tf(num,den);>>[gm,pm,wcg,wcp]=margin(g0);>>margin(g0);>>figure(1)ffisBodeDiagramGm=ffisBodeDiagramGm=6.44dB(atS.Mrad/sec),Pm=9.35deg(at6.13rati/sec)50oo-5Q1Gw11D",蓿1Gw11D",蓿Frequencyfract/sec)52-2〔BillsillwELICL-11D幅值裕