基于Python的音频三边定位算法：分贝响度转距离的实现疑问

问题背景 #

我正在用Python开发声源定位算法，通过3个已知位置的麦克风定位高声声源。思路是把每个麦克风作为圆心，用音频响度估算的距离作为半径画圆，交集区域即为声源位置。

目前遇到的问题是：用PyAudio录制音频计算出的响度是相对分贝值（例如远麦克风测得-30dB，近麦克风测得-15dB），无法直接用于三边定位。需要将这些相对分贝值转换为以米为单位的绝对距离。已知“距离每增加一倍，声级降低6分贝”，但不清楚具体实现方式。

当前计算响度的代码如下：

def _loudness(chunk):
    '''Calculate and return volume of input samples
    Input chunk should be a numpy array of samples for analysis, as
    returned by sound card.  Sound card should be in 16-bit mono mode.
    Return value is measured in dB, will be from 0dB (maximum
    loudness) down to -80dB (no sound).  Typical very loud sound will
    be -1dB, typical silence is -36dB.
    '''
    data = np.array(chunk, dtype=float) / 32768.0
    ms = math.sqrt(np.sum(data ** 2.0) / len(data))
    if ms < 10e-8: ms = 10e-8
    return 10.0 * math.log(ms, 10.0)

解决方案 #

1. 声级与距离的核心关系 #

声音在自由空间传播遵循平方反比定律，声级（dB）与距离的数学关系为：
[
L_2 = L_1 - 20 \log_{10}\left(\frac{r_2}{r_1}\right)
]
其中：

• (L_1) = 距离声源(r_1)处的声级（dB）
• (L_2) = 距离声源(r_2)处的声级（dB）

当(r_2 = 2r_1)时，代入公式可得到(L_2 = L_1 - 6)，这正是“距离加倍声级降6dB”的理论依据。

2. 校准参考点的必要性 #

你的代码输出的是相对声级（相对于麦克风能捕捉到的最大音量），因此必须先进行一次环境校准：

• 将声源放置在已知距离（例如1米）处，记录此时麦克风测得的声级(L_{ref})
• 这个参考值是将相对分贝转换为绝对距离的关键

3. 分贝转距离的代码实现 #

基于参考点，可推导出距离计算公式：
[
r = r_{ref} \times 10^{\frac{L_{ref} - L}{20}}
]
其中：

• (r_{ref}) = 校准用的参考距离（默认1米）
• (L_{ref}) = 参考距离处测得的声级（dB）
• (L) = 当前麦克风测得的声级（dB）

在现有代码基础上新增转换函数即可：

import math
import numpy as np

def _loudness(chunk):
    '''Calculate and return volume of input samples
    Input chunk should be a numpy array of samples for analysis, as
    returned by sound card.  Sound card should be in 16-bit mono mode.
    Return value is measured in dB, will be from 0dB (maximum
    loudness) down to -80dB (no sound).  Typical very loud sound will
    be -1dB, typical silence is -36dB.
    '''
    data = np.array(chunk, dtype=float) / 32768.0
    ms = math.sqrt(np.sum(data ** 2.0) / len(data))
    if ms < 10e-8: ms = 10e-8
    return 10.0 * math.log(ms, 10.0)

def db_to_distance(db, ref_db, ref_distance=1.0):
    '''Convert relative dB value to absolute distance in meters
    Args:
        db: 当前麦克风测得的声级（dB）
        ref_db: 参考距离处的声级（dB）
        ref_distance: 校准用的参考距离，默认1米
    Returns:
        声源到麦克风的距离（米）
    '''
    # 限制dB值范围，避免计算异常
    if db < -80:
        db = -80
    ratio = (ref_db - db) / 20.0
    return ref_distance * math.pow(10, ratio)

4. 实际使用的注意事项 #

• 校准环境一致性：环境中的反射、背景噪音会影响声级测量，校准步骤必须和实际定位在同一空间完成
• 麦克风同步：确保三个麦克风的录制完全同步，否则不同时间点的响度差异会导致定位误差
• 交集区域处理：由于测量误差，三个圆通常不会精确交于一点，建议用最小二乘法拟合最优位置，或取交集区域的几何中心作为声源位置

内容的提问来源于stack exchange，提问作者Jevon Mao