métodos numéricos - juanrodriguezprieto.files.wordpress.com · ejemplo de la regla del trapecio f...
TRANSCRIPT
![Page 1: Métodos Numéricos - juanrodriguezprieto.files.wordpress.com · Ejemplo de la regla del trapecio f ( x ) = 0.2 + 25 x - 200 x 2 + 675 x 3 - 900 x 4 + 400 x 5 Aproxime la integral](https://reader030.vdocumento.com/reader030/viewer/2022021522/5c5d427309d3f2d72f8cb1b7/html5/thumbnails/1.jpg)
Métodos Numéricos
Juan Manuel Rodríguez Prieto
I.M., M.Sc., Ph.D.
![Page 2: Métodos Numéricos - juanrodriguezprieto.files.wordpress.com · Ejemplo de la regla del trapecio f ( x ) = 0.2 + 25 x - 200 x 2 + 675 x 3 - 900 x 4 + 400 x 5 Aproxime la integral](https://reader030.vdocumento.com/reader030/viewer/2022021522/5c5d427309d3f2d72f8cb1b7/html5/thumbnails/2.jpg)
Integración numérica
![Page 3: Métodos Numéricos - juanrodriguezprieto.files.wordpress.com · Ejemplo de la regla del trapecio f ( x ) = 0.2 + 25 x - 200 x 2 + 675 x 3 - 900 x 4 + 400 x 5 Aproxime la integral](https://reader030.vdocumento.com/reader030/viewer/2022021522/5c5d427309d3f2d72f8cb1b7/html5/thumbnails/3.jpg)
Integración numérica
• Objetivo: aproximar el valor de la integral
• Limitaciones de la integración analítica
• Las expresiones analíticas de son
desconocidas
• tiene una integral analítica complicada o
desconocida
I = f (x)dxa
b
ò
f (x)
f (x)
![Page 4: Métodos Numéricos - juanrodriguezprieto.files.wordpress.com · Ejemplo de la regla del trapecio f ( x ) = 0.2 + 25 x - 200 x 2 + 675 x 3 - 900 x 4 + 400 x 5 Aproxime la integral](https://reader030.vdocumento.com/reader030/viewer/2022021522/5c5d427309d3f2d72f8cb1b7/html5/thumbnails/4.jpg)
Integración numérica
de Newton-Cotes
• Métodos que remplazan una función complicada
o datos tabulados por un polinomio de
aproximación que es fácil de integral
• Donde es un polinomio de la forma
I = f (x)dxa
b
ò @ fn(x)dxa
b
ò
fn(x)
![Page 5: Métodos Numéricos - juanrodriguezprieto.files.wordpress.com · Ejemplo de la regla del trapecio f ( x ) = 0.2 + 25 x - 200 x 2 + 675 x 3 - 900 x 4 + 400 x 5 Aproxime la integral](https://reader030.vdocumento.com/reader030/viewer/2022021522/5c5d427309d3f2d72f8cb1b7/html5/thumbnails/5.jpg)
Integración numérica
de Newton-Cotes
Aproximación de una integral mediante el área
bajo una línea recta.
![Page 6: Métodos Numéricos - juanrodriguezprieto.files.wordpress.com · Ejemplo de la regla del trapecio f ( x ) = 0.2 + 25 x - 200 x 2 + 675 x 3 - 900 x 4 + 400 x 5 Aproxime la integral](https://reader030.vdocumento.com/reader030/viewer/2022021522/5c5d427309d3f2d72f8cb1b7/html5/thumbnails/6.jpg)
Integración numérica
de Newton-Cotes
Aproximación de una integral mediante el área
bajo una parábola.
![Page 7: Métodos Numéricos - juanrodriguezprieto.files.wordpress.com · Ejemplo de la regla del trapecio f ( x ) = 0.2 + 25 x - 200 x 2 + 675 x 3 - 900 x 4 + 400 x 5 Aproxime la integral](https://reader030.vdocumento.com/reader030/viewer/2022021522/5c5d427309d3f2d72f8cb1b7/html5/thumbnails/7.jpg)
La regla del trapecio
La regla del trapecio es la primera de las formulas de integración de Newton-
Cotes. Corresponde al caso donde el polinomio es de primer grado (línea
recta).
La recta que pasa por los puntos y esta dada por:
El área bajo la línea recta f1(x) es una aproximación de la integral de f(x) entre
los límites a y b
I = f (x)dxa
b
ò @ a0 + a1xdxa
b
ò(a, f (a)) (b, f (b))
a0 + a1x = f (a)+f (b)- f (a)
b- a(x- a)
I @ f (a)+f (b)- f (a)
b- a(x- a)dx
a
b
ò = (b- a)f (a)+ f (b)
2
![Page 8: Métodos Numéricos - juanrodriguezprieto.files.wordpress.com · Ejemplo de la regla del trapecio f ( x ) = 0.2 + 25 x - 200 x 2 + 675 x 3 - 900 x 4 + 400 x 5 Aproxime la integral](https://reader030.vdocumento.com/reader030/viewer/2022021522/5c5d427309d3f2d72f8cb1b7/html5/thumbnails/8.jpg)
La regla del trapecio
I = f (x)dxa
b
ò @ (b- a)f (a)+ f (b)
2
Geométricamente, la regla del trapecio equivale a aproximar el área
bajo la curva f(x), como el área del trapecio que se forma al unir los
puntos (a, f (a)) (b, f (b))
![Page 9: Métodos Numéricos - juanrodriguezprieto.files.wordpress.com · Ejemplo de la regla del trapecio f ( x ) = 0.2 + 25 x - 200 x 2 + 675 x 3 - 900 x 4 + 400 x 5 Aproxime la integral](https://reader030.vdocumento.com/reader030/viewer/2022021522/5c5d427309d3f2d72f8cb1b7/html5/thumbnails/9.jpg)
Error de la regla del trapecio
Et = -1
12f ''(x)(b- a)3
Cuando usamos la integral bajo un segmento de línea recta para
aproximar la integral bajo una curva, se tiene un error que puede ser
importante. Una estimación al error de truncamiento para una sola
aplicación de la regla del trapecio es
Donde esta en algún lugar del intervalo de a a bx
![Page 10: Métodos Numéricos - juanrodriguezprieto.files.wordpress.com · Ejemplo de la regla del trapecio f ( x ) = 0.2 + 25 x - 200 x 2 + 675 x 3 - 900 x 4 + 400 x 5 Aproxime la integral](https://reader030.vdocumento.com/reader030/viewer/2022021522/5c5d427309d3f2d72f8cb1b7/html5/thumbnails/10.jpg)
Ejemplo de la regla del trapecio
f (x) = 0.2 + 25x- 200x2 + 675x3 - 900x4 + 400x5
Aproxime la integral de la curva
Entre a=0 y b=0.8 usando la regla del trapecio
f (a = 0) = 0.2
f (b = 0.8) = 0.232
b- a= 0.8
I @ (b- a)f (a)+ f (b)
2= 0.8
0.2 + 0.232
2= 0.1728
![Page 11: Métodos Numéricos - juanrodriguezprieto.files.wordpress.com · Ejemplo de la regla del trapecio f ( x ) = 0.2 + 25 x - 200 x 2 + 675 x 3 - 900 x 4 + 400 x 5 Aproxime la integral](https://reader030.vdocumento.com/reader030/viewer/2022021522/5c5d427309d3f2d72f8cb1b7/html5/thumbnails/11.jpg)
Ejemplo de la regla del trapecio
f (x) = 0.2 + 25x- 200x2 + 675x3 - 900x4 + 400x5
Aproxime la integral de la curva
![Page 12: Métodos Numéricos - juanrodriguezprieto.files.wordpress.com · Ejemplo de la regla del trapecio f ( x ) = 0.2 + 25 x - 200 x 2 + 675 x 3 - 900 x 4 + 400 x 5 Aproxime la integral](https://reader030.vdocumento.com/reader030/viewer/2022021522/5c5d427309d3f2d72f8cb1b7/html5/thumbnails/12.jpg)
Ejemplo de la regla del trapecio
f (x) = 0.2 + 25x- 200x2 + 675x3 - 900x4 + 400x5
Aproxime la integral de la curva
El área bajo la curva es mucho mayor que el área debajo de la
aproximación lineal, de acuerdo a la gráfica mostrada en la
diapositiva anterior.
¿Cuanto es el error debido al aproximar la integral del polinomio
de grado usando la regla del trapecio?
I = f (x)dx0
0.8
ò = (0.2 + 25x- 200x2 + 675x3 - 900x4 + 400x5 )dx0
0.8
ò
= 1.640533
![Page 13: Métodos Numéricos - juanrodriguezprieto.files.wordpress.com · Ejemplo de la regla del trapecio f ( x ) = 0.2 + 25 x - 200 x 2 + 675 x 3 - 900 x 4 + 400 x 5 Aproxime la integral](https://reader030.vdocumento.com/reader030/viewer/2022021522/5c5d427309d3f2d72f8cb1b7/html5/thumbnails/13.jpg)
Ejemplo de la regla del trapecio
El error porcentual es de
E(%) =1.640533 - 0.1728
1.640533*100 = 89.5%
¿Cómo se puede disminuir el error?
Se puede dividir el intervalo de interés, en más
intervalos. En otras palabras aplicar varias veces
la regla del trapecio en el intervalo de interés
![Page 14: Métodos Numéricos - juanrodriguezprieto.files.wordpress.com · Ejemplo de la regla del trapecio f ( x ) = 0.2 + 25 x - 200 x 2 + 675 x 3 - 900 x 4 + 400 x 5 Aproxime la integral](https://reader030.vdocumento.com/reader030/viewer/2022021522/5c5d427309d3f2d72f8cb1b7/html5/thumbnails/14.jpg)
La regla del trapecio de aplicación
múltipleUna forma de mejorar la precisión
de la regla del trapecio consiste en
dividir el intervalo de integración
de a a b en varios segmentos, y
aplicar el método a cada uno de
ellos. La área asociada a cada uno
de los intervalos se suman después
para obtener la integral en todo el
intervalo. Las ecuaciones
resultantes se llaman fórmulas de
integración, de aplicación múltiple
o compuesta.
Vamos a dividir el intervalo de
interés en n segmentos del mismo
ancho, es decir tendremos n+1
puntos igualmente espaciados.
![Page 15: Métodos Numéricos - juanrodriguezprieto.files.wordpress.com · Ejemplo de la regla del trapecio f ( x ) = 0.2 + 25 x - 200 x 2 + 675 x 3 - 900 x 4 + 400 x 5 Aproxime la integral](https://reader030.vdocumento.com/reader030/viewer/2022021522/5c5d427309d3f2d72f8cb1b7/html5/thumbnails/15.jpg)
La regla del trapecio de aplicación
múltiple
Aplicando la regla del trapecio a cada una de las integrales
![Page 16: Métodos Numéricos - juanrodriguezprieto.files.wordpress.com · Ejemplo de la regla del trapecio f ( x ) = 0.2 + 25 x - 200 x 2 + 675 x 3 - 900 x 4 + 400 x 5 Aproxime la integral](https://reader030.vdocumento.com/reader030/viewer/2022021522/5c5d427309d3f2d72f8cb1b7/html5/thumbnails/16.jpg)
La regla del trapecio de aplicación
múltiple
Simplificando, se obtiene
I = f (x)dxa
b
ò @h
2f (a)+ 2 f (xi )
i=1
n-1
å + f (b)é
ëê
ù
ûú
I = f (x)dxa
b
ò @ (b- a)
f (a) + 2 f (xi )i=1
n-1
å + f (b)
2n
é
ë
êêêê
ù
û
úúúú
![Page 17: Métodos Numéricos - juanrodriguezprieto.files.wordpress.com · Ejemplo de la regla del trapecio f ( x ) = 0.2 + 25 x - 200 x 2 + 675 x 3 - 900 x 4 + 400 x 5 Aproxime la integral](https://reader030.vdocumento.com/reader030/viewer/2022021522/5c5d427309d3f2d72f8cb1b7/html5/thumbnails/17.jpg)
Erro de la regla del trapecio de
aplicación múltiple
Et = -(b- a)3
12n3f ''(xi )
i=1
n-1
å
Simplificando, se obtiene
Ea =(b- a)3
12n2
f ''(xi )i=1
n-1
å
n
=(b- a)3
12n2f ''
Así, si se duplica el número de segmentos, el error de truncamiento se divide
entre 4
![Page 18: Métodos Numéricos - juanrodriguezprieto.files.wordpress.com · Ejemplo de la regla del trapecio f ( x ) = 0.2 + 25 x - 200 x 2 + 675 x 3 - 900 x 4 + 400 x 5 Aproxime la integral](https://reader030.vdocumento.com/reader030/viewer/2022021522/5c5d427309d3f2d72f8cb1b7/html5/thumbnails/18.jpg)
Ejemplo de la regla del trapecio
f (x) = 0.2 + 25x- 200x2 + 675x3 - 900x4 + 400x5
Aproxime la integral de la curva
Entre a=0 y b=0.8 usando la regla del trapecio de aplicación múltiple
Usaremos inicialmente n=2, lo cual da un h=0.4f (a = 0) = 0.2
f (x1 = 0.4) = 2.456
f (b = 0.8) = 0.232b- a
n= 0.4
I = f (x)dxa
b
ò @h
2f (a)+ 2 f (x1)+ f (b)[ ]
=0.4
20.2 + 2 *2.456 + 0.232[ ]
= 1.0688
![Page 19: Métodos Numéricos - juanrodriguezprieto.files.wordpress.com · Ejemplo de la regla del trapecio f ( x ) = 0.2 + 25 x - 200 x 2 + 675 x 3 - 900 x 4 + 400 x 5 Aproxime la integral](https://reader030.vdocumento.com/reader030/viewer/2022021522/5c5d427309d3f2d72f8cb1b7/html5/thumbnails/19.jpg)
Ejemplo de la regla del trapecio
f (x) = 0.2 + 25x- 200x2 + 675x3 - 900x4 + 400x5
Aproxime la integral de la curva
Entre a=0 y b=0.8 usando la regla del trapecio de aplicación múltiple
Usaremos inicialmente n=2, lo cual da un h=0.4
![Page 20: Métodos Numéricos - juanrodriguezprieto.files.wordpress.com · Ejemplo de la regla del trapecio f ( x ) = 0.2 + 25 x - 200 x 2 + 675 x 3 - 900 x 4 + 400 x 5 Aproxime la integral](https://reader030.vdocumento.com/reader030/viewer/2022021522/5c5d427309d3f2d72f8cb1b7/html5/thumbnails/20.jpg)
Ejemplo de la regla del trapecio
f (x) = 0.2 + 25x- 200x2 + 675x3 - 900x4 + 400x5
Aproxime la integral de la curva
Entre a=0 y b=0.8 usando la regla del trapecio de aplicación múltiple
Usaremos inicialmente n=2, lo cual da un h=0.4
Utilizando dos divisiones del intervalo el error ha disminuido de un
89.5% a un 34.9%
E(%) =1.640533 -1.0688
1.640533*100 = 34.9%
![Page 21: Métodos Numéricos - juanrodriguezprieto.files.wordpress.com · Ejemplo de la regla del trapecio f ( x ) = 0.2 + 25 x - 200 x 2 + 675 x 3 - 900 x 4 + 400 x 5 Aproxime la integral](https://reader030.vdocumento.com/reader030/viewer/2022021522/5c5d427309d3f2d72f8cb1b7/html5/thumbnails/21.jpg)
Ejemplo de la regla del trapecio
f (x) = 0.2 + 25x- 200x2 + 675x3 - 900x4 + 400x5
Aproxime la integral de la curva
Entre a=0 y b=0.8 usando la regla del trapecio de aplicación múltiple
Usaremos n=3, lo cual da un h=0.2667
f (a = 0) = 0.2
f (x1 = 0.2667) = 1.4327
f (x2 = 0.5333) = 3.4872
f (b = 0.8) = 0.232
b- a
n= 0.2667
I = f (x)dxa
b
ò @h
2f (a)+ 2 f (x1) + 2 f (x2 )+ f (b)[ ]
=0.4
20.2 + 2 *1.4327 + 2 * 3.4872 + 0.232[ ]
= 1.3695
![Page 22: Métodos Numéricos - juanrodriguezprieto.files.wordpress.com · Ejemplo de la regla del trapecio f ( x ) = 0.2 + 25 x - 200 x 2 + 675 x 3 - 900 x 4 + 400 x 5 Aproxime la integral](https://reader030.vdocumento.com/reader030/viewer/2022021522/5c5d427309d3f2d72f8cb1b7/html5/thumbnails/22.jpg)
Ejemplo de la regla del trapecio
f (x) = 0.2 + 25x- 200x2 + 675x3 - 900x4 + 400x5
Aproxime la integral de la curva
Entre a=0 y b=0.8 usando la regla del trapecio de aplicación múltiple
Usaremos inicialmente n=3, lo cual da un h=0.2667
![Page 23: Métodos Numéricos - juanrodriguezprieto.files.wordpress.com · Ejemplo de la regla del trapecio f ( x ) = 0.2 + 25 x - 200 x 2 + 675 x 3 - 900 x 4 + 400 x 5 Aproxime la integral](https://reader030.vdocumento.com/reader030/viewer/2022021522/5c5d427309d3f2d72f8cb1b7/html5/thumbnails/23.jpg)
Ejemplo de la regla del trapecio
f (x) = 0.2 + 25x- 200x2 + 675x3 - 900x4 + 400x5
Aproxime la integral de la curva
Entre a=0 y b=0.8 usando la regla del trapecio de aplicación múltiple
Usaremos inicialmente n=3, lo cual da un h=0.2667
Utilizando tres divisiones del intervalo el error ha disminuido de un
89.5% a un 16.5%
E(%) =1.640533 -1.3695
1.640533*100 = 16.5%
![Page 24: Métodos Numéricos - juanrodriguezprieto.files.wordpress.com · Ejemplo de la regla del trapecio f ( x ) = 0.2 + 25 x - 200 x 2 + 675 x 3 - 900 x 4 + 400 x 5 Aproxime la integral](https://reader030.vdocumento.com/reader030/viewer/2022021522/5c5d427309d3f2d72f8cb1b7/html5/thumbnails/24.jpg)
Ejemplo de la regla del trapecio
múltiple en Matlab
n = 5;a = 0;b = 0.8;x = linspace(a,b,n+1);x2 = x.*x;x3 = x2.*x;x4 = x3.*x;x5 = x4.*x;y = 0.2+25*x-200*x2+675*x3-900*x4+400*x5;
integral = y(1);for i = 2:n
integral = integral + 2*y(i);endintegral = integral + y(n+1);h = (b-a)/n;integral = integral*h/2;integral
![Page 25: Métodos Numéricos - juanrodriguezprieto.files.wordpress.com · Ejemplo de la regla del trapecio f ( x ) = 0.2 + 25 x - 200 x 2 + 675 x 3 - 900 x 4 + 400 x 5 Aproxime la integral](https://reader030.vdocumento.com/reader030/viewer/2022021522/5c5d427309d3f2d72f8cb1b7/html5/thumbnails/25.jpg)
Ejemplo de la regla del trapecio
múltiple en Matlab
n h Integral
2 0,4000 1,0688
3 0,2667 1,3696
4 0,2000 1,4848
5 0,1600 1,5399
6 0,1333 1,5703
7 0,1143 1,5887
8 0,1000 1,6008
9 0,0889 1,6091
10 0,0800 1,6150
11 0,0727 1,6195
12 0,0667 1,6228
13 0,0615 1,6254
14 0,0571 1,6275
15 0,0533 1,6292
Resultados obtenidos para diferentes
valores de n
A medida que n incrementa, el valor de
la integral que se obtiene usando la regla
del trapecio múltiple se aproxima a la
solución analítica
for n= 2:15a = 0;b = 0.8;x = linspace(a,b,n+1);x2 = x.*x;x3 = x2.*x;x4 = x3.*x;x5 = x4.*x;y = 0.2+25*x-200*x2+675*x3-900*x4+400*x5;
integral = y(1);for i = 2:n
integral = integral + 2*y(i);endintegral = integral + y(n+1);h = (b-a)/n;integral = integral*h/2;resultado(n-1,:)= [n,h, integral];end
![Page 26: Métodos Numéricos - juanrodriguezprieto.files.wordpress.com · Ejemplo de la regla del trapecio f ( x ) = 0.2 + 25 x - 200 x 2 + 675 x 3 - 900 x 4 + 400 x 5 Aproxime la integral](https://reader030.vdocumento.com/reader030/viewer/2022021522/5c5d427309d3f2d72f8cb1b7/html5/thumbnails/26.jpg)
Regla de Simpson 1/3Además de aplicar la regla del trapecio con una
segmentación fina, otra forma de obtener una
estimación más exacta de una integral consiste en
usar polinomios de grados superior para unir los
puntos.
Por ejemplo, otro punto entre la mitad entre f(a)
y f(b), los tres puntos se pueden unir con una
parábola. Si hay dos puntos igualmente
espaciados entre f(a) y f(b), los cuatro puntos se
pueden unir con mediante un polinomio de
tercer grado.
Las formulas que resultan de tomar las integrales
bajo esos polinomios se conocen como reglas de
Simpson
f (x0 ) = a0 + a1x0 + a2x0
2
f (x1) = a0 + a1x1 + a2x1
2
f (x2 ) = a0 + a1x2 + a2x2
2
1 x0 x0
2
1 x1 x1
2
1 x2 x2
2
é
ë
êêêê
ù
û
úúúú
a0
a1
a2
é
ë
êêê
ù
û
úúú
=
f (x0 )
f (x1)
f (x2 )
é
ë
êêê
ù
û
úúú
a0
a1
a2
é
ë
êêê
ù
û
úúú
=
-(x1x2 ) / (x0x1 + x0x2 - x1x2 - x0
2 ) -(x0x2 ) / (x0x1 - x0x2 + x1x2 - x1
2 ) (x0x1) / (x0x1 - x0x2 - x1x2 + x2
2 )
(x1 + x2 ) / (x0x1 + x0x2 - x1x2 - x0
2 ) (x0 + x2 ) / (x0x1 - x0x2 + x1x2 - x1
2 ) -(x0 + x1) / (x0x1 - x0x2 - x1x2 + x2
2 )
-1/ (x0x1 + x0x2 - x1x2 - x0
2 ) -1/ (x0x1 - x0x2 + x1x2 - x1
2 ) 1 / (x0x1 - x0x2 - x1x2 + x2
2 )
é
ë
êêêê
ù
û
úúúú
f (x0 )
f (x1)
f (x2 )
é
ë
êêê
ù
û
úúú
![Page 27: Métodos Numéricos - juanrodriguezprieto.files.wordpress.com · Ejemplo de la regla del trapecio f ( x ) = 0.2 + 25 x - 200 x 2 + 675 x 3 - 900 x 4 + 400 x 5 Aproxime la integral](https://reader030.vdocumento.com/reader030/viewer/2022021522/5c5d427309d3f2d72f8cb1b7/html5/thumbnails/27.jpg)
I = f (x)dxa
b
ò @ f2(x)dxa
b
ò =h
3f (x0 )+ 4 f (x1)+ f (x2 )[ ]
Regla de Simpson 1/3
I @ (b- a)f (x0 )+ 4 f (x1)+ f (x2 )
6
é
ëêù
ûú
I @ Ancho por altura promedio
![Page 28: Métodos Numéricos - juanrodriguezprieto.files.wordpress.com · Ejemplo de la regla del trapecio f ( x ) = 0.2 + 25 x - 200 x 2 + 675 x 3 - 900 x 4 + 400 x 5 Aproxime la integral](https://reader030.vdocumento.com/reader030/viewer/2022021522/5c5d427309d3f2d72f8cb1b7/html5/thumbnails/28.jpg)
Et = -1
90h5 f 4 (x)
Regla de Simpson 1/3
Error
Et = -(b- a)5
2880f 4 (x)
Así, la regla de Simpson 1/3 es más exacta que la regla del trapecio. Además,
es más exacta de lo esperado, porque en lugar de ser el error proporcional a la
tercera derivada, el error es proporcional a la cuarta derivada. En otras
palabras, da resultados exactos para polinomios cúbicos aun cuando se
obtenga un parábola.
![Page 29: Métodos Numéricos - juanrodriguezprieto.files.wordpress.com · Ejemplo de la regla del trapecio f ( x ) = 0.2 + 25 x - 200 x 2 + 675 x 3 - 900 x 4 + 400 x 5 Aproxime la integral](https://reader030.vdocumento.com/reader030/viewer/2022021522/5c5d427309d3f2d72f8cb1b7/html5/thumbnails/29.jpg)
Aplicación de la regla de Simpson
1/3
f (x) = 0.2 + 25x- 200x2 + 675x3 - 900x4 + 400x5
Aproxime el valor de la integral de la siguiente función
Usando la regla de Simpson 1/3
f (a = 0) = 0.2
f (x1 = 0.4) = 2.456
f (b = 0.8) = 0.232
I @ 0.80.2 + 4 *2.456 + 0.232
6
é
ëêù
ûú= 1.367467
![Page 30: Métodos Numéricos - juanrodriguezprieto.files.wordpress.com · Ejemplo de la regla del trapecio f ( x ) = 0.2 + 25 x - 200 x 2 + 675 x 3 - 900 x 4 + 400 x 5 Aproxime la integral](https://reader030.vdocumento.com/reader030/viewer/2022021522/5c5d427309d3f2d72f8cb1b7/html5/thumbnails/30.jpg)
Aplicación de la regla de Simpson
1/3
Error
f (x) = 0.2 + 25x- 200x2 + 675x3 - 900x4 + 400x5
Aproxime el valor de la integral de la siguiente función
Et =1.640533 -1.367467( )
1.640533100 = 16.6%
Es aproximadamente 5 veces más precisa que una sola aplicación de la regla del
trapecio
![Page 31: Métodos Numéricos - juanrodriguezprieto.files.wordpress.com · Ejemplo de la regla del trapecio f ( x ) = 0.2 + 25 x - 200 x 2 + 675 x 3 - 900 x 4 + 400 x 5 Aproxime la integral](https://reader030.vdocumento.com/reader030/viewer/2022021522/5c5d427309d3f2d72f8cb1b7/html5/thumbnails/31.jpg)
Regla de Simpson 1/3 de aplicación
múltiple
La regla de Simpson se mejora al dividir el intervalo de integración en
varios segmentos de un mismo tamaño.
Se debe utilizar un número par de segmentos para implementar el método
f (x) = 0.2 + 25x- 200x2 + 675x3 - 900x4 + 400x5
Aproxime el valor de la integral de la siguiente función
Usando la regla de Simpson 1/3 de aplicación múltiple
f (a = 0) = 0.2
f (x1 = 0.2) = 1.288
f (x2 = 0.4) = 2.456
f (x3 = 0.6) = 3.464
f (b = 0.8) = 0.232
I @ 0.80.2 + 4(1.288 + 3.464)+ 2(2.456)+ 0,232
12
é
ëêù
ûú= 1.623467
![Page 32: Métodos Numéricos - juanrodriguezprieto.files.wordpress.com · Ejemplo de la regla del trapecio f ( x ) = 0.2 + 25 x - 200 x 2 + 675 x 3 - 900 x 4 + 400 x 5 Aproxime la integral](https://reader030.vdocumento.com/reader030/viewer/2022021522/5c5d427309d3f2d72f8cb1b7/html5/thumbnails/32.jpg)
Regla de Simpson3/8
De manera similar a la obtención de la regla del trapecio y Simpson 1/3,
es posible ajustar un polinomio de Lagrange de tercer grado a cuatro
puntos e integrarlo
I @3h
8f (x0 )+ 3 f (x1)+ 3 f (x2 )+ f (x3)( )
Et = -(b- a)5
6480f 4 (x)
La regla 3/8 es más exacta que la regla de 1/3.
Por lo general, se prefiere la regla de Simpson 1/3, ya que alcanza una
exactitud de tercer orden con tres puntos en lugar de los cuatro puntos
requeridos en la versión 3/8. No obstante, la regla de 3/8 es útil cuando
el número de segmentos es impar.
![Page 33: Métodos Numéricos - juanrodriguezprieto.files.wordpress.com · Ejemplo de la regla del trapecio f ( x ) = 0.2 + 25 x - 200 x 2 + 675 x 3 - 900 x 4 + 400 x 5 Aproxime la integral](https://reader030.vdocumento.com/reader030/viewer/2022021522/5c5d427309d3f2d72f8cb1b7/html5/thumbnails/33.jpg)
Regla de Simpson3/8
f (x) = 0.2 + 25x- 200x2 + 675x3 - 900x4 + 400x5
Aproxime el valor de la integral de la siguiente función
Usando la regla de Simpson 3/8
f (a = 0) = 0.2
f (x1 = 0.2667) = 1.432724
f (x2 = 0.5333) = 3.487177
f (b = 0.8) = 0.232
I @ 0.80.2 + 3(1.432774 + 3.487177)+ 0.232
8
é
ëêù
ûú= 1.519170
Et =1.640533 -1.519170
1.640533100 = 7.4%
![Page 34: Métodos Numéricos - juanrodriguezprieto.files.wordpress.com · Ejemplo de la regla del trapecio f ( x ) = 0.2 + 25 x - 200 x 2 + 675 x 3 - 900 x 4 + 400 x 5 Aproxime la integral](https://reader030.vdocumento.com/reader030/viewer/2022021522/5c5d427309d3f2d72f8cb1b7/html5/thumbnails/34.jpg)
Tarea
8 + 4 cos(x)dx0
p /2
ò
Evalué la integral siguiente:
(a) De forma analítica
(b) Con una aplicación de la regla de trapecio
(c) Con aplicación múltiple de la regla de trapecio n=3
(d) Con una aplicación de la regla de Simpson 1/3
(e) Con una aplicación de la regla de Simpson 3/8