Université Abdelmalek Essaadi
Ecole Nationale des Sciences Appliquées de Tétouan
Département Sciences des Technologies Industrielles
Filière : Génie Mécanique
Compte rendu des travaux pratiques qui comportent sur :
Les Méthodes itératives
Réalisé par : Encadré par :
Bilal Yosra
M. Boukour Mustapha
El hafidi Hind
ANNEE UNIVERSITAIRE : 2021-2022 1
Table des matières
Remerciement ....................................................................................................................................3
Introduction........................................................................................................................................4
LA METHODE DE JACOBI .....................................................................................................................5
LA METHODE DE GAUSS-SEIDEL ...........................................................................................................6
LA METHODE DE RELAXATION ..............................................................................................................7
LA METHODE D ’EULER ........................................................................................................................8
PROBLEME DE CAUCHY .......................................................................................................................9
Objectif Général de la TP 1: ............................................................................................................... 10
I. CODE ASSOCIE A LA METHODE DE JACOBI : .................................................................................. 11
II. CODE ASSOCIE A LA METHODE DE GAUSS-SEIDEL : ........................................................................ 13
III. CODE ASSOCIE A LA METHODE DE RELAXATION : ....................................................................... 14
Conclusion .................................................................................................................................... 15
Objectif générale du TP 2 : ............................................................................................................... 16
Exercice 1 : La méthode d’Euler explicite ................................................................................. 17
Exercice 2 : Problème de Cauchy ............................................................................................. 18
Exercice 3 : Résolution numérique de l’équation différentielle ............................................... 19
ANNEE UNIVERSITAIRE : 2021-2022 2
Remerciement
Nous souhaitons adresser nos remerciements les plus sincères à Monsieur
M.BOUKOUR qui, en tant que professeur encadrant, s’est toujours montré à
l’écoute et très disponible tout au long de la réalisation de ces travaux pratiques.
Ainsi Nous lui remercions pour son aide et pour le temps qu’il a bien voulu nous
consacrer.
Nos remerciements s’étendent également aux membres du groupe pour leurs
efforts complices qui ont contribués au succès de ces Tps.
Ces travaux pratiques étaient une chance pour nous d’assimiler les données
théoriques et développer notre aspect pratique.
ANNEE UNIVERSITAIRE : 2021-2022 3
Introduction
ANNEE UNIVERSITAIRE : 2021-2022 4
LA METHODE DE JACOBI
ANNEE UNIVERSITAIRE : 2021-2022 5
LA METHODE DE GAUSS-SEIDEL
ANNEE UNIVERSITAIRE : 2021-2022 6
LA METHODE DE RELAXATION
ANNEE UNIVERSITAIRE : 2021-2022 7
LA METHODE D’EULER
ANNEE UNIVERSITAIRE : 2021-2022 8
PROBLEME DE CAUCHY
ANNEE UNIVERSITAIRE : 2021-2022 9
Objectif Général de la TP 1:
Sur les Méthodes itératives de résolution des systèmes linéaires
Ecrire les codes des trois m´méthodes itératives vues dans le cours avec le Tests
d’arrêt des itérations
Variables d’entrées: (n, A, b, X(0), ε)
Résultats obtenus: k, X(k)
Pour la m´méthode de relaxation de paramètre ω, les variables d’entrée sont :
Variables d’entrées: (ω, n, A, b, X(0), ε)
Sur des exemples simples, on comprend la convergence ou la divergence d’une telle
m´méthode itérative.
Sur des exemples extraits de la fiche des TDS, on compare numériquement la vitesse
de convergence des trois m´méthodes
On applique les m´ethodes it´eratives pour des matrices de tailles élèves creuses
quand trouve souvent dans l’exploitation de la m´ethode des ´el´ements finis.
A titre d’exemple :
ANNEE UNIVERSITAIRE : 2021-2022 10
I. CODE ASSOCIE A LA METHODE DE JACOBI :
Exécution du programme :
On introduit alors le programme qui permet à l’utilisateur d’entrer une matrice de son choix
ANNEE UNIVERSITAIRE : 2021-2022 11
format compact
disp('random number');
n = input('Donner les dimentions de la matrice\t');
A = input('Donner la matrice de la form of [A00 A01 ... ; A10 A11
....]\t');
disp(A);
%A = [10 1; 2 10];
b = input('Donner la matrice colonne B du solution de la form
[value1; value2]');
disp(b);
%b = [11 12];
X = input('Donner la matrice colonne X des solutions initiales de la
form [value1; value2]');
disp(X);
%X = [0 0];
tol = input('Donner la marge d''erreur epsilon');
[result, itr ] = Jacobi(n, A, b, X, tol, 0);
if (itr < 99)
disp('resultat final au sens de Jacobi');
disp(transpose(result));
disp('nombre des iteration');
disp(itr);
else
disp('Slution impossible au sens de Jacobi');
end
[result, itr] = gauss_Seiden(n, A, b, X, tol, 0);
if (itr < 99)
disp('resultat final au sens de Gauss-Seidel');
disp(transpose(result));
disp('nombre des iteration');
disp(itr);
else
disp('Slution impossible au sens de Gauss-Seidel');
end
[result, itr] = Relaxation(n, A, b, X, tol, 0);
if (itr < 99)
disp('resultat final au sens de Relaxation');
disp(transpose(result));
disp('nombre des iteration');
disp(itr);
else
disp('Slution impossible au sens de Relaxation');
end
ANNEE UNIVERSITAIRE : 2021-2022 12
II. CODE ASSOCIE A LA METHODE DE GAUSS-SEIDEL :
EXECUTION DU PROGRAMME
ANNEE UNIVERSITAIRE : 2021-2022 13
III. CODE ASSOCIE A LA METHODE DE RELAXATION :
n = input('Enter number of equations, n: ');
A = zeros(n,n+1);
x1 = zeros(1,n);
A=[4 2 3 8; 3 -5 2 -14; -2 3 8 27];
x1 = [0 0 0];
tol = input('Enter tolerance, tol: ');
m = input('Enter maximum number of iterations, m: ');
w = input('Enter the parameter w (omega): ');
k = 1;
while k <= m
err = 0;
for i = 1 : n
s = 0;
for j = 1 : n
s = s-A(i,j)*x1(j);
end
s = w*(s+A(i,n+1))/A(i,i);
if abs(s) > err
err = abs(s);
end
x1(i) = x1(i)+s;
end
if err <= tol
break;
else
k = k+1;
end
end
fprintf('The solution vector after %d iterations is :\n', k);
for i = 1 : n
fprintf(' %11.8f \n', x1(i));
end
EXECUTION DU PROGRAMME
ANNEE UNIVERSITAIRE : 2021-2022 14
Conclusion
Apres l’exécution pour plusieurs matrices, on constate qu’il se peut que la méthode de Jacobi
converge tandis que la méthode de Gauss Seidel diverge et pour certains exemples en trouve
l’inverse. C’est-à-dire gauss Seidel converge mais la méthode de Jacobi diverge.
ANNEE UNIVERSITAIRE : 2021-2022 15
Objectif générale du TP 2 :
Sur les Méthodes itératives de équations différentielles
Sur des exemples concrets des équations différentielles, on applique les méthodes d’Euler pour
trouver une solution approchée avec une précision demandée.
Exercice 1 : Tester la méthode d’Euler explicite sur
Pour trois valeurs de N, N = 3, N = 10, N = 50.
Tracer sur le même graphique les trois courbes ainsi que la courbe de la fonction exponentielle.
Exercice 2 : On considère le problème de Cauchy :
Exercice 3 : Résoudre numériquement l’équation différentielle
Par la méthode d’Euler explicite et implicite pour h = 0.1, h = 0.05 et h = 0.01 et
représenter ainsi les solutions approchées.
ANNEE UNIVERSITAIRE : 2021-2022 16
Exercice 1 : La méthode d’Euler explicite
clc, clear
a=0;
b=1;
y0=1;
f = @(x,y) y;
x=0:1/100:1 ;
[x1,y1] = EulerExplicite (a, b, y0, 3, f);
[x2,y2] = EulerExplicite (a, b, y0, 10, f);
[x3,y3] = EulerExplicite (a, b, y0, 50, f);
plot(x1,y1,'+y',x2,y2,'+b',x3,y3,'+r',x,exp(x),'y','LineWidth',2);
legend({'n=3','n=10','n=50','exp'},'Location','northwest');
function [x,y] = EulerExplicite (t0, t, y0, N, f)
h = (t - t0 ) / N;
x = t0:h:t ;
y = zeros (1,N);
y(1) = y0;
for i = 1 : N
y(i+1) = y(i) + h * f(x(i),y(i));
end
end
ANNEE UNIVERSITAIRE : 2021-2022 17
Exercice 2 : Problème de Cauchy
La solution théorique du problème (PC) :
∀ 𝑥 ∈ [𝑎, 𝑏] , 𝑦(𝑥) = 𝑦0 + ∫𝑓(𝑡, 𝑦(𝑡)) 𝑑𝑡.
Pour h = 0.1 et N = 30, les valeurs trouvées sont :
L’erreur globale est égale à 166.0525
ANNEE UNIVERSITAIRE : 2021-2022 18
Exercice 3 : Résolution numérique de l’équation différentielle
Programme du méthode d’Euler Explicite:
function [t,y] = Euler(f,tmin,tmax,Nint,y0) % Méthode d'Euler
h = (tmax-tmin)/Nint % h valeur du pas
t = linspace(tmin,tmax,Nint+1) ; % vecteur de t discrétisé t=[tmin,tmax]
y(1) = y0 ; % initialisation : y(1)=y(t0) = y0
for n = 2 : Nint+1
y(n) =y(n-1) + h*feval(f,t(n-1),y(n-1)) ;
% Calcul d'Euler
end
end
Les solutions approchées pour h= 0.1
ANNEE UNIVERSITAIRE : 2021-2022 19
Les solutions approchées pour h = 0,05
Code du méthode d’Euler Implicite:
function [x,y_i] = Myeulerimp(x0,xn,y0,h,f) x = x0 : h :xn;
n=length(x) y_i=zeros(1,n)+y0; y_i(1)=y0;
Erreur_i=1
while Erreur_i>=0.001 E_i=y_i(n);
for i=1:n-1 y_i(i+1)=y_i(i)+h*f(x(i+1),y_i(i+1)); end
Erreur_i = abs(E_i - y_i(n)); end
end
ANNEE UNIVERSITAIRE : 2021-2022 20
La partie commande :
ANNEE UNIVERSITAIRE : 2021-2022 21
