MÉTODO DE INTEGRACIÓN DE ROMBERG
Sea el valor de la integral que aproxima a , mediante una partición
de subintervalos de longitud , usando la regla del trapecio. Entonces,
donde es el error de truncamiento que se comete al aplicar la regla trapecial.
El método de extrapolación de Richardson combina dos aproximaciones de
integración numérica, para obtener un tercer valor más exacto.
El algoritmo más eficiente dentro de este método, se llama Integración de Romberg,
la cual es una fórmula recursiva.
Supongamos que tenemos dos aproximaciones: e , con subintervalos h1 y h2
respectivamente.
Se ha visto que el error que se comete con la regla del trapecio para n subintervalos está
dado por las siguientes fórmulas:
donde ƒ′′(s) 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 (que es el algoritmo de Romberg), tenemos:
Esta fórmula es solo una parte del algoritmo de Romberg. Para entender el método, es
conveniente pensar que se trabaja en niveles de aproximación. En un primer nivel (que
llamamos 0), es cuando aplicamos la regla del Trapecio, y para poder usar la fórmula
anterior, debemos de duplicar cada vez el número de subintervalos: así, podemos comenzar
con un subintervalo, luego con dos, cuatro, ocho, etc., hasta donde se desee.
Posteriormente, pasamos al segundo nivel de aproximación (el 1), que es donde se usa la
fórmula anterior, tomando las parejas contiguas de aproximación del nivel anterior, y que
corresponden cuando .
Después pasamos al tercer nivel de aproximación (el 2), pero aquí cambia la fórmula de
Romberg, y así sucesivamente hasta el último nivel, que se alcanza cuando solo contamos
con una pareja del nivel anterior.
Desde luego, el número de niveles de aproximación que se alcanzan, depende de las
aproximaciones que se hicieron en el nivel 0. En general, si en el primer nivel, iniciamos con
n aproximaciones, entonces alcanzaremos a llegar hasta el nivel de aproximación n.
Hacemos un diagrama para explicar un poco más lo anterior.
Ejemplo 1.
Usar el algoritmo de Romberg, para aproximar la integral
usando segmentos de longitud .
Solución.
Primero calculamos las integrales del nivel 1, usando la regla del trapecio para las
longitudes de segmentos indicadas:
Con estos datos, tenemos:
Ahora pasamos al segundo nivel de aproximación donde usaremos la fórmula que se dedujo
anteriormente:
donde es la integral menos exacta (la que usa menos subintervalos) e es la
más exacta (la que usa el doble de subintervalos).
En un diagrama vemos lo siguiente:
Para avanzar al siguiente nivel, debemos conocer la fórmula correspondiente. De forma
similar a la deducción de la fórmula,
se puede ver que la fórmula para el siguiente nivel de aproximación (nivel 3) queda como
sigue:
dónde:
es la integral más exacta
es la integral menos
exacta
En el siguiente nivel (nivel 4) se tiene la fórmula
En el ejemplo anterior, obtenemos la aproximación en el nivel 3 como sigue:
Así, podemos concluir que el valor de la aproximación, obtenido con el método de Romberg
en el ejemplo 1, es:
Ejemplo 2.
Usar el algoritmo de Romberg para aproximar la integral:
Agregando a la tabla anterior donde .
Solución.
Calculamos con la regla del trapecio:
Tenemos entonces la siguiente tabla:
De donde concluimos que la aproximación buscada es:
Ejemplo 3.
Aproximar la siguiente integral:
usando el método de Romberg con segmentos de longitud
, , ,
Solución.
Igual que arriba, primero usamos la regla del trapecio (con los valores de h indicados) para
llenar el nivel 1. Tenemos entonces que:
A continuación, usamos las fórmulas de Romberg para cada nivel y obtenemos la siguiente
tabla:
De donde concluimos que la aproximación buscada es:
Podemos escribir una fórmula general para calcular las aproximaciones en cada uno de los
niveles como sigue:
ALGORITMO DE INTEGRACIÓN DE ROMBERG
El método, tomando en cada reiteración una división a la mitad del intervalo, se define de
forma recursiva así:
La extrapolación a orden cero R(n,0) es equivalente a la Regla del trapecio con n+2
puntos. La de orden uno R(n,1) es equivalente a la Regla de Simpson con n+2 puntos.
El proceso se suele realizar hasta que se cumpla alguna de las condiciones, según sea la
necesidad:
Error absoluto: ∈𝑎= |𝑅 (𝑛, 𝑚) − 𝑅 (𝑛, 𝑚 − 1) | ≤ c
𝑅 (𝑛, 𝑚) −𝑅 (𝑛, 𝑚−1)
Error relativo: ∈𝑟= | |≤c
𝑅 (𝑛, 𝑚)
Error relativo porcentual: ∈𝑟𝑝 = |𝑅 (𝑛, 𝑚) −𝑅 (𝑛, 𝑚−1) | * 100 ≤ c
𝑅 (𝑛, 𝑚)
donde c es una cota suficiente de error, definida en forma absoluta, relativa o porcentual
según corresponda.
Ejemplo 1.
Aplicar el algoritmo de integración de Romberg a la integral:
tomando
Solución.
En este caso no sabemos exactamente cuantas aproximaciones debemos
hacer con la regla del trapecio. Así que para comenzar hacemos los cálculos
correspondientes a uno, dos, cuatro y ocho subintervalos:
Con estos datos, podemos hacer los cálculos desde el nivel 1 hasta el nivel 3.
Tenemos la siguiente tabla:
Haciendo los cálculos de los errores, nos damos cuenta que efectivamente la
aproximación se obtiene hasta el nivel 3, donde ∈𝑟𝑝= 0.008%
Por lo tanto, concluimos que la aproximación buscada es:
