帮我写一个用python求频谱图分贝的有效值

好的，以下是一个用Python计算频谱图分贝有效值的示例代码：

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

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

# 计算音频信号的频谱图
freqs, _, spec = plt.specgram(audio, Fs=fs, NFFT=1024, cmap=plt.cm.gist_heat)

# 计算频谱图每个点对应的分贝有效值
spec_db = 20 * np.log10(spec / np.max(spec))

# 显示频谱图
plt.imshow(spec_db, cmap=plt.cm.gist_heat, origin='lower', aspect='auto',
           extent=[freqs[0], freqs[-1], 0, audio.shape[0] / fs])
plt.xlabel('Frequency (Hz)')
plt.ylabel('Time (s)')
plt.colorbar()
plt.show()

这段代码中，我们首先使用`scipy.io.wavfile.read()`函数读取音频文件。然后，使用`plt.specgram()`函数计算音频信号的频谱图。该函数的参数中，`Fs`为采样率，`NFFT`为每个窗口的采样点数，`cmap`为色彩映射。

接下来，我们计算频谱图每个点对应的分贝有效值，即利用20*log10(spec/max(spec))的公式，将每个点表示为分贝数，其中`spec`为频谱图矩阵。

最后，使用`plt.imshow()`函数显示频谱图，并使用`plt.colorbar()`函数添加一个色条。