¿Existen atajos para sistemas numéricamente aproximados de ecuaciones diferenciales ordinarias cuando son autónomos?

10

Los algoritmos existentes para resolver ODE manejan funciones , dondeyRn. Pero en muchos sistemas físicos, la ecuación diferencial es autónoma, entoncesdyreyret=F(y,t)yRnorte,yRn, con latfuera. Con este supuesto simplificador, ¿qué mejoras se pueden ver en los métodos numéricos existentes? Por ejemplo, sin=1, el problema se convierte ent=dyreyret=F(y)yRnortetnorte=1 y pasamos a una clase completamente diferente de algoritmos para integrar integrales unidimensionales. Paran>1, la mejora máxima posible es reducir la dimensión deyen 1, porque el caso dependiente del tiempo se puede simular agregandotay, cambiando el dominio deydeRnaRn+1.t=reyF(y)norte>1ytyyRnorteRnorte+1

postergar el trabajo real
fuente

Respuestas:

2

Yo diría que una mejora significativa es que en el alcance de los enfoques de paso de tiempo, donde propaga ynorteynorte+1=U(ynorte) utilizando un mapa de solución U , puede determinar el propagador (o al menos partes de ) una vez y luego reutilícela en cada paso del tiempo.

Por ejemplo, en el caso lineal, tendría ty=UNAy , donde UNA es una matriz. El operador de solución U(y)=Exp(UNAΔt)y consiste principalmente en una matriz exponencial. Para sistemas autónomos, esta costosa evaluación exponencial de matriz se requiere solo una vez para la propagación completa, en contraste con un sistema dependiente del tiempo, en el que debe realizar esta evaluación en cada paso de tiempo.

Para los sistemas no lineales no es tan fácil, pero dependiendo del algoritmo se pueden reutilizar ciertas evaluaciones costosas.

carlosvalderrama
fuente