FFT for the given data to find the dominant frequency

Question

Naveen K S 2019-5-3

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/460152-fft-for-the-given-data-to-find-the-dominant-frequency

评论： Star Strider 2022-6-3

采纳的回答： Star Strider

forcecoefficients1.xlsx

在 MATLAB Online 中打开

hello everyone,

I have a data file containing a variable called "Lift coefficient" collected over a period of 1.25 seconds with time steps of 1e-5. I am performing the FFT but i'm not getting the desirable result. I am using the following code. Please find the data file also. Thankyou in advance

x=Cl;
plot(t,x)
t=VarName1;
y=fft(x);
f = (0:length(y)-1)*(10000)/length(y);
plot(f,abs(y))

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Star Strider 2019-5-3

2
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/460152-fft-for-the-given-data-to-find-the-dominant-frequency#answer_373524

编辑：Star Strider 2019-5-3

在 MATLAB Online 中打开

Try this:

[D,S] = xlsread('forcecoefficients1.xlsx');
t = D(:,1);
Cl = D(:,2);
Ts = mean(diff(t));
Fs = 1/Ts;
Fn = Fs/2;
L = numel(t);
Clm = Cl-mean(Cl);                                  % Subtract Mean (Mean = 0 Hz)
FCl = fft(Clm)/L;
Fv = linspace(0, 1, fix(L/2)+1)*Fn;
Iv = 1:numel(Fv);
[pks,locs] = findpeaks(abs(FCl(Iv))*2, 'MinPeakHeight',0.01);
figure
plot(Fv, abs(FCl(Iv))*2)
grid
xlim([0 1E3])
text(Fv(locs), pks, sprintf('\\leftarrow Dominant Frequency = %.2f', Fv(locs)), 'HorizontalAlignment','left')

The dominant frequency appears to be 226.4.

EDIT —

Added plot image:

FFT for the given data to find the dominant frequency- 2019 05 03.png

2 个评论
显示无隐藏无

Guanting SU 2019-10-29

This is absoluately helpful!

Star Strider 2019-10-29

@Guanting SU — Thank you!

请先登录，再进行评论。

Answer 2

Image Analyst 2019-5-3

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/460152-fft-for-the-given-data-to-find-the-dominant-frequency#answer_373525

编辑：Image Analyst 2019-5-3

在 MATLAB Online 中打开

Is this homework?

Here's a start.

clc;    % Clear the command window.
close all;  % Close all figures (except those of imtool.)
clear;  % Erase all existing variables. Or clearvars if you want.
workspace;  % Make sure the workspace panel is showing.
format long g;
format compact;
fontSize = 20;
data = readtable('forcecoefficients1.xlsx');
t = data.x_Time;
x = data.Cl;
% Get rid of garbage first point.
x(1) = x(2);
subplot(2, 1, 1);
plot(t, x, 'b-')
title('x vs. t', 'FontSize', fontSize);
xlabel('t', 'FontSize', fontSize);
ylabel('x', 'FontSize', fontSize);
grid on;
y=fftshift(fft(x - mean(x)));
f = (0:length(y)-1)*(10000)/length(y);
subplot(2, 1, 2);
plot(abs(y), 'b-', 'LineWidth', 2)
grid on;
title('Spectral Power', 'FontSize', fontSize);
xlabel('frequency', 'FontSize', fontSize);
ylabel('Power', 'FontSize', fontSize);
% Enlarge figure to full screen.
set(gcf, 'Units', 'Normalized', 'OuterPosition', [0, 0.04, 1, 0.96]);

It looks like there are two main spectral ranges with power in them:

See the two spikes close to the origin? You might want to set xlim() to zoom in and see them. Be sure to properly set up the units for the frequency axis. Now, what frequencies do you see the two spikes at? I don't want to do everything for you since I strongly suspect this is your homework and I don't want to get you in trouble for copying.