Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Scribd Logo
Importer

Bienvenue sur Scribd !

  • 139 vues

    Méthodes Numériques-TP

    Transféré par

    Hind
    Ce compte rendu décrit trois méthodes itératives pour résoudre des systèmes linéaires: la méthode de Jacobi, la méthode de Gauss-Seidel et la méthode de relaxation. Il présente également deux méthodes numériques pour résoudre des équations différentielles: la méthode d'Euler explicite et implicite. Des codes sont fournis pour illustrer la mise en œuvre de ces méthodes.

    Droits d'auteur :

    © All Rights Reserved
    Téléchargez comme PDF, TXT ou lisez en ligne sur Scribd

    Méthodes Numériques-TP

    Transféré par

    Hind
    139 vues21 pages
    Ce compte rendu décrit trois méthodes itératives pour résoudre des systèmes linéaires: la méthode de Jacobi, la méthode de Gauss-Seidel et la méthode de relaxation. Il présente également deux méthodes numériques pour résoudre des équations différentielles: la méthode d'Euler explicite et implicite. Des codes sont fournis pour illustrer la mise en œuvre de ces méthodes.

    Description originale:

    Traveaux sur les méthodes itératives

    Copyright

    © © All Rights Reserved

    Formats disponibles

    PDF, TXT ou lisez en ligne sur Scribd

    Avez-vous trouvé ce document utile ?

    Ce contenu est-il inapproprié ?

    Ce compte rendu décrit trois méthodes itératives pour résoudre des systèmes linéaires: la méthode de Jacobi, la méthode de Gauss-Seidel et la méthode de relaxation. Il présente également deux méthodes numériques pour résoudre des équations différentielles: la méthode d'Euler explicite et implicite. Des codes sont fournis pour illustrer la mise en œuvre de ces méthodes.

    Droits d'auteur :

    © All Rights Reserved
    Téléchargez comme PDF, TXT ou lisez en ligne sur Scribd
    Télécharger au format pdf ou txt
    139 vues21 pages

    Méthodes Numériques-TP

    Transféré par

    Hind
    Ce compte rendu décrit trois méthodes itératives pour résoudre des systèmes linéaires: la méthode de Jacobi, la méthode de Gauss-Seidel et la méthode de relaxation. Il présente également deux méthodes numériques pour résoudre des équations différentielles: la méthode d'Euler explicite et implicite. Des codes sont fournis pour illustrer la mise en œuvre de ces méthodes.

    Droits d'auteur :

    © All Rights Reserved
    Téléchargez comme PDF, TXT ou lisez en ligne sur Scribd
    Télécharger au format pdf ou txt
    Vous êtes sur la page 1sur 21

    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

    Vous aimerez peut-être aussi