1
.TH obs_gram 1 "April 1993" "Scilab Group" "Scilab Function"
4
obs_gram - observability gramian
7
Go=obs_gram(A,C [,dom])
13
: real matrices (of appropriate dimensions)
16
: string (\fV"d'\fR or \fV"c"\fR (default value))
21
Observability gramian of the pair \fV(A,C)\fR or linear
22
system \fVsl\fR (\fVsyslin\fR list).
23
\fVdom\fR is the domain which can be
26
: continuous system (default)
35
Go = | e C'C e dt or Go = > A' C'C A
41
.LA Go = \int_0^{\infty} e^{A't}C'Ce^{At}dt
42
.LA \qquad Go = \sum_{0}^{\infty} {A'}^k C'C{A}^k
46
A=-diag(1:3);C=rand(2,3);
47
Go=obs_gram(A,C,'c'); // <=> w=syslin('c',A,[],C); Go=obs_gram(w);
48
norm(Go*A+A'*Go+C'*C,1)
49
norm(lyap(A,-C'*C,'c')-Go,1)
50
A=A/4; Go=obs_gram(A,C,'d'); //discrete time case
51
norm(lyap(A,-C'*C,'d')-Go,1)
54
ctr_gram, obsvss, obsv_mat, lyap