Semana 12

METODOS NUMERICOS
Facultad de Ingeniería Electrónica
semestre 2024-II

Mag Jose Luis Acuña Guillermo 1

 Mag Jose Luis Acuña Guillermo

CUADRATURA
GAUSSIANA

2
 Mag Jose Luis Acuña Guillermo

Todas las fórmulas previas para integración

numérica estaban basadas en valores
uniformemente espaciados: esto significa
que los valores estaban predeterminados.
Así, con una fórmula de tres términos hay
tres parámetros, donde los coeficientes
(factores de ponderación) se aplican a
cada uno de los valores de función. Una
fórmula con tres parámetros corresponde a
un polinomio de segundo grado, uno
menos que el número de parámetros.

3
 Mag Jose Luis Acuña Guillermo

Gauss observó que si se elimina el

requisito de que la función sea evaluada
en valores predeterminados, una
fórmula de tres términos contiene seis
parámetros (los tres valores son ahora
incógnitas, más los tres pesos) y debe
corresponder a un polinomio de
interpolación de grado 5. Las fórmulas
basadas en este principio se denominan
fórmulas de cuadratura gaussiana.
Pueden aplicarse sólo cuando se conoce
f(x) explícitamente, de modo que
pueden evaluarse en cualquier valor
deseado de x.
4
 Mag Jose Luis Acuña Guillermo

Introduciendo una función peso wC[a,b]

con w(x)>0, x(a,b), nos planteamos el
problema de encontrar (n+1) coeficientes
{c0, c1,c2,….,cn} tales que

 w ( x ) f ( x ) dx  ci f ( xi )
b

a
i =1 (1)

eligiendo los puntos {x0, x1,x2,….,xn}[a,b]

que maximicen el orden del método.

5
 Mag Jose Luis Acuña Guillermo

POLINOMIOS DE LEGENDRE

Los primeros cinco polinomios de Legendre son :

P0 ( x ) = 1
P1 ( x ) = x

P2 ( x ) = ( 3x 2 − 1)
1
2
P3 ( x ) = ( 5 x 3 − 3x )
1
2
P4 ( x ) = ( 35 x 4 − 30 x 2 + 3)
1
8
6
 Mag Jose Luis Acuña Guillermo

El polinomio de Legendre de grado n se puede obtener por

medio de la fórmula de

n (
x − 1)
1 d
Pn ( x ) = n 2 n

2 n ! dx

O bien a partir de la fórmula recursiva:

( n + 1) Pn+1 ( x ) − ( 2n + 1) xPn ( x ) + nPn −1 ( x ) = 0

7
 Mag Jose Luis Acuña Guillermo

P2 ( x ) = 3x 2 – 1 /2

P4 ( x ) = 35x 4 – 30x 2 + 3 /8

P3 ( x ) = 5x 3 – 3x /2

POLINOMIOS DE LEGENDRE 8
 Mag Jose Luis Acuña Guillermo

LOS POLINOMIOS DE CHEBYSHEV

T0 ( x ) = 1
T1 ( x ) = x
T2 ( x ) = 2 x − 1
2

T3 ( x ) = 4 x − 3x
3

T4 ( x ) = 8 x − 8 x + 1
4 2

T5 ( x ) = 16 x − 20 x + 5 x
5 3

9
 Mag Jose Luis Acuña Guillermo

El polinomio de Chebyshev de grado

(n+1) se puede obtener por medio de la
fórmula de recurrencia:

Tn +1 ( x ) = 2 xTn ( x ) − Tn −1

10
 Mag Jose Luis Acuña Guillermo

Polinomios de Chebyshev
11
 Mag Jose Luis Acuña Guillermo

Polinomios ortogonales

Polinomios de Intervalo Peso w() Polinomio genérico

1 dk 2
[-1,1] 1 ( )
k
Legendre Pk ( x ) = k k
x −1
2 k ! dx

1
Tchebychev [-1,1] Tk ( x ) = cos ( k arccos ( x ) )
1− x 2

−x
Laguerre [0,+) e Lk ( x ) = e
d k −x k
dx
x
k
e x ( )
− x2 d k  − x2 
H k ( x ) = ( −1) e x2
k
k 
e 
Hermite (-,+) e dx  

12
 Mag Jose Luis Acuña Guillermo

Los primeros polinomios de Legendre, Tchebychev, Laguerre y Hermite

Polinomios K=0 K=1 K=2 K=3

3x 2 − 1 x ( 5 x 2 − 3)
Pk ( x ) 1 x 2 2
Tk ( x ) 1 x 2 x2 −1 x ( 4 x 2 − 3)

Lk ( x ) 1 1− x 2 − 4x + x 2 − x3 + 9 x2 − 18x + 6

Hk ( x) 1 2x 4 x2 − 2 (
4 x 2 x2 − 3 )
13
 Mag Jose Luis Acuña Guillermo

Teorema

Si f  C
2n+ 2
( a, b) entonces existe   ( a, b ) tal
que

f(
b 2n+ 2)
()
n


