Modélisation et Calcul Scientifique 2 - MT 01
M1 SPMA/TI
Feuille de TD n◦ 1 : Intégration numérique
Exercice 1 : Application de la méthode des trapèzes
Z 5/2
Estimer f (t)dt à partir des données suivantes :
0
x 0 1/2 1 3/2 2 5/2
f (x) 3/2 2 2 1.6364 1.2500 0.9564
en utilisant la méthode des trapèzes composée.
Exercice 2 : Ordre d’une formule d’intégration
Z 1
˜ 4 1 2 4 1
Quel est l’ordre de la formule I(f ) = f − − f (0) + f qui approche f (t)dt ?
3 2 3 3 2 −1
Exercice 3 : Méthodes de Newton-Cotes : poids, ordre et formules d’erreur
On considère les formules obtenues à l’aide des polynômes d’interpolation pour des points également
répartis sur [−1, 1], c’est-à-dire
2
xi = −1 + i pour 0 ≤ i ≤ n.
n
Ce sont les méthodes de Newton-Cotes.
3.1. Montrer qu’il existe des poids λi , 0 ≤ i ≤ n uniques tels que la méthode de Newton-Cotes à
n + 1 points soit d’ordre n.
3.2. Calculer les poids de la méthode de Simpson, qui est la méthode de Newton-Cotes pour n = 2
sur [−1, 1].
3.3. Montrer les formules de symétrie suivantes : xn−i = −xi , Ln−i (x) = Li (−x) et λn−i = λi , pour
tout 0 ≤ i ≤ n.
3.4. En déduire que si n est pair, la méthode de Newton-Cotes à n + 1 points est d’ordre n + 1.
3.5. Calculer les noyaux de Péano sur [−1, 1] pour la méthode du point milieu et la méthode des
trapèzes (n = 1) et vérifier qu’ils sont de signe constant sur [−1, 1]. Montrer que
1
(a) |E(f )| ≤ ||f 00 ||∞ , pour la méthode du point milieu,
3
2
(b) |E(f )| ≤ ||f 00 ||∞ , pour la méthode des trapèzes (n = 1).
3
2.6. En déduire les formules d’erreur pour les méthodes composées correspondantes.
Exercice 4 : Ordre et erreur d’une formule d’intégration
Z 1
1
Soit f ∈ C ([0, 1]). On considère la formule de quadrature élémentaire pour approcher f (t)dt :
0
˜ ) = λ0 f (0) + λ1 f (ξ) + λ2 f (0), où ξ ∈]0, 1[ et λ0 , λ1 , λ2 sont des réels.
I(f 0
Z 1
On pose E(f ) = f (t)dt − (λ0 f (0) + λ1 f (ξ) + λ2 f 0 (0)).
0
4.1. Déterminer les paramètres λ0 , λ1 , λ2 et ξ pour que le formule de quadrature soit exacte si f est
un polynôme de degré inférieur ou égal à 3.
4.2. Les paramètres ξ, λ0 , λ1 , λ2 ainsi déterminés, calculer E(x → x4 ) et en déduire l’ordre de la
méthode. Déterminer le noyau de Péano K associé à la méthode.
4.3. Vérifier avec scilab que K est toujours positif. En déduire une expression de l’erreur E(f ) lorsque
f ∈ C 4 ([0, 1]).
4.4. A l’aide d’un changement de variable, construire une méthode de quadrature élémentaire sur
un intervalle [a, b] et donner la valeur de l’erreur pour la formule composée.
1
2
Exercice 5 : Calcul de méthodes de Gauss
5.1. Calculer les méthodes de Gauss associées au poids ω = 1 sur [−1, 1] dans le cas de 1, 2, puis 3
points. √
5.2. Calculer les méthodes de Gauss associées au poids ω = x sur [0, 1] dans le cas de 1 point.
Exercice 6 : Méthodes de Gauss ”fermées”
On reprend les méthodes de Gauss sur [−1, 1] en imposant les points extrêmes de l’intervalle comme
points de la formule de quadrature, à savoir x0 = −1 et xn = 1. Ces méthodes sont appelées méthodes
de Gauss-Lobatto-Legendre.
6.1. La méthode à n + 1 points est d’ordre 2n − 1. Pourquoi ?
6.2. Calculer les méthodes obtenues pour le poids ω = 1 sur [−1, 1] dans le cas de 2 (n = 1), puis 3
(n = 2) et 4 (n = 3) points.
Exercice 7 : Calcul d’intégrale généralisée
Z +∞
Estimer cos t e−t dt en fonction de A ≥ 1. En déduire un programme qui permet de calculer
A
Z +∞
−t
cos t e dt avec une précision ε fixée.
0
On pourra utiliser la méthode de Simpson dont une estimation de l’erreur de la formule composée
Z b
h4
pour calculer f (t)dt est donnée par |ESimpson (f )| ≤ ||f (4) ||∞ (b − a).
a 2880
Exercice 8 : Calcul d’intégrale généralisée
8.1. Montrer que la formule de Gauss à 1 point associée au poids ω(x) = x−β pour 0 < β < 1 sur
Z 1
−β 1 1−β
l’intervalle [0, 1] vaut f (t)t dt ' f .
0 1−β 2−β Z
ε
dt
8.2. En déduire une approximation de l’intégrale √ en utilisant la méthode du point
0 t(1 + t)
milieu, puis la méthode de Gauss calculée ci-dessus. Quelle est la valeur exacte de cette intégrale ?
