Mna1 T05-1
Mna1 T05-1
Mna1 T05-1
Integración numérica
1 Introducción
4 Cuadratura de Gauss
Cuadratura de Gauss-Legendre
Cuadratura de Gauss-Chebyshev
Cuadratura de Gauss-Laguerre
Cuadratura de Gauss-Hermite
5 Integración múltiple
Nodos equiespaciados
Nodos no equiespaciados
Introducción
Aproximación de la integral
Longitud de arco Z 1 p
L= 1 + cos2 xdx
0
Función de Bessel Z π
1
Jn (x) = cos (x sin θ − nθ)dθ
x 0
Función error Z x
2 2
erf (x) = √ e−t dt
π 0
z 0 3 6 9 12 15 18 21 24 27 30
f (z) 190 141 104 77.5 54.7 42.5 31.5 23.3 17.3 12.8 9.5
El cálculo de la fuerza total, es decir, de una integral definida cuyo integrando no posee
una expresión analı́tica, es imprescindible para un correcto diseño del mástil.
¿Qué podemos hacer si nuestra función no tiene una expresión explı́cita y solo conocemos
la función en algunos puntos?
Métodos Numéricos Aplicados I PER 3362 5 / 63
Introducción
INTEGRACIÓN NUMÉRICA
N O D O S E Q U I E S PA C I A D O S N O D O S N O E Q U I E S PA C I A D O S
Newton-Cotes Cuadratura de Gauss
I N T E G R A C I Ó N S I M P L E Y M Ú LT I P L E
Objetivos
○ Conocer la expresión general de las técnicas de cuadratura
○ Comprender e implementar las fórmulas derivadas de las expresiones de
Newton-Cotes, cerradas y abiertas
○ Comprender e implementar los casos particulares de la cuadratura de Gauss
○ Aplicar la integración numérica sobre integración múltiple
Métodos Numéricos Aplicados I PER 3362 6 / 63
Introducción
Cuadratura numérica
La técnica de cuadratura numérica consiste en calcular una integral a partir de unos
valores discretos, obtenidos del valor del integrando en los diferentes nodos o en valores
tabulados.
En general, la cuadratura numérica consiste en obtener
Z b Xn
I= f (x) dx ≈ ai f (xi ).
a i=0
Cuadratura numérica
La técnica de cuadratura numérica consiste en calcular una integral a partir de unos
valores discretos, obtenidos del valor del integrando en los diferentes nodos o en valores
tabulados.
En general, la cuadratura numérica consiste en obtener
Z b Xn
I= f (x) dx ≈ ai f (xi ).
a i=0
0.6 0.6
0.4 0.4
0.2 0.2
0 0
0 0.5 1 1.5 2 0 0.5 1 1.5 2
x x
Métodos Numéricos Aplicados I PER 3362 7 / 63
Introducción
Cuadratura numérica
La técnica de cuadratura numérica consiste en calcular una integral a partir de unos
valores discretos, obtenidos del valor del integrando en los diferentes nodos o en valores
tabulados.
En general, la cuadratura numérica consiste en obtener
Z b Xn
I= f (x) dx ≈ ai f (xi ).
a i=0
f (n+1) (ξ(x)) Qn
Error cometido en f (x) ≈ ln (x): (x) = i=0 (x − xi ).
(n + 1)!
Por tanto: n
X Z b
I≈ f (xi ) Li (x) dx,
i=0 a
Por tanto: n
X Z b
I≈ f (xi ) Li (x) dx,
i=0 a
Xn Z b
𝑥𝑖
I ≈𝑥 f (xi ) Li (x) dx,
𝑖+1 a
i=0
𝑎 𝑥 b−a
Nodos equiespaciados: x ∈ [a, b], h = n
, {x0 = a, x1 , x2 , . . . , xn = b}
𝑎 𝑎+ℎ 𝑎 + 2ℎ 𝑏−ℎ 𝑏
Simpson 3/8:
3h5 (iv)
3h 2a + b a + 2b
I≈ f (a) + 3f + 3f + f (b) , =− f (ξ)
8 3 3 80
Milne:
2h 3a + b a+b a + 3b
I≈ 7f (a) + 32f + 12f + 32f + 7f (b) ,
45 4 2 4
8h7 (vi)
=− f (ξ)
945
Métodos Numéricos Aplicados I PER 3362 11 / 63
Contenidos
1 Introducción
4 Cuadratura de Gauss
5 Integración múltiple
Desarrollo
Puntos (n = 1): x0 = a, x1 = b
Z b Z b
I ≈ f (a) L0 (x) dx + f (b) L1 (x) dx
a a
Z b Z b
x−b x−a
= f (a) dx + f (b) dx
a a−b b−a a
f (a) b f (b) b
Z Z
= − (x − b) dx + (x − a) dx
h a h a
h
= (f (a) + f (b)) .
2
Método de Trapecios
h
I≈ (f (a) + f (b))
2
3
Error: = − h12 f 00 (ξ(x))
Métodos Numéricos Aplicados I PER 3362 13 / 63
Fórmulas cerradas de Newton-Cotes
Método de Trapecios
0.4 0.4
0.2 0.2
0 0
0 1 2 0 1 2
x x
n−1
hX h
I≈ (f (xi ) + f (xi+1 )) = (f (x0 ) + 2f (x1 ) + 2f (x2 ) + · · · + 2f (xn−1 ) + f (xn ))
2 i=0 2
2
Error: = − h12 (b − a)f 00 (ξ(x))
Métodos Numéricos Aplicados I PER 3362 14 / 63
Fórmulas cerradas de Newton-Cotes
Método de Trapecios
Trapecios.m
function I=Trapecios(f,a,b,n)
% fórmula de Trapecios compuesta
h=(b-a)/n;
x=a:h:b;
pesos=[1 2*ones(1,n-1) 1];
I=h/2*sum(pesos.*f(x));
end
1 Introducción
4 Cuadratura de Gauss
5 Integración múltiple
Desarrollo
a+b
Puntos (n = 2): x0 = a, x1 = 2
, x2 = b
Z b Z b Z b
a+b
I ≈ f (a) L0 (x) dx + f L1 (x) dx + f (b) L2 (x) dx
a 2 a a
h a+b
= f (a) + 4f + f (b)
3 2
Método de Simpson
h a+b
I≈ f (a) + 4f + f (b)
3 2
h5 (iv)
Error: = − f (ξ(x))
90
h4
Error: = − (b − a)f (iv) (ξ(x))
180
Simpson.m
function I=Simpson(f,a,b,n)
% fórmula de Simpson compuesta
h=(b-a)/n; x=a:h:b;
pesos=ones(1,n+1); pesos(2:2:n)=4; pesos(3:2:n-1)=2;
I=h/3*sum(pesos.*f(x));
end
Métodos Numéricos Aplicados I PER 3362 18 / 63
Fórmulas cerradas de Newton-Cotes
Método de Simpson
Simpson.m
function I=Simpson(f,a,b,n)
% fórmula de Simpson compuesta
h=(b-a)/n;
x=a:h:b;
pesos=ones(1,n+1);
pesos(2:2:n)=4; pesos(3:2:n-1)=2;
I=h/3*sum(pesos.*f(x));
end
Rπ
Ejemplo 2. Calcula I = 02 sin(x)e−x dx, utilizando los métodos de Trapecios y
Simpson tomando 4 y 8 subintervalos. Para cada caso, calcula el error cometido
−π
sabiendo que el resultado analı́tico es 1−e2 2 .
Trapecios Simpson T S
n=4 0.380591 0.395839 0.015469 0.000221
n=8 0.392183 0.396047 0.003877 0.000013
n
X Z b
I≈ f (xi ) Li (x) dx
i=0 a
𝑎 𝑎+ℎ 𝑎 + 2ℎ 𝑏−ℎ 𝑏
1 Introducción
4 Cuadratura de Gauss
5 Integración múltiple
Desarrollo
a+b
Puntos (n = 0): x−1 = a, x0 = 2
, x1 = b
Z b
a+b
I ≈ f (x0 ) dx = (b − a)f (x0 ) = (b − a)f
a 2
h3 00
Error: = − f (ξ(x))
3
x−1 ←- x0 ,→ x1
n/2
X
I ≈ (x2i+1 − x2i−1 )f (x2i )
i=0
= 2h (f (x0 ) + f (x2 ) + · · · + f (xn−2 ) + f (xn ))
b−a 2 00
Error: = 6
h f (ξ)
Métodos Numéricos Aplicados I PER 3362 27 / 63
Fórmulas abiertas de Newton-Cotes
Método de Punto Medio
PuntoMedio.m
function I=PuntoMedio(f,a,b,n)
% fórmula de punto medio compuesta
h=(b-a)/(n+2);
x=a+h:h:b-h;
I=2*h*sum(f(x(1:2:end)));
end
𝐹
x 0.00 1.52 3.04 4.56 6.08 7.60 9.12
F (x) 0.00 40.04 57.83 62.28 46.71 53.38 22.24 𝛼
𝐹
x 0.00 1.52 3.04 4.56 6.08 7.60 9.12
F (x) 0.00 40.04 57.83 62.28 46.71 53.38 22.24 𝛼
40
20
1.52
xi = i, i = −1, 0, . . . , 13, 10
2 0
-1.52 0 1.52 3.04 4.56 6.08 7.6 9.12 10.64
x
𝐹
x 0.00 1.52 3.04 4.56 6.08 7.60 9.12
F (x) 0.00 40.04 57.83 62.28 46.71 53.38 22.24 𝛼
40
20
1.52
xi = i, i = −1, 0, . . . , 13, 10
2 0
-1.52 0 1.52 3.04 4.56 6.08 7.6 9.12 10.64
x
𝐹
x 0.00 1.52 3.04 4.56 6.08 7.60 9.12
F (x) 0.00 40.04 57.83 62.28 46.71 53.38 22.24 𝛼
40
20
1.52
xi = i, i = −1, 0, . . . , 13, 10
2 0
-1.52 0 1.52 3.04 4.56 6.08 7.6 9.12 10.64
x
Cuadratura de Gauss
Objetivo
Obtención de los nodos x1 , x2 , · · · , xn y de los coeficientes c1 , c2 , . . . , cn que minimizan
el error obtenido en la aproximación:
Z b Xn
w(x)f (x) dx ≈ ci f (xi ),
a i=1
hpi , pj i = 0, i 6= j.
Fórmulas de cuadratura
w(x) = 1, [a, b] = [−1, 1] ú Cuadratura de Gauss-Legendre
w(x) = √ 1 2 , [a, b] = [−1, 1] ú Cuadratura de Gauss-Chebyshev
1−x
Error
b
f (2n) (ξ(x))
Z
= p2n (x)w(x) dx, ξ ∈ (a, b).
(2n)! a
1 Introducción
4 Cuadratura de Gauss
Cuadratura de Gauss-Legendre
Cuadratura de Gauss-Chebyshev
Cuadratura de Gauss-Laguerre
Cuadratura de Gauss-Hermite
5 Integración múltiple
Z 1 n
X
f (x) dx ≈ ci f (xi )
−1 i=1
ú w(x) = 1
ú [a, b] = [−1, 1]
ú Polinomios pi (x):
p0 (x) = 1,
p1 (x) = x,
1
pk+1 (x) = k+1
[(2k + 1)xpk (x) − kpk−1 (x)] , k = 1, 2, . . . , n − 1.
ú Raı́ces de pn (x):
1 1 4k − 1
xi = 1 − 2 + 3 cos π
8n 8n 4n + 2
ú Coeficientes ci :
2
ci =
(1 − x2i )(p0n (xi ))2
i 1 2 3 4 5
xi -0.577350 0.577350 - - -
n=2
ci 1.000000 1.000000 - - -
xi 0.000000 -0.774597 0.774597 - -
n=3
ci 0.888889 0.555556 0.555556 - -
xi -0.339981 -0.861136 0.339981 0.861136 -
n=4
ci 0.652145 0.347855 0.652145 0.347855 -
xi 0.000000 -0.538469 -0.906180 0.538469 0.906180
n=5
ci 0.568889 0.478629 0.236927 0.478629 0.236927
Tabla: Nodos y coeficientes de la cuadratura de Gauss-Legendre
= 0.109364.
Métodos Numéricos Aplicados I PER 3362 37 / 63
Contenidos
1 Introducción
4 Cuadratura de Gauss
Cuadratura de Gauss-Legendre
Cuadratura de Gauss-Chebyshev
Cuadratura de Gauss-Laguerre
Cuadratura de Gauss-Hermite
5 Integración múltiple
1
ú w(x) = √
1−x2
ú [a, b] = [−1, 1]
ú Polinomios pi (x):
p0 (x) = 1,
p1 (x) = x,
pk+1 (x) = 2xpk−1 (x) − pk−2 , k = 2, 3, . . . , n.
ú Raı́ces de pn (x):
2k − 1
xi = cos π
2n
ú Coeficientes ci :
π
ci =
n
Aproximamos la integral:
Z 1 n n
f (x) X πX
√ dx ≈ ci f (xi ) = f (xi )
−1 1 − x2 i=1
n i=1
i 1 2 3 4 5
n=2 xi -0.707107 0.707107 - - -
n=3 xi -0.866025 0.000000 0.866025 - -
n=4 xi -0.923880 -0.382683 0.382683 0.923880 -
n=5 xi -0.951057 -0.587785 0.000000 0.587785 0.951057
Tabla: Nodos de la cuadratura de Gauss-Chebyshev
π −0.951057
+ e−0.587785 + e0 + e0.587785 + e0.951057 = 3.977463
I≈ e
5
Métodos Numéricos Aplicados I PER 3362 41 / 63
Contenidos
1 Introducción
4 Cuadratura de Gauss
Cuadratura de Gauss-Legendre
Cuadratura de Gauss-Chebyshev
Cuadratura de Gauss-Laguerre
Cuadratura de Gauss-Hermite
5 Integración múltiple
Z +∞ n
X
e−x f (x) dx ≈ ci f (xi )
0 i=1
ú w(x) = e−x
ú [a, b] = [0, +∞)
ú Polinomios pi (x):
p0 (x) = 1,
p1 (x) = 1 − x,
pk+2 (x) = (2k + 3 − x)pk+1 (x) − (k + 1)2 pk (x), k = 0, 1, . . . , n − 2.
ú Raı́ces de pn (x):
2 2
j0i −2 + j0k
xi = 1+ ,
4hn 48h2n
siendo hn = n + 21 y j0i la raı́z i-ésima de la función J0 (x), es decir, la función de
Bessel de primera especie y orden.
ú Coeficientes ci :
(n!)2 xi
ci = 2 .
pn+1 (xi )
Métodos Numéricos Aplicados I PER 3362 43 / 63
Cuadratura de Gauss-Laguerre
R +∞
Ejemplo 8. Calcula la integral I = 0
e−10x sin(x) dx utilizando la cuadratura
de Gauss-Laguerre para n = 3
R +∞
El dominio de la integral está en [0, +∞), pero no es de la forma I = 0
e−x f (x) dx
ú Cambio de variable:
y = 10x ⇒ dy = 10 dx
entonces Z +∞ Z +∞
1 y
I= e−10x sin(x) dx = e−y sin dy.
0 10 0 10
Polinomios de Laguerre hasta n = 4:
p0 (x) = 1,
p1 (x) = 1 − x,
p2 (x) = x2 − 4x + 2,
p3 (x) = −x3 + 9x2 − 18x + 6,
p4 (x) = x4 − 16x3 + 72x2 − 96x + 24,
1 Introducción
4 Cuadratura de Gauss
Cuadratura de Gauss-Legendre
Cuadratura de Gauss-Chebyshev
Cuadratura de Gauss-Laguerre
Cuadratura de Gauss-Hermite
5 Integración múltiple
Z +∞ n
2 X
e−x f (x) dx ≈ ci f (xi )
−∞ i=1
2
ú w(x) = e−x
ú [a, b] = (−∞, +∞)
ú Polinomios pi (x):
p0 (x) = 1,
p1 (x) = 2x,
pk+2 (x) = 2xpk+1 (x) − 2(k + 1)pk (x), k = 0, 1, . . . , n − 2.
ú Coeficientes ci : √
2n−1 n! π
ci = 2
.
n2 Hn−1 (xi )
R +∞ 2
Ejemplo 9. Calcula la integral I = −∞
e−4x |x| dx utilizando la cuadratura de
Gauss-Hermite para n = 4
Cambio de variable:
y = 2x ⇒ dy = 2 dx
entonces Z +∞ Z +∞
2 1 2
I= e−4x |x| dx = e−y |y| dy.
−∞ 4 −∞
p0 (x) = 1,
p1 (x) = 2x,
p2 (x) = 4x2 − 2,
p3 (x) = 8x3 − 12x,
p4 (x) = 16x4 − 48x2 + 12.
R +∞ 2
Ejemplo 9. Calcula la integral I = −∞
e−4x |x| dx utilizando la cuadratura de
Gauss-Hermite para n = 4
Calculamos las raı́ces de p4 (x):
>> xi=roots([16 0 -48 0 12])
xi =
-1.650680123885785
1.650680123885786
-0.524647623275290
0.524647623275290
Obtenemos los coeficientes:
>> syms x
>> p3=8*x.∧3-12*x;
>> ci=2∧(4-1)*factorial(4)*sqrt(pi)/4∧2./(double(subs(p3,x,xi)).∧2);
El valor de la integral para n = 4 será:
f=@(y) abs(y);
>> I=1/4*sum(ci.*f(xi))
I =
0.278258872775874
Integración múltiple
1
0.9
0.8
0.5
2
f : R → R, (x, y) ∈ [a, b] × [c, d] 0.7
0
2
0.5 0.6
0
0
y -2 -0.5
x
Objetivo
Aproximación de integrales dobles utilizando:
Nodos equiespaciados (método de Trapecios)
Nodos no equiespaciados (fórmulas de cuadratura de Gauss)
1 Introducción
4 Cuadratura de Gauss
5 Integración múltiple
Nodos equiespaciados
Nodos no equiespaciados
Rd
Si nombramos g(x) = c
f (x, y) dy, entonces
Z b
I= g(x) dx.
a
d−c
b−a
Con la notación h = , k=
, aproximamos la integral como
n m
Z b n−1 n−1 Z d Z d
hX hX
I= g(x) dx ≈ g(xi ) + g(xi+1 ) = f (xi , y) dy + f (xi+1 , y) dy
a 2 i=0 2 i=0 c c
n−1
( m−1 m−1
)
hX k X k X
≈ f (xi , yj ) + f (xi , yj+1 ) + f (xi , yj ) + f (xi , yj+1 ) .
2 i=0 2 j=0 2 j=0
b−a d−c
siendo h = n
, k= m
.
𝑦𝑚 1 2 2 2 1
𝑦𝑚−1 2 4 4 4 2
𝑦2 2 4 4 4 2
2 4 4 4 2
𝑦1
1 2 2 2 1
𝑦0
𝑥0 𝑥1 𝑥2 𝑥𝑛−1 𝑥𝑛
1 Introducción
4 Cuadratura de Gauss
5 Integración múltiple
Nodos equiespaciados
Nodos no equiespaciados
Enunciado
Calcula la integral
s
x ∂f
2
∂f
2
I= (x, y) + (x, y) dA,
∂x ∂y
R
x2 + y 2 + f 2 (x, y) = 9,
y el recinto de integración es
R = {(x, y) ∈ R2 : 0 ≤ x ≤ 1, 0 ≤ y ≤ 1}.
Solución
p
La función f que aparece en el integrando es f (x, y) = 9 − x2 − y 2 , y sus derivadas
parciales
∂f x ∂f y
(x, y) = − p , (x, y) = − p ,
∂x 9−x −y 2 2 ∂y 9 − x2 − y 2
de modo que el integrando es
s 2 2 s
∂f ∂f x2 + y 2
(x, y) + (x, y) = .
∂x ∂y 9 − x2 − y 2
El recinto de integración es un cuadrado [0, 1] × [0, 1], por lo que la integral a resolver es
Z 1Z 1s
x2 + y 2
I= dy dx.
0 0 9 − x2 − y 2
El resultado es
I ≈ 0.267814255559730.
u = 2x − 1 ⇒ du = 2 dx, v = 2y − 1 ⇒ dv = 2 dy.
Por tanto,
s r
Z 1 Z 1 Z 1 Z 1
x2 + y 2 1 u2 + 2u + v 2 + 2v + 2
I= dy dx = dv du
0 0 9 − x2 − y 2 4 −1 −1 34 − u2 − 2u − v 2 − 2v
u1 = v1 = −0.339981, c1 = 0.652145;
u2 = v2 = −0.861136, c2 = 0.347855;
u3 = v3 = −u1 , c3 = c1 ;
u4 = v4 = −u2 , c4 = c2 .
De este modo,
s
4 4
1 XX u2i + 2ui + vj2 + 2vj + 2
I≈ ci cj ≈ 0.267770529696778.
4 i=1 j=1 34 − u2i − 2ui − vj2 − 2vj
Métodos Numéricos Aplicados I PER 3362 61 / 63
Para finalizar...
i Lecciones magistrales
L Material complementario: A fondo
Bibliografı́a recomendada
TEST DE APRENDIZAJE!!