14
14
<refentry xmlns="http://docbook.org/ns/docbook" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:svg="http://www.w3.org/2000/svg" xmlns:ns5="http://www.w3.org/1999/xhtml" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:db="http://docbook.org/ns/docbook" version="5.0-subset Scilab" xml:id="ode_discrete" xml:lang="ja">
16
<refname>ode_discrete</refname>
17
<refpurpose>常微分方程式ソルバ,
23
<synopsis>y=ode("discrete",y0,k0,kvect,f)</synopsis>
31
<para>実数のベクトルまたは行列 (初期条件).</para>
37
<para>実数スカラー (初期時間).</para>
43
<para>外部 (関数または文字列またはリスト).</para>
49
<para>整数 (初期時間).</para>
63
(最初の引数を<literal>"discrete"</literal>とする )
64
この構文により,<literal>ode</literal>は
65
初期状態量 <literal>y(k0)</literal> から再帰的に
66
<literal>y(k+1)=f(k,y(k))</literal>を計算し,
67
<literal>kvect</literal>の中にある<literal>k</literal>について
68
<literal>y(k)</literal>を返します.
69
<literal>kvect(1)</literal> は<literal>k0</literal>以上である
72
<para>他の引数およびその他のオプションは
73
<literal>ode</literal>と同じです.
74
<link linkend="ode">ode</link>のヘルプを参照ください.
79
<programlisting role="example"><![CDATA[
16
<refname>ode_discrete</refname>
17
<refpurpose>常微分方程式ソルバ,
23
<synopsis>y=ode("discrete",y0,k0,kvect,f)</synopsis>
31
<para>実数のベクトルまたは行列 (初期条件).</para>
37
<para>実数スカラー (初期時間).</para>
43
<para>外部 (関数または文字列またはリスト).</para>
49
<para>整数 (初期時間).</para>
63
(最初の引数を<literal>"discrete"</literal>とする )
64
この構文により,<literal>ode</literal>は
65
初期状態量 <literal>y(k0)</literal> から再帰的に
66
<literal>y(k+1)=f(k,y(k))</literal>を計算し,
67
<literal>kvect</literal>の中にある<literal>k</literal>について
68
<literal>y(k)</literal>を返します.
69
<literal>kvect(1)</literal> は<literal>k0</literal>以上である
72
<para>他の引数およびその他のオプションは
73
<literal>ode</literal>と同じです.
74
<link linkend="ode">ode</link>のヘルプを参照ください.
79
<programlisting role="example"><![CDATA[
80
80
y1=[1;2;3]; deff("yp=a_function(k,y)","yp=A*y+B*u(k)")
81
81
A=diag([0.2,0.5,0.9]); B=[1;1;1];u=1:10;n=5;
82
82
y=ode("discrete",y1,1,1:n,a_function);
85
85
// ここで, [y3,y5,y7,y9] で y を評価します
86
86
y=ode("discrete",y1,1,3:2:9,a_function)
87
87
]]></programlisting>
89
<refsection role="see also">
91
<simplelist type="inline">
93
<link linkend="ode">ode</link>
89
<refsection role="see also">
91
<simplelist type="inline">
93
<link linkend="ode">ode</link>