example7_4a.m

来自「这是一个MATLAB6。5 辅助控制系统设计与仿真例题的所有源码」· M 代码 · 共 13 行

% 计算超前校正器的传递函数
k0=1000;n1=1;d1=conv([1 0],[1 2]);
sope=tf(k0*n1,d1);
[mag,phase,w]=bode(sope);
gama=45;[mu,pu]=bode(sope,w);
gam=gama*pi/180;
alfa=(1-sin(gam))/(1+sin(gam));
adb=20*log10(mu);am=10*log10(alfa);
ca=adb+am;wc=spline(adb,w,am);
T=1/(wc*sqrt(alfa));
alfat=alfa*T;
Gc=tf([T 1],[alfat 1])