如何绘制python中谐波的总和？

import numpy as np
import matplotlib.pyplot as plt

# 定义谐波函数
def harmonic(x, n, amplitude, frequency):
    return amplitude * np.sin(n * frequency * x)

# 设置参数
x = np.linspace(0, 2*np.pi, 1000)  # x轴范围
n_harmonics = 5  # 谐波的数量
amplitude = 1  # 谐波的振幅
frequency = 1  # 谐波的频率

# 计算谐波的总和
y = np.zeros_like(x)
for n in range(1, n_harmonics+1):
    y += harmonic(x, n, amplitude, frequency)

# 绘制谐波的总和
plt.plot(x, y)
plt.xlabel('x')
plt.ylabel('y')
plt.title('Sum of Harmonics')
plt.show()