1 of 27

Integración de Romberg�

2 of 27

Métodos de aproximación numérica

¿Problemática?

Necesitamos integrar una función continua complicada o una función tabulada donde los valores de x y f(x) están dados como un conjunto discreto de puntos.

Usamos métodos de integración numéricos para obtener la integral aproximada

Formas cerradas de las fórmulas de integración de Newton-Cotes

Integración de Romberg

Regla de Simpson 1/3

Regla de Simpson 3/8

Regla del trapecio

3 of 27

  • Es una técnica de Integración de ecuaciones

Integración de Romberg

  • Se basa en la extrapolación de Richardson que es un método es un método que combina dos estimaciones numéricas de la integral para obtener un tercer valor más exacto.
  • Es una técnica recursiva, se basa en aplicaciones sucesivas de la regla del trapecio
  • Se utiliza para generar una estimación/aproximación de la integral dentro de una tolerancia de error preespecificada

4 of 27

O(h2)

n = 1

I1,1

  • Técnica recursiva … combina dos estimaciones/aproximaciones numéricas de la integral para obtener un tercer valor más exacto.

  • Esas dos estimaciones o aproximaciones se obtienen de aplicar la regla del trapecio

Integración de Romberg

5 of 27

O(h2)

n = 1

I1,1

n = 2

I2,1

  • Técnica recursiva … combina dos estimaciones/aproximaciones numéricas de la integral para obtener un tercer valor más exacto.

  • Esas dos estimaciones o aproximaciones se obtienen de aplicar la regla del trapecio

Integración de Romberg

6 of 27

O(h2)

O(h4)

n = 1

I1,1

I1,2

n = 2

I2,1

  • Técnica recursiva … combina dos estimaciones/aproximaciones numéricas de la integral para obtener un tercer valor más exacto.

  • Esas dos estimaciones o aproximaciones se obtienen de aplicar la regla del trapecio

Integración de Romberg

7 of 27

O(h2)

O(h4)

O(h6)

n = 1

I1,1

I1,2

n = 2

I2,1

  • Técnica recursiva … combina dos estimaciones/aproximaciones numéricas de la integral para obtener un tercer valor más exacto.

  • Esas dos estimaciones o aproximaciones se obtienen de aplicar la regla del trapecio

Integración de Romberg

8 of 27

O(h2)

O(h4)

O(h6)

n = 1

I1,1

I1,2

n = 2

I2,1

n = 4

I3,1

  • Técnica recursiva … combina dos estimaciones/aproximaciones numéricas de la integral para obtener un tercer valor más exacto.

  • Esas dos estimaciones o aproximaciones se obtienen de aplicar la regla del trapecio

Integración de Romberg

9 of 27

O(h2)

O(h4)

O(h6)

n = 1

I1,1

I1,2

n = 2

I2,1

n = 4

I3,1

I2,2

  • Técnica recursiva … combina dos estimaciones/aproximaciones numéricas de la integral para obtener un tercer valor más exacto.

  • Esas dos estimaciones o aproximaciones se obtienen de aplicar la regla del trapecio

Integración de Romberg

10 of 27

O(h2)

O(h4)

O(h6)

n = 1

I1,1

I1,2

I1,3

n = 2

I2,1

n = 4

I3,1

I2,2

  • Técnica recursiva … combina dos estimaciones/aproximaciones numéricas de la integral para obtener un tercer valor más exacto.

  • Esas dos estimaciones o aproximaciones se obtienen de aplicar la regla del trapecio

Integración de Romberg

11 of 27

O(h2)

O(h4)

O(h6)

O(h8)

n = 1

I1,1

I1,2

I1,3

n = 2

I2,1

n = 4

I3,1

I2,2

  • Técnica recursiva … combina dos estimaciones/aproximaciones numéricas de la integral para obtener un tercer valor más exacto.

  • Esas dos estimaciones o aproximaciones se obtienen de aplicar la regla del trapecio

Integración de Romberg

12 of 27

O(h2)

O(h4)

O(h6)

O(h8)

n = 1

I1,1

I1,2

I1,3

n = 2

I2,1

n = 4

I3,1

I2,2

n = 8

I4,1

  • Técnica recursiva … combina dos estimaciones/aproximaciones numéricas de la integral para obtener un tercer valor más exacto.

  • Esas dos estimaciones o aproximaciones se obtienen de aplicar la regla del trapecio

Integración de Romberg

13 of 27

O(h2)

O(h4)

O(h6)

O(h8)

n = 1

I1,1

I1,2

I1,3

n = 2

I2,1

n = 4

I3,1

I2,2

n = 8

I4,1

I3,2

  • Técnica recursiva … combina dos estimaciones/aproximaciones numéricas de la integral para obtener un tercer valor más exacto.

  • Esas dos estimaciones o aproximaciones se obtienen de aplicar la regla del trapecio

Integración de Romberg

14 of 27

O(h2)

O(h4)

O(h6)

O(h8)

n = 1

I1,1

I1,2

I1,3

n = 2

I2,1

n = 4

I3,1

I2,2

n = 8

I4,1

I3,2

I2,3

  • Técnica recursiva … combina dos estimaciones/aproximaciones numéricas de la integral para obtener un tercer valor más exacto.

  • Esas dos estimaciones o aproximaciones se obtienen de aplicar la regla del trapecio

Integración de Romberg

15 of 27

O(h2)

O(h4)

O(h6)

O(h8)

n = 1

I1,1

I1,2

I1,3

I1,4

n = 2

I2,1

n = 4

I3,1

I2,2

n = 8

I4,1

I3,2

I2,3

