Integración Por Cuadratura de Gauss

Profesor Alejandro Sandoval Ramos
INTEGRACIÓN NUMÉRICA POR EL MÉTODO DE CUADRATURA DE GAUSS
En la figura se ilustra la curva de la función f(x) que se debe integrar entre los límites a y b. En el
inciso a se muestra como se integraría usando el método trapezoidal, se observa que una
importante cantidad de área no es tomada en cuenta por el trapezoide de vértices (𝑥0 , 𝐴, 𝐵, 𝑥1 ) por
lo cual la aproximación a la integral por este método nos genera un error significativo.
Nieves & Domínguez (2014)
El inciso b nos muestra el método de Gauss con dos puntos para aproximar la integral de la función
𝑓(𝑥), en dicho método en lugar de tomar los dos puntos A y B en los extremos del intervalo se
escogen dos puntos interiores C y D, se traza una línea recta que pase por estos dos puntos y se
extiende hasta los extremos del intervalo para formar el trapezoide sombreado. Como se puede
observar en la figura parte del trapezoide queda por encima de la curva y parte por abajo, si los
puntos C y D se escogen adecuadamente se pueden llegar a igualar las dos zonas de modo que el
área del trapezoide sea igual al área bajo la curva.
El método de Cuadratura de Gauss puede extenderse a tres o mas puntos, de manera general el
algoritmo tiene la forma:
𝑏 𝑛
𝑏−𝑎 𝑏−𝑎 𝑎+𝑏
∫ 𝑓(𝑥)𝑑𝑥 ≈ ∑ 𝑤𝑖 𝑓 ( 𝑧𝑖 + )
2 2 2
𝑎 𝑖=1
𝐷𝑜𝑛𝑑𝑒
𝑎 𝑒𝑠 𝑒𝑙 𝑙í𝑚𝑖𝑡𝑒 𝑖𝑛𝑓𝑒𝑟𝑖𝑜𝑟 𝑑𝑒 𝑖𝑛𝑡𝑒𝑔𝑟𝑎𝑐𝑖ó𝑛
𝑏 𝑒𝑠 𝑒𝑙 𝑙í𝑚𝑖𝑡𝑒 𝑠𝑢𝑝𝑒𝑟𝑖𝑜𝑟 𝑑𝑒 𝑖𝑛𝑡𝑒𝑔𝑟𝑎𝑐𝑖ó𝑛
𝑤𝑖 𝑒𝑠 𝑒𝑙 𝑣𝑎𝑙𝑜𝑟 𝑑𝑒𝑙 𝑐𝑜𝑒𝑓𝑖𝑐𝑖𝑒𝑛𝑡𝑒 𝑑𝑒 𝑙𝑎 𝑐𝑢𝑎𝑑𝑟𝑎𝑡𝑢𝑟𝑎 𝑑𝑒 𝐺𝑎𝑢𝑠𝑠
𝑧𝑖 𝑒𝑠 𝑒𝑙 𝑣𝑎𝑙𝑜𝑟 𝑑𝑒 𝑙𝑎𝑠 𝑎𝑏𝑠𝑐𝑖𝑠𝑎𝑠 𝑑𝑒 𝑙𝑎 𝑐𝑢𝑎𝑑𝑟𝑎𝑡𝑢𝑟𝑎 𝑑𝑒 𝐺𝑎𝑢𝑠𝑠
1
Academia de Matemáticas Aplicadas ESIQIE
Los valores de los coeficientes y las abscisas se muestran en la siguiente tabla y varían de acuerdo
al número de puntos que se emplean.
Nieves & Domínguez (2014)
Ejemplo 1. Calcule el valor de la siguiente integral usando el método de la cuadratura de Gauss con
2 puntos.
1.2
2
∫ 𝑒 𝑥 𝑑𝑥
0.2
𝑏 𝑛
∫ 𝑓(𝑥)𝑑𝑥 ≈ ∑ 𝑤𝑖 𝑓 ( 𝑧𝑖 + )
2 2 2
𝑎 𝑖=1
Solución
Como n=2
𝑏 2
∫ 𝑓(𝑥)𝑑𝑥 ≈ ∑ 𝑤𝑖 ∙ 𝑓 ( 𝑧𝑖 + )
2 2 2
𝑎 𝑖=1
𝑏
𝑏−𝑎 𝑏−𝑎 𝑎+𝑏 𝑏−𝑎 𝑎+𝑏
∫ 𝑓(𝑥)𝑑𝑥 ≈ (𝑤1 ∙ 𝑓 ( 𝑧1 + ) + 𝑤2 ∙ 𝑓 ( 𝑧2 + ))
2 2 2 2 2
𝑎
Para la cuadratura de Gauss de dos puntos los valores de w y z son
𝑤1 = 𝑤2 = 1; 𝑧1 = −0.5773502692; 𝑧2 = 0. 5773502692
Profesor Alejandro Sandoval Ramos
Se realizan la sustitución de estos valores obteniendo
1.2 − 0.2 1.2 − 0.2 0.2 + 1.2

