Campo matemático de ecuaciones diferenciales ordinarias numéricas
En el campo matemático de ecuaciones diferenciales ordinarias numéricas , un integrador geométrico es un método numérico que preserva las propiedades geométricas del flujo exacto de una ecuación diferencial.
Ejemplo de péndulo
Podemos motivar el estudio de los integradores geométricos considerando el movimiento de un péndulo .
Supongamos que tenemos un péndulo cuyo cuerpo tiene masa y cuya varilla no tiene masa y su longitud es . Supongamos que la aceleración debida a la gravedad es . Denotemos por el desplazamiento angular de la varilla desde la vertical y por el momento del péndulo. El hamiltoniano del sistema, la suma de sus energías cinética y potencial , es
Es natural tomar el espacio de configuración de todos como el círculo unitario , de modo que se encuentra en el cilindro . Sin embargo, tomaremos , simplemente porque el espacio -es más fácil de representar gráficamente. Definamos
y . Experimentemos usando algunos métodos numéricos simples para integrar este sistema. Como es habitual, seleccionamos un tamaño de paso constante, , y para un entero no negativo arbitrario escribimos . Usamos los siguientes métodos.
(Nótese que el método de Euler simpléctico trata q mediante el método de Euler explícito e implícito).
La observación de que es constante a lo largo de las curvas solución de las ecuaciones de Hamilton nos permite describir las trayectorias exactas del sistema: son las curvas de nivel de . Trazamos, en , las trayectorias exactas y las soluciones numéricas del sistema. Para los métodos de Euler explícito e implícito tomamos , y z 0 = (0,5, 0) y (1,5, 0) respectivamente; para los otros dos métodos tomamos , y z 0 = (0, 0,7), (0, 1,4) y (0, 2,1).
El método de Euler explícito (o implícito) se desarrolla en espiral desde el origen (o hacia el origen). Los otros dos métodos muestran el comportamiento cualitativo correcto, ya que la regla del punto medio implícito coincide con la solución exacta en mayor medida que el método de Euler simpléctico.
Recuerde que el flujo exacto de un sistema hamiltoniano con un grado de libertad preserva el área, en el sentido de que
Para todos .
Esta fórmula se puede verificar fácilmente a mano. En nuestro ejemplo del péndulo, vemos que el flujo numérico del método explícito de Euler no preserva el área; es decir,
Se puede realizar un cálculo similar para el método de Euler implícito, donde el determinante es
Sin embargo, el método de Euler simpléctico preserva el área:
Por lo tanto , la regla del punto medio implícita tiene propiedades geométricas similares.
En resumen: el ejemplo del péndulo muestra que, además de que los métodos de Euler explícitos e implícitos no son buenas opciones para resolver el problema, el método de Euler simpléctico y la regla del punto medio implícita concuerdan bien con el flujo exacto del sistema, y la regla del punto medio concuerda más estrechamente. Además, estos dos últimos métodos preservan el área, al igual que el flujo exacto; son dos ejemplos de integradores geométricos (de hecho, simplécticos ).
Método de marco móvil
El método del marco móvil se puede utilizar para construir métodos numéricos que preserven las simetrías de Lie de la EDO. Los métodos existentes, como el de Runge-Kutta, se pueden modificar utilizando el método del marco móvil para producir versiones invariantes. [1]
^ Pilwon Kim (2006), "Invariantización de esquemas numéricos utilizando marcos móviles"
Lectura adicional
Hairer, Ernst; Lubich, Christian; Wanner, Gerhard (2002). Integración numérica geométrica: algoritmos que preservan la estructura para ecuaciones diferenciales ordinarias . Springer-Verlag. ISBN 3-540-43003-2.
Leimkuhler, Ben; Reich, Sebastian (2005). Simulación de dinámica hamiltoniana . Cambridge University Press. ISBN 0-521-77290-7.
Budd, CJ; Piggott, MD (2003). "Integración geométrica y sus aplicaciones". Manual de análisis numérico. Vol. 11. Elsevier. págs. 35–139. doi :10.1016/S1570-8659(02)11002-7. ISBN 9780444512475.
Kim, Pilwon (2007). "Invariantización de esquemas numéricos utilizando marcos móviles". BIT Numerical Mathematics. Vol. 47, núm. 3. Springer. págs. 525–546. doi :10.1007/s10543-007-0138-8.