Python中的n阶傅里叶级数曲线拟合

Question

我一直在寻找一种在Python中编码片段的方法，它可以计算任何n阶傅里叶级数曲线拟合。为了计算傅里叶级数曲线拟合的某个阶数，比如3阶很简单，但是在阶数n可变的情况下这样做仍然不可行。也许有人做过，但我的搜索还没找到。我想知道是否有人可以提供帮助。谢谢。

Answer 1

公式是

n-th cos_coeff = (2/T)*integral(-T/2,T/2, f(t)*cos(n*t*2*pi/2)dt)
n-th sin coeff = (2/T)*integral(-T/2,T/2, f(t)*sin(n*t*2*pi/2)dt)

检查scypi和scipy.integrate，了解有关集成的详细信息。

这应该是

cos_coeff(f, T, N) = (2/T)*quad(lambda t: f(t)*cos(N*t*2*math.pi/2),-T/2,T/2)

（虽未测试）

我不熟悉离散傅里叶变换，但您也可以从中计算所述系数。校验 http://docs.scipy.org/doc/scipy/reference/tutorial/fftpack.html