IdentifiantMot de passe
Loading...
Mot de passe oublié ?Je m'inscris ! (gratuit)
Navigation

Inscrivez-vous gratuitement
pour pouvoir participer, suivre les réponses en temps réel, voter pour les messages, poser vos propres questions et recevoir la newsletter

MATLAB Discussion :

programme pivot de Gauss


Sujet :

MATLAB

Vue hybride

Message précédent Message précédent   Message suivant Message suivant
  1. #1
    Futur Membre du Club
    Inscrit en
    Février 2009
    Messages
    4
    Détails du profil
    Informations forums :
    Inscription : Février 2009
    Messages : 4
    Par défaut programme pivot de Gauss
    salam,

    je veux un programme élimination de gauss pivot total
    j'ai un algorithme et je sais pas comment le programmer sous matlab
    merci

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    11
    12
    13
    14
    15
    16
    17
    18
    19
    20
    21
    22
    23
    24
    25
    26
    27
    28
    29
    30
    31
    32
    33
    34
    35
    36
    37
    38
    39
    40
    41
    42
    43
    44
    45
    46
    47
    48
    49
    50
    51
    52
    A:tableau(20,20)reel, i,j,k,n:entier; s :reel;
    b,x: tableau(20)reel; rik : reel;
    AA,BB,amax:reel;i0:entier; (* Gauss Pivot Partiel*)
    j0:entier;Inc:tableau(20)reel; (*Gauss P.T*)
    Début 
    (*Lecture des données *)
        lire(n); lire(b(i),i=1,n); lire((a(i,j),j=1,n)i=1,n);
    (*Initialisation de Inc à  1, 2 , 3, ....., n *)
    pour i =1 à n faire inc(i) := i
                           fait;
        Pour k =1 à n-1
         faire  (*recherche du pivot *) 
             amax := a(k,k); i0 := k;j0:=k;
              pour i=k à n faire
                                       pour j =k à n faire si abs(a(i,j))>abs(amax)
                                                                           alors amax := a(i,j); j0:=j;i0:=i
                                                                      fsi
                                                              fait;
                                       fait;
     
          (*permutation des lignes i0 et k dans A et B *)
              pour j =K à n
                   faire aa := a(i0,j); a(i0,j) := a(k,j); 
                            a(k,j) := aa 
                    fait;
              bb := b(i0); 
              b(i0) := b(k); 
              b(k) := bb;
          (*permutation des colonnes j0 et k dans A et inc*)
              pour  i =k à n faire aa := a(i,j0); a(i,j0):=a(i,k);
                                              a(i,k):=aa
                                      fait;
              bb:=inc(j0); inc(j0):= inc(k); inc(k) := bb
    (* triangularisation de Gauss*)          
                pour i =k+1 à n faire rik := A(i,k)/A(k,k);
                                                   b(i) := b(i)-rik*b(k);   
                                                  pour j = k+1 à n
                                                     faire  A(i,j):= A(i,j)-rik*A(k,j)  fait
                                           fait
        fait;
    (* La remontée *) 
      x(inc(n))=b(n)/A(n,n);
    pour i =n-1 à 1 pas -1
      faire  s:=0; pour j =i+1 à n 
                              faire s := s+A(i,j)*x(inc(j))  fait;
                  x(inc(i)):= (b(i)-s)/A(i,i)
      fait;
     (* Impression des résultats *)
    pour i = 1 à n 
         faire écrire(i,x(i))  fait
     
    FIN

  2. #2
    Modérateur

    Homme Profil pro
    Ingénieur en calculs scientifiques
    Inscrit en
    Août 2007
    Messages
    4 639
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : Royaume-Uni

    Informations professionnelles :
    Activité : Ingénieur en calculs scientifiques

    Informations forums :
    Inscription : Août 2007
    Messages : 4 639
    Par défaut
    Bonjour,

    qu'as-tu fait? essayer de faire? Ou bloques-tu?
    Pour une bonne utilisation des balises code c'est ici!
    Petit guide du voyageur MATLABien : Le forum La faq Les tutoriels Les sources


    La nature est un livre écrit en langage mathématique. Galilée.

  3. #3
    Futur Membre du Club
    Inscrit en
    Février 2009
    Messages
    4
    Détails du profil
    Informations forums :
    Inscription : Février 2009
    Messages : 4
    Par défaut
    j'ai fais ça

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    11
    12
    13
    14
    15
    16
    17
    18
    19
    20
    21
    22
    23
    24
    25
    26
    27
    28
    29
    30
    31
    32
    33
    34
    35
    36
    37
    38
    39
    40
    41
    42
    43
    44
    45
    46
    47
    48
    49
    50
    51
    52
    53
    54
    55
    56
    57
    58
    59
    60
    61
    62
    63
    64
    65
    66
    67
    68
    69
    70
    71
    72
    73
    74
    75
    76
    77
    78
    79
    80
    81
    82
    83
    %A:tableau(20,20)reel,
    %i,j,k,n:entier;
    %s :reel;
    %b,x: tableau(20)reel;
    %rik : reel;
    %AA,BB,amax:reel;
    %i0:entier; (* Gauss Pivot Partiel*)
    %j0:entier;
    %Inc:tableau(20)reel; (*Gauss P.T*)
    clear all
    clc
    %(*Lecture des donnees *)
       a=input('entre la matrice A: ');
       b=input('entre le vecteur B: ');
    %(*Initialisation de Inc a 1, 2 , 3, ....., n *)
    n=length(a);
    for i=1:n
        inc(i) = i;
    end
        for k=1:n-1
            %("rechercher du pivot")
            amax =a(k,k);
            i0=k;
            j0=k;
            for i=k:n
                for j=k:n 
                    if abs(a(i,j))>abs(amax)
                        amax=a(i,j);
                        j0=j;
                        i0=i;
                    end
                end
            end
     
            %      (*permutation des lignes i0 et k dans A et B *)
            for j=k:n
                aa=a(i0,j);
                a(i0,j)=a(k,j);
                a(k,j)=aa;
            end
                     bb = b(i0); 
              b(i0) = b(k); 
              b(k) = bb;
     
            % (*permutation des colonnes j0 et k dans A et inc*)
            for i=k:n 
                aa=a(i,j0);
                a(i,j0)=a(i,k);
                a(i,k)=aa;
            end
            bb=inc(j0);
            inc(j0)=inc(k);
            inc(k)=bb;
            %(* triangularisation de Gauss*)          
            for  i =k+1:n 
                rik=a(i,k)/a(k,k);
                b(i)=b(i)-rik*b(k);
                for j=k+1:n
                    a(i,j)=a(i,j)-rik*a(k,j);
                end
            end
        end
        %(* La remontée *) 
        x(inc(n))=b(n)/a(n,n);
        for i=n-1:1
            s=0
            for j=i+1:n
                s=s-a(i,j)*x(inc(j))  
            end
        end
        % (* Impression des résultats *)
        for i=1:n
            disp(x(1))
     
        end
     
             disp('[==========LES MATRICES(A)===========]')
    a
    disp('[==========LES VECTEURS(B)===========]')
     
    disp('[==========LES SOLUTIONS(X)==========]')
    x
    disp('[================[Fin]===============]')

  4. #4
    Modérateur

    Homme Profil pro
    Ingénieur en calculs scientifiques
    Inscrit en
    Août 2007
    Messages
    4 639
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : Royaume-Uni

    Informations professionnelles :
    Activité : Ingénieur en calculs scientifiques

    Informations forums :
    Inscription : Août 2007
    Messages : 4 639
    Par défaut
    Et? Qu'est ce qui ne fonctionne pas? As-tu un message d'erreur?
    Pour une bonne utilisation des balises code c'est ici!
    Petit guide du voyageur MATLABien : Le forum La faq Les tutoriels Les sources


    La nature est un livre écrit en langage mathématique. Galilée.

  5. #5
    Futur Membre du Club
    Inscrit en
    Février 2009
    Messages
    4
    Détails du profil
    Informations forums :
    Inscription : Février 2009
    Messages : 4
    Par défaut
    tu peux me donné la solution svp

  6. #6
    Modérateur

    Homme Profil pro
    Ingénieur en calculs scientifiques
    Inscrit en
    Août 2007
    Messages
    4 639
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : Royaume-Uni

    Informations professionnelles :
    Activité : Ingénieur en calculs scientifiques

    Informations forums :
    Inscription : Août 2007
    Messages : 4 639
    Par défaut
    Nous ne sommes pas là pour faire les éxos à ta place. Donc tu nous a montré ce que tu as fait, c'est bien, mais pour avancer, il faut que tu nous dises pourquoi tu penses que ton programme ne fonctionne pas.
    Pour une bonne utilisation des balises code c'est ici!
    Petit guide du voyageur MATLABien : Le forum La faq Les tutoriels Les sources


    La nature est un livre écrit en langage mathématique. Galilée.

Discussions similaires

  1. [Débutant] Code Pivot de Gauss
    Par didier35 dans le forum MATLAB
    Réponses: 1
    Dernier message: 27/03/2014, 22h27
  2. Résolution pivot de Gauss
    Par haraigo dans le forum C
    Réponses: 7
    Dernier message: 16/05/2008, 23h08
  3. Explication de pivot de gauss
    Par bilou_2007 dans le forum Mathématiques
    Réponses: 6
    Dernier message: 01/03/2007, 22h33
  4. [Pivot de Gauss] probleme si pivot nul
    Par jmjmjm dans le forum Mathématiques
    Réponses: 33
    Dernier message: 02/02/2007, 15h47
  5. [LG]Matrice et pivot de Gauss
    Par Loopingus dans le forum Langage
    Réponses: 3
    Dernier message: 16/03/2005, 17h26

Partager

Partager
  • Envoyer la discussion sur Viadeo
  • Envoyer la discussion sur Twitter
  • Envoyer la discussion sur Google
  • Envoyer la discussion sur Facebook
  • Envoyer la discussion sur Digg
  • Envoyer la discussion sur Delicious
  • Envoyer la discussion sur MySpace
  • Envoyer la discussion sur Yahoo