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 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51
| clear all
close all
A=rdir('D:\Travail\DOSSIER_TEST\TEST_BOUCLE\test0*\*.txt');
for p=1:size(A,1)
fid=fopen(A(p).name,'r');
%% Traitement du fichier
N = 15; % N = Nombre de ligne à passer
for n=1:N;
fgetl(fid);
end;
%fréquence échantillonage
f_ech = 10000;
%tmax = temps d'enregistrement d'un événement
tmax = 0.5002;
%axe de temps
t = [0:1/f_ech:tmax-1/f_ech];
t=t';
% Lecture des données qui se trouvent sur deux colonnes
donnee=fscanf(fiche,'%g %g',[2 inf]);
% Mettre le vecteur de donnees dans le bon ordre
donnee=donnee';
voie1 = donnee(:,1);
voie2 = donnee(:,2);
voie1_fft = fft(voie1);
%fréquence pas
freq = [0:1/tmax:f_ech-1/tmax];
freq = freq';
%spectre d'énergie
spi = abs(voie1_fft).^2; %voie1_fft.*conj(voie1_fft)
figure
plot(freq, spi);
fclose(fid);
end |
Partager