10: Integración Numérica de ODEs
- Última actualización
- 30 oct 2022
- Guardar como PDF
( \newcommand{\kernel}{\mathrm{null}\,}\)
Este artículo describe los métodos numéricos para resolver el problema del valor inicial, que es un tipo estándar de problema que aparece en muchos campos de la física. Supongamos que tenemos un sistema cuyo estado en el momentot es descrito por un vector→y(t), que obedece a la ecuación diferencial ordinaria de primer orden (ODE) para la forma:
d→ydt=→F(→y(t),t).
Aquí,→F se da alguna función vectorizada, cuyas entradas son (i) el estado instantáneo→y(t) y (ii) el tiempo actualt. Entonces, dado un tiempo inicialt0 y un estado inicial→y(t0), el objetivo es encontrar→y(t) para tiempos posteriores.
Conceptualmente, el problema del valor inicial es distinto del problema de resolver una ODE discutido en el artículo sobre ecuaciones de diferencia finita. Ahí, se nos dio un par de límites con ciertas condiciones límite, y el objetivo era encontrar la solución entre los dos límites. En este caso, se nos da el estado en un momento inicialt0, y nuestro objetivo es encontrar→y(t) para algún conjunto de tiempos futurost>t0. Esto a veces se denomina “integración” de la ODE, porque la solución tiene la forma
→y(t)=→y(t0)+∫tt0dt′→F(→y(t′),t′).
Sin embargo, a diferencia de la integración numérica ordinaria (es decir, la computación de una integral definida), el valor del integrando no se conoce de antemano, debido a la→F dependencia de lo desconocido→y(t).
