如何在 python 中的 FFT 上的箱之间插值

Question

我有四个频率峰值，这是我在应用FFT后得到的。现在我想知道这些频率峰值的精确值。有不同的插值方法。如何使用python对四个频率峰值使用这种插值方法

k=index of the max (possibly local) magnitude of an DFT
X[i]=bin "i" of an DFT 
|X[i]|=magnitude of DFT at bin "i“
k'=the interpolated bin location
Quadradic Method
y1=|X[k - 1]|
y2=|X[k]|
y3=|X[k + 1]|
d=(y3 - y1) / (2 * (2 * y2 - y1 - y3))
k'=k + d
Thanks in advance

Answer 1

如果你想要一个非常精确的估计，那么使用Sinc插值而不是抛物线，也许是连续近似来定位极值。