Discussion :

[Newenda]

Bonjour,

J'ai des fichiers (.raw) du type :

FFT Raw Data
============

Data Filename : C:\Documents and Settings\guillaume.chambon\Bureau\MAINSANT\TEST_2011_02_24_OSC_58\rawdata0011.raw
Creation Date : Thursday, February 24, 2011 10:48:22
Version : 1

General Heading : none
Start Of Test : Thursday, February 24, 2011 10:42:45

Param File XRef : {9BC5F336-4391-4B38-B186-260850B7E120}
Param File : none
Param Desc : none
Owner ID : none

Measuring System
- Title : PP 60 rugueux corrig
- Type : PARALLEL PLATES
- C1 : 1.768400e+04 Pa/Nm
- C2 : 3.000000e+01 1/rad
- C3 : 1.000000e+00 mm
- C4 : 5.100000e-06 Kgm^2

Air Bearing Inertia : 1.570000e-05 Kgm^2

Test Run : 1
Temperature Run : 1
Temperature Index : 1
Data Point : 11

Point Raw Torque Uncor Raw Disp Raw Disp Force(g)
0 -2.257303e-04 0.000000e+00 -1.020748e-05 -1.096230e+02
1 -2.460893e-04 2.041757e-05 1.021052e-05 -1.096230e+02

[.....]

2040 -6.068914e-05 9.764924e-06 4.416929e-07 -1.113550e+02
2041 -8.146503e-05 1.952985e-05 1.020705e-05 -1.113550e+02
2042 -1.022132e-04 7.101763e-06 -2.220601e-06 -1.113550e+02
2043 -1.229267e-04 1.509125e-05 5.769315e-06 -1.113550e+02
2044 -1.435987e-04 1.864213e-05 9.320630e-06 -1.113550e+02
2045 -1.642221e-04 3.550881e-06 -5.770182e-06 -1.113550e+02
2046 -1.847898e-04 2.219301e-05 1.287238e-05 -1.113550e+02
2047 -2.052947e-04 4.438602e-06 -4.881595e-06 -1.113550e+02

FFT Successful

Raw Torque Amp : 2.047464e+03
Raw Ang Disp Amp : 3.863534e-05 rad
Raw Phase Angle (rad) : 1.765852e-01 rad
Raw Phase Angle (deg) : 1.011759e+01 deg
G* : 5.176815e+04 Pa
G' : 5.096325e+04 Pa
G\\ : 9.093305e+03 Pa
Eta* : 4.773697e+03 Pa
Phase Angle : 1.011674e+01 deg
Stress : 1.999916e+01 Pa
Frequency : 1.725949e+00 Hz
Peak Strain : 3.863534e-04

***** EOF *****

Avec un header en tête et un header de fin. Au milieu les données (ici 2048 mais cela peut varier suivant les fichiers).

Question: j'aimerais pour tout fichier indépendamment du nombre de ligne de data, retrouver la valeur de frequency (ici = 1.725949e+00 Hz) en rouge.

Y'a t il un moyen simple, je suis paumé.

Merci d'avance !

[Jerome Briot]

Deux solutions :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
fid = fopen('rawdata0011.raw','rt');
 
str = fgetl(fid);
 
while ~strncmp(str,'Frequency :',11)
    str = fgetl(fid);    
end
 
fclose(fid);
 
F = sscanf(str,'Frequency : %f');

ou

Code :

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

1
2
3
4
5
X = textread('rawdata0011.raw','%s','delimiter','\n');
 
idx = strncmp(X,'Frequency :',11);
 
F = sscanf(X{idx},'Frequency : %f');

Ces deux codes supposent que la ligne commençant par Frequency est toujours présente. Dans le cas contraire, il faudra le prévoir.

[Newenda]

parfait merci !