b


a
w ( x ) f ( x ) dx =
i =0
ci f ( xi ) +
( 2n + 2 ) !  a
w ( x ) Pn+12dx (*)

donde las abscisas  x0 , x1, x2 ,...., xn  son las

raíces del polinomio ortogonal Pn+1 en el
intervalo [a,b] respecto al peso w(x) y
cc0 , c1, c2 ,...., cn son los coeficientes .
14
 Mag Jose Luis Acuña Guillermo

Ejemplo 1

Para aproximar, mediante una fórmula de Gauss

de tres puntos, la integral

 −1
e − x2
dx,

vamos a utilizar la fórmula de Gauss-Legendre y

Gauss-Tchebychev.

15
 Mag Jose Luis Acuña Guillermo

a) Las raíces del polinomio de Legendre P3 son

3 3
x0 = − , x1 =0 y x2 =
5 5
y los coeficientes valen
5 8 5
c0 = , c1 = y c2 =
9 9 9
por lo tanto, la aproximación obtenida es:
1 − x2 − x02 − x12 − x22
−1 e dx  c0e + c1e + c2e

1
 ( 5  0.548812 + 8 + 5  0.548812 )
9
 1.498679
16
 Mag Jose Luis Acuña Guillermo

b) Cuadratura de Gauss Chebyshev

1
f ( x) n

 dx   wi f ( xi )
−1 1− x 2
i =1

En este caso los pesos correspondientes son


wi = , i = 1, 2,..., n
n
17
 Mag Jose Luis Acuña Guillermo

es decir, la fórmula de Gauss-Chebyshev de n puntos es:

 n
Gn ( f ) 
n
 f (x )
i =1
i

 ( 2i − 1)  
xi = cos   , i = 1, 2,..., n
 2n 

18
 Mag Jose Luis Acuña Guillermo

Las raíces del polinomio de Tchebychev T3 son

3 3
x0 = − , x1 =0 y x2 =
2 2
y los coeficientes valen ahora


c0 = c1 = c2 =
3

− x2
dx =  w ( x )e − x2
1 1
−1 e
−1
dx

19
 ( )
Mag Jose Luis Acuña Guillermo

1 1 2
= 1− x e 2 −x
dx
−1
1− x 2

2 2 2
2 − x0 2 − x1 2 − x2
 c0 1 − x e0
+ c1 1 − x e 1
+ c2 1 − x e
2


 ( 0.236183 + 1 + 0.236183)
3
 1.541858

20
 Mag Jose Luis Acuña Guillermo

Para la fórmula de Gauss-Legendre se puede

utilizar, de forma sencilla, la expresión (*) con
n=2 para el estudio del error. Puesto que la
derivada sexta de la función − x 2 es
f ( x) = e

f ( vi )
( x ) = 8 ( −15 + 90 x 2
− 60 x + 8 x e
4 6
) − x2
, x
Puede comprobarse que

( vi ) ( vi )
max f ( x) =−f ( 0 ) = 120
−1 x 1

21
Mag Jose Luis Acuña Guillermo

22
 Mag Jose Luis Acuña Guillermo

Por lo tanto, una del error cometido en la

aproximación de la integral viene dada por

( )
2
120 1 x 5 x − 3
2 2
R( f )   dx
6! −1 4
1 2 1
= . =  0.047619
6 7 21
El valor buscado es
1
−1
− x2
e dx = 1.4936482....
f(
b 2n+ 2)
( )
n


b


a
w ( x ) f ( x ) dx =
i =0
ci f ( xi ) +
( 2n + 2 ) ! 
a
w ( x ) Pn+12dx

23
 Mag Jose Luis Acuña Guillermo

(Regla de Gauss-Legendre con 2 nodos). Si f es continua

en [-1,1], entonces

 −1   1 
-1 f ( x ) dx  G2 ( f ) = f  3  +  + E2 ( f )
1
f
 3
tiene grado de precisión n=3 y si fC4[-1,1] entonces

f ( ) (c)
4
E2 ( f ) =
135

para algún punto c[-1,1].

24
 Mag Jose Luis Acuña Guillermo

Ejemplo 2. Vamos a usar la cuadratura gaussiana

con 2 nodos para aproximar

dx
−1 x + 2 = Ln ( 2 ) − Ln (1)  1.09861
1

Solución dx
−1 x + 2  G2 ( f ) = f ( −0.57735) + f ( 0.57735)
1

1 1
= +
−0.57735 + 2 0.57735 + 2

= 0.70291 + 0.38800 = 1.09091,

Regla Trapecio
T ( f , 2 ) = f ( −1) + f (1) = 1 + 0.33333 = 1.33333,

1+ 2 +1 3
S ( f ,1) = f ( −1) + 4 f ( 0 ) + f (1) = = 1.11111.
Regla Simpson

3
25
 Mag Jose Luis Acuña Guillermo

La regla general de Gauss-Legendre (cuadratura

gaussiana) con N nodos es exacta para polinomios
de grado ≤ (2N-1) y su fórmula de cuadratura es

