Discussion :

[Adjen]

Bonjour et meilleurs voeux pour cette année 2013 !

J'ai bricolé un bout de script qui me permet d'ouvrir simultanément plusieurs fichiers qui se présentent sous la forme d'un entête de quelques lignes suivi d'un tableau 9 colonnes et environ 10k lignes. Je conserve pour chaque fichier la deuxième colonne que je range dans une matrice que j'appelle p. J'obtiens alors une matrice 10k lignes et autant de colonnes que de fichiers ouverts.

J'aimerais maintenant obtenir l'indice des valeurs de p comprises entre 1.98 et 2.02 pour chaque fichier, donc pour chaque colonne, sachant que je ne connais absolument pas le nombre d'éléments se trouvant dans cette fenêtre et qu'il est susceptible d'être différent d'un fichier à l'autre.

Code :

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

1
2
3
for i = 1:numel(filename)
    ind_2bar(:,i) = find(2-2*0.01<p(:,i) & p(:,i)<2+2*0.01);
end

me renvoit :

Code :

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

Subscripted assignment dimension mismatch.

Simple problème de syntaxe j'imagine mais j'avoue bloqué complètement.

Merci pour votre aide.

[Invité]

Bonjour, et bonne année également

Chaque fichier ayant j'imagine un nombre différent de valeurs qui correspondent à ce critère, tu ne pourras pas stocker ces indices dans un simple tableau ayant un nombre de colonnes fixe.
Il te faut utiliser un tableau de cellules. Voir aussi ce lien.

[FLB]

Il y a bien autant de ligne dans tous les fichiers?

je ne connais absolument pas le nombre d'éléments se trouvant dans cette fenêtre et qu'il est susceptible d'être différent d'un fichier à l'autre

Je te conseille d'utiliser une structure pour stocker l'information au lieu d'un tableau :

Code :

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

ind_2bar{i} = find(2-2*0.01<p(:,i) & p(:,i)<2+2*0.01);

Seconde solution tu stockes non pas le résultat du find, mais celui du test logique

Code :

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

ind_2bar(:,i) = 2-2*0.01<p(:,i) & p(:,i)<2+2*0.01;

Car dans ce cas le résultat est toujours de même taille.

[Adjen]

Je n'avais pas pensé en effet qu'un simple tableau ne pouvait pas supporter des colonnes au nombre d'éléments différent. J'ai opté pour la solution du tableau de cellules du coup.

Code :

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

1
2
3
for i = 1:numel(filename)
    ind_2bar{i} = num2cell(find(2-2*0.01<p(:,i) & p(:,i)<2+2*0.01));
end

fonctionne impec !

Merci à vous deux !

[Invité]

Cela peut dépendre de ce que tu souhaites faire par la suite, mais l'utilisation de num2cell est inutile ici.

[Adjen]

En effet. Je l'avais ajouté suite au message d'erreur suivant :

Code :

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

Cell contents reference from a non-cell array object.

Certainement autre chose dans mon code que j'ai réglé sans même m'en rendre compte.

[Invité]

Tu devais avoir gardé en mémoire ta variable ind_2bar sous forme de simple tableau numérique.
Cell contents assignment to a non-cell array object.

[Adjen]

C'est fort possible oui. Je n'avais lancé qu'une partie du script à ce moment là. Merci encore !