Metodo de Romberg

Instituto Tecnolgico de La Paz
Ingeniera Bioqumica
Programacin y Mtodos numricos
Alumnos:
-Higuera Brquez Carolina
-Orozco Flores Ramn
-Pia Espinoza Jennie
Docente:
Ramn de Jess Armenta Machado
Mtodo de Integracin de Romberg
La Paz B.C.S., a 14 de noviembre de 2013
Sea
el valor de la integral que aproxima a
una particin de subintervalos de longitud

trapecio. Entonces,
donde
, mediante
y usando la regla del
es el error de truncamiento que se comete al aplicar la regla.
El mtodo de extrapolacin de Richardson combina dos aproximaciones

de integracin numrica, para obtener un tercer valor ms exacto.
de
El algoritmo ms eficiente dentro de ste mtodo, se llama Integracin

Romberg, la
cual
es
una
frmula
recursiva.
Supongamos que tenemos dos aproximaciones:
Se puede demostrar que el error que se comete con la regla del trapecio
para n subintervalos est dado por las siguientes frmulas:
donde
es un promedio de la doble derivada entre ciertos valores que
pertenecen a cada uno de los subintervalos.
Ahora bien, si suponemos que el valor de
es constante, entonces:
Sustituyendo esto ltimo en nuestra primera igualdad, tenemos que:
De aqu podemos despejar
En el caso especial cuando

tenemos:
(que es el algoritmo de Romberg),
Esta frmula es solo una parte del algoritmo de Romberg. Para entender el
mtodo, es conveniente pensar que se trabaja en niveles de aproximacin. En
un primer nivel, es cuando aplicamos la regla del Trapecio, y para poder usar la
frmula anterior, debemos de duplicar cada vez el nmero de subintervalos:
as, podemos comenzar con un subintervalos, luego con dos, cuatro, ocho, etc.,
hasta donde se desee.
Posteriormente, pasamos al segundo nivel de aproximacin, que es donde se

usa la frmula anterior, tomando las parejas contiguas de aproximacin del
nivel anterior, y que corresponden cuando
Despus pasamos al nivel tres de aproximacin, pero aqu cambia la frmula

de Romberg, y as sucesivamente hasta el ltimo nivel, que se alcanza cuando
solo contamos con una pareja del nivel anterior.
Desde luego, el nmero de niveles de aproximacin que se alcanzan, depende
de las aproximaciones que se hicieron en el nivel 1. En general, si en el primer
nivel, iniciamos con n aproximaciones, entonces alcanzaremos a llegar hasta
el nivel de aproximacin n.
Esquema de integracin de Romberg
0
1
2
3
4
R (J, 0)
Regla del
trapecio
R (0,0)
R (1,0)
R (2,0)
R (3,0)
R (4,0)
R (J, 1)
Regla de
Simpson
R
R
R
R
(1,1)
(2,1)
(3,1)
(4,1)
R (J, 2)
Regla de
Boole
R (J, 3)
Tercera
mejora
R (J, 4)
Cuarta
mejora
R (2,2)
R (3,2)
R (4,2)
R (3,3)
R (4,3)
R (4,4)
Se define la sucesin {R (J, K): J K} j= 0 de formulas de cuadratura para

