Discussion :

[Xeher]

Bonjour,

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

Code :

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

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 :

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

struct(p1{1}).data

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

Code :

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

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

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

Code :

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

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

[Jerome Briot]

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.

[Xeher]

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 :

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

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 :

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

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 :

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

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.

[Jerome Briot]

Un peu trop de eval pour moi dans ce code

Voici une alternative :

Code :

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

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.

[Xeher]

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);

[Jerome Briot]

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.

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 :

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

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

[Xeher]

Merci, ça fonctionne !

Je commence à y voir un peu plus clair