結論
可制御な線形システムを座標変換しても、可制御だということを述べる。
つまり:
線形システム\begin{eqnarray}\dot{x}(t) &=& Ax(t)+Bu(t) \\ y(t)&=&Cx(t) + Du(t) \end{eqnarray}
の状態 \(x(t)\) を別の状態 \(z(t)\) に \(x(t) = Tz(t)\) (\(T\) は正則行列)で座標変換したとき、
システムの可制御性は保存する。
※状態空間表現についてはこちらをご覧ください。
※可制御性についてはこちらをご覧ください。
わかること:
解析しやすくするため、システムを座標変換したとき、
「急に制御できなくなる事態は発生しない」と保証できる。
例
モータの角度、角速度が計測できる直流モータの状態空間表現は
\begin{eqnarray}&\frac{d}{dt}\begin{pmatrix} x_1 (t) \\ x_2 (t)\end{pmatrix}&=\begin{pmatrix} 0 & 1 \\ 0 & -\frac{1}{T_m}\end{pmatrix}\begin{pmatrix} x_1 (t) \\ x_2 (t)\end{pmatrix}+\begin{pmatrix} 0 \\ \frac{1}{T_m K_E}\end{pmatrix}u(t) \\ &\begin{pmatrix} y_1 (t) \\ y_2 (t)\end{pmatrix}&= \begin{pmatrix} c_{\theta} & 0 \\ 0 & c_{\omega}\end{pmatrix}\begin{pmatrix} x_1 (t) \\ x_2 (t)\end{pmatrix}\end{eqnarray}とかける。※詳細はこちらをご覧ください。
このシステムは \(T_m, K_E \neq 0,\infty\) であるため
\[{\rm rank} \begin{pmatrix}B & AB\end{pmatrix} = {\rm rank}\begin{pmatrix}0 & \frac{1}{T_m K_E} \\ \frac{1}{T_m K_E} & -\frac{1}{T_m^2 K_E} \end{pmatrix} = 2\]
より可制御である。
次にこのモータのモデルを座標変換する。ここでは
\[\begin{pmatrix} x_1 (t) \\ x_2 (t)\end{pmatrix} = \underbrace{\begin{pmatrix} 1 & 1 \\ 0 & -\frac{1}{T_m}\end{pmatrix}}_{\text{\(T\)}} \begin{pmatrix} z_1 (t) \\ z_2 (t)\end{pmatrix}\]
とする。
※この正則行列 \(T\) は適切に選ぶ。さまざまな方法が提案されている。
座標変換を代入・整理・計算すると、モータのモデルは以下のように書き直すことができる。
\begin{eqnarray}&\frac{d}{dt}\begin{pmatrix} z_1 (t) \\ z_2 (t)\end{pmatrix}&=\begin{pmatrix} 0 & 0 \\ 0 & -\frac{1}{T_m}\end{pmatrix}\begin{pmatrix} z_1 (t) \\ z_2 (t)\end{pmatrix}+\begin{pmatrix} \frac{1}{K_E} \\ -\frac{1}{K_E}\end{pmatrix}u(t) \\ &\begin{pmatrix} y_1 (t) \\ y_2 (t)\end{pmatrix}&= \begin{pmatrix} c_{\theta} & c_{\theta} \\ 0 & -\frac{c_{\omega}}{T_m}\end{pmatrix}\begin{pmatrix} z_1 (t) \\ z_2 (t)\end{pmatrix}\end{eqnarray}
この変換後のシステムの可制御性を判定する。
\[{\rm rank} \begin{pmatrix}B & AB\end{pmatrix} = {\rm rank}\begin{pmatrix}\frac{1}{K_E} & 0 \\ -\frac{1}{K_E} & \frac{1}{T_m K_E} \end{pmatrix} = 2\]
よって、変換後のシステムも可制御。
座標変換しても、(同じモノを表現しているため)可制御という性質は変わらない。
証明
「元のシステムが可制御ならば、座標変換後のシステムは可制御」であることを示す。
つまり:\begin{eqnarray}\dot{x}(t) &=& Ax(t)+Bu(t) \\ y(t)&=&Cx(t) + Du(t) \end{eqnarray}について
\begin{eqnarray}{\rm rank} \begin{pmatrix}B & AB & \cdots & A^{n-1}B \end{pmatrix}=n \tag{Assumption}\end{eqnarray}が成立するとき、座標変換後のシステム
\begin{eqnarray}\dot{z}(t) &=& \underbrace{T^{-1}AT}_{\text{状態 \(z\) での行列 \(A\)}} z(t)+\underbrace{T^{-1}B}_{\text{状態 \(z\) での行列 \(B\)}}u(t) \\ y(t)&=&\underbrace{CT}_{\text{状態 \(z\) での行列 \(C\)}} z(t) + \underbrace{D}_{\text{状態 \(z\) での行列 \(D\)}}u(t) \end{eqnarray}の可制御の条件
\begin{eqnarray}{\rm rank} \begin{pmatrix}T^{-1}B & (T^{-1}AT)(T^{-1}B) & \cdots & (T^{-1}AT)^{n-1}(T^{-1}B) \end{pmatrix}=n \tag{Conclusion}\end{eqnarray}
が成り立つことを示す。
(Conclusion) の左辺をさらに式変形すると
\begin{eqnarray}{\rm rank} \begin{pmatrix}T^{-1}B & T^{-1}AB & \cdots & T^{-1}A^{n-1}B \end{pmatrix} \end{eqnarray}
ランクの性質より
\begin{eqnarray}{\rm rank}\left(\underbrace{T^{-1}}_{\text{変換行列より正則.}} \underbrace{\begin{pmatrix}B & AB & \cdots & A^{n-1}B \end{pmatrix}}_{\text{(Assumption) より rank \(=n\).}}\right)=n \end{eqnarray}
よって座標変換後の可制御性行列のランクも \(n\) となる。
したがって可制御なシステムを座標変換しても、可制御のまま。
参考文献
[1] 梶原 “線形システム制御入門,” コロナ社, 2000.
北海道大学大学院情報科学研究科修士課程修了。
機械メーカにて開発業務に従事したのち、フリーのエンジニア・講師として活動中。