ode_discrete

パラメータ

y0

実数のベクトルまたは行列 (初期条件).

t0

実数スカラー (初期時間).

f

外部 (関数または文字列またはリスト).

k0

整数 (初期時間).

kvect

整数ベクトル.

説明

(最初の引数を"discrete"とする ) この構文により,odeは 初期状態量 y(k0) から再帰的に y(k+1)=f(k,y(k))を計算し, kvectの中にあるkについて y(k)を返します. kvect(1) はk0以上である 必要があります.

他の引数およびその他のオプションは odeと同じです. odeのヘルプを参照ください.

例

y1=[1;2;3]; deff("yp=a_function(k,y)","yp=A*y+B*u(k)")
A=diag([0.2,0.5,0.9]); B=[1;1;1];u=1:10;n=5;
y=ode("discrete",y1,1,1:n,a_function);
y(:,2)-(A*y1+B*u(1))

// ここで, [y3,y5,y7,y9] で y を評価します
y=ode("discrete",y1,1,3:2:9,a_function)

参照

• ode — 常微分方程式ソルバ