TP méthode numérique TP3 :L’intégration numérique
But de TP
1- Utiliser les fonctions de manipulation 1 et 2 ci-dessous pour calculassions
Calculez les approximations de l’intégrale :
π
𝐼(𝑓) = ∫ sin (x2) 𝑑𝑥
0
en utilisant les méthodes du point milieu, du trapèze et de Simpson..
-calculer l’intégrale, pour n =5 et 10 sous-intervalles.
-Sachant que la valeur exacte est 0,7726
- comparer le résultat obtenu avec la valeur exacte
2- utiliser le programme de manipulation 3 pour écrire le code de méthode
de la fonction Simpson
MÉTHODE DES RECTANGLES
Dans cette méthode, on calcule l’intégrale numérique en réalisant une somme de
surfaces de rectangles. Nous subdivisons l'intervalle [a,b]en n sous-intervalles dont
les bornes sont xi. Les longueurs de ces sous intervalles sont .
hi =xi+1 - xi Nous construisons les rectangles dont les côtés sont hi et f(xi).
b n
∫ f ( x ) dx ≈ h ∑ f ( xi)
a i=0
1.2 Méthode du point milieu
L’idée de base de cette méthode, consiste à subdiviser l’intervalle [𝑎, 𝑏] en 𝑛 sous-
intervalles [𝑥𝑘, 𝑥𝑘+1]. Dans le cas où les sous-intervalles sont équidistants, on
écrira h = (𝑏−𝑎)/𝑛.
[Tapez un texte] Page 1
TP méthode numérique TP3 :L’intégration numérique
Ainsi, le schéma numérique de cette méthode s’écrira comme :
b
∫ f ( x ) dx ≈ h ∑ f ( xi) xi=( a+ h2 )+i∗h
n
a i=0
Manipulation .1
Ouvrir un fichier script et écrire le programme suivant de la fonction Matlab de
méthode point milieu avec le nome point.m
function result_integration = midpoint(f, a, b, n)
h = (b-a)/n;
result = 0;
for i = 0:(n-1)
result = result + f((a + h/2) + i*h);
end
result_integration = h*result;
end
2-Méthode de trapèze
Pour cette règle, nous ferons la même mise en place que pour la règle du milieu.
b−a
Nous allons décomposer l'intervalle [a,b] en n sous-intervalles de largeur, h= n .
Ensuite, sur chaque sous-intervalle, nous allons approximer la fonction avec une
ligne droite qui est égale aux valeurs de la fonction à chaque extrémité de
l'intervalle.
Chacun de ces morceaux est un trapèze et comme nous pouvons le constater, cette
méthode peut donnée une très bonne valeur d'approximation de l'integration
[Tapez un texte] Page 2
TP méthode numérique TP3 :L’intégration numérique
La zone du trapèze dans l'intervalle [ xi-1,xi]est donnée par,
En faisant une petite simplification, nous arrivons à la règle générale des trapèzes.
b
∫ f ( x ) dx=h¿]
a
Manipulation 2
1- Ouvrir un fichier script et écrire le programme dessous de la fonction Matlab
de méthode trapèze avec le nome trapez.m
function integral = trapez(f, a, b, n)
h = (b-a)/n;
result = 0.5*f(a) + 0.5*f(b);
for i = 1:(n-1)
result = result + f(a + i*h);
end
integral = h*result;
end
3 Méthode de Simpson
C'est la dernière méthode que nous allons examiner et dans ce cas, nous diviserons à
nouveau l'intervalle [a,b] en n sous-intervalles. Cependant, contrairement aux deux
méthodes précédentes, nous devons exiger que n soit pair. La raison de ceci sera
évidente dans un peu. La largeur de chaque sous-intervalle est,
[Tapez un texte] Page 3
TP méthode numérique TP3 :L’intégration numérique
Pour la règle de Simpson, nous allons approximer la fonction avec un quadratique et
nous allons exiger que le quadratique soit d'accord avec trois des points de nos sous-
intervalles
Manipulation 3
1
On se propose de calculer l’intégrale définie f = ∫ sin (x) . e x , Ecrire un rogramme
0
qui calcule cette intégrale en utilisant la méthode du simpson et en divisant
l’intervalle d’intégration en 10
f=inline('sin(x).*exp(x)');
a=0;
b=1;
n=10;
h=(b‐a)/n;
x=[a:h:b];
y=f(x);
sigma1=0;
k=n/2;
for i=1:k
sigma1=sigma1+y(2*i‐1);
end
sigma2=0;
for i=1:k
sigma2=sigma2+y(2*i);
end
Is=(h/3)*(y(1)+y(n+1)+2*sigma1+4*sigma2)
% Utilisation de la fonction predefinie Matlab: QUAD
Ism=quad(f,0,1)
[Tapez un texte] Page 4
