Université de Savoie L3 MASS

Méthodes numériques TD 3

Exercice 1 : On veut approcher l’intégrale

Z 1/2
I= exp(−x2 ) dx
0

1. Donner une valeur approchée de I en utilisant la méthode

• rectangle à gauche
• rectangle à droite
• point milieu
• trapèze
• Simpson
2. Utiliser les méthodes des rectangles à gauche, des trapèzes et de Simpson composites avec 3 sous-
intervalles pour approcher I.
3. Donner dans les trois cas un majorant de l’erreur.
4. En combien de sous-intervalles faut-il découper [0, 1/2] dans les méthodes des rectangles à gauche
composite et des trapèzes composite pour avoir un majorant du même ordre de grandeur que celui
obtenu avec la méthode de Simpson pour trois sous-intervalles.

Exercice 2 : On considère l’intégrale Z 3

I= log(x) dx
1
1. Calculer la valeur exacte de I
2. Donner une valeur approchée de I en utilisant la méthode des trapèzes composite avec 4 sous-
intervalles. Comparer avec la valeur exacte.
3. Combien faut-il prendre de sous-intervalles pour assurer une erreur inférieure à 10−3 ?

Exercice 3 : Soit f une fonction C ∞ (R, R).

1. On considère l’approximation
Z 1     
2 1 1
f (x) dx ' 2f − − f (0) + 2f
−1 3 2 2
Donner le degré de précision de cette formule de quadrature.

1
 b−a
2. On subdivise l’intervalle [a, b] en J sous-intervalles [aj , aj+1 ] avec aj = a + jh et h = . À l’aide
J
d’un changement de variable en déduire une formule de quadrature pour approcher l’intégrale
Z aj+1
f (x) dx
aj

Rb
puis une formule de quadrature composite pour l’intégrale a f (x) dx.

Exercice 4 :
R1
1. Soit f : [0, 1] → R de classe C 1 , on pose I[0,1] (f ) = 0 f (t) dt et
         
˜ 1 2 0 1 0 2
I(f ) = α f +f +β f −f
3 3 3 3

Déterminer α et β pour que l’approximation de I par I˜ soit de degré de précision le plus élevé
possible et préciser alors le degré.

2. On pose
f 0 (0)
 
˜
I(f ) = α f (0) + + βf (c)
2
Déterminer α, β et c ∈ [0, 1] pour que l’approximation de I par I˜ soit de degré de précision le plus
élevé possible et préciser alors le degré.

3. Déterminer des constantes α, β telles que l’approximation

Z 1
xf (x) dx ' αf (0) + βf (1)
0

soit exacte pour f polynôme de degré aussi élevé que possible.

2
 Corrigés

Exercice 1 :
2
1. On pose f (x) = e−x

• rectangle à gauche
f (0) 1
I˜ = =
2 2
• rectangle à droite
f (1/2)
I˜ = ≈ 0, 3894
2
• point milieu
f (1/4)
I˜ = ≈ 0, 469707
2
• trapèze   
1 1 1
I˜ = × f (0) + f ≈ 0, 4447
2 2 2
• Simpson     
1 1 1
I˜ = f (0) + 4f +f ≈ 0, 461371
6×2 4 2
1
2. On prend J = 3 sous-intervalles sur [0, 1/2], ce qui revient à poser h = et les points
6
1 1 1
a0 = 0, a1 = , a2 = , a3 =
6 3 2
• rectangles à gauche
3
X
I˜3 = h f (aj−1 )
j=1
    
1 1 1
= f (0) + f +f
6 6 3
≈ 0, 477907

• trapèzes
3
X f (aj−1 ) + f (aj )
I˜3 = h
2
j=1
 
f (a0 ) + f (a3 )
= h f (a1 ) + f (a2 ) +
2
≈ 0, 459474

3
 • Simpson composite
           
˜ h f (0) + f (1/2) 1 1 1 1 5
I3 = +f +f +2 f +f +f
3 2 6 3 12 4 12
≈ 0, 461282

3. Grâce un logiciel de calcul formel (ou la calculatrice) on a

2
f (x) = e−x
2
f 0 (x) = −2xe−x
2
f 00 (x) = (4x2 − 2)e−x
 
3 −x2
f (3) (x) = −8x x2 − e
2
2
f (4) (x) = (16x4 − 4x2 + 12)e−x

et les majorants
  
1 
(4) (2)
M1 = f 0

M4 = |f (0)| = 12, M2 = |f (0)| = 2, ≈ 0, 778801
2 

Ceci nous permet de majorer l’erreur commise avec les trois méthodes composites

• rectangles à gauche
(b − a)2
E3 (f ) ≤ M1 ≈ 0, 03245
2J
• trapèzes
(b − a)3
E3 (f ) ≤ M2 ≈ 0, 0023
12J 2
• Simpson composite
(b − a)5
E3 (f ) ≤ M4 ≈ 1, 6.10−6
2880J 4
4. Soit JR le nombre nécessaire de sous-intervalles nécessaire pour obtenir une précision de 10−6 avec
la méthode des rectangles, JT pour la méthode des trapèzes.

