Octave/Matlab: l'Ajout de nouveaux éléments d'un vecteur

Ayant un vecteur x et j'ai du ajouter un élément (newElem) .

Quelle est la différence entre 

x(end+1) = newElem;

et 

x = [x newElem];

?

InformationsquelleAutor URL87 | 2013-04-24
  1. 91

    x(end+1) = newElem est un peu plus robuste.

    x = [x newElem] ne fonctionnera que si x est une ligne de vecteur, si elle est un vecteur colonne x = [x; newElem] doit être utilisé. x(end+1) = newElem, cependant, travaille pour de ligne et de colonne de vecteurs.

    En général cependant, la croissance de vecteurs doivent être évités. Si vous faites cela, beaucoup de choses, qu'il peut apporter à votre code vers le bas à un rampement. Pensez-y: la croissance d'une matrice s'agit d'allouer un nouvel espace, de tout recopier de plus, l'ajout de ce nouvel élément, et le nettoyage de l'ancien mess...Tout à fait une perte de temps si vous saviez la bonne taille à l'avance 🙂

    • Aussi, pour la deuxième méthode, x doit être initialisé en premier!
    • pas de problème! J'ai fait de même, sans doute en même temps. @Dan, c'est vrai, mais que la question visée "ayant un vecteur x (de taille n)", je suppose n à être non-nul, et le vecteur d'être initialisé déjà 🙂
    • Encore, à l'aide de end plutôt que de certaines variables n (qui pourrait être un global pour tout ce que vous savez!) est le plus universel, robuste-pas-coût moyen d'aller
    • x doit être initialisé pour les deux méthodes.
    • Je ne le pense pas, je n'ai pas de MATLAB pour vérifier, mais dans l'octave vous pouvez certainement le faire x(end+1)=... sans x cours d'initialisation et je suis presque certain que vous pouvez le faire dans MATLAB ainsi
    • donne Undefined fonction ou une variable "x".

    InformationsquelleAutor ThijsW
  2. 27

    Juste pour ajouter à @ThijsW réponse, il y a une augmentation considérable de la vitesse de l'avantage de la première méthode par rapport à la méthode de concaténation:

    big = 1e5;
tic;
x = rand(big,1);
toc

x = zeros(big,1);
tic;
for ii = 1:big
    x(ii) = rand;
end
toc

x = []; 
tic; 
for ii = 1:big
    x(end+1) = rand; 
end; 
toc 

x = []; 
tic; 
for ii = 1:big
    x = [x rand]; 
end; 
toc

   Elapsed time is 0.004611 seconds.
   Elapsed time is 0.016448 seconds.
   Elapsed time is 0.034107 seconds.
   Elapsed time is 12.341434 seconds.

    J'ai eu ces moments de course dans 2012b cependant quand j'ai couru le même code sur le même ordinateur dans matlab 2010a-je obtenir

    Elapsed time is 0.003044 seconds.
Elapsed time is 0.009947 seconds.
Elapsed time is 12.013875 seconds.
Elapsed time is 12.165593 seconds.

    Donc je suppose que la vitesse de l'avantage ne s'applique qu'aux plus récentes versions de Matlab

    • +1, Modifié pour ajouter l'évidence. Je vais le tester de nouveau sur un "vrai" PROCESSEUR (je suis sur cette merde pas fiable bon à rien APU chose maintenant...)
    • Là, tous mieux maintenant 🙂
    • de même pour moi, je reçois 0.028 pour la 3ème option et 8.909 pour la dernière
    • Je pense aussi que le JIT optimisation pour le x(end+1) cas est assez récente (R2012a ou alors...). Je vagualey souviens avoir lu quelque chose comme ça dans certains changenotes à smoe point. J'ai aussi obtenir des résultats très différents sur mon APU/Matlab 2010, mais je ne sais pas si c'est dû à l'Matlab version ou l'APU...
    • et ThiijsW voir ma récente de comparaison entre les anciennes et nouvelles matlab
    • Ah 🙂 OK, MAINTENANT c'est mieux 🙂

    InformationsquelleAutor Dan
  3. 4

    Comme mentionné précédemment, l'utilisation de x(end+1) = newElem a l'avantage qu'il permet de concaténer votre vecteur par un scalaire, peu importe si votre vecteur est transposée ou pas. Par conséquent, il est plus robuste pour l'ajout de scalaires.

    Cependant, ce ne doit pas être oublié, c'est que x = [x newElem] travaillera également lorsque vous essayez d'ajouter de multiples éléments à la fois. En outre, cette généralise un peu plus, naturellement, pour le cas où vous voulez concaténer des matrices. M = [M M1 M2 M3]

    Dans l'ensemble, si vous voulez une solution qui vous permet de concaténer votre vecteur x avec newElem qui peut ou peut ne pas être un scalaire, cela devrait faire l'affaire:

     x(end+(1:numel(newElem)))=newElem
    • Je pense que votre dernier exemple devrait être: x(end+1:fin+longueur(newElem)) = newElem
    • Merci de trouver le bug, j'ai mis à jour la réponse pour résoudre le problème.
    • Avec mon Matlab2011b, il y avait aussi une drastique (~50x) amélioration de la vitesse sur le vecteur de concaténation avec cette méthode contre la a=[a b] méthode.
    • le [a b], la méthode est la concaténation. Je suppose que vous dire que l'extension de la matrice est plus rapide que la concaténation?
    InformationsquelleAutor Dennis Jaheruddin