Chapitre 2

CHAPITRE 2
RESOLUTION DES EQUATIONS NON LINEAIRES
1. Introduction
La résolution des problèmes de l’ingénieur débouche souvent sur deux types de recherche
de solutions d’une équation à une variable :
a.) recherche d’une racine de l’équation f ( x )  0 où est une fonction transcendante
ou numérique de x.
Les méthodes numériques que nous présenterons dans ce chapitre conduiront (sous
certaines conditions) à l’approximation d’une racine de l’équation f( x)0.
Pour l’ingénieur, la recherche des racines complexes de f ( x )  0 est relativement peu
courante sauf en commande des processus. C’est pourquoi la méthode dite de Newton est
généralement présentée pour rechercher les racines complexes de :
f x  0 (1)
b.) Recherche de plusieurs ou toutes les racines de Pn ( x )  0 où Pn ( x ) est un
polynôme de degré n en x . Ce problème est un cas particulier du premier ; il pourrait

donc être résolu par la méthode précédente.
2. Méthode de la bissection (dichotomie)
La méthode de la bissection repose sur l’idée toute simple qu’en général de part et d’autre
d’une solution de l’équation (1) f  x   0 une fonction continue f ( x ) change de signe
et passe du positif au négatif ou vice versa.
Supposons qu’il y ait effectivement un changement de signe autour d’une racine r de
f x.
f x
r
x1 xm x2
Analyse numérique pour Ingénieur : Par EGBOHO M. Franck/EPAC-UAC Page 20

Soit  x1 , x2  un intervalle ayant un changement de signe de f  x  c à d :
f  x1   f  x 2   0
(2)
On pose :
x  x2
xm  1
2
le point milieu de l’intervalle  x1 , x2  . Il s’agit alors de déterminer entre les intervalles
 x1 , xm  et  xm , x2  celui qui possède encore un changement de signe et la racine se

trouvera forcement dans cet intervalle d’où l’algorithme suivant :
Algorithme (algorithme de la bissection)
1. Etant donné un intervalle  x1 , x2  pour lequel f ( x ) possède un changement de

signe ;
2. Etant donné  le critère d’arrêt et N le nombre maximal d’itérations ;
x1  x 2
3. Poser : xm 
2
x 2  x1
4. Si  :
2 xm
 Convergence atteinte ;
 Ecrire la racine xm ;
 Ecrire f xm ;
 Arrêt.
     
5. Ecrire x1 , x2 , xm , f x1 , f x2 , f xm ;
   
6. Si f x1  f xm  0 alors x2  xm ;
7. Si f  xm   f  x2   0 alors x1  xm ;
8. Si le nombre maxi. d’itérations N est atteint :
 convergence non atteinte en N itérations ;
 Arrêt.
9. Retour à l’étape 3.

Remarque
1. L’expression :
x 2  x1

2 xm
est une approximation de l’erreur relative. En effet, à l’étape 3 de l’algorithme, la racine
recherchée est soit dans l’intervalle  x1 , xm  ou  xm , x2  qui sont toutes deux de
longueur :
x2  x1
2
ce qui constitue une borne supérieure de l’erreur absolue. Ainsi, en divisant par xm ,
on obtient une approximation assez fiable de l’erreur relative.

2. Prendre garde au cas où la racine recherchée est 0 car il y a risque de division par 0
au cours de l’évaluation de l’erreur relative. Ce cas est toutefois rare en pratique.
3. Il est parfois utile d’introduire un test d’arrêt sur la valeur de f  x  qui doit tendre
également vers 0.
Algorithme pratique de la méthode de bissection

Données : x1 , x2 {bornes de l’intervalle}
E {précision souhaitée}
f(x) {fonction}
Calculs :
 x2  x1 
 ln 
N  partie entière  E  1
 ln 2 
 
y1  f  x1 
pour i  1 à N

  x1  x2  ; y  f x
 xm 
2 m  m

 si y1 ym 0 alors x2  xm
 sinon x1  xm ; y1  ym


écrire xm
Exercice :
En utilisant la méthode de la bissection, chercher la racine carrée de 2 en prenant
l’intervalle 1 , 2 pour une précision de E = 10-3.
Solution
Calculons d’abord le nombre d’itérations nécessaires pour obtenir la précision de
E  10 3 .
Nous avons :
 x2  x1   21 
ln  ln 3 
 
N  INT  E  1  INT  10  1
 ln 2   ln 2 
   
 INT ( 10 ,966 )  10
x  2  x 2  2  0; f ( x )  x 2  2  0
y1  f  x1   f  1   1  2  1
y2  f  x 2   f  2   4  2  2
On commence les calculs itératifs
1ère itération
x  x2 1  2
xm  1   1,5
2 2
ym  f  xm   1,5 2  2  0 ,25
y1  ym 0 donc on sélectionne l'intervalle  x1 , xm  (en prenant
x2  xm ) et ainsi de suite.
Les résultats obtenus en fonction du nombre d’itérations sont :

