Selection d'une colonne dans une structure

Bonjour,

j'ai importé plusieurs fichiers p11_f_rep1_strain.txt à p16_f_rep1_strain.txt avec les commandes suivantes :

Code:

---

1 2 3 4 5 6

numfiles = 6; p1 = cell(1, numfiles); for k = 1:numfiles myfilename = sprintf('p1%d_f_rep1_strain.txt', k); p1{k} = importdata(myfilename); end

---

Je cherche maintenant a accéder aux données pour effectuer des opérations.
J'arrive à récupérer les données de la matrice du premier fichier p11_f_rep1_strain.txt avec la commande :

Code:

---

struct(p1{1}).data

---

Et j'arrive à récupérer les données à une colonne et une ligne donnée par la commande :

Code:

---

struct(p1{1}).data(1,4)

---

Mais je n'arrive pas à accéder à une colonne, j'ai tenté

Code:

---

struct(p1{1}).data(:,4)

---

Le problème est que j'avais testé avec des fichiers légers l'importation et je n'avais pas besoin d'utiliser la commande struc pour récupérer les données car matlab ne me créé pas de structure pour les fichiers légers. J'avais donc simplement fait p1{1}(:,4).
Mais cela ne fonctionne pas pour de gros fichiers de données.

J'avais également essayé d'importer un fichier à la fois. Cela marche car je n'ai alors qu'un "data" et je n'ai alors qu'à faire data(:,4).
Mais comme j'ai plusieurs fichier, je voulais réduire le nombre de lignes de mon modèle.

Quelqu'un a-t-il une idée ?

Merci d'avance

impordata n'est peut être pas la meilleure solution.

Montre nous le contenu d'un des fichiers texte (ou attache le au message) et dis nous quelles données tu souhaites lire à l'intérieur.

J'ai mis un fichier de données en pièce jointe.

Lorsque j'ai étudié ce fichier tout seul, je renommais les colonnes ainsi :

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17

val=data;   %% Les colonnes sont renommées pour faciliter la suite pt=val(:,1); tps=val(:,2); for i=1:14 eval(['def_a_R' num2str(i) '=val(:,4+6*(i-1));']); eval(['def_b_R' num2str(i) '=val(:,4+6*(i-1));']); eval(['def_c_R' num2str(i) '=val(:,5+6*(i-1));']); eval(['angle_R' num2str(i) '=val(:,6+6*(i-1));']); eval(['Spmax_R' num2str(i) '=val(:,7+6*(i-1));']); eval(['Seq_R' num2str(i) '=val(:,8+6*(i-1));']); end top_tour=val(:,87); torque=val(:,88); vitesse=val(:,89);

---

Puis je traçais les données des colonnes Seq :

Code:

---

1 2 3 4

for i=1:14 figure (i) eval(['plot(tps,Seq_R' num2str(i) ');']) end

---

Et je cherchais la valeur maximale des colonnes Seq :

Code:

---

1 2 3

for i=1:14 eval(['Seq_max_R' num2str(i) '=max(val(:,8+6*(i-1)));']) end

---

Mon objectif est d'importer plusieurs fichiers en même temps pour avoir le maximum des colonnes Seq pour chacun de mes fichiers.
Je souhaite ensuite tracer ces évolutions.

Un peu trop de eval pour moi dans ce code :aie:

Voici une alternative :

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27

fmt = ['%f' repmat(' %f',1,88)]; fid = fopen('p11_f_rep1_strain.txt','r');   data = textscan(fid,fmt,'collectoutput',true,'headerlines',5);   fclose(fid);   pt = data{1}(:,1); tps = data{1}(:,2);   i = 1:14; def_a_R = data{1}(:,4+6*(i-1)); def_b_R = data{1}(:,4+6*(i-1)); def_c_R = data{1}(:,5+6*(i-1)); angle_R = data{1}(:,6+6*(i-1)); Spmax_R = data{1}(:,7+6*(i-1)); Seq_R = data{1}(:,8+6*(i-1));   Seq_max_R = max(Seq_R);   top_tour = data{1}(:,87); torque = data{1}(:,88);   for i = 1:14 figure(i) plot(tps,Seq_R(:,i')); end

---

Par contre, dans ton cas, l'utilisation des variables intermédiaires et la quantité de valeurs à tracer par vecteur consomment pas mal de mémoire.

Désolé pour cette réponse tardive,
Merci pour l'alternative au eval, parc contre je n'arrive pas à trouver ce que signifie la première ligne avec la commande repmat.

Cela ne me dérange pas de consommer de la mémoire, je dois tracer ces données brutes pour détecter une éventuelle erreur.

Si je fais l'import de deux autres fichiers après le fclose, je vais pouvoir renommé les colonnes de chaque fichier par une boucle ?
pt{j}=data{j}(:,1);

Citation:

---

Envoyé par Xeher

je n'arrive pas à trouver ce que signifie la première ligne avec la commande repmat.

---

Pour un fichier à N colonnes avec des valeurs réelles, cela permet de répéter automatique N fois le format %f.

Citation:

---

Envoyé par Xeher

Si je fais l'import de deux autres fichiers après le fclose, je vais pouvoir renommé les colonnes de chaque fichier par une boucle ?
pt{j}=data{j}(:,1);

---

Non, tu ne dois pas changer d'indice pour data :

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22

numfiles = 6;   fmt = ['%f' repmat(' %f',1,88)];   for k = 1:numfiles   myfilename = sprintf('p1%d_f_rep1_strain.txt', k);   fid = fopen(myfilename,'r');   data = textscan(fid,fmt,'collectoutput',true,'headerlines',5);   fclose(fid);   pt{k} = data{1}(:,1); tps{k} = data{1}(:,2);   ...   clear data   end

---

Merci, ça fonctionne !

Je commence à y voir un peu plus clair :D