ex5_1.m

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

G=tf(1,[conv([1,1],[1,2]) 0]);
hold off; rlocus(G);  axis('square')
[xx,yy]=rlasymp(G); hold on; plot(xx,yy,':')
[K,P]=rlocfind(G)

hold off; t=0:0.2:15; Y=[];
for k=[0.1:0.1:1, 2:5]
   Gk=feedback(k*G,1); 
   y=step(Gk,t); Y=[Y y];
end
plot(t,Y)