使用python生成宽带语谱图的代码

生成宽带语谱图的常用方法是使用短时傅里叶变换（STFT），将音频信号分解成一系列短时窗口的频谱，并将这些频谱进行叠加得到语谱图。下面是使用Python实现生成宽带语谱图的代码示例：

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

# 读取音频文件
fs, audio = wavfile.read('audio.wav')

# 设置STFT参数
nperseg = int(fs * 0.025)  # 窗口长度
noverlap = int(fs * 0.01)  # 窗口重叠长度
nfft = 2 ** int(np.ceil(np.log2(nperseg)))  # FFT长度

# 计算STFT
f, t, Zxx = signal.stft(audio, fs, nperseg=nperseg, noverlap=noverlap, nfft=nfft)

# 将幅度谱转换为分贝
Zxx_db = 20 * np.log10(np.abs(Zxx))

# 绘制语谱图
plt.pcolormesh(t, f, Zxx_db, cmap='jet')
plt.ylabel('Frequency [Hz]')
plt.xlabel('Time [sec]')
plt.show()

其中，`audio.wav`是输入的音频文件名，`fs`为采样率，`audio`为音频数据。`nperseg`、`noverlap`和`nfft`分别为STFT的参数，`signal.stft()`函数用于计算STFT。最后将幅度谱转换为分贝，并用`plt.pcolormesh()`函数绘制语谱图。