Integración de Romberg�
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
Integración de Romberg
| O(h2) | | | | … |
n = 1 | I1,1 | | | | |
| | ||||
| | | |||
| | | | ||
| | | | | |
| | ||||
Integración de Romberg
| O(h2) | | | | … |
n = 1 | I1,1 | | | | |
n = 2 | I2,1 | ||||
| | | |||
| | | | ||
… | | | | | |
| | ||||
Integración de Romberg
| O(h2) | O(h4) | | | … |
n = 1 | I1,1 | I1,2 | | | |
n = 2 | I2,1 | ||||
| | | |||
| | | | ||
… | | | | | |
| | ||||
Integración de Romberg
| O(h2) | O(h4) | O(h6) | | … |
n = 1 | I1,1 | I1,2 | | | |
n = 2 | I2,1 | ||||
| | | |||
| | | | ||
… | | | | | |
| | ||||
Integración de Romberg
| O(h2) | O(h4) | O(h6) | | … |
n = 1 | I1,1 | I1,2 | | | |
n = 2 | I2,1 | ||||
n = 4 | I3,1 | | |||
| | | | ||
… | | | | | |
| | ||||
Integración de Romberg
| O(h2) | O(h4) | O(h6) | | … |
n = 1 | I1,1 | I1,2 | | | |
n = 2 | I2,1 | ||||
n = 4 | I3,1 | I2,2 | |||
| | | | ||
… | | | | | |
| | ||||
Integración de Romberg
| O(h2) | O(h4) | O(h6) | | … |
n = 1 | I1,1 | I1,2 | I1,3 | | |
n = 2 | I2,1 | ||||
n = 4 | I3,1 | I2,2 | |||
| | | | ||
… | | | | | |
| | ||||
Integración de Romberg
| 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 | |||
| | | | ||
… | | | | | |
| | ||||
Integración de Romberg
| 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 | | | ||
… | | | | | |
| | ||||
Integración de Romberg
| 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 | | ||
… | | | | | |
| | ||||
Integración de Romberg
| 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 | ||
… | | | | | |
| | ||||
Integración de Romberg
| 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?
Integración de Romberg
| 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?
Integración de Romberg
| 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
¿Cómo combinamos dos estimaciones?
Forma general de la Integración de Romberg:
Integración de Romberg
| 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
¿Cómo combinamos dos estimaciones?
Forma general de la Integración de Romberg:
CASO ESPECIAL:
dividiendo a la mitad sucesivamente
Integración de Romberg
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
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:
Integración de Romberg
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:
Integración de Romberg
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:
Integración de Romberg
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
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
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
Formas cerradas de las fórmulas de integración de Newton-Cotes
Comentarios finales
Comentario para el caso de integrales múltiples (dobles o triples…)
EJEMPLO
se pueden calcular como integrales iteradas
1° saplican métodos a la primera dimensión manteniendo constantes los valores de la segunda dimensión.
2° se aplica e el método para integrar la segunda dimensión.
Ejemplo aplicando la regla del trapecio con dos segmentos