在信号处理理论中，方波可近似表示为多个正弦波的叠加。事实上，任意一个方波 信号都可以使用傅里叶变换为多个正弦波表示。利用 numpy 和 matplotlib 在坐标系中绘制方 波的无穷级数表示。请尝试调节正弦波的个数、幅度以及周期，尽可能使方波边缘平滑。方 波无穷级数表达式

方波的无穷级数表达式为：

$$f(x) = \frac{4}{\pi} \sum_{n=1}^{\infty} \frac{\sin((2n-1)\omega x)}{2n-1}$$

其中，$\omega$ 为方波的角频率，$x$ 为时间变量。为了绘制方波的无穷级数表示，我们可以选取一定范围内的时间点，计算出对应的函数值，然后在坐标系中绘制出来。以下是 Python 代码的实现：

import numpy as np
import matplotlib.pyplot as plt

# 方波的无穷级数表示
def square_wave(x, n):
    result = np.zeros_like(x)
    for i in range(1, n+1):
        coef = 1 / (2*i - 1)
        term = coef * np.sin((2*i - 1) * x)
        result += term
    return 4 / np.pi * result

# 绘制方波的无穷级数表示
x = np.linspace(-np.pi, np.pi, 1000)
y1 = square_wave(x, 1)
y2 = square_wave(x, 3)
y3 = square_wave(x, 5)
y4 = square_wave(x, 10)

plt.plot(x, y1, label="1 term")
plt.plot(x, y2, label="3 terms")
plt.plot(x, y3, label="5 terms")
plt.plot(x, y4, label="10 terms")
plt.legend()
plt.show()

运行上述代码，即可在坐标系中绘制出方波的无穷级数表示。在代码中，我们定义了 `square_wave()` 函数来计算方波的函数值，然后根据所选取的正弦波个数绘制出相应的图像。你可以尝试调节正弦波的个数、幅度以及周期，以获得更加平滑的方波表示。