Chapitre 4 Intégration1
Chapitre 4 Intégration1
Chapitre 4 Intégration1
1- Description du problème
f fonction mal connue mais ne disposant pas de singularité sur [a; b].
1
Exemple :𝑓 (𝑥 ) = intégrable sur [0; 1] mais possède une singularité en 0.
√𝑥
Lorsqu’on connait une primitive de f (notée ici F) sur [a; b], on peut calculer directement I .
𝑏
𝐼 = ∫ 𝑓 (𝑥) 𝑑𝑥 = 𝐹 (𝑏) − 𝐹(𝑎)
𝑎
Problème :
La plupart des fonctions f ne disposent pas d’expressions analytique pour leurs primitives même dans le cas
de fonctions s’écrivant très simplement.
1 2 1
Exemples :∫0 𝑒 −𝑥 𝑑𝑥 ; ∫0 cos(𝑥 2 ) 𝑑𝑥
𝑏 𝑛
𝐼 (𝑓 ) = ∫ 𝑓 (𝑥 )𝑑𝑥 ≈ 𝐼𝑛 (𝑓 ) = ∑ 𝛼𝑖 𝑓 (𝑥𝑖 )
𝑎 𝑖=0
Définition 2 : Une formule de quadrature est dite exacte sur un ensemble V si l‘erreur de
l’approximation𝐸𝑛 (𝑓 ) = 𝐼(𝑓 ) − 𝐼𝑛 (𝑓 ) = 0
Définition 3 : Une formule de quadrature est dite de degré de précision n si elle est exacte pour
Formules globales:
Nous avons vu au chapitre précédent que nous pouvons approcher une fonction f par un polynôme
𝑏
d’interpolation 𝑃𝑛 (𝑥 ) = ∑𝑛𝑖=0 𝐿𝑖 (𝑥)𝑓 (𝑥𝑖 ) .L’idée va donc être de remplacer le calcul de 𝐼(𝑓 )par∫𝑎 𝑃𝑛 (𝑥 )𝑑𝑥.
Cela se justifie pour les deux raisons suivantes:
L’intégration de polynômes est très simple, et ne nécessite que les quatre opérations élémentaires.
Dans la pratique, assez souvent, on ne connait pas la forme analytique de f, mais seulement sa valeur
en certains points𝑥 𝑘 .
Définition 1 : On appelle formule de quadrature de type interpolation, la formule obtenue par le calcul
𝑏 𝑏
de∫𝑎 𝑃𝑛 (𝑥 )𝑑𝑥pour approcher∫𝑎 𝑓 (𝑥 )𝑑𝑥 , 𝑃𝑛 (𝑥 )étant un polynôme d’interpolation de f.
Proposition 1 : Une formule de quadrature à (n + 1) points est de degré au moins n, si et seulement si elle
est de type interpolation à n + 1 points.
2.1-Méthode de Newton-Cotes: plus généralement, les formules d’intégration qui peuvent s’écrire sous la
𝑥 +𝑛ℎ
forme∫𝑥 0 𝑓 (𝑥) 𝑑𝑥 ≈ ∑𝑛𝑖=0 𝐴𝑖 𝑓 (𝑥0 + 𝑖ℎ)sont connues sous le nom de « formule de Newton-Cotes ».
0
Les constantes 𝐻𝑖 sont appelées coefficients de Cotes, elles vérifient :∑𝑛𝑖=0 𝐻𝑖 = 1 ; 𝐻𝑖 = 𝐻𝑖−1
D’où
𝑏 𝑏
𝐼 (𝑓 ) = ∫𝑎 𝑓 (𝑥 )𝑑𝑥 = 𝐼 (𝑃1 ) = ∫𝑎 𝑃1 (𝑥 )𝑑𝑥
1 𝑏
ℎ
= ∑ 𝑓 (𝑥𝑘 ) ∫ 𝐿𝑘 (𝑥 )𝑑𝑥 = = (𝑓(𝑎) + 𝑓 (𝑏))
𝑎 2
𝑘=0
𝑏−𝑎 𝑥𝑛 − 𝑥0
= (𝑓(𝑎) + 𝑓 (𝑏))𝑎𝑣𝑒𝑐 ℎ =
2 𝑛
Remarque 1 :
Le trapèze composé revient à intégrer une interpolation linéaire par morceaux de la fonction.
Notons que
(𝑥 − 𝑥1 )(𝑥 − 𝑏) 𝑏
𝑏−𝑎
𝐿0 (𝑥 ) = ⇒ ∫ 𝐿0 (𝑥 )𝑑𝑥 =
(𝑎 − 𝑥1 )(𝑎 − 𝑏) 𝑎 6
(𝑥 − 𝑎)(𝑥 − 𝑏) 𝑏
4(𝑏 − 𝑎)
𝐿1 (𝑥 ) = ⇒ ∫ 𝐿1 (𝑥 )𝑑𝑥 =
(𝑥1 − 𝑎)(𝑥1 − 𝑏) 𝑎 6
(𝑥 − 𝑎)(𝑥 − 𝑥1 ) 𝑏
𝑏−𝑎
𝐿2 (𝑥 ) = ⇒ ∫ 𝐿0 (𝑥 )𝑑𝑥 =
(𝑏 − 𝑥1 )(𝑏 − 𝑎) 𝑎 6
Remarques
- Cette formule a un degré de précision de 3.
- Comme pour le trapèze, il est préférable de subdiviser en sous intervalle
𝑥 𝑥 𝑥
On trouve∫𝑥 𝑛 𝑓 (𝑥 )𝑑𝑥 = ∫𝑥 1 𝑓 (𝑥 )𝑑𝑥 + ∫𝑥 2 𝑓 (𝑥 )𝑑𝑥 + ⋯ +
0 0 1
𝑥𝑛 𝑓(𝑎)+𝑓(𝑏)
∫𝑥𝑛−1 𝑓 (𝑥 )𝑑𝑥 ≈ ℎ [ 2
+ 𝑓 (𝑥1 ) + 𝑓 (𝑥2 ) … + 𝑓 (𝑥𝑛−1 )] ≈
𝑏−𝑎 𝑓(𝑎)+𝑓(𝑏)
[ + ∑𝑛−1
𝑖=1 𝑓 (𝑥𝑖 )]
𝑛 2
On trouve
𝑥𝑛 𝑥1 𝑥2 𝑥𝑛
∫ 𝑓 (𝑥) 𝑑𝑥 = ∫ 𝑓 (𝑥) 𝑑𝑥 + ∫ 𝑓 (𝑥) 𝑑𝑥 + ⋯ + ∫ 𝑓 (𝑥 )𝑑𝑥
𝑥0 𝑥0 𝑥1 𝑥𝑛−1
ℎ
≈ [𝑓 (𝑥0 ) + 𝑓 (𝑥𝑛 ) + 4(𝑓 (𝑥1 ) + 𝑓 (𝑥3 ) + 𝑓 (𝑥5 ) + ⋯
3
+ 𝑓 (𝑥𝑛−1 )) + 2(𝑓 (𝑥2 ) + 𝑓 (𝑥4 ) + 𝑓 (𝑥6 ) + ⋯
+ 𝑓 (𝑥𝑛−2 ))]
ℎ ℎ𝑛
≈ 3 [𝑦0 + 𝑦𝑛 + 2 ∑𝑖 𝑝𝑎𝑖𝑟 𝑦𝑖 + 4 ∑𝑖 𝑖𝑚𝑝𝑎𝑖𝑟 𝑦𝑖 ] = {𝑓(𝑎) + 𝑓 (𝑏) + 2 ∑𝑘−1 𝑘−1
𝑖=1 𝑓2𝑘 + 4 ∑𝑖=0 𝑓2𝑘+1
3
Sur chaque triplet de sous intervalles, la courbe est remplacée par une cubique.
Pour chaque 4uplet de valeurs :
3ℎ
𝑆𝑖 = (𝑓 + 3𝑓1 + 3𝑓2 + 𝑓3 )
8 0
La règle composite de Simpson 3/8
Pour trouver l’intégrale I(f) est :
𝑏
3ℎ
∫ 𝑓(𝑥) 𝑑𝑥 = (𝑓 + 3𝑓1 + 3𝑓2 + 2𝑓3 + ⋯ + 2𝑓𝑛−3 + 3𝑓𝑛−2 + 3𝑓𝑛−1 + 𝑓𝑛 )
𝑎 8 0
Exercice 1
On lance une fusée verticalement du sol et l’on mesure pendant les premières 80 secondes l’accélération 𝛾:
t (en s) 0 10 20 30 40 50 60 70 80
𝛾(en m/s ) 30 31.63 33.44 35.47 37.75 40.33 43.29 46.70
2
50.67
Calculer la vitesse V de la fusée à l’instant t = 80s, par la méthode des trapèzes puis par Simpson.
80 ℎ 10
- 𝑇 = ∫0 𝛾(𝑠)𝑑𝑠 = 2 [𝛾(𝑡0 ) + 𝛾(𝑡8 ) + 2 ∑𝑛−1
𝑖=1 𝛾 (𝑡𝑖 )] = (30 + 50. .67 + 2(31.63 + 33.44 + ⋯ +
2
46.70)) = 3089, 5𝑚/𝑠.
ℎ
𝑆 = (𝛾 (𝑡0 ) + 𝛾 (𝑡8 ) + 4(𝛾(𝑡1 ) + 𝛾(𝑡3 ) + 𝛾(𝑡5 ) + 𝛾(7)) + 2(𝛾(𝑡2 ) + 𝛾(𝑡4 ) + 𝛾(𝑡6 )))
3
10
= [(30 + 50. .67 + 4(31.63 + 35.47 + 40.33 + 46.70) + 2
3
× (33.44 + 37.75 + 46.70)] = 3087, 2𝑚/𝑠.
Exercice 6
ℎ 𝜋
𝑇= (𝑓(𝑎) + 𝑓 (𝑏)) = (0 + 1) = 0.785
2 4
Soit S l’approximation de I par la méthode de Simpson. Celle-ci s’écrit
(𝑏 − 𝑎 ) 𝑎+𝑏 𝜋
𝑆= (𝑓(𝑎) + 4𝑓 ( ) + 𝑓(𝑏)) = (0 + 4 ∗ 0.707107 + 1)1 − 2027
6 2 12
ℎ 𝜋1
𝑆 = (𝑦0 + 𝑦4 + 4(𝑦1 + 𝑦3 ) + 2𝑦2 ) = [(0 + 1 + 4(0.38. . . + 0.92. . . ) + 2 × 0.707. . . )]
3 83
= 1.000135
𝜋
Les points d’appui donnés dans cet exercice correspondent à la fonction sin x. Et𝐼 = ∫0 𝑓 (𝑥) 𝑑𝑥 = 1On
2
constate donc
que l’approximation de I donnée par la méthode de Simpson est meilleure que celle par les trapèzes, puisque
|S − I| = 0.000135 et |T − I| = 0.012884.
2.6- Méthode de Newton (n=3) :
𝑏−𝑎
𝑥0 = 𝑎 ; 𝑥1 = 𝑥0 + ℎ, 𝑥2 = 𝑥0 + 2ℎ; 𝑥3 = 𝑥0 + 3ℎ = 𝑏; 𝑎𝑣𝑒𝑐 ℎ =
3
𝑏
L’approximation à l’ordre 3 de la fonctionnelle 𝐼 (𝑓 ) = ∫𝑎 𝑓 (𝑥) 𝑑𝑥 est donnée par ∑𝑛=3
𝑖=0 𝐴𝑖 𝑓 (𝑥𝑖 ) avec 𝐴𝑖 =
𝑥 3
𝐼 (𝐿𝑖 ) = ∫𝑥 3 𝐿𝑖 (𝑥) 𝑑𝑥 où 𝐿𝑖 est le polynôme de Lagrange au point𝑥𝑖 . Après calcul, il vient :𝐴0 = 8 ℎ; 𝐴1 =
0
9 9 3
ℎ; 𝐴2 = 8 ℎ; 𝐴3 = 8 ℎ.
8
Et alors
𝑥3 3
∫𝑥0 𝑓 (𝑥) 𝑑𝑥 ≈ 8 ℎ[𝑓(𝑎) + 3𝑓 (𝑏) + 3𝑓 (𝑐 ) + 𝑓 (𝑑 )]
Ou encore
𝑏
3
∫ 𝑓 (𝑥) 𝑑𝑥 ≈ ℎ[𝑓(𝑥0 ) + 3𝑓 (𝑥𝑛 ) + 3(𝑓 (𝑥1 ) + 𝑓 (𝑥3 )]
𝑎 8
2𝑎+𝑏 𝑎+2𝑏
Où 𝑐 = 𝑒𝑡 𝑑 =
3 3
Théorème : Soit f une fonction n fois continument dérivable sur [a,b]telle que 𝑓 𝑛+1 existe sur ]a,b[.
Si les valeurs de f aux points 𝑥0 , 𝑥1 , … , 𝑥𝑛 appartenant à [a ;b] sont connues, et ∑𝑛𝑖=0 𝐴𝑖 𝑓 (𝑥𝑖 ) est
𝑏
l’approximation d’ordre n de ∫𝑎 𝑓 (𝑥) 𝑑𝑥 alors,
𝑏 𝑛 𝑏 𝑛
𝑀𝑛+1
|∫ 𝑓 (𝑥 )𝑑𝑥 − ∑ 𝐴𝑖 𝑓 (𝑥𝑖 )| ≤ ∫ |∏(𝑥 − 𝑥𝑖 )| 𝑑𝑥
𝑎 (𝑛 + 1)! 𝑎
𝑖=0 𝑖=0
Formule générale :
𝑏 𝑛−1
𝑏 − 𝑎 𝑓 (𝑎 ) + 𝑓 ( 𝑏 ) (𝑏 − 𝑎 ) 3
|∫ 𝑓 (𝑥 )𝑑𝑥 − [ + ∑ 𝑓 (𝑥𝑖 )]| ≤ 𝑀2
𝑎 𝑛 2 12𝑛2
𝑖=1
Il vient pour appliquer la méthode des trapèzes au calcule d’une intégrale avec une erreur inférieur à 𝜀, on
(𝑏−𝑎)3
choisit de n de trapèze grâce à la majoration :𝑀2 ≤𝜀
12𝑛2
Remarque 1:
• La formule du trapèze a un degré de précision de 1 (l’erreur sera nulle pour les polynômes de degré 0 et 1).
Et pour déterminer le nombre (suffisant) de sous intervalles partiels de l’intervalle d’intégration tel que
(𝑏−𝑎)5
l’erreur soit inférieur à 𝜀, il suffit de faire la majoration :𝑀4 ≤𝜀.
180𝑛4
Remarque 2:
-Cette méthode a un degré de précision de 3.Si nous doublons le nombre de sous intervalles alors
- Simpson 1/3 composée revient à intégrer une interpolation quadratique par morceaux de la fonction .
- La formule est d’ordre 4.
Erreur dans la formule des Simpson 3/8
𝑏
3ℎ 3 (𝑏 − 𝑎 ) 5
|∫ 𝑓 (𝑥 )𝑑𝑥 − (𝑓(𝑎) + 3𝑓 (𝑥1 ) + 3𝑓 (𝑥2 ) + 𝑓(𝑏))| ≤ 𝑀4
𝑎 8 2560
Et pour déterminer le nombre (suffisant) de sous intervalles partiels de l’intervalle d’intégration tel que
(𝑏−𝑎)5
l’erreur soit inférieur à 𝜀, il suffit de faire la majoration :𝑀4 ≤𝜀.
80𝑛4
Remarque : Cette méthode aura le même degré de précision et le même ordre que Simpson 1/3
Exemple :
1
2 𝑎+𝑏 1 1
𝐼 = ∫ 𝑒 −𝑥 𝑑𝑥 ; 𝑎 = 0 ; 𝑏 = 1; 𝑐 = = ; 𝑓 (0) = 1; 𝑓 (1) = 0.36788; 𝑓 ( ) = 0.7788
2 2 2
0
1- Trapèze : ℎ = 𝑏 − 𝑎 = 1
6
ℎ 1
𝐼 (𝑓 ) = ∫ 𝑓 (𝑥 )𝑑𝑥 = (𝑓(𝑎) + 𝑓(𝑏)) = (1 + 0.36788) = 6480.6839
0 2 2
(𝑏−𝑎)
Simpson 1/3(n = 2) h= 2 = 3 ,
(𝑏 − 𝑎 ) 1
𝐼 (𝑓 ) = (𝑓(𝑎) + 4𝑓 (𝑐 ) + 𝑓 (𝑏)) (0 + 4 ∗ 0.7788 + 0.36788) = 0.74718
6 6
4- Changement de variable affine
Souvent on définit d’abord une formule de quadrature sur l’intervalle [0;1] ou sur l’intervalle [-1;1] et puis
on la généralise à l’intervalle [xi ; xi+1 ] par un changement de variable affine.
Soit 𝑥𝜖[a ; b] et soit𝑦𝜖 [𝑐; 𝑑], on cherche une transformation 𝑦 = 𝑔 (𝑥) qui envoie l’intervalle [a;b] dans
l’intervalle [c;d] ainsi
𝑑 𝑏
∫𝑐 𝑔 (𝑦)𝑑𝑦 = ∫𝑎 𝑓(𝑔(𝑥 ))𝑔′(𝑥)𝑑𝑥
Si 𝑔′(𝑥) est une constante, i.e. si g est une transformation affine𝑔(𝑥 ) = 𝑚𝑥 + 𝑞, alors
𝑑 𝑏
∫ 𝑔 (𝑦)𝑑𝑦 = 𝑚 ∫ 𝑓 (𝑚𝑥 + 𝑞 )𝑑𝑥
𝑐 𝑎
On obtient
𝑑−𝑐 𝑐𝑏 − 𝑎𝑑
𝑚= ;𝑞=
𝑏−𝑎 𝑏−𝑎
Par conséquent
𝑑−𝑐 𝑐𝑏−𝑎𝑑
𝑦 = 𝑏−𝑎 𝑥 + ,d’où
𝑏−𝑎
𝑑
𝑑−𝑐 𝑏 𝑑−𝑐 𝑐𝑏 − 𝑎𝑑
∫ 𝑓 (𝑦)𝑑𝑦 = ∫ 𝑓( 𝑥+ ) 𝑑𝑥
𝑐 𝑏−𝑎 𝑎 𝑏−𝑎 𝑏−𝑎
Exemple
Transformer l’intervalle [0;1] dans l’intervalle [𝑥𝑖 ; 𝑥𝑖+1 ] par un changement de variable affine.
On a𝑦 = (𝑥𝑖+1 − 𝑥𝑖 )𝑥 + 𝑥𝑖 et
𝑥𝑖+1 1
∫ 𝑓(𝑦)𝑑𝑦 = (𝑥𝑖+1 − 𝑥𝑖 ) ∫ 𝑓((𝑥𝑖+1 − 𝑥𝑖 )𝑥 + 𝑥𝑖 )𝑑𝑥
𝑥𝑖 0
Transformer l’intervalle [-1;1] dans l’intervalle [𝑥𝑖 ; 𝑥𝑖+1 ] par un changement de lorsque x=1/2 alors 𝑦 =
𝑥𝑖+1 −𝑥𝑖
etc.variable affine
2
𝑥𝑖+1 1
𝑥𝑖+1 − 𝑥𝑖 𝑥𝑖+1 − 𝑥𝑖
∫ 𝑓 (𝑦)𝑑𝑦 = ( ) ∫ 𝑓 (𝑥𝑖 + (1 + 𝑥)) 𝑑𝑥
𝑥𝑖 2 0 2
5- Formules de Gauss
On se pose la question suivante : comment choisir au mieux les points d'intégration 𝑥𝑖 pour que la formule de
quadrature soit de degré de précision le plus élevé possible ? Le problème revient donc à trouver à la fois les
(𝑛)
poids 𝛼𝑖 , 𝑖 = 0, . . . , 𝑛 et les points 𝑥𝑖 , 𝑖 = 0, . . . , 𝑛 de la formule de quadrature :
𝑏 𝑛
(𝑛)
∫ 𝑓(𝑥 )𝑑𝑥 ≃ ∑ 𝛼𝑖 𝑓(𝑥𝑖 )
𝑎 𝑖=0
On a donc 2n + 2 inconnues à déterminer ! On cherche alors une formule de quadrature exacte sur 𝑃2𝑛+1 , ce
qui donne les 2n + 2 équations :
𝑏 𝑛
(𝑛)
∫ 𝑥 𝑘 𝑑𝑥 = ∑ 𝛼𝑖 𝑥𝑖𝑘 ; 𝑘 = 0,1, … ,2𝑛 + 1
𝑎 𝑖=0
Remarque 1. Le degré de précision de la formule proposée peut-il être 2n + 2? La réponse est non. En effet
; soit le polynôme de degré 2n + 2,
𝑄(𝑥) = (𝑥 − 𝑥0 )2 (𝑥 − 𝑥1 )2 . . . (𝑥 − 𝑥𝑛 )2 .
𝑏
Ce polynôme est strictement positif (𝑠𝑖 𝑥 ≠ 𝑥𝑖 , 𝑖 = 0, . . . , 𝑛) donc ∫𝑎 𝑄(𝑥) 𝑑𝑥 > 0 alors que
∑𝑛𝑖=0 𝛼𝑖(𝑛) 𝑄(𝑥𝑖 ) = 0
Théorème 1. Les formules de type Gauss3 sont stables et convergentes pour toute fonction continue.
Remarque 2
1. Stable signifie les poids intervenant dans la formule de quadrature sont strictement positifs.
(𝑛) 𝑏
2. Convergente signifie lim ∑𝑛𝑖=0 𝛼𝑖 𝑓(𝑥𝑖 ) = ∫𝑎 𝑓(𝑥) 𝑑𝑥
𝑛→∞
3. Les formules de type Gauss sont souvent utiles pour les intégrales impropres convergentes.
Déterminons les poids 𝛼0 , 𝛼1 et les points d'intégration 𝑥0 et 𝑥1 pour que la formule proposée soit exacte sur
P3 . Les quatre équations sont les suivantes :
1
𝑓 ≡ 1, R(f) = 0 ⇔ ∫ 1 𝑑𝑥 = 2 = 𝛼0 + 𝛼1
−1
1
𝑓(𝑥 ) = 𝑥, R(f) = 0 ⇔ ∫ 𝑥 𝑑𝑥 = 0 = 𝛼0 𝑥0 + 𝛼1 𝑥1
−1
1
2
𝑓(𝑥 ) = 𝑥 , R(f) = 0 ⇔ ∫ 𝑥 2 𝑑𝑥 =
2
= 𝛼0 𝑥02 + 𝛼1 𝑥12
−1 3
1
𝑓(𝑥) = 𝑥 3 , R(f) = 0 ⇔ ∫ 𝑥 3 𝑑𝑥 = 0 = 𝛼0 𝑥03 + 𝛼1 𝑥13
−1
D'où,
𝛼0 𝑥0 = −𝛼1 𝑥1
{
𝛼0 𝑥1 = (𝑥12 − 𝑥02 )
De la dernière équation on tire :
2
Soit 𝛼1 = 0, 𝛼0 = 0𝑒𝑡 𝑥0 = 0 ce qui est impossible puisque 𝛼0 𝑥02 + 𝛼1 𝑥12 = 3.
1
1 5 1 5 1 1
∫ 𝑓(𝑥) 𝑑𝑥 = 𝑓 (−1) + 𝑓 (− ) + 𝑓 ( ) + 𝑓(1) + 𝑅(𝑓 )
−1 6 6 √5 6 √5 6
Exercice 1
1
On veut calculer l'intégrale∫0 𝑐𝑜𝑠(𝑥 2 )𝑑𝑥 dx avec une précision de 10−8. Déterminer le nombre
de points nécessaires pour l'évaluer avec cette précision par la méthode des trapèzes et par
la méthode de Simpson.