clear%清空内存中的矩阵常量和变量
K=9;
num=1;
den=conv([1 0],[1 2]);
s1=tf(K*num,den);
%前向通路的传递函数
%从截止频率出发,进行校正
%*********************************
wc=4;%截止频率不小于4
L=bode(s1,wc);
Lwc=20*log10(L);%求出校正曲线wc等于20rad/s处的幅值
a=10^(0.1*Lwc);%计算a
T=1/(wc*sqrt(a));%计算T
%*********************************
%以上从截止频率出发进行校正
nc=[T,1];
dc=[a*T,1];
sysc=tf(nc,dc);%求校正后的传递函数