ex411.m

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

num=1000;den=conv([1,0],conv([0.1,1],[0.001,1]));
G=tf(num,den);                %未校正系统的开环传递函数
[Gm,Pm,Wcg,Wcp]=margin(G);  %1) 未校正系统的频域响应参数，计算需要相角欲度Pm(γ)
w=0.1:0.1:10000;              %确定频率的取值范围和频率采样的间隔值
[mag,phase]=bode(G,w);magdb=20*log10(mag);  %对数幅频、相频特性向量，计算需要magdb
%(对数幅频响应值)
phim1=45;deta=8;              % 2）(4.1.9)式设置系统设计参数phim1( )；deta( )
phim=phim1-Pm+deta;                        % (4.1.9)式相位超前角φ 
bita=(1-sin(phim*pi/180))/(1+sin(phim*pi/180));        %求出β值
n=find(magdb+10*log10(1/bita)<=0.0001);       %找出满足magdb+10*log10(1/bita)<=0.0001
%式的magdb向量所有的下标值
wc=n(1);      %通常magdb(1)+10*log10(1/bita)>0.0001;取 的第1项为 (比实际频率值大
%10倍)，这是因为w=0.1:0.1:10000，而下标向量 的值=1:1:100000
w1=(wc/10)*sqrt(bita);w2=(wc/10)/sqrt(bita);               % 3)  (4.1.11)和(4.1.12)式
numc=[1/w1,1];denc=[1/w2,1];                           % (4.1.13)式，取K=1 
Gc=tf(numc,denc);                                     %校正器的传递函数                       
GmdB=20*log10(Gm);                         %4) 下面程序校验系统校正后的系统指标
GcG=Gc*G;[Gmc,Pmc,wcgc,wcpc]=margin(GcG);        %GcG是校正后系统的开环传递函数
GmcdB=20*log10(Gmc);                       
disp('未校正系统的开环传递函数和频域响应参数：h，γ，wc')
G,[GmdB,Pm,Wcp],
disp('校正装置传递函数和校正后系统开环传递函数')
Gc,GcG,
disp('校正后系统的频域响应参数：h，γ，wc')
[GmcdB,Pmc,wcpc],
disp('校正装置的参数T和β值：T，β')
T=1/w1;[T,bita],
bode(G,GcG);figure(2);margin(GcG)