MATLAB在系统综合校正中的应用

%% 
% MATLAB在系统综合校正中的应用
% eg.某单位负反馈系统校正前开环传递函数为G1(s)=100/(s*(0.04*s+1)*(0.01*s+1)),
% 校正后开环传递函数为G2(s)=100*(0.5*s+1)/(s*(5*s+1)*(0.04*s+1)*(0.01*s+1)),
% 利用MATLAB求校正前后相位裕度，校正前后系统是否稳定？
num=[100];  %校正前
den=conv([0.04 1 0],[0.01 1]);
sys=tf(num,den)
[gm,pm,wcg,wcp]=margin(sys)  %gm为幅值裕量，wcg为幅值裕量处的频率值；
% pm为相位裕量，wcp为其频率值
margin(sys)  %开环伯德图，由图知接近临界稳定，稳定储备差
grid

 

 num=[50 100]; %校正后
den1=conv([5 1 0],[0.04,1])
den=conv(den1,[0.01 1])
sys=tf(num,den)
[gm,pm,wcg,wcp]=margin(sys)
margin(sys)  %开环伯德图，由图知，校正后系统稳定，且稳定裕度有较大提高
grid

 

 

%% 
% MATLAB在系统综合校正中的应用
% eg.某单位负反馈系统校正前开环传递函数为G1(s)=100/(s*(0.04*s+1)*(0.01*s+1)),
% 校正后开环传递函数为G2(s)=100*(0.5*s+1)/(s*(5*s+1)*(0.04*s+1)*(0.01*s+1)),
% 利用MATLAB求校正前后相位裕度，校正前后系统是否稳定？
num=[100];  %校正前
den=conv([0.04 1 0],[0.01 1]);
sys=tf(num,den)
[gm,pm,wcg,wcp]=margin(sys)  %gm为幅值裕量，wcg为幅值裕量处的频率值；
% pm为相位裕量，wcp为其频率值
margin(sys)  %开环伯德图，由图知接近临界稳定，稳定储备差
grid

num=[50 100]; %校正后
den1=conv([5 1 0],[0.04,1])
den=conv(den1,[0.01 1])
sys=tf(num,den)
[gm,pm,wcg,wcp]=margin(sys)
margin(sys)  %开环伯德图，由图知，校正后系统稳定，且稳定裕度有较大提高
grid