example7_2a.m

<MATLAB 7辅助控制系统设计与仿真>书籍内的源代码.

%求解校正器的传递函数
kc=5;
s_1=-3.7518+9.7854i;
nk1=2;
dk1=conv(conv([1 0],[0.25 1]),[0.1 1]);
ngv=polyval(nk1,s_1);dgv=polyval(dk1,s_1);
g=ngv/dgv;zetag=angle(g);zetag_d=zetag*180/pi;
mg=abs(g);ms=abs(s_1);
zetas=angle(s_1);zetas_d=zetas*180/pi;
tz=(sin(zetas)-kc*mg*sin(zetag-zetas))/(kc*mg*ms*sin(zetag));
tp=-(kc*mg*sin(zetas)+sin(zetag+zetas))/(ms*sin(zetag));
nk=[tz,1];dk=[tp,1];
Gc=tf(nk,dk)