ode_discrete

Parâmetros

y0

vetor de reais ou matriz (condições iniciais).

t0

escalar real (tempo inicial).

f

função externa i.e. função, string ou lista.

k0

inteiro (tempo inicial).

kvect

vetor de inteiros.

Descrição

Com esta sintaxe (primeiro argumento igual a "discrete") ode, computa recursivamente y(k+1)=f(k,y(k)) de um estado inicial y(k0) e retorna y(k) para k em kvect. kvect(1) deve ser maior do que ou igual a k0.

Outros argumentos e opções são os mesmos que para ode, veja a página de ajuda de ode.

Exemplos

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))

// Agora, y é avaliado em [y3,y5,y7,y9]
y = ode("discrete", y1, 1, 3:2:9, a_function)

Ver Também

• ode — Solucionador de equações diferenciais ordinárias