绘制脉冲响应wav.文件的能量衰变曲线

绘制脉冲响应wav文件的能量衰减曲线，需要进行以下步骤：

1. 读取脉冲响应wav文件，获取其采样率和波形数据。

2. 对波形数据进行快速傅里叶变换（FFT），得到频域幅度谱。

3. 计算频域幅度谱每个频率的能量值，即幅度的平方。

4. 将能量值按照频率从低到高排序，计算累积能量，并归一化到总能量。

5. 绘制累积能量曲线，横坐标为频率，纵坐标为归一化的累积能量。

下面是一个Python代码示例：

import numpy as np
import matplotlib.pyplot as plt
import scipy.io.wavfile as wavfile

# 读取wav文件
sample_rate, data = wavfile.read("pulse_response.wav")

# 进行快速傅里叶变换
fft_data = np.fft.fft(data)
freq = np.fft.fftfreq(len(data), 1.0 / sample_rate)
freq = freq[:len(data) // 2]  # 取频率为正的部分
fft_data = np.abs(fft_data[:len(data) // 2])  # 取幅度谱的值

# 计算能量值
energy = fft_data ** 2

# 按照频率从低到高排序，计算累积能量
sort_idx = np.argsort(freq)
cum_energy = np.cumsum(energy[sort_idx])
cum_energy /= cum_energy[-1]

# 绘制能量衰减曲线
plt.plot(freq[sort_idx], cum_energy)
plt.xlabel("Frequency (Hz)")
plt.ylabel("Cumulative Energy")
plt.title("Energy Decay Curve")
plt.show()

运行以上代码，即可绘制出脉冲响应wav文件的能量衰减曲线。注意，这里的能量衰减曲线是归一化的，因此纵坐标的最大值为1。