Matlab 实现系统校正

function  revise = bode(w)%w:校正后的截止频率
num1 = 40;
den1 = conv([0.2 1 0],[0.0625 1]);
sys1 = tf(num1,den1);
% [mag,phase,17.2]=bode(sys); 
figure('name','sys1未校正系统')%创建一个窗口
margin(sys1);
%
[mag]=bode(sys1,w);
MAG = 20*log10(mag);
%
a = 10^(-MAG/10);%超前
Tc = 1/(w*sqrt(a));
num2 = conv([40],[a*Tc 1]);
den2 = conv([0.01251 0.26255 1 0],[Tc 1]);
sys2 = tf(num2,den2);
figure('name','超前校正系统')%创建一个窗口
margin(sys2);
%
b = 10^(-MAG/20);%滞后
Tz = 1/(0.1*w*b);
numh = conv([40],[b*Tz 1]);
denh = conv([0.01251 0.26255 1 0],[Tz 1]);
sysh = tf(numh,denh);
figure('name','滞后校正系统')%创建一个窗口
margin(sysh);
%
revise = [a Tc b Tz];

未校正系统伯德图