Fiche de TD N°2

November 1, 2024

1 Exercice 1
On considère la fonction suivante :

f (x) = x3 − 4x2 + 6x − 2

Nous souhaitons approximer sa dérivée première f ′ (x) en x = 2 en utilisant les

différences finies avec un pas h = 0.1.

2 Solution:
1. Méthodes des différences finies
Pour approximer f ′ (x) en utilisant les différences finies, nous pouvons utiliser
trois méthodes :
1. **Méthode de la différence finie avant** :
f (x + h) − f (x)
f ′ (x) ≈
h
2. **Méthode de la différence finie arrière** :
f (x) − f (x − h)
f ′ (x) ≈
h
3. **Méthode de la différence finie centrée** :
f (x + h) − f (x − h)
f ′ (x) ≈
2h
2. Calculs des valeurs de la fonction
Pour calculer les valeurs de f (x), f (x + h) et f (x − h), nous avons :

f (2) = 23 − 4 · 22 + 6 · 2 − 2 = 8 − 16 + 12 − 2 = 2

f (2.1) = (2.1)3 − 4 · (2.1)2 + 6 · 2.1 − 2 = 9.261 − 17.64 + 12.6 − 2 = 2.221

f (1.9) = (1.9)3 − 4 · (1.9)2 + 6 · 1.9 − 2 = 6.859 − 14.44 + 11.4 − 2 = 1.819
3. Approximation de la dérivée en x = 2

1
 Nous appliquons les méthodes de différences finies aux valeurs calculées :
- **Différence finie avant :**
f (2.1) − f (2) 2.221 − 2
f ′ (2) ≈ = = 2.21
h 0.1
- **Différence finie arrière :**
f (2) − f (1.9) 2 − 1.819
f ′ (2) ≈ = = 1.81
h 0.1
- **Différence finie centrée :**
f (2.1) − f (1.9) 2.221 − 1.819
f ′ (2) ≈ = = 2.01
2h 0.2
Résultats
Les approximations obtenues pour f ′ (2) sont :
- **Différence finie avant** : f ′ (2) ≈ 2.21 - **Différence finie arrière** :
f (2) ≈ 1.81 - **Différence finie centrée** : f ′ (2) ≈ 2.01
′

La méthode des différences finies centrée offre souvent une meilleure précision
pour l’approximation de la dérivée.

3 Exercice 2 de différentiation numérique par la

méthode des différences finies
On considère la fonction suivante :

g(x) = sin(x) + x2

Nous souhaitons approximer sa dérivée première g ′ (x) en x = π

4 en utilisant les
différences finies avec un pas h = 0.01.

3.1 solution
1. Méthodes des différences finies
Nous utiliserons les trois méthodes de différences finies pour approximer
g ′ (x) :
1. **Méthode de la différence finie avant** :
g(x + h) − g(x)
g ′ (x) ≈
h
2. **Méthode de la différence finie arrière** :
g(x) − g(x − h)
g ′ (x) ≈
h
3. **Méthode de la différence finie centrée** :
g(x + h) − g(x − h)
g ′ (x) ≈
2h

2
 2. Calculs des valeurs de la fonction
π
Pour calculer les valeurs de g(x), g(x + h) et g(x − h) avec x = 4 et h = 0.01,
nous obtenons :

 π 2 √
π π 2  π 2
g = sin + + = ≈ 0.7071 + 0.6169 = 1.3240
4 4 4 2 4
π  π  π 2
g + 0.01 = sin + 0.01 + + 0.01 ≈ 0.7141 + 0.6365 = 1.3506
4 4 4
π  π  π 2
g − 0.01 = sin − 0.01 + − 0.01 ≈ 0.7000 + 0.5974 = 1.2974
4 4 4
3. Approximation de la dérivée en x = π4
Appliquons maintenant les méthodes de différences finies aux valeurs cal-
culées :
- **Différence finie avant :**
  g π + 0.01
− g π
1.3506 − 1.3240
′ π 4 4
g ≈ = = 2.66
4 h 0.01
- **Différence finie arrière :**
  g π
− g π − 0.01
1.3240 − 1.2974
′ π 4 4
g ≈ = = 2.66
4 h 0.01
- **Différence finie centrée :**
 π  g π + 0.01
− g π
− 0.01


1.3506 − 1.2974
g′ ≈ 4 4
= = 2.66
4 2h 0.02
Résultats
Les approximations obtenues pour g
′ π4 sont les mêmes dans ce cas :

