Discussion :

[jonnyx]

Bonjour,
Je viens vers vous car j’ai besoin d’importer un fichier texte dans Matlab et je suis un peu perdu avec toutes les fonctions qui existent… Je pense pouvoir réussir à parvenir à mes fins tout seul mais j’ai peur de faire du ‘bricolage’ en n’utilisant pas la bonne fonction et me retrouver avec beaucoup de lignes de code pour pas grand-chose !
J’ai donc un fichier texte qui a :
‘n’ lignes de données inutiles
‘N’ lignes de données utiles
‘n2’ lignes de données inutiles
‘N’ lignes de données utiles
‘n3 ‘….
…
Je souhaiterais avoir un fichier sans espace et sans saut de ligne sous la forme
‘n’ lignes de données utiles
‘n’ lignes de données utiles
….
Le nombre de lignes de données inutiles varie d’un bloc à l’autre (n1 , n2), je ne peux donc pas le passer en paramètre (type Headerline par exemple).
En revanche : je connais le nombre N de lignes utiles pour chaque bloc de données utiles. Ce nombre est spécifié dans l’entête (données inutiles) et je l’ai récupéré grâce à un bout de code, me disant que ça pourrait servir.
Je sais que la première ligne chaque bloc de données utiles commence par ‘1’ (premier point de mesure) et la dernière ligne commence par ‘N’ .
Je vous serais vraiment très reconnaissant de votre aide, je ne suis pas du tout à l’aise avec la lecture/modification de fichiers.
Ci-joints deux fichiers ‘file.txt’ fichier d’origine et ‘file2.txt’ fichier que je souhaiterais avoir.

Je vous remercie

[Jerome Briot]

Une solution :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
fid = fopen('file.txt', 'r');
 
X = textscan(fid, '%s', 'delimiter', '\n');
 
fclose(fid);
 
idx = cellfun(@(x) ~isempty(x) && isstrprop(x(1), 'digit') , X{1}, 'UniformOutput', true);
 
M = cellfun(@(x) str2num(x), X{1}(idx));
 
M = reshape([M{:}], 4, []).';

Ce qui donne :

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
 
M =
 
    1.0000         0    1.0000    1.0000
    2.0000    0.0000    1.0000    1.0000
    3.0000    0.0000    1.0000    1.0000
    4.0000    0.0000    1.0000    1.0000
    5.0000    0.0001    1.0000    1.0000
    6.0000    0.0001    1.0000    1.0000
    7.0000    0.0001    1.0000    1.0000
    8.0000    0.0001    1.0000    1.0000
    9.0000    0.0001    1.0000    1.0000
   10.0000    0.0001    1.0000    1.0000
    1.0000         0    1.0000    1.0000
    2.0000    0.0000    1.0000    1.0000
    3.0000    0.0000    1.0000    1.0000
    4.0000    0.0000    1.0000    1.0000
    5.0000    0.0001    1.0000    1.0000
    6.0000    0.0001    1.0000    1.0000
    7.0000    0.0001    1.0000    1.0000
    8.0000    0.0001    1.0000    1.0000
    9.0000    0.0001    1.0000    1.0000
   10.0000    0.0001    1.0000    1.0000

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
>> M(:,2)
 
ans =
 
   1.0e-03 *
 
         0
    0.0153
    0.0305
    0.0458
    0.0610
    0.0763
    0.0916
    0.1068
    0.1221
    0.1373
         0
    0.0153
    0.0305
    0.0458
    0.0610
    0.0763
    0.0916
    0.1068
    0.1221
    0.1373

Il suffit ensuite de chercher les valeurs égalent à 1 dans la première colonne de M pour séparer chaque bloc :

Code :

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

1
2
3
4
5
6
>> find(M(:,1)==1)
 
ans =
 
     1
    11