Ecuaciones diferenciales ordinarias

Ejemplo:

dy 3t
=t e −2 y
dt
y ( 0 )=0
La solución analítica

y=1/5*t.*exp(3*t)-1/25*exp(3*t)+1/25*exp(-2*t);

Método de Euler

y n+ 1= y n +hf (t n , y n )
Donde

f ( t n , y n )=t n e3 t −2 y n
n

La ecuación anterior se podía resolver con diferentes pasos de tiempo

Por ejemplo h=0.1 (la línea azul s la solución analítica y la línea roja es la solución numérica)

3.5

2.5

1.5

0.5

0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Refinando el paso de tiempo (h=0.01) se acerca más a la solución analítica (la línea roja que esta
pegada con la línea azul).

3.5

2.5

1.5

0.5

0
0 0.2 0.4 0.6 0.8 1 1.2 1.4

Método de Taylor

Es posible generar un serie de Taylor de ”y”

y n+ 1= y ( t n +h )

y n= y ( t n )

Al momento de obtener la serie de Taylor

d y ( t n ) h2 d2 y ( t n ) h3 d 3 y ( t n ) h 4 d 4 y ( t n )
y n+ 1= y ( t n ) + h + + + +…
dt 2 ! d t2 3! d t 3 4 ! d t4

Se tiene que evaluar las derivadas.

La primera derivada ya la tenemos

d y (t n)
=f (t n , y n )
dt

La segunda derivada
 ( )
2
d y (t n ) d d y (tn) d
2
= = f (t n , y n )
dt dt d t dt

El termino nos quedaría como

d ∂ f ( t n , y n ) ∂ f ( t n , y n ) dy
f (t n , yn )= +
dt ∂t ∂y dt

Entonces

d ∂ f ( t n , y n ) ∂ f ( tn , y n )
f (t n , yn)= + f (t n , y n )
dt ∂t ∂y

Para nuestro ejemplo utilizaríamos el método de Taylor de segundo orden.

2
d y ( t n ) h2 d y ( t n )
y n+ 1= y ( t n ) + h +
dt 2 ! d t2

Que nos viene quedando como:

3tn
f ( t n , y n )=t n e −2 y n

Para el siguiente término

d ∂ f ( t n , y n ) ∂ f ( tn , y n )
f (t n , yn)= + f (t n , y n )
dt ∂t ∂y
∂ f (t n , yn ) 3tn 3tn 3 tn
=3 t n e + e =e ( 3 t n +1 )
∂t
∂ f (t n , yn )
=−2
∂y
Por lo que

d
f ( t n , y n ) =e ( 3 t n +1 ) −2 [ t n e −2 y n ]=e ( t n +1 ) +4 y n
3tn 3t n 3t n

dt
Entonces

h2 3 t
y n+ 1= y ( t n ) + h [ t n e −2 y n ] +
3t n

2!
[ e ( t n+ 1 ) +4 y n ]
n
 3.5

3 solucion analitica
solucion numerica Euler
2.5 solucion Taylor 2do orden

1.5

0.5

0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1

Segundo ejemplo

()
2
y y
y'= −
t t
y ( 1 )=1
Evaluar en el rango 1 ≤t ≤1.2

Solución analítica

La solución propuesta es.

n
A y =B
Donde

dy A ' (t ) B '(t ) 1−n

+ y= y
dt nA (t ) nA(t )
Entonces

1−n=2 n=−1
A ' (t ) 1
= A ( t )=t
A (t ) t
 B ' (t) 1
= B ( t ) =C+lnt
A(t ) t 2

Entonces la solución analítica queda como

t y −1=C+lnt
Utilizando la condición inicial

y ( 1 )=1
1=C
Entonces

t
y=
1+lnt

Con el método de Euler

y n+ 1= y n +hf (t n , y n )
Que en este ejemplo

( )
2
yn yn
f ( t n , y n )= −
tn tn

[ ( )]
2
y y
y n+ 1= y n +h n − n
tn tn

Con el método de Taylor de Segundo orden

2
d y ( t n ) h2 d y ( t n )
y n+ 1= y ( t n ) + h +
dt 2 ! d t2

Que nos viene quedando como:

( )
2
y y
f ( t n , y n )= n − n
tn tn

Para el siguiente término

 d ∂ f ( t n , y n ) ∂ f ( tn , y n )
f (t n , yn)= + f (t n , y n )
dt ∂t ∂y
∂ f (t n , yn) − yn yn
2
= 2
+2 3
∂t tn tn

∂ f (t n , yn ) 1 yn
= −2 2
∂y tn tn

Por lo que

( )( )
d − yn y n2 1 y y n y n2
f ( t n , y n ) = 2 +2 3 + −2 n2 −
dt tn tn tn tn t n t n2

Entonces desarrollando el producto

d − yn y n2 y n y n2 y n2 y n3
f ( t , y ) = 2 +2 3 + 2 − 3 −2 3 +2 4
dt n n tn tn t n tn tn tn

simplificando

d − y n2 y n3
f ( t n , y n ) = 3 +2 4
dt tn tn

Entonces

[ ( )] [ ]
2 2 2 3
y y h − yn y
y n+ 1= y ( t n ) + h n − n + +2 n4
tn tn 2 ! tn 3
tn

Resultados (h=0.05)
 1.018

1.016 solucion analitica

solucion numerica Euler
1.014 solucion Taylor 2do orden

1.012

1.01

1.008

1.006

1.004

1.002

1
1 1.05 1.1 1.15 1.2 1.25

