Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 78 vistas

    Ejemplos de Integración Numerica - Metodos Numericos

    Cargado por

    Angel Banda
    1) El documento describe el cálculo numérico de la fuerza y la posición de la línea de acción en un mástil de velero mediante las reglas de integración de Newton-Cotes y Simpson. 2) Se aplican las reglas trapezoidal compuesta y Simpson 1/3 compuesta para calcular la fuerza total y la posición de la línea de acción para un mástil de 30 metros de altura, obteniendo resultados de 1402.7 N y 13.72 m respectivamente con la trapezoidal, y 1586.2 N y 12.022 m con Simpson 1

    Copyright:

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

    Ejemplos de Integración Numerica - Metodos Numericos

    Cargado por

    Angel Banda
    78 vistas6 páginas
    1) El documento describe el cálculo numérico de la fuerza y la posición de la línea de acción en un mástil de velero mediante las reglas de integración de Newton-Cotes y Simpson. 2) Se aplican las reglas trapezoidal compuesta y Simpson 1/3 compuesta para calcular la fuerza total y la posición de la línea de acción para un mástil de 30 metros de altura, obteniendo resultados de 1402.7 N y 13.72 m respectivamente con la trapezoidal, y 1586.2 N y 12.022 m con Simpson 1

    Descripción original:

    integrales

    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?

    1) El documento describe el cálculo numérico de la fuerza y la posición de la línea de acción en un mástil de velero mediante las reglas de integración de Newton-Cotes y Simpson. 2) Se aplican las reglas trapezoidal compuesta y Simpson 1/3 compuesta para calcular la fuerza total y la posición de la línea de acción para un mástil de 30 metros de altura, obteniendo resultados de 1402.7 N y 13.72 m respectivamente con la trapezoidal, y 1586.2 N y 12.022 m con Simpson 1

    Copyright:

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

    Ejemplos de Integración Numerica - Metodos Numericos

    Cargado por

    Angel Banda
    1) El documento describe el cálculo numérico de la fuerza y la posición de la línea de acción en un mástil de velero mediante las reglas de integración de Newton-Cotes y Simpson. 2) Se aplican las reglas trapezoidal compuesta y Simpson 1/3 compuesta para calcular la fuerza total y la posición de la línea de acción para un mástil de 30 metros de altura, obteniendo resultados de 1402.7 N y 13.72 m respectivamente con la trapezoidal, y 1586.2 N y 12.022 m con Simpson 1

    Copyright:

    © All Rights Reserved
    Descargue como PDF, TXT o lea en línea desde Scribd
    Descargar como pdf o txt
    Está en la página 1de 6

    Ejemplos Integración numérica: Newton-Cotes

    P-19.10 La fuerza en el mástil de un velero se puede representar mediante la siguiente función:


    𝑧 𝑧
    𝑓 (𝑧) = 200( )𝑒 −2𝐻
    5+𝑧
    donde 𝑧 = la elevación sobre la cubierta y 𝐻 = la altura del mástil. La fuerza total 𝐹 ejercida sobre el mástil se
    puede determinar integrando esta función sobre la altura del mástil:
    𝐻
    𝐹 = ∫ 𝑓 (𝑧)𝑑𝑧
    0

    La línea de acción también se puede determinar mediante integración:


    𝐻
    ∫0 𝑧𝑓 (𝑧)𝑑𝑧
    𝑑= 𝐻
    ∫0 𝑓 (𝑧)𝑑𝑧

    a. Utilice la regla trapezoidal compuesta para calcular 𝐹 y 𝑑 para el caso donde 𝐻 = 30 (n = 6).
    𝑏 𝑓(𝑥0)+2 ∑𝑛−1
    𝑖=1 𝑓(𝑥𝑖 )+𝑓(𝑥𝑛 )
    Trapecio compuesto, ∫𝑎 𝑓(𝑥)𝑑𝑥 ≈ (𝑏 − 𝑎) 2𝑛

    30 𝑓(𝑧0 )+2 ∑6−1


    𝑖=1 𝑓(𝑧𝑖 )+𝑓(𝑧6 ) 30
    Fuerza resultante, 𝐹 = ∫0 𝑓 (𝑧)𝑑𝑧 ≈ (30 − 0) = (𝑓(𝑧0 ) +
    2∗6 12
    30−0
    2 ∑5𝑖=1 𝑓 (𝑧𝑖 ) + 𝑓 (𝑧𝑛 )), 𝑧𝑖 = 0 + 𝑖 =𝑖∗5
    6
    30
    5
    𝐹 = ∫ 𝑓(𝑧)𝑑𝑧 ≈ (𝑓 (0) + 2(𝑓(5) + 𝑓 (10) + 𝑓 (15) + 𝑓 (20) + 𝑓 (25)) + 𝑓(30))
    0 2
    30
    5
    𝐹 = ∫ 𝑓 (𝑧)𝑑𝑧 ≈ (0 + 2(71.653 + 68.456 + 55.182 + 42.176 + 31.479) + 23.200)
    0 2
    30
    𝐹 = ∫ 𝑓(𝑧)𝑑𝑧 ≈ 1402.7
    0
    𝐻 30
    ∫0 𝑧𝑓(𝑧)𝑑𝑧 ∫0 𝑧𝑓(𝑧)𝑑𝑧
    Posición de la línea de acción, 𝑑 = 𝐻 = 1402.7
    ∫0 𝑓(𝑧)𝑑𝑧

    30
    Cálculo de la integral, ∫0 𝑧𝑓 (𝑧)𝑑𝑧

    30 𝑧0 𝑓(𝑧0 )+2 ∑6−1


    𝑖=1 𝑧𝑖 𝑓(𝑧𝑖 )+𝑧6 𝑓(𝑧6 ) 30−0
    ∫0 𝑧𝑓(𝑧)𝑑𝑧 ≈ (30 − 0) 2∗6
    , 𝑧𝑖 = 0 + 𝑖 6
    =𝑖∗5
    30
    5
    ∫ 𝑧𝑓(𝑧)𝑑𝑧 ≈ (0𝑓(0) + 2(5𝑓(5) + 10𝑓(10) + 15𝑓(15) + 20𝑓(20) + 25𝑓(25)) + 30𝑓(30))
    0 2
    30
    ∫ 𝑧𝑓(𝑧)𝑑𝑧 ≈ 2.5 ∗ (0 + 2(5 ∗ 71.653 + 10 ∗ 68.456 + 15 ∗ 55.182 + 20 ∗ 42.176 + 25 ∗ 31.479) + 30 ∗ 23.200)
    0

    30
    ∫ 𝑧𝑓(𝑧)𝑑𝑧 ≈ 19245.250
    0

    𝐻
    ∫0 𝑧𝑓 (𝑧)𝑑𝑧 19245.250
    →𝑑= 𝐻 = = 13.720
    ∫0 𝑓 (𝑧)𝑑𝑧 1402.7

    b. Repita (a), pero use la regla compuesta 1∕3 de Simpson.


    𝑏 𝑓(𝑥0)+4 ∑𝑛−1 𝑛−2
    𝑖=1,3,5 𝑓(𝑥𝑖 )+2 ∑𝑗=2,4 𝑓(𝑥𝑗)+𝑓(𝑥𝑛 )
    Simpson 1/3, compuesto, ∫𝑎 𝑓 (𝑥)𝑑𝑥 ≈ (𝑏 − 𝑎) 3𝑛

    30 𝑓(𝑧0 )+4 ∑6−1 6−2


    𝑖=1,3,5 𝑓(𝑧𝑖 )+2 ∑𝑗=2,4 𝑓(𝑧𝑗 )+𝑓(𝑧6 ) 30−0
    Fuerza resultante, 𝐹 = ∫0 𝑓(𝑧)𝑑𝑧 ≈ (30 − 0) 3∗6
    , 𝑧𝑖 = 0 + 𝑖 6
    =𝑖∗5

    𝐹 ≈ (1.6667)(𝑓(𝑧0 ) + 4(𝑓(𝑧1 ) + 𝑓(𝑧3 ) + 𝑓(𝑧5 )) + 2(𝑓(𝑧2 ) + 𝑓(𝑧4 )) + 𝑓(𝑧6 ))

    𝐹 ≈ (1.6667)(0 + 4(71.653 + 68.456 + 55.182) + 2(42.176 + 31.479) + 23.200)


    30
    → 𝐹 = ∫ 𝑓 (𝑧)𝑑𝑧 ≈ 1586.2
    0
    𝐻 30
    ∫0 𝑧𝑓(𝑧)𝑑𝑧 ∫0 𝑧𝑓(𝑧)𝑑𝑧
    Posición de la línea de acción, 𝑑 = 𝐻 = 1586.2
    ∫0 𝑓(𝑧)𝑑𝑧

    30
    Cálculo de la integral, ∫0 𝑧𝑓 (𝑧)𝑑𝑧

    30 𝑧0 𝑓(𝑧0 )+4 ∑6−1 6−2


    𝑖=1,3,5 𝑧𝑖 𝑓(𝑧𝑖 )+2 ∑𝑖=2,4 𝑧𝑖 𝑓(𝑧𝑖 )+𝑧6 𝑓(𝑧6 ) 30−0
    ∫0 𝑧𝑓(𝑧)𝑑𝑧 ≈ (30 − 0) 3∗6
    , 𝑧𝑖 = 0 + 𝑖 6
    =𝑖∗5
    30
    5
    ∫ 𝑧𝑓(𝑧)𝑑𝑧 ≈ (0𝑓(0) + 4(5𝑓(5) + 10𝑓(10) + 15𝑓(15)) + 2(20𝑓(20) + 25𝑓(25) + 30𝑓(30))
    0 3
    30
    ∫ 𝑧𝑓(𝑧)𝑑𝑧 ≈ (1.667) ∗ (0 + 4(5 ∗ 71.653 + 10 ∗ 68.456 + 15 ∗ 55.182) + 2(20 ∗ 42.176 + 25 ∗ 31.479) + 30 ∗ 23.200)
    0

    30
    ∫ 𝑧𝑓(𝑧)𝑑𝑧 ≈ 19069.1631
    0

    𝐻
    ∫0 𝑧𝑓(𝑧)𝑑𝑧 19069.1631
    →𝑑= 𝐻 = = 12.022
    ∫0 𝑓(𝑧)𝑑𝑧 1586.2
    Ejemplos Integración numérica: Cuadratura de Gauss-Legendre

    a. Calculo analítico,
    8
    I = ∫ (−0.005𝑥 4 + 0.86𝑥 3 − 4.2𝑥 2 + 6.3𝑥 + 2)𝑑𝑥 = −0.001𝑥 5 + 0.215𝑥 4 − 1.4𝑥 3 + 3.15𝑥 2 + 2𝑥|80
    0

    𝐼 = 348.672

    c. Cuadratura de Gauss
    𝑏 𝑁
    𝑏−𝑎 𝑎+𝑏 𝑏−𝑎
    𝐼 = ∫ 𝑓 (𝑥)𝑑𝑥 ≈ ∑ 𝑤𝑁,𝑘 𝑓( + 𝑥𝑁,𝑘 )
    𝑎 2 2 2
    𝑘=1

    𝑓 (𝑥) = −0.005𝑥 + 0.86𝑥3 − 4.2𝑥2 + 6.3𝑥 + 2


    4

    3
    8
    8−0
    𝐼 = ∫ 𝑓(𝑥)𝑑𝑥 ≈ ∑ 𝑤3,𝑘 𝑓(4 + 4𝑥𝑁,𝑘 )
    0 2
    𝑘=1

    𝐼 ≈ 4(𝑤3,1 𝑓(4 + 4𝑥3,1 ) + 𝑤3,2 𝑓(4 + 4𝑥3,2 ) + 𝑤3,3 𝑓(4 + 4𝑥3,3 ))

    𝐼 ≈ 4(0.5556𝑓(4 − 4 ∗ 0.7746) + 0.8889𝑓(4 + 4 ∗ 0) + 0.5556𝑓(4 + 4 ∗ 0.7746))

    𝐼 ≈ 4 ∗ (0.5556 ∗ 4.8929590 + 0.8889 ∗ 13.760000 + 0.5556 ∗ 129.99436) = 348.699

    d. Mediante la función de intg de Scilab

    b. Cuadratura de Gauss
    𝑓 (𝑥) = 𝑥𝑒 2𝑥
    2
    3
    3−0 3 3
    𝐼 = ∫ 𝑓(𝑥)𝑑𝑥 ≈ ∑ 𝑤2,𝑘 𝑓 ( + 𝑥2,𝑘 ) = 1.5(𝑤2,1 𝑓(1.5 + 1.5𝑥2,1) + 𝑤2,2 𝑓(1.5 + 1.5𝑥2,2))
    0 2 2 2
    𝑘=1
    𝐼 ≈ 1.5(1 ∗ 𝑓(1.5 − 1.5 ∗ 0.5774) + 1 ∗ 𝑓(1.5 + 1.5 ∗ 0.5774))

    𝐼 ≈ 406.29

    c. Mediante la función de intg de Scilab

    a. Cuadratura de Gauss-Legendre con dos nodos


    2 1.5 −x2
    erf(1.5)= ∫ e dx
    √π 0
    2
    𝑓(𝑥) = 𝑒 −𝑥
    2
    2 1.5 2 1.5 − 0
    𝐼= ∫ 𝑓(𝑥)𝑑𝑥 ≈ ∑ 𝑤2,𝑘 𝑓(0.75 + 0.75𝑥2,𝑘 )
    √π 0 √π 2 𝑘=1

    2
    𝐼≈ 0.75(𝑤2,1 𝑓(0.75 + 0.75𝑥2,1 ) + 𝑤2,2 𝑓(0.75 + 0.75𝑥2,2 ))
    √π
    2
    𝐼≈ 0.75(1 ∗ 𝑓(0.75 − 0.75 ∗ 0.5774) + 1 ∗ (0.75 + 0.75 ∗ 0.5774))
    √π
    2
    𝐼≈ 0.75(1 ∗ 0.9044 + 1 ∗ 0.2467) = 0.8463
    √𝜋

    b. Cuadratura de Gauss-Legendre con tres nodos


    2
    𝑓(𝑥) = 𝑒 −𝑥
    3
    1.5
    2 2 1.5 − 0
    𝐼= ∫ 𝑓(𝑥)𝑑𝑥 ≈ ∑ 𝑤3,𝑘 𝑓(0.75 + 0.75𝑥3,𝑘 )
    √𝜋 0 √𝜋 2 𝑘=1

    2
    𝐼≈ 0.75(𝑤3,1 𝑓(0.75 + 0.75𝑥3,1 ) + 𝑤3,2 𝑓(0.75 + 0.75𝑥3,2 ) + 𝑤3,3 𝑓(0.75 + 0.75𝑥3,3 ))
    √𝜋
    2
    𝐼≈ 0.75(0.5556𝑓(0.75 − 0.75 ∗ 0.7746) + 0.8889𝑓(0.75 + 0.75 ∗ 0) + 0.5556𝑓(0.75 + 0.75 ∗ 0.7746))
    √𝜋
    2
    𝐼≈ 0.75 ∗ (0.5556 ∗ 0.9718241 + 0.8889 ∗ 0.5697828 + 0.5556 ∗ 0.1700933) = 0.9656
    √𝜋

    c. Error relativo con base a erf de Scilab

    Scilab → erf(1.5) = 0.9661


    |erf(1.5) − erfGL−2 | |0.9661 − 0.8463|
    Er(N: 2) = 100 = % = 12.4 %
    erf(1.5) 0.9661
    |erf(1.5) − erfGL−3 | |0.9661 − 0.9656|
    Er(N: 3) = 100 = 100 = 0.05 %
    erf(1.5) 0.9661

    b. Cuadratura de Gauss-Legendre con dos nodos


    𝑧
    𝑓 (𝑧) = 200( )𝑒 −2𝑧/30
    5+𝑧
    2
    30
    30 − 0
    𝐼 = ∫ 𝑓(𝑧)𝑑𝑧 ≈ ∑ 𝑤2,𝑘 𝑓(15 + 15𝑧2,𝑘 )
    0 2
    𝑘=1

    𝐼 ≈ 15(𝑤2,1 𝑓(15 + 15𝑧2,1 ) + 𝑤2,2 𝑓(15 + 15𝑧2,2 ))

    𝐼 ≈ 15(1 ∗ 𝑓(15 − 15 ∗ 0.5774) + 1 ∗ (15 + 15 ∗ 0.5774))


    30
    𝑧 −2𝑧/30
    ∫ 200( )𝑒 𝑑𝑧 ≈ 15 ∗ (1 ∗ 73.273 + 1 ∗ 34.273) = 1613.2
    0 5+𝑧

    c. Mediante la función de intg de Scilab

    Utilice el método de cuadratura de Gauss de dos puntos para estimar el valor promedio de la
    2
    siguiente función, f(x) = 2 entre a = 1, b = 5.
    1+x
    1 5 2
    ̅
    f[1,5] = 5−1 ∫1 dx
    1+x2

    2
    𝑓 (𝑥) =
    1 + x2
    2
    5
    5−1
    𝐼 = ∫ 𝑓(𝑥)𝑑𝑥 ≈ ∑ 𝑤2,𝑘 𝑓(3 + 2𝑥2,𝑘 )
    1 2
    𝑘=1

    𝐼 ≈ 2(𝑤2,1 𝑓(3 + 2𝑥2,1 ) + 𝑤2,2 𝑓(3 + 2𝑥2,2 ))

    𝐼 ≈ 2(1 ∗ 𝑓(3 − 2 ∗ 0.5774) + 1 ∗ (3 + 2 ∗ 0.5774))

    𝐼 ≈ 2 ∗ (1 ∗ 0.4540539 + 1 ∗ 0.1095149) = 1.1271


    1 5 2 1
    f̅[1,5] = 5−1 ∫1 1+x2
    dx ≈ 4 1.1271 = 0.2818

    Pesos 𝒘𝒊 nodos 𝒙𝒊 de Gauss-Legendre

    También podría gustarte