ex3_23.m

来自「反馈控制系统分析与设计,虽然有点抽象」· M 代码 · 共 9 行

G=tf([2,1],[conv([1,1],[1,2]),0]);
Gc=tf([5,6],[1,0]); H=tf(1,[0.01,1]);
[S,T]=sensifcn(G,Gc,H);  G1=Gc*S;
w=logspace(-2,2);
mag_a=abs(polyval(G1.den{1},w))./abs(polyval(G1.num{1},w));
semilogx(w,20*log10(mag_a));
mag_m=abs(polyval(T.den{1},w))./abs(polyval(T.num{1},w));
figure; semilogx(w,20*log10(mag_m));