intstable.m

控制系统计算机辅助设计——MATLAB语言与应用(源代码）

function key=intstable(G,Gc,H)
GG=minreal(feedback(G*Gc,H)); Go=H*G*Gc; Go1=minreal(Go); p=eig(GG);
z0=eig(Go); z1=eig(Go1); zz=setdiff(z0,z1); % 找出开环对销极点
if (G.Ts>1), % 离散系统判定
   key=any(abs(p)>1); if key==0, key=2*any(abs(zz)>1); end
else, % 连续系统判定
   key=any(real(p)>0); if key==0, key=2*any(real(zz)>0); end
end