串联滞后-超前校正及matlab代码实现

想想叫啥名

已于 2024-04-23 13:23:39 修改

阅读量978

点赞数 5

文章标签：滞后-超前校正自动控制原理

于 2024-04-23 13:07:26 首次发布

本文链接：https://blog.csdn.net/qq_35239859/article/details/138121303

版权

syms b_1;
K = 100; 
s = tf('s'); 
G0 = K/(s*(s/10+1)*(s/60+1)); % 原系统
gamma = 60; %期望相角裕度
% h = 10; %期望幅值裕度
wc = 10; % 期望截止频率，将该处设为校正点
lb = 0.1; 
ub = 1000; %3个10倍频程
[mag,phase,w]=bode(G0);
gamma0_wc=180+spline(w,phase,wc);
phim = gamma - gamma0_wc + 6;  %设置校正装置最大相角
a = (1+sind(phim))/(1-sind(phim)); % 求解超前校正参数a
T=1/(wc*sqrt(a));
Gclead=tf([a*T,1],[T,1]);
Glead = Gclead*G0; % 超前校正后系统传函
[mag_c,phase_c,w_c]=bode(Glead);
mag_wm=20*log10(spline(w_c,mag_c,wc));% 超前校正后系统wc处幅值
solve(-20*log10(b_1)==mag_wm,b_1);% 利用校正点处20lgb的对称性求解b
b=double(ans); %利用double将原本ans为syms类型转化为数值型
T=1/(0.1*wc*b);
Gclag=tf([b*T,1],[T,1]);
Gc = Gclag*Gclead; % 滞后超前校正函数
G = Gc*G0; % 校正后系统
bode(G0,Gc,G,{lb,ub});
grid on