ex412.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));          %α与β互为倒数
n=find(magdb+10*log10(alpha)<=0.0001);
%找出满足magdb+10*log10(alpha)<=0.0001式的magdb向量所有的下标值
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)  %绘出校正后的阶跃响应图