Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 152 vistas

    Integración Simpson y Trapecio Metodos Numericos

    Cargado por

    Jeancarlos Bonilla
    Este documento describe varios métodos numéricos para la integración, incluidas las reglas de Newton-Cotes, Trapecio y Simpson. La regla del trapecio aproxima la función a integrar con una línea recta entre los puntos inicial y final. La regla de Simpson 1/3 usa un polinomio cuadrático. Dividiendo el intervalo en más subintervalos y aplicando Simpson 1/3 de forma repetida aumenta la precisión. El documento también discute el cálculo del error asociado con estos métodos.

    Copyright:

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

    Integración Simpson y Trapecio Metodos Numericos

    Cargado por

    Jeancarlos Bonilla
    152 vistas25 páginas
    Este documento describe varios métodos numéricos para la integración, incluidas las reglas de Newton-Cotes, Trapecio y Simpson. La regla del trapecio aproxima la función a integrar con una línea recta entre los puntos inicial y final. La regla de Simpson 1/3 usa un polinomio cuadrático. Dividiendo el intervalo en más subintervalos y aplicando Simpson 1/3 de forma repetida aumenta la precisión. El documento también discute el cálculo del error asociado con estos métodos.

    Descripción original:

    De la materia de metodos numericos

    Título original

    Integración Simpson y Trapecio metodos numericos

    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 varios métodos numéricos para la integración, incluidas las reglas de Newton-Cotes, Trapecio y Simpson. La regla del trapecio aproxima la función a integrar con una línea recta entre los puntos inicial y final. La regla de Simpson 1/3 usa un polinomio cuadrático. Dividiendo el intervalo en más subintervalos y aplicando Simpson 1/3 de forma repetida aumenta la precisión. El documento también discute el cálculo del error asociado con estos métodos.

    Copyright:

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

    Integración Simpson y Trapecio Metodos Numericos

    Cargado por

    Jeancarlos Bonilla
    Este documento describe varios métodos numéricos para la integración, incluidas las reglas de Newton-Cotes, Trapecio y Simpson. La regla del trapecio aproxima la función a integrar con una línea recta entre los puntos inicial y final. La regla de Simpson 1/3 usa un polinomio cuadrático. Dividiendo el intervalo en más subintervalos y aplicando Simpson 1/3 de forma repetida aumenta la precisión. El documento también discute el cálculo del error asociado con estos métodos.

    Copyright:

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

    Universidad Politécnica Salesiana

    Sede Cuenca

    Curso: Métodos Numéricos

    ACTIVIDAD: Algoritmos de Newton-Cotes. Método de


    integración de Romberg.
    Introducción
    Las fórmulas de Newton-Cotes son las formas de
    integración numérica más empleadas. Las mismas
    se basan en la sustitución de la función que se
    desea integrar o de los datos tabulados
    (mediciones) por un polinomio de grado
    determinado, mucho más fácil de integrar.
    b b

    I =∫ f ( x)dx≈∫ f n (x)dx ;
    a a
    n n−1
    f n ( x)=an x +an−1 x +⋯+a1 x +a0
    Introducción

    a) polinomio de grado 1; b) polinomio de grado 2


    Introducción

    a) forma cerrada; b) forma abierta


    Regla del trapecio
    b b

    I =∫ f ( x)dx≈∫ f 1 ( x)dx ;
    a a

    f 1 ( x)=a1 x +a0

    Se tienen las coordenadas de los puntos extremos del


    intervalo de integración: (a, f(a)), (b, f(b)). Y por ellos se
    traza una recta, cuya ecuación es:
    Regla del trapecio

    I≈∫
    a
    [ f (a)+
    f (b)−f (a)
    b−a ]
    ( x−a) dx ;

    f (a)+ f (b)
    I =(b−a)
    2
    Fórmula para Regla del Trapecio.
    Ejemplo 1: Calcular, aplicando la regla del
    trapecio:
    4

    ∫ (1−e −2x
    )dx
    0
    0 −2×4
    f (a)+ f (b) (1−e )+(1−e )
    I =(b−a) =(4−0)
    2 2
    1
    1− 8
    e 1
    I=4 =2(1− )≈2×0.999665≈1.999339
    2 2980.958
    Ejemplo 1:
    Error de la Regla del Trapecio

    1 '' 3
    E t =− f (ξ)(b−a) ; ξ∈(a , b)
    12
    Et es el error verdadero de truncamiento. Sin embargo, es
    indeterminado el valor de ξ.
    Error aproximado

    1 ̄ '' 3
    E a=− ( f ) ( x)(b−a) .
    12

    f ( x)=1−e
    −2x ' −2x ' ' −2x
    ; f (x)=2 e ; f ( x)=−4 e ;
    Error de la Regla del Trapecio
    b
    ''
    ∫f ( x)dx
    ̄ '' a
    ( f ) ( x)=
    b−a
    Para el ejemplo 1:
    4
    1
    ∫ (−4 e
    0
    −2x
    )dx=2 8 −1
    e ( )
    ' ' 1 1
    ( ̄f ) (x)= 2 8 −1
    4 e ( )
    Error de la Regla del Trapecio
    Para el ejemplo 1:

    1 1 1
    E a =−
    12 2 e(8 ) 3
    −1 (4) =2.665772

    4
    1 1
    2(e )
    ∫ (1−e −2x
    )dx=4+ 8
    −1 =3.500168
    0

    E v =I v −I a =3.500168−1.999339=1.500829
    1.500829
    Er= ≈42.9 ‰
    3.500168
    Reglas de Simpson

    Cuando el intervalo de integración se divide en varios


    subintervalos, se obtienen métodos de integración de
    mayor precisión. En tales casos se utilizan, para sustituir
    la función integrando o los valores tabulados, polinomios
    de grado n, en dependencia de la cantidad de puntos
    interiores que se tenga. Las fórmulas que se obtienen de
    tomar integrales bajo esos polinomios se llaman fórmulas
    de Simpson.
    Regla de Simpson 1/3
    b b

    I =∫ f ( x)dx≈∫ f 2 (x) dx ;
    a a

    x0+ x2 b−a x 2−x 0


    a=x 0 ; b=x 2 ; x 1= →h= ≡
    2 2 2
    f 2 ( x): polinomio de Lagrange de 2do grado

    x2

    I =∫
    x0
    [ ( x−x 1 )( x−x 2 )
    ( x 0 −x 1 )( x 0 −x 2 )
    f ( x 0 )+
    ( x−x 0 )( x−x 2 )
    ( x 1−x 0 )( x 1 −x 2 )
    f ( x 1 )+
    ( x−x 0 )( x−x 1 )
    ( x 2−x 0 )( x 2 −x 1 ) ]
    f ( x 2 ) dx
    Regla de Simpson 1/3

    h
    I = [ f ( x 0 )+4f ( x 1 )+f ( x 2 ) ]
    3
    5
    (b−a) (4 )
    E t =− f (ξ); ξ∈(a , b)
    2880
    5
    (b−a) (4 )
    E a=− ̄
    ( f ) (x)
    2880
    b
    (4)
    ∫f ( x)dx
    (4) a
    ̄
    ( f ) (x)=
    b−a
    Ejemplo 2: Calcular, aplicando la regla de
    Simpson 1/3:
    4

    ∫ (1−e −2x
    )dx
    0
    h 2
    I=
    3
    [ f ( x 0 )+4f ( x 1 )+f ( x 2 ]
    ) =
    3
    [ (1−e
    0
    )+4 (1−e
    −2×2
    )+(1−e
    −2×4
    )]

    I≈3.284268
    5
    (b−a) ̄ (4)
    E a=− ( f ) ( x)
    2880
    4

    ∫ (−16 e −2x
    )dx
    (4) 0 1 −8
    ( ̄f ) (x)= = ×8(e −1)
    4 4
    ̄ (4)
    ( f ) ( x)≈−1.9993291

    5
    (4)
    E a=− (−1.9993291)=0.71087257
    2880
    4
    1 1
    ∫ (1−e
    0
    −2x
    )dx=4 +
    (
    2 e 8 )
    −1 =3.500168

    E v =I v −I a =3.500168−3.284268=0.215899743
    0.215899743
    Er= ≈6.17 ‰
    3.500168
    Aplicación múltiple de los algoritmos de
    Newton-Cotes
    La precisión de las fórmulas de Newton-Cotes mejora
    considerablemente si el intervalo (a, b) se divide en n
    subintervalos de igual longitud, definiendo el paso como:

    b−a
    h=
    n
    Aplicación múltiple de Simpson 1/3
    b x2 x4 xn

    I =∫ f ( x)dx=∫ f ( x)dx +∫ f ( x)dx +⋯+ ∫ f ( x)dx .


    a x0 x2 x n−2

    Al sustituir la regla de Simpson 1/3 en cada integral


    se obtiene:

    f ( x 0 )+4f ( x 1 )+f ( x 2 ) f ( x 2 )+4f ( x 3 )+f ( x 4 ) f ( x n−2 )+4f ( x n−1 )+ f ( x n )


    I ≈2h +2h +⋯+2h
    6 6 6

    (b−a)5 (4)
    E a =− 4
    ( ̄f ) ( x)
    180 n
    Ejemplo 3: Calcular, aplicando la regla de
    Simpson 1/3 con n = 4:
    4
    4−0
    ∫ (1−e −2x
    )dx ; h=
    4
    =1
    0

    (0,1) (1,2) (2,3) (3,4)


    x 0=0 ; x 1=1; x 2=2 ; x 3 =3 ; x 4 =4
    2 4

    I=∫ f 2 ( x)dx+∫ f 2 (x)dx


    0 2
    f (x 0 )+4f ( x 1 )+ f (x 2 ) f ( x 2 )+4f ( x 3 )+f ( x 4 )
    I ≈2h +2h
    6 6

    I≈2 [
    f (0)+ 4f (1)+f (2)
    6
    +2 ] [
    f (2)+4f (3)+ f (4)
    6 ]
    −2 −4
    f (0)=0 ; f (1)=1−e ; f (2)=1−e ;
    −6 −8
    f (3)=1−e ; f (4)=1−e
    I≈3.470592
    E v =I v −I a =3.500168−3.470592=0.029576
    Regla de Simpson 3/8

    b b

    I =∫ f ( x)dx≈∫ f 3 (x) dx ;
    a a

    3h b−a
    I = [ f ( x 0 )+3f (x 1 )+3f ( x 2 )+f ( x 3 ) ] ; h=
    8 3
    5
    (b−a) (4)
    E t =− f (ξ); ξ∈(a , b)
    6480
    Integración de Romberg

    Esta es una técnica para obtener integración numérica


    más eficiente. Los métodos que están agrupados bajo
    esta técnica usan dos estimaciones de una integral para
    calcular una tercera estimación, de mayor exactitud. A
    estas técnicas se les conoce como extrapolación de
    Richardson.

    -Estudiar el epígrafe 22.2 “Integración de Romberg”, libro


    de texto básico, desde la página 649 a la 655.
    Conclusiones
    Deberes
    - Estudiar, en el libro de texto básico, el capítulo 21
    “Fórmulas de integración de Newton-Cotes”, epígrafes 21.1,
    21.2; páginas de la 619 a la 640.

    - Resolver el problema 21.2, literales a), b), d) y e) de la


    página 646 del libro de texto básico.

    También podría gustarte