¿Cuál es el criterio de terminación de las iteraciones?

  • Técnica recursiva … combina dos estimaciones/aproximaciones numéricas de la integral para obtener un tercer valor más exacto.

  • Esas dos estimaciones o aproximaciones se obtienen de aplicar la regla del trapecio

Integración de Romberg

16 of 27

O(h2)

O(h4)

O(h6)

O(h8)

n = 1

I1,1

I1,2

I1,3

I1,4

n = 2

I2,1

n = 4

I3,1

I2,2

n = 8

I4,1

I3,2

I2,3

Las filas j indican integrales considerando c/vez con mas segmentos

Las columnas k, indican el nivel de la integración

Estimación del error relativo porcentual:

¿Cuál es el criterio de terminación de las iteraciones?

  • Técnica recursiva … combina dos estimaciones/aproximaciones numéricas de la integral para obtener un tercer valor más exacto.

  • Esas dos estimaciones o aproximaciones se obtienen de aplicar la regla del trapecio

Integración de Romberg

17 of 27

O(h2)

O(h4)

O(h6)

O(h8)

n = 1

I1,1

I1,2

I1,3

I1,4

n = 2

I2,1

n = 4

I3,1

I2,2

n = 8

I4,1

I3,2

I2,3

  • Técnica recursiva … combina dos estimaciones/aproximaciones numéricas de la integral para obtener un tercer valor más exacto.

  • Esas dos estimaciones o aproximaciones se obtienen de aplicar la regla del trapecio

Las filas j indican integrales considerando c/vez con mas segmentos

Las columnas k, indican el nivel de la integración

¿Cómo combinamos dos estimaciones?

Forma general de la Integración de Romberg:

Integración de Romberg

18 of 27

O(h2)

O(h4)

O(h6)

O(h8)

n = 1

I1,1

I1,2

I1,3

I1,4

n = 2

I2,1

n = 4

I3,1

I2,2

n = 8

I4,1

I3,2

I2,3

  • Técnica recursiva … combina dos estimaciones/aproximaciones numéricas de la integral para obtener un tercer valor más exacto.

  • Esas dos estimaciones o aproximaciones se obtienen de aplicar la regla del trapecio

Las filas j indican integrales considerando c/vez con mas segmentos

Las columnas k, indican el nivel de la integración

¿Cómo combinamos dos estimaciones?

Forma general de la Integración de Romberg:

CASO ESPECIAL:

dividiendo a la mitad sucesivamente

Integración de Romberg

19 of 27

Aproximación de la integral con O(h6)

Aproximación de la integral con O(h8)

Aproximación de la integral con O(h4)

Los coeficientes suman 1

Representan factores de ponderación

que, conforme aumenta la exactitud, dan un peso relativamente mayor a la mejor estimación de la integral.

¿Cómo combinamos dos estimaciones?

Forma general de la Integración de Romberg:

CASO ESPECIAL:

dividiendo a la mitad sucesivamente

20 of 27

EJEMPLO

Mejorar la estimación de la integral de la siguiente función entre 0 y 0.8 utilizando los resultados de estimaciones simples y múltiples de la regla del trapecio:

  • Si se combinan las estimaciones con 1 y 2 segmentos resulta:

Integración de Romberg

21 of 27

EJEMPLO

Mejorar la estimación de la integral de la siguiente función entre 0 y 0.8 utilizando los resultados de estimaciones simples y múltiples de la regla del trapecio:

  • Si se combinan las estimaciones con 1 y 2 segmentos resulta:
  • También podemos utilizar las estimaciones con 2 y 4 segmentos para obtener:

Integración de Romberg

22 of 27

EJEMPLO

Mejorar la estimación de la integral de la siguiente función entre 0 y 0.8 utilizando los resultados de estimaciones simples y múltiples de la regla del trapecio:

  • Si se combinan las estimaciones con 1 y 2 segmentos resulta:
  • También podemos utilizar las estimaciones con 2 y 4 segmentos para obtener:

Integración de Romberg

23 of 27

EJEMPLO

Mejorar la estimación de la integral de la siguiente función entre 0 y 0.8 utilizando los resultados de dos estimaciones de la integral O(h4):

Combinan las estimaciones O(h4) y calcular una integral O(h6) :

Integración de Romberg

24 of 27

EJEMPLO

Proceso iterativo de aproximación:

Estimación del error relativo porcentual:

¿Cuál es el criterio de terminación de las iteraciones?

Integración de Romberg

25 of 27

Métodos de aproximación numérica

¿Problemática?

Necesitamos integrar una función continua complicada o una función tabulada donde los valores de x y f(x) están dados como un conjunto discreto de puntos.

Usamos métodos de integración numéricos para obtener la integral aproximada

Formas cerradas de las fórmulas de integración de Newton-Cotes

Integración de Romberg

Regla de Simpson 1/3

Regla de Simpson 3/8

Regla del trapecio

26 of 27

Formas cerradas de las fórmulas de integración de Newton-Cotes

Comentarios finales

  • La integración de Romberg es más eficientes que las reglas del trapecio y de Simpson.

  • La reglas de Simpson, las integración de Romberg y la cuadratura de Gauss son más eficientes y exactas que la regla del trapecio.

  • Para la regla del Trapecio y las de Simpson, de aplicación múltiple, conforme número de segmentos aumenta (n) disminuye el error. Pero también, para grandes valores de n, el error vuelve a aumentar debido a que los errores de redondeo comienzan a dominar.

27 of 27

Comentario para el caso de integrales múltiples (dobles o triples…)

EJEMPLO

se pueden calcular como integrales iteradas

saplican métodos a la primera dimensión manteniendo constantes los valores de la segunda dimensión.

se aplica e el método para integrar la segunda dimensión.

Ejemplo aplicando la regla del trapecio con dos segmentos