Analyse numérique

Chapitre 1 : Zéros d’une fonction non linéaire

Dr. Imed MAHFOUDHI

mahfoudhi.imed5@[Link]

11 septembre 2024
1/27
Imed MAHFOUDHI 2023-2024 ENIM
 PLAN

1 Chapitre I : Zéros d’une fonction non linéaire

Introduction
Méthode de dichotomie (Bissection)
Méthode du point fixe (approximations successives)
Méthode de Newton

2/27
Imed MAHFOUDHI 2023-2024 ENIM
 Chapitre I : Zéros d’une fonction non linéaire

Introduction

L’objet essentiel de ce chapitre est la résolution numérique des

équatiions non linéaires. Étant donné une fonction

f :I⊂R→R (1)
Une valeur α ∈ I, telle que f (α) = 0 est appelée racine ou zéro
de la fonction f .
En générale, il n’existe pas des formules qui permettent le
calcul des valeurs exactes des racines de la fonction f . Pour
cette raison, on s’intéresse à chercher l’approximation des
racines de f à l’aide des méthodes numériques.
Les méthodes employées pour approcher une racine α de f sont
en générale des méthodes itératives : elles consistent à
construire une suite (xn ), telle que

lim xn = α.
n→∞

3/27
Imed MAHFOUDHI 2023-2024 ENIM
 Chapitre I : Zéros d’une fonction non linéaire

Introduction

Definition
Une valeur α est une racine séparée ou localisée de la fonction f
s’il existe un intervalle J ⊂ I, tel que α est l’unique racine f
dans l’intervalle J .

Definition
Soit xn une suite convergente vers la solution α de l’équation
f (x) = 0. On pose en = xn − α. La suite xn est dite convergente
d’ordre p ≥ 1 s’il existe une constante finie C > 0, telle que
|en+1 |
lim = C.
n→∞ |en |p

- La constante C est appelée facteur de convergence ou

constante asymptotique de l’erreur.
- Si p = 1, on dit que la convergence est linéaire.
4/27 - Si p = 2, on dit que la convergence est quadratique.
Imed MAHFOUDHI 2023-2024 ENIM
 Chapitre I : Zéros d’une fonction non linéaire

Introduction

Remarque
1) Soit xn une suite convergente vers α avec un ordre p et une
constante asymptotique d’erreur C et yn une suite
convergente vers α avec un ordre q et une constante
asymptotique D. Alors si (p > q) ou et (p = q et C < D), on
dit que la suite xn converge plus vite que yn .
2) La convergence des méthodes pour la détermination des
racines d’une équation non linéaire dépend en générale de
la donnée initiale x0 . Le plus souvent, on ne sait établir que
des résultats de convergence locale, c’est à dire valable
seulement pour une valeur x0 appartenant à un certain
voisinage de la racine α pour tout choix de x0 ∈ I sont
dites globalement convergentes vers α.

5/27
Imed MAHFOUDHI 2023-2024 ENIM
 Chapitre I : Zéros d’une fonction non linéaire

Méthode de dichotomie (Bissection)

Méthode de dichotomie (Bissection)

6/27
Imed MAHFOUDHI 2023-2024 ENIM
 Chapitre I : Zéros d’une fonction non linéaire

Méthode de dichotomie (Bissection)

La méthode de dichotomie est basée sur le théorème des

valeurs intermédaires.
Théorème des valeurs intermédaires
Soit f une fonction continue sur I = [a, b], telle que f (a)f (b) < 0.
Alors il existe α ∈ [a, b] tel que f (α) = 0.

La condition f (a)f (b) ≤ 0 signifie que f (a) et f (b) sont de signes

opposés (où que l’un des deux est nul). L’hypothèse de
continuité est essentielle !

7/27
Imed MAHFOUDHI 2023-2024 ENIM
 Chapitre I : Zéros d’une fonction non linéaire

Méthode de dichotomie (Bissection)

Ce théorème affirme qu’il existe au moins une solution de

l’équation (f (x) = 0) dans l’intervalle [a, b]. Pour le rendre
effectif est trouver une solution (approchée) de l’équation
(f (x) = 0), il s’agit maintenant de l’appliquer sur un intervalle
suffisamment petit. On va voir que cela permet d’obtenir un
solution de l’équation (f (x) = 0) comme la limite d’une suite.
Étant donné une fonction continue

