blob: b0c3f5898143c461c7bd41fa33d0e45b01d2d8ac (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
|
s=%s ;//convert to state space
TFcont=syslin ('c',2/(s^3+6*s^2+11*s+6))
SScont=tf2ss (TFcont)
[Ac ,Bc ,U, ind ]=canon( SScont( 2 ) , SScont( 3 ) )
disp(Ac,"A=")
disp(Bc,"B=")
C=[1 0 0]
p=cont_mat(Ac,Bc)
disp (p," controllability matrix=");
d=det(p)
if d==0
printf ("matrix is singular, so the system is uncontrollable");
else
printf ("system is controllable ");
end
g= obsv_mat (Ac,C);
disp (g," Observability Matrix=");
i= det(g)
if i ==0
printf ("matrix is singular, so the system is unobservable");
else
printf (" system is observable ");
end
|
