di44.m

详细介绍了Mathlab的使用方法与技巧

num=7;den=conv([1,0],[1/6,1]);
den=conv(den,[1/2,1]);
G=tf(num,den);
w=0.01:0.01:1000;
[mag,phase]=bode(G,w);
%[Gm,Pm,Wcg,Wcp]=margin(G);    %未校正系统的频域响应参数，计算需要相角欲度Pm(γ)
gamma_cas=40;delta=6;gamma_l=gamma_cas+delta;
n=find(180+phase-(gamma_l)<=0.1);
wgamma_l=n(1)/100;                %n(1)的值与实际ω'c频率相差100倍     
[mag,phase]=bode(G,wgamma_l);
rr=-20*log10(mag);beta=10^(rr/20);
w2=wgamma_l/10;w1=beta*w2;
numc=[1/w2,1];denc=[1/w1,1];Gc=tf(numc,denc),GcG=Gc*G
margin(GcG),beta