Université UMBB Faculté de Technologie
Département ISE 1ère année Master Electrotechnique
Module : TP Méthodes numériques Enseignant : N. Cherrat
TP° 04 Transport et méthode de différences finie
Exercice 1
On considère l’équation du transport suivant (EDP)
𝜕𝑢 𝜕𝑢
+ =0
𝜕𝑡 𝜕𝑥
{
2
𝑢(𝑥, 𝑡 = 0) = 𝑒 −𝑥
1) Montre que l’algorithme de la méthode de différence finie est écrire sous la forme :
𝑗+1 𝑗 𝑗 𝑗
𝑢𝑖 = 𝑢𝑖 − 𝑟(𝑢𝑖+1 − 𝑢𝑖 )
2) Calculer la valeur de 𝑟 pour Δ𝑥 = 0.4 et Δ𝑡 = 0.2
3) En utilisant MATLAB, calculer les valeurs de 𝑢01 , 𝑢11 , 𝑢21 , 𝑢31 et 𝑢41
4) En utilisant MATLAB, calculer les valeurs de 𝑢02 , 𝑢12 , 𝑢22 , 𝑢32 .
5) Conclure !
Exercice 2
On considère l’équation du transport suivant (EDP)
𝜕𝑢 𝜕𝑢
+5 = 0 ∀𝑥 ∈ [0 5], ∀𝑡 ∈ [0 5]
𝜕𝑡 𝜕𝑥
𝑢(𝑥, 𝑡 = 0) = 5
{ 𝑢(𝑥 = 5, 𝑡) = 1
Avec l’algorithme de la méthode de différence finie est écrire sous la forme :
𝑗+1 𝑗 𝑗 𝑗
𝑢𝑖 = 𝑢𝑖 − 5(𝑢𝑖+1 − 𝑢𝑖 )
1) Définît les commandes suivantes : linspace, length, zeros
2) Faire la discrétisation du temps en 6 valeurs en utilisant linspace.
3) Faire la discrétisation de l’axe 𝑥 en 6 valeurs en utilisant linspace.
4) Calculer les valeurs de 𝑢01 , 𝑢11 𝑢21 , 𝑢31 , et 𝑢41 / 𝑢02 , 𝑢12 𝑢22 , 𝑢32 et 𝑢42 / 𝑢03 , 𝑢13 𝑢23 , 𝑢33 et
𝑢43 . en utilisant les fonctions "for" et "while"
5) Tracer graphiquement ces solutions. Conclure !
Exercice supplémentaire :
Soit l’équation non linéaire : 𝑓(𝑥) = 𝑥 2 − 𝑥 − 2 = 0
1) Tracer le graphe 𝑦 = 𝑓(𝑥) sur un intervalle 𝑥 = −10: 0.001: 10 tel qu’il vous permet
de localiser la solution de l’équation.
Université UMBB Faculté de Technologie
Département ISE 1ère année Master Electrotechnique
Module : TP Méthodes numériques Enseignant : N. Cherrat
2) Il est à noter que, les solutions exactes de cette équation sont 𝑥1 = −1 et 𝑥2 = 2
et on veut trouver la première racine 𝑥1 de cette équation en utilisant :
a) La méthode de point fixe :
Trouver des formes possibles de la fonction 𝑔(𝑥) tel que 𝑓(𝑥) = 0 → 𝑔(𝑥) = 𝑥.
Quelle est la fonction qui vérifie le théorème de la méthode de point fixe, sur
l'intervalle [0, 3] ?
Ecrire un programme Matlab qui donne la solution de cette équation. Prendre
𝜀 = 0.01 et 𝑥0 = 0 puis 𝑥0 = 3. Conclure !
b) La méthode de Newton-Raphson
Ecrire un programme Matlab qui donne la solution de cette équation. Prendre
𝜀 = 0.01 et 𝑥0 = 2 puis 𝑥0 = 3. Conclure !
4) Comparer les résultats des différentes méthodes implémentées.