CN113395638B - 一种基于等效源法的室内声场扬声器重放方法 - Google Patents
一种基于等效源法的室内声场扬声器重放方法 Download PDFInfo
- Publication number
- CN113395638B CN113395638B CN202110573110.7A CN202110573110A CN113395638B CN 113395638 B CN113395638 B CN 113395638B CN 202110573110 A CN202110573110 A CN 202110573110A CN 113395638 B CN113395638 B CN 113395638B
- Authority
- CN
- China
- Prior art keywords
- loudspeaker
- sound field
- sound
- equivalent source
- microphone
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
- H04R5/00—Stereophonic arrangements
- H04R5/02—Spatial or constructional arrangements of loudspeakers
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04S—STEREOPHONIC SYSTEMS
- H04S7/00—Indicating arrangements; Control arrangements, e.g. balance control
- H04S7/30—Control circuits for electronic adaptation of the sound field
- H04S7/302—Electronic adaptation of stereophonic sound system to listener position or orientation
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
- H04R2205/00—Details of stereophonic arrangements covered by H04R5/00 but not provided for in any of its subgroups
- H04R2205/024—Positioning of loudspeaker enclosures for spatial sound reproduction
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Signal Processing (AREA)
- Circuit For Audible Band Transducer (AREA)
Abstract
本发明涉及一种基于等效源法的室内声场扬声器重放方法,是一种测量目标重放区域边界处的房间脉冲响应,利用FFT得到频域内声传递函数,利用等效源法插值重放区域内部扬声器到预测点的声传递函数,利用声场逆运算,求得扬声器阵列的驱动函数。声场逆运算通常为病态的,具有不适定性,使用正则化方法确保求解的稳定性。正则化方法选择Tikhonov正则化。本发明使用等效源法插值扬声器阵列声传递函数进行声场重放,较目前使用的声压匹配方法,增加了扬声器阵列与传声器阵列排布的自由度,更为经济实用,可操作度高,计算速度快,重放效果精度高。在封闭空间声场重放中,可以取得良好的效果。
Description
技术领域
本发明属于空间声场重放技术领域,涉及一种基于等效源法的室内声场扬声器重放方法。
背景技术
声场重放(Sound Field Reproduction,SFR)旨在重现一个完整的声学环境或者合成一个期望的声学场景。空间声场重放技术是一种通过耳机或者多个扬声器为听者提供空间立体声体验的方法。使用耳机重放声场是比较便捷的方法,但是如果长期佩戴耳机会引起耳朵的不适,耳机重放在空间定位方面也存在不足。因此,对于多个听者,实现空间声的重放,需要借助扬声器重放的方式。为了避免扬声器阵列重放过程中复杂的计算,“一种多区域声场重放系统中扬声器布放优化的迭代方法”,其发明专利公开号为CN1l2565972A,为中国专利,提出了一种基于声压匹配的扬声器声场重放方法。但是此方法的控制点只限于传声器的测量点,重放误差大。因此,本设计采用等效源法对扬声器传递函数进行插值计算,降低重放误差。
发明内容
要解决的技术问题
为了避免现有技术的不足之处,本发明提出一种基于等效源法的室内声场扬声器重放方法。
技术方案
一种基于等效源法的室内声场扬声器重放方法,其特征在于步骤如下:
步骤1、确定扬声器阵列的位置:将扬声器个数为L的扬声器阵列布置在所需的期望声场边界或在房间中贴墙摆放;
步骤2、设置目标声场区域:将个数为M的传声器置于目标声场区域边界,测量每个扬声器到每个传声器的房间脉冲响应,利用傅里叶变换,得到传声器个点处扬声器阵列的传递函数G+;在测量频率范围内,传声器位置上脉冲声源产生的峰值声压级至少高于相应频段内背景噪声45dB;
所述目标声场区域为房间中听众期望的空间声场位置;
步骤3:设置目标区域内的重放声压pdes,为已录制的声场环境声压或为合成的期望声场场景声压;
步骤4:将目标区域内按照个数为K网格化分布,设置等效声源的位置,个数记为N;第n个等效源rn到第m个传声器rm间的自由空间传递函数为:
传声器rm点处的期望声压pdes为N个等效源辐射声压之和:
其中wn为等效声源的权值系数,矩阵形式为:
pdes=Gn,mWn;
步骤5:使用Tikhonov正则化方法求解pdes=Gn,mWn,得到等效源权值系数Wn;
步骤6:预测点k处的声压为:
在目标重放声场区域内预测点k处的由扬声器阵列产生的声压表示由等效声源表示:
矩阵形式为:
Gk=Gn,kWn
插值后的扬声器到目标声场重放区域内的声传递函数为:
步骤7:重放区域内的预先测量的或合成的期望声压为pdes,利用Tikhonov正则化求解声学逆运算:
计算扬声器的权重系数W,以此扬声器的权重系数,作为扬声器的驱动信号,从而实现扬声器声场重放。
所述扬声器阵列为随意排布的扬声器阵列。
所述扬声器阵列包括但不限于:线阵列,平面阵列,环形阵列或球形阵列。
所述等效声源的个数N大于传声器测量点个数M。
有益效果
本发明提出的一种基于等效源法的室内声场扬声器重放方法,是一种测量目标重放区域边界处的房间脉冲响应,利用FFT得到频域内声传递函数,利用等效源法插值重放区域内部扬声器到预测点的声传递函数,利用声场逆运算,求得扬声器阵列的驱动函数。声场逆运算通常为病态的,具有不适定性,使用正则化方法确保求解的稳定性。正则化方法选择Tikhonov正则化。
本发明使用等效源法插值扬声器阵列声传递函数进行声场重放,较目前使用的声压匹配方法,增加了扬声器阵列与传声器阵列排布的自由度,更为经济实用,可操作度高,计算速度快,重放效果精度高。在封闭空间声场重放中,可以取得良好的效果。
附图说明
图1:某封闭空间中二维平面扬声器、传声器与等效源的分布示意图
图2:声传递函数测量流程图
图3:本发明声场重放方法的流程图
具体实施方式
现结合实施例、附图对本发明作进一步描述:
步骤1:确定扬声器阵列的位置。扬声器阵列的排布形式可以根据需要为线阵列,平面阵列,环形阵列,球形阵列等,也可以是随意排布的扬声器阵列。扬声器的个数为L。设定扬声器阵列在房间中的摆放位置。可以贴墙摆放,也可以根据需求,将扬声器阵列布置在所需的期望声场边界。
步骤2:设置目标声场区域。目标声场区域通常为房间中听众期望的空间声场位置。在目标声场区域边界放置传声器,测量每个扬声器到每个传声器的房间脉冲响应。在测量频率范围内,传声器位置上脉冲声源产生的峰值声压级应至少高于相应频段内背景噪声45dB。传声器的个数为M。利用傅里叶变换,得到传声器个点处扬声器阵列的传递函数G+。
步骤3:设置目标区域内的重放声压pdes,一般为已录制的声场环境声压,也可为合成的期望声场场景声压。
步骤4:设置目标区域内的预测点位置,可取网格化分布,个数为K。设置等效声源的位置,个数记为N。第n个等效源rn到第m个传声器rm间的自由空间传递函数为其中为波数,f为频率,c为空气中的声速,一般为340m/s。传声器rm点处的期望声压pdes为N个等效源辐射声压之和,则
其中wn为等效声源的权值系数。将上式写成矩阵形式为:
pdes=Gn,mWn (2)
步骤5:通常情况下,等效声源的个数N大于传声器测量点个数M,则式(2)为欠定线性方程组。使用Tikhonov正则化方法求解,求得等效源权值系数Wn。
步骤6:预测点k处的声压为:
在目标重放声场区域内预测点k处的由扬声器阵列产生的声压可以表示由等效声源表示,则:
矩阵形式为:
Gk=Gn,kWn (5)
插值后的扬声器到目标声场重放区域内的声传递函数为:
步骤7:重放区域内的预先测量的或合成的期望声压为pdes,利用Tikhonov正则化求解声学逆运算:
求得扬声器的权重系数W。此扬声器的权重系数,可作为扬声器的驱动信号,从而实现扬声器声场重放。
现结合实例、附图详细描述本发明的声场重放方法:假设在某封闭空间中,为简明与清晰描述,以二维平面声场重放为例(三维范例也可很容易的扩展)。
封闭空间中二维平面扬声器、传声器与等效源的分布如附图1所示。房间长宽高为4m,4m,3m。在高为1.5m平面内,以中心点为坐标原点。环形扬声器阵列的半径为1m。环形等效源半径为1.2m,个数为45。方形声场重放区域长宽为1m,1m。
步骤1:测量每个扬声器到所有传声器位置的房间脉冲响应。在测量频率范围内,传声器位置上脉冲声源产生的峰值声压级应至少高于相应频段内背景噪声45dB。进行傅里叶变换计算每个扬声器分别到每个传声器的声传递函数,记为G+。
步骤2:设定重放区域内预测点k的个数为100,呈网格离散均匀分布。
步骤3:利用等效源法,计算扬声器到预测点处的声传递函数。
Gk=Gn,kWn (8)
扬声器到预测点k处的声传递函数可以由等效源表示。Gn,k为等效源到k点处的传递函数,Wn为等效源的权值系数。
步骤5:重放区域内的预先测量的或合成的期望声压为pdes,利用Tikhonov正则化求解声学逆运算:
求得扬声器的权重系数。从而实现扬声器声场重放。
Claims (4)
1.一种基于等效源法的室内声场扬声器重放方法,其特征在于步骤如下:
步骤1、确定扬声器阵列的位置:将扬声器个数为L的扬声器阵列布置在所需的期望声场边界或在房间中贴墙摆放;
步骤2、设置目标声场区域:将个数为M的传声器置于目标声场区域边界,测量每个扬声器到每个传声器的房间脉冲响应,利用傅里叶变换,得到传声器个点处扬声器阵列的传递函数G+;在测量频率范围内,传声器位置上脉冲声源产生的峰值声压级至少高于相应频段内背景噪声45dB;
所述目标声场区域为房间中听众期望的空间声场位置;
步骤3:设置目标区域内的重放声压pdes,为已录制的声场环境声压或为合成的期望声场场景声压;
步骤4:将目标区域内按照个数为K网格化分布,设置等效声源的位置,个数记为N;第n个等效源rn到第m个传声器rm间的自由空间传递函数为:
传声器rm点处的期望声压pdes为N个等效源辐射声压之和:
其中wn为等效声源的权值系数,矩阵形式为:
pdes=Gn,mWn;
步骤5:使用Tikhonov正则化方法求解pdes=Gn,mWn,得到等效源权值系数Wn;
步骤6:预测点k处的声压为:
在目标重放声场区域内预测点k处的由扬声器阵列产生的声压表示由等效声源表示:
矩阵形式为:
Gk=Gn,kWn
插值后的扬声器到目标声场重放区域内的声传递函数为:
步骤7:重放区域内的预先测量的或合成的期望声压为pdes,利用Tikhonov正则化求解声学逆运算:
计算扬声器的权重系数W,以此扬声器的权重系数,作为扬声器的驱动信号,从而实现扬声器声场重放。
2.根据权利要求1所述基于等效源法的室内声场扬声器重放方法,其特征在于:所述扬声器阵列为随意排布的扬声器阵列。
3.根据权利要求1所述基于等效源法的室内声场扬声器重放方法,其特征在于:所述扬声器阵列包括但不限于:线阵列,平面阵列,环形阵列或球形阵列。
4.根据权利要求1所述基于等效源法的室内声场扬声器重放方法,其特征在于:所述等效声源的个数N大于传声器测量点个数M。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110573110.7A CN113395638B (zh) | 2021-05-25 | 2021-05-25 | 一种基于等效源法的室内声场扬声器重放方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110573110.7A CN113395638B (zh) | 2021-05-25 | 2021-05-25 | 一种基于等效源法的室内声场扬声器重放方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113395638A CN113395638A (zh) | 2021-09-14 |
CN113395638B true CN113395638B (zh) | 2022-07-26 |
Family
ID=77619045
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110573110.7A Active CN113395638B (zh) | 2021-05-25 | 2021-05-25 | 一种基于等效源法的室内声场扬声器重放方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113395638B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114390402B (zh) * | 2022-01-04 | 2024-04-26 | 杭州老板电器股份有限公司 | 一种吸油烟机的音频注入控制方法、装置及吸油烟机 |
CN115038010B (zh) * | 2022-04-26 | 2023-12-19 | 苏州清听声学科技有限公司 | 一种基于扬声器阵列的声场重建控制方法及系统 |
CN116684784B (zh) * | 2023-06-29 | 2024-03-12 | 中国科学院声学研究所 | 一种基于参量阵扬声器阵列的声重放方法及系统 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102121847A (zh) * | 2010-12-16 | 2011-07-13 | 合肥工业大学 | 一种瞬态声场重建方法 |
CN105181121A (zh) * | 2015-05-29 | 2015-12-23 | 合肥工业大学 | 采用加权迭代等效源法的高精度近场声全息算法 |
CN107566969A (zh) * | 2017-07-20 | 2018-01-09 | 西北工业大学 | 一种封闭环境内部低频声场重构方法 |
CN108632709A (zh) * | 2018-04-19 | 2018-10-09 | 北京工业大学 | 一种沉浸式宽带3d声场重放方法 |
CN109883532A (zh) * | 2019-03-14 | 2019-06-14 | 哈尔滨工程大学 | 一种声源识别与声场预报方法 |
CN110598238A (zh) * | 2019-07-16 | 2019-12-20 | 西北工业大学 | 一种舱室结构振动响应再现方法 |
CN110765673A (zh) * | 2019-12-30 | 2020-02-07 | 上海索辰信息科技有限公司 | 基于等效源法的噪声声场重建方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2592845A1 (en) * | 2011-11-11 | 2013-05-15 | Thomson Licensing | Method and Apparatus for processing signals of a spherical microphone array on a rigid sphere used for generating an Ambisonics representation of the sound field |
-
2021
- 2021-05-25 CN CN202110573110.7A patent/CN113395638B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102121847A (zh) * | 2010-12-16 | 2011-07-13 | 合肥工业大学 | 一种瞬态声场重建方法 |
CN105181121A (zh) * | 2015-05-29 | 2015-12-23 | 合肥工业大学 | 采用加权迭代等效源法的高精度近场声全息算法 |
CN107566969A (zh) * | 2017-07-20 | 2018-01-09 | 西北工业大学 | 一种封闭环境内部低频声场重构方法 |
CN108632709A (zh) * | 2018-04-19 | 2018-10-09 | 北京工业大学 | 一种沉浸式宽带3d声场重放方法 |
CN109883532A (zh) * | 2019-03-14 | 2019-06-14 | 哈尔滨工程大学 | 一种声源识别与声场预报方法 |
CN110598238A (zh) * | 2019-07-16 | 2019-12-20 | 西北工业大学 | 一种舱室结构振动响应再现方法 |
CN110765673A (zh) * | 2019-12-30 | 2020-02-07 | 上海索辰信息科技有限公司 | 基于等效源法的噪声声场重建方法 |
Non-Patent Citations (3)
Title |
---|
《基于等效源法插值的声场重放方法》;杜博凯,曾向阳,洪 汐;《声学学报》;20220228;第41卷(第1期);第103-107页 * |
基于声源叠加原理的大型号角扬声器声辐射特性分析;曾向阳,薛旭瑞,王海涛;《西北工业大学学报》;20131231;第31卷(第6期);第854-857页 * |
封闭空间声场重构的多层等效源法;刘延善,曾向阳;《声学学报》;20200531;第45卷(第3期);第367-376页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113395638A (zh) | 2021-09-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113395638B (zh) | 一种基于等效源法的室内声场扬声器重放方法 | |
JP6069368B2 (ja) | 組み合わせ又はハイブリッド制御方法を適用する方法 | |
Bernschütz | A spherical far field HRIR/HRTF compilation of the Neumann KU 100 | |
US20050223877A1 (en) | Sound system and method for creating a sound event based on a modeled sound field | |
JP5024792B2 (ja) | 全方位周波数指向性音響装置 | |
Sakamoto et al. | Sound-space recording and binaural presentation system based on a 252-channel microphone array | |
US20150180433A1 (en) | Sound processing apparatus, sound processing method, and program | |
JPWO2018008395A1 (ja) | 音場形成装置および方法、並びにプログラム | |
CN109891503A (zh) | 声学场景回放方法和装置 | |
JP6143571B2 (ja) | 音像定位装置 | |
JP2000092578A (ja) | スピーカ装置 | |
CN109196884A (zh) | 声音再现系统 | |
JP6613078B2 (ja) | 信号処理装置及びその制御方法 | |
JP2006507727A (ja) | オーディオ信号を再生するためのオーディオ再生システムおよび方法 | |
CN104581541A (zh) | 一种可定位多媒体视听装置及其控制方法 | |
CN108632709B (zh) | 一种沉浸式宽带3d声场重放方法 | |
JP2008118559A (ja) | 3次元音場再生装置 | |
CN115038010B (zh) | 一种基于扬声器阵列的声场重建控制方法及系统 | |
JPWO2019208285A1 (ja) | 音像再現装置、音像再現方法及び音像再現プログラム | |
JP2013009112A (ja) | 収音再生装置、プログラム及び収音再生方法 | |
JP5331553B2 (ja) | マルチチャンネル音響におけるインパルス応答測定システム、残響音生成装置、及び残響音生成システム | |
CN110099351B (zh) | 一种声场回放方法、装置和系统 | |
JP5941373B2 (ja) | スピーカアレイ駆動装置およびスピーカアレイ駆動方法 | |
JP2013106172A (ja) | 指向性スピーカ装置 | |
Gao et al. | Multizone sound reproduction with adaptive control of scattering effects |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |