program3_5.m

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

%MATLAB Program3_5
A=[-2 2 -1;0 -2 0;1 -4 0];
B=[0;0;1];
C=[1 -1 1];
D=[0];
%判断系统的能控性
Qc=ctrb(A,B);
n=rank(Qc);
L=length(A);
if n ==L
    str='系统是状态完全能控的'
else 
    str='系统状态不完全能控，对系统按能控性进行结构性分解'
    [Abar Bbar Cbar T K]=ctrbf(A,B,C)
end

%判断系统的能观察性
Q0=obsv(A,C);
r=rank(Q0);
L=length(A);
if r ==L
    str='系统是状态完全能观测'
else 
    str='系统不完全能观测，对系统按能观测进行结构分解'
    [Abar Bbar Cbar T K]=obsvf(A,B,C)
end