ex415.m

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

num=1600;den=conv([1,0],conv([1,2],[1,40]));
G=tf(num,den);   %未校正系统的传递函数
[h,gamma,wg,wc]=margin(G);h=20*log10(h);                %未校正系统的频域响应参数
w=0.001:0.001:100;        %设置计算频率点
[mag,phase]=bode(G,w);                          %未校正系统的幅频与相频特性值向量
disp('未校正系统的参数：h,wc,γ');
[ h,wc, gamma],
gamma1=40;delta=6;             %设计要求的相角裕度 ，(4.1.32)式中的 设置为 
phim=gamma1-gamma+delta;                             % (4.1.32)式
alpha=(1+sin(phim*pi/180))/(1-sin(phim*pi/180));             %(4.1.33)式
magdb=20*log10(mag);                         %用分贝值表示的未校正系统的幅值裕度
n=find(magdb+10*log10(alpha)<=0.0001);     %找出满足magdb+10*log10(alpha)<=0.0001式%的magdb向量所有的下标值
wc=n(1);wcc=wc/1000;                     % 的值与wc相差1000倍
w3=wcc/sqrt(alpha);w4=sqrt(alpha)*wcc;                  %(4.1.34)式
numc1=[1/w3,1];denc1=[1/w4,1];                        %(4.1.35)式
Gc1=tf(numc1,denc1);                                 %超前校正部分的传递函数
w1=wcc/10;w2=w1/alpha;                              %取 =  ； = 
numc2=[1/w1,1];denc2=[1/w2,1];                        %(4.1.36)式
Gc2=tf(numc2,denc2);                                 %滞后校正部分的传递函数
Gc12=Gc1*Gc2;                                  %串联超前-滞后校正网络的传递函数
GcG=Gc12*G;                                               %校正后系统的传递函数
[Gmc,Pmc,wcgc,wcpc]=margin(GcG);
GmcdB=20*log10(Gmc);
disp ('超前校正部分的传递函数'),Gc1,
disp ('滞后校正部分的传递函数'),Gc2,
disp ('串联超前-滞后校正网络的传递函数'),Gc12,
disp('校正后系统的开环传递函数'),GcG,
disp('校正后系统的性能参数：h,wc,γ及α值'),[GmcdB,wcpc,Pmc,alpha],
bode(G,GcG)