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 — 常微分方程式ソルバ