f : I = [a, b] → R.
Cette méthode consiste à construire des intervalles emboités
In = [an , bn ], tels que In+1 ⊂ In ⊂ In−1 ⊂ · · · ⊂I où la longueur
de l’intervalle In+1 est la moitié de celle de l’intervalle In en
concervant la propriété f (an )f (bn ) ≤ 0.
Voici comment construire une suite d’intervalles emboı̂tés, dont
la longueur tend vers 0, et contenant chacun une solution de
l’équation (f (x) = 0).
8/27 On part d’une fonction f : [a, b] → R continue, avec a < b, et
Imed MAHFOUDHI 2023-2024 ENIM
 Chapitre I : Zéros d’une fonction non linéaire

Méthode de dichotomie (Bissection)

• Si f (a)f ( a+b a+b

2 ) ≤ 0, alors il existe c ∈ [a, 2 ] tel que f (c) = 0.
• Si f (a)f ( a+b a+b
2 ) > 0, cela implique que f ( 2 )f (b) ≤ 0, donc il
existe c ∈ [ a+b2 , b] tel que f (c) = 0.

Nous avons obtenu un intervalle de longueur moitié dans

lequel l’équation (f (x) = 0) admet une solution. On itère alors
le procédé pour diviser de nouveau l’intervalle en deux. Voici
le processus complet :
9/27
Imed MAHFOUDHI 2023-2024 ENIM
 Chapitre I : Zéros d’une fonction non linéaire

Méthode de dichotomie (Bissection)

• Au rang 0 :On pose a0 = a, b0 = b. Il existe une solution x0

de l’équation (f (x) = 0) dans l’intervalle [a0 , b0 ].
• Au rang 1 :
- Si f (a0 )f ( a0 +b
2 ) ≤ 0, alors on pose a1 = a0 et b1 = x1 =
0 a0 +b0
2 ,
a0 +b0
- sinon on pose a1 = x1 = 2 et b1 = b.
- Dans les deux cas, il existe une solution x1 de l’équation
(f (x) = 0) dans l’intervalle [a1 , b1 ].
• ···
• Au rang n : supposons construit un intervalle [an , bn ], de
longueur b−a2n , et contenant une solution xn de l’équation
(f (x) = 0). Alors :
- Si f (an )f ( an +b
2 ) ≤ 0, alors on pose an+1 = an et
n

an +bn
bn+1 = xn = 2 ,
- sinon on pose an+1 = xn = an +b 2
n
et bn+1 = bn .
- Dans les deux cas, il existe une solution xn+1 de l’équation
(f (x) = 0) dans l’intervalle [an+1 , bn+1 ].
10/27
Imed MAHFOUDHI 2023-2024 ENIM
 Chapitre I : Zéros d’une fonction non linéaire

Méthode de dichotomie (Bissection)

E XEMPLE

Exemple
Appliquer la méthode de la dichotomie pour résoudre f (x) = 0
avec f (x) = x2 − 2.
1 Prendre a0 = 1, b0 = 2 et calculer an , bn pour n = 1, 2, 3.
√
2 Trouver le nombre d’itérations n tel que |an − 2| ≤ 10−10 .

11/27
Imed MAHFOUDHI 2023-2024 ENIM
 Chapitre I : Zéros d’une fonction non linéaire

Méthode de dichotomie (Bissection)

Exercice N :1
Montrer que la suite xn est convergente vers α la racine de l’équation
f (x) = 0.

Proposition
Les suites an , bn et xn sont construites de façon que

[an , bn ] ⊂ [an−1 , bn−1 ] ⊂ · · · ⊂ [a0 , b0 ]

bn−1 − an−1
bn − an = tel que f (an )f (bn ) ≤ 0, ∀n ≥ 1
2
et
bn + an
xn = , ∀n ≥ 1
2
Il est clair que an et bn sont respectivement croissante et décroissante
et on a |bn − an | = 21 |bn−1 − an−1 | = 2−n |b − a| → 0, quand n →
+∞.
12/27
Imed MAHFOUDHI 2023-2024 ENIM
 Chapitre I : Zéros d’une fonction non linéaire

Méthode du point fixe (approximations successives)

Méthode du point fixe (approximations

successives)

13/27
Imed MAHFOUDHI 2023-2024 ENIM
 Chapitre I : Zéros d’une fonction non linéaire

Méthode du point fixe (approximations successives)

Definition
(Point fixe) Soit f une fonction de R −→ R. On dit que x0 est un
point fixe de f si et seulement si f (x0 ) = x0 .

Definition
(zéro d’une fonction) Soit f une fonction de R −→ R et x0 est un
zéro de f si et seulement si x0 est point fixe de g(x) = f (x) + x.

