潮汐数据的FFT - 在MATLAB中识别Mf分量的问题

时间:2015-02-09 13:03:10

标签: matlab fft frequency frequency-analysis

我在MATLAB中运行一些潮汐数据的FFT。我正在使用FFT函数:

y = tide_data
Fs = 12*24; % Sampling frequency (number of samples per day)
L = numel(y(:,ii)); % Length of time series
NFFT = 2^nextpow2(L); % Next power of 2 from length of y
Y = fft(y(:,ii),NFFT)/L;
f = Fs/2*linspace(0,1,NFFT/2+1);
tt = (1./f);
figure (1)
semilogx(tt,2*abs(Y(1:NFFT/2+1)), 'k');

我成功地拾取了大部分主要潮汐物种,主要是主峰位于12.42小时(M2潮汐成分),但我似乎无法获得每两周一次的成分(Mf),甚至虽然在时域中显而易见。它应该具有大约327小时或14天的频率。请参阅下面的图片链接。有没有人对这是为什么有任何建议?

enter image description here

0 个答案:

没有答案