-
**Différence finie avant** : g ′ π4 ≈ 2.66 -
**Différence finie arrière** :

g 4 ≈ 2.66 - **Différence finie centrée** : g ′ π4 ≈ 2.66
′ π

La méthode des différences finies centrée est généralement plus précise, mais
ici toutes les méthodes donnent des résultats identiques.

4 Exercice : Equations différentielles ordinaires

On considère le problème suivant :
dy
= x + y, y(0) = 1
dx
Nous souhaitons approximer la solution y(x) aux points x = 0.1 et x = 0.2 en
utilisant les méthodes suivantes avec un pas h = 0.1 :
1. La méthode d’Euler. 2. La méthode de Runge-Kutta d’ordre 4.

3
4.1 Solution
1. Méthode d’Euler
La méthode d’Euler est donnée par la formule suivante :

yn+1 = yn + h · f (xn , yn )

où f (x, y) = x + y, x0 = 0, et y0 = 1.
Itérations
1. **Première itération** (x0 = 0, y0 = 1) :

y1 = y0 + h · f (x0 , y0 ) = 1 + 0.1 · (0 + 1) = 1 + 0.1 = 1.1

2. **Deuxième itération** (x1 = 0.1, y1 = 1.1) :

y2 = y1 + h · f (x1 , y1 ) = 1.1 + 0.1 · (0.1 + 1.1) = 1.1 + 0.1 · 1.2 = 1.1 + 0.12 = 1.22

Donc, l’approximation de y(0.2) avec la méthode d’Euler est y ≈ 1.22.

2. Méthode de Runge-Kutta d’ordre 4
La méthode de Runge-Kutta d’ordre 4 est définie par les étapes suivantes
:
1. Calculer k1 = h · f (xn , yn ). 2. Calculer k2 = h · f xn + h2 , yn + k21 . 3.
Calculer k3 = h · f xn + h2 , yn + k22 . 4. Calculer k4 = h · f (xn + h, yn + k3 ). 5.

Mettre à jour yn+1 = yn + k1 +2k2 +2k

6
3 +k4
.
Itérations
1. **Première itération** (x0 = 0, y0 = 1) :

k1 = h · f (x0 , y0 ) = 0.1 · (0 + 1) = 0.1

 
h k1
k2 = h · f x0 + , y0 + = 0.1 · (0.05 + 1.05) = 0.1 · 1.1 = 0.11
2 2
 
h k2
k3 = h · f x0 + , y0 + = 0.1 · (0.05 + 1.055) = 0.1 · 1.105 = 0.1105
2 2
k4 = h · f (x0 + h, y0 + k3 ) = 0.1 · (0.1 + 1.1105) = 0.1 · 1.2105 = 0.12105
k1 + 2k2 + 2k3 + k4 0.1 + 2 · 0.11 + 2 · 0.1105 + 0.12105
y1 = y0 + = 1+ ≈ 1.1103
6 6
2. **Deuxième itération** (x1 = 0.1, y1 ≈ 1.1103) :

k1 = h · f (x1 , y1 ) = 0.1 · (0.1 + 1.1103) = 0.1 · 1.2103 = 0.12103

 
h k1
k2 = h·f x1 + , y1 + = 0.1·(0.15 + 1.170815) = 0.1·1.320815 = 0.1320815
2 2
 
h k2
k3 = h·f x1 + , y1 + = 0.1·(0.15 + 1.17634075) = 0.1·1.32634075 = 0.132634075
2 2
k4 = h·f (x1 +h, y1 +k3 ) = 0.1·(0.2+1.242934075) = 0.1·1.442934075 = 0.1442934075

4
 k1 + 2k2 + 2k3 + k4 0.12103 + 2 · 0.1320815 + 2 · 0.132634075 + 0.1442934075
y2 = y1 + = 1.1103+ ≈ 1.2428
6 6
Ainsi, l’approximation de y(0.2) avec la méthode de Runge-Kutta d’ordre 4
est y ≈ 1.2428.
Conclusion
Les approximations de y(0.2) sont : - **Méthode d’Euler** : y ≈ 1.22 -
**Méthode de Runge-Kutta d’ordre 4** : y ≈ 1.2428
La méthode de Runge-Kutta d’ordre 4 fournit une approximation plus précise
par rapport à la méthode d’Euler.