UNIVERSITE
ABDELMALEK
ESSAADI
ﺍاﻟﺳﻌﺩدﻱي ﺍاﻟﻣﺎﻟﻙك ﻋﺑﺩد ﺟﺎﻣﻌﺔ
FACULTE
DES
SCIENCES
ﺍاﻟﻌﻠﻭوﻡم
ﻛﻠﻳﯾﺔ
TETOUAN
ﻥن ﺗﻁطﻭوﺍا
TSE- M1- Semestre 2
Année : 2019-2020
TD Méthodes Numériques
Série N°2
Exercice 1 :
On veut calculer l’unique racine positive r de l’équation f(x) = 0 par la méthode de point fixe
où
f(x)=ex −x−2 à l’intervalle [1,2]
1. Trouver l’équation de la forme g(x) = x qui se ramène à g(r) = r
2. Calculer les itérés (4 décimal près) en partant depuis x0 = 1
3. Écrire l’algorithme de point fixe, ainsi que le code Matlab correspondant.
Exercice 2 :
x
On considère l’équation : 3cosx = 10 .
1. En appliquant la fonction log (décimal) de chaque côté, écrire cette équation sous la
forme g(x) = x.
2. Vérifier que g est contractante sur l’intervalle [0, π/6], de rapport k = log(e)/ 3
3. Calculer les 3 premières itérées de x0 = 0.2, en précisant à chaque étape une majoration
de l’erreur commise.
Exercice 3 :
Calculer les zéros des fonctions suivantes, à l'aide des deux algorithmes Point fixe et Newton :
1. x3− x − 1 = 0
2. x − 0, 2 sin x − 0, 5 = 0
3. x−tanx=0
4. e-x − cos x = 0.
Donner le code Matlab correspondant, puis comparer.
Exercice 4 :
Écrire la méthode de Newton pour la résolution du système suivant :
x2+ 2xy = 0,
xy + 1 = 0.
Exercice 5 :
On désire résoudre le système d’équations non linéaires ci-dessous par la méthode de Newton-
Raphson :
f(x,y)=x2+y2-3=0
g(x,y)=xy-1=0
PR.
ZUGARI
ASMAA
1
1. Exprimer la matrice Jacobienne J(x,y) en fonction de f et g.
2. Donner le code Matlab de résolution de ce système, en prenant comme vecteur initial
(x0, y0) = (0.5, 1.5). En déduire le vecteur solution (x*, y*), ainsi que le nombre
d’itérations.
Exercice 6 :
Résoudre le problème de Cauchy suivant :
𝒅𝒚
= −𝒚 𝒕 − 𝟓𝒆*𝒕 − 𝒔𝒊𝒏 𝟓𝒕 ,
𝒚 𝟎 = 𝟏
𝒅𝒙
Par la méthode d’Euler implicite, combinée avec la méthode de Newton-Raphson, pour
𝑡𝜖 0, 3 , discrétisé en N=100 points équidistants.
1. Donner le schéma d’Euler implicite pour cette équation différentielle.
2. En posant x=yn+1 et x0=yn , donner l’équation non linéaire G(x) à résoudre par la
méthode de Newton-Raphson.
3. Donner le code Matlab permettant de calculer y(t). On fixera epsilon =1e-6
4. Tracer y(t) sous Matlab pout 𝑡𝜖 0, 3 ,
PR.
ZUGARI
ASMAA
2