Algorithme de décomposition de Cholesky

Version imprimable

14/03/2006, 12h59
mounstein

Algorithme de décomposition de Cholesky

Bonjour,

alors voilà j'ai un problème pour traduire mon algorithme de décomposition de Cholesky en pascal !!! :cry:

Le problème c'est comment traduire les indices de lignes et de colonnes de la matrice !!
Dois-je utiliser un tableau avec 2 dimensions ou y a t-il une autre astuce ??

Merci de répondre
14/03/2006, 13h02
M.Dlb

Re: cholesky

Citation:

Envoyé par mounstein

dois je utiliser un tableau avec 2 dimension ou y a t-il une autre astuce??

Non y'a pas d'autre astuce ;)
14/03/2006, 13h16
mounstein

:(

est ce que tu peux m'orienter stp?? :oops:
14/03/2006, 16h09
mounstein

reponse

est ce qu'il ya quelqu'un pour m'aider svp??

j'ai mis mon algorithme sous langage pascal!! la compilation est faite mais le resultat est etrange!! j'arrive pas a voire la matrice L que je devais voire en final!
14/03/2006, 17h11
Eric Sigoillot

Oui et ? On a oublié notre boule de cristal, donc sans code, sans résultat, sans rien en fait, on ne fera pas grand chose...

@++

oui c vrai

oui c'est vrai desolée voici le code!

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
program cholesky;
 
const
   n=3;
   m=3;
 
var
 
   a:array [1..n,1..m] of real;
   r:array [1..n,1..m] of real;
   i,j,k:integer;
   s1:real;
   s2:real;
 
begin
      s1:=0;
      s2:=0;
      writeln('entrez votre matrice A');
       for i:=1 to n do
        for j:=1 to m do
         readln(a[i,j]);
 
       repeat
         begin
 
           for k:=1 to i-1 do
           s1:=s1+sqrt(r[i,i]);
 
           r[i,i]:=sqrt(a[i,i]-s1);
 
           if r[i,i]=0 then
             writeln('matrice non reguliere')
 
           while (j>1) and (j<=n) do
             for k:=1 to i-1 do
             s2:=s2+(r[i,k]*r[j,k]);
 
             r[i,i]:=(a[i,j]-s2)/r[i,i];
 
           end;
 
      until i=n
end.

:oops:

14/03/2006, 19h27
Loceka

Faire des boucles, c'est bien, mais incrémenter/décrémenter les indices, c'est mieux ! :P

En gros ton repeat ne sert à rien si tu ne fait pas varier "i", de même ton while tourne dans le vide vu que "j" ne varie pas...
14/03/2006, 22h55
mounstein

heu..

j'ai pas bien compris ou dois je incrementer et decrementer?

votre extrait de code
Code:

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 repeat begin for k:=1 to i-1 do s1:=s1+sqrt(r[i,i]); r[i,i]:=sqrt(a[i,i]-s1); if r[i,i]=0 then writeln('matrice non reguliere') while (j>1) and (j<=n) do for k:=1 to i-1 do s2:=s2+(r[i,k]*r[j,k]); r[i,i]:=(a[i,j]-s2)/r[i,i]; end; until i=n
pose les problèmes:
1- ( sans importance mais tout de même)
avec repeat begin end until, la paire begin end est suppreflue
2- la boucle se résume à
repeat .... until i=n
mais i n'est pas incrementé dans la boucle et son initialisation est plus que douteuse car elle résulte de la fin d'un index for i=1 to n

Je n'ai pas regardéle code mais à prioris cela doit être quelquechose du type
Code:

1 2 3 4 5 6 i:=0; repeat inc(i); ..... { liste des choses à faire } until i=n;