MÉTODO DEL TRAPECIO (Autoguardado) PDF
MÉTODO DEL TRAPECIO (Autoguardado) PDF
MÉTODO DEL TRAPECIO (Autoguardado) PDF
De Chiapas
Facultad de Ingeniería, Campus I.
Métodos Numéricos
𝟐𝒅𝒐 “C”
Fecha: 03/05/2018
1
MÉTODO DEL TRAPECIO
LA REGLA DE TRAPECIO.
La regla del trapecio es la primera de las fórmulas de integración cerradas de
Newton-Cotes. Corresponde al caso donde el polinomio de la ecuación es de
primer orden.
𝑏 𝑏
𝐼 = ∫ 𝑓(𝑥 )𝑑𝑥 ≅ ∫ 𝑓1 (𝑥 )𝑑𝑥
𝑎 𝑎
Tal que
𝑓 (𝑏 ) − 𝑓 (𝑎 )
𝑓1 (𝑥 ) = 𝑓 (𝑎) + (𝑥 − 𝑎 )
𝑏−𝑎
El área bajo esta línea recta es una aproximación de la integral de (𝑓(𝑥)) entre los
límites a y b:
𝑏
𝑓 (𝑏) − 𝑓(𝑎
𝐼 = ∫ [𝑓(𝑎) + (𝑥 − 𝑎)] 𝑑𝑥
𝑎 𝑏−𝑎
𝒇(𝒂) + 𝒇(𝒃)
𝑰 = (𝒃 − 𝒂)
𝟐
2
Trapecio simple
Es un método de integración que se efectúa de dos maneras, de forma
tabular, este a su vez se presenta, condicionada, o no condicionada y cuando existe
la presencia de una función, es un método que trabaja únicamente con dos puntos
estos son encontrados en la tabla cuando se realiza de forma tabular, y cuando se
presenta una función “a” vendría siendo el límite inferior y “b” el superior, carece de
números de iteraciones, son formas de cálculo directo.
En la gráfica se observa que el área sombreada entre los puntos f(a) y f(b) de la
función dada f(x) al evaluar la misma proporciona el valor exacto del área, mientras
que el trapecio formado por el intervalo a y b con longitud igual, el valor de la misma,
es igual a la integral calculada mediante la regla del trapecio.
𝒇(𝒂) + 𝒇(𝒃)
𝑰 = (𝒃 − 𝒂)
𝟐
ALGORITMO
• Paso 1: Identificar los límites del área bajo la curva que se desea calcular.
• Paso 2: Evaluar los límites en la función.
• Paso 2.1: si no se cuenta con una función seleccionar los datos
directamente de la tabulación.
• Paso 3: Sustituir en la fórmula de la regla del trapecio simple.
• Paso 4: dar solución al problema.
• Paso 5: En caso necesario calcular el error.
3
FORMULA PARA CALCULAR EL ERROR
𝟏
𝑬𝒕 = − 𝒇"(𝝃)(𝒃 − 𝒂)𝟑
𝟏𝟐
ALGORITMO
• Paso 1: Divida el intervalo [a, b] en subintervalos de su longitud h.
• Paso 2: Aproxime en cada subintervalo la función f(x) por una recta.
• Paso 3: Aproxime el área bajo la curva f en el intervalo [a, b] mediante la suma
de las áreas de los trapecios.
𝑏−𝑎
ℎ=
𝑛
𝒙𝟏 𝒙𝟐 𝒙𝒏
𝑰 = ∫ 𝒇(𝒙)𝒅𝒙 + ∫ 𝒇(𝒙)𝒅𝒙 + ⋯ + ∫ 𝒇(𝒙)𝒅𝒙
𝒙𝟎 𝒙𝟏 𝒙𝒏−𝟏
4
Y agrupando términos:
𝑛−1
ℎ
𝐼 = [𝑓(𝑥0 + 2 ∑ 𝑓 (𝑥1 ) + 𝑓(𝑥𝑛 )]
2
𝑖=1
para calcular el error en esta formula en caso de no tener la parte analitica es:
(𝑏 − 𝑎)3
𝐸𝑎 = 2
𝑓̅ ´´
12𝑛
5
EJEMPLOS:
Aplicación simple de la regla del trapecio.
Desde a=0 hatas b=0.8, se sabe que el area bajo la curva calculando mediante la
manera analitica es 1.640533 aplicando la regla siemple del trapecio.
Paso 1: Identificar los límites del área bajo la curva que se desea calcular.
a=0 b = 0.8
6
paso 3: Sustituir en la fórmula de la regla del trapecio simple.
𝒇(𝒂) + 𝒇(𝒃)
𝑰 = (𝒃 − 𝒂)
𝟐
Sustituyendo se tiene:
𝟎. 𝟐 + 𝟎. 𝟐𝟑𝟐
𝐈 = (𝟎. 𝟖 − 𝟎) = 𝟎. 𝟏𝟕𝟐𝟖
𝟐
La cual presenta un error de:
𝑬𝒕 = 1.640533 – 0.1728 = 1.467733
Que corresponde a un error relativo porcentual, de 𝐸𝑡 =89.5 %. La razón de que
sea tan grande este error es evidente en la gráfica. Observe que el área bajo la
línea recta no toma en cuenta una porción significativa de la integral que está por
encima de la línea.
Habrá situaciones donde tal vez no conozcan previamente el valor verdadero. Por
lo tanto se requerirá de una estimación de error aproximado. Y para obtener dicha
estimación se calculara la segunda derivada de la función original.
Función original
𝒇(𝒙) = 𝟎. 𝟐 + 𝟐𝟓𝒙 − 𝟐𝟎𝟎𝒙𝟐 + 𝟔𝟕𝟓𝒙𝟑 − 𝟗𝟎𝟎𝒙𝟒 + 𝟒𝟎𝟎𝒙𝟓
Primera derivada:
F´(x)= 25 + 400x+2025x^2-3600x^3+2000x^4
Segunda derivada:
f´´(x)= 400 + 4050x - 10800x^2 + 2000x^3
Quedando la función anterior.
Se calcula el valor promedio de la segunda derivada
0.8
∫ (−400 + 4050𝑥 − 10800𝑥 2 + 8000𝑥 3 )
𝑓´´(𝑥 ) = 0 = −60
0.8 − 0
7
Que se sustituye en la ecuación Para calcular el error.
𝟏
𝑬𝒕 = − 𝒇"(𝝃)(𝒃 − 𝒂)𝟑
𝟏𝟐
Que dando:
𝟏
𝑬𝒂 = − 𝟏𝟐 (−𝟔𝟎)(𝟎. 𝟖 − 𝟎)𝟑 = 𝟐. 𝟓𝟔 𝑼𝟐
Ejemplo 2:
Evalué la integral siguiente:
𝝅/𝟐
∫ (𝟔 + 𝟑𝒄𝒐𝒔𝒙)𝒅𝒙
𝟎
𝝅/𝟐
∫ (𝟔 + 𝟑𝒄𝒐𝒔𝒙)𝒅𝒙
𝟎
𝝅 𝝅
[𝟔 ( ) + 𝟑𝒔𝒆𝒏 ( ) − [𝟔(𝟎) + 𝟑𝒔𝒆𝒏(𝟎)] ] = 12.42477 𝑼𝟐
𝟐 𝟐
8
Intervalos: a= 0 b= 𝝅/𝟐
𝒇(𝒂) + 𝒇(𝒃)
𝑰 = (𝒃 − 𝒂)
𝟐
𝝅 𝟗 + 𝟔 𝟏𝟓𝝅
𝑰 = ( − 𝟎) = = 𝟏𝟏. 𝟕𝟖𝟎𝟗𝟕 𝑼𝟐
𝟐 𝟐 𝟒
𝑰 = 𝟏𝟏. 𝟕𝟖𝟎𝟗𝟕 𝑼𝟐
Margen de error:
𝟏𝟐.𝟒𝟐𝟒𝟕𝟕−𝟏𝟏.𝟕𝟖𝟎𝟗𝟕
𝜺𝒕 = [ ] ∗ 𝟏𝟎𝟎 = 𝟓. 𝟏𝟖%
𝟏𝟐.𝟒𝟐𝟒𝟕𝟕
9
𝝅 𝟗+𝟐(𝟖.𝟏𝟐𝟏𝟑𝟐)+𝟔
𝑰 = ( − 𝟎) = 𝟏𝟐. 𝟐𝟔𝟖𝟗𝟓 𝑼𝟐
𝟐 𝟐(𝟐)
𝟏𝟐.𝟒𝟐𝟒𝟕𝟕−𝟏𝟐.𝟐𝟔𝟖𝟗𝟓
Margen de error: 𝜺𝒕 = [ ] ∗ 𝟏𝟎𝟎 = 𝟏. 𝟐𝟓%
𝟏𝟐.𝟒𝟐𝟒𝟕𝟕
n=4
x f(x)
0 9
𝝅 8.77164
𝟖
𝝅 8.12132
𝟒
𝟏𝟕𝝅 7.123819
𝟒𝟓
𝝅 6
𝟐
𝝅 𝟗+𝟐(𝟖.𝟕𝟕𝟏𝟔𝟒+𝟖.𝟏𝟐𝟏𝟑𝟐+𝟕.𝟏𝟐𝟑𝟖𝟏𝟗)+𝟔
𝑰 = ( − 𝟎) = 𝟏𝟐. 𝟑𝟕𝟔𝟔𝟎𝟔 𝑼𝟐
𝟐 𝟐(𝟒)
Margen de error.
𝟏𝟐.𝟒𝟐𝟒𝟕𝟕−𝟏𝟐.𝟑𝟕𝟔𝟔𝟎𝟔
𝜺𝒕 = [ ] ∗ 𝟏𝟎𝟎 = 𝟎. 𝟑𝟖𝟕%
𝟏𝟐.𝟒𝟐𝟒𝟕𝟕
10
Calcular con el metodo de la regla de trapecio multiple con n=20.
n = 20
11
12