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:
© Copyright 2024