Problem using fft();

2 次查看（过去 30 天）
atharva aalok
atharva aalok 2021-10-21
评论： atharva aalok 2021-10-21
Ran in:
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?

请先登录，再回答此问题。

采纳的回答

KSSV
KSSV 2021-10-21
编辑：KSSV 2021-10-21
Remove the mean from the data and then use fft.
Say A is your data, (array). For fft use:
A = A-mean(A) ;
  1 个评论
atharva aalok
atharva aalok 2021-10-21
Yeah lol I just realised a second after posting thanks!

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Parametric Spectral Estimation 的更多信息

标签

产品


版本

R2021a

Community Treasure Hunt

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

Start Hunting!

Translated by