Extraction de données dans un fichier texte

Bonjour à tous,

Merci par avance de l'aide que vous pourrez m'offrir.

J'aimerais extraire des données d'un fichier texte de 150 000 lignes. Je sais c'est une question posée souvent.

J'ai déjà réussi à extraire certaines informations, les informations suivantes se reproduise plusieurs avec des valeurs différentes:

Code:

---

1 2 3 4

EVODRV: SAVE INFORMATION RELATED TO TIME 8.6400E-04 E+8 S ( 1.0000E+00 DAY) ON LCM DIRECTORY NAMED 'DEPL-DAT0002'. -> FINAL BURNUP AT TIME = 1.000000E+00 DAYS FUEL BURNUP = 3.309378E+01 MW*D/TONNE NEUTRON EXPOSURE = 2.273049E-02 N/KB

---

Dans cette partie je voulais récupérer les valeurs de BURNUP et de TIME, j'ai donc exécuté le script suivant:

Code:

---

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

fid=fopen('pincell_uo2_burnup.result','rt'); tline = fgetl(fid); BURNUP = [0]; TIME = [0]; while ischar(tline) a = strread(tline,'%s'); if ~isempty(a) && strcmp(a{1},'FUEL') && strcmp(a{2},'BURNUP') BURNUP = [BURNUP; str2double(a{4})]; end   if ~isempty(a) && strcmp(a{1},'->') && strcmp(a{3},'BURNUP') TIME = [TIME; str2double(a{7})]; end tline = fgetl(fid); end fclose(fid); BURNUP TIME

---

J'ai réussi grâce à d'autres posts sur ce même forum. Je vous en remercie d'ailleurs.


---Dorénavant j'aimerais récupérer la valeur du KINF:

Code:

---

++ TRACKING CALLED= 6 TIMES FINAL KINF=1.362294E+00 FINAL KEFF=1.000000E+00 B2= 6.16706E-03 PRECISION= 1.31E-07

---

Je ne peux pas procéder de la même manière car il n'y a pas d'espace entre le KINF et sa valeur.

---J'aimerais également récupérer les valeurs des "NUMBER DENSITY", le problème étant que le script énoncé plus haut ne travail que ligne par ligne, or ici il me faudrait, par exemple, uniquement la valeur pour l'U235 (qui se trouve deux lignes après).

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13

CROSS SECTION OF MERGED/CONDENSED ISOTOPE 'O16 0001':   NUMBER DENSITY = 4.6455E-02     CROSS SECTION OF MERGED/CONDENSED ISOTOPE 'U235 0001':   NUMBER DENSITY = 1.1745E-03     CROSS SECTION OF MERGED/CONDENSED ISOTOPE 'U238 0001':   NUMBER DENSITY = 2.2051E-02

---

Merci beaucoup.

Je viens de trouver pour mon 2nd problème.

En fait, dès qu'il m'a trouvé:

Code:

---

CROSS SECTION OF MERGED/CONDENSED ISOTOPE 'U235 0001':

---

Je dois lui mettre, dans ma boucle if, de sauter deux fois de ligne et ensuite de me renvoyer la valeur cherchée.

Le problème néanmoins vient de la recherche du:

Code:

---

'U235

---

Code:

---

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

fid=fopen('pincell_uo2_burnup.result','rt'); tline = fgetl(fid); U235 = []; while ischar(tline) a = strread(tline,'%s'); if ~isempty(a) && strcmp(a{1},'CROSS') && strcmp(a{6},'U235') tline = fgetl(fid); tline = fgetl(fid); b = strread(tline,'%s'); U235 = [U235; str2double(b{4})]; end tline = fgetl(fid); end fclose(fid); U235

---

Cela ne fonctionne pas et si je modifie par:

Code:

---

strcmp(a{6},''U235')

---

Il y a un problème de syntaxe forcément à cause du << ' >>







Encore un petit edit de mon monologue. (Si ca peut servir à d'autres tant mieux)

Je viens de réussir pour le premier problème, pour récupérer la valeur du KINF, c'est du bidouillage mais ca fonctionne:

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 28 29 30

fid=fopen('pincell_uo2_burnup.result','rt'); tline = fgetl(fid); BURNUP = [0]; TIME = [0]; KINF = []; while ischar(tline)   if strfind(tline, 'FUEL BURNUP = '); a = strread(tline,'%s'); BURNUP = [BURNUP; str2double(a{4})]; end   if strfind(tline, '-> FINAL BURNUP AT TIME = '); a = strread(tline,'%s'); TIME = [TIME; str2double(a{7})]; end   if strfind(tline, 'TIMES FINAL KINF='); a = strread(tline,'%s','delimiter', '='); b = strread(a{3},'%s'); KINF= [KINF; str2double(b{1})]; end   tline = fgetl(fid); end fclose(fid); BURNUP TIME KINF plot(BURNUP,KINF)

---

Bonjour,

il manque encore un guillemet :

Code:

---

strcmp(a{6},'''U235')

---

Merci!

C'est vraiment aussi simple que cela ^^

Une version souvent plus efficace consiste à utiliser TEXTREAD et STRNCMP comme ceci :

Code:

---

1 2 3 4 5 6 7

X = textread('pincell_uo2_burnup.result','%s','delimiter','\n');   str = 'NUMBER DENSITY =';   idx = strncmp(X,str,numel(str));   strread(char(X(idx)).',[str '%f'])

---