第1个回答 2013-06-08
k=1/0.01;
G0=zpk([],[0 -2],k);
[h0,r,wx,wc]=margin(G0)
wm=4.4
L=bode(G0,wm)
Lwc=20*log10(L)
a=10^(-0.01*Lwc)
T=1/(wm*sqrt(a));
phi=asin((a-1)/(a+1))
Gc=(1/a)*tf([a*T 1],[T 1]);
Gc=a*Gc;
G=Gc*G0;
bode(G,'r',G0,'b--');
grid;
[h,r,wx,wc]=margin(G)