TPMN Grupo 7
TPMN Grupo 7
TPMN Grupo 7
Facultad Politcnica
INTEGRANTES:
Jonathan Dario Bernal Fernandez
Elber Andrs Rojas Diarte
Tema:
Unidad 6: Solucin Numrica de Ecuaciones Diferenciales
6.1 Mtodo de Heun con un solo corrector, con a = 1/2
6.2 Mtodo del Punto medio con a = 1
6.3 Mtodo de Raiston, con a = 2/3
6.4 Mtodo de Runge-Kutta clsico
Ao 2016
RUNGE-KUTTA
El mtodo de Runge Kutta es un mtodo numrico de resolucin de
ecuaciones diferenciales que surge como una mejora del mtodo de Euler. El
mtodo de Euler se puede considerar como un mtodo de Runge Kutta de
primer orden, el de Heun, es un mtodo de Runge Kutta de orden dos.
Los mtodos de Runge-Kutta logran la exactitud del procedimiento de una
serie de Taylor sin requerir el clculo de derivadas superiores. Existen
muchas variaciones, pero todas se pueden denotar en la forma generalizada
de la ecuacin
yi + 1 = yi + F(xi,yi,h)h
Donde F(xi,yi,h) se conoce como la funcin incremento la cual puede
interpretarse como una pendiente representativa en el intervalo. La funcin
incremento se escribe en forma general como:
F = a1k1 + a2k2 +.+ ankn
Donde las a son constantes y las k son:
k1 = f(xi,yi)
k2 = f(xi + p1h,yi + q11k1h)
k3 = f(xi + p2h,yi + q21k1h + q22k2h)
kn = f(xi + pnh,yi + q2n-1k1h + qn-1,2k2h + . + qn-1,n-1kn-1h)
Donde las p y q son constantes.
Como cada k es una evaluacin funcional, esta recurrencia hace que los
mtodos Runge-Kutta sean eficientes para la programacin. Existen varios
tipos de mtodos Runge-Kutta al emplear diferentes nmeros de trminos
en la funcin incremento como la especificada por n.
n = 1, es el mtodo de Euler. Una vez se elige n, se evalan las a, p y q al
igualar la funcin incremento a los trminos en la serie de expansin de
Taylor. La versin de segundo orden para la ecuacin en su forma
generalizada es:
Donde:
Los valores de a1, a2, p1 y q11 son evaluados al igualar el trmino de segundo
orden de la ecuacin dada con la expansin de la serie de Taylor.
Desarrollando tres ecuaciones para evaluar las cuatro incgnitas:
Como se tienen tres ecuaciones con cuatro incgnitas se tiene que suponer
el valor de una de ellas. Suponiendo que se especific un valor para a2, se
puede resolver de manera simultnea el sistema de ecuaciones obtenido:
Como se puede elegir un nmero infinito de valores para a2, hay un nmero
infinito de mtodos Runge-Kutta de segundo orden.
a2 = 1/2: Mtodo de Heun con un solo corrector, donde:
de
ecuaciones.
Una
versin
ampliamente
usada
es:
Ejemplo
dy
yx 2 1.2 y
dx
Donde:
y(0)=1
h = 0.25
Solucin
y i 1 y i k 2 h
k 1 f(xi, y i )
k 2 = f(x i
1
1
h , y i k 1 h)
2
2
Primera iteracin
k 1 f(x 0 , y 0 ) f (0 , 1)
k1 (1)(0) 2 1.2(1)
k1 1.2
1
k 2 f ( x0 h
2
y0
1
k1 h )
2
1
k 2 f (0 (0.25)
2
1
1 (1.2)(0.25))
2
,
k 2 f (0.125,0.85)
k 2 0.85(0.125) 2 1.2(0.85)
k 2 1.006718
y1 1 (1.006718)0.25
y1 0.748320
Segunda iteracin
x1 x 0 h
x1 0 0.25
x1 0.25
k 1 f(x 1 , y1 ) f (0.25 , 0.748320)
k 2 f (0.25
1
(0.25)
2
1
0.748320 (0.851432)(0.25))
2
,
k 2 f (0.375,0.641891)
k 2 0.641891(0.375) 2 1.2(0.641891)
k 2 0.680003
y 2 0.748320 ( 0.680003)0.25
y 2 0.578319
Tercera iteracin
x 2 x1 h
x 2 0.25 0.25
x 2 0.5
k 1 f(x 2 , y 2 ) f (0.5,0.578319)
k1 (0.578319)(0.5) 2 1.2(0.578319)
k1 0.549403
k 2 f ( x2
1
h
2
y2
1
k1 h )
2
1
k 2 f (0.5 (0.25)
2
1
0.578319 (0.549403)(0.25))
2
,
k 2 f (0.625,0.509643)
k 2 0.509643(0.625) 2 1.2(0.509643)
k 2 0.4125
y 3 0.578319 (0.4125)0.25
y 3 0.4752
Cuarta iteracin
x3 x 2 h
x3 0.5 0.25
x3 0.75
k 1 f(x 3 , y 3 ) f (0.75,0.4752)
k1 (0.4752)( 0.75) 2 1.2(0.4752)
k1 0.3029
1
k 2 f ( x3 h
2
y3
1
k1 h )
2
k 2 f (0.75
1
(0.25)
2
0.4752
,
k 2 f (0.875,0.4373)
k 2 0.4373(0.875) 2 1.2(0.4373)
k 2 0.1900
y 4 0.4752 ( 0.1900)0.25
1
( 0.3029)(0.25))
2
y 4 0.4277
x4 x3 h
x 4 0.75 0.25
x4 1
Vectores solucin
X
y
0
1
0.25
0.7483
0.5
0.5783
0.75
0.4752
1
0.4277
Se resuelve el mismo problema anterior pero esta vez mediante el uso del
mtodo Runge kutta de tercer grado, de valor inicial, en el intervalo de x=0
a x=1.
dy
yx 2 1.2 y
dx
Donde:
y(0)=1
h = 0.25
Solucin.
En el mtodo de Runge kutta de tercer orden se utilizan las siguientes
formulas:
1
y i 1 y i (k 1 4k 2 k 3 )h
6
k 1 f(xi, y i )
k 2 = f(x i
1
1
h , y i k 1 h)
2
2
k 3 f(x i h , y i k 1 h 2k 2 h)
Primera iteracin
k1 f(x 0 , y 0 ) f (0 , 1)
k1 (1)(0) 2 1.2(1)
k1 1.2
k 2 f ( x0
1
h
2
y0
,
1
k1 h )
2
1
k 2 f (0 (0.25)
2
1
1 (1.2)(0.25))
2
,
k 2 f (0.125,0.85)
k 2 0.85(0.125) 2 1.2(0.85)
k 2 1.0067
k 3 f(x o h , y o k 1 h 2k 2 h)
k 3 f (0 (0.25), (1) (1.2)(0.25) 2(1.0067)(0.25))
k 3 f (0.25,0.7966)
k 3 0.7966(0.25) 2 1.2(0.7966)
k 3 0.9062
1
y 1 y 0 (k 1 4k 2 k 3 )h
6
y1 0.7445
Segunda iteracin
x1 x 0 h
x1 0 0.25
x1 0.25
k 1 f(x 1 , y1 ) f (0.25 , 0.7445)
k 2 f ( x1
1
h
2
y1
1
k1 h )
2
k 2 f (0.25
1
(0.25)
2
0.7445
1
( 0.8469)( 0.25))
2
k 2 f (0.375,0.6386)
k 2 0.6386(0.375) 2 1.2(0.6386)
k 2 0.6765
k 3 f(x 1 h , y1 k 1 h 2k 2 h)
k 3 f (0.25 (0.25), (0.7445) (0.8469)(0.25) 2(0.6765)( 0.25))
k 3 f (0.5,0.6178)
k 3 0.6178(0.5) 2 1.2(0.6178)
k 3 0.5870
1
y 2 y1 (k 1 4k 2 k 3 )h
6
y 2 0.5720
Tercera iteracin
x 2 x1 h
x 2 0.25 0.25
x 2 0.5
k 1 f(x 2 , y 2 ) f (0.5,0.5720)
k1 (0.5720)(0.5) 2 1.2(0.5720)
k1 0.5434
k 2 f ( x2
1
h
2
y2
1
k1 h )
2
k 2 f ( 0 .5
1
(0.25)
2
0.5720
1
( 0.5434)(0.25))
2
k 2 f (0.625,0.5041)
k 2 0.5041(0.625) 2 1.2(0.5041)
k 2 0.4080
k 3 f(x 2 h , y 2 k 1 h 2k 2 h)
k 3 f (0.5 (0.25), (0.5720) (0.5434)(0.25) 2(0.4080)(0.25))
k 3 f (0.75,0.5038)
k 3 0.5038(0.75) 2 1.2(0.5038)
k 3 0.3212
1
y 3 y 2 (k 1 4k 2 k 3 )h
6
y 3 0.4679
Cuarta iteracin
x3 x 2 h
x3 0.5 0.25
x3 0.75
k 1 f(x 3 , y 3 ) f (0.75,0.4679)
k1 (0.4679)( 0.75) 2 1.2(0.4679)
k1 0.2986
k 2 f ( x3
1
h
2
y3
1
k1 h )
2
k 2 f (0.75
1
(0.25)
2
0.4679
,
k 2 f (0.875,0.4306)
k 2 0.4306(0.875) 2 1.2(0.4306)
k 2 0.1871
1
( 0.2983)( 0.25))
2
k 3 f(x 3 h , y 3 k 1 h 2k 2 h)
k 3 f (0.75 (0.25), (0.4679) (0.2983)(0.25) 2(0.1871)(0.25))
k 3 f (1,0.4489)
k 3 0.4489(1) 2 1.2(0.4489)
k 3 0.0898
1
y 4 y 3 (k 1 4k 2 k 3 )h
6
y 4 0.4206
x 4 x3 h
x 4 0.75 0.25
x4 1
Vectores solucin
X
y
0
1
0.25
0.7445
0.5
0.5720
0.75
0.4679
1
0.4206
x1 x0 h ; x2 x1 h ; etc
Del mismo modo, se calcula el valor aproximado de la solucin, y2, en el punto
x2 = x1 + h:
Con muy buena precisin, sin que, para ello, sea necesario que los h
k1 h f ( x1 , y1 )
El procedimiento
consta de los siguientes pasos:
x1 sean
= x0 muy
+ h, pequeos.
se calculan los
siguientes nmeros:
h
k
k
k 1 hhf (fx( x0 , y, y0 ) 1 )
2
1 2
h f ( x0 ,2yun
Parak 2 calcular
aproximado
de la solucin y1 en el
0 valor
h)
k
h
k
2
2
k
h
f
(
x
3
1
h
k ,
punto
k3 h f ( x0 , y0 22 )
2
2
y1
k h f ( x h, y k3 )
k 44 h f ( x0 1h, y0 k31)
1 1
KK
(k1 (2k
k 2 2k23 k2k 4
) 2k3
0 0
6 6 1
k4 )
y entonces se toma:
y 2 y1
K0
y1 y0 K 0
Y as, sucesivamente
Ejemplo
y(0) = 2 ,
con h = 0.2
Ejercicios Propuestos:
1. Considerar el problema de valor inicial:
3.
Dada
la
y(0) = 1
ecuacin
la
siguiente
tabla:
Bibliografa
Chapra, S y Canale, R (2003). Mtodos numricos para ingenieros. Editorial Macgraw Hill.
Cuarta edicin. Mxico.
http://matema.ujaen.es/jnavas/web_modelos/problemas/ptema7.pdf
http://aprendeenlinea.udea.edu.co/lms/moodle/mod/page/view.php?id=24538
http://www.frsn.utn.edu.ar/gie/an/mnedo/ejercicios%20propuestos.pdf
http://mathstools.com/section/main/Metodos_de_Runge_Kutta?lang=es