使用传声器计算SPL

---

前言

如何将传声器测量的数值转换成与声压计相同的数值？首先传声器又俗称麦克风，功能是将声能转化成机械能，然后再将机械能转换成电能。然后声卡再将模拟信号装换成数值信号，balabala~~~
无论他如何转换都是线性的，所以我们只需要使用校准音源校准我们的计算的数值，然后再做一个加权。
怕有些童鞋不知道SPL计算公式，还是写一下公式吧。

• $$\text{SPL=20*}lo{g}_{10}(\text{prefp0})$$

---

一、使用到的工具

• 标准声源，一般是1000Hz 94dB和1000Hz 114dB
• 标准声压计
• 传声器和功放

  

二、测量标准声源。

1、验证标准声压计与标准声源。

其实这一步可以去掉，主要是我想验证一下两个标准器件。

• 94dB
  
• 114dB
  

2、使用传声器录制标准音源。

这里使用的软件是Audacity；注意将电脑里面的音量设置里面的mic录制音量，与增益适当减小，还有功放里的增益。使传声器录制114dB时录制的波形超出了记录范围。如我录制的就刚刚好。

三、频谱分析

使用Audacity自带频谱分析工具；选择一段音频片段，然后点击工具栏的 分析/频谱分析进行频谱分析。

• 94dB
  

  峰值1000Hz 处的SPL为-20.2dB

• 114dB
  

  峰值1000Hz 处的SPL为-0.6dB

  可以看到两个相差的绝对值为20dB左右。所以如果要得到测量的 94/114 只需要在计算的基础上加上一个补偿值。

四、python实现

import librosa as lb
from librosa import display
import soundfile as sf
import numpy as np
import matplotlib.pyplot as plt

path = "音频路径"

data, rate = sf.read(path)

#进行短时傅里叶变换
freqs, times, mags = lb.reassigned_spectrogram(y=data, sr=rate)

mags_db = lb.amplitude_to_db(mags, ref=1.0)

# 输出第11个片段的声压值， 61.5是计算的补偿值，这个值是用94/114减去你计算的值
print(np.round(np.max(mags_db[:, 10]), 1) + 61.5) 

fig, ax = plt.subplots()

img = display.specshow(mags_db, x_axis="s", y_axis="log", sr=rate, ax=ax)

fig.colorbar(img, ax=ax, format="%+2.f dB")

plt.show()

---