python中求无穷的傅立叶级数

Python中可以使用`scipy`库中的`fftpack`模块来求解傅立叶级数。下面是一个求解正弦函数的傅立叶级数的例子：

import numpy as np
from scipy.fftpack import fft, ifft
import matplotlib.pyplot as plt

# 定义正弦函数
def f(x):
    return np.sin(x)

# 定义采样点数和采样区间
N = 64
T = 2 * np.pi
# 生成采样点
x = np.linspace(0, T, N, endpoint=False)
# 计算采样点的函数值
y = f(x)

# 进行傅立叶变换
yf = fft(y)
# 计算频率
xf = np.linspace(0.0, 1.0/(2.0*T), N//2)

# 绘制原始信号和频域图像
fig, axs = plt.subplots(2, 1)
axs[0].plot(x, y)
axs[0].set_title('Original Signal')
axs[1].plot(xf, 2.0/N * np.abs(yf[:N//2]))
axs[1].set_title('Frequency Domain')

plt.show()

这个程序会生成一个正弦函数，并对其进行采样。然后，使用`fft`函数进行傅立叶变换，并计算频率。最后，使用`matplotlib`库绘制原始信号和频域图像。运行该程序会得到一个包含两个子图的窗口，第一个子图是原始信号的图像，第二个子图是频域图像的图像。