Discussion :

[ngossi]

Bonjour à tous,
Voici une fonction "test"(ci dessous) que j'ai écris sous MATLAB. En fait l'idée est d'avoir la même réponse entre le produit simple de MATLAB et mon programme que j'ai écris en dessous. Mais il se trouve que les résultats ne sont pas les mêmes.
SVP je souhaiterais avoir une aide de votre part pour savoir où exactement se trouve mon erreur
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
%Mon programme
function test
clc
 
text = 'Matrice source'
M = [29      95      60      60      64      52      61      63;...
123     85      119     114     118     103     95      103;...
65      79      82      106     115     95      97      73;...
71      79      99      82      119     122     96      106;...
106     86      90      148     132     132     116     162;...
93      99      79      103     123     113     115     117;...
141     155     88      102     108     95      118     103;...
134     160     93      104     97      102     105     128]
 
M = M*M %produit simple dans MATLAB
 
%programme similaire qui devrait produire le même résultat que celui de MATLAB
T = double([]);
for i = 1:8
    for j = 1:8
        T(i,j)=0;
        for k=1:8
            T(i,j) = T(i,j)+M(i,k)*M(k,j);
        end
    end
end
 
T
 
end

[Invité]

Bonjour,

Ton calcul est bon, mais il faudrait peut-être mettre la ligne

Code :

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

M = M*M %produit simple dans MATLAB

qui change M après le calcul de T ou mettre le résultat dans une autre variable que M.
Une petite remarque: l'initialisation de T peut se faire sous la forme T = zeros(8); tu n'as alors plus besoin de T(i,j)=0;

[ngossi]

En faisant T = zeros(8) en lieu et place de T(i,j)=0; j'obtiens un résultat satisfaisant et du coup, le résultats sont les même.
Merci beaucoup pour cette aide.

[Invité]

Non attention je parle de l'initialisation: à la place de T = double([]);. Ce qui donne:

Code :

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

1
2
3
4
5
6
7
8
T = zeros(8);
for i = 1:8
    for j = 1:8
        for k=1:8
            T(i,j) = T(i,j)+M(i,k)*M(k,j);
        end
    end
end

Voire en simplifiant un peu:

Code :

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

1
2
3
4
5
6
T = zeros(8);
for i = 1:8
    for j = 1:8
        T(i,j) = sum(M(i,:).*M(:,j).');
    end
end

Je te conseille de regarder ce point de la