`decic`

Encontra condições iniciais consistentes para uma equação diferencial totalmente implícita.

Categoria: Mathematics > Numerical Integration and Differential Equations

Sintaxe

[y0mod, yp0mod] = decic(odefun, t0, y0, fixed_y0, yp0, fixed_yp0)

Ajusta as componentes livres de y0 e yp0 de modo que odefun(t0, y0mod, yp0mod) = 0.

[y0mod, yp0mod] = decic(___, options)

Aplica também tolerâncias de odeset.

Descrição

Um solver implícito como ode15i precisa de valores iniciais que já satisfaçam a equação. decic os calcula mantendo as componentes marcadas como fixas e resolvendo para as demais. O número de componentes livres deve coincidir com o número de equações.

Se a equação não for consistente para nenhuma escolha dos valores livres, a função informa que não convergiu.

Argumentos de entrada

odefun: Handle res = odefun(t, y, yp) do sistema implícito f(t, y, y') = 0.
t0: Instante inicial.
y0, yp0: Palpites iniciais para a solução e sua derivada.
fixed_y0, fixed_yp0: Vetores de 1 e 0 que marcam quais componentes devem ser mantidas (1) ou podem ser alteradas (0).

Argumentos de saída

y0mod, yp0mod: Valores iniciais consistentes, próximos dos palpites, prontos para ode15i.

Exemplo

y' + y = 0 com y(0) = 1

f = @(t, y, yp) yp + y;
[y0, yp0] = decic(f, 0, 1, 1, 0, 0)

Saída

y0 = 1
yp0 = -1

Veja também

ode15i odeset