obs_gram

obs_gram — observability gramian

Calling sequence

Go=obs_gram(A,C [,dom])  
Go=obs_gram(sl)  

Parameters

A,C : real matrices (of appropriate dimensions)
dom : string ("d' or "c" (default value))
sl : syslin list

Description

Observability gramian of the pair (A,C) or linear system sl (syslin list). dom is the domain which can be

"c": continuous system (default)
"d": discrete system

Examples



A=-diag(1:3);C=rand(2,3);
Go=obs_gram(A,C,'c');     // <=> w=syslin('c',A,[],C); Go=obs_gram(w);
norm(Go*A+A'*Go+C'*C,1)
norm(lyap(A,-C'*C,'c')-Go,1)
A=A/4; Go=obs_gram(A,C,'d');    //discrete time case
norm(lyap(A,-C'*C,'d')-Go,1)
 
  

See also

ctr_gram, obsvss, obsv_mat, lyap