In control theory, visible state of a system
観測可能性とは、 システム の内部状態を外部出力の知識からどれだけ正確に推測できるか を示す尺度である。 制御理論 において、線形システムの 観測可能性と 制御可能性は数学的に 双対関係 にある。
観測可能性の概念は、ハンガリー系アメリカ人エンジニアの ルドルフ・E・カルマン によって線形動的システムに対して導入されました。 [1] [2] 出力の測定からシステムの状態を推定するように設計された動的システムは、 カルマンフィルタ などのそのシステムの 状態観測器 と呼ばれます。
意味 状態空間表現 でモデル化された物理システムを考えてみましょう 。 状態ベクトルと制御ベクトル のあらゆる可能な変化に対して、出力(物理的には、これは一般的に センサーによって得られる情報に相当 し ます)からの情報のみを用いて現在の状態を推定できる場合 、システムは観測可能であると言われます。言い換えれば、システムの出力からシステム全体の動作を判断できるということです。一方、システムが観測可能でない場合、出力の測定だけでは区別できない状態軌跡が存在します。
線形時間不変システム 状態空間表現における時間不変線形システム の場合 、システムが観測可能かどうかを確認するための便利なテストがあります。 状態変数( MIMO システムの詳細については 状態空間 を参照)を持つ SISO システムを考えます。 n {\displaystyle n}
x ˙ ( t ) = A x ( t ) + B u ( t ) {\displaystyle {\dot {\mathbf {x} }}(t)=\mathbf {A} \mathbf {x} (t)+\mathbf {B} \mathbf {u} (t)} y ( t ) = C x ( t ) + D u ( t ) {\displaystyle \mathbf {y} (t)=\mathbf {C} \mathbf {x} (t)+\mathbf {D} \mathbf {u} (t)}
可観測性マトリックス 観測可能性行列 の 列 ランクが 、次のように定義される
とき、かつその場合に限り、
O = [ C C A C A 2 ⋮ C A n − 1 ] {\displaystyle {\mathcal {O}}={\begin{bmatrix}C\\CA\\CA^{2}\\\vdots \\CA^{n-1}\end{bmatrix}}} が に等しい 場合、システムは観測可能である。この検定の根拠は、 列が線形独立であれば、各 状態変数は出力変数の線形結合を通して観測可能であるということである。観測可能性は、連続時間 状態観測器 の設計における必要十分条件である 。 n {\displaystyle n} n {\displaystyle n} n {\displaystyle n} y {\displaystyle y}
可観測性指数 線形時間不変離散システムの観測 可能性指数 は、次式を満たす最小の自然数である。 ここで、 v {\displaystyle v} rank ( O v ) = rank ( O v + 1 ) {\displaystyle {\text{rank}}{({\mathcal {O}}_{v})}={\text{rank}}{({\mathcal {O}}_{v+1})}}
O v = [ C C A C A 2 ⋮ C A v − 1 ] . {\displaystyle {\mathcal {O}}_{v}={\begin{bmatrix}C\\CA\\CA^{2}\\\vdots \\CA^{v-1}\end{bmatrix}}.}
観測不可能な部分空間 線形システムの 観測不可能な部分空間は 、 [3] で与えられる線形写像の核である。 N {\displaystyle N} G {\displaystyle G}
G : R n → C ( R ; R n ) x ( 0 ) ↦ C e A t x ( 0 ) {\displaystyle {\begin{aligned}G\colon \mathbb {R} ^{n}&\rightarrow {\mathcal {C}}(\mathbb {R} ;\mathbb {R} ^{n})\\x(0)&\mapsto Ce^{At}x(0)\end{aligned}}}
ここでは から まで の連続関数の集合である 。 は次のようにも書ける [3] C ( R ; R n ) {\displaystyle {\mathcal {C}}(\mathbb {R} ;\mathbb {R} ^{n})} R {\displaystyle \mathbb {R} } R n {\displaystyle \mathbb {R} ^{n}} N {\displaystyle N}
N = ⋂ k = 0 n − 1 ker ( C A k ) = ker O {\displaystyle N=\bigcap _{k=0}^{n-1}\ker(CA^{k})=\ker {\mathcal {O}}} システムが観測可能であるのは の場合のみなので 、システムが観測可能であるのは が ゼロ部分空間である場合のみである。 rank ( O ) = n {\displaystyle \operatorname {rank} ({\mathcal {O}})=n} N {\displaystyle N}
観測不可能な部分空間には以下の性質がある: [3]
N ⊂ K e ( C ) {\displaystyle N\subset Ke(C)} A ( N ) ⊂ N {\displaystyle A(N)\subset N} N = ⋃ { S ⊂ R n ∣ S ⊂ K e ( C ) , A ( S ) ⊂ N } {\displaystyle N=\bigcup \{S\subset R^{n}\mid S\subset Ke(C),A(S)\subset N\}}
検出可能性 観測可能性よりもやや弱い概念が 検出可能性 である。観測不可能な状態がすべて安定している場合、システムは検出可能である。 [4]
検出可能性条件はセンサーネットワーク の文脈において重要である 。 [5] [6]
機能的観測可能性 機能的観測 可能性は、完全な状態観測が不可能な場合(測定信号の欠如のため)に古典的な観測可能性の概念を拡張した性質であり、代わりに出力からの情報のみを使用して 線形関数 を推定できる条件を確立する。 [7] 正式には、(通常は低次元の) 行列 ( ただし)が与えられた 場合、システムが機能的に観測可能であるのは、 [8] z ( t ) = F x ( t ) {\displaystyle \mathbf {z} (t)=\mathbf {F} \mathbf {x} (t)} r × n {\displaystyle r\times n} F {\displaystyle \mathbf {F} } r ≤ n {\displaystyle r\leq n}
rank [ O F ] = rank O . {\displaystyle \operatorname {rank} {\begin{bmatrix}{\mathcal {O}}\\\mathbf {F} \end{bmatrix}}=\operatorname {rank} {\mathcal {O}}.} 機能的観測可能性は、機能的観測器(ダルーアッシュ観測器 [9] とも呼ばれる )が漸近的に推定するように設計できる 十分条件と必要条件を決定する重要な概念である 。特定の条件下では、機能的観測可能性と出力制御可能性は数学的に 双対関係 にある。 [10] z ( t ) {\displaystyle \mathbf {z} (t)}
線形時間変動システム 連続 線形 時変システム を考える
x ˙ ( t ) = A ( t ) x ( t ) + B ( t ) u ( t ) {\displaystyle {\dot {\mathbf {x} }}(t)=A(t)\mathbf {x} (t)+B(t)\mathbf {u} (t)\,} y ( t ) = C ( t ) x ( t ) . {\displaystyle \mathbf {y} (t)=C(t)\mathbf {x} (t).\,} 行列 、 、 が与えられ、入力と出力と が全てについて与えられているとすると、 によって 定義 さ れる の 零 空間 に ある加法定数ベクトル内で を 決定することが可能である。 A {\displaystyle A} B {\displaystyle B} C {\displaystyle C} u {\displaystyle u} y {\displaystyle y} t ∈ [ t 0 , t 1 ] ; {\displaystyle t\in [t_{0},t_{1}];} x ( t 0 ) {\displaystyle x(t_{0})} M ( t 0 , t 1 ) {\displaystyle M(t_{0},t_{1})}
M ( t 0 , t 1 ) = ∫ t 0 t 1 φ ( t , t 0 ) T C ( t ) T C ( t ) φ ( t , t 0 ) d t {\displaystyle M(t_{0},t_{1})=\int _{t_{0}}^{t_{1}}\varphi (t,t_{0})^{T}C(t)^{T}C(t)\varphi (t,t_{0})\,dt} ここで、 は 状態遷移行列 です。 φ {\displaystyle \varphi }
が特異 でない 場合には 、 を一意に決定することができます 。実際、 が の零空間にある 場合には 、 の初期状態と の初期状態を区別することはできません 。 x ( t 0 ) {\displaystyle x(t_{0})} M ( t 0 , t 1 ) {\displaystyle M(t_{0},t_{1})} x 1 {\displaystyle x_{1}} x 2 {\displaystyle x_{2}} x 1 − x 2 {\displaystyle x_{1}-x_{2}} M ( t 0 , t 1 ) {\displaystyle M(t_{0},t_{1})}
上記のように定義された行列には 次の特性があることに注意してください。 M {\displaystyle M}
M ( t 0 , t 1 ) {\displaystyle M(t_{0},t_{1})} 対称的 である M ( t 0 , t 1 ) {\displaystyle M(t_{0},t_{1})} は 、 t 1 ≥ t 0 {\displaystyle t_{1}\geq t_{0}} M ( t 0 , t 1 ) {\displaystyle M(t_{0},t_{1})} 線形 行列微分方程式を満たす d d t M ( t , t 1 ) = − A ( t ) T M ( t , t 1 ) − M ( t , t 1 ) A ( t ) − C ( t ) T C ( t ) , M ( t 1 , t 1 ) = 0 {\displaystyle {\frac {d}{dt}}M(t,t_{1})=-A(t)^{T}M(t,t_{1})-M(t,t_{1})A(t)-C(t)^{T}C(t),\;M(t_{1},t_{1})=0} M ( t 0 , t 1 ) {\displaystyle M(t_{0},t_{1})} 方程式を満たす M ( t 0 , t 1 ) = M ( t 0 , t ) + φ ( t , t 0 ) T M ( t , t 1 ) φ ( t , t 0 ) {\displaystyle M(t_{0},t_{1})=M(t_{0},t)+\varphi (t,t_{0})^{T}M(t,t_{1})\varphi (t,t_{0})} [11]
観測可能性行列の一般化 システムが で観測可能であるのは、 行列が特異でないような に 区間 が存在する場合のみです 。 [ t 0 , t 1 ] {\displaystyle [t_{0},t_{1}]} [ t 0 , t 1 ] {\displaystyle [t_{0},t_{1}]} R {\displaystyle \mathbb {R} } M ( t 0 , t 1 ) {\displaystyle M(t_{0},t_{1})}
が解析的である場合 、系は区間[ , ]で観測可能であり 、かつ正の整数 k が存在し 、 [12] A ( t ) , C ( t ) {\displaystyle A(t),C(t)} t 0 {\displaystyle t_{0}} t 1 {\displaystyle t_{1}} t ¯ ∈ [ t 0 , t 1 ] {\displaystyle {\bar {t}}\in [t_{0},t_{1}]}
rank [ N 0 ( t ¯ ) N 1 ( t ¯ ) ⋮ N k ( t ¯ ) ] = n , {\displaystyle \operatorname {rank} {\begin{bmatrix}&N_{0}({\bar {t}})&\\&N_{1}({\bar {t}})&\\&\vdots &\\&N_{k}({\bar {t}})&\end{bmatrix}}=n,} ここで 、および は再帰的に次のように定義される。 N 0 ( t ) := C ( t ) {\displaystyle N_{0}(t):=C(t)} N i ( t ) {\displaystyle N_{i}(t)}
N i + 1 ( t ) := N i ( t ) A ( t ) + d d t N i ( t ) , i = 0 , … , k − 1 {\displaystyle N_{i+1}(t):=N_{i}(t)A(t)+{\frac {\mathrm {d} }{\mathrm {d} t}}N_{i}(t),\ i=0,\ldots ,k-1}
例 行列 が解析的に変化するシステムを考える ( − ∞ , ∞ ) {\displaystyle (-\infty ,\infty )}
A ( t ) = [ t 1 0 0 t 3 0 0 0 t 2 ] , C ( t ) = [ 1 0 1 ] . {\displaystyle A(t)={\begin{bmatrix}t&1&0\\0&t^{3}&0\\0&0&t^{2}\end{bmatrix}},\,C(t)={\begin{bmatrix}1&0&1\end{bmatrix}}.}
すると 、この行列の階数は 3 なので、このシステムは のすべての非自明な区間で観測可能です 。 [ N 0 ( 0 ) N 1 ( 0 ) N 2 ( 0 ) ] = [ 1 0 1 0 1 0 1 0 0 ] {\displaystyle {\begin{bmatrix}N_{0}(0)\\N_{1}(0)\\N_{2}(0)\end{bmatrix}}={\begin{bmatrix}1&0&1\\0&1&0\\1&0&0\end{bmatrix}}} R {\displaystyle \mathbb {R} }
非線形システム システム が与えられます 。 ここで、 状態ベクトル、 入力ベクトル、および 出力ベクトル は、滑らかなベクトル場になります。 x ˙ = f ( x ) + ∑ j = 1 m g j ( x ) u j {\displaystyle {\dot {x}}=f(x)+\sum _{j=1}^{m}g_{j}(x)u_{j}} y i = h i ( x ) , i ∈ p {\displaystyle y_{i}=h_{i}(x),i\in p} x ∈ R n {\displaystyle x\in \mathbb {R} ^{n}} u ∈ R m {\displaystyle u\in \mathbb {R} ^{m}} y ∈ R p {\displaystyle y\in \mathbb {R} ^{p}} f , g , h {\displaystyle f,g,h}
観測空間をすべての繰り返し リー導関数を 含む空間と定義すると、 のときのみ、 系は で観測可能となる。 ここで O s {\displaystyle {\mathcal {O}}_{s}} x 0 {\displaystyle x_{0}} dim ( d O s ( x 0 ) ) = n {\displaystyle \dim(d{\mathcal {O}}_{s}(x_{0}))=n}
d O s ( x 0 ) = span ( d h 1 ( x 0 ) , … , d h p ( x 0 ) , d L v i L v i − 1 , … , L v 1 h j ( x 0 ) ) , j ∈ p , k = 1 , 2 , … . {\displaystyle d{\mathcal {O}}_{s}(x_{0})=\operatorname {span} (dh_{1}(x_{0}),\ldots ,dh_{p}(x_{0}),dL_{v_{i}}L_{v_{i-1}},\ldots ,L_{v_{1}}h_{j}(x_{0})),\ j\in p,k=1,2,\ldots .} [13] 非線形動的システムにおける観測可能性の初期の基準は、グリフィスとクマール [14] 、コウ、エリオットとターン [15] 、そしてシン [16] によって発見されました。
非線形時間変動システムに対する観測可能性基準も存在する。 [17]
静的システムと一般位相空間 観測可能性は、定常状態システム(通常は代数方程式と不等式で定義されるシステム)や、より一般的には、 における集合に対しても特徴付けられる 。 [18] [19]動的システムの場合、 カルマンフィルタ やその他の観測器 の挙動を予測するために観測可能性基準が使用されるのと同様に、 における集合の観測可能性基準は、 データ調整 やその他の静的推定器の挙動を予測するために使用される 。非線形システムの場合、観測可能性は個々の変数に対して特徴付けられるだけでなく、全体的な挙動だけでなく局所的な推定器の挙動に対しても特徴付けられる。 R n {\displaystyle \mathbb {R} ^{n}} R n {\displaystyle \mathbb {R} ^{n}}
参照
参考文献 ^ Kalman, RE (1960). 「制御システムの一般理論について」. IFAC Proceedings Volumes . 1 : 491– 502. doi :10.1016/S1474-6670(17)70094-8. ^ Kalman, RE (1963). 「線形動的システムの数学的記述」. Journal of the Society for Industrial and Applied Mathematics, Series A: Control . 1 (2): 152– 192. doi :10.1137/0301010. ^ abc ソンタグ, ED, 「数学的制御理論」, 応用数学テキスト, 1998 ^ 「制御可能性と観測可能性」 (PDF) 。 2019年6月10日時点の オリジナル (PDF)からアーカイブ 。 2024年5月19日 閲覧。 ^ Li, W.; Wei, G.; Ho, DWC; Ding, D. (2018年11月). 「センサーネットワークのための重み付け均一検出可能性」. IEEE Transactions on Neural Networks and Learning Systems . 29 (11): 5790– 5796. Bibcode :2018ITNNL..29.5790L. doi :10.1109/TNNLS.2018.2817244. PMID 29993845. S2CID 51615852. ^ Li, W.; Wang, Z.; Ho, DWC; Wei, G. (2019). 「カルマンコンセンサスフィルタリング問題における誤差共分散の有界性について」. IEEE Transactions on Automatic Control . 65 (6): 2654– 2661. doi :10.1109/TAC.2019.2942826. S2CID 204196474. ^ Fernando, T.; Trinh, HM; Jennings, L. (2010). 「機能的観測可能性と最小次数線形機能オブザーバーの設計」 IEEE Transactions on Automatic Control . 55 (5): 1268– 1273. doi :10.1109/TAC.2010.2042005. ^ Jennings, L.; Fernando, T.; Trinh, HM (2011). 「固有空間の観点から見た機能的観測可能性の存在条件」. IEEE Transactions on Automatic Control . 56 (12): 2957– 2961. doi :10.1109/TAC.2011.2162793. ^ Darouach, M. (2002). 「線形システムのための機能オブザーバーの存在と設計」. IEEE Transactions on Automatic Control . 45 (5): 940– 943. doi :10.1109/9.855563. ^ Montanari, AN; Duan, C.; Motter, AE (2025). 「ネットワークにおける目標制御と推定における制御性と観測性の双対性」. IEEE Transactions on Automatic Control . doi :10.1109/TAC.2025.3552001. ^ Brockett, Roger W. (1970). 有限次元線形システム . John Wiley & Sons. ISBN 978-0-471-10585-5 。 ^ Eduardo D. Sontag、「数理制御理論:決定論的有限次元システム」 ^ 教授による非線形システム理論の講義ノート博士。 D.ジェルツェマ 教授JMAScherpen と教授アジュヴァン・デル・シャフト。 ^ グリフィス, EW; クマール, KSP (1971). 「非線形システムの観測可能性について: I」. 数学解析応用ジャーナル . 35 : 135–147 . doi :10.1016/0022-247X(71)90241-1. ^ Kou, Shauying R.; Elliott, David L.; Tarn, Tzyh Jong (1973). 「非線形システムの観測可能性」. 情報制御 . 22 : 89–99 . doi : 10.1016/S0019-9958(73)90508-1 . ^ Singh, Sahjendra N. (1975). 「測定不可能な入力を持つ非線形システムにおける観測可能性」. International Journal of Systems Science . 6 (8): 723– 732. doi :10.1080/00207727508941856. ^ Martinelli, Agostino (2022). 「観測可能性ランク条件の時間変動非線形システムへの拡張」. IEEE Transactions on Automatic Control . 67 (9): 5002– 5008. Bibcode :2022ITAC...67.5002M. doi :10.1109/TAC.2022.3180771. ISSN 0018-9286. S2CID 251957578. ^ Stanley, GM; Mah, RSH (1981). 「プロセスデータ推定における観測可能性と冗長性」 (PDF) . 化学工学科学 . 36 (2): 259– 272. Bibcode :1981ChEnS..36..259S. doi :10.1016/0009-2509(81)85004-X. ^ Stanley, GM; Mah, RSH (1981). 「プロセスネットワークにおける観測可能性と冗長性の分類」 (PDF) . 化学工学科学 . 36 (12): 1941– 1954. doi :10.1016/0009-2509(81)80034-6.
外部リンク