ex6_3.m

来自「matlab 程序 专用于解答多项式 和 自动控制的各类习题 从我们老师那里」· M 代码 · 共 13 行

% Gc=(1+b.Tbs)/(1+Ts)*(1+a*Tas)/(1+Tas)
b=0.5;
Tb=8;
a=2;
Ta=0.5;
num=conv([b*Tb,1],[a*Ta,1]);
den=conv([Tb,1],[Ta,1]);
g=tf(num,den);

[x0,y0,w]=bode(g);
[x,y]=bd_asymp(g,w);
subplot(211), semilogx(w,20*log10(x0(:)),x,y)
subplot(212), semilogx(w,y0(:))