MATLAB 校正系统程序.docx

clear; G1=tf(20,[1,0.4,0]) %系统的开环传递函数 G=feedback(G1,1,-1) %原系统的闭环传递函数 figure(1) step(G) %原系统的单位阶跃响应 figure(2) rlocus(G1) %原系统的根轨迹 figure(3) margin(G1) %原系统的bode图 [Gm0,Pm0,Wg0,Wc0]=margin(G1) deta=5 %补偿角 Phim=45-Pm0+deta %校正装置最大超前角 a=sin(Phim*pi/180) alpha=(1+a)/(1-a) Wm=Wc0*alpha^(1/4) T=1/(Wm*sqrt(alpha)); alpha_T=alpha*T Gc=tf([alpha_T,1],[T,1]) %校正装置 Ga=G1*Gc %校正后的传递函数 G11=feedback(G1,1); %原系统闭环传递函数 Gaa=feedback(Ga,1); %校正后系统闭环传递函数 figure(4) rlocus(Ga)%校正后系统的根轨迹 legend(Ga) figure(5) margin(Ga) %校正后系统的频域特性 [Gm,Pm,Wcm,Wpm]=margin(Ga) figure(6) step(G11,-.k,Gaa)%原系统与校正后系统的时域特性的比较 legend(G1,