如何创建振幅和相位变化的合成余弦波？

import numpy as np
import matplotlib.pyplot as plt

def create_cosine_wave(amplitude, phase_shift, frequency, duration, sampling_rate):
    # 计算采样点数
    num_samples = int(duration * sampling_rate)
    
    # 计算时间间隔
    time_interval = 1.0 / sampling_rate
    
    # 创建时间序列
    time = np.arange(0, duration, time_interval)
    
    # 计算每个采样点的角度
    angles = 2 * np.pi * frequency * time + phase_shift
    
    # 计算每个采样点的振幅
    amplitudes = amplitude * np.cos(angles)
    
    return amplitudes

# 设置参数
amplitude = 1.0  # 振幅
phase_shift = np.pi/2  # 相位变化
frequency = 1.0  # 频率
duration = 1.0  # 持续时间（秒）
sampling_rate = 44100  # 采样率

# 创建合成余弦波
wave = create_cosine_wave(amplitude, phase_shift, frequency, duration, sampling_rate)

# 绘制波形图
plt.plot(wave)
plt.xlabel('Time')
plt.ylabel('Amplitude')
plt.title('Synthesized Cosine Wave')
plt.show()