FFT on low sample count signal

Question

Jona Gladines 2019-1-8

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/438738-fft-on-low-sample-count-signal

评论： Star Strider 2019-1-14

I have sampled data of a slowly variating signal, about 1.85e-4 Hz . However since the variation is so slow I only have about 1.5 periods of data sampled at a rate about 28 times higher, so every 188 seconds I take a sample of the signal for a total of 9600 seconds. I want to do a fft analysis on this data to find the most important spectral component.

This gives a good approximation of the data I have:

Fs = 0.0053;           % Sampling frequency                    
T = 1/Fs;              % Sampling period       
L = 9800;              % Length of signal
t=0:T:L;               % Time Vector 
freq=1/5400;
X = (0.06*sin(2*pi*freq*t)+15.3)+(0.08*(rand(size(t))-0.5));
figure;
plot(t,X)
title('Signal')
xlabel('t (seconds)')
ylabel('X(t)')

I tried an FFT on this, but could nog get any good results. The highest spectral component is in the same order of magnitute as the original signal, but not exactly what I hoped for. I'm only interested in the frequency of the component.

Y = fft(X);
f=linspace(0, fs/2, N/2);
figure;
plot(f,abs(Y(1:N/2))) 
title('Single-Sided Amplitude Spectrum of X(t)')
xlabel('f (Hz)')
xlim([0 0.001]);
ylim([0 5]);

is there any way to improve without increasing the number of effective signal samples??

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

请先登录，再进行评论。

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

Answer 1

Star Strider 2019-1-8

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/438738-fft-on-low-sample-count-signal#answer_355486

Without seeing your signal, it is not possible to suggest a specif approach. One option may be doiing a nonlinear regression on your signal, using the approach in Curve fitting to a sinusoidal function (link).

4 个评论
显示 2更早的评论隐藏 2更早的评论

Star Strider 2019-1-11

在 MATLAB Online 中打开

Here is one approach:

The Code —

D = load('data.mat');
dates = D.dates;
p_mean = D.p_mean;
G = findgroups(dates);                                                      % Unique Time Values
p_meanu = splitapply(@mean, p_mean, G);                                     % Take Mean Of Values For Each ‘date’ Element
N = numel(p_meanu);                                                         % Number Of Samples
datesv = linspace(min(dates), max(dates), N);                               % Create Uniformly-Sampled Vector For Resampling Interpolation
Ts = datenum(mean(diff(datesv)));                                           % Sampling Interval (Days)
Fs = 1/Ts;                                                                  % Sampling Frequency (Samples/Day)
Fn = Fs/2;                                                                  % Nyquist Frequency (Samples/Day)
p_meanur = resample(p_meanu, datesv);                                       % Resample To Uniformly-Sampled Time Vector
NFFT = 2*nextpow2(N);                                                       % Length Of Fourier Transform
p_meanurm = p_meanur - mean(p_meanur);                                      % Subtreact Mean
FTp_meanur = fft(p_meanurm, 2^NFFT)/N;                                      % Fourier Transform
Fv = linspace(0, 1, fix(numel(FTp_meanur)/2)+1)*Fn;                         % Frequency Vector
Iv = 1:numel(Fv);                                                           % Correspoinding Index Vector
[pks, loc] = findpeaks(abs(FTp_meanur(Iv)), 'MinPeakHeight',0.0075);        % Fourier Transform Peak Am,plitudes & Locations
figure
plot(Fv, abs(FTp_meanur(Iv)))
hold on
plot(Fv(loc), pks, '^r', 'MarkerFaceColor','r')
grid
pklblc = compose('Amp = %6.4f\nPerd = %6.4f', [pks; 1./Fv(loc)]');
text(Fv(loc), pks, pklblc, 'HorizontalAlignment','left', 'VerticalAlignment','bottom')
ylim([0  0.04])

The Plot —

FFT on low sample count signal - 2019 01 11.png

The approach is straightforward. The Fourier Transfomr prefers uniformly-sampled signals, so I created a uniform vector of dates, and then used the resample function to resample it. (This is preferable to simple interpolation, because the resample funciton uses an anti-aliasing filter to remove unwanted frequencies from the the resampled signal.) The rest is relatively straightforward fft code, with a longer Fourier transform to provide improved frequency resolution. The plot presents the amplitude and the calculated period for each peak.

I did my best to comment-document it, so it should be relatively easy for you to experiment with it.

Jona Gladines 2019-1-14

Thanks,

This is indeed what I was looking for. I didn't think about resampling the data.

Jona

Star Strider 2019-1-14

As always, my pleasure.

Resampling to uniform sampling intervals is necessary in order to get reliable results from the fft.

请先登录，再进行评论。

FFT on low sample count signal

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

采纳的回答

4 个评论
显示 2更早的评论隐藏 2更早的评论

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

FFT on low sample count signal

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

采纳的回答

4 个评论 显示 2更早的评论隐藏 2更早的评论

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

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

4 个评论
显示 2更早的评论隐藏 2更早的评论