Cuadratura Gaussiana

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 3

Cuadratura Gaussiana

La cuadratura Gaussiana consiste en seleccionar de manera óptima los nodos 𝑥1 , 𝑥2 , … , 𝑥𝑛 y las


constantes 𝑐1 , 𝑐2 , … , 𝑐𝑛 de tal manera que:
𝑏 𝑛

∫ 𝑓(𝑥)𝑑𝑥 = ∑ 𝑐𝑖 𝑓(𝑥𝑖 )
𝑎 𝑖=1

Para lograr la mayor exactitud de grado de precisión necesitamos un polinomio de grado 2𝑛 − 1.

Por ejemplo, para determinar la fórmula de cuadratura cuando 𝑛 = 2, es decir:


𝑏
∫ 𝑓(𝑥)𝑑𝑥 = 𝑐1 𝑓(𝑥1 ) + 𝑐2 𝑓 (𝑥2 )
𝑎

Utilizando un polinomio de grado 3, tenemos:


𝑏 𝑏
∫ 𝑓 (𝑥)𝑑𝑥 = ∫ (𝑎0 + 𝑎1 𝑥 + 𝑎2 𝑥 2 + 𝑎3 𝑥 3 )𝑑𝑥
𝑎 𝑎

𝑏 𝑏 𝑏 𝑏
= 𝑎0 ∫ 𝑑𝑥 + 𝑎1 ∫ 𝑥𝑑𝑥 + 𝑎2 ∫ 𝑥 2 𝑑𝑥 + 𝑎3 ∫ 𝑥 3 𝑑𝑥
𝑎 𝑎 𝑎 𝑎

Entonces la fórmula de la cuadratura exacta es exacta para 1, 𝑥, 𝑥 2 , 𝑥 3 . Por lo que se genera el


siguiente sistema utilizando el intervalo [−1,1]
1
∫ 𝑑𝑥 = 𝑐1 + 𝑐2 → 𝑐1 + 𝑐2 = 2
−1

1
∫ 𝑥𝑑𝑥 = 𝑐1 𝑥1 + 𝑐2 𝑥2 → 𝑐1 𝑥1 + 𝑥2 𝑐2 = 0
−1
1
2
∫ 𝑥 2 𝑑𝑥 = 𝑐1 𝑥12 + 𝑐2 𝑥22 → 𝑐1 𝑥12 + 𝑐2 𝑥22 =
−1 3
1
∫ 𝑥 3 𝑑𝑥 = 𝑐1 𝑥13 + 𝑐2 𝑥23 → 𝑐1 𝑥13 + 𝑐2 𝑥23 = 0
−1

1
Resolviendo el sistema de 4 × 4 nos queda como solución única 𝑐1 = 𝑐2 = 1, 𝑥1 = √3,
1
𝑥2 = −√3

Por lo tanto, la fórmula de cuadratura para 𝑛 = 2 es:


1
√3 −√3
∫ 𝑓(𝑥)𝑑𝑥 = 𝑓 ( ) + 𝑓 ( )
−1 3 3

Para un intervalo arbitrario [𝑎, 𝑏] utilizaremos el cambio de variable:


2𝑥 − 𝑎 − 𝑏
𝑢=
𝑏−𝑎
𝑢 (𝑏 − 𝑎) + 𝑎 + 𝑏
𝑥=
2
𝑏 1
𝑢 (𝑏 − 𝑎) + 𝑎 + 𝑏 𝑏 − 𝑎
∫ 𝑓 (𝑥)𝑑𝑥 = ∫ 𝑓 ( ) 𝑑𝑢
𝑎 −1 2 2

Ejemplo: Aproxime la siguiente integral usando la cuadratura gaussiana con 𝑛 = 2


1.5
∫ 𝑥 2 ln(𝑥)𝑑𝑥
1

Solución:

El cambio de variable
𝑢 (𝑏 − 𝑎) + 𝑎 + 𝑏
𝑥=
2

0.5𝑢 + 2.5
𝑥=
2

1.5
2
0.5𝑢 + 2.5 2
1
0.5𝑢 + 2.5 0.5
∫ 𝑥 ln(𝑥)𝑑𝑥 = ∫ ( ) ln ( ) 𝑑𝑢
1 −1 2 2 2
2
√3 √3
0.5 ( 3 ) + 2.5 0.5 ( 3 ) + 2.5
0.5
= ln
2 2 2
[( ) ( )
2
−√3 −√3
0.5 ( ) + 2.5 0.5 ( ) + 2.5
3 3
+ ln
2 2
( ) ( )]
= 0.192268
1.5
Real ∫1 𝑥 2 ln(𝑥)𝑑𝑥 = 0.192259
Ejemplo: Considere la fórmula de cuadratura:
1
∫ 𝑒 𝑥 𝑓(𝑥)𝑑𝑥 = 𝐴𝑓 (0) + 𝐵𝑓 ′ (0) + 𝐶𝑓(1)
0

Encuentre el valor de las constantes 𝐴, 𝐵 y 𝐶 para que la fórmula tenga el mayor grado de
exactitud. Determinar el grado.

Solución:

Sea 𝑃 (𝑥) = 𝑥 𝑘 con 𝑘 ∈ ℤ+

Si 𝑘 = 0 → 𝑃(𝑥) = 1
1
∫ 𝑒 𝑥 𝑑𝑥 = 𝐴(1) + 𝐵 (0) + 𝐶 (1) → 𝐴 + 𝐶 = 𝑒 − 1
0

Si 𝑘 = 1 → 𝑃(𝑥) = 𝑥
1
∫ 𝑒 𝑥 𝑥𝑑𝑥 = 𝐴(0) + 𝐵(1) + 𝐶 (1) → 𝐵 + 𝐶 = 1
0

Si 𝑘 = 2 → 𝑃(𝑥) = 𝑥 2
1
∫ 𝑒 𝑥 𝑥 2 𝑑𝑥 = 𝐴(0) + 𝐵(0) + 𝐶 (1) → 𝐶 = 𝑒 − 2
0

𝐴+𝐶 =𝑒−1
{ 𝐵+𝐶 =1
𝐶 = 𝑒−2

Entonces 𝐴 = 1, 𝐵 = 3 − 𝑒 y 𝐶 = 𝑒 − 2

Si 𝑘 = 3 → 𝑃(𝑥) = 𝑥 3
1
∫ 𝑒 𝑥 𝑥 3 𝑑𝑥 = 𝐴(0) + 𝐵(0) + 𝐶 (1) → 𝐶 = −2𝑒 + 6
0

Que contradice el procedimiento anterior, de manera que 𝑘 = 2

También podría gustarte