Métodos de Runge Kutta

dy
=f (t , y)
dt
Desarrollando por serie de Taylor

d y ( t n ) h2 d2 y ( t n ) h3 d 3 y ( t n ) h 4 d 4 y ( t n )
y n+ 1= y ( t n ) + h + + + +…
dt 2 ! d t2 3! d t 3 4 ! d t4

Método de Runge Kutta de segundo orden.

2
d y ( t n ) h2 d y ( t n )
y n+ 1= y ( t n ) + h +
dt 2 ! d t2

En el método de Runge Kutta se tendría que:

h
y n+ 1= y ( t n ) +
2
[ f ( t n , y n ) + f (t n +ah , y n +kh)]

Desarrollando por series de Taylor la función f ¿ )

 [
2 2 2
∂ f (t n, yn) ∂ f (t n, yn) 1 ∂ f ( tn , y n ) ∂ f (t n, yn) ∂ f ( tn , y n
f ( t n +ah , y n +kh )=f ( t n , y n ) + ah+ kh+ ( ah )2+2 ( ah )( kh ) +
∂t ∂y 2! ∂t 2
∂t∂ y ∂ y2

Sustituyendo en la ecuación (truncando hasta los términos de la primer derivada de la función f)

y n+ 1= y ( t n ) +
h
2 [
2 f (t n, yn )+
∂ f (t n , yn)
∂t
ah+
∂ f (t n , yn)
∂y
kh ]
Simplificando se tiene que

[
h ∂ f ( t n , y n) ∂ f (t n , yn )
]
2
y n+ 1= y ( t n ) + h f ( t n , y n ) + a+ k
2 ∂t ∂y

Y comparándolo con la ecuación de Taylor de segundo orden

[ ]
2
d y ( tn ) h d y ( tn )
2
h ∂ f ( tn , yn ) ∂ f ( t n , y n )
2
y n+ 1= y ( t n ) + h + = y n+1 = y ( n)
t +h f ( n n) 2 !
t , y + + f (t n , y n )
dt 2 ! d t2 ∂t ∂y
Para que fuesen iguales

a=1
k =f (t n , y n)

Entonces el método de Runge-Kutta de segundo orden queda como:

h
y n+ 1= y ( t n ) +
2
[ f ( t n , y n ) + f (t n +h , y n +kh)]
Donde

k =f (t n , y n)

Y Regularmente se presenta como

h
y n+ 1= y ( t n ) + [ k +k ]
2 1 2
Donde

k 1=f (t n , y n)

k 2=f (t n +h , y n +k 1 h)

Para el ejemplo anterior nos da los siguientes resultados

1.018

1.016 solucion analitica

solucion numerica Euler
1.014 solucion Taylor 2do orden
solucion Runge Kutta 2do orden
1.012

1.01

1.008

1.006

1.004

1.002

1
1 1.05 1.1 1.15 1.2 1.25

Regularmente se utiliza método de Runge Kutta de cuarto orden.

h
y n+ 1= y ( t n ) + [ k +2 k 2 +2 k 3 +k 4 ]
6 1
Donde

k 1=f (t n , y n)

(h
k 2=f t n + , y n +
2
k1 h
2 )
 (h
k 3=f t n+ , y n+
2
k2 h
2 )
k 4=f ( t n+ h , y n + k 3 h )

1.015

solucion analitica
solucion Runge Kutta 4to. orden

1.01

1.005

1
1 1.05 1.1 1.15 1.2 1.25

Método de Adams Bashforth

Si se tiene los valores de y en los pasos anteriores, estos se pueden utilizar

y n+ 1= y ( t n ) + h [ α f ( t n , y n )+ β f (t n−1 , y n−1) ]

Obtener la serie de Taylor de f (t n−1 , y n−1)

f ( t n−1 , y n−1 ) =f (t n −h , y n−Δ y )

Δ y = y n− y n−1

De tal forma que

∂ f (t n, yn) ∂ f (t n , yn)
f ( t n−h , y n −Δ y ) =f ( t n , y n ) −h −Δ y +…
∂t ∂y
Esta otra expresión quedaría como
 2
d y (t n) h2 d y (t n)
y n−1= y ( t n−h )= y ( t n )−h + −…
dt 2 ! d t2

Por lo que
2 2
d y ( t n ) h2 d y ( t n ) h d y (t n)
2
Δ y = y n− y n−1=h
dt
−
2! d t2
+ …=h f ( t n , y n ) −
2! dt2
+…

sustitituyendo

∂ f (t n, yn) ∂ f ( tn , yn )
f ( t n−h , y n −Δ y ) =f ( t n , y n ) −h −h f ( t n , y n ) +…
∂t ∂y

Entonces

[ [
y n+ 1= y ( t n ) + h α f ( t n , y n ) + β f ( t n , y n ) −h
∂ f ( tn , y n )
∂t
−h f ( t n , y n )
∂ f (t n , yn)
∂y ]]
Comparando con la expresión de Taylor de segundo orden

y n+ 1= y ( t n ) + h
d y ( t n ) h2 d2 y ( t n )
dt
+
2 ! d t2
= y n+1= y ( t n ) +h f ( t n , y n ) +
2! ∂t [
h2 ∂ f ( t n , y n ) ∂ f ( t n , y n )
+
∂y ]
f (t n , y n )

Para que fuesen iguales

α + β=1
−1
β=
2
3
α=
2

y n+ 1= y ( t n ) + h
[ 3 1
f (t , y )− f (t , y )
2 n n 2 n−1 n−1 ]