Problem using fft();

3 visualizaciones (últimos 30 días)
atharva aalok
atharva aalok el 21 de Oct. de 2021
Comentada: atharva aalok el 21 de Oct. de 2021
I have a periodic signal with TimePeriod~2.3 that is f~.4347
But I am getting another peak at f = 0.
load('myData.mat');
figure(1);
plot(tSol, omegaSol);
% This is periodic with TimePeriod~2.3 or f~.4347
Fs = length(tSol) / (max(tSol)-min(tSol)); % Sampling frequency
T = 1/Fs; % Sampling period
L = Fs * (max(tSol) - min(tSol)); % Length of signal
t = (0:L-1)*T; % Time vector
Y = fft(omegaSol);
P2 = abs(Y/L);
P1 = P2(1:L/2+1);
P1(2:end-1) = 2*P1(2:end-1);
f = Fs*(0:(L/2))/L;
limit_f = f < 1;
f = f(limit_f);
P1 = P1(limit_f);
figure(2);
plot(f,P1)
title('Single-Sided Amplitude Spectrum of X(t)')
xlabel('f (Hz)')
ylabel('|P1(f)|')
Why Is there another peak at f ~ 0?

Respuesta aceptada

KSSV
KSSV el 21 de Oct. de 2021
Editada: KSSV el 21 de Oct. de 2021
Remove the mean from the data and then use fft.
Say A is your data, (array). For fft use:
A = A-mean(A) ;
  1 comentario
atharva aalok
atharva aalok el 21 de Oct. de 2021
Yeah lol I just realised a second after posting thanks!

Iniciar sesión para comentar.

Más respuestas (0)

Categorías

Más información sobre Parametric Spectral Estimation en Help Center y File Exchange.

Productos


Versión

R2021a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by