aproximar la integral de f (x) en [a, b] de la siguiente manera
R (J, 0) = T (J) para J 0, que es la regla recursiva del trapecio
R (J, 1) = S (J) para J 1, que es la regla recursiva de Simpson
R (J, 2) = B (J) para J 2, que es la regla recursiva de Boole
Las formulas {R (J, 0)} de partida se usan para generar las primeras mejores
{R (J, 1)} que, a su vez, se usan para generar las segundas mejoras {R (J, 2)};
R (J, 1) =
41 R ( J ,0 )R(J 1,0)
411
R (J, 2) =
4 R ( J ,1 )R (J 1, 1)
4 21
para J 1
para J 2
La regla general para construir recursivamente las mejores es
R (J, K) =
4 R ( J , K 1 )R( J 1, K 1)
4 K 1
para J K
Ejemplo:
Calcule el rea bajo la curva de una funcin f (x) = xx en [1, 2] para
n=10 aplicando el mtodo de Romberg.
Se calcular el primer nivel en donde de realizar la regla del trapecio
dependiendo de las particiones que indique cada seccin.
Rn
J1
h
2
R Trapecio
ba
n
0.1
0.05
2.008553692
ba
21
0.5
0.25
10.04276846
ba
2
2
0.25
0.125
5.021384231
ba
23
0.125
0.0625
2.510692115
ba
24
0.0625
0.03125
1.255346058
ba
25
0.03125
0.015625
0.627673028
ba
6
2
0.015625
0.0078125
0.313836514
ba
27
0.0078125
0.00390625
0.156918257
R Simpson
ba
8
2
0.00390625
0.001953125
0.078459128
ba
29
0.001953125
0.000976562
5
0.039229564
10
ba
210
0.000976562
5
0.000488281
25
0.019614782
R Trapecio:
x ex
1
1
2
2
1
[(1e ) (2e )] =
2.008553692
x ex
1
x ex
1
0.5
2
x ex
1
[(1e ) (2e )] =
x ex
1
1.255346058
[(1e1) (2e2)] =
x ex
1
0.0078125
2
[(1e1) (2e2)] =
0.156918257
0.25
2
2
1
[(1e ) (2e )] =
x ex
1
0.00390625
2
[(1e1) (2e2)] =
0.078459128
0.125
2
2
1
[(1e ) (2e )] =
2.510692115
5.021384231
0.015625
2
0.313836514
10.04276846
x ex
x ex
1
0.001953125
2
[(1e1) (2e2)] =
0.039229564
0.0625
2
[(1e1) (2e2)] =
x ex
1
0.019614782
0.0009765625
2
[(1e1) (2e2)] =
x ex
1
0.03125
2
[(1e1) (2e2)] =
0.627673028
Esquema de integracin de Romberg
R (J, 0)
Regla del
trapecio
R (J, 1)
Regla de
Simpson
2.0085536
92
10.042768
46
5.0213842
31
2.5106921
15
1.2553460
58
0.6276730
28
0.3138365
14
0.1569182
57
0.0784591
28
0.0392295
64
0.0196147
82
12.987507
38
3.2809226
56
1.6737944
4
1.2504735
24
0.6270594
1
0.3137598
75
0.1569086
78
0.0784679
3
0.0392294
14
0.0196147
63
1
2
3
4
5
6
7
8
9
10
R (J, 2) Regla de
Boole
R (J, 3)
Tercera
mejora
R (J, 4)
Cuarta
mejora
R Simpson:
K
R (J, K) =
A1 =
4 R ( J , K 1 )R( J 1, K 1)
4 K 1
41 ( 10.04276846 )2.008553692
411
= 12.98750738
A6 =
46 ( 0.313836514 )0.627673028
4 61
0.313759875
A2 =
2
4 ( 5.021384231 )10.04276846
421
A7 =
= 3.280922656
A3 =
43 ( 2.510692115 ) 5.021384231
4 31
A8 =
A9 =
= 1.250473524
A5 =
4 8 ( 0.078459128 )0.156918257
4 81
4 9 ( 0. 039229564 ) 0.078459128
4 91
A10 =
410 ( 0.019614782 )0.039229564

410 1
0.019614763
Para la tercera columna ser, la formula de Romberg quedar as:

A
42 A n+1 A n
4 21
A21 =
42 ( 3.347589488 )6.695178973
421
A22 =
4 2 () 2
4 1
4 2 () 2
4 1
A24 =
4 2 () 2
4 1
A25 =
2
4 () 2
4 1
0.039229414
= 0.62705941
0.07846793
4 4 (1.255346058 ) 2.510692115
4 4 1
45 ( 0.627673028 )1.255346058
45 1
0.156908678
= 1.67379444
A4 =
7
4 ( 0.156918257 ) 0.313836514
4 71
A26 =
2
4 () 2
4 1
4 2 () 2
4 1
A28 =
4 2 () 2
4 1
A29 =
4 2 () 2
4 1
A210 =
2
4 () 2
4 1

Metodo de Romberg

Cargado por

Copyright:

Formatos disponibles

Metodo de Romberg

Cargado por

Información del documento

Descripción original:

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Metodo de Romberg

Cargado por

Copyright:

Formatos disponibles

Instituto Tecnolgico de La Paz

Mtodo de Integracin de Romberg

La Paz B.C.S., a 14 de noviembre de 2013

el valor de la integral que aproxima a

una particin de subintervalos de longitud

es el error de truncamiento que se comete al aplicar la regla.

El mtodo de extrapolacin de Richardson combina dos aproximaciones

El algoritmo ms eficiente dentro de ste mtodo, se llama Integracin

Supongamos que tenemos dos aproximaciones:

Sustituyendo esto ltimo en nuestra primera igualdad, tenemos que:

De aqu podemos despejar

En el caso especial cuando

(que es el algoritmo de Romberg),

Posteriormente, pasamos al segundo nivel de aproximacin, que es donde se

Despus pasamos al nivel tres de aproximacin, pero aqu cambia la frmula

Se define la sucesin {R (J, K): J K} j= 0 de formulas de cuadratura para

La regla general para construir recursivamente las mejores es

Esquema de integracin de Romberg

410 ( 0.019614782 )0.039229564

Para la tercera columna ser, la formula de Romberg quedar as:

También podría gustarte