制御工学ブログ

制御工学の研究者を20年やっている国立大学教員が制御工学の基礎から専門まで広く説明します。記事内では、動画やMATLABコードを交えながらわかりやすく解説する方針です。制御工学チャンネル(YouTube,動画ポータル)を運営しています。

システムの可制御性・可観測性について

この記事ではシステムの可制御性・可観測性についてまとめます。まず,可制御性について述べたいと思います。なお,可制御性・可観測性について説明した動画へのリンクは最下部に置いています。

なお、状態フィードバック制御の全体像は次の記事でまとめています。

状態フィードバック制御・状態方程式に基づく制御のまとめ

状態方程式表現

状態方程式の

状態方程式の記号

まず制御対象の状態方程式表現は上図のようになっています。まず状態方程式において、制御入力は u、制御出力は  y です。状態ベクトル x と表現していて、その要素数 n 個あります。

このようなシステムの次数のことをnという記号を用いて表していて、 n 次のシステムに対して A 行列は  n 次の正方行列です。 Bは仮に1入力システムの場合には  n次の縦長ベクトルで、入力数が多数ある場合は行列の形になります。

 C n 次の横ベクトルで、仮に出力数が多数ある場合には、行列の形になります。

システムの可制御性

まず可制御から説明をします。可制御性は、入力によって自由に状態を動かせるかを表現する指標であり、制御工学分野の重要事項の一つです。まず、可制御か可安定であれば状態フィードバックで安定が可能です。例えば、制御対象に対して、 u=- Kx という状態フィードバック制御を施す場合、これによって

\begin{equation} \dot x (t) = (A-BK)x(t) \end{equation}

で与えられる自律システムが得られます。この  A-BK固有値の実部が負であれば自律システムは安定になります。もし、システムが可制御もしくは可安定であるならば、そのようなことを実現する Kが存在するといえます。

ここで可制御であれば可安定であるというのは成り立ちますが、逆に可安定であっても
場合によっては可制御にならないようなケースがあります。

例えばもともと Aが安定で、 Bが零ベクトルというケースでは、 Kを適当に与えてもシステムの安定性は保たれますが可制御ではなく、極を任意に配置するといったこともできません。

可制御であれば、適切なフィードバックゲイン Kを与えることで (A-BK) の極を任意に配置することができます。

可制御性行列と可制御性判定の図

可制御性の判定

図に示すように、ここで可制御性を調べるためには、可制御性行列Uc というものを用います。この行列 U_c B AB、…と並べた行列になっています。仮にこの行列のランクが nであれば、このシステムは可制御になります。

次数ごとの可制御性行列

まず、 n = 1 の場合は状態数 1 であり、 U_cは次の形になります。

\begin{equation} U_c = B\end{equation}

このとき、 Bが非零であれば  n = 1のシステムは可制御となります。

次に、 n = 2 の場合です。状態数は2であり、 U_cは次の形になります。

可制御性行列02

最後に、 n = 4 の場合です。 U_cは以下で与えられます。

可制御性行列03

入力数が1の場合は、行列サイズが  n \times nになります。そのランクが  n であれば可制御で、そうでなければ不可制御です。

次に、具体例を見ていきましょう。

可制御性判別の例を3つ表示する図

可制御性判別の例

例1: まず  n が2のシステムとしてA, B行列が図の左のように与えられている2入力のシステムを考えます。 n = 2の2次システムですので、可制御性行列 U_cのランクが2であることから可制御になります。入力数が複数になると行列が正方でなく横長になります。

 

例2: 次の例としては  n =3としてA, Bが図の中央のように与えられている例を考えます。このケースでは、システムの次数は 3 です。ここで、可制御性行列のランクが 3 ありますので可制御となります。

なお、この場合の A, B行列の対は可制御正準系と呼ばれる形式になっています。


例3:  n =3のシステムとして A, Bが図の右のように与えられているケースを考えます。最後にこの3番目のケースですが、こちらは  U_c のランクは2になります。そのため n =3であるのに対し、ランクが2であるため不可制御となります。

状態方程式可観測性

次に可観測について説明します。可観測もしくは可検出であれば制御対象の全状態を状態オブザーバーによって推定することが可能です。ここで先ほどの可制御の場合と同様に可観測であれば、可検出であるという性質が成り立ちます。その逆に可検出であっても必ずしも可観測は成り立たないといえます。

可観測性を判別するための方法を示した図

可観測性の判別

ここで可観測性を調べるためには、可観測性行列を調べれば良いです。この可観測
行列は  C,CA…と縦に並べた行列になります。

次数ごとの可観測性行列

まず、 n = 1 の場合は状態数 1 であり、 U_oは次の形になります。

可観測性行列01

このとき、 Cが非零であれば可観測になります。

次に、 n = 2 の場合です。状態数は2であり、 U_oは次の形になります。

可観測性行列02

最後に、 n = 3 の場合です。状態数は3であり、 U_oは次の形になります。

可観測性行列03

この行列のランクが nであれば、可観測  n 未満であれば不可観測となります。ここでは示しませんが、可制御なシステムで可制御正準形があるように、可観測なシステムでは可観測正準形があります。

可制御性行列の場合と同様に、出力数が1であれば正方行列となりますが、出力数が2以上であると縦長の行列になります。

可制御・可観測に関する動画

最後に、可観測性に関する動画を示して終わりたいと思います。

youtu.be

youtu.be

 

関連ページ(可制御・可観測)

>>制御工学チャンネル:500本以上の制御動画ポータルサイト Control Engineering Channel - 01-Prof.Okajima (control-theory.com)

前の記事:

 

blog.control-theory.com

自己紹介

岡島 寛 (熊本大学工学部情報電気工学科准教授)

制御工学の研究をしています。モデル誤差抑制補償器,状態推定,量子化制御など

研究室HP

岡島研究室(システム制御 control-theory.com)

制御動画ポータルサイト

制御工学チャンネル(伝達関数・状態方程式・MATLABなど)

電気動画ポータルサイト

電気電子チャンネル (半導体・電気・電子工作など)

YouTube 

制御工学チャンネル (制御YouTubeチャンネル)

 

本記事をお読みいただきありがとうございます。役に立った、と思われましたら、ブックマーク・シェア等のアクションをしていただければ嬉しいです。