chapter7_11.m

《MATLAB在控制系统中的应用》实例程序

ng0=[1];dg0=conv([1,0],conv([1,1],[1,2])); 
t=[0:0.01:30];w=logspace(-2,2); 
KK=20;delta=0.05; 
g0=tf(ng0,dg0);b1=feedback(g0,1); 
[gm,pm,wcg,wcp]=margin(KK*g0) 
[pos,tr,ts,tp]=stepchar(b1,delta) 
Pm=45;wc=1.5+0.2; 
[ng1,dg1]=lead5(ng0,dg0,1,Pm,wc,w); 
gc1=tf(ng1,dg1) 
g01=g0*gc1; 
[ng2,dg2]=lag2(g01.num{1},g01.den{1},w,KK,Pm); 
gc2=tf(ng2,dg2) 
g02=KK*g01*gc2 
b2=feedback(g02,1); 
step(b1,t);grid on,hold on;step(b2,t);hold off; 
figure,bode(KK*g0,w);grid on,hold on;bode(g02,w);hold off; 
[gm,pm,wcg,wcp]=margin(g02) 
[pos,tr,ts,tp]=stepchar(b2,delta)