issteady.m

以《随机控制理论导论》为基础

function C=issteady(den);

% Distinguishes the steady condition of the system using Astrom's theory;
% den - denominator coeff. array of TF;
% Written by Paul Wong, Mar 25th, 2008

% Important!!!
% Dealing with numeric input only.


C=1;
A(1,:)=den;
n=length(den);
for k=2:n
    if A(k-1,1)<=0
        C=0;
        break;
    end
    for i=1:n-k+1
        A(k,i)=A(k-1,i)-(A(k-1,n-k+2)/A(k-1,1))*A(k-1,n+3-k-i);
    end
end