ctrb

状態空間モデルの可制御性

構文

Co = ctrb(A,B)

Co = ctrb(sys)

説明

動的システムを駆動する制御信号を有限時間内に任意の状態に適用することが可能な場合、その動的システムは "可制御" と言われます。この特性は "可到達性" とも呼ばれます。ctrb は、可制御性行列を状態行列から、または状態空間モデルから計算します。この行列を使用して、可制御性を判定することができます。

たとえば、Nx 個の状態、Ny 個の出力、および Nu 個の入力をもつ連続時間状態空間モデルを考えます。

$\begin{array}{l} \dot{x} = A x + B u \\ y = C x + D u \end{array}$

ここで、x、u、および y はそれぞれ状態、入力、および出力を表し、A、B、C、および D は以下のサイズをもつ状態空間行列です。

A は Nx 行 Nx 列の実数値行列または複素数値行列。
B は Nx 行 Nu 列の実数値行列または複素数値行列。
C は Ny 行 Nx 列の実数値行列または複素数値行列。
D は Ny 行 Nu 列の実数値行列または複素数値行列。

ctrb により生成される可制御性行列 $C o = [\begin{matrix} B & A B & A^{2} B & \dots & A^{n - 1} B \end{matrix}]$ がフルランクの場合、すなわち、ランクが状態空間モデルの状態数と等しい場合、システムは可制御です。可制御性行列 Co は、Nx 行 Nxu 列です。例については、SISO 状態空間モデルの可制御性を参照してください。

Co = ctrb(A,B) は、状態行列 A と入力から状態への行列 B を使用して、可制御性行列 Co を返します。Co がフルランクの場合、すなわち、Co のランクが状態数と等しい場合、システムは可制御です。

例

Co = ctrb(sys) は、状態空間モデル sys の可制御性行列を返します。この構文は、次と等価です。

Co = ctrb(sys.A,sys.B);

例

すべて折りたたむ

システムの可制御性のチェック

ライブスクリプトを開く

行列 A および B を定義します。

A = [1  1;
     4 -2];
B = [1 -1;
     1 -1];

可制御性行列を計算します。

Co = ctrb(A,B);

非可制御状態の数を特定します。

unco = length(A) - rank(Co)

unco = 
1

非可制御状態は Co がフルランク 2 ではないことを示します。したがって、システムは可制御ではありません。

SISO 状態空間モデルの可制御性

ライブスクリプトを開く

この例では、2 つの状態をもつ次の SISO 状態空間モデルを考えます。

$A = [\begin{array}{cc} - 1.5 & - 2 \\ 1 & 0 \end{array}] B = [\begin{array}{c} 0.5 \\ 0 \end{array}] C = [\begin{array}{cc} 0 & 1 \end{array}] D = 1$

次の状態空間行列で定義される SISO 状態空間モデルを作成します。

A = [-1.5,-2;1,0];
B = [0.5;0];
C = [0,1];
D = 1;
sys = ss(A,B,C,D);

可制御性行列を計算し、ランクを求めます。

Co = ctrb(sys)

Co = 2×2

    0.5000   -0.7500
         0    0.5000

可制御性行列のサイズは、行列 A および行列 B のサイズに依存します。たとえば、行列 A が Nx 行 Nx 列の行列であり、行列 B が Nx 行 Nu 列の行列である場合、結果の行列 Co は Nx 行 Nxu 列になります。ここで、Nx は状態の数、Nu は入力の数です。

rank(Co)

ans = 
2

可制御性行列 Co のランクが状態の数と等しいため、システム sys は可制御です。

あるいは、行列 A および行列 B のみを使用して可制御性行列を求めることもできます。

Co = ctrb(sys.A,sys.B);
rank(Co)

ans = 
2

入力引数

すべて折りたたむ

`A` — 状態行列
`Nx` 行 `Nx` 列の行列

状態行列。Nx 行 Nx 列の行列として指定します。ここで Nx は状態の数です。

`B` — 入力から状態への行列
`Nx` 行 `Nu` 列の行列

入力から状態への行列。Nx 行 Nu 列の行列として指定します。ここで Nx は状態の数、Nu は入力の数です。

`sys` — 状態空間モデルまたは状態空間モデル配列
`ss` モデルオブジェクト | `genss` モデルオブジェクト | `uss` モデルオブジェクト

状態空間モデルまたは状態空間モデル配列。以下として指定します。

状態空間 (ss) モデルオブジェクト (入力 A、B、C および D が数値行列の場合、または別のモデルオブジェクトタイプから変換する場合)。
一般化状態空間モデル (genss) オブジェクト (1 つ以上の行列 A、B、C および D に realp パラメーターや一般化行列 (genmat) などの調整可能なパラメーターが含まれている場合)。関数は調整可能なパラメーターの現在の値を使用します。
不確かさをもつ状態空間モデル (uss) オブジェクト (1 つ以上の入力 A、B、C および D に不確かさをもつ行列が含まれている場合)。関数は不確かさをもつパラメーターのノミナル値を使用します。不確かさをもつモデルを使用するには Robust Control Toolbox™ ソフトウェアが必要です。

出力引数

すべて折りたたむ

`Co` — 可制御性行列
配列

可制御性行列。配列として返されます。sys の実体と、対応する配列は次のとおりです。

Nx 個の状態と Nu 個の入力をもつ単一の状態空間モデルである場合、結果の配列 Co は Nx 行 Nxu 列になります。
状態空間モデルの配列 sys(:,:,j1,...,jN) である場合、Co は N+2 個の次元をもつ配列、すなわち Co(:,:,j1,...,jN) になります。

参照

[1] Paige, C. C. "Properties of Numerical Algorithms Related to Computing Controllability." IEEE Transactions on Automatic Control. Vol. 26, Number 1, 1981, pp. 130-138.