program3_4.m

此为matlab的代码。包括数值o45算法及o23算法

%MATLAB Program3_4
A=[4 1 0 0;0 4 1 0;0 0 4 1;0 0 0 4];
B=[0 0;1 2;0 0;2 1];
C=[1 0 2 0;2 0 4 2];

%判断系统的能控性
Qc=ctrb(A,B);
n=rank(Qc);
L=length(A);
if n ==L
    str='系统是状态完全能控的'
else 
    str='系统是状态不完全能控的'
end

%判断系统的能观察性
Q0=obsv(A,C);
r=rank(Q0);
L=length(A);
if r ==L
    str='系统是状态完全能观测'
else 
    str='系统是状态不完全能观测'
end