𝐼≈ (1 ∙ 𝑓 ( (−0.5773502692) + )+1
2 2 2
1.2 − 0.2 0.2 + 1.2
∙𝑓( (0.5773502692) + ))
2 2
𝐼 ≈ 0.5(1 ∙ 𝑓(0.41132) + 1 ∙ 𝑓(0.98868))
Los valores obtenidos 0.41132 y 0.98868 se utilizan para evaluar la función que es
2
𝑓(𝑥) = 𝑒 𝑥
2
𝑓(0.41132) = 𝑒 0.41132
𝑓(0.41132) = 1.18434
2
𝑓(0.98868) = 𝑒 0.98868
𝑓(0.98868) = 2.65777
Sustituimos los valores de la función evaluada obteniendo
𝐼 ≈ 0.5(1 ∙ 1.18434 + 1 ∙ 2.65777)

𝐼 ≈ 1.92106
Ejemplo 2. Calcule el cambio de entalpia especifica que sufre el aire al pasar de 430°C a 100°C
𝑇2
∆𝐻 = ∫ 𝐶𝑝(𝑇)𝑑𝑇
𝑇1
La capacidad calorífica a presión constante del aire en este rango de temperatura está dada por
𝐶𝑝 = 0.02894 + 0.4147 × 10−5 𝑇 + 0.3191 × 10−8 𝑇 2 − 1.965 × 10−12 𝑇 3

𝐷𝑜𝑛𝑑𝑒
𝑇 = °𝐶
𝐾𝐽
𝐶𝑝 = 𝑚𝑜𝑙 °𝐶
Use para su cálculo la cuadratura de Gauss de tres puntos
Solución
Aplicamos la formula de cuadratura de Gauss para n=3

100 3
∆𝐻 = ∫ 𝐶𝑝(𝑇)𝑑𝑇 ≈ ∑ 𝑤𝑖 ∙ 𝑓 ( 𝑧𝑖 + )
2 2 2
400 𝑖=1
3
Academia de Matemáticas Aplicadas ESIQIE
𝑏−𝑎 𝑏−𝑎 𝑎+𝑏 𝑏−𝑎 𝑎+𝑏 𝑏−𝑎 𝑎+𝑏

∆𝐻 ≈ (𝑤1 ∙ 𝑓 ( 𝑧1 + ) + 𝑤2 ∙ 𝑓 ( 𝑧2 + ) +𝑤3 ∙ 𝑓 ( 𝑧3 + ))
2 2 2 2 2 2 2
Para la cuadratura de Gauss de tres puntos los valores de w y z son
𝑤1 = 𝑤3 = 0.55555; 𝑤2 = 0.88888; 𝑧1 = −0.7745966692; 𝑧2 = 0; 𝑧3 = 0.7745966692

Se realizan la sustitución de estos valores
100 − 430 100 − 430 430 + 100

