xiandai3.m

来自「该程序可以用于所以现代控制理论里的关于系统的能观性、能控性等判断」· M 代码 · 共 15 行

A1=[0 1 0;0 0 1;-6 -11 -6];B1=[0;0;1];C1=[4 5 1];%N1=[C1;C1*A1;C1*A1*A1];  不能观
ob1=obsv(A1,C1);N1=size(A1);n1=N1(1);rs1=rank(ob1);
if rs1==n1
    disp('系统1是可观的')
else
    disp('系统1是不可观的')
end
A2=[0 1 0 0;0 0 0 2;0 0 0 1;0 -2 5 0];B2=[0;1;0;0];C2=[1 0 0 0];
ob2=obsv(A2,C2);N2=size(A2);n2=N2(1);rs2=rank(ob2);
if rs2==n2
    disp('系统2是可观的')
else
    disp('系统2是不可观的')
end
%N2=[C2;C2*A2;C2*A2*A2;C2*A2*A2*A2];n2=rank(N2)  %能观