如何绘制两个不同频率的正弦波的和？

import numpy as np
import matplotlib.pyplot as plt

# 定义参数
f1 = 1   # 第一个正弦波频率
A1 = 1   # 第一个正弦波振幅
phi1 = 0 # 第一个正弦波相位

f2 = 2   # 第二个正弦波频率
A2 = 0.5 # 第二个正弦波振幅
phi2 = np.pi/2 # 第二个正弦波相位

fs = 100 # 采样率，每秒采样点数目

# 计算时间轴上的采样点
t = np.arange(0, 1, 1/fs)

# 计算两个正弦波的波形
y1 = A1 * np.sin(2*np.pi*f1*t + phi1)
y2 = A2 * np.sin(2*np.pi*f2*t + phi2)

# 计算两个正弦波的和
y = y1 + y2

# 绘制波形图
plt.plot(t, y1, color='blue', label='f1=1Hz')
plt.plot(t, y2, color='red', label='f2=2Hz')
plt.plot(t, y, color='green', label='f1+f2')
plt.legend()

# 添加标题和坐标轴标签
plt.title('Sum of Two Sinusoidal Waves')
plt.xlabel('Time')
plt.ylabel('Amplitude')

# 显示图像
plt.show()