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

Question

我一直在寻找一种方法，用Python编写代码来计算任意阶数的傅里叶级数曲线拟合。计算某个特定阶数的傅里叶级数曲线拟合，比如三阶的，还是比较简单的。但是如果要让这个阶数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)

你可以查看一下scipy和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