14/27
Imed MAHFOUDHI 2023-2024 ENIM
 Chapitre I : Zéros d’une fonction non linéaire

Méthode du point fixe (approximations successives)

Exemple
Trouver les équivalences entre les points fixes et les zéros
des fonctions suivantes :
▶ f1 (x) = x2 − x − 3
▶ f2 (x) = x3 − 3x
9 −x
▶ f3 (x) = 10 e
▶ f4 (x) = sin(x) − x

15/27
Imed MAHFOUDHI 2023-2024 ENIM
 Chapitre I : Zéros d’une fonction non linéaire

Méthode du point fixe (approximations successives)

Proposition
1 Soit f : [a, b] → R une fonction continue telle que
f ([a, b]) ⊂ [a, b]. Alors la fonction f admet au moins un
point fixe dans [a, b]. Soit g(x) = f (x) − x pour tout
x ∈ I = [a, b].
2 La fonction g est continue sur I et on a

g(a) = f (a)−a ≥ 0 et g(b) = f (b)−b ≤ 0 ⇒ g(a)g(b) ≤ 0.

D’après le théorème des valeurs intermédiaires, il

existe α ∈ I, tel que g(α) = 0. Ainsi, on a démontré
qu’il existe α ∈ I tel que f (α) = α.

Principe de la méthode Cette methode consiste à écrire

f (x) = 0 sous la forme g(x) = x. La recherche d’une solution de
l’ équation f (x) = 0 est équivalent à un point fixe de g.
16/27
Imed MAHFOUDHI 2023-2024 ENIM
 Chapitre I : Zéros d’une fonction non linéaire

Méthode du point fixe (approximations successives)

Exemple
Soit f (x) = x3 + x − 1, on peut choisir la fonction g définie
par g(x) = 1 − x3 .

Definition
On dit que la fonction f est contractante (où f est une
contraction) s’il existe k ∈ [0, 1[ tel que

|f (x) − f (y)| ≤ k|x − y|, ∀x, y ∈ [a, b].

Le réel k est appelé rapport de contraction

17/27
Imed MAHFOUDHI 2023-2024 ENIM
 Chapitre I : Zéros d’une fonction non linéaire

Méthode du point fixe (approximations successives)

Fonction k-lipschitzienne
Soit f : [a, b] −→ R est une fonction différentiable sur [a, b] inclu
dans R. Si sup |f ′ (x)| ≤ k pour tout x dans [a, b], alors f est
x∈[a,b]
k-lipschitzienne sur [a, b].

Théorème du point fixe

Soit f : [a, b] −→ R est une fonction contractante telle que

f ([a, b]) ⊂ [a, b]. Alors f admet un unique point fixe α ∈ [a, b]. De
plus, pour tout point initial x0 ∈ [a, b] la suite (xn )n≥0 définie
par xn+1 = f (xn ) converge vers α ∈ [a, b].

18/27
Imed MAHFOUDHI 2023-2024 ENIM
 Chapitre I : Zéros d’une fonction non linéaire

Méthode du point fixe (approximations successives)

C OROLLAIRE

▶ Le résultat du Théorème 18 reste valable si on remplace

l’hypothèse ”f est contractante ” par : f de classe C1 sur
[a, b] et supa≤x≤b |f ′ (x)| < 1.
▶ Soit f une fonction réelle admettant un point fixe α.
Supposons que f est continûment différentialle sur un
voisinage de α. Alors
(i) Si |f ′ (α)| < 1, la méthode du point fixe est localement
convergente, i.e., il existe un voisinage V de α tel que, pour
tout choix de x0 ∈ V, la suite (xn ) définie par xn+1 = f (xn )
converge vers α.
(ii) Si |f ′ (α)| > 1, la méthode du point fixe est divergente pour
tout choix de x0 .

19/27
Imed MAHFOUDHI 2023-2024 ENIM
 Chapitre I : Zéros d’une fonction non linéaire

Méthode du point fixe (approximations successives)

E XERCICE
Exercice N :2
On se propose de résoudre numériquement l’équation

(E1 ) : f (x) = x3 + x − 1

(1) Montrer que (E1 ) admet une solution réelle unique α et

que α ∈]0, 1[.
(2) Montrer que l’équation (E1 ) est équivalente à l’équation

(E2 ) : gi (x) = x

pour i = 1, 2, 3 on a,

1
g1 (x) = x3 + 2x − 1, g2 (x) = 1 − x3 et g3 (x) = .
1 + x2
20/27
Imed MAHFOUDHI 2023-2024 ENIM
 Chapitre I : Zéros d’une fonction non linéaire

Méthode du point fixe (approximations successives)

E XERCICE

(3) Etudier dans chacun des trois cas la convergence de la

suite (xn )n∈N définie par la méthode de point fixe

x0 ∈ I = [0, 1]
xn+1 = g(xn ).

pour la recherche de α.
(4) Dans le cas où il y a une convergence, donner un intervalle
I0 ⊂ I tel que la méthode converge pour tout choix x0 ∈ I0 .

21/27
Imed MAHFOUDHI 2023-2024 ENIM
 Chapitre I : Zéros d’une fonction non linéaire

Méthode du point fixe (approximations successives)

E XERCICE

Exercice N :3
Soit la fonction f définie par f (x) = 2x3 − x − 2. On se propose de
trouver les racines réelles de f par la (Méthode de point fixe).
(1) Montrer que que f posséde une seule racine réelle α ∈ [1, 2]
(2) Etudier la convergence des trois méthodes itératives
suivantes : x0 ∈ [1, 2] donné et
a) xn+1 = 2x3n − 2.
2
b) xn+1 = 2x2n −1
.
xn

