Find the frequency for signals by using FFT

Question

Mohanned Al Gharawi 2021-4-6

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/794432-find-the-frequency-for-signals-by-using-fft

评论： Star Strider 2021-4-7

signals.xlsx

Hello every body,

Thank in advance. I have several signals. They have been created by using a collections of images for one week. Anyhow, the siganls are available in the attached excel sheet. The interval time between each reading is 1-minute. My question is how we could the FFT to calculate the domenant frequecny (or the frequency) for each signal.

Thank you again.

3 个评论
显示 1更早的评论隐藏 1更早的评论

Mohanned Al Gharawi 2021-4-7

Thank you for your answer, but I got confused!!!

That's my code:

s=data;

t = 1:length(s);

Ts = 1/(1*60); % Sampling Frequency in Hz

Fs = 1/Ts ; % Sampling Frequency

Fn = Fs/2; % Nyquist Frequency

L = length(s);

fts = fft(s)/L; % Normalised Fourier Transform

Fv = linspace(0, 1, fix(L/2)+1)*Fn; % Frequency Vector

Iv = 1:length(Fv); % Index Vector

amp = abs(fts(Iv))*2; % Spectrum Amplitude

phs = unwrap(angle(fts(

Mohanned Al Gharawi 2021-4-7

I got zero value for the frequency!!!

请先登录，再进行评论。

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

Answer 1

Star Strider 2021-4-7

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/794432-find-the-frequency-for-signals-by-using-fft#answer_668872

在 MATLAB Online 中打开

One approach (I recognise my code from some previous Answer, so I’ll supply the rest of it):

T1 = readtable('signals.xlsx', 'VariableNamingRule','preserve');
t = T1.('time (minutes)');
L = numel(t);
s = T1{:,2:9};
Ts = mean(diff(t))*60;                                      % Sampling Interval in Seconds
Fs = 1/Ts;                                                  % Sampling Frequency in Hz
Fn = Fs/2;                                                  % Nyquist Frequency
s_mean = mean(s);
fts = fft(s - s_mean)/L;                                    % Normalised Fourier Transform
Fv = linspace(0, 1, fix(L/2)+1)*Fn;                         % Frequency Vector
Iv = 1:length(Fv);                                          % Index Vector
amp = abs(fts(Iv))*2;                                       % Spectrum Amplitude
phs = unwrap(angle(fts));                                   % Spectrum Phase
figure
subplot(2,1,1)
plot(Fv, abs(fts(Iv,:))*2)
xlim([0 2]*1E-4)
grid
subplot(2,1,2)
plot(Fv, phs(Iv,:))
xlim([0 2]*1E-4)
grid