Université Abou Bekr Belkaïd – Tlemcen
Faculté de Technologie
Département de Télécommunication
Licence 2 Année universitaire 2019/2020
TD n°2 : Intégration numérique
Exercice 1 : Question de cours Du chapitre I au chapitre II
1. Soit f une fonction intégrable sur le segment [a;b] (a < b). Retrouver la formule des trapèzes
T1 en faisant une interpolation (de f ) à 2 points.
2. Retrouver la formule de Simpson en faisant une interpolation en 3 points dont les abscisses
a+b
sont respectivement a, b et .
2
Exercice 2 :
On lance une fusée verticalement du sol et l’on mesure pendant les premières 80 secondes
l’accélération γ :
Calculer la vitesse V de la fusée à l’instant t = 80s, par la méthode des rectangles, des trapèzes
puis par Simpson.
Exercice 3 :
Une tige métallique (une antenne) se déforme lorsqu'elle est soumise à une charge (un poids).
Les mesures expérimentales ont donné les résultats pour la force de torsion s (qu'on appelle la
contrainte) en fonction du taux de déformation e
Les mesures sont prises jusqu’à la rupture de la tige (dans ce cas, e = 0,25).
1. L’aire sous la courbe de la fonction s(e) représente l’énergie nécessaire pour provoquer la
0,25
rupture de l'antenne et elle s'écrit : ∫ s( e)de
0
Calculer cette aire par la méthode de votre choix.
Exercice 4 :
Trouver le nombre n de subdivisions nécessaires de l’intervalle d’intégration [−π, π], pour
π
−3
évaluer à 0.5 10 prés, grâce à la méthode de Simpson, l’intégrale ∫ cos(x )dx
−π
Même question avec la méthode des trapèzes. Comparer les résultats.
Exercice 5 :
L’objectif de cet exercice pratique est d’implémenter et de comparer deux méthodes
b
d'intégration numérique pour calculer I = ∫ f (x)dx
a
1. Écrire une fonction Matlab/Scilab d’entête
function [S] = Int_Simpson(a,b,f,n)
qui utilise la méthode de Simpson avec n intervalles de même longueur pour obtenir une
approximation S de I.
2
2. Tester la fonction Int_Simpson sur l’exemple suivant : I = ∫ e x dx =e 2
−1.
0
Vous devez observer la convergence vers e 2
− 1 lorsque n → ∞.
3. Écrire une deuxième fonction Matlab/Scilab d’entête
function [T] = Int_Trapezes(a,b,f,n)
qui utilise la méthode des trapèzes avec n intervalles de même longueur pour obtenir une
approximation T de I.
2
4. Tester la fonction Int_Trapezes sur l’exemple : I = ∫ e x dx =e 2
−1.
0
2
Vous devez observer la convergence vers e − 1 lorsque n → ∞.
5. Comparer la vitesse de convergence des deux méthodes.
Culture G. : Les noms qui ont marqués l'histoire des intégrales sont : Leibniz ,
Riemann et Lebesgue ! Pères de théorie à l'origine du signal d'aujourd'hui !