Scribd Logo
Importer

Bienvenue sur Scribd !

  • 55 vues

    Chap III Intégration Numérique

    Transféré par

    Youba Bekakria
    Ce document présente deux méthodes pour calculer une intégrale numériquement : la méthode des trapèzes et la méthode de Simpson. Il donne les formules pour ces méthodes et un exemple numérique pour comparer la précision des deux méthodes.

    Droits d'auteur :

    © All Rights Reserved
    Téléchargez comme PDF, TXT ou lisez en ligne sur Scribd

    Chap III Intégration Numérique

    Transféré par

    Youba Bekakria
    55 vues3 pages
    Ce document présente deux méthodes pour calculer une intégrale numériquement : la méthode des trapèzes et la méthode de Simpson. Il donne les formules pour ces méthodes et un exemple numérique pour comparer la précision des deux méthodes.

    Titre original

    Chap III Intégration numérique

    Copyright

    © © All Rights Reserved

    Formats disponibles

    PDF, TXT ou lisez en ligne sur Scribd

    Avez-vous trouvé ce document utile ?

    Ce contenu est-il inapproprié ?

    Ce document présente deux méthodes pour calculer une intégrale numériquement : la méthode des trapèzes et la méthode de Simpson. Il donne les formules pour ces méthodes et un exemple numérique pour comparer la précision des deux méthodes.

    Droits d'auteur :

    © All Rights Reserved
    Téléchargez comme PDF, TXT ou lisez en ligne sur Scribd
    Télécharger au format pdf ou txt
    55 vues3 pages

    Chap III Intégration Numérique

    Transféré par

    Youba Bekakria
    Ce document présente deux méthodes pour calculer une intégrale numériquement : la méthode des trapèzes et la méthode de Simpson. Il donne les formules pour ces méthodes et un exemple numérique pour comparer la précision des deux méthodes.

    Droits d'auteur :

    © All Rights Reserved
    Téléchargez comme PDF, TXT ou lisez en ligne sur Scribd
    Télécharger au format pdf ou txt
    sur 3

    III.

    Intégration numérique

    Soit 𝑓𝑓 une fonction définie dans un intervalle [𝑎𝑎 ,ܾ] et soit à calculer l’aire inclus dans cet
    intervalle compris entre la fonction 𝑓et l’axe des 𝑥 noté 𝐼𝐼 DSSHOpHLQWpJUDOHGHf TXLGRQQHUDXQH
    SULPLWLYHQRWpH)GpILQLHalors SDU:

    𝑏𝑏
    𝐼𝐼 = � 𝑓𝑓(𝑥𝑥)𝑑𝑑𝑑𝑑 = 𝐹𝐹 (𝑏𝑏) − 𝐹𝐹(𝑎𝑎)
    𝑎𝑎

    Or souvent la primitive F est inconnue ou que la fonction f , n’est connue que pour un nombre
    r fini de points
    dans ce cas, on calcul une valeur approximative de l’intégrale sous l'une des méthodes suivante :

    D
    ∗ Méthode des trapèzes généralisée
    Soient 𝑥𝑥0 , 𝑥𝑥1 , … . 𝑥𝑥𝑛𝑛 , (n+1) points équidistants dans un intervalle [𝑎𝑎 , 𝑏𝑏] ou
    𝑏𝑏−𝑎𝑎
    𝑥𝑥0 = 𝑎𝑎 et 𝑥𝑥𝑛𝑛 = 𝑏𝑏 alors 𝑥𝑥𝑖𝑖 = 𝑥𝑥0 + 𝑖𝑖ℎ avec ℎ = appelé le pas et n le nombre de sous

    I
    𝑛𝑛
    intervalles égaux formés par ces points [𝑥𝑥𝑖𝑖 , 𝑥𝑥𝑖𝑖+1 ] ; 𝑖𝑖 = 0, … , 𝑛𝑛 − 1
    SA
    𝑑𝑑 𝑛−1
    ℎ
    ≅ ( 𝑓(𝑎) + 𝑓(𝑏) + 2 � 𝑓(𝑥𝑖𝑖 ))
    𝑎𝑎 2 𝑖𝑖=1

    Cas particulier : pour n = 1


    N
    𝑏𝑏 ℎ 𝑏𝑏−𝑎𝑎
    𝐼𝐼 = ∫𝑎𝑎 𝑓𝑓 (𝑥𝑥)𝑑𝑑𝑑𝑑 ≅ ( 𝑓𝑓(𝑎𝑎) + 𝑓𝑓 (𝑏𝑏) ) = (𝑓𝑓(𝑎𝑎) + 𝑓𝑓(𝑏𝑏))
    2 2𝑛𝑛
    BE

    ∗ Méthode de Simpson généralisée

    L’intervalle [𝑎𝑎 , 𝑏𝑏] est divisé en 2n sous intervalles égaux alors :

    𝑛𝑛� 𝑛𝑛� −1
    𝑏𝑏 2 2
    ℎ
    𝐼𝐼 = � 𝑓𝑓 (𝑥𝑥)𝑑𝑑𝑑𝑑 ≅ ( 𝑓𝑓(𝑎𝑎) + 𝑓𝑓(𝑏𝑏) + 4 � 𝑓𝑓(𝑥𝑥2𝑖𝑖−1 ) + 2 � 𝑓𝑓(𝑥𝑥2𝑖𝑖 ))
    𝑎𝑎 3
    𝑖𝑖=1 𝑖𝑖=1

    Cas particulier : pour n = 2

    𝑏𝑏 ℎ 𝑏𝑏−𝑎𝑎
    𝐼𝐼 = ∫𝑎𝑎 𝑓𝑓 (𝑥𝑥 )𝑑𝑑𝑥𝑥 ≅ � 𝑓𝑓(𝑎𝑎) + 𝑓𝑓(𝑏𝑏) + 4𝑓𝑓(𝑥𝑥1 )� = � 𝑓𝑓(𝑎𝑎) + 𝑓𝑓(𝑏𝑏) + 4𝑓𝑓(𝑥𝑥1 )�
    3 6
    Evaluation de l’erreur commise sur les 2 méthodes par rapport à l’intégrale exacte

    ξ = | 𝐼𝐼 (𝑓𝑓 ) − 𝐼𝐼𝐼𝐼(𝑓𝑓 )|= intégral exact - intégral approché


    L’estimation par rapport à l’erreur théorique pour les trapèzes n = 1

    ξ ≤ ξ theo
    ( 𝑏𝑏 – 𝑎𝑎 )3 𝑀𝑀2
    ξ theo = M 2 = max|𝑓𝑓 ′′ (𝑥𝑥)|
    12.𝑛𝑛2
    𝑥𝑥 ∈ [𝑎𝑎, 𝑏𝑏]

    L’estimation par rapport à l’erreur théorique pour Simpson n = 2

    D
    ( 𝑏𝑏 – 𝑎𝑎 )5 𝑀𝑀4
    ξ theo = M 4 = max|𝑓𝑓 ⁗ (𝑥𝑥)|
    180.𝑛𝑛4
    R

    𝑥𝑥 ∈ [𝑎𝑎, 𝑏𝑏]

    I
    Application numérique
    SA
    On se propose d’approximer l’intégrale :

    2 1
    ∫1 𝑑𝑑𝑑𝑑 = [ln(𝑥𝑥 + 1)]12 = ln(3) − ln(2) = 0,40547
    𝑥𝑥+1
    P
    N
    Par les méthodes d’intégrations des trapèzes et de Simpson

    ∗ Formule des trapèzes généralisée


    𝑏𝑏 ℎ
    BE

    𝑛𝑛−1
    𝐼𝐼 = ∫𝑎𝑎 𝑓𝑓 (𝑥𝑥 )𝑑𝑑𝑑𝑑 ≅ ( 𝑓𝑓(𝑎𝑎) + 𝑓𝑓 (𝑏𝑏) + 2 ∑𝑖𝑖=1 𝑓𝑓 (𝑥𝑥𝑖𝑖 ) ) ; 𝑥𝑥𝑖𝑖 = 𝑥𝑥0 + 𝑖𝑖ℎ
    2

    𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃 𝑛𝑛 = 1
    2 𝑑𝑑𝑑𝑑 ℎ 𝑏𝑏−𝑎𝑎
    𝐼𝐼𝐼𝐼(𝑓𝑓) = ∫1 = �𝑓𝑓(𝑎𝑎) + 𝑓𝑓(𝑏𝑏)� = �𝑓𝑓(1) + 𝑓𝑓(2)� = 0.41667
    𝑥𝑥+1 2 2𝑛𝑛

    𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃 𝑛𝑛 = 2
    2
    𝑑𝑑𝑥𝑥 𝑏𝑏 − 𝑎𝑎
    𝐼𝐼𝑛𝑛(𝑓𝑓) = � = �𝑓𝑓(𝑎𝑎) + 2[𝑓𝑓(𝑥𝑥1 )] + 𝑓𝑓(𝑏𝑏)�
    1 𝑥𝑥 + 1 2𝑛𝑛
    3
    𝑥𝑥1 =
    2
    2
    𝑑𝑑𝑥𝑥 𝑏𝑏 − 𝑎𝑎
    � = �𝑓𝑓(1) + 2[𝑓𝑓(3/2)] + 𝑓𝑓(2)� = 0,40833
    1 𝑥𝑥 + 1 4
    ∗ Formule de Simpson généralisée (n pair )
    𝑛𝑛� 𝑛𝑛� −1
    𝑏𝑏 2 2
    ℎ
    𝐼𝐼 = � 𝑓𝑓(𝑥𝑥)𝑑𝑑𝑑𝑑 ≅ ( 𝑓𝑓(𝑎𝑎) + 𝑓𝑓(𝑏𝑏) + 4 � 𝑓𝑓(𝑥𝑥2𝑖𝑖−1 ) + 2 � 𝑓𝑓(𝑥𝑥2𝑖𝑖 ))
    𝑎𝑎 3
    𝑖𝑖=1 𝑖𝑖=1

    𝑃𝑃𝑃𝑃𝑃𝑃𝑃𝑃 𝑛𝑛 = 2

    2 𝑑𝑑𝑑𝑑 𝑏𝑏−𝑎𝑎 3
    𝐼𝐼𝐼𝐼(𝑓𝑓 ) = ∫1 = (𝑓𝑓(𝑎𝑎) + 𝑓𝑓(𝑏𝑏) + 4[𝑓𝑓 (𝑥𝑥1 )]) ; 𝑥𝑥1 =
    𝑥𝑥+1 3𝑛𝑛 2

    2
    𝑑𝑑𝑑𝑑 1
    � = (𝑓𝑓(1) + 𝑓𝑓(2) + 4[𝑓𝑓(3/2)]) = 0,40555
    1 𝑥𝑥 + 1 6

    D
    pour n = 4

    In(f)=∫₁² 1/(x+1)dx = 1/8.[1/2+1/3+2.(f(5/4))+f(3/2))+f(7/4))]= 0.40619 Trapèze


    In(f)=∫₁² 1/(x+1)dx = 1/12.[1/2+1/3+4.(f(5/4))+f(7/4))+2.f(3/2)]= 0.40547 Simpson

    I
    SA
    n Méthode des trapèzes Méthode de Simpson

    1 0,41667
    2 0,40833 0,40555
    4 0,40619 0,40547
    N
    BE

    Vous aimerez peut-être aussi