stabdetc.m

来自「用MATLAB的脚本语言M语言编写的」· M 代码 · 共 16 行

function [V,ss]=stabdetc(A,X)
n=length(A); [nc,nr]=size(X); ss=[];
if nc==n & nr==1, A0=A; B0=X;
elseif nc==1 & nr==n, A0=A'; B0=X';
else, 
   V=[]; error('uncompetible (A,B)');
end
C=ctrb(A0,B0); nB=rank(C);
if nB==n, V=1; 
else
   [Ac,Bc,Cc]=ctrbf(A0,B0,ones(1,n)); 
   Anc=Ac(1:n-nB,1:n-nB); ee=eig(Anc);
   ii=find(real(ee)>=0); ss=ee(ii);
   V=~(length(ss)>0);
end