TP3 Methode Integration
TP3 Methode Integration
TP3 Methode Integration
TP N°3
Intégration numérique
I. Introduction
Lorsque le calcul analytique est difficile (voir impossible) le calcul numérique de l’intégral d’une fonction
s’impose. L’une des manières les plus naturelles pour calculer l’aire sous une courbe est de l’approximer par
une aire facilement calculable.
(b − a )
b
Pour évaluer l’intégrale ∫ f ( x)dx
a
on divise l’intervalle a ≤ x ≤ b en (n) sous intervalles de longueur ∆x =
n
b x1 x2 b n xk
alors : ∫ f ( x)dx = ∫ f ( x)dx + ∫ f ( x)dx + ... + ∫ f ( x)dx = ∑ ∫ f ( x)dx
a a x1 xn−1 k =1 xk −1
L’intégrale de chaque sous-intervalle peut être approximé par la surface sous chaque intervalle.
(b − a )
∆x = ;
n
Pour k = 1 à (n − 1)
x = a + ∆x × k ;
s = s + f ( x);
fin pour
( f (a) + f (b))
T = ∆x × + ∆x × s;
2
III. 2 Application
1. Implémenter la formule d’intégration numérique de Trapèze sous forme d’une fonction Matlab:
function T=MethTrap(f, a, b, N)
2. Appliquer la fonction produite pour trouver l’intégrale de la fonction suivante avec N=30 :
2
∫ x + x cos( x ) dx
1
3. Afficher l’erreur d’approximation après le calcule de la valeur exacte.
function S=MethSimp(f, a, b, N)
Avec les arguments d’entrée :
a. f : la fonction à Intégrer (à entrer comme un string - entre guillemet ‘ ’ ).
b. a,b : l’intervalle d’intégration.
c. N : nombre de division.
Et les arguments de sortie :
d. S : la valeur approximative de l’intégrale.
2. Appliquer la fonction produite pour trouver l’intégrale de la fonction suivante avec N=20 :
4
∫ xe
x
dx
2
3. Afficher l’erreur d’approximation après le calcule de la valeur exacte.