Integración de Romberg

Universidad de Santiago de Chile
Facultad de Ingeniería
Departamento de Ingeniería Mecánica
Cálculo Numérico
Integración Numérica
INTEGRACIÓN NUMÉRICA
Sabemos que el problema de determinar el área bajo la curva de una función
f(x) conduce al concepto de integral, el cual conocemos de la forma:
𝑏
𝑓 π‘₯ 𝑑π‘₯
π‘Ž
El problema β€œse vuelve aún peor” cuando dicha integral es muy difícil de
resolver o, definitivamente, no se puede mediante métodos analíticos. Por
ejemplo, encuentre la solución para la siguiente integral:
1.2
1.1
(ex βˆ’ 1)sin π‘₯
𝑑π‘₯
2
tan π‘₯ βˆ’ π‘₯
INTEGRACIÓN NUMÉRICA
En este capítulo estudiaremos como resolver este tipo de integrales mediante
aproximaciones numéricas. Para ello estudiaremos los siguiente métodos:
β€’ Método del Trapecio: simple, corregido y compuesto
β€’ Método de Simpson: simple y compuesto
β€’ Método de Romberg
β€’ Cuadratura Gaussiana
INTEGRACIÓN NUMÉRICA
Representación Gráfica Método del Trapecio
Suponga que aproximamos una función mediante una línea recta entre los
puntos a y b.
1 panel entre a y b
INTEGRACIÓN NUMÉRICA
Si integramos al área bajo la recta con la que aproximamos f(x) obtenemos el
siguiente resultado.
Método del Trapecio
𝑏
𝑇=
π‘Ž
π‘βˆ’π‘Ž
𝑓 π‘₯ 𝑑π‘₯ =
𝑓 π‘Ž + 𝑓(𝑏) + 𝐸𝑇
2
Donde:
𝐸𝑇 = βˆ’
1
12
𝑏 βˆ’ π‘Ž 3 𝑓 β€²β€² (πœ‰); πœ‰ ∈ (π‘Ž, 𝑏)
es el error de cuadratura o integración.
INTEGRACIÓN NUMÉRICA
Una mejora al método del trapecio es la siguiente:
Método del Trapecio Corregido
𝑏
𝑇𝑐 =
π‘Ž
Donde:
π‘βˆ’π‘Ž
π‘βˆ’π‘Ž
𝑓 π‘₯ =
𝑓 π‘Ž + 𝑓(𝑏) βˆ’
2
12
𝐸𝑇𝐢 = βˆ’
1 π‘βˆ’π‘Ž 5
90
2
2
𝑓′ π‘Ž + 𝑓′ 𝑏
𝑓 β€²β€² (πœ‰); πœ‰ ∈ (π‘Ž, 𝑏)
es el error de cuadratura o integración.
+ 𝐸𝑇𝐢
INTEGRACIÓN NUMÉRICA
Representación Gráfica Método del Trapecio Compuesto
Aumenta el número de paneles con el fin de disminuir el error cometido.
INTEGRACIÓN NUMÉRICA
Método del Trapecio Compuesto
Con n número de paneles, podemos escribir el método del trapecio de la
siguiente forma:
𝑇 𝑗 =
π‘₯𝑛
π‘₯0
Donde:
y
β„Ž
𝑓 π‘₯ 𝑑π‘₯ =
𝑓 π‘₯0 + 2
2
π‘›βˆ’1
𝑓(π‘₯𝑖 ) + 𝑓(π‘₯𝑛 ) + 𝐸𝑇𝑗
𝑖=1
π‘₯𝑛 βˆ’ π‘₯0
β„Ž=
= π‘₯𝑖+1 βˆ’ π‘₯𝑖
𝑛
𝐸𝑇𝑗 = βˆ’
1 β„Ž3
12
𝑛2
es el error de cuadratura integración.
𝑓 β€²β€² (πœ‰); πœ‰ ∈ (π‘Ž, 𝑏)
INTEGRACIÓN NUMÉRICA
Si integramos al área bajo la curva cuadrática con la que aproximamos f(x)
obtenemos el siguiente resultado
Método de Simpson
𝑆=
π‘₯2
π‘₯0
Donde:
β„Ž
𝑓 π‘₯ 𝑑π‘₯ = 𝑓 π‘₯0 + 4𝑓(π‘₯1 ) + 𝑓(π‘₯2 ) + 𝐸𝑆
3
π‘₯2 βˆ’ π‘₯0
β„Ž=
= π‘₯2 βˆ’ π‘₯1 = π‘₯1 βˆ’ π‘₯0
2
y
𝐸𝑆 = βˆ’
1
90
β„Ž5 𝑓 𝑖𝑣 (πœ‰); πœ‰ ∈ (π‘Ž, 𝑏)
es el error de cuadratura o integración.
INTEGRACIÓN NUMÉRICA
Representación Gráfica Método de Simpson
Suponga que aproximamos una función mediante una curva cuadrática entre
los puntos a y b, tomando un punto intermedio.
INTEGRACIÓN NUMÉRICA
Realizando el mismo procedimiento utilizado para derivar la fórmula del
trapecio compuesto podemos obtener una fórmula para Simpson compuesto.
Método de Simpson Compuesto
𝑆 𝑗 =
π‘₯ 2π‘š
𝑓 π‘₯ 𝑑π‘₯
π‘₯0
β„Ž
=
𝑓 π‘₯0 + 2
3
π‘š βˆ’1
π‘š βˆ’1
𝑓 π‘₯2𝑖 + 4
𝑖=1
𝑓 π‘₯2π‘–βˆ’1 + 𝑓(π‘₯2π‘š ) + 𝐸𝑆𝑗
𝑖=1
π‘₯2π‘š βˆ’ π‘₯0
β„Ž=
= π‘₯𝑖+1 βˆ’ π‘₯𝑖
2π‘š
Donde:
π‘₯𝑖 = π‘₯0 + π‘–β„Ž
Y
𝐸𝑆𝑗 = βˆ’
1
90
β„Ž5
π‘š
𝑖𝑣
𝑓
(πœ‰);
𝑖=1
es el error de cuadratura o integración.
πœ‰ ∈ (π‘Ž, 𝑏)
INTEGRACIÓN NUMÉRICA
Retomemos el ejemplo del comienzo
1.2
1.1
Otros ejemplos ver Mathcad.
(ex βˆ’ 1)sin π‘₯
𝑑π‘₯
2
tan π‘₯ βˆ’ π‘₯
INTEGRACIÓN NUMÉRICA
Integración de Romberg
La fórmula de Romberg se puede escribir de la siguiente manera:
1
𝑗 β‰₯ π‘˜; 𝑅 𝑗, π‘˜ = 𝑅 𝑗, π‘˜ βˆ’ 1 βˆ’ π‘˜
(𝑅 𝑗, π‘˜ βˆ’ 1 βˆ’ 𝑅(𝑗 βˆ’ 1, π‘˜ βˆ’ 1))
4 βˆ’1
La idea es encontrar un método de integración que en forma iterativa
disminuya el error de integración
¿Qué indican los parámetros J y k?
j: Indica el número de paneles, es decir, la cantidad de subintervalos con los
que se desea aproximar una integral.
k: Indica el orden del método, por ejemplo, para k = 0 se está utilizando el
método del trapecio y para k = 1 se está utilizando el método de Simpson.
INTEGRACIÓN NUMÉRICA
Integración de Romberg
Del esquema de Romberg antes descrito se puede obtener lo siguiente:
π‘˜ = 0 βˆ’ π‘‡π‘Ÿπ‘Žπ‘π‘’π‘π‘–π‘œ βˆ’ 𝑅 𝑗, 0 = 𝑇 𝑗 , 𝑗 β‰₯ 0; 𝑁 = 2𝑗 β†’ π‘œ(β„Ž2 )
π‘˜ = 1 βˆ’ π‘†π‘–π‘šπ‘π‘ π‘œπ‘› βˆ’ 𝑅 𝑗, 1 = 𝑆 𝑗 , 𝑗 β‰₯ 1; 𝑁 = 2𝑗 β†’ π‘œ(β„Ž4 )
π‘˜ = 2 βˆ’ π΅π‘œπ‘œπ‘™π‘’ βˆ’ 𝑅 𝑗, 2 = 𝐡 𝑗 , 𝑗 β‰₯ 2; 𝑁 = 2𝑗 β†’ π‘œ(β„Ž6 )
A medida que aumentamos el parámetro β€œk” aumentamos el orden del error
INTEGRACIÓN NUMÉRICA
Integración de Romberg
¿Cómo opera Romberg?
Error
o(h2)
o(h4)
o(h6)
Trapecio
Simpson
Boole
R(J,1)
R(J,2)
Índice J
N° Paneles (2^J)
R(J,0)
0
1
T(0)=R(0,0)
1
2
T(1)=R(1,0)
R(1,1)
2
4
T(3)=R(3,0)
R(2,1)
R(2,2)
3
8
T(4)=R(4,0)
R(3,1)
R(3,2)
𝐸 = |𝑅 𝑗, 𝑗 βˆ’ 𝑅(𝑗 βˆ’ 1, 𝑗 βˆ’ 1)| < πœ€
Ejemplo:
o(h8)
R(J,3)
R(3,3)
INTEGRACIÓN NUMÉRICA
Cuadratura Gaussiana
En general, podemos escribir las fórmulas de cuadratura de la forma:
𝑛
𝑏
𝑓 π‘₯ 𝑑π‘₯ =
π‘Ž
𝑀𝑗 𝑓 π‘₯𝑗 + 𝐸
𝑗 =0
Donde wj se llaman pesos de cuadratura y xj nodos o puntos de precisión.
Por ejemplo, para el método del trapecio:
π‘βˆ’π‘Ž
𝑀𝑗 =
; 𝑗 = 1,2
2
π‘₯1 = π‘Ž y π‘₯2 = 𝑏
1
𝐸=βˆ’
𝑏 βˆ’ π‘Ž 3 𝑓 β€²β€² (πœ‰)
12
INTEGRACIÓN NUMÉRICA
Cuadratura Gaussiana
La cuadratura de Gauss escoge los puntos de precisión y los pesos de forma
conveniente con el objeto que el error de integración sea lo mas cercano a
cero.
Representación Gráfica:
INTEGRACIÓN NUMÉRICA
Cuadratura Gaussiana
La fórmula de Gauss que estudiaremos es la fórmula de Legendre, que se
escribe de la siguiente forma:
𝑛
1
𝑓 π‘₯ 𝑑π‘₯ =
βˆ’1
𝑀𝑗 𝑓 π‘₯𝑗 + 𝐸
𝑗 =0
Los pesos y puntos de precisión se obtienen de la siguiente tabla:
INTEGRACIÓN NUMÉRICA
Límites de integración de las fórmulas de cuadratura Gaussiana
Es posible notar que la fórmula de cuadratura Gaussiana como la hemos
definido sólo es aplicable a problemas con límites de integración -1 y 1.
Sin embargo, es posible hacer una cambio de variables de la forma:
𝑦 = π‘šπ‘₯ + 𝑛
Al encontrar las constantes m y n, se obtiene:
π‘βˆ’π‘Ž π‘₯+𝑏+π‘Ž
π‘βˆ’π‘Ž π‘₯
𝑦=
; 𝑑𝑦 =
2
2
Donde a y b son los nuevos límites de integración.
Ejemplo:
INTEGRACIÓN NUMÉRICA
Cotas de Error
De las fórmulas de error mencionadas para cada uno de los métodos
estudiados, se pueden obtener cotas de los errores y, con ello, poder definir el
espaciamiento h.
Consideremos el error para el método del trapecio:
Suponga que |f’’(x)| ≀ K para a ≀ x ≀ b, de aquí es posible escribir lo siguiente
1 3
|𝐸𝑇 | ≀ βˆ’ β„Ž 𝐾
12
Ejemplo:
INTEGRACIÓN NUMÉRICA
Exactitud de las fórmulas de integración
Consideremos la fórmula de error del trapecio
1
𝐸𝑇 = βˆ’
𝑏 βˆ’ π‘Ž 3 𝑓 β€²β€² (π‘₯)
12
Es posible observar que para polinomios de grado 1 la fórmula es cero para
todo x en [a,b].
Por lo tanto, podemos decir que la fórmula del trapecio es exacta hasta
polinomios de grado 1.
INTEGRACIÓN NUMÉRICA
Exactitud de las fórmulas de integración
Consideremos ahora la fórmula del error para el método de Simpson:
1 5 𝑖𝑣
𝐸𝑆 = βˆ’ β„Ž 𝑓 (π‘₯)
90
Es posible notar que para polinomios de grado 3 la fórmula del error es cero
para cualquier x en [a,b].
Por lo tanto, podemos decir que la fórmula de Simpson es exacta hasta
polinomios de grado 3.
En los casos de cuadratura Gaussiana, una fórmula de n sumandos es exacta
hasta polinomios de grado 2n+1.
INTEGRACIÓN NUMÉRICA
Integración Múltiple
Cada uno de los métodos anteriormente vistos se pueden emplear para
realizar integraciones múltiples.
Sabemos que una integral múltiple se resuelve en forma iterativa, es decir:
𝑏
𝑑
𝑏
𝑑
𝑓 π‘₯, 𝑦 𝑑π‘₯𝑑𝑦 =
π‘Ž
𝑐
𝑓 π‘₯, 𝑦 𝑑π‘₯ 𝑑𝑦
π‘Ž
𝑐
Lo cual también sirve para integrales de mayores dimensiones.
La idea es resolver en forma iterada cada integral aplicando los métodos
vistos.
Ejemplo: