计算音频 4000hz 峰值能量

要计算音频在4000Hz处的峰值能量，可以先对音频信号进行分帧处理，然后对每一帧进行加窗和快速傅里叶变换（FFT），得到每一帧的频谱图。然后在频谱图上找到在4000Hz处的峰值，计算峰值能量值即可。

以下是一个简单的示例代码：

import numpy as np
import librosa

# 读取音频文件
y, sr = librosa.load('audio.wav', sr=None)

# 分帧处理
frame_length = int(0.025 * sr)  # 帧长
frame_step = int(0.01 * sr)  # 帧移
frames = librosa.util.frame(y, frame_length=frame_length, hop_length=frame_step).T

# 加窗和FFT
window = np.hanning(frame_length)
frames *= window
mag_frames = np.abs(np.fft.rfft(frames, axis=1))

# 计算峰值能量
peak_freqs = np.argmax(mag_frames, axis=1)
peak_vals = np.max(mag_frames, axis=1)
peak_energy = peak_vals ** 2

# 计算在4000Hz处的峰值能量
freqs = np.fft.rfftfreq(frame_length, d=1./sr)
idx_4000 = np.argmin(np.abs(freqs - 4000))
idx_peak_4000 = np.argmin(np.abs(peak_freqs - idx_4000))
energy_4000 = peak_energy[idx_peak_4000]

print(energy_4000)

其中，`librosa`是一个用于音频处理的库，可以用来读取音频文件和进行分帧处理。以上代码仅供参考，具体计算方法可以根据实际需要进行调整。