Bonjour, cliquez-ici pour vous inscrire et participer au forum.
  • Login:



+ Répondre à la discussion
Affichage des résultats 1 à 1 sur 1

Methode de Jacobi(Matlab)

  1. nubnab

    Date d'inscription
    mars 2008
    Messages
    1

    Methode de Jacobi(Matlab)

    bonjour

    je ne sais pas si je suis dans la bonne rubrique pour poster ce message mais bon je vais essayer


    je sollicite votre aide pour ecrire 2 fonctions sous matlab:
    j'aimerais aussi sans trop vous en demandez, de m'expliquer les solutions obtenus.

    1. ecrire une fonction [r]=Rhoj(A) qui donne le rayon spectrale* de la matrice J associée à la méthode itérative de Jacobi(On pourra utiliser la fonction de Matlab eigs)
    (*c'est le module de la plus grande valeur propre de J)

    Rappel de la méthode de Jacobi

    Pour une matrice A dont les éléments diagonaux sont nuls, on effectue la decomposition A=D - (E+F), avec D la matrice diagonale associée à A,

    E l'opposée de la matrice triangulaire strictement inferieure à A, et F l'opposée de la matrice triangulaire strictement superieure associée à A.

    on pose alors

    J= D-¹(E+F)=I- D-¹A

    (D-¹:matrice inverse de D)

    resoudre le systeme lineaire AX=b par la méthode de Jacobi, c'est trouver lalimite de la suite Xk definie par, Xo arbitraire et

    Xk+1=JXk + D-¹b

    la matrice J se calculant de façon rapide et exacte à l'aide de A et D


    2eme fonction


    ecrire une methode [X,err,it,prec]=[Jacobi(A,b,itm,tol) qui resout le systeme lineaire AX=b par la methode de Jacobi


    les arguments d'entrée et de sortie sont:

    input

    A matrice carré de taille n,n
    b matice de taille n,p
    itm ||AX - b ||/||b|| recherché
    (les deux en normes 2)



    output

    X matrice n,p solution

    err
    0 si inversion reussie
    1 si nb iteration depassé
    2 si rayon spectrale trop grand


    it nombre d'iterations effectuées

    prec ||AX - b ||/||b|| obtenue
    (les deux en normes 2)

    -----

     


    • Publicité







Sur le même thème :





 

Discussions similaires

  1. Identité de Jacobi
    Par Rasta--Rocket dans le forum Mathématiques du supérieur
    Réponses: 3
    Dernier message: 06/01/2008, 14h15
  2. MATLAB - Les intégrateurs ODE pour la méthode des lignes
    Par UterP dans le forum Mathématiques du supérieur
    Réponses: 0
    Dernier message: 13/05/2007, 12h37
  3. Méthode de quasi-newton avec matlab
    Par rogeg dans le forum Mathématiques du supérieur
    Réponses: 1
    Dernier message: 12/01/2007, 19h12
  4. Champ vectoriel - Identité de Jacobi
    Par Lévesque dans le forum Mathématiques du supérieur
    Réponses: 21
    Dernier message: 26/04/2006, 10h11
  5. avantage , limite des methodes de jacobi et de gauss jordan( urgent)
    Par fullcrum dans le forum Mathématiques du supérieur
    Réponses: 3
    Dernier message: 26/03/2006, 10h53