Discussion :

[linda05]

Bonjour à tous;
J’ai deux questions :

1. Si j’ai un vecteur

Code :

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

 V=[2 ; 5 ;7]

,je cherche un vecteur VV tel que :

Code :

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

VV=[2; 2 ; 2 ;5 ; 5 ; 5 ; 7 ; 7 ; 7] ?

2. Si j’ai une matrice A carrée d’ordre 5:

Code :

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

A=[33 44 55 22 11;11 22 44 66 44;33 99 88 11 77;22 44 33 44 77 ;44 22 33 11 11]

et deux matrices B et C contenant les indices lignes et colonnes de A :

Code :

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

1
2
 B=[3 2 1 ;2 3 1 ;4 2 1 ;2 4 1 ;4 3 1 ;3 4 1]
C=[5 4 ;5 3 ;5 2]

Je cherche un vecteur X tel que :

Code :

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

X=[A(C(1,end),B(1,1)) ; A(C(1,end),B(2,1)) ; A(C(2,end),B(3,1)) ; A(C(2,end),B(4,1)) ; A(C(3,end),B(5,1)) ; A(C(3,end),B(6,1)) ]

ie:

Code :

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

 X=[33;44;11;99;66;44]??

merci

[Gooby]

Bonjour,

Pour être franc, je n'ai pas compris ce que tu souhaites/ce que tu attends de nous. Quelles sont les questions exactement?

Code :

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

VV = [V(1);V(1);V(1);V(2);V(2);V(2);V(3);V(3);V(3)

?

[linda05]

Merci pour votre reponse:

VV = [V(1);V(1);V(1);V(2);V(2);V(2);V(3);V(3);V(3)

mais si je doubler les éléments par exemple 10 fois ,donc je trouve ton idée long non?

ma question c'est comment trouver les vecteurs VV et X d'une façon optimisée?
merci

[Gooby]

Ok, dans ce cas tu peux utiliser une boucle pour VV

Code :

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

1
2
3
4
 
for i = 1:length(V)
    VV(1+(i-1)*3:i*3,1) = V(i);
end

Je regarde plus tard pour la seconde question.

PS: tu peux entrer en paramètre le nombre de répétitions de ta variable =>

Code :

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

1
2
3
4
5
parameter = 3;
for i = 1:length(V)
    VV(1+(i-1)*parameter:i*parameter,1) = V(i);
end

[linda05]

Merci,
mais est ce qu'il possible sans boucle?

[Dombrai]

salut,

Code :

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

1
2
3
 
Temp = repmat(V',nb_repetitions,1);
VV = Temp(:);

[Invité]

Bonjour,

Réplication entrelacée

[linda05]

Merci Dombri;
salut Winjerome et merci bien:
donc d’après vous j'ai ecrit ce code pour X:

Code :

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

1
2
3
rows=kron(C( :,end),ones(2,1))
cols=B(:,1)
X=A( sub2ind(size(A),rows,cols) )

est ce que vous avez mieux que ça?

J’ai une autre question:
Je cherche un moyen de récupérer les indices de lignes de C et B corresponds à la première valeur max de X
Donc ds l’exemple :
La valeur max est :99
Les indices de lignes corresponds :
Pour B c’est 4
Pour C c’est 2 ??


merci

[Paul-Ceria]

1 :

Code :

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

1
2
V=[2 ; 5 ;7];
VV=kron(V, ones(3,1))

2 :

Code :

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

1
2
e=[kron(C(:,end),ones(2,1)),B(:,1)];
X=diag(A(e(:,1),e(:,2)))

3 :

Code :

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

1
2
3
[c,b]=find(A==max(A(:)));
indice_C=find(C(:,end)==c)
indice_B=find(B(:,1)==b)