13
13
<refentry xmlns="http://docbook.org/ns/docbook" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:svg="http://www.w3.org/2000/svg" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:db="http://docbook.org/ns/docbook" version="5.0-subset Scilab" xml:lang="ja" xml:id="findBDK">
15
<refname>findBDK</refname>
16
<refpurpose>カルマンゲインおよび離散時間システムの B, D システム行列</refpurpose>
20
<synopsis>[B,D,K] = findBDK(S,N,L,R,A,C,METH,JOB,NSMPL,TOL,PRINTW)
21
[B,D,RCND] = findBDK(S,N,L,R,A,C,METH,JOB)
22
[B,D,K,Q,Ry,S,RCND] = findBDK(S,N,L,R,A,C,METH,JOB,NSMPL,TOL,PRINTW)
31
<para>整数, ブロックハンケル行列のブロック行の数</para>
49
<para>行列, findR をコールすることにより計算される結合されたブロックハンケル行列の
69
<para>整数, 使用する方法に関するオプション</para>
74
<para> 過去の入力と出力を指定するMOESP法;</para>
80
<para> N4SID 法;</para>
92
<para>計算するシステム行列の種類を指定するオプション:</para>
103
<para> BおよびD行列を計算.</para>
116
共分散行列およびカルマン予測ゲインを計算する際に使用される
118
このパラメータは,共分散行列とカルマン予測ゲインが必要ない
119
場合には,不要です.NSMPL = 0の場合, K, Q, Ry, および S は計算されません.
128
行列のランクを推定する際に使用される許容誤差.
129
TOL > 0 の場合, TOL に指定した値が, 条件数の逆数の下界として
131
デフォルト: prod(size(matrix))*epsilon_machine ただし,
132
epsilon_machine は相対機械精度です.
139
<para>整数, 警告メッセージの出力に関するスイッチ</para>
144
<para>= 1: 警告メッセージを出力;</para>
150
<para>= 0: 警告メッセージを出力しない.</para>
162
<para>状態空間実現 SYS = (A,B,C,D) (syslin オブジェクト)を計算</para>
168
<para>カルマン予測ゲイン K (NSMPL > 0の場合)</para>
186
<para>状態出力相互共分散</para>
194
ランク決定, 最小二乗またはリカッチ方程式の解に関わる行列の条件数の逆数が含まれます.
203
システムの次数,行列AおよびCおよび関連する
204
結合ブロックハンケル行列のR因子を指定すると,
205
部分空間同定法 (MOESP または N4SID) により,
206
離散時間システムのシステム行列BおよびDとカルマンゲインを返します.
210
<para>[B,D,K] = findBDK(S,N,L,R,A,C,METH,JOB,NSMPL,TOL,PRINTW) は,
211
システム行列 B (JOB = 1の場合),
212
システム行列 BおよびD (JOB = 2の場合),
213
カルマン予測ゲイン K (NSMPL > 0の場合) を計算します.
15
<refname>findBDK</refname>
16
<refpurpose>カルマンゲインおよび離散時間システムの B, D システム行列</refpurpose>
20
<synopsis>[B,D,K] = findBDK(S,N,L,R,A,C,METH,JOB,NSMPL,TOL,PRINTW)
21
[B,D,RCND] = findBDK(S,N,L,R,A,C,METH,JOB)
22
[B,D,K,Q,Ry,S,RCND] = findBDK(S,N,L,R,A,C,METH,JOB,NSMPL,TOL,PRINTW)
31
<para>整数, ブロックハンケル行列のブロック行の数</para>
49
<para>行列, findR をコールすることにより計算される結合されたブロックハンケル行列の
69
<para>整数, 使用する方法に関するオプション</para>
74
<para> 過去の入力と出力を指定するMOESP法;</para>
80
<para> N4SID 法;</para>
92
<para>計算するシステム行列の種類を指定するオプション:</para>
103
<para> BおよびD行列を計算.</para>
116
共分散行列およびカルマン予測ゲインを計算する際に使用される
118
このパラメータは,共分散行列とカルマン予測ゲインが必要ない
119
場合には,不要です.NSMPL = 0の場合, K, Q, Ry, および S は計算されません.
128
行列のランクを推定する際に使用される許容誤差.
129
TOL > 0 の場合, TOL に指定した値が, 条件数の逆数の下界として
131
デフォルト: prod(size(matrix))*epsilon_machine ただし,
132
epsilon_machine は相対機械精度です.
139
<para>整数, 警告メッセージの出力に関するスイッチ</para>
144
<para>= 1: 警告メッセージを出力;</para>
150
<para>= 0: 警告メッセージを出力しない.</para>
162
<para>状態空間実現 SYS = (A,B,C,D) (syslin オブジェクト)を計算</para>
168
<para>カルマン予測ゲイン K (NSMPL > 0の場合)</para>
186
<para>状態出力相互共分散</para>
194
ランク決定, 最小二乗またはリカッチ方程式の解に関わる行列の条件数の逆数が含まれます.
203
システムの次数,行列AおよびCおよび関連する
204
結合ブロックハンケル行列のR因子を指定すると,
205
部分空間同定法 (MOESP または N4SID) により,
206
離散時間システムのシステム行列BおよびDとカルマンゲインを返します.
216
<programlisting role=""><![CDATA[
210
<para>[B,D,K] = findBDK(S,N,L,R,A,C,METH,JOB,NSMPL,TOL,PRINTW) は,
211
システム行列 B (JOB = 1の場合),
212
システム行列 BおよびD (JOB = 2の場合),
213
カルマン予測ゲイン K (NSMPL > 0の場合) を計算します.
216
<programlisting role=""><![CDATA[
217
217
x(k+1) = Ax(k) + Bu(k) + Ke(k), k >= 1,
218
218
y(k) = Cx(k) + Du(k) + e(k),
219
219
]]></programlisting>
221
ただし, x(k) および y(k) はそれぞれ長さ N および Lのベクトルです.
225
<para>[B,D,RCND] = findBDK(S,N,L,R,A,C,METH,JOB) は,
226
長さ 4 のベクトル RCND も返します.
227
このベクトルには,ランク決定に関連する条件数の逆数が含まれます.
231
<para>[B,D,K,Q,Ry,S,RCND] = findBDK(S,N,L,R,A,C,METH,JOB,NSMPL,TOL,PRINTW) は,
233
状態量, 出力, 出力状態 (相互)共分散行列 Q, Ry, および S
234
(カルマンゲイン計算で使用される)も返します.
235
RCNDベクトルには,ランク決定, 最小二乗またはリカッチ方程式の解に関連する
242
適当な引数 METH および JOBDで定義されます.
243
findR ではMETH = 1かつJOBD = 1,findBDKではMETH = 1を使用しますす.
244
FINDR では METH = 1, findBDK では METH = 2 を使用することができます.
247
出力引数の数は変換しますが,次のように入力引数に対応します,
249
<programlisting role=""><![CDATA[
221
ただし, x(k) および y(k) はそれぞれ長さ N および Lのベクトルです.
225
<para>[B,D,RCND] = findBDK(S,N,L,R,A,C,METH,JOB) は,
226
長さ 4 のベクトル RCND も返します.
227
このベクトルには,ランク決定に関連する条件数の逆数が含まれます.
231
<para>[B,D,K,Q,Ry,S,RCND] = findBDK(S,N,L,R,A,C,METH,JOB,NSMPL,TOL,PRINTW) は,
233
状態量, 出力, 出力状態 (相互)共分散行列 Q, Ry, および S
234
(カルマンゲイン計算で使用される)も返します.
235
RCNDベクトルには,ランク決定, 最小二乗またはリカッチ方程式の解に関連する
242
適当な引数 METH および JOBDで定義されます.
243
findR ではMETH = 1かつJOBD = 1,findBDKではMETH = 1を使用しますす.
244
FINDR では METH = 1, findBDK では METH = 2 を使用することができます.
247
出力引数の数は変換しますが,次のように入力引数に対応します,
249
<programlisting role=""><![CDATA[
250
250
B = findBDK(S,N,L,R,A,C,METH,1) または
251
251
[B,D] = findBDK(S,N,L,R,A,C,METH,2) または
252
252
[B,D,RCND] = findBDK(S,N,L,R,A,C,METH,2)
253
253
]]></programlisting>
257
<programlisting role="example"><![CDATA[
257
<programlisting role="example"><![CDATA[
258
258
//指定した線形システムからデータを生成
259
259
A = [ 0.5, 0.1,-0.1, 0.2;
260
260
0.1, 0, -0.1,-0.1;