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