自动控制课程设计报告书

1设计目的…………………………………………………………22设计容与条件…………………………………………………22.1设计容……………………………………………………22.2设计条件……………………………………………………23滞后校正特性及设计一般步骤…………………………………23.1滞后特性校正………………………………………………23.2滞后校正设计一般步骤……………………………………34校正系统分析……………………………………………………34.1校正参数确定………………………………………………34.2校正前后系统特征根及图像………………………………64.3函数动态性能指标及其图像………………………………104.4系统校正前后根轨迹及其图像……………………………114.5Nyquist图…………………………………………………124.6Bode图………………………………………………………155设计心得体会……………………………………………………176设计主要参考文献………………………………………………18串联滞后校正装置设计1、设计目的：1)了解控制系统设计的一般方法、步骤。2)掌握对系统进行稳定性分析、稳态误差分析以及动态特性分析的方法。3)掌握利用MATLAB对控制理论容进行分析和研究的技能。4)提高分析问题解决问题的能力。2、设计容与条件：2.1设计容：1)阅读有关资料。2)对系统进行稳定性分析、稳态误差分析以及动态特性分析。3)绘制根轨迹图、Bode图、Nyquist图。4)设计校正系统，满足工作要求。2.2设计条件：已知单位负反馈系统的开环传递函数0KG(S)S(0.0625S1)(0.2S1)，试用频率法设计串联滞后校正装置，使系统的相位裕度050，静态速度误差系数1vK40s，增益欲度17dB。3、滞后校正特性及设计一般步骤：3.1滞后特性校正：滞后校正就是在前向通道中串联传递函数为)(sGc的校正装置来校正控制系统，)(sGc的表达式如下所示。1,11)(aTsaTssGc其中，参数a、T可调。滞后校正的高频段是负增益，因此，滞后校正对系统中高频噪声有削弱作用，增强了抗干扰能力。可以利用滞后校正的这一低通滤波所造成的高频衰减特性，降低系统的截止频率，提高系统的相位裕度，以改善系统的暂态性能。滞后校正的基本原理是利用滞后网络的高频幅值衰减特性使系统截止频率下降，从而使系统获得足够的相位裕度。或者，是利用滞后网络的低通滤波特性，使低频信号有较高的增益，从而提高了系统的稳态精度。可以说，滞后校正在保持暂态性能不变的基础上，提高开环增益。也可以等价地说滞后校正可以补偿因开环增益提高而发生的暂态性能的变化。3.2滞后校正设计一般步骤：1）根据稳态误差要求，确定开环增益K。2）利用确定好的开环增益K，绘制未经校正系统的伯德图，确定未校正系的幅值穿越频率1c，相位裕量和增益裕量GM。3）根据相角裕度要求'，确定已校正系统的截止频率2c。在2c处，原系统的相位为12~5180'2c，其中，'为期望的相角裕量，12~5为附加滞后校正环节，在2c处引起的相角滞后。4）根据下式确定滞后环节的参数b和T0lg20'cLb（1）'1.01cbT（2）利用（1）式确定的b，由式（2）计算出滞后环节的参数T。如果参数T过大，难以实现，可将式（2）中的0.1适当加大，通常在0.1~0.25之间选取，而2c的估计值相应地在14~6围确定。5）校验校正后系统的相角裕度和幅值裕度。如果矫正后的系统满足要求，校正工作结束；如果校正后的系统不能满足指标要求，跳回第2步，重复进行以后的步骤，直到满足给出的指标要求为止。6）根据滞后校正环节的参数，确定滞后校正环节的实现。可以用硬件实现，设计滞后校正网络；也可以考虑用软件实现编制相应的程序算法。4、校正系统分析：4.1校正参数确定：已知单位负反馈系统的开环传递函数G（s）=)12.0)(10625.0(0SSSK，试用频率法设计串联滞后校正装置，使系统的相位裕度050，静态速度误差系数140sKv，增益欲度dBKg17lg20。计算如下：（1）40)12.0)(10625.0(lim000SSSSKKKsv40)102.0)(10625.0(lim00SSKs4010K400KK所以校正前开环传递函数为)102.0)(10625.0(40)(SSSsG（2）由1)2.0(1)0625.0(4022ccccA可得13.121crad/s,相角裕度78.14。（3）由50'，取附加补偿角为-62.0arctan0625.0arctan90令180，180902.0arctan0625.0arctan'''2.0arctan0625.0arctan90cccjG5014~60'56180'cjGsradc/4.2'(4)滞后装置的传递函数)1()1(TSbTssG，'101~511cbT-20lgb=20lg4.240,所以b=0.053'101~511cbT4.2'cT=39.3~78.61，取T=48所以滞后校正传递函数为)481()2.5441()(SSsG矫正后传递函数为)481()544.21()12.0)(10625.0(40)(SSSSSsG（5）、用MATLAB求出系统校正后的幅值裕度和相角裕度,并画出Bode图：在MATLAB中输入：得到的结果是Gm=8.9369Pm=50.6177Wcp=8.4961Wcg=1.9919由以上数据可以看出，在串联了一个滞后校正装置后系统稳定,γ﹦50.6177°，满足0503、利用MATLAB函数求系统校正前与校正后的特征根，并判断系统是否稳定，为什么？校正前：num=[40];den=[0.01250.262510]g=tf(num,den);sys=feedback(g,1);pzmap(g);den=[0.01250.2625140]t=tf(num,den);pzmap(t);[p,z]=pzmap(g);den=sys.den{1};r=roots(den);disp(r)语段执行结果如下：-23.41871.2094+11.6267i1.2094-11.6267i开环零极点分布图：由开环零极点分布图可知，系统有极点在右半平面，故系统开环不稳定。校正后：程序如下所示：Clearall;num=40;den=conv([0.062510],[0.21]);Gs=tf(num,den);n1=[2.5441];d1=[481];Gc=tf(n1,d1);sope=Gs*Gc;Go=feedback(sope,1,-1)eig(Go)语段执行结果如下：Transferfunction:101.8s+40----------------------------------------------0.6s^4+12.61s^3+48.26s^2+102.8s+40eig(Go)ans=-16.8327-1.8507+2.1709i-1.8507-2.1709i-0.4867计算数据表明，特征根中无实部为正的根，所以闭环系统是稳定的。4、利用MATLAB作出系统校正前与校正后的单位脉冲响应曲线，单位阶跃响应曲线，单位斜坡响应的曲线图，分析三种曲线关系。求出系统动态性能指标的值，以及稳态误差的值，并分析其有何变化？校正前：程序如下：num=40;den=conv([0.062510],[0.21]);g=tf(num,den);sys=feedback(g,1,-1);[y,t]=step(sys);C=dcgain(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.98*C&&y(s)1.02*Cs=s-1;endts=t(s)figure(1);step(sys);figure(2);impulse(sys)figure(3);[num2,den2]=tfdata(sys,'v');step(num2,[den2,0])ess=1-dcgain(sys)结果如下：tp=1.9421max_overshoot=876.8570tr=0.1142ts=1.9801ess=0校正后：程序如下所示num=40;den=conv([0.062510],[0.21]);Gs=tf(num,den);n1=[2.5441];d1=[481];Gc=tf(n1,d1);sope=Gs*Gc;Go=feedback(sope,1,-1)[y,t]=step(Go);C=dcgain(Go);[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.98*C&&y(s)1.02*Cs=s-1;endts=t(s)figure(1);step(Go)figure(2);impulse(Go)figure(3);[num2,den2]=tfdata(Go,'v');step(num2,[den2,0])ess=1-dcgain(Go)Transferfunction:101.8s+40----------------------------------------------0.6s^4+12.61s^3+48.26s^2+102.8s+40tp=1.4405max_overshoot=23.6780tr=0.5602ts=5.5220ess=0校正前阶跃响应校正后阶跃响应校正前脉冲响应校正后脉冲响应校正前斜坡响应校正后斜坡响应三条曲线关系：斜坡响应曲线的导数是阶跃，阶跃响应曲线的导数是冲激响应曲线。%，st,pt的值都变小，rt的值变大，sse的值不变。5.绘制系统的根轨迹图，并求出其分离点、汇合点及与虚轴交点的坐标和相应点的增益K*值，得出系统稳定时增益K*的变化围。校正前：程序如下所示clearn=1;z=[];p=[0,-3,-5];[num,den]=zp2tf(z,p,n);rlocus(num,den);title('rootlocus');校正前系统的根轨迹图分离点d=-1.24,K=K*=0.205与虚轴交点±3.76i,K=K*=3.2无汇合点校正后：程序如下所示：clear;num=1;den=conv([0.510],[0.1671]);Gs=tf(num,den);n1=[8.41];d1=[35.011];Gc=tf(n1,d1);G=Gs*Gc;k=0:0.05:200;%给定k的围rlocus(G,k)%绘制给定k的围下的跟轨迹[k,poles]=rlocfind(G)%交互式的选取根轨迹上的增益，这里用于选取其临界稳定值Selectapointinthegraphicswindowselected_point=-0.0071-3.3043ik=30.3047poles=-7.8579-0.0191+3.3094i-0.0191-3.3094i-0.1204校正后的根轨迹图：分离点d=-2.29,K=K*=0.6与虚