Université des Sciences et de la Technologie d’Oran MB
Département Electrotechnique
Licence 2
Niveau : Semestre 4Module : Méthodes Numériques
TP n°2 : Intégration numérique
But du TP : Le but de ce TP est l’implémentation d’un algorithme de calcul d’une intégrale approchée
par les méthodes ( Trapèzes et Simpson d’ordre 4). Puis d’étudier la variation de l’erreur lorsque l’on
augmente le nombre de sous intervalles.
Théorie brève :
x 0− xn
Soit une fonction f à intégrer sur l’intervalle [x0 xn] avec un pas de discrétisationh= et n (le
n
nombre de sous intervalles). On peut estimer cette intégrale par plusieurs méthodes dont celles des
Trapèzes et celles de Simpson.
Méthode des trapèzes pour un pas constant
[ ]
n −1
I tr =h /2∗ f ( x 0 ) + f ( x n ) +2∗∑ f ( x i )
i =1
Méthode de Simpson d’ordre 4
[ ]
n /2−1 n/ 2
I S=h/3∗ f ( x 0 ) + f ( xn ) + 2∗ ∑ f ( x 2 i ) +4∗∑ f ( x i 2−1)
i=1 i=1
Préparation théorique : (obligatoire et individuelle. Doit être remise au début de chaque séance de TP)
Soit l’intégrale à calculer
1 1
I =∫ f ( x ) dx=∫ sin ( x ) cos ( x ) dx
2
−1 −1
1- Calculer I par la méthode analytique
2- En utilisant un nombre de sous intervalles n=10, estimer la valeur approchée de I par la méthode
des trapèzes puis par la méthode de Simpson d’ordre 4
3- Comparer les résultats
4- Elaborer l’algorithme de calcul de I par la méthode des trapèzes et la méthode de Simpson
d’ordre 4 (il est recommandé de vous servir de l’instruction sum) qui :
a. Déclare
i. Le nombre de sous intervalles n=10
ii. La valeur analytique de I
b. Détermine le pas de discrétisation h
c. Remplit un tableau nommé x avec les abscisses des points d’appui
d. Remplit un tableau nommé f avec les ordonnées des points d’appui
e. Calcul l’intégrale de f par la méthode des Trapèzes
f. Calcul l’intégrale de f par la méthode de Simpson d’ordre 4
g. Calcule et affiche les erreurs d’intégration des deux méthodes
5- En déduire le programme Matlab correspondant.
Travail demandé :
1. Ouvrer une fenêtre d’éditeur de script et écrire le programme Matlab préparé dans la
question N°5
2. Exécuter le programme et corriger les erreurs s’il y en a
3. Sur votre compte rendu reporter les résultats :
a. Dans un tableau relever les erreurs de calcul pour chaque méthode
b. Faite varier n=10,20, 30,40 et relever à chaque fois les erreurs de calcul
c. Commenter les résultats