%線性系統 H(S)=(s+alph)/(s^3+10s^2+27s+18)，當alph分別取－1，0，1時，
%判別系統的可控性和可觀性，并求出相應的狀態方程。
clc
clear
more on
for alph=[-1:1]
   alph
   num=[1,alph];
   den=[1 10 27 18];
   [a,b,c,d]=tf2ss(num,den)
   cam=ctrb(a,b)
   rcam=rank(cam)
   oam=obsv(a,c)
   roam=rank(oam)
end
more off