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

我正在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天。请参阅下面图片的链接。有没有人对为什么会这样有什么建议

---

什么是Mf？你能重新措辞这个问题吗？什么是你想要的，什么是你看不到的？你能发布代码来创建这些图片吗？MWew不可能告诉您做错了什么，因为您没有向我们展示您是如何做的，以及所需的行为。那么，您希望输出结果是什么呢？你的问题是这样表达的，你期待着一个领域专家来拯救这一天。不幸的是，并非所有人都有潮汐数据方面的经验，因此请重新表述您的问题，以显示频域中的预期输入和预期输出是什么，这将更适合像我这样的信号处理专家。很抱歉，希望我能通过更具体的描述来改进我所追求的？你希望看到什么样的振幅？如果它很小，你试过对数图吗？您的样本窗口是否与潮汐周期的自然频率一致？如果没有，那么您应该在执行FFT之前对数据加窗-您可能会丢失矩形窗口的旁瓣或噪声地板中的低振幅数据。