TD2-AN-2018

Université Claude Bernard, Lyon 1 Licence Sciences & Technologies
43, boulevard du 11 novembre 1918 Spécialité : Mathématiques

69622 Villeurbanne cedex, France Intro. analyse numérique L2- Printemps 2018
Série d’exercices no 2/5

Calculs approchés d’intégrales
Exercice 1. Partiel 2017- Exercice - Parties 2 et 3
1. Partie 2. (6 points) (30 min)

On pose f (x) = 8x4 − 8x2 + 1.
Z 1
(a) (1 point) Justifier que l’intégrale I = f (x)dx est bien définie et la calculer.
−1
(b) (1 point) Rappeler la formule de quadrature par la méthode de Simpson pour une fonc-
tion f sur un intervalle [a, b] où a, b ∈ R et a < b.
(c) (1 point) Donner la valeur de l’approximation de cette intégrale obtenue par la méthode
de Simpson.
(d) (1 point) Calculer l’erreur exacte en se servant du calcul explicite de I.
(e) (1 point) Rappeler la formule de l’erreur de quadrature pour cette méthode. Donnez une
estimation de l’erreur pour la fonction f de cette partie, et la comparer avec l’erreur
calculée dans la question précédente.
(f) (1 point) Est-ce que cela confirme ce que l’on savait sur l’ordre de cette méthode ?
2. Partie 3. (4 points) (30 min)
(a) (2 points) On note z1 et z2 deux réels non nuls de ] − 1, 1[, tels que z1 < z2 et α, β, γ
trois coefficients réels.
Déterminer z1 , z2 , α, β et γ tels que l’expression
G(f ) = αf (z1 ) + βf (0) + γf (z2 ),
soit une formule d’intégration numérique exacte pour les polynômes de degré inférieur
ou égal à 4 sur [−1, 1].
(b) (2 points) Cette méthode est-elle d’ordre 5 ? D’ordre 6 ?
1
Exercice 2. Simpson
Soit p la fonction polynôme de la variable réelle x définie par
35 4 15 2
p(x) = x − x + 3.
16 2
1. Calculer l’intégrale
Z 2
p(x)dx.
−2
2. Donner la valeur de l’approximation de cette intégrale obtenue par la méthode de Simpson

appliquée aux intervalles successifs [−2, 0] et [0, 2]. Donner une estimation de l’erreur.
Z 2
3. On note M (ϕ) l’approximation d’une intégrale de la forme ϕ(x)dx :
−2
8 4 8
M (ϕ) = ϕ(−1) − ϕ(0) + ϕ(1).
3 3 3
Z 2
(a) Montrer que M (ϕ) = ϕ(x)dx si ϕ est un polynôme de degré inférieur ou égal à 3.
−2
Z 2
(b) Comparer M (p) et p(x)dx.
−2
4. On note z1 < z2 deux nombres réels distincts de [−2, 2] et α, β et γ trois nombres réels.
(a) Déterminer z1 , z2 , α, β et γ de telle sorte que si l’on pose
N (ϕ) = αϕ(z1 ) + βϕ(0) + γϕ(z2 ),
nous ayons
Z 2
N (ϕ) = ϕ(x)dx si ϕ est un polynôme de degré inférieur ou égal à 4.
−2
(b) Vérifier les résultats trouvés à l’aide du polynôme p.
Exercice 3. Newton-Cotes
Soit f : [−1, 1] → R une fonction continue.
Calculer les coefficients des formules suivantes
Z 1
1 1
1. f (x)dx ' λ0 f (− ) + λ1 f ( ), exacte pour f ∈ P1 .
−1 3 3
Z 1
3 1 1 3
2. f (x)dx ' λ0 f (− ) + λ1 f (− ) + λ2 f ( ) + λ3 f ( ), exacte pour f ∈ P3 .
−1 5 5 5 5
3. Construire la formule composée correspondante à la première question avec des nœuds
équidistants et trouver l’expression de l’erreur commise.
2
Exercice 4. Formules de Gauss-Lobatto
Dans ce type de formule de quadrature, les nœuds extrêmes sont imposés et égaux aux extrémités
de l’intervalle d’intégration.
Pour [a, b] = [−1, 1], la formule de Gauss-Lobatto à deux nœuds intérieurs approchant l’intégrale
Z 1
I(f ) = f (x)dx,
−1
s’écrit en tenant compte de la symétrie :
J(f ) = λ (f (−1) + f (1)) + µ (f (−α) + f (α)).
1. Calculer les valeurs de λ et µ pour avoir la formule d’ordre la plus élevée possible.
2. Quel est cet ordre ?
Exercice 5. Méthode de Simpson

Soient ϕ une fonction continue sur R et α ∈ R∗+ un réel fixé. On considère la fonction ψ : [0, α] →
R définie pour tout t ∈ [0, α] par
Z t
t
ψ(t) = ϕ(s)ds − (ϕ(−t) + 4ϕ(0) + ϕ(t)) − kt2 ,
−t 3
où k est un réel défini par la condition ψ(α) = 0. On suppose que ϕ ∈ C 4 ([−α, α]).
1. Calculer pour tout t ∈ [0, α], ψ 0 (t), ψ 00 (t) et ψ (3) (t).
2. En utilisant le théorème de Rolle, successivement pour ψ, ψ 0 et ψ 00 montrer qu’il existe
θ ∈]0, α[ tel que ψ (3) (θ) = 0.
3. En déduire qu’il existe η ∈] − α, α[ tel que k s’exprime en fonction de ϕ(4) (η).
4. Retrouver ainsi la formule d’intégration de Simpson pour ϕ sur un intervalle [−α, α] :
Z α
α α5
ϕ(s)ds = (ϕ(−α) + 4ϕ(0) + ϕ(α)) − ϕ(4) (η), η ∈] − α, α[.
−α 3 90
5. On considère une fonction f ∈ C (4) ([a, b]) (a, b ∈ R, a < b). Soient xi = a+ih, i = 0, ..., n
(b − a)
où h = . En utilisant un changement de variable et la question précédente, retrouver
n
la méthode d’intégration de Simpson pour f sur [a, b], à savoir
Z b n−1
hX xi + xi+1
f (x)dx = f (xi ) + 4f ( ) + f (xi+1 ) + e(f ),
a 6 i=0 2
b − a 4 (4)
où e(f ) = − h f (c) avec c ∈ [a, b].
2880
3
Exercice 6. Méthode des trapèzes-Hermite
1. Soient ϕ ∈ C 4 ([−1, 1]) et p3 le polynôme d’Hermite de de ϕ de degré ≤ 3 vérifiant

p3 (−1) = ϕ(−1), p03 (−1) = ϕ0 (−1),
p3 (1) = ϕ(1), p03 (1) = ϕ0 (1).
(a) Déterminer l’expression du polynôme p3
(b) En déduire la formulation d’intégration numérique élémentaire suivante
Z 1
1 2
ϕ(s)ds = ϕ(−1) + ϕ(1) + (ϕ0 (−1) − ϕ0 (1)) + ϕ(4) (η), avec η ∈ [−1, 1].
−1 3 45
2. En déduire la formule d’intégration numérique composée
Z b
h2 0 0 h4
f (x)dx = Tn (f ) + ((f (a) − f (b)) + (b − a)f (4) (c),
a 12 720
b−a
où f ∈ C 4 ([a, b]), h = , n ∈ N∗ , c ∈ [a, b] et Tn (f ) est la valeur obtenue par la
n
méthode des trapèzes composés avec le pas h.

TD2-AN-2018

Transféré par

Droits d'auteur :

Formats disponibles

TD2-AN-2018

Transféré par

Informations du document

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

TD2-AN-2018

Transféré par

Droits d'auteur :

Formats disponibles

Université Claude Bernard, Lyon 1 Licence Sciences & Technologies

43, boulevard du 11 novembre 1918 Spécialité : Mathématiques

Série d’exercices no 2/5

Exercice 1. Partiel 2017- Exercice - Parties 2 et 3

1. Partie 2. (6 points) (30 min)

2. Donner la valeur de l’approximation de cette intégrale obtenue par la méthode de Simpson

(b) Vérifier les résultats trouvés à l’aide du polynôme p.

s’écrit en tenant compte de la symétrie :

J(f ) = λ (f (−1) + f (1)) + µ (f (−α) + f (α)).

Exercice 5. Méthode de Simpson

1. Soient ϕ ∈ C 4 ([−1, 1]) et p3 le polynôme d’Hermite de de ϕ de degré ≤ 3 vérifiant

Vous aimerez peut-être aussi