超前滞后电气与电子测量技术.pdf

超前滞后 matlab 程序1： %超前校正 pos=input(Type %OS ); % 输入期望超调量 Tp=input(Type peak time ); % 期望的峰值时间 Kv=input(Type value of Kv ); % 期望的Kv numg=[10]; % 定义传递函数的分子 deng=poly([0 -36 -100]); % 定义传递函数的分母deng=poly([0 -36 -100]) G=tf(numg,deng); s=tf([1 0],1); % 创建传递函数,s sG=s*G; sG=minreal(sG); % 零极点对消 K=dcgain(Kv/sG); % 计算增益K. G(s) G=zpk(K*G) % 将传递函数转换为零极点模型 z=(-log(pos/100))/(sqrt(pi^2+log(pos/100)^2)); % 计算所需的阻尼系数 Pm=atan(2*z/(sqrt(-2*z^2+sqrt(1+4*z^4))))*(180/pi); % 计算相角裕度 wn=pi/(Tp*sqrt(1-z^2)); % 计算所需固有频率 wBW=wn*sqrt((1-2*z^2)+sqrt(4*z^4-4*z^2+2)); % 选取带宽 w=0.01:0.5:1000; % 设置频率范围 [M,P]=bode(G,w); % 获取Bode 数据 [Gm,Pm,Wcg,Wcp]=margin(G); % 找到当前相角裕度 Pmreq=atan(2*z/(sqrt(-2*z^2+sqrt(1+4*z^4))))*(180/pi); % 计算所需相角裕度 Pmreqc=Pmreq+10; % 加一个10 度的修正量 Pc=Pmreqc-Pm; % 计算超前校正所需的相角 % 设计超前校正 beta=(1-sin(Pc*pi/180))/(1+sin(Pc*pi/180)); % 找到β magpc=1/sqrt(beta); % 找到校正装置的最大幅值 for k=1:1:length(M); % 找到未校正系统幅值为1/magpc 的频率 % 这个频率是新的相位裕量频率 if M(k)-(1/magpc)=0; % 找最大幅值 wmax=w(k); % 最大幅值处的频率 break % 停止循环 end end % 计算超前校正的零极点以及增益 zc=wmax*sqrt(beta); % 找到超前校正器的最低截止频率 pc=zc/beta; % 找到超前校正器的最高截止频率 Kc=1/beta; % 计算超前校正器的增益 Gc(s) Gc=tf(Kc*[1 zc],[1 pc]); Gc=zpk(Gc) Ge(s)=G(s)Gc(s) Ge=G*Gc sGe=s*Ge; sGe=minreal(sGe); % 零极点对消 Kv=dcgain(sGe)