Indexation en vecteur

Question

Quelqu'un peut-il expliquer pourquoi dans l'exemple A le résultat est le vecteur 1x6 (ce qui est logique et était ce que je m'attendais) alors que dans l'exemple B le résultat est le vecteur 1x4?

Dans l'exemple B si je prédéfinis la taille de newvec pour être un vecteur 1x6, le résultat est le vecteur 1x6 correct. Juste ne pas comprendre ce qui se passe.

Exemple A

vec = [0 2 3 0 5 0]  %1x6 
newvec(vec == 0) = 1 %produces a 1 x 6 vector 

Exemple B

vec = [0 2 3 0 5 3]  %1 x 6 
newvec(vec == 0) = 1 %produces a 1 x 4 vector 

Answer 1

EDIT: Merci à Dan pour avoir souligné l'erreur.

vec == 0 produit une matrice logique. Maintenant, si la taille de newvec n'est pas définie alors newvec(vec==0) =1 créera toujours un nouveau vecteur de ligne newvec.

Les étapes suivantes décideront de la taille de cette newvec

temp1 = logical_matrix(:); 
trailing_zeros = number of zeros after the occurrence of last `true` value in temp1 
size(newvec) = [1, numel(logical_matrix) - trailing_zeros] ; 

Lorsque vous prédéfinir la taille de newvec alors newvec(v==0) = 1 ne peut pas changer la taille de newvec, il n'attribue 1 aux indices correspondants pour lesquels vec==0 est true.

Answer 2

Pourquoi, nous ne le saurons jamais. C'est un choix de conception par The MathWorks, et à mon humble avis, il viole en effet le principle of least surprise.

Il sent un peu de premature optimization (en termes d'utilisation de RAM) et/ou que les internes de MATLAB appellent find avant l'affectation, ou quelque chose de similaire. Mais c'est n'importe qui devine, vraiment.

Une façon de faire ce que vous voulez plus commune, est la suivante:

newvec = (vec == 0); 

Answer 3

Si newvec n'existe pas lorsque vous appelez la deuxième ligne, ne fait Matlab aussi grand qu'il a besoin de tenir les index vous définissez sur 1.

Qu'est-ce que vous êtes en train de faire est:

newvec([1 4 6]) = 1; ou

newvec([1 4]) = 1;

De même, si vec était en fait une etc. 2D/3D matrice, newvec sortira comme 1 x N, où N est quel que soit le nombre d'index le plus élevé, et ne conserve pas la forme de la matrice d'origine.

Donc, dans ces cas, vous devez soit:

1) Préallouer newvec être la taille de vec premier.

2) Est-ce newvec = (vec==0) au lieu