exp3_4.m

来自「一个很好的例子,是一本报书籍上的,大家好好看看啊 ,一定不回失望的啊」· M 代码 · 共 17 行

%线性系统 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