ANÁLISIS NUMÉRICO

Mag. Carlos Alberto Ardila Albarracín

BLOQUE 3. INTEGRACIÓN NUMÉRICA

3.7. MÉTODO DE ROMBERG
CURSO DE ANÁLISIS NUMÉRICO. BLOQUE 3

MÉTODO DE ROMBERG
La integración de Romberg es una técnica diseñada para obtener
integrales numéricas de funciones de manera eficiente

Se basa en la técnica de:

combinar dos estimaciones con la regla del trapecio cuyo error de estimación es O(h2)
para obtener una nueva estimación de O(h 4) [es decir, con menor error]

Posteriormente,

combinar dos estimaciones con la regla del trapecio con error de estimación O(h4)
para obtener una nueva estimación de O(h6) [es decir, con menor error]

2
CURSO DE ANÁLISIS NUMÉRICO. BLOQUE 3

MÉTODO DE ROMBERG
La integración de Romberg es una técnica diseñada para obtener
integrales numéricas de funciones de manera eficiente

Posteriormente,

combinar dos estimaciones con la regla del trapecio con error de estimación O(h6)
para obtener una nueva estimación de O(h8) [es decir, con menor error]

Hasta (aunque se podría seguir)…

combinar dos estimaciones con la regla del trapecio con error de estimación O(h8)
para obtener una nueva estimación de O(h10) [es decir, con menor error]

3
CURSO DE ANÁLISIS NUMÉRICO. BLOQUE 3

MÉTODO DE ROMBERG
Se comienza con:
combinar dos estimaciones con la
regla del trapecio cuyo error de
estimación es O(h2)

Para que veamos la efectividad del método,

se usarán estimaciones con regla del Trapecio
con 1, 2, 4, 8, y 16 segmentos
(es decir k = 5 aproximaciones)
para estimar:

 /3

2
sen ( x)dx
0

4
CURSO DE ANÁLISIS NUMÉRICO. BLOQUE 3

MÉTODO DE ROMBERG
Ahora se continúa de la siguiente manera:
combinar dos estimaciones con la regla del trapecio cuyo error de estimación es O(h2)
para obtener una nueva estimación de O(h 4) [es decir, con menor error]

Se calcula así:
I(1,2) =(4/3)I(2,1) - (1/3)I(1,1)

Se calcula así:
I(2,2) =(4/3)I(3,1) - (1/3)I(2,1)

Se calcula así:
I(3,2) =(4/3)I(4,1) - (1/3)I(3,1)

Se calcula así:
I(4,2) =(4/3)I(5,1) - (1/3)I(4,1)

Noten que vamos llenando parte de una matriz (K × K):

Columna 1: Aproximaciones del método del Trapecio
Columna 2: Aproximaciones con las fórmulas Romberg O(h4)

5
CURSO DE ANÁLISIS NUMÉRICO. BLOQUE 3

MÉTODO DE ROMBERG
Posteriormente,
combinar dos estimaciones con la regla del trapecio con error de estimación O(h4)
para obtener una nueva estimación de O(h6) [es decir, con menor error]

Se calcula así:
I(1,3) =(16/15)I(2,2) - (1/15)I(1,2)

Se calcula así:
I(2,3) =(16/15)I(3,2) - (1/15)I(2,2)

Se calcula así:
I(2,3) =(16/15)I(4,2) - (1/15)I(3,2)

Noten que seguimos llenando parte de una matriz (K × K):

Columna 1: Aproximaciones del método del Trapecio
Columna 2: Aproximaciones con las fórmulas Romberg O(h4)
Columna 3: Aproximaciones con las fórmulas Romberg O(h6)
6
CURSO DE ANÁLISIS NUMÉRICO. BLOQUE 3

MÉTODO DE ROMBERG
Posteriormente,
combinar dos estimaciones con la regla del trapecio con error de estimación O(h6)
para obtener una nueva estimación de O(h8) [es decir, con menor error]

Se calcula así:
I(1,4) =(64/63)I(2,3) - (1/63)I(1,3)

Se calcula así:
I(2,4) =(64/63)I(3,3) - (1/63)I(2,3)

Noten que seguimos llenando parte de una matriz (K × K):

Columna 4: Aproximaciones con las fórmulas Romberg O(h8)

7
CURSO DE ANÁLISIS NUMÉRICO. BLOQUE 3

MÉTODO DE ROMBERG
Hasta (aunque se podría seguir)…
combinar dos estimaciones con la regla del trapecio con error de estimación O(h8)
para obtener una nueva estimación de O(h10) [es decir, con menor error]

Se calcula así:
I(1,5) =(256/255)I(2,4) - (1/255)I(1,4)

Terminamos de llenar parte de una matriz (K × K):

Columna 5: Aproximación con fórmula Romberg O(h10)

8
CURSO DE ANÁLISIS NUMÉRICO. BLOQUE 3

MÉTODO DE ROMBERG
Calculamos el error…

I (1, k )  I (1, k  1)
Ea  *100%
I (1, k )

Para el ejemplo que venimos tratando:

I (1,5)  I (1, 4) 0.30709242  0.30709241

Ea  *100%  *100%  6.1493*10 6%
I (1,5) 0.30709242

Logrado usando solamente un total de 31 segmentos

a partir de las 5 aproximaciones con el Método del Trapecio

---------- FIN DEL DOCUMENTO

9