di42.m

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

num=10;den=[0.2,1,0,0];G=tf(num,den);         %未校正系统的传递函数
[Gm,Pm,Wcg,Wcp]=margin(G);    %未校正系统的频域响应参数，计算需要相角欲度Pm(γ)
w=0.1:1:10000;                %频率采样点的确定，注意起始频率为0.1
[mag,phase]=bode(G,w);magdb=20*log10(mag);  %对数幅频、相频特性值，计算需要magdb(对数幅频响应值)
phim1=35;deta=18;              %(4.1.9)式设置系统设计参数phim1(γ')；deta(△)
phim=(phim1-Pm+deta)/2;                        %(4.1.9)式相位超前角φ,取φ/2 
alpha=(1+sin(phim*pi/180))/(1-sin(phim*pi/180));  %α与β互为倒数
%w=1:1:100000;               %w在这里作为magdb向量的下标，比实际频率值小0.1
n=find(magdb+10*log10(alpha)<=0.0001); %找出满足magdb+10*log10(alpha)<=0.0001式的magdb向量所有的下标值(w值)
wc=n(1)+0.1;                         %通常magdb(1)+10*log10(1/bita)>0.0001;故取n的第1项为wc+0.1
w1=wc/sqrt(alpha);w2=sqrt(alpha)*wc;        %(4.1.11)和(4.1.12)式，注意：α与β互为倒数
numc=(1/alpha)*[1/w1,1];denc=[1/w2,1];      %(4.1.13) 式，其中K=1/α
GG=tf(numc,denc);Gc=GG*GG;                     %Gc是2个校正网络的串联
GcG=(alpha)^2*Gc*G;                             % α^2是校正网络增加的放大器
disp('显示单级校正网络传递函数、2级串联校正网络传递函数及α、T值,')
T=1/w2;GG,Gc,[alpha,T],margin(GcG)
figure(2);sys=feedback(GcG,1);step(sys)