如何在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插值代替抛物线，也许用逐次近似来定位极值。