Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Romberg

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 13

MÉTODO DE INTEGRACION DE ROMBERG

INTEGRACIÓN DE ROMBERG
Sea el valor de la integral que aproxima a
, mediante una partición de subintervalos de

longitud y usando la regla del trapecio. Entonces,

donde es el error de truncamiento que se comete al


aplicar la regla.
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 éste método, se llama
Integración de Romberg , la cual es una fórmula recursiva.

Supongamos que tenemos dos aproximaciomnes : e


Se puede demostrar que el error que se comete con la regla
del trapecio para n subintervalos está dado por las siguientes
fórmulas:

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 (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, 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,
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 nivel tres de aproximación, 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 1. 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ómula


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:

donde:
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 concluímos 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 concluímos 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
Los coeficientes en cada una de las fórmulas en el
método de Romberg, deben sumar 1.
Así se tiene la siguiente fórmula recursiva:

donde:
es la integral más exacta
es la integral menos exacta

y el indice k indica el nivel de integración o de


aproximación. Por ejemplo, digamos que ,
entonces tenemos:
que es nuestra fórmula del nivel 2 de aproximación.
Como todo proceso iterativo, éste se detiene cuando se
obtiene una aproximación suficientemente buena. En este
caso se pide que:

donde es la cota suficiente.

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 hasta el nivel 4.
Tenemos la siguiente tabla:

Haciendo los cálculos de los errores, nos damos cuenta que


efectivamente la aproximación se obtiene hasta el nivel 4,
donde .
Por lo tanto, concluímos que la aproximación buscada es:

Ir al Indice
Msc. Pilar Nuñez Blas
Doc. Facultad de Ciencias de la UNASAM

También podría gustarte