xiandai4.m

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

A1=[2 0 0;0 4 1;0 0 4];B1=[1;0;1];C1=[1 1 0];co1=ctrb(A1,B1);m1=rank(co1);
N1=size(A1);n1=N1(1);%能控
if m1==n1
    disp('系统1是能控的')
else
    disp('系统1是不能控的')
end
a=poly(A1);   %  syms r;I=eye(n1);expand(det(r*I-A1));
T1=fliplr(co1)*[1 0 0;-10 1 0;32 -10 1];
A=inv(T1)*A1*T1
B=inv(T1)*B1
C=C1*T1