Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 33 vistas

    Tare2 - Toapanta Gilson

    Cargado por

    Jilson Toapanta
    Este documento describe y compara varios métodos numéricos para aproximar el cálculo de integrales definidas. El método del trapecio divide el área bajo la curva en trapecios y aproxima el área total como la suma de las áreas de los trapecios. El método de Boole usa cinco puntos para ajustar un polinomio de cuarto orden e integrarlo. El método de Gauss-Legendre determina coeficientes óptimos para aproximar la integral como una suma ponderada de los valores de la función en dos puntos.

    Copyright:

    © All Rights Reserved
    Descargue como PDF, TXT o lea en línea desde Scribd

    Tare2 - Toapanta Gilson

    Cargado por

    Jilson Toapanta
    33 vistas10 páginas
    Este documento describe y compara varios métodos numéricos para aproximar el cálculo de integrales definidas. El método del trapecio divide el área bajo la curva en trapecios y aproxima el área total como la suma de las áreas de los trapecios. El método de Boole usa cinco puntos para ajustar un polinomio de cuarto orden e integrarlo. El método de Gauss-Legendre determina coeficientes óptimos para aproximar la integral como una suma ponderada de los valores de la función en dos puntos.

    Título original

    Tare2_Toapanta Gilson

    Derechos de autor

    © © All Rights Reserved

    Formatos disponibles

    PDF, TXT o lea en línea desde Scribd

    ¿Le pareció útil este documento?

    ¿Este contenido es inapropiado?

    Este documento describe y compara varios métodos numéricos para aproximar el cálculo de integrales definidas. El método del trapecio divide el área bajo la curva en trapecios y aproxima el área total como la suma de las áreas de los trapecios. El método de Boole usa cinco puntos para ajustar un polinomio de cuarto orden e integrarlo. El método de Gauss-Legendre determina coeficientes óptimos para aproximar la integral como una suma ponderada de los valores de la función en dos puntos.

    Copyright:

    © All Rights Reserved
    Descargue como PDF, TXT o lea en línea desde Scribd
    Descargar como pdf o txt
    33 vistas10 páginas

    Tare2 - Toapanta Gilson

    Cargado por

    Jilson Toapanta
    Este documento describe y compara varios métodos numéricos para aproximar el cálculo de integrales definidas. El método del trapecio divide el área bajo la curva en trapecios y aproxima el área total como la suma de las áreas de los trapecios. El método de Boole usa cinco puntos para ajustar un polinomio de cuarto orden e integrarlo. El método de Gauss-Legendre determina coeficientes óptimos para aproximar la integral como una suma ponderada de los valores de la función en dos puntos.

    Copyright:

    © All Rights Reserved
    Descargue como PDF, TXT o lea en línea desde Scribd
    Descargar como pdf o txt
    de 10

    MÉTODO DEL TRAPECIO

    Para resolver una integral definida existen métodos numéricos que nos permiten obtener
    aproximaciones cercanas al resultado real. El método del trapecio es un método numérico
    cuya estrategia es representar el área bajo la curva a través de un trapecio, en este caso se
    presentará un solo trapecio bajo la curva.[1]

    Figura 1. Método del trapecio con n=1.

    Observando la gráfica, se obtiene el área aproximada de la curva mediante:


    𝑏 (𝑓(𝑏) + 𝑓(𝑎))
    ∫ 𝑓(𝑥)𝑑𝑥 ≅ (𝑏 − 𝑎) ∙
    𝑎 2

    Sin embargo, es posible esta mejorar esta aproximación tomando dos trapecios, es decir, n es
    igual a 2, y para ello se divide el intervalo de integración [a, b], a la mitad para tener dos
    trapecios [1].

    Figura 2. Método del trapecio con n=2.

    De esta manera se tiene que el área bajo la curva está dada por dos trapecios:
    𝑏 (𝑓(𝑥1 ) + 𝑓(𝑎)) (𝑓(𝑏) + 𝑓(𝑥1 ))
    ∫ 𝑓(𝑥)𝑑𝑥 ≅ (𝑥1 − 𝑎) ∙ + (𝑏 − 𝑥1 ) ∙
    𝑎 2 2
    Siguiendo la regularidad vista anteriormente, si se toma n trapecios la aproximación de la
    integral queda de la siguiente manera:
    𝑏
    ℎ
    ∫ 𝑓(𝑥)𝑑𝑥 ≅ [𝑓(𝑎) + 2𝑓(𝑥1 ) + 2𝑓(𝑥2 ) + ⋯ + 2𝑓(𝑥𝑛−1 ) + 𝑓(𝑏)]
    𝑎 2

    𝑏 𝑛−1
    ℎ
    ∫ 𝑓(𝑥)𝑑𝑥 ≅ [𝑓(𝑎) + 2 ∑ 𝑓(𝑥𝑖 ) + 𝑓(𝑏)]
    𝑎 2
    𝑖=1

    Donde h será igual a:


    𝑏−𝑎
    ℎ=
    𝑛

    y los valores de cada punto será:


    𝑥1 = 𝑎 + ℎ
    𝑥2 = 𝑎 + 2ℎ

    𝑥𝑛−1 = 𝑎 + (𝑛 − 1) ∙ ℎ
    MÉTODO DE BOOLE
    La regla de Regla de Boole también denominada en algunos libros como método de Milne
    es aquella en la que utilizamos 5 puntos de una función determinada para la aproximación
    del área bajo la función, tal y como se muestra en la figura [3]:

    Figuar 3. Método de Boole.

    De manera similar al método de Simpson 1/3 es posible ajustar un polinomio de Lagrange de cuarto orden de
    5 puntos e integrarlo [2]:
    𝑏 𝑏
    ∫ 𝑓(𝑥)𝑑𝑥 ≅ ∫ 𝑓4 (𝑥)𝑑𝑥
    𝑎 𝑎
    𝑏 (𝑏 − 𝑎)
    ∫ 𝑓(𝑥)𝑑𝑥 ≅ [7𝑓(𝑎) + 32𝑓(𝑥1 ) + 12𝑓(𝑥2 ) + 32𝑓(𝑥3 ) + 7𝑓(𝑏)]
    𝑎 90
    𝑏
    2ℎ
    ∫ 𝑓(𝑥)𝑑𝑥 ≅ [7𝑓(𝑎) + 32𝑓(𝑥1 ) + 12𝑓(𝑥2 ) + 32𝑓(𝑥3 ) + 7𝑓(𝑏)]
    𝑎 45

    Donde h es igual a:
    𝑏−𝑎
    ℎ=
    4

    y los valores de cada punto será:


    𝑥1 = 𝑎 + ℎ

    𝑥2 = 𝑎 + 2ℎ
    Por otro lado se puede obtener una fórmula de Boole compuesta en grupos de 5 puntos, esto
    2ℎ
    se obtiene al realizar una sumatoria de la formula anteriormente vista, [7𝑓 (𝑥0 ) +
    45
    32𝑓 (𝑥1 ) + 12𝑓 (𝑥2 ) + 32𝑓 (𝑥3 ) + 7𝑓 (𝑥4 )], por lo que se tiene [4]:
    𝑏
    2
    ∫ 𝑓(𝑥)𝑑𝑥 ≈ ℎ(7𝑓(𝑥0 ) + 32𝑓(𝑥1 ) + 12𝑓(𝑥2 ) + 32𝑓(𝑥3 ) + 7𝑓(𝑥4 ))
    𝑎 45
    2
    + ℎ(7𝑓(𝑥4 ) + 32𝑓(𝑥5 ) + 12𝑓(𝑥6 ) + 32𝑓(𝑥7 ) + 7𝑓(𝑥8 )) + ⋯
    45
    2
    + ℎ(7𝑓(𝑥4𝑛−4 ) + 32𝑓(𝑥4𝑛−3 ) + 12𝑓(𝑥4𝑛−2 ) + 32𝑓(𝑥4𝑛−1 ) + 7𝑓(𝑥4𝑛 ))
    45

    Observando la formula se tiene dos elementos extremos que denominaremos como:

    𝐸 = 𝑓 (𝑥0 ) + 𝑓 (𝑥4𝑛 )

    Además, se tiene puntos impares que podremos denotar como:


    𝐼𝑖𝑚𝑝 = 𝑓(𝑥1 ) + 𝑓(𝑥3 ) + 𝑓(𝑥5 ) + 𝑓(𝑥7 ) + ⋯ + 𝑓(𝑥4𝑛−3 ) + 𝑓(𝑥4𝑛−1 )

    De igual forma se tiene elementos pares por lo que se le denominará como:


    𝐼𝑝𝑎𝑟 = 𝑓(𝑥2 ) + 𝑓(𝑥6 ) + 𝑓(𝑥10 ) + ⋯ + 𝑓(𝑥4𝑛−2 )

    Y por último los elementos sobrantes se les denominará como extremos acumulables:
    𝐼𝑒𝑥𝑡_𝑎𝑐𝑢 = 𝑓(𝑥4 ) + 𝑓(𝑥8 ) + ⋯ + 𝑓(𝑥4𝑛−4 )

    De forma más reducida se tiene que la fórmula de Boole compuesta es igual a:


    𝑏
    2
    ∫ 𝑓(𝑥)𝑑𝑥 ≈ ℎ(7𝐸 + 32𝐼𝑖𝑚𝑝 + 12𝐼𝑝𝑎𝑟 + 14𝐼𝑒𝑥𝑡_𝑎𝑐𝑢 )
    𝑎 45
    MÉTODO GAUSS-LEGENDRE DE DOS PUNTOS HASTA 6 PUNTOS

    Figura 4. Método de Gauss-Legendre.

    Las fórmulas particulares de cuadratura de Gauss se denominan fórmulas de Gauss-


    Legendre. El objetivo de la cuadratura de Gauss es determinar los coeficientes de una
    ecuación de la forma:
    𝐼 ≅ 𝑐0𝑓 (𝑥0 ) + 𝑐1 𝑓(𝑥1 )
    .

    Donde 𝑐0 y 𝑐1 son valores a calcular. Debido a que los puntos 𝑥1 y 𝑥0 no se encuentran en


    los extremos, estos valores también se convierten en variables a determinar. Para determinar
    estos valores se supone una serie de condiciones, donde se supone que la ecuación
    𝑐0 𝑓(𝑥0 ) + 𝑐1 𝑓(𝑥1 ) se ajusta con exactitud:
    A la integral de una constante:
    1
    𝑐0 𝑓(𝑥0 ) + 𝑐1𝑓 (𝑥1 ) = ∫ 1𝑑𝑥
    −1

    𝑐0 𝑓(𝑥0 ) + 𝑐1 𝑓(𝑥1 ) = 2
    𝑐0 + 𝑐1 = 2
    A la integral de una función lineal:

    1
    𝑐0 𝑓(𝑥0 ) + 𝑐1 𝑓(𝑥1 ) = ∫ 𝑥𝑑𝑥
    −1

    𝑐0 𝑥0 + 𝑐1 𝑥1 = 0
    A la integral de una función parabólica:
    1
    𝑐0𝑓 (𝑥0 ) + 𝑐1 𝑓(𝑥1 ) = ∫ 𝑥 2 𝑑𝑥
    −1

    2
    𝑐0 𝑥0 2 + 𝑐1 𝑥1 2 =
    3
    Y a la integral de una función cúbica:

    1
    𝑐0𝑓 (𝑥0 ) + 𝑐1 𝑓(𝑥1 ) = ∫ 𝑥 3 𝑑𝑥
    −1

    𝑐0 𝑥0 3 + 𝑐1 𝑥1 3 = 0

    Se tiene un sistema de ecuaciones:


    𝑐0 + 𝑐1 = 2
    𝑐0 𝑥0 + 𝑐1 𝑥1 = 0
    2
    𝑐0 𝑥0 2 + 𝑐1 𝑥1 2 =
    3
    𝑐0 𝑥0 3 + 𝑐1 𝑥1 3 = 0
    De donde se despeja 𝑐0 y 𝑥0 de las dos primeras ecuaciones:
    𝑐0 + 𝑐1 = 2
    𝑐0 = 2 − 𝑐1
    𝑐0 𝑥0 + 𝑐1 𝑥1 = 0
    𝑐1𝑥1
    𝑥0 = −
    𝑐0
    𝑐1𝑥1
    𝑥0 = −
    2 − 𝑐1
    A continuación se reemplaza 𝑐0 y 𝑥0 en la cuarta ecuación:
    𝑐0 𝑥0 3 + 𝑐1 𝑥1 3 = 0
    𝑐1𝑥1 3
    (2 − 𝑐1) (− ) + 𝑐1𝑥1 3 = 0
    2 − 𝑐1
    𝑐1 3 𝑥1 3
    − + 𝑐1𝑥1 3 = 0
    (2 − 𝑐1 )2
    𝑐1 2
    𝑐1𝑥1 3 (− + 1) = 0
    (2 − 𝑐1)2

    𝑐1 2
    (− + 1) = 0
    (2 − 𝑐1)2

    𝑐1 2
    =1
    (2 − 𝑐1)2

    𝑐1 2 = (2 − 𝑐1 )2
    𝑐12 = 4 − 4𝑐1 + 𝑐1 2
    0 = 4 − 4𝑐1
    4 = 4𝑐1
    𝑐1 = 1
    Se obtiene 𝑐1 , y se reemplaza en 𝑐0 = 2 − 𝑐1 :
    𝑐0 = 2 − 1
    𝑐0 = 1
    𝑐 𝑥
    1 1
    En la ecuación 𝑥0 = − 2−𝑐 se reemplaza los valores de 𝑐1 y 𝑐0:
    1

    𝑐1𝑥1
    𝑥0 = −
    2 − 𝑐1
    (1)𝑥1
    𝑥0 = −
    2−1
    𝑥0 = −𝑥1
    2
    En la ecuación 𝑐0 𝑥0 2 + 𝑐1 𝑥1 2 = 3 se reemplaza los valores de 𝑐1 , 𝑐0 Y 𝑥0 :

    2
    𝑐0 𝑥0 2 + 𝑐1 𝑥1 2 =
    3
    2
    𝑥0 2 + 𝑥1 2 =
    3
    2
    (−𝑥1 )2 + 𝑥1 2 =
    3
    2
    2𝑥1 2 =
    3
    1
    𝑥1 2 =
    3
    1
    𝑥1 = √
    3

    1
    𝑥1 =
    √3
    Reemplazando en 𝑥0 = −𝑥1 se tiene:
    𝑥0 = −𝑥1
    1
    𝑥0 = −
    √3
    Por ultimo reemplazamos los valores en:
    𝐼 ≅ 𝑐0𝑓 (𝑥0 ) + 𝑐1 𝑓(𝑥1 )
    1 1
    𝐼 ≅ (1)𝑓 (− ) + (1)𝑓 ( )
    √3 √3

    Observemos que el límite de integración se tomó en el intervalo [-1,1]. Para utilizar, en el


    caso que los límites de integración sean diferentes, se los cambia con un cambio de variable,
    𝑥 = 𝑎0 + 𝑎1 𝑥𝑑 . Entonces se toma el valor de a como límite inferior y 𝑥𝑑 = −1:
    𝑎 = 𝑎0 + 𝑎1 (−1)
    Y para el límite superior se toma x=b, y 𝑥𝑑 = 1:
    𝑏 = 𝑎0 + 𝑎1 (1)
    Al resolver las ecuaciones anteriores se obtiene que:
    𝑏+𝑎 𝑏−𝑎
    𝑎0 = 𝑦 𝑎1 =
    2 2
    Estos valores sustituimos en 𝑥 = 𝑎0 + 𝑎1 𝑥𝑑 :
    𝑏+𝑎 𝑏−𝑎
    𝑥= + ∙ 𝑥𝑑
    2 2

    Al derivar la ecuación se tiene:


    𝑏−𝑎
    𝑑𝑥 = ∙ 𝑑𝑥𝑑
    2
    Método Gauss-Legendre con más de dos puntos
    Aparte de la fórmula de dos puntos descrita en la sección anterior, se pueden desarrollar
    versiones con más puntos en la forma general:
    )
    𝐼 ≅ 𝑐0𝑓 (𝑥0 ) + 𝑐1 𝑓(𝑥1 ) + ⋯ + 𝑐𝑛−1 𝑓(𝑥𝑛−1
    donde n es el número de puntos. A continuación, se presenta una tabla con los
    valores de c y hasta de 6 puntos.

    Tabla 1. Fórmulas de Gauss-Legendre.

    Para obtener el error que produce al utilizar este método, se usará:


    𝑉𝑟𝑒𝑎𝑙 − 𝑉𝑚𝑒𝑑𝑖𝑑𝑜
    𝜀%𝑟 = ∗ 100%
    𝑉𝑟𝑒𝑎𝑙

    También podría gustarte