状態量初期値および離散時間システムのシステム行列BとD
[[x0] [,B [,D]] [,V] [,rcnd]] = findBD(jobx0,comuse [,job],A [,B],C [,D],Y [,U,tol,printw,ldwork])
状態量初期値を計算するかどうかを指定する整数オプション:
1 : 状態量初期値 x0を計算;
2 : 状態量初期値を計算しない (x0 はゼロとする).
システム行列BおよびDを計算または使用するかどうかを指定する整数オプション:
1 : jobによる指定に基づき行列BおよびDを計算;
2 : jobによる指定に基づき行列BおよびDを使用;
3 : 行列BおよびDを計算も使用もしない.
システム行列のどちらを計算または使用するかを指定する整数オプション:
1 : システム行列Bのみを計算/使用 (D はゼロであるとする);
2 : 行列BおよびDを計算/使用.
jobx0 = 2 かつ comuse = 2, または comuse = 3 の場合, job を指定できません.
指定したシステムの状態行列
オプション, 指定したシステムの入力行列
指定したシステムの出力行列
オプション, 指定したシステムの直達項
t行l列の出力データ行列です. Yのj列目には,連続した時間増加のj番目の出力要素である t 個の値が含まれます.
t行m列の入力データ系列行列 (jobx0 = 1 かつ comuse = 2, または comuse = 1 の時に入力). Uのj列目には,連続する時間増分のj番目の入力要素の値が含まれます.
オプションで,行列のランクを推定する際に使用される許容誤差. tol > 0 の場合, tol に指定した値が, 条件数の逆数の下界として 使用されます. 推定される条件数が 1/tol より小さいm行n列の行列は フルランクであると判定されます. デフォルト: m*n*epsilon_machine ただし, epsilon_machine は相対機械精度です.
オプション, 警告メッセージの出力に関するスイッチ
1: 警告メッセージを出力;
0: 警告メッセージを出力しない.
デフォルト: printw = 0.
(オプション) ワークスペースの大きさ. デフォルト: 計算式 LDWORK = MAX( 必要なワークスペースの大きさの最小値, 2*CSIZE/3, CSIZE - ( m + l )*t - 2*n*( n + m + l ) - l*m ) ただし CSIZE はキャッシュの大きさ(単位:倍精度ワード)です.
状態ベクトル初期値
システム入力行列
システム直達行列
Aを実数Schur形式に縮小させるn行n列直交行列 (jobx0 = 1 または comuse = 1 の時に出力).
(オプション) ランク決定に関連する行列の条件数の逆数.
findBD 関数は, SLICOTルーチン IB01CD により, 状態量初期値および離散時間システムにおけるシステム行列 B と D を推定します.
注意: 上の例の各行では, 最後に パラメータ tol, printw, ldwork を含めることもできます.
FINDBD は状態量初期値および離散時間システムn システム行列 Br と Dr を推定します. この際,システム行列 A, C そして場合によっては B, D および システムにおける入力および出力の軌道を指定します.
モデルの構造は以下のようになります :
where x(k) は,(時刻 kにおける) n 次元状態ベクトルで,
u(k) はn次元入力ベクトル,
y(k) は l 次元出力ベクトル,
そして A, B, C, および D は適当な次元の実数行列です.
n行m列のシステム入力行列は, jobx0 = 1 かつ comuse = 2 の時に 入力パラメータ, comuse = 1 の時に出力パラメータとなります.
The l行m列のシステム行列 Dは, jobx0 = 1, comuse = 2 かつ job = 2 の時に 入力パラメータ, comuse = 1 かつ job = 2 の時に出力パラメータとなります.
n次ベクトルである状態量初期値の推定値 x(0) は, jobx0 = 1 の時, および jobx0 = 2 かつ comuse <= 2 の時 に出力パラメータとなり, この場合, 0 に設定されます.
ldwork が指定されたが, 必要なワークスペースの大きさの最小値よりも 小さい場合, かわりに最小値が使用されます.
// 指定した線形システムからデータを生成 A = [ 0.5, 0.1,-0.1, 0.2; 0.1, 0, -0.1,-0.1; -0.4,-0.6,-0.7,-0.1; 0.8, 0, -0.6,-0.6]; B = [0.8;0.1;1;-1]; C = [1 2 -1 0]; SYS=syslin(0.1,A,B,C); nsmp=100; U=prbs_a(nsmp,nsmp/5); Y=(flts(U,SYS)+0.3*rand(1,nsmp,'normal')); // Rを計算 S=15;L=1; [R,N,SVAL] = findR(S,Y',U'); N=3; METH=3;TOL=-1; [A,C] = findAC(S,N,L,R,METH,TOL); [X0,B,D] = findBD(1,1,2,A,C,Y',U') SYS1=syslin(1,A,B,C,D,X0); Y1=flts(U,SYS1); clf();plot2d((1:nsmp)',[Y',Y1']) | ![]() | ![]() |