Le document présente trois exercices sur la résolution numérique d'équations non linéaires. Le premier exercice propose trois méthodes de point fixe pour résoudre une équation et les valide avec Matlab. Le deuxième exercice détermine les points fixes et leur domaine de convergence pour une fonction. Le troisième exercice dérive et implémente l'algorithme de Newton pour le calcul de racines carrées.
0 évaluation0% ont trouvé ce document utile (0 vote)
150 vues6 pages
Le document présente trois exercices sur la résolution numérique d'équations non linéaires. Le premier exercice propose trois méthodes de point fixe pour résoudre une équation et les valide avec Matlab. Le deuxième exercice détermine les points fixes et leur domaine de convergence pour une fonction. Le troisième exercice dérive et implémente l'algorithme de Newton pour le calcul de racines carrées.
Le document présente trois exercices sur la résolution numérique d'équations non linéaires. Le premier exercice propose trois méthodes de point fixe pour résoudre une équation et les valide avec Matlab. Le deuxième exercice détermine les points fixes et leur domaine de convergence pour une fonction. Le troisième exercice dérive et implémente l'algorithme de Newton pour le calcul de racines carrées.
Le document présente trois exercices sur la résolution numérique d'équations non linéaires. Le premier exercice propose trois méthodes de point fixe pour résoudre une équation et les valide avec Matlab. Le deuxième exercice détermine les points fixes et leur domaine de convergence pour une fonction. Le troisième exercice dérive et implémente l'algorithme de Newton pour le calcul de racines carrées.
Prof: BAKKALI Abderrahmane TD3: Equations non linéaires Exercice 1:
On cherche à résoudre l'équation
qui possède les deux racines et ainsi une troisième racine située près de . On vous propose les méthodes des points fixes suivantes pour obtenir .
1) Lesquelles, parmi ces trois méthodes des points fixes, sont susceptibles de converger vers ? (Ne pas calculer les itérations).
2) Retrouver les racines de l'équation par l'utilisation de la fonction de Matlab
pointfixe(g,x0, tol, nmax) (voir dernier chapitre du cours).
Exercice 2:
1) Obtenir tous les points fixes de la fonction
2) Déterminer pour chaque point fixe trouvé en (1) les valeurs de pour lesquelles ces points fixes sont attractifs (c-à-d la méthode de point fixe susceptible de converger).
Exercice 3:
1) A l'aide de la méthode de Newton-Raphson, montrer comment on obtient l'algorithme
pour le calcul de .
2) Obtenir un algorithme similaire pour calculer .
Filière SMP-S6-parcours mécanique Année universitaire: 2020/2021 Faculté des Sciences de Tétouan Prof: BAKKALI Abderrahmane 3) On pose N=2. Pour calculer , écrire un programme sous Matlab, utilisant l'algorithme donné en (1), avec critère de convergence . Vérifier bien que le programme converge vers . Filière SMP-S6-parcours mécanique Année universitaire: 2020/2021 Faculté des Sciences de Tétouan Prof: BAKKALI Abderrahmane Solution: TD3 Exercice 1:
1) La condition de convergence correspond à .
on a
La condition est vérifiée pour les trois méthodes des points fixes. Donc les trois méthodes des points fixes proposées sont susceptibles de converger vers . 2) >> g=@(x) -sqrt(exp(x)./3); >> x0=-0.43;tol=10^-15; >> [zero,erreur,niter] = pointfixe(g,x0,tol,50)
zero =
-4.5896e-001
erreur =
3.3307e-016
niter =
23
On voit bien que la fonction de Matlab utilisant la méthode de point fixe, pour la fonction , converge vers la première racine de . Pour converger vers on a choisit proche de la première racine.
Retrouver les 2 autres racines en changeant le point initial. Refaire la même chose pour et
Exercice 2:
1) on a
ce qui implique
donc, soit ou
les points fixes de la fonction sont:
Filière SMP-S6-parcours mécanique Année universitaire: 2020/2021 Faculté des Sciences de Tétouan Prof: BAKKALI Abderrahmane et
2) la condition de convergence est
pour ,
En appliquant la condition de convergence
ce qui donne .
Pour que la méthode de point fixe converge vers 0, doit être compris dans ]-1 1[.
pour ,
En appliquant la condition de convergence
Ce qui donne , d'où
Pour que la méthode de point fixe converge vers , doit être compris dans ]1 3[.
Exercice 3:
1) Pour calculer , on va considérer une fonction dont la racine est
Il suffit de poser .
La méthode Newton-Raphson recherchant le zéro de est donnée par
En remplaçant et dans l'équation précédente, et après des
simplifications on trouve
2) Il suffit de poser et refaire les mêmes étapes faites dans la question (1) pour arriver à Filière SMP-S6-parcours mécanique Année universitaire: 2020/2021 Faculté des Sciences de Tétouan Prof: BAKKALI Abderrahmane 3) Le programme sous Matlab: function [zeros,erreur,niter] = newton(N,x0,tol,nmax) zeros=x0; for niter=1:nmax x=zeros; zeros=(x+N/x)/2; erreur=abs(zeros-x); if erreur<tol return end end disp('il n y a de convergence avec le nombre d''iterration indiqué')
Ce programme est écrit dans un fichier m-file et enregistré portant le nom de la fonction (newton).
Les entrées de la fonction sont:
N: Le nombre dont on veut calculer la racine au carré.
x0: Initiation de programme, doit être différent de zéros.
tol: Critère d'arrêt ou erreur toléré.
nmax: Nombre d'itération.
Les arguments de sortie de la fonction sont:
zéros: La racine au carré de N.
erreur: Erreur estimée.
niter: Nombre d'itération effectué jusqu'à la convergence.
Exécution de la fonction:
La fonction est exécutée dans la fenêtre window command:
