Discussion :

[Wel Kol]

Bonjour,

J'ai un problème au niveau du stockage de mes résultats en utilisant la boucle for. A la fin du calcul, il n'y a que la dernier calcul qui est stocké. Ci -après mon code :

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
% A est une matrice contenant les composantes de vecteurs
 
nump=size(indice,1);
 
for k=2:nump
 
% Vecteur de reférence
 
v1=[A(indice(1),4),A(indice(1),5)];
 
 
v2 =[A(indice(k),4),A(indice(k),5)];
 
p=dot(vp1,vp2);
 
if p=o
   vp=v2;
else
    vo=v2;
end
 
end

Quelqu'un pourrait m'expliquer comment stocker mes résultats à chaque itération?
Merci d'avance

[thecrazydonut]

Par exemple

Code :

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

1
2
3
4
5
6
 
a=NaN(1,N);
for i=1:N
    a_temp=rand(1,1);
    a(i)=a_temp;
end

À chaque itération, un nouveau a_temp est calculé et est stocké ensuite dans la matrice a. Ne pas oublier la première étape de préallocation.

[Wel Kol]

Au fait j'ai une matrice A de dimension n*2. Chaque ligne correspond à un vecteur. J'aimerais calculer le produit scalaire du premier vecteur (première ligne) par rapport aux autres vecteurs ( 2ème ligne jusqu'à la nième ligne) et stocker les résultats. Comment vais je coder cela avec Matlab?

[Jerome Briot]

La réponse a été donnée par thecrazydonut. Montre nous ce que tu as codé (même si c'est faux)

Sinon, la ligne suivante :

Code :

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

v1 = [A(indice(1),4),A(indice(1),5)];

peut aussi s'écrire :

Code :

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

v1 = A(indice(1),4:5);

[Wel Kol]

Ci après mon code après quelques modifications faites suite à vos suggestions :

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
A est une matrice n*2
indice est une matrice contenant des indices
v22=zeros;
nump=size(indice,1);
 
for k=2:nump
 
% Vecteur de reférence
 
v1=[A(indice(1), 4:5)];
 
Les autres vecteurs
 
v2 =[A(indice(k),4:5)]; % 
 
% v22(k) =v2;   Quand je mets cela, pour pouvoir stocker à chaque itération, ça ne marche pas; Le message d'erreur suivant apparaît "In an assignment  A(I) = B, the number of elements in B and I must be the same. "
 
p=dot(vp1,vp2);
 
if p=o
   vp=v2;
else
    vo=v2;
end
 
end

A la fin du calcul, j'obtiens que les résultats du dernier calcul. Pourriez-vous me dire ce qui fait que je n'arrive pas à stocker mes résultats?

[Invité]

Bonjour,

v22(k) ne représente qu'un seul élément, alors que v2 =[A(indice(k),4:5)]; possède deux valeurs. Tu ne peux stocker deux valeurs à la place qu'une seule.

Concernant cette ligne : v22=zeros;, je te conseille de lire Qu'est-ce que la préallocation de mémoire ?.
Une fois fait... ton tableau v22 sera au final de taille nump x 2 : ta préallocation se fera donc ainsi :

Code :

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

v22 = zeros(nump, 2);

puis le stockage :

Code :

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

v22(k,:) = v2;

les deux points signifiant toutes la ligne, ici v22(k,1) et v22(k,2)