ex6_20.m

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

num=[3,5,3,2]; wn=30;
den=conv([1,3,2],conv([1,1,3],[1,0.4,1]));
G=tf(num,den); Gs=std_tf(1,wn,4);
Gc=inv(G)*feedback(Gs,1,1); zpk(Gc)

KK=[0.8:0.1:1.5];
for i=1:length(KK)
   Gc1=KK(i)*Gc; G_c=feedback(G*Gc1,1);
   step(G_c); hold on
end