∆𝐻 ≈ (0.55555 ∙ 𝑓 ( (−0.7745966692) + ) + 0.88888
2 2 2
100 − 430 430 + 100
∙𝑓( (0) + ) + 0.55555
2 2
100 − 430 430 + 100
∙𝑓( (0.7745966692) + ))
2 2
∆𝐻 ≈ −165(0.55555 ∙ 𝑓(392.80845) + 0.88888 ∙ 𝑓(265) + 0.55555 ∙ 𝑓(137.19155))
Se evalúa la función 𝐶𝑝(𝑇) con los valores obtenidos
𝐶𝑝(𝑇) = 0.02894 + 0.4147 × 10−5 𝑇 + 0.3191 × 10−8 𝑇 2 − 1.965 × 10−12 𝑇 3
𝐶𝑝(𝑇 = 392.80845) = 0.02894 + 0.4147 × 10−5 𝑇 + 0.3191 × 10−8 𝑇 2 − 1.965 × 10−12 𝑇 3

𝐶𝑝(𝑇 = 392.80845) = 0.03094
𝐶𝑝(𝑇 = 265) = 0.03023
𝐶𝑝(𝑇 = 137.19155) = 0.02956
Sustituimos los valores de la función evaluada en la ecuación de la ∆𝐻
∆𝐻 ≈ −165(0.55555 ∙ 0.03094 + 0.88888 ∙ 0.03023 + 0.55555 ∙ 0.02956)

∆𝐻 ≈ −165(0.55555 ∙ 0.03094 + 0.88888 ∙ 0.03023 + 0.55555 ∙ 0.02956)
∆𝐻 ≈ −9.97951
REFERENCIAS
Nieves, A. & Domínguez, F. (2014) Métodos Numéricos aplicados a la Ingeniería. Cuarta Edición.
Editorial Patria.
Márquez, R. & Vázquez, G. (2020) Problemario de Métodos Numéricos. Academia de Matemáticas
Aplicadas-ESIQIE-IPN.

Integración Por Cuadratura de Gauss

Cargado por

Copyright:

Formatos disponibles

Integración Por Cuadratura de Gauss

Cargado por

Información del documento

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Integración Por Cuadratura de Gauss

Cargado por

Copyright:

Formatos disponibles

Profesor Alejandro Sandoval Ramos

INTEGRACIÓN NUMÉRICA POR EL MÉTODO DE CUADRATURA DE GAUSS

Nieves & Domínguez (2014)

Nieves & Domínguez (2014)

Para la cuadratura de Gauss de dos puntos los valores de w y z son

Se realizan la sustitución de estos valores obteniendo

1.2 − 0.2 1.2 − 0.2 0.2 + 1.2

𝐼 ≈ 0.5(1 ∙ 𝑓(0.41132) + 1 ∙ 𝑓(0.98868))

𝐼 ≈ 0.5(1 ∙ 1.18434 + 1 ∙ 2.65777)

𝐶𝑝 = 0.02894 + 0.4147 × 10−5 𝑇 + 0.3191 × 10−8 𝑇 2 − 1.965 × 10−12 𝑇 3

Use para su cálculo la cuadratura de Gauss de tres puntos

Aplicamos la formula de cuadratura de Gauss para n=3

𝑏−𝑎 𝑏−𝑎 𝑎+𝑏 𝑏−𝑎 𝑎+𝑏 𝑏−𝑎 𝑎+𝑏

Para la cuadratura de Gauss de tres puntos los valores de w y z son

𝑤1 = 𝑤3 = 0.55555; 𝑤2 = 0.88888; 𝑧1 = −0.7745966692; 𝑧2 = 0; 𝑧3 = 0.7745966692

100 − 430 100 − 430 430 + 100

∆𝐻 ≈ −165(0.55555 ∙ 𝑓(392.80845) + 0.88888 ∙ 𝑓(265) + 0.55555 ∙ 𝑓(137.19155))

Se evalúa la función 𝐶𝑝(𝑇) con los valores obtenidos

𝐶𝑝(𝑇) = 0.02894 + 0.4147 × 10−5 𝑇 + 0.3191 × 10−8 𝑇 2 − 1.965 × 10−12 𝑇 3

𝐶𝑝(𝑇 = 392.80845) = 0.02894 + 0.4147 × 10−5 𝑇 + 0.3191 × 10−8 𝑇 2 − 1.965 × 10−12 𝑇 3

∆𝐻 ≈ −165(0.55555 ∙ 0.03094 + 0.88888 ∙ 0.03023 + 0.55555 ∙ 0.02956)

También podría gustarte