ex422a0.m

来自「详细介绍了Mathlab的使用方法与技巧」· M 代码 · 共 11 行

num=10;den=conv([1,1],conv([1,3],[1,5]));G=tf(num,den);
for Km=0:0.1:10000
Gc=Km;GcG=feedback(Gc*G,1);
[num,den]=tfdata(GcG,'v');
p=roots(den);pr=real(p);prm=max(pr);
pr0=find(prm>=-0.001);n=length(pr0);
if n>=1
    break
end;end
step(GcG,0:0.001:3);Km