Discussion :

[Cypette]

Bonjour,

J'ai grandement besoin de votre aide dans le cadre d'un travail où je suis bloqué

Mon but:

Créer une matrice 15x15 avec des valeurs T(i,j)

A la fin je dois avoir la somme des valeurs de chaque ligne égale aux différentes valeurs du vecteur O (somme des valeurs de la première ligne = O(1)- somme des valeurs de la 2eme ligne = O(2)). Je dois faire la même chose avec la somme des éléments de chaque colonne de la matrice et le vecteur D.

Mes formules pour calculer les T(i,j):

A(i) = 1/(somme de 1 à 15(B(j).D(j)))
B(j) = 1/(somme de 1 à 15(A(i).O(i)))

On commence par calculer A en posant B=1

Enfin, T(i,j) = A(i).B(j).O(i).D(j)

Mon code ne fonctionne pas sans que je sache pourquoi...

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
function []= OD()
tic
O=[191;43;110;78;49;372;398;193;442;137;199;53;72;90;527];
D=[363 46 133 79 70 236 194 90 310 127 168 70 57 130 881];
 
n=15;
Results=[];
T=ones(12);
A=ones(n,1);
B=ones(1,n);
A1=[];
 
A1=1/dot(B,D);
A=A1*ones(n,1);
B1=[];
B=[];
T1=zeros(n,1);
T2=zeros(1,n);
dif1=ones(n,1);
dif2=ones(1,n);
dif11=2;
dif22=2;
 
while (dif11>1 && dif22>1)
 
B1=1/dot(A,O);
B=B1*ones(1,n)
A1=1/dot(B,D);
A=A1*ones(n,1)
 
for i=1:n
    for j=1:n
        if i==j
          T(i,j)=0;
        else
        T(i,j)=A(i)*B(j)*O(i)*D(j);
        end
    end
end
 
for l=1:n
for k=1:n
    T1(l)=T(l,k)+T1(l);
end
end
 
for p=1:n
for q=1:n
    T2(p)=T(q,p)+T2(q);
end
end  
dif1=O-T1
dif2=D-T2
 
dif11=sum(dif1,1);
dif22=sum(dif2,2);
end
 
clearvars n A B A1 B1 O D
save('Results','T','T1','T2')
toc
end

Merci beaucoup pour votre aide

Cypette

[Gooby]

Bonjour,

Mon code ne fonctionne pas sans que je sache pourquoi...

Peux tu être un peu plus précise?

Je n'ai pas très bien compris ce passage,

Mes formules pour calculer les T(i,j):

A(i) = 1/(somme de 1 à 15(B(j).D(j)))
B(j) = 1/(somme de 1 à 15(A(i).O(i)))

Juste pour être bien sûr, les "." sont bien ici des multiplications?

[Cypette]

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
function []= OD()
tic
O=[191;43;110;78;49;372;398;193;442;137;199;53;72;90;527];
D=[363 46 133 79 70 236 194 90 310 127 168 70 57 130 881];
 
n=15;
Results=[];
A=ones(n,1);
B=ones(n,1);
T1=zeros(n,1);
T2=zeros(1,n);
dif1=2*ones(n,1);
dif2=2*ones(1,n);
 
 
while (dif1(1)>1 ||dif1(1)<-1|| dif1(2)<-1 || dif1(2)>1|| dif1(3)>1 || dif1(3)<-1|| dif1(4)>1 || dif1(4)<-1 || dif1(5)>1 || dif1(5)<-1|| dif1(6)>1 || dif1(6)<-1|| dif1(7)>1 || dif1(7)<-1|| dif1(8)>1 || dif1(8)<-1|| dif1(9)>1 || dif1(9)<-1|| dif1(10)>1 || dif1(10)<-1|| dif1(11)>1 || dif1(11)<-1|| dif1(12)>1 || dif1(12)<-1|| dif1(13)>1 || dif1(13)<-1|| dif1(14)>1 || dif1(14)<-1|| dif1(15)>1 || dif1(15)<-1 ||dif2(1)>1 || dif2(1)<-1|| dif2(2)>1 || dif2(2)<-1|| dif2(3)>1 || dif2(3)<-1|| dif2(4)>1 || dif2(4)<-1 || dif2(5)>1 || dif2(5)<-1|| dif2(6)>1 || dif2(6)<-1|| dif2(7)>1 || dif2(7)<-1|| dif2(8)>1 || dif2(8)<-1|| dif2(9)>1 || dif2(9)<-1|| dif2(10)>1 || dif2(10)<-1|| dif2(11)>1 || dif2(11)<-1|| dif2(12)>1 || dif2(12)<-1|| dif2(13)>1 || dif2(13)<-1|| dif2(14)>1 || dif2(14)<-1|| dif2(15)>1 || dif2(15)<-1)
B=B'
A=1./(B.*D);
A=A'
B=1./(A.*O);
 
for i=1:n
    for j=1:n
        if i==j
        T(i,j)=0;
        else
        T(i,j)=A(i)*B(j)*O(i)*D(j);
        end
    end
end
 
for l=1:n
for k=1:n
    T1(l)=T(l,k)+T1(l);
end
end
 
for p=1:n
for q=1:n
    T2(p)=T(q,p)+T2(q);
end
end  
dif1=O-T1;
dif2=D-T2;
 
end
 
clearvars n A B A1 B1 O D
save('Results','T','T1','T2')
toc
end

Voici une autre ébauche de code (cela ne marche toujours pas puisque le code fini par me sortir des valeurs qui sont des NaN)

En effet désolé du manque de précisions: les . sont des multiplier.

[Gooby]

Tu peux dans un premier temps essayer de lancer ton code pas à pas afin de déterminer quel ligne ne te renvoie pas le résultat attendu.

Il y a des parties dans ton code que je n'ai pas bien compris notamment celle ci:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
for p=1:n
for q=1:n
    T2(p)=T(q,p)+T2(q);
end
end  
dif1=O-T1
dif2=D-T2
 
dif11=sum(dif1,1);
dif22=sum(dif2,2);
end

Et pourquoi utiliser une boucle while?