13
c) xn+1 = 1 + 2 .

22/27
Imed MAHFOUDHI 2023-2024 ENIM
 Chapitre I : Zéros d’une fonction non linéaire

Méthode de Newton

Méthode de Newton

23/27
Imed MAHFOUDHI 2023-2024 ENIM
 Chapitre I : Zéros d’une fonction non linéaire

Méthode de Newton

Ètant donné une fonction f de classe C2 , pour déterminer la

racine de l’équation f (x) = 0, la méthode de Newton consiste à
construire à aprtir de x0 une suite (xn ) définie par
f (xn )
xn+1 = xn − .
f ′ (xn )
Interpretation géometrique : xn+1 est l’intersection de la
tangente à la courbe Cf de la fonction f au point (xn , f (xn )) avec
l’axe des abscisses. En fait, l’équation de la tangente T à Cf en
(xn , f (xn )) est donnée par T : y = f ′ (xn )(x − xn ) + f (xn ).
Remarque : Il faut que f ′ (xn ) ̸= 0.
La méthode de Newton est une méthode de point fixe pour la
f (x)
fonction ψ(x) = x − f ′ (x) .

f (x)
ψ(x) = x − f ′ (x) si f ′ (x) ̸= 0
f (x ) (2)
xn+1 = xn − f ′ (xnn )
24/27
Alors on a : f (α) = 0 si et seulement
Imed MAHFOUDHI 2023-2024 si α est un point
ENIM fixe de ψ.
 Chapitre I : Zéros d’une fonction non linéaire

Méthode de Newton

Théorème de Newton
Si f ∈ C2 , f (α) = 0 et f ′ (α) ̸= 0 alors il existe δ telle que, si
|x0 − α| ≤ δ la suite définie dans (2) converge vers α. De plus, la
xn+1 − α f ′′ (α)
convergence est quadratique : lim =
n→∞ (xn − α)2 2f ′ (α)

Remarque :
▶ Si f ′ (α) ̸= 0 alors ψ ′ (α) = 0,la convergence est quadratique.
▶ Si f ′ (α) = 0, la convergence est linéaire. On considère alors
la méthode de Newton modifiée :
f (x )
xn+1 = xn − p f ′ (xnn ) n ≥ 0 avec p la multiplicité de α.

25/27
Imed MAHFOUDHI 2023-2024 ENIM
 Chapitre I : Zéros d’une fonction non linéaire

Méthode de Newton

E XERCICE
Exercice N :4
Soit α une racine double de la fonction f tel que

f (α) = f ′ (α) = 0 et f (2) (α) ̸= 0.

En tenant compte du fait qu’on peut écrire la fonction f comme

f (x) = (x − α)2 h(x) où h(α) ̸= 0.

1 Vérifier que la méthode de Newton pour l’approximation

de la racine α est seulement d’ordre 1.
On considère la méthode de Newton modifiée suivante :
f (xk )
xk+1 = xk − 2 ′ .
f (xk )

2 Vérifier que cette méthode est d’ordre deux si l’on veut

26/27 approcher α.
Imed MAHFOUDHI 2023-2024 ENIM
 Chapitre I : Zéros d’une fonction non linéaire

Méthode de Newton

Merci pour votre attention

27/27
Imed MAHFOUDHI 2023-2024 ENIM