Université M’hammed Bougarra Boumerdes :
Faculté des hydrocarbures et de la chimie :
Compte rendu de TP
méthodes numérique
Fait par :
SALIM Haithem Ammar
BELAMRI Abderraouf
MOSETEFA Saadeddine
LAMU 19
But du TP :
Implémentation d’un code sous Matlab calculant la racine d’une équation non linéaire par la
méthode de la Dichotomie.
[Link]éambule :
Soit la fonction f (x) = e−√x − 2 + x2 définie sur R+.
a. Ecrire une fonction Matlab qui reçoit un argument x et renvoie la valeur de f (x) :
b. Tracer f (x) avec fplot en appelant la fonction précédente dans un autre script :
c. Localiser la racine de l’équation f (x) = 0 dans un intervalle utilisable par la Dichotomie :
[1 , 2]
d. Trouver cette racine avec la fonction fzero :
2. Algorithme de la Dichotomie :
Supposant que [a, b] est l’intervalle contenant la racines de l’équation f (x) = 0 et satisfaisant la
condition d’application de la Dichotomie f (a). f (b) < 0. On recherche la racine r avec une
précision e (au sens étroit) suivant l’algorithme suivant :
lire a, b, e
r
a+b
2
tant que
b−a
2
> e faire
si f (a). f (r) < 0 alors
br
sinon
a r
fin si
r
a+b
2
fin tant que
afficher r
a. Ecrire le code Matlab correspondant à l’algorithme précédent :
f=@(x) exp(-sqrt(x))-2+x^2
a=1; %début de l'intervalle
b=2; %fin de l'intervalle
k=0; %début d'initiale de k
while k<17 %17 itération
xm=(a+b)/2;
if
f(xm)==0
disp(xm)
elseif f(a)*f(xm)<0
b=xm;
else
a=xm;
end
k=k+1;%agmenter k avec +1 chaque itération
end
b. Résoudre l’équation f (x) = 0 avec le code obtenu sur l’intervalle précédent à la
cinquième décimale près :
c. Ajouter un test de l’applicabilité de la Dichotomie avant le calcul :
d. Afficher le nombre d’itérations effectuées pour obtenir le résultat souhaité :
𝒃−𝒂
e. Tracer un graphe de l’erreur absolue e = en fonction des itérations :
𝟐