matlabT型失真矫正,用matlab的guide写了个滞后校正的

不知道有没有用

s &＃61; tf(&＃39;s&＃39;);

gs &＃61; input(&＃39;请输入原开环传递函数\n 例如: 1/(s*(0.1*s&＃43;1)*(0.05*s&＃43;1)) \ngs &＃61;&＃39;);

% 稳态误差

N &＃61; input(&＃39;请输入给定信号表达式\n 例如: 1/s^2 \n N &＃61; &＃39;);

E &＃61;s*N/(1&＃43;gs);

G1 &＃61; E;

clear s

syms s

[num,den] &＃61; tfdata(G1);%开始转换类型 tf->sym

Na &＃61; size(den);

for i&＃61;1:Na(1)

for k&＃61;1:Na(2)

Num&＃61;poly2sym(num{i,k},s);

Den&＃61;poly2sym(den{i,k},s);

sym_G1(i,k)&＃61;Num/Den;

end

end

assignin(&＃39;base&＃39;,&＃39;sym_G1&＃39;,sym_G1);

E &＃61; sym_G1;

ess2 &＃61; limit(E,s,0,&＃39;right&＃39;);

ess &＃61; input(&＃39;请输入期望的稳态误差 ess &＃61; &＃39;);

K &＃61; double(ess2)/ess ;%求得满足稳态误差的最小开环增益

clear s

s &＃61; tf(&＃39;s&＃39;);

disp([&＃39;最小开环增益K &＃61; &＃39;,num2str(K)])

K &＃61; input(&＃39;请输入开环增益K &＃61; &＃39;);

gs&＃61;K*gs%更新系统传递函数

figure(1)

margin(gs);

[gm0 pm0 wg0 wp0] &＃61; margin(gs);

gm0&＃61;20*log10(gm0);

disp([&＃39;原系统幅值裕度Lg &＃61; &＃39;,num2str(gm0),&＃39;dB &＃39;,&＃39;相角裕度γ &＃61; &＃39;,num2str(pm0)])

lg &＃61; input(&＃39;请输入目标幅值裕度Lg &＃61; &＃39;);

gama &＃61; input(&＃39;请输入目标相角裕度γ &＃61; &＃39;);

disp([&＃39;选用滞后校正 &＃39;])

% 截止频率

gama1&＃61;gama&＃43;5;

[mu pu w]&＃61;bode(gs);

wc&＃61;spline(pu,w,(gama1-180));%相对于pu的w的值的规律插值得到相对于(gama1-180)的值

disp([&＃39;截止频率wc &＃61; &＃39;,num2str(wc)]);

% 计算β

na&＃61;polyval(gs.num{1},j*wc);

da&＃61;polyval(gs.den{1},j*wc);

g&＃61;na/da;

g1&＃61;abs(g);

h&＃61;20*log10(g1);

beta&＃61;10^(h/20);

disp([&＃39;β &＃61; &＃39;,num2str(beta)]);

% 计算校正环节参数

T&＃61;10/wc

disp(&＃39;校正环节&＃xff1a;&＃39;)

gc&＃61;(T*s&＃43;1)/(beta*T*s&＃43;1)

disp(&＃39;校正后系统传递函数&＃xff1a;&＃39;)

g&＃61;gs*gc

disp(&＃39;校正后系统传伯德图&＃xff1a;&＃39;)

figure(2)

bode(g);

[gm pm wg wp]&＃61;margin(g);

disp([&＃39;校正后系统幅值裕度Lg &＃61; &＃39;,num2str(gm),&＃39;dB &＃39;,&＃39;相角裕度γ &＃61; &＃39;,num2str(pm)])

%阶跃响应

sys1&＃61;feedback(gs,1);

sys2&＃61;feedback(g,1);

figure(3)

hold on

step(sys1)

step(sys2)

hold off

figure(4)

hold on

bode(gs)

bode(g)

hold off

程序写的比较烂

滞后校正参数计算.mlappinstall

111.09KB

MLAPPINSTALL

28次下载