可観測性グラミアン
Go=obs_gram(A,C [,dom]) Go=obs_gram(sl)
(適当な次元の)実数行列
文字列 ("d' または "c" (デフォルト値))
"d'
"c"
syslin リスト
syslin
対(A,C)または 線形システムsl (syslin リスト)の 可観測性グラミアン. dom は領域で,以下の値とすることができます
(A,C)
sl
dom
連続システム (デフォルト)
離散システム
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'); //離散時間の場合 norm(lyap(A,-C'*C,'d')-Go,1)