• pour la méthode des rectangles on doit assurer

(b − a)2
M1 ≤ 10−6
2JR

(b − a)2
=⇒ JR ≥ M1 106
2
on trouve JR ≥ 97351.
• pour la méthode des trapèzes
JT ≥ 145

Exercice 2 :

4
 1. Une intégration par parties donne
Z 3 Z 3
3
I= 1 × log(x) dx = [x log x]1 − 1 dx = 3 log 3 − 2 ≈ 1, 2958
1 1

2. la fonction x 7→ log x est concave donc au-dessus de toutes ses cordes. On en déduit que la valeur
approchée de l’intégrale par la méthode des trapèzes sera toujours inférieure, quelque soit le nombre
de sous-intervalles, à la valeur exacte de I calculée à la question précédente.
Si on divise l’intervalle [1, 3] en J = 4 sous-intervalles cela revient à prendre h = 1/2 et

a0 = 1, a1 = 3/2, a2 = 2, a3 = 5/2, a4 = 3

La méthode des trapèzes composite donne

    
˜ f (1) + f (3) 3 5
I(f ) = h +h f + f (2) + f ≈ 1, 2821.
2 2 2

3. Pour obtenir un précision de 10−3 on doit diviser l’intervalle en JT sous-intervalles avec

(b − a)3
M2 < 10−3 =⇒ J ≥ 26
12J 2
où on a calculé M2 = 1.

Exercice 3 :
1. La méthode est d’ordre de précision égal à 3.

2. On introduit le changement de variable ϕj qui envoie [−1, 1] sur [aj , aj+1 ]

aj + aj+1 aj+1 − aj
ϕj (t) = +t
2 2
En faisant ce changement de variable sous l’intégrale on obtient
Z aj+1
aj+1 − aj 1
Z
f (x) dx = f ◦ ϕj (t) dt
aj 2 −1
aj+1 − aj
' × 23 (2f (φ(−1/2)) − f (φ(0)) + 2f (φ(1/2)))
2
On en déduit la formule de quadrature pour l’intervalle [aj , aj+1 ]
Z aj+1       
aj+1 − aj aj+1 + 3aj aj+1 + aj 3aj+1 + aj
f (x) dx ' 2f −f + 2f
aj 3 4 2 4
      
aj+1 − aj h h 3h
= 2f aj + − f aj + + 2f aj +
3 4 2 4
Enfin en sommant on revient à l’intégrale sur l’intervalle [a, b]
Z b J−1       
hX h h 3h
f (x) dx ' 2f aj + − f aj + + 2f aj +
a 3 4 2 4
j=0

5
Exercice 4 :

1. Pour déterminer quel ordre de précision maximal on peut atteindre avec cette méthode de quadra-
ture, on teste la méthode sur les monômes xi pour en déduire des conditions sur les coefficients α,
β

• on teste la méthode contre f (x) = 1, la méthode est exacte pour f (et donc par linéarité pour
les constantes) ssi
1
1 = 2α ⇐⇒ α =
2
la méthode est donc au moins d’ordre 0 si α = 1/2. Dans la suite on suppose que α = 1/2
• on teste le monôme f (x) = x
 
1 1 2
=α + + β(0 − 0)
2 3 3

cette égalité n’introduit pas de condition supplémentaire sur les coefficients, la méthode est
donc au moins d’ordre 1
• on teste f (x) = x2 , la méthode est d’ordre au moins 2 si la relation suivante est satisfaite
   
1 1 1 4 2 4 1
= + +β − ⇐⇒ β = −
3 2 9 9 3 3 12

• avec α = 1/2 et β = −1/12 on vérifie que la méthode est exacte pour le monôme x3 et donc
que la méthode est d’ordre au moins 3
• on teste maintenant f (x) = x4
   
1 1 1 16 4 1 8 31
6= + − − =
5 2 81 81 12 27 27 162

=⇒ si α 6= 1/2 alors la méthode n’est même pas d’ordre 0, si α = 1/2 et β 6= −1/12 alors la
méthode est d’ordre 1, enfin si α = 1/2, β = −1/12 alors la méthode est d’ordre 3.

2. on reprend la même démarche pour trouver des relations entre α, β et c

• de degré au moins 0 ssi α + β = 1

• de degré au moins 1 ssi
α 1
+ βc =
2 2
• de degré au moins 2 ssi
1
βc2 =
3
la résolution de ce système (attention non-linéaire !) donne

1
c=


2



4

β=
 3
1


 α=−


3

6
 On en déduit donc que c’est pour ces valeurs de α, β, c que l’ordre de la méthode est maximal. Il
est de plus exactement égale à 2 car en testant f (x) = x3 , on a
Z 1    3
3 1 1 1 1 4 1
x dx = =6 =− 0+ ×0 + ×
0 4 6 3 2 3 2

3. idem, on trouve que la méthode est d’ordre de précision maximal pour α = 1/6, β = 1/3. L’ordre
est alors exactement égal à 1.