Campus Livre
ode23t
Solver de regra do trapézio para equações moderadamente rígidas e sistemas algébrico-diferenciais.
Categoria: Mathematics > Numerical Integration and Differential Equations
Sintaxe
[t, y] = ode23t(odefun, tspan, y0)Integra y' = odefun(t, y).
[t, y] = ode23t(odefun, tspan, y0, options)Aplica também opções de odeset, incluindo uma matriz 'Mass'.
Descrição
A regra do trapézio implícita, com um interpolante livre, não amortece oscilações, o que faz de ode23t o solver rígido preferido quando a solução é oscilatória e se quer preservar sua amplitude. Aceita matriz de massa singular, então também resolve equações algébrico-diferenciais de índice 1.
Para amortecer fortemente transientes rígidos, ode15s ou ode23tb são melhores.
Argumentos de entrada
odefun- Handle
dydt = odefun(t, y). tspan[t0 tf]ou um vetor de instantes de saída.y0- Vetor de valores iniciais.
Argumentos de saída
t, y- Instantes e solução, com o mesmo formato de
ode45.
Exemplo
Decaimento rumo a um cosseno
f = @(t, y) -50*(y - cos(t));
[t, y] = ode23t(f, [0 1], 0);
fprintf('%.3f\n', y(end))Saída
0.557