2eme Année physique Méthodes numériques et programmation
Série TP N=4 (Solution)
Intégration numérique
1 Exercice
Soient deux fonction : f (x) = sin(x) et g(x) = x
1+x2
Utilisant la méthode de point milieu donnez une valeur numérique approchée de :
Z π
I1 = f (x)dx
−π
Z 3
I2 = g(x)dx
0
Solution 1
1. f (x) = sin(x)
b+a π−π
Ipm 1 = (b − a)f ( ) = (π − (−π))sin( ) = 2πsin(0) = 0
2 2
Alors que : Z 0 Z π
I1 = f (x)dx + f (x)dx
−π 0
1 0−π
Ipm = (0 − (−π))sin( ) = −π
2
2 π+0
Ipm = (π − 0)sin(( )) = π
2
Alors :
1 2
I1 = |Ipm | + |Ipm | = 2π
Remarque
1. Si la valeur de la fonction est négative dans l'interval d'intégration. L'intégral est l'opposé de
l'aire calculer pour la surface entre l'axe des abcisses et la courbe de la fonction.
2. Si la fonction prend des valeurs positives et négatives alternativement, on calacule séparément
l'intégrale de la fonction dans l'intervalle positive et negative, et on fait la somme avec la valeur
absolue.
2. g(x) = x
1+x2
3+0 6 18
Ipm 2 = (3 − 0)g( )=3× = = 1.3846
2 13 13
2 Exercice
π
Déterminer par la méthode des trapèzes puis par celle de Simpson : f (x)dx, sur la base du
R 2
0
tableau suivant :
π π 3π π
x 0 8 4 8 2
f (x) 0 0.382683 0.707107 0.923880 1
Ces points d'appui sont ceux donnant sin(x), comparer alors les résultats obtenus avec la valeur exacte.
http://rahab.e-monsite.com 1 2016-2017
2eme Année physique Méthodes numériques et programmation
Solution 2
Z π
2
f (x)dx
0
1- la méthode des trapèzes :
n = 4 le pas h = xn −x n
0
= π8
h P4
It = 2 k=1 (f (xk−1 + f (xk ))
= h2 [(f (x0 ) + f (x1 )) + (f (x1 ) + f (x2 )) + (f (x2 ) + f (x3 )) + (f (x3 ) + f (x4 ))]
π
= 16 [(0 + 0.382683) + (0.382683 + 0.707107) + (0.707107 + 0.923880) + (0.923880 + 1)]
= 0.987116
Le programme Matlab :
>> I=trapeze(0,pi/2,4)
I =
0.98711590069363
xn −x0
2- La méthode de Simpson : n = 2, Le pas h = n = π
4
4
hX π
Is (f ) = (f (xk−1 ) + 4f (x¯k ) + (f (xk )) = [(0+4(0.382683)+0.707107)+(0.707107+4(0.923880)+1)]
6 k=1 24
Is = 1.000135
Le programme Matlab :
>> I=Simpson(0,pi/2,4)
I =
1.00013466065011
π
Et on a 02 sin(x)dx = 1, alors :
R
|I − Is | = |1 − 1.000135| = 0.000135 et |I − It | = |1 − 0.987116| = 0.012884
On constate donc que l'approximation de I donnée par la méthode de Simpson est meilleure que celle
par les trapèzes.
3 Exercice
On lance une fusée verticalement du sol et l'on mesure pendant les premières 80 secondes l'accélé-
ration γ :
t(ens) 0 10 20 30 40 50 60 70 80
γ(enm/s2 ) 30 31.63 33.44 35.47 37.75 40.33 43.29 46.70 50.67
Calcule la vitesse V de la fusée à l'instant t = 80s, par la méthode des trapçèes puis par celle Simpson.
Solution 3
On sait que l'acceleration γ est la dérivée de la vitesse V , donc :
Z 80
V (t = 80s) = γ(t)dt
0
http://rahab.e-monsite.com 2 2016-2017
2eme Année physique Méthodes numériques et programmation
A. Méthode des trapèzes
1. Méthode de trapèze de Matlab :
>> x=[0 10 20 30 40 50 60 70 80]
x =
0 10 20 30 40 50 60 70 80
>> y=[30 31.63 33.44 35.47 37.75 40.33 43.29 46.70 50.67]
y =
30.0000 31.6300 33.4400 35.4700 37.7500 40.3300 43.2900 46.7000 50.6700
>> I=trapz(x,y)
I =
3.0895e+003
2. Programme Matlab de notre cours :
function I=trapeze(a,b,n)
h=(b-a)/n;
I=0;
x=[a:h:b];
F= [30 31.63 33.44 35.47 37.75 40.33 43.29 46.70 50.67];
for i=1:n
I = I+h*(F(i)+F(i+1))/2;
end
On a n = 8, l'appel du programme sera ainsi :
>> I=trapeze(0,80,8)
I =
3.0895e+003
- La vitesse du fusée estimée par la méthode de trapèse est alors de : 3089.5 m/s.
B. Méthode de Simpson
>> I=Simpson(0,80,8)
I =
3.087166666666667e+003
- La vitesse du fusée estimée par la méthode de simpson est alors de : 3087.16 m/s.
4 Exercice
Rπ
Calculer à l'aide de la méthode des trapèzes l'intégrale I = 0 sin x2 dx avec un nombre de points
d'appui n = 5 puis n = 10.
http://rahab.e-monsite.com 3 2016-2017
2eme Année physique Méthodes numériques et programmation
Solution 4
1. Avec n = 5 :
function I=trapeze(a,b,n)
h=(b-a)/n;
I=0;
x=[a:h:b];
F=sin(x.^2);
for i=1:n
I = I+h*(F(i)+F(i+1))/2;
end
et l'exécution :
>> I=trapeze(0,pi,5)
I =
0.50443075989740
>> I=trapeze(0,pi,10)
I =
0.72238098855789
- Alors que la valeur 'exacte' est approximativement 0, 772651. Avec ce pas plus petit l'approxima-
tion numérique est meilleure.
5 Exercice
Approximer par la formule de point milieu décomposant l'intervalle d'intégration en dix parties,
l'intégrale : Z 2
√
xdx
1
Solution
On a : n = 10 alors : h = b−a
n = 2−1
10 = 0.1
Le vecteur x est :
x = [1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2]
Le vecteur x̄ :
xb = [1.05 1.15 1.25 1.35 1.45 1.55 1.65 1.75 1.85 1.95]
L'inégrale
√: √ √ √ √ √ √ √ √ √
I = 0.1( 1.05 + 1.15 + 1.25 + 1.35 + 1.45 + 1.55 + 1.65 + 1.75 + 1.85 + 1.95) = 1.2190
Le programme Matlab :
function I=PointMilieuComposite(a,b,n)
h=(b-a)/n;
I=0;
x=[a:h:b];
for j=1:n
http://rahab.e-monsite.com 4 2016-2017
2eme Année physique Méthodes numériques et programmation
xb=(x(j)+x(j+1))/2;
I=I+h*sqrt(xb);
end
L'exécution :
>> I=PointMilieuComposite(1,2,10)
I =
1.2190
6 Exercice
Évaluer à l'aide de la méthode de Simpson, avec 20 subdivision de l'intervalle d'intégration l'inté-
grale : Z π
cos(x)dx
−π
Solution
Si on fait l'intégrale sur l'intervalle [−π, π], on obtient I = 1.3951e − 016 ' 0 :
>> I=Simpson(-pi,pi,20)
I =
1.3951e-016
On sait que la fonction cos(x) s'annulle à −π
2 et π
2 et change de signe alors :
Z π Z −π Z π Z π
2 2
cos(x)dx = | cos(x)dx| + | −π
cos(x)dx| + | cos(x)dx|
π
−π −π 2 2
on obtient alors :
Figure 1 La fonction f (x) = cos(x) dans [−π, π]
http://rahab.e-monsite.com 5 2016-2017
2eme Année physique Méthodes numériques et programmation
>> I=Simpson(-pi,-pi/2,20)
I =
-1.0000
>> I=Simpson(-pi/2,pi/2,20)
I =
2.0000
>> I=Simpson(pi/2,pi,20)
I =
-1.0000
Donc :
I = | − 1| + |2| + | − 1| = 4
http://rahab.e-monsite.com 6 2016-2017
