p200.m

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

alpha=1.2; T=1; alpha0=1.5:0.5:5;
G=tf([alpha*T,1]/alpha,[T,1]);
[x,y,w1]=nyquist(G,{.01,10});  [m,p,w2]=bode(G,{.01,10});
x=x(:); y=y(:); m=m(:); p=p(:);
for i=1:length(alpha0)
   G1=tf([alpha0(i)*T,1]/alpha0(i),[T,1]);
   [xx,yy]=nyquist(G1,w1); [mm,pp]=bode(G1,w2);
   x=[x,xx(:)]; y=[y,yy(:)]; 
   m=[m,mm(:)]; p=[p,pp(:)];
end
subplot(211); semilogx(w2,20*log10(m))
subplot(212); semilogx(w2,p); 
figure, plot([x,x],[y -y]), 
axis([0,1,-0.4,0.4]); axis('square')