一种基于等效源法的室内声场扬声器重放方法

1.本发明属于空间声场重放技术领域，涉及一种基于等效源法的室内声场扬声器重放方法。


背景技术：

2.声场重放(sound field reproduction，sfr)旨在重现一个完整的声学环境或者合成一个期望的声学场景。空间声场重放技术是一种通过耳机或者多个扬声器为听者提供空间立体声体验的方法。使用耳机重放声场是比较便捷的方法，但是如果长期佩戴耳机会引起耳朵的不适，耳机重放在空间定位方面也存在不足。因此，对于多个听者，实现空间声的重放，需要借助扬声器重放的方式。为了避免扬声器阵列重放过程中复杂的计算，“一种多区域声场重放系统中扬声器布放优化的迭代方法”，其发明专利公开号为cn1l2565972a，为中国专利，提出了一种基于声压匹配的扬声器声场重放方法。但是此方法的控制点只限于传声器的测量点，重放误差大。因此，本设计采用等效源法对扬声器传递函数进行插值计算，降低重放误差。


技术实现要素：

3.要解决的技术问题
4.为了避免现有技术的不足之处，本发明提出一种基于等效源法的室内声场扬声器重放方法。
5.技术方案
6.一种基于等效源法的室内声场扬声器重放方法，其特征在于步骤如下：
7.步骤1、确定扬声器阵列的位置：将扬声器个数为l的扬声器阵列布置在所需的期望声场边界或在房间中贴墙摆放；
8.步骤2、设置目标声场区域：将个数为m的传声器置于目标声场区域边界，测量每个扬声器到每个传声器的房间脉冲响应，利用傅里叶变换，得到传声器个点处扬声器阵列的传递函数g
+
；在测量频率范围内，传声器位置上脉冲声源产生的峰值声压级至少高于相应频段内背景噪声45db；
9.所述目标声场区域为房间中听众期望的空间声场位置；
10.步骤3：设置目标区域内的重放声压p
des
，为已录制的声场环境声压或为合成的期望声场场景声压；
11.步骤4：将目标区域内按照个数为k网格化分布，设置等效声源的位置，个数记为n；第n个等效源r
n
到第m个传声器r
m
间的自由空间传递函数为：
[0012][0013]
其中为波数，f为频率，c为空气中的声速为340m/s；
[0014]
传声器r
m
点处的期望声压p
des
为n个等效源辐射声压之和：
[0015][0016]
其中w
n
为等效声源的权值系数，矩阵形式为：
[0017]
p
des
＝g
n，m
w
n
；
[0018]
步骤5：使用tikhonov正则化方法求解p
des
＝g
n，m
w
n
，得到等效源权值系数w
n
；
[0019]
步骤6：预测点k处的声压为：
[0020][0021]
在目标重放声场区域内预测点k处的由扬声器阵列产生的声压表示由等效声源表示：
[0022][0023]
矩阵形式为：
[0024]
g
k
＝g
n，k
w
n
[0025]
插值后的扬声器到目标声场重放区域内的声传递函数为：
[0026][0027]
步骤7：重放区域内的预先测量的或合成的期望声压为p
des
，利用tikhonov正则化求解声学逆运算：
[0028][0029]
计算扬声器的权重系数w，以此扬声器的权重系数，作为扬声器的驱动信号，从而实现扬声器声场重放。
[0030]
所述扬声器阵列为随意排布的扬声器阵列。
[0031]
所述扬声器阵列包括但不限于：线阵列，平面阵列，环形阵列或球形阵列。
[0032]
所述等效声源的个数n大于传声器测量点个数m。
[0033]
有益效果
[0034]
本发明提出的一种基于等效源法的室内声场扬声器重放方法，是一种测量目标重放区域边界处的房间脉冲响应，利用fft得到频域内声传递函数，利用等效源法插值重放区域内部扬声器到预测点的声传递函数，利用声场逆运算，求得扬声器阵列的驱动函数。声场逆运算通常为病态的，具有不适定性，使用正则化方法确保求解的稳定性。正则化方法选择tikhonov正则化。
[0035]
本发明使用等效源法插值扬声器阵列声传递函数进行声场重放，较目前使用的声压匹配方法，增加了扬声器阵列与传声器阵列排布的自由度，更为经济实用，可操作度高，计算速度快，重放效果精度高。在封闭空间声场重放中，可以取得良好的效果。
附图说明
[0036]
图1：某封闭空间中二维平面扬声器、传声器与等效源的分布示意图
[0037]
图2：声传递函数测量流程图
[0038]
图3：本发明声场重放方法的流程图
具体实施方式
[0039]
现结合实施例、附图对本发明作进一步描述：
[0040]
步骤1：确定扬声器阵列的位置。扬声器阵列的排布形式可以根据需要为线阵列，平面阵列，环形阵列，球形阵列等，也可以是随意排布的扬声器阵列。扬声器的个数为l。设定扬声器阵列在房间中的摆放位置。可以贴墙摆放，也可以根据需求，将扬声器阵列布置在所需的期望声场边界。
[0041]
步骤2：设置目标声场区域。目标声场区域通常为房间中听众期望的空间声场位置。在目标声场区域边界放置传声器，测量每个扬声器到每个传声器的房间脉冲响应。在测量频率范围内，传声器位置上脉冲声源产生的峰值声压级应至少高于相应频段内背景噪声45db。传声器的个数为m。利用傅里叶变换，得到传声器个点处扬声器阵列的传递函数g
+
。
[0042]
步骤3：设置目标区域内的重放声压p
des
，一般为已录制的声场环境声压，也可为合成的期望声场场景声压。
[0043]
步骤4：设置目标区域内的预测点位置，可取网格化分布，个数为k。设置等效声源的位置，个数记为n。第n个等效源r
n
到第m个传声器r
m
间的自由空间传递函数为其中为波数，f为频率，c为空气中的声速，一般为340m/s。传声器r
m
点处的期望声压p
des
为n个等效源辐射声压之和，则
[0044][0045]
其中w
n
为等效声源的权值系数。将上式写成矩阵形式为：
[0046]
p
des
＝g
n，m
w
n
ꢀꢀꢀ
(2)
[0047]
步骤5：通常情况下，等效声源的个数n大于传声器测量点个数m，则式(2)为欠定线性方程组。使用tikhonov正则化方法求解，求得等效源权值系数w
n
。
[0048]
步骤6：预测点k处的声压为：
[0049][0050]
在目标重放声场区域内预测点k处的由扬声器阵列产生的声压可以表示由等效声源表示，则：
[0051][0052]
矩阵形式为：
[0053]
g
k
＝g
n，k
w
n
ꢀꢀꢀ
(5)
[0054]
插值后的扬声器到目标声场重放区域内的声传递函数为：
[0055][0056]
步骤7：重放区域内的预先测量的或合成的期望声压为p
des
，利用tikhonov正则化求解声学逆运算：
[0057][0058]
求得扬声器的权重系数w。此扬声器的权重系数，可作为扬声器的驱动信号，从而实现扬声器声场重放。
[0059]
现结合实例、附图详细描述本发明的声场重放方法：假设在某封闭空间中，为简明与清晰描述，以二维平面声场重放为例(三维范例也可很容易的扩展)。
[0060]
封闭空间中二维平面扬声器、传声器与等效源的分布如附图1所示。房间长宽高为4m，4m，3m。在高为1.5m平面内，以中心点为坐标原点。环形扬声器阵列的半径为1m。环形等效源半径为1.2m，个数为45。方形声场重放区域长宽为1m，1m。
[0061]
步骤1：测量每个扬声器到所有传声器位置的房间脉冲响应。在测量频率范围内，传声器位置上脉冲声源产生的峰值声压级应至少高于相应频段内背景噪声45db。进行傅里叶变换计算每个扬声器分别到每个传声器的声传递函数，记为g
+
。
[0062]
步骤2：设定重放区域内预测点k的个数为100，呈网格离散均匀分布。
[0063]
步骤3：利用等效源法，计算扬声器到预测点处的声传递函数。
[0064]
g
k
＝g
n，k
w
n
ꢀꢀꢀ
(8)
[0065]
扬声器到预测点k处的声传递函数可以由等效源表示。g
n，k
为等效源到k点处的传递函数，w
n
为等效源的权值系数。
[0066]
步骤4：将步骤3所得的声传递函数g
k
与传声器测量的声传递函数g
+
进行插值，得到扬声器到重放区域的整体声传递函数
[0067]
步骤5：重放区域内的预先测量的或合成的期望声压为p
des
，利用tikhonov正则化求解声学逆运算：
[0068][0069]
求得扬声器的权重系数。从而实现扬声器声场重放。