nk.m

该程序可以用于所以现代控制理论里的关于系统的能观性、能控性等判断

function  nk(A,B)
co=ctrb(A,B);N=size(A);n=N(1);rs=rank(co);   %  rs为系统的控制矩阵的秩，co为控制矩阵,n为系统维数
if rs==n                                 %   N1=size(A1);n1=N1(1);
    disp('系统是可控的')
else
    disp('系统是不可控的')
end