我想使用numpy和scipy实现xy数据集的傅立叶级数开发。
首先,我想将我的数据与前8个余弦拟合,并仅绘制第一个谐波。因此,我编写了以下两个函数定义:
# fourier series defintions
tau = 0.045
def fourier8(x, a1, a2, a3, a4, a5, a6, a7, a8):
return a1 * np.cos(1 * np.pi / tau * x) + \
a2 * np.cos(2 * np.pi / tau * x) + \
a3 * np.cos(3 * np.pi / tau * x) + \
a4 * np.cos(4 * np.pi / tau * x) + \
a5 * np.cos(5 * np.pi / tau * …Run Code Online (Sandbox Code Playgroud)