isproper

R2026b

動的システムモデルがプロパーであるかどうかの判定

構文

B = isproper(sys) B = isproper(sys,"elem") [B,sysr,info] = isproper(sys) [B,sysr,info] = isproper(sys,info)

説明

B = isproper(sys) は、動的システムモデル sys がプロパーな場合に 1 (true) の論理値を返し、それ以外の場合は 0 (false) の論理値を返します。

プロパーなモデルは相対次数 ≤ 0 であり、因果性があります。SISO 伝達関数および零点-極-ゲインモデルは、分子係数の次数が分母係数の次数以下である場合 (すなわち、それらが少なくとも零点と同じ個数の極をもっている場合) にプロパーになります。SISO エントリがすべてプロパーであれば、MIMO 伝達関数はプロパーです。通常の状態空間モデル (行列 E をもたない状態空間モデル) は常にプロパーです。可逆の行列 E をもつ記述子状態空間モデルは常にプロパーです。特異な (可逆でない) 行列 E をもつ記述子状態空間モデルは、モデルが少なくとも零点と同じ個数の極をもっている場合にプロパーです。

sys がモデル配列の場合、配列内のすべてのモデルがプロパーであれば B は 1 になります。

B = isproper(sys,"elem") はモデル配列 sys 内の各モデルをチェックし、sys と同じサイズの logical 配列を返します。logical 配列は sys のどのモデルがプロパーであるかを示します。dsys がモデル配列ではなく単一のモデルである場合、"elem" 引数は効果がありません。

[B,sysr,info] = isproper(sys) は単一のモデル sys を受け取り、同等のモデル sysr も返します。sys がプロパーでない場合は、sysr = sys になります。sys が、非可逆行列 E をもつプロパーな記述子状態空間モデルである場合、sysr はそれより少ない状態 (低い次数) になり、非特異行列 E をもちます。info 構造体には、変換が反映された射影行列が含まれます。射影行列の詳細については、Spectral Projector Matricesを参照してください。info の詳細については、dss2ss と sminreal を参照してください。

[B,sysr,info] = isproper(sys,info) は、変換を通じて伝播する情報が含まれる info 行列または構造体も受け取ります。最後の出力として返される info 構造体には、更新された初期条件と、変換を反映して更新された射影行列が含まれます。

例

すべて折りたたむ

モデルがプロパーかどうかの検証

ライブスクリプトを開く

SISO 連続時間伝達関数 $H_{1} = s$ を作成します。

H1 = tf([1 0],1);

H1 がプロパーかどうかをチェックします。

B1 = isproper(H1)

B1 = logical
   0

SISO 伝達関数は、分子の次数が分母の次数以下である場合にプロパーになります。すなわち、伝達関数が少なくとも零点と同じ個数の極をもっている場合です。H1 は 1 つの零点をもち、極をもたないので、isproper コマンドは 0 を返します。

次に、極と零点を 1 つずつもつ伝達関数 $H_{2} = s / (s + 1)$ を作成します。

H2 = tf([1 0],[1 1]);

H2 がプロパーかどうかをチェックします。

B2 = isproper(H2)

B2 = logical
   1

H2 がもつ極と零点の数は同じなので、isproper は 1 を返します。

等価の低次元化されたモデルの計算

ライブスクリプトを開く

状態空間モデルを組み合わせると、必要以上の状態を含む結果が生成されることがあります。isproper を使用して、等価の低次元化されたモデルを計算します。

H1 = ss(tf([1 1],[1 2 5]));
H2 = ss(tf([1 7],[1]));
H = H1*H2;
size(H)

State-space model with 1 outputs, 1 inputs, and 4 states.

H はプロパーで可約です。isproper は低次元化されたモデルを返します。

[isprop,Hr] = isproper(H);
size(Hr)

State-space model with 1 outputs, 1 inputs, and 2 states.

H と Hr は、ボード線図が示すように等価です。

bodeplot(H,Hr,'r--')
legend('original','reduced')

MATLAB figure

参照

[1] Varga, Andràs. "Computation of irreducible generalized state-space realizations." Kybernetika 26.2 (1990): 89-106.

バージョン履歴

R2006a より前に導入

参考

ss | dss