How to extract correct FFT?

Question

Samuel 2012-4-23

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/36300-how-to-extract-correct-fft

Hello all,

I thank you guys all in advance for your help!

I am currently trying to plot an FFT plot of the strain gauge data I collected from the oscilloscope, and currently have a routine as defined below:

clear all; close all; clc
t = load('T0011.CSV');
time=t(1:125000,1);
gage=t(1:125000,2); % load files and assign to variables
totaltime=1; % total time captured by scope
figure(1)
plot(time,gage);
xlabel('time(s)')
ylabel('signal(volts)')
title('{\bf MTS data}')
legend('gage');
%
%compute frequency
%
fs=length(gage)/(gage(length(time))+abs(time(1)))
m=length(t);
n=m;
% n=pow2(nextpow2(m));
y=fft(gage)/n;
f=(0:n-1)*(fs/n);
figure(2)
semilogx(f(1:n/2),abs(y(1:n/2)));
title('Frequency (Logarithmic x-axis, Linear y-axis)','fontsize',14)
xlabel('Frequency (Hz)')
ylabel('|y(f)|')
legend('gage1');

Alright, the code below gives the bottom two figures: Signal Data http://imageshack.us/photo/my-images/109/fftb.jpg/ FFT Data http://imageshack.us/photo/my-images/833/fft2.jpg/ If you take a look, the first mode natural frequency is around 200 Hz, which is clearly not correct- going through the peak analysis reveals that the dominant frequency is around 40Hz.

I have noticed that even though the function FFT has been developed and made available, going through the different programming causes different frequency plots to come up (with different values), and wanted to make sure that my approach is correct.

Any help and/or advice will be greatly appreciated!

Sam

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

请先登录，再进行评论。

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

Answer 1

Dr. Seis 2012-4-23

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/36300-how-to-extract-correct-fft#answer_45440

在 MATLAB Online 中打开

Just to eliminate a few things, try this:

dt = time(2)-time(1);
Fs = 1/dt;
M = length(gage);
N = pow2(nextpow2(M));
Nyquist = Fs/2;
df = Fs/N;
f = fftshift(-Nyquist:df:Nyquist-df);
Y = fft(gage)*dt; % I would use dt or 1/Fs for normalization

2 个评论
显示无隐藏无

Samuel 2012-4-23

Excellent. Thanks so much for the prompt and accurate solution!

One more question to follow up in plotting it. I went ahead with the following code and made the plot command as shown above. Specifically:

plot(f(1:N/2),abs(Y(1:N/2)));

However, this command gives a straight line at 0 and 0 elsewhere.

The same command in semilogx command gives the correct plot, however.

Should anything be changed to plot it linearly?

thanks again.

sam

Dr. Seis 2012-4-23

I can't imagine why that would be... I would think that simply replacing "semilogx" with "plot" wouldn't cause weird behavior. Can you upload another jpg onto imageshack?

请先登录，再进行评论。