Discussion :

[Djeu Kuru]

Bonjour,

J'aimerai réaliser un fft sur mon signal (en pièce jointe)
J'ai suivis un tutoriel mais je ne comprends pas très bien comment obternir un graph de densité spéctral de puissance.
Pourriez-vous m'aider ?

merci beaucoup pour votre aide,
DK

[Djeu Kuru]

Voici mon essai sur mon signal :

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
 
Fs = 1;                    % Sampling frequency
T = 1/Fs;                     % Sample time
L = length(I2(:,2));                     % Length of signal
t = (0:L-1)*T;                % Time vector
y = I2(:,2);
 
NFFT = 2^nextpow2(L); % Next power of 2 from length of y
Y = fft(y,NFFT)/L;
f = Fs/2*linspace(0,1,NFFT/2+1);
 
% Plot single-sided amplitude spectrum.
plot(f,2*abs(Y(1:NFFT/2+1))) 
axis([0 0.08 0 max(2*abs(Y(1:NFFT/2+1)))+10])
title('Single-Sided Amplitude Spectrum of y(t)')
xlabel('Frequency (Hz)')
ylabel('|Y(f)|')

Pensez-vous que c'est correct ?
Merci beaucoup