制御工学ブログ

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

状態空間モデルと応答

この記事では状態空間モデルと制御系の応答についてまとめます。関連動画は以下のものです。また、状態方程式について説明した関連記事リンクは最下部に置いています。

 

youtu.be

執筆者情報:岡島 寛 (熊本大学工学部情報電気工学科准教授,Web, YouTube)約20年教員をやっています。モデル誤差抑制補償器,状態推定,量子化制御など制御工学の研究をしています。

状態空間モデルと次数

状態空間表現の全体像

制御対象の動特性を表現する方法として状態方程式があります。説明は次の記事「状態方程式の概要」で行っております。次数 nが大きいとより複雑な動特性を表現することができます。ここでは、 n=1,  n=2,  n=3の3パターンの制御対象を説明します。次の図が n = 1のケースです。

n = 1の場合

 Aはスカラであり、 B, Cも1入力である場合にはスカラになります。教科書的な説明では係数として整数を用いることが多いですが、実際は小数点以下の値を持つような実数が係数として与えられることになります。他方、状態が何等かの物理量で出力も物理量の場合には C C = 1となることもあります。

次に、 n = 2の場合を考えます。

n = 2の場合

状態数は 2であり、行列 Aは2×2の行列となります。これらを当てはめると行列ベクトル形式の状態空間表現が求まりますが、個々の状態量に分解して x_1および x_2のそれぞれのダイナミクスを書き下すこともできます。

 n = 3の場合は以下の通りです。

n = 3の場合

状態数は 3であり、行列 Aのサイズは3×3となっています。

状態空間モデルの入出力数

次に、入力数が 2の場合を示します。

入力数が2の場合

このとき、行列 Bが3行2列の行列となっています。また、入力はスカラではなくベクトルとして2行1列の行列となります。

最後に、2入力2出力の多入出力系の状態空間モデルを示します。

2入力2出力系

行列 Cが2行3列の行列になっていることがわかります。出力 yの要素を y_1, y_2とするとそれらが縦に並ぶ2行1列の形式で出力信号が表現されます。

数値シミュレーション例

数値シミュレーション条件

パラメータ値が異なる9種類の制御対象について、応答の比較を行います。まず、インパルス応答を比較します。全ての図において、対象1を基準として与え、それと対象2,3、対象4,5、対象6,7、対象8,9を比較してます。

インパルス応答

次に、ステップ応答を示します。

ステップ応答

ここでは、意識的に対象の極の変更、零点の変更を行っています。それぞれ、同じ入力信号を印加しても応答波形が全く違うことがわかります。対象のパラメータ A,B,Cが適切に与えられていないと、対象の特性を表現できません。また、対象ごとに制御のしやすさやしにくさは変わってきます。

極は零点については以下の記事を参考にしてください。

 

blog.control-theory.com

blog.control-theory.com

blog.control-theory.com

関連記事

以下は状態方程式について説明している記事になります。

blog.control-theory.com

blog.control-theory.com

本記事は以上です。役に立った、と思われましたら、ブックマーク・シェア等のアクションをしていただければ嬉しいです。

 

制御工学の研究者紹介:制御工学チャンネル:500本以上の制御動画ポータルサイト - 日本の制御研究者 (control-theory.com)

制御の書籍紹介:制御工学チャンネル:500本以上の制御動画ポータルサイト - 制御工学関連の書籍 (control-theory.com)

制御の書籍購入:Amazon.co.jp: 制御工学Amazonアソシエイト

制御キットの購入:Amazon.co.jp: 制御キットAmazonアソシエイト