自动控制原理课程设计__频率法设计串联滞后超前校正装置

目录一、设计任务.......................................................1二、设计要求........................................................1三、设计步骤.........................................................11.未校正前系统的性能分析.........................................11.1开环增益0K...........................................11.2校正前系统的各种波形图................................21.3由图可知校正前系统的频域性能指标......................41.4特征根................................................41.5判断系统稳定性........................................51.6分析三种曲线的关系....................................51.7求出系统校正前动态性能指标及稳态误差..................51.8绘制系统校正前的根轨迹图..............................51.9绘制系统校正前的Nyquist图............................62.校正后的系统的性能分析.........................................72.1滞后超前校正...........................................72.2校正后系统的各种波形图................................82.3由图可知校正后系统的频域性能指标......................102.4特征根................................................112.5判断系统稳定性........................................112.6分析三种曲线的关系....................................112.7求出系统校正后动态性能指标及稳态误差..................112.8绘制系统校正后的根轨迹图和Nyquist图..................11四、心得体会.......................................................13五、主要参考文献....................................................131一、设计任务已知单位负反馈系统的开环传递函数0()(0.11)(0.011)KGSSSS，试用频率法设计串联滞后——超前校正装置。（1）使系统的相位裕度045（2）静态速度误差系数250/vKrads（3）幅值穿越频率30/Crads二、设计要求（1）首先，根据给定的性能指标选择合适的校正方式对原系统进行校正，使其满足工作要求。要求程序执行的结果中有校正装置传递函数和校正后系统开环传递函数，校正装置的参数T，等的值。（2）利用MATLAB函数求出校正前与校正后系统的特征根，并判断其系统是否稳定，分析原因。（3）利用MATLAB作出系统校正前与校正后的单位脉冲响应曲线，单位阶跃响应曲线，单位斜坡响应曲线，分析这三种曲线的关系，求出系统校正前与校正后的动态性能指标σ%、tr、tp、ts以及稳态误差的值，并分析其有何变化。（4）绘制系统校正前与校正后的根轨迹图，并求其分离点、汇合点及与虚轴交点的坐标和相应点的增益K值，得出系统稳定时增益K的变化范围。绘制系统校正前与校正后的Nyquist图，判断系统的稳定性，并说明理由。（5）绘制系统校正前与校正后的Bode图，计算系统的幅值裕量，相位裕量，幅值穿越频率和相位穿越频率。判断系统的稳定性，并说明理由。三、设计步骤开环传递函数0()(0.11)(0.011)KGSSSS1、未校正前系统的性能分析1.1开环增益2已知系统中只有一个积分环节，所以属于I型系统由静态速度误差系数250/vKrads可选取vK=600rad/ssradKSSSKSSHSSGKssV/600)101.0)(11.0(lim)()(lim0000（1—1）开环传递函数为)101.0)(11.0(600)(SSSSG1.2通过MATLAB绘制出校正前系统的bode图和校正前系统的单位阶跃响应图MATALAB程序为：cleark=600;n1=1;d1=conv(conv([10],[0.11]),[0.011]);s1=tf(k*n1,d1);figure(1);sys=feedback(s1,1);step(sys)c=dcgain(sys);[y,t]=step(sys);[max_y,k]=max(y);tp=t(k)max_overshoot=100*(max_y-c)/cr1=1;while(y(r1)0.1*c)r1=r1+1;endr2=1;while(y(r2)0.9*c)r2=r2+1;endtr=t(r2)-t(r1)s=length(t);whiley(s)0.98c&&y(s)1.02*cs=s-1;3endts=t(s)figure(2);margin(s1);holdonfigure(3);sys=feedback(s1,1);impulse(sys)figure(4);step(k*n1,[d1,0])ess=1-dcgain(sys)图1-1校正前系统的bode图图1-2校正前系统的单位阶跃响应4图1-3校正前系统的单位脉冲响应图1-4校正前系统的单位斜坡响应1.3由图得校正前系统的频域性能指标幅值裕度hL=-14.7dB；穿越频率x=31.6/rads;相角裕度r=-26.8度；剪切频率c=69.8/rads。1.4特征根51.5判断系统稳定性（1）由图1可以看出，0c之前对数幅频渐近特性曲线所对应的相频特性曲线穿越了。（2）由特征根可以看出，有根在右半平面，因此系统不稳定，按本题要求，需要进行串联滞后超前校正。1.6分析校正前单位脉冲响应曲线，单位阶跃响应曲线，单位斜坡响应曲线，这三种曲线的关系单位斜坡响应的一次导数是阶跃响应曲线，阶跃响应的一次导数是冲击响应。1.7求出系统校正前动态性能指标σ%、tr、tp、ts以及稳态误差的值%max_overshoot=560.9737tr=0.0145tp=0.1489ts=0.1998ess=01.8绘制系统校正前的根轨迹图，并求其分离点、汇合点及与虚轴交点的坐标和相应点的增益K值，得出系统稳定时增益K的变化范围。程序：cleark=600;n1=1;d1=conv(conv([10],[0.11]),[0.011]);s1=tf(k*n1,d1);k=0:0.05:200;figure(1);rlocus(s1,k)figure(2);nyquist(s1)[k,poles]=rlocfind(s1)Selectapointinthegraphicswindow6selected_point=-97.7488-1.5528ik=0.0391poles=1.0e+002*-1.0248-0.0376+0.1466i-0.0376-0.1466i图1-5校正前系统的根轨迹分离点-10与虚轴的交点为-97.74881.5528i，当取-97.7488-1.5528i点时，k=(00.0391)1.9绘制系统校正前的Nyquist图，判断系统的稳定性，并说明理由。图1-6校正前系统的耐奎斯特曲线7因为系统的耐奎斯特曲线顺时针包围（-1，j0）点1圈，所以R=-1，没有实部为正的极点所以P=0，Z=P-R=1，闭环系统不稳定。2、校正后的系统的性能分析2.1滞后超前校正题目中要求30/Crads,取C=35rad/s，过C处作一斜率为-20dB/dec的直线作为期望特性的中频段。图2-1为使校正后系统的开环增益不低于250rad/s，期望特性的低频段应与未校正系统特性一致。而未校正系统的低频段斜率与期望特性的中频段斜率同为-20dB/dec，即两线段平行，为此，需在期望特性的中频段与低频段之间用一斜率为-40dB/dec的直线作连接线。连接线与中频段特性相交的转折频率2距C不宜太近，否则难于保证系统相角裕度的要求。现按10~52CC的原则选取75352为使校正装置不过于复杂，期望特性的高频段与未校正系统特性一致。由于未校正系统高频段特性的斜率是-60dB/dec，故期望特性中频段与高频段之间也应有斜率为-40dB/dec的直线作为连接线。用未校正系统的特性Lo减去期望特性，就得到串联校正装置的对数幅频特性Lc,-20db/dec-40db/dec-40db/dec100Wc0=69.810-20db/dec-40db/decWc=35-60db/decW2=8-20db/dec+20db/dec8它表明，应在系统中串联相位滞后-超前校正装置。其传递函数为:STSTSaTSbTSGC12121111)(（1—2）式中：143.071122bT1.01011aTba1由上图可以写出aKlg20)lg(lg20lg2010351015.17a058.0b47.22T006.01T因此，串联滞后-超前校正装置的传递函数为SSSSSGC006.0147.211.01143.01校正后系统的开环传递函数为101.01006.0147.21143.0600)()(0SSSSSSGSGSGC（1—3）2.2通过MATLAB绘制出校正后系统的bode图和校正后系统的单位阶跃响应图MATALAB程序为：clearn1=600;d1=conv(conv([10],[0.11]),[0.011]);s1=tf(n1,d1);s2=tf([0.1431],[2.471]);s3=tf([0.11],[0.0061]);sope=s1*s2*s3;figure(1);margin(sope);holdonfigure(2);sys=feedback(sope,1);step(sys)[y,t]=step(sys);c=dcgain(sys);[max_y,k]=max(y);9tp=t(k)max_overshoot=100*(max_y-c)/cr1=1;while(y(r1)0.1*c)r1=r1+1;endr2=1;while(y(r2)0.9*c)r2=r2+1;endtr=t(r2)-t(r1)s=length(t);whiley(s)0.98*c&&y(s)1.02*cs=s-1;endts=t(s)figure(3);sys=feedback(s1,1);impulse(sys)figure(4);step(k*n1,[d1,0])ess=1-dcgain(sys)图2-2校正后系统的bode图10图2-3校正后系统的单位阶跃响应图图2-4校正后系统的单位脉冲响应图图2-5校正后系统的单位斜坡响应图2.3由图得校正前系统的频域性能指标幅值裕度hL=16.7dB；穿越频率x=122/rads