我有四个频率峰值，这是在应用 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

提前致谢