programme pivot de Gauss

Version imprimable

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:

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

16/11/2010, 22h56
Jean Dumoncel

Bonjour,

qu'as-tu fait? essayer de faire? Ou bloques-tu?

j'ai fais ça

Code:

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]===============]')

16/11/2010, 23h03
Jean Dumoncel

Et? Qu'est ce qui ne fonctionne pas? As-tu un message d'erreur?
16/11/2010, 23h10
azertysq123

tu peux me donné la solution svp
16/11/2010, 23h21
Jean Dumoncel

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.
16/11/2010, 23h28
azertysq123

oui ce ne fonctionne pas
et je ne sais pas c'est quoi erreur
17/11/2010, 18h39
FR119492

Salut!

Citation:

ne fonctionne pas

Qu'est-ce que ça veut dire? Ton ordinateur s'éteint-il tout seul? Y a-t-il de la fumée qui en sort? Y a-t-il un message qui s'affiche à l'écran, et, si oui, lequel?
Jean-Marc Blanc