x1 x2 xm ym
N
1 1 2 1,5 0,25
2 1 1,5 1,25 - 0,4375
3 1,25 1,5 1,375 - 10,94
4 1,375 1,5 1,4375 0,06641
5 1,375 1,4375 1,4062 - 0,02246
6 1,4062 1,4375 1,4219 0,02173
7 1,4141 1,4219 1,4141 - 0,00043
8 1,4141 1,4219 1,4180 0,01064
9 1,4141 1,4180 1,4160 0,00510
10 1,4141 1,4160 1,4150 0,00233
Exercice :
Montrer que la fonction f ( x )  x 3  x 2  3x  3 possède un zéro dans l’intervalle

1, 2 . Trouver sa valeur approchée en 10 itérations.
Solution
f ( 1 )  4,0; f ( 2 )  3,0
f ( 1 )  f ( 2 )  4  3  12 0
Alors, la fonction f ( x ) donnée possède un zéro dans  1, 2  .

Recherchons cette solution par la méthode de la bissection :
1 2
Le point milieu de  1, 2  est : xm   1,5
2
f  xm   f ( 1,5 )  1,875
f ( 1,5 )  f ( 2 ) 0  l’intervalle  1,5 , 2  possède encore un changement de signe,

ce qui n’est pas le cas pour l’intervalle  1 , 1,5  puisque f(1)  f(1,5) 0 . Ensuite on
calcule le point milieu xm de  1,5 , 2  et ainsi de suite :

Erreur absolue
n x1 x2 xm f ( x1 ) f ( x2 ) f ( xm ) liée à xm
1
x  x1
2 2
1 1,0 2,0 1,5 -4,0 3,0 -1,875 0,5
2 1,5 2,0 1,75 -1,875 3,0 0,17187 0,25
3 1,5 1,75 1,625 -1,875 0,17187 -0,94335 0,125
4 1,625 1,75 1,6875 -0,94335 0,17187 -0,40942 0,0625
5 1,6875 1,75 1,71875 -0,40942 0,17187 -0,12478 0,03125
6 1,71875 1,75 1,734375 -0,12478 0,17187 0,022029 0,015625
7 1,71875 1,734375 1,72656 -0,12478 0,022029 -0,05175 0,007812
8 1,72656 1,734375 1,73046 -0,05175 0,022029 -0,014972 0,0039075
9 1,73046 1,734375 1,7324175 -0,014972 0,022029 0,1994533 0,0019525
10 1,73046 1,7324175 1,731438 -0,014972 0,1994533 0,190645 0,0019575
On remarque aisément que la longueur de l’intervalle entourant la racine est divisée par 2
à chaque itération. Cette constatation permet de déterminer à l’avance le nombre d’itérations
nécessaires pour obtenir une certaine erreur absolue ∆r sur la racine r.
Soit par exemple L  x2  x1 la longueur de l’intervalle de départ, après une itération le
L
nouvel intervalle est de longueur et après N itérations la longueur de l’intervalle est :
2
L
2N
Pour connaître la valeur de N nécessaires pour avoir :
L
r
2N
il suffit de résoudre cette inéquation en fonction de N pour trouver la condition
suivante :
 L
ln  
N  r 
ln 2

Chapitre 2

Transféré par

Droits d'auteur :

Formats disponibles

Chapitre 2

Transféré par

Informations du document

Titre original

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

Chapitre 2

Transféré par

Droits d'auteur :

Formats disponibles

CHAPITRE 2

RESOLUTION DES EQUATIONS NON LINEAIRES

b.) Recherche de plusieurs ou toutes les racines de Pn ( x )  0 où Pn ( x ) est un

polynôme de degré n en x . Ce problème est un cas particulier du premier ; il pourrait

Analyse numérique pour Ingénieur : Par EGBOHO M. Franck/EPAC-UAC Page 20

 x1 , xm  et  xm , x2  celui qui possède encore un changement de signe et la racine se

1. Etant donné un intervalle  x1 , x2  pour lequel f ( x ) possède un changement de

Analyse numérique pour Ingénieur : Par EGBOHO M. Franck/EPAC-UAC Page 21

recherchée est soit dans l’intervalle  x1 , xm  ou  xm , x2  qui sont toutes deux de

on obtient une approximation assez fiable de l’erreur relative.

Algorithme pratique de la méthode de bissection

Analyse numérique pour Ingénieur : Par EGBOHO M. Franck/EPAC-UAC Page 22

Analyse numérique pour Ingénieur : Par EGBOHO M. Franck/EPAC-UAC Page 23

Montrer que la fonction f ( x )  x 3  x 2  3x  3 possède un zéro dans l’intervalle

Alors, la fonction f ( x ) donnée possède un zéro dans  1, 2  .

f ( 1,5 )  f ( 2 ) 0  l’intervalle  1,5 , 2  possède encore un changement de signe,

Analyse numérique pour Ingénieur : Par EGBOHO M. Franck/EPAC-UAC Page 24

Soit par exemple L  x2  x1 la longueur de l’intervalle de départ, après une itération le

Analyse numérique pour Ingénieur : Par EGBOHO M. Franck/EPAC-UAC Page 25

Vous aimerez peut-être aussi