ex6_3.m

反馈控制系统设计与分析（matlab程序）

G=tf(100,[0.04,1,0]); wc=100; gamma=[30,40,50,60,70,80,90];
f1=figure; f2=figure; m=[]; p=[]; y=[];
w=logspace(-1,3); t=0:0.005:0.2;
for i=1:length(gamma)
   Gc=leadlagc(G,wc,gamma(i),10);
   G_o=Gc*G; G_c=feedback(G_o,1);
   [m1,p1]=bode(G_o,w); y1=step(G_c,t);
   m=[m,m1(:)]; p=[p,p1(:)]; y=[y,y1];
end
figure(f1), 
subplot(211), semilogx(w,20*log10(m))
subplot(212), semilogx(w,p)
figure(f2), plot(t,y)

Gc

G=tf(100,[0.04,1,0]); gamma=60; wc=[50,100,200,300,500];
f1=figure; f2=figure; m=[]; p=[]; y=[];
w=logspace(-1,3); t=0:0.001:0.1;
for i=1:length(wc)
   Gc=leadlagc(G,wc(i),gamma,10);
   G_o=Gc*G; G_c=feedback(G_o,1);
   [m1,p1]=bode(G_o,w); y1=step(G_c,t);
   m=[m,m1(:)]; p=[p,p1(:)]; y=[y,y1];
end
figure(f1), 
subplot(211), semilogx(w,20*log10(m))
subplot(212), semilogx(w,p)
figure(f2), plot(t,y)