sorder

パラメータ

meth

使用する手法を定義するための整数オプション:

=

1 : 過去の入力および出力を用いるMOESP法;

=

2 : N4SID 法.

alg

入出力データから構築される 結合されたブロックハンケル行列の三角分解を計算するアルゴリズムを 定義するための整数オプション:

=

1 : 相関行列のコレスキーアルゴリズム;

=

2 : 高速 QR アルゴリズム;

=

3 : 標準 QR アルゴリズム.

jobd

行列BおよびDを後でMOESP法により計算するかどうかを指定する 整数オプション:

=

1 : 行列 B および D を後でMOESP法により計算します;

=

2 : 行列 B および D はMOESP法により計算されません.

このパラメータは mesh = 2 には関連しません.

batch

連続的なデータ処理を行うかどうか, また,連続処理の際に,カレントのデータブロックを先頭ブロック, 中間ブロック,または最終ブロックとするかどうかを以下のように 指定する整数オプション:

=

1 : 連続データ処理で先頭ブロック;

=

2 : 連続データ処理で中間ブロック;

=

3 : 連続データ処理で最終ブロック;

=

4 : 1ブロックのみ (連続データ処理をしない).

conct

連続データ処理における連続するデータブロックを1回の試行に配置する かどうかを以下のように指定する整数オプション:

=

1 : カレントのデータブロックの処理は, 前のデータブロックと連続し,また, 次のブロックも継続処理されます.;

=

2 : カレントのデータブロックと前後のデータブロックの間は接続されません.

このパラメータは, batch = 4 には使用されません.

s

処理されるブロックハンケル行列の入力/出力のブロック行数. s > 0

Y

the t行l列の出力データ系列行列. Yのj列には, j番目の出力要素のtの時系列の値が含まれます.

U

(オプション) t行m列の入力データ系列行列. Column j of Uのj列には, j番目の入力要素のtの時系列の値が含まれます.デフォルト: U = [].

tol

(オプション) 許容誤差を含む長さ2のベクトル: tol(1) - 行列のランクを推定する際に使用される許容誤差. tol(1) > 0の場合, 条件数の逆数の下界として,指定したtol(1)の値が使用されます; m行n列の行列で,推定された条件数が 1/tol(1) より小さい場合,フルランクとみなされます. tol(1) <= 0 の場合, デフォルト値 m*n*epsilon_machineが使用されます. ただし,epsilon_machine は相対機械精度です.

tol(2) - システム次数の推定値を定義する際に使用される許容誤差. tol(2) >= 0 の場合, tol(2)以上の最後の特異値の添字が推定値となります. (tol(2)よち小さな特異値はゼロとみなされます.) tol(2) = 0の場合,デフォルト値, tol(2) = s*epsilon_machine*sval(1), が内部的に計算され,使用されます. ただし,sval(1)は最大特異値,epsilon_machineは相対機械精度です. tol(2) < 0 の場合, 推定値は次の要素への指数的な差異が最大となる特異値の添字となります.

デフォルト: tol(1:2) = [0,-1].

printw

(オプション) 警告メッセージの出力に関するスイッチ.

=

1: 警告メッセージを出力します;

=

0: 警告メッセージを出力しません.

デフォルト: printw = 0.

ldwork

(オプション) ワークスペースの大きさ. デフォルト : 以下の式で計算

nr = 2*( m + l )*s
LDWORK = ( t - 2*s + 3 + 64 )*nr
if ( CSIZE > MAX( nr*nr + t*( m + l ) + 16, 2*nr ) ) then
   LDWORK = MIN( LDWORK, CSIZE - nr*nr - t*( m + l ) - 16 )
else
   LDWORK = MIN( LDWORK, MAX( 2*nr, CSIZE/2 ) )
end if

LDWORK = MAX( minimum workspace size needed, LDWORK ) ただし,CSIZE は倍精度ワード単位のキャッシュの大きさです.

最小ワークスペースの大きさよりも小さなLDWORK が指定された場合, 替わりに最小値が使用されます.

Ri

(オプション) batch = 2 または 3の場合, Rの 2*(m+l)*s行2*(m+l)*s列の部分 (alg <> 2の場合は,上三角) に 連続データ処理におけるこの mexfile の前のコールで計算された (上三角)行列 R を含む必要があります. conct = 1の場合, Rには同じく前のコールで設定された列が1列追加されます.

alg = 2の場合, R は m+l+1 個の前のコールで設定された追加の列を有します.

このパラメータは, batch = 1 または batch = 4の場合は使用されません.

Ro

batch = 3 または 4の場合, Rの2*(m+l)*s行2*(m+l)*s列 の部分には, 結合されたブロックハンケル行列のQR分解に基づく 上三角分解 Rが含まれ, システム行列を計算する際により詳細な情報が必要となります. batch = 1 または 2の場合, R にはこのmexfileの次のコールで必要とされる 中間結果が含まれます. batch = 1 または 2, かつ conct = 1の場合, R には列が追加され, リターンの前に値が設定されます. batch = 1 または 2 で alg = 2の場合,R には m+l+1個の列が追加され, 同じくリターンの前に値が設定されます.

n

システムの次数.

sval

(オプション) システムの次数を推定する際に使用される特異値.

rcnd

(オプション) meth = 2の場合, 長さ2のベクトルで, ランク決定や最小二乗解で使用される 行列の条件数の逆数を含みます.

説明

sorder - SLICOT の IB01ADルーチンにより離散時間システムの次数を計算する関数.

ブロックが1個の場合 (データ系列 Y, U): [R,n,sval,rcnd] = sorder(meth,alg,jobd,4,conct,s,Y,U);

ブロックがf個の場合 (データ系列 Yj, Uj, j = 1 : f):

R = sorder(meth,alg,jobd,1,conct,s,Y1,U1);
for j = 2 : f - 1
   R = sorder(meth,alg,jobd,2,conct,s,Yj,Uj,tol,printw,ldwork,R)
end
[R,n,sval,rcnd] = sorder(meth,alg,jobd,3,conct,s,Yf,Uf,tol);

sorder は,線形時不変システムの行列を推定するために コレスキーまたは(高速)QR分解および 部分空間同定法 (MOESP およびN4SID)による 入出力データの前処理を行い, 離散時間実現の次数を推定します.

モデルの構造は以下のようになります :

x(k+1) = Ax(k) + Bu(k) + w(k),   k >= 1,
y(k)   = Cx(k) + Du(k) + e(k),

ただし, x(k) は(時刻kにおける) n次元状態ベクトル,

u(k) はm次元入力ベクトル,

y(k) はl次元出力ベクトル,

w(k) はn次元状態外乱ベクトル,

e(k) はl次元出力外乱ベクトル,

として A, B, C, および D は適当な次元の実数行列です.

コメント

1. コレスキーまたは高速QRアルゴリズムは,(大きなデータブロックの場合) QRアルゴリズムよりも はるかに高速ですが, 修正行列 H'*H が正定でない場合には使用できません. このような場合, 十分なワークスペースが確保され, batch = 4 の時に コードは自動的に QR アルゴリズムに切り替えます.

2. ldwork が指定されたが, ワークスペースの大きさの最小値よりも小さい場合, 替わりに最小値が使用されます.

参照

• findBD — 状態量初期値および離散時間システムのシステム行列BとD
• sident — 離散時間状態空間実現およびカルマンゲイン