GN ( f ) = wN ,1 f ( xN ,1 ) + wN ,2 f ( xN ,2 ) + .... + wN ,N f ( xN ,N ) (**)
Los nodos xN,k y los pesos wN,k que hay usar están
tabulados en una tabla, así como la forma de los
términos del error EN(f); estos términos pueden usarse
para estimar la precisión del método de integración de
Gauss-Legendre. Los nodos, son de hecho, las raíces
de los polinomios de Legendre. Para la regla de
Gauss-Legendre con 3 nodos son
− ( 0.6 ) ( 0.6 )1/ 2
1/ 2
, 0,
y los pesos correspondientes son 5 / 9, 8 / 9 y 5 / 9.
26
 Mag Jose Luis Acuña Guillermo

1 N

 f ( x )dx =  wN , k (x N , k )+ E N ( f )
−1 k =1
N Nodos, x N , k Pesos, wN ,k Error, E N ( f )
2 0.5773502692 1.0000000000 f ( ) (c )
4
T 0.5773502692 1.0000000000
135
A 3 ±0.7745966692 0.5555555556 f ( ) (c )
6

0.0000000000 0.8888888888
15 750
B
4 ±0.8611363116 0.3478548451 f ( ) (c )
8

L ±0.3399810436 0.6521451549
3 472875
5 ±0.9061798459 0.2369268851 f ( ) (c )
10
A ±0.5384693101 0.4786286705
0.5688888888 1237 732 650
0.0000000000

27
Mag Jose Luis Acuña Guillermo

28
 Mag Jose Luis Acuña Guillermo

1 N

 f ( x ) dx =  wN , k (x N , k )+ E N ( f )
−1 k =1

T
N Nodos, x N , k Pesos, wN ,k Error, E N (f )
6 ±0.9324695142 0.1713244924 f ( ) (c ) 213 (6!)
12 4

±0.6612093865 0.3607615730
A (12!)
3
13!
±0.2386191861 0.4679139346
7 ±0.9491079123 0.1294849662 f ( ) (c ) 215 (7!)
14 4

B ±0.7415311856 0.2797053915
(14!)
3
±0.4058451514 0.3818300505 15!
0.4179591837
L 0.0000000000
(16 )
8 ±0.9602898565 0.1012285363 f (c ) 2 17
(8!)
4

±0.7966664774 0.2223810345
A (16! )
3
17!
±0.5255324009 0.3137066459
±0.1834346425 0.3626837834

29
29
 Mag Jose Luis Acuña Guillermo

(Regla de Gauss-Legendre con 3 nodos). Si f es

continua en [-1,1], entonces

( )
5 f − 3 / 5 + 8 f ( 0) + 5 f − 3 / 5 ( )+E (f)
-1 f ( x ) dx  G2 ( f ) =
1
3
9

y existe algún c[-1,1] tal que

f ( ) (c)
6
E3 ( f ) =
15750

30
 Mag Jose Luis Acuña Guillermo

Teorema (Traslación del método de Gauss-Legendre).

Supongamos que tenemos los nodos {xN,k} y los
pesos {wN,k} necesarios para aplicar la regla de
Gauss-Legendre con N nodos en [-1,1]. Entonces,
para aplicar el método de Gauss-legendre en un
intervalo [a,b], se puede usar el cambio de variable
a +b b-a b−a
t= + x con dt = dx (a)
2 2 2
y la relación
a+b b–a  b–a
a ( ) −1  2 + 2 x  2 dt
b 1
f t dt = f (b)
proporciona la fórmula de cuadratura
b–a N a+b b–a 
f ( t )dt = 
b
a 2 k =1
wN ,k f 
 2
+
2
x N ,k 
 (c)
31
 Mag Jose Luis Acuña Guillermo

Ejemplo Vamos a usar el método de Gauss-

Legendre con tres nodos para aproximar
dt
= Ln ( 5 ) –Ln (1)  1.609438
5
1 t
Solución En este caso a=1 y b=5, así que la fórmula
5 dt 5 –1 3  a +b b – a 
1 t
= G3 ( f ) = 
2 k =1
w3,k f 
 2
+
2
x3,k 

 
 


( 1/ 2
)
= 2 w3,1 f 3 − 2 ( 0.6 ) + w3,2 f ( 3 + 0 ) + w3,3 f 3 + 2 ( 0.6 ) (
1/ 2
) 

 5 8 5 
 9 9 9 
3.446359 + 2.666667 + 0.099096
= ( 2) = 1.602694
9
32
 Mag Jose Luis Acuña Guillermo

CUADRATURA DE GAUSS-LAGUERRE

Los pesos wi y los correspondientes ceros de algunos de los primeros

polinomios de Laguerre se presentan en la siguiente tabla:

 n

 f ( x )dx   wi f ( xi )
−x
e
0 i =1

33
Mag Jose Luis Acuña Guillermo

34
 GRACIAS

MAG JOSE LUIS ACUÑA GUILLERMO 35