%確定系統的可控性和可觀測性
A=[-2 3;3 -2];
B=[1 1.2;1 3];
C=[2 1;1 -2];
D=0;
n=length(A);
OB=obsv(A,C)
rOB=rank(OB)%可觀測性判斷
if rOB==n
  disp('System is observable')
elseif rOB<n
  disp('System is unobservable')
end
CO=ctrb(A,B)
rCO=rank(CO)%可控性判斷
if rCO==n
  disp('System is controllable')
elseif rCO<n
  disp('System is uncontrollable')
end