计算音频 各频率 峰值能量

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

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

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

# 计算每个频率的峰值能量
freqs = np.fft.rfftfreq(frame_length, d=1./sr)
freq_energy = np.zeros_like(freqs)
for i in range(len(freqs)):
    idx = np.where(peak_freqs == i)[0]
    if len(idx) > 0:
        freq_energy[i] = np.sum(peak_energy[idx])

print(freq_energy)

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