xiandai2.m

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

A1=[1 1 0;0 1 0;0 1 1];B1=[0;1;0];%  M1=[B1 A1*B1 A1*A1*B1];
co1=ctrb(A1,B1);N1=size(A1);n1=N1(1);rs1=rank(co1); %不能控  n1=rank(M1)
if rs1==n1
    disp('系统1是可控的')
else
    disp('系统1是不可控的')
end
A2=[1 3 2;0 2 0;0 1 2];B2=[2 1;1 1;-1 -1];  % M2=[B2 A2*B2 A2*A2*B2];n2=rank(M2);
co2=ctrb(A2,B2);N2=size(A2);n2=N2(1);rs2=rank(co2);  %能控
if rs2==n2
    disp('系统2是可控的')
else
    disp('系统2是不可控的')
end

A3=[0 1 0 0;0 0 -1 0;0 0 0 1;0 0 5 0];B3=[0;1;0;-2];
% M3=[B3 A3*B3 A3*A3*B3 A3*A3*A3*B3];n3=rank(M3);
co3=ctrb(A3,B3);N3=size(A3);n3=N3(1);rs3=rank(co3); %能控
if rs3==n3
    disp('系统3是可控的')
else
    disp('系统3是不可控的')
end
A4=[0 1 0 0;3 0 0 2;0 0 0 1;0 -2 0 0];B4=[0;1;0;0];
% M4=[B4 A4*B4 A4*A4*B4 A4*A4*A4*B4];n4=rank(M4); %不能控
co4=ctrb(A4,B4);N4=size(A4);n4=N4(1);rs4=rank(co4);
if rs4==n4
    disp('系统4是可控的')
else
    disp('系统4是不可控的')
end