Tronc Commun PC Université Abdelmalek Essaâdi
Module : Analyse Numérique Faculté des Sciences de Tétouan
Pr. Y. El Yazidi 2024-2025
TD3 Interpolation, intégration et dérivation numérique
Exercice 1. On veut interpoler la fonction f (x) = ln(x) par un polynôme aux points x0 = 1, x1 = 2, x2 = 3,
x3 = 4 et x4 = 5.
1. Trouver une expression algébrique de ce polynôme en résolvant un système linéaire.
2. Estimer la valeur de f (3.2) avec le polynôme trouvé puis calculer l’erreur absolue.
Exercice 2. Dans une expérience on a mesuré toutes les 10 secondes la vitesse (en m/s) d’écoulement de l’eau
dans une conduite cylindrique. On a calculé à l’aide de ces données la table de différences divisées suivante :
ti V (ti ) δ[ti−1 , ti ] δ[ti−2 , ti−1 , ti ] δ[ti−3 , ti−2 , ti−1 , ti ]
0 2, 00
10 1, 89
20 1, 72
30 1, 44
1. Compléter la table.
2. Trouver l’approximation de la vitesse V (en m/s) à t̄ = 15s avec un polynôme de Newton de degré 2.
3. Retrouver une autre approximation en utilisant un polynôme de Newton de degré 3.
4. Estimer l’erreur commise sur la vitesse calculée en (2) sachant que Π4 (t) ≈ Π4 (t3 ).
5. En utilisant la formule des différences finies centrées trouver une approximation de l’accélération d’écoulement
de l’eau γ à t = 20s.
∫ 1
e−x dx avec n = 2. Déduire une
2
Exercice 3. 1. Par la méthode des Trapèze calculé une approximation de
0
estimation d’erreur d’intégration numérique.
∫ 1
2. En utilisant la formule de quadrature de Simpson calculer une approximation de sin3 (x)dx avec n = 4.
0
Déduire une estimation d’erreur d’intégration numérique.
Exercice 4. L’objectif de cet exercice est la construction d’une méthode d’intégration numérique. Pour cela on
procède d’abord par une interpolation. Soit f une fonction définie sur R.
1. Chercher le polynôme d’interpolation P de f aux points − 21 , 0 et 12 .
∫ 1
2. Intégrer P sur l’intervalle [−1, 1], puis trouver une formule de quadrature pour l’intégral f (x)dx.
−1
∫ b
3. Déduire une formule de quadrature pour l’intégral f (x)dx en utilisant le changement de variable suivant
a
b+a b−a
t= 2 + 2 x.
∫ 1
e−x ,
2
4. Application, supposons que f (x) = trouver une approximation de f (x)dx.
0
5. Comparer le résultat trouvé avec la question 1 de l’exercice 1.
Exercice 5. Soit l’approximation
∫ x0 +h ( ( ))
h 2h
f (x) dx ≈ f (x0 ) + 3f x0 +
x0 4 3
( )
1. Obtenir un développement de Taylor de f x0 + 2h3 jusqu’à l’ordre 4 et donner une nouvelle expression du
terme de droite.
2. Obtenir un développement de Taylor à l’ordre 4 du terme de gauche.
3. Soustraire les expressions obtenues en a) et en b) pour obtenir le premier terme de l’erreur. En déduire
l’ordre de la méthode proposée.
4. Quel est le degré de précision de cette méthode.
Exercice 6. On considère le θ-schéma
( ) ( )
′ f (x + h) − f (x) f (x) − f (x − h)
f (x) ≈ (1 − θ) +θ = Appθ (h).
h h
Montrer que les deux premiers termes de l’erreur associée au θ-schéma (Appθ (h)) sont donnés par :
(2θ − 1) ′′ h2
hf (x) − f (3) (x),
2 6
et en déduire l’ordre de précision du θ-schéma en fonction de θ.
Exercice 7. 1. Résoudre par la méthode d’Euler sur R l’équation différentielle suivante :
y ′ + y = 2 sin x avec y(0) = 0
2. Résoudre numériquement par la méthode de Runge-Kutta d’ordre 2 sur R l’équation différentielle suivante :
y ′ − y = (x + 1)ex avec y(0) = 1
3. Utiliser la méthode d’Euler et celles de Runge-Kutta d’ordre 4 pour résoudre numériquement sur R l’équation
suivante :
y ′ + y = x − ex + cos x avec y(0) = 0
Exercice 8. Soit le problème de Cauchy suivant :
{
y ′ (t) = y(t) − 2t
y(t)
y(0) = 1
1. Calculer la solution approximative de cette équation en t = 0.2 à l’aide de méthode de Runge-Kutta d’ordre
2 avec un pas h = 0.2.
√
2. Sachant que la solution exacte est y(t) = 2t + 1, comparer le résultat obtenu avec la solution exacte.
Exercice 9. Soit le problème de Cauchy suivant :
{
y ′ (t) = e−t − 2y(t), t ∈ [0, 1]
y(0) = 1
1. Montrer que la fonction f définie par f (t, y) = e−t − 2y est Lipschitzienne par rapport à y uniformément
par rapport à t, et donner la constante de Lipschitz correspondante.
2. Montrer que ce problème admet une solution unique.
3. Donner la solution exacte de ce problème ainsi que y(0.2).
4. Appliquer la méthode d’Euler à ce problème, et donner l’approximation de y(0.35) obtenue à l’aide d’un pas
de discrétisation h = 0.1.
5. Comparer le résultat obtenu avec la solution exacte.
