Comment générer toutes les paires de deux vecteurs dans MATLAB à l'aide de vectorisées code?

Plus d'une fois maintenant, j'ai besoin de générer de toutes les paires possibles de deux vecteurs dans MATLAB ce que je fais avec des boucles qui prennent juste quelques lignes de code, c'est à dire

vec1 = 1:4;
vec2 = 1:3;
i = 0;
pairs = zeros([4*3 2]);
for val1 = vec1
    for val2 = vec2
         i = i + 1;
         pairs(i,1) = val1;
         pairs(i,2) = val2;
    end
end

Génère ...

1 1
1 2
1 3
2 1
2 2
2 3
3 1
3 2
3 3
4 1 
4 2
4 3

Il doit y avoir une meilleure façon de faire ce qui est plus MATLAB'esque?

n.b. nchoosek ne fait pas de paires inversées ce dont j'ai besoin (c'est à dire 2 1 ainsi que 1 2), je ne peux pas juste l'inverse et ajouter le nchoosek sortie parce que le symétrique paires seront inclus deux fois.

InformationsquelleAutor Brendan | 2011-09-16
  1. 91

    Essayer

    [p,q] = meshgrid(vec1, vec2);
pairs = [p(:) q(:)];

    Voir le MESHGRID documentation. Bien que ce n'est pas exactement ce que la fonction est pour, mais si vous plissez les yeux lors qu'il est drôle, ce que vous demandez est exactement ce qu'il fait.

    • Que c'est génial. Je savais que les deux meshgrid et serialising de vecteurs 2d, mais n'a jamais envisagé de les utiliser de cette façon.
    InformationsquelleAutor Lambdageek
  2. 12

    Vous pouvez utiliser

    a = 1:4;
b = 1:3;
result = combvec(a,b);
result = result'
    • À noter que cela nécessite le Réseau de Neurones boîte à outils. Pour ceux qui en ont, il semble être la meilleure solution.
    • Comme cela mieux que la solution retenue, comme il le généralise à de multiples vecteurs.
    • Une réponse tardive je sais, mais il vaut la peine de futurs visiteurs sachant que c'est très lent par rapport à la solution retenue pour les grands vecteurs.
    InformationsquelleAutor efirvida
  3. 3

    Une autre solution pour collection:

    [idx2, idx1] = find(true(numel(vec2),numel(vec1)));
pairs = [reshape(vec1(idx1), [], 1), reshape(vec2(idx2), [], 1)];
    • +1: j'ai supprimé ma solution, j'aime vôtre mieux. J'ai pris la liberté de corriger ça, cependant, pour obtenir les combinaisons de vec1 et vec2.
    • Pas besoin de remodeler. Les Indices de find doit toujours être vecteurs colonnes, ainsi que vec1(idx1) etc.
    • Je pense que le remodelage est nécessaire. Étant donné que idx1 est un vecteur (n'a pas d'importance si de ligne ou de colonne), si vec1 est un vecteur ligne puis vec1(idx1) sera également un vecteur ligne. Il en va de même pour vec2.
    • OK, vous avez raison. Restaurée. Merci.
    InformationsquelleAutor yuk
  4. 3

    Vous pourriez le faire en reproduisant les matrices à l'aide de repmat, puis en tournant le résultat dans un vecteur colonne à l'aide de reshape.

    a = 1:4;
b = 1:3;
c = reshape( repmat(a, numel(b), 1), numel(a) * numel(b), 1 );
d = repmat(b(:), length(a), 1);
e = [c d]

e =

     1     1
     1     2
     1     3
     2     1
     2     2
     2     3
     3     1
     3     2
     3     3
     4     1
     4     2
     4     3

    Bien sûr, vous pouvez vous débarrasser de toutes les variables intermédiaires à partir de l'exemple ci-dessus.

    InformationsquelleAutor Praetorian
  5. 2

    Vous pouvez utiliser le bon vieux opérations matricielles, par exemple, dans 

    x = [3,2,1];
y = [11,22,33,44,55];
v = [(ones(length(y),1) * x)(:), (ones(length(x), 1) * y)'(:)]

    Edit: c'est l'Octave de la syntaxe, MATLAB ressemblera à ceci:

    x = [3,2,1];
y = [11,22,33,44,55];
A = ones(length(y),1) * x;
B = (ones(length(x), 1) * y)';
v = [A(:) B(:)]

    dans les deux cas, le résultat sera

    v =
 3    11
 3    22
 3    33
 3    44
 3    55
 2    11
 2    22
 2    33
 2    44
 2    55
 1    11
 1    22
 1    33
 1    44
 1    55
    • Ce n'est pas valable MATLAB syntaxe.
    • merci, mon mauvais. J'ai corrigé la syntaxe.
    InformationsquelleAutor jihor
  6. 1

    Ici un plus de MATLAB'esque façon de trouver des combinaisons. Celui-ci peut également être facilement étendu à plus de 2 vecteurs (et aussi non-combinaisons numériques):

    v1 =   1:  1:  3;
v2 =  11: 11: 44;
v3 = 111:111:555;

dimensions = cellfun(@numel, {v1,v2,v3});

[i1,i2,i3] = ind2sub(dimensions, 1:prod(dimensions));

combinations = [v1(i1); v2(i2); v3(i3)]'
    InformationsquelleAutor pAtrick
  7. 1

    À partir de la version R2015a, vous pouvez le faire à l'aide de repelem et repmat:

    >> vec1 = 1:4;
>> vec2 = 1:3;
>> pairs = [repelem(vec1(:), numel(vec2)) ...
            repmat(vec2(:), numel(vec1), 1)]

pairs =

     1     1
     1     2
     1     3
     2     1
     2     2
     2     3
     3     1
     3     2
     3     3
     4     1
     4     2
     4     3

    Ce type de solution évite supplémentaires variables intermédiaires requis par certains des autres solutions (comme celles basées sur les meshgrid) qui pourrait conduire à des problèmes de mémoire pour les plus grands vecteurs.

    InformationsquelleAutor gnovice
  8. 0

    Ce que vous cherchez est la produit cartésien

    cartprod est la fonction qui l'implémente. Vous pouvez le trouver dans le linéaire algèbre paquet de sorte que vous devez faire:

       >> pkg install -forge linear-algebra
   >> pkg load linear-algebra 
   >> sortrows(cartprod(1:4,1:3))                                            
    ans =                                                                                           
       1   1                                                                  
       1   2                                                                  
       1   3                                                                  
       2   1                                                                  
       2   2                                                                  
       2   3                                                                  
       3   1                                                                  
       3   2                                                                  
       3   3                                                                  
       4   1                                                                  
       4   2                                                                  
       4   3
    InformationsquelleAutor elviejo79