CN113949983A - 一种听音区域声效恢复方法及装置 - Google Patents

一种听音区域声效恢复方法及装置 Download PDF

Info

Publication number
CN113949983A
CN113949983A CN202110575028.8A CN202110575028A CN113949983A CN 113949983 A CN113949983 A CN 113949983A CN 202110575028 A CN202110575028 A CN 202110575028A CN 113949983 A CN113949983 A CN 113949983A
Authority
CN
China
Prior art keywords
sound
zone
preset
generated
sound source
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.)
Granted
Application number
CN202110575028.8A
Other languages
English (en)
Other versions
CN113949983B (zh
Inventor
王松
张聪
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wuhan Polytechnic University
Original Assignee
Wuhan Polytechnic University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Wuhan Polytechnic University filed Critical Wuhan Polytechnic University
Priority to CN202110575028.8A priority Critical patent/CN113949983B/zh
Publication of CN113949983A publication Critical patent/CN113949983A/zh
Application granted granted Critical
Publication of CN113949983B publication Critical patent/CN113949983B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04SSTEREOPHONIC SYSTEMS 
    • H04S7/00Indicating arrangements; Control arrangements, e.g. balance control
    • H04S7/30Control circuits for electronic adaptation of the sound field
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04SSTEREOPHONIC SYSTEMS 
    • H04S7/00Indicating arrangements; Control arrangements, e.g. balance control
    • H04S7/30Control circuits for electronic adaptation of the sound field
    • H04S7/302Electronic adaptation of stereophonic sound system to listener position or orientation

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • Signal Processing (AREA)
  • Circuit For Audible Band Transducer (AREA)

Abstract

本发明公开了一种听音区域声效恢复方法及装置,包括:获取原始声源在预设坐标系中的声源位置坐标,和原始声源对应的频域信号;获取扬声器阵列中各扬声器在预设坐标系中的扬声器位置坐标;根据预设中心坐标确定扬声器阵列对应的声区,获取声区内各声音采样点的采样点坐标;根据各声音采样点的采样点坐标、声源位置坐标和频域信号确定原始声源在声区内产生的径向质点速度;根据各扬声器的扬声器位置坐标和各声音采样点的采样点坐标确定扬声器阵列在声区内产生的径向质点速度;根据原始声源和扬声器阵列在声区内产生的径向质点速度确定扬声器阵列的控制信号;根据控制信号控制所述扬声器阵列完成听音区域声效恢复,提高了听音区域声效恢复的准确度。

Description

一种听音区域声效恢复方法及装置
技术领域
本发明涉及声场重建技术领域,尤其涉及一种听音区域声效恢复方法及装置。
背景技术
目前,声音对比度控制方法(Acoustic Contrast Control method)通过最大化明声区范围内和暗声区范围内的声音对比度实现声场控制,该方法可以获得最大的明声区和暗声区对比度,但是在明声区范围内重建的声场和原始声源产生的声场的声压误差较大。已有声压匹配方法(Pressure Matching method)在控制区域内匹配原始声场和重建声场,可以获得最小的声压重建误差,但是明声区和暗声区范围内的声音对比度较小。声音对比度和声压重建误差都是声场控制比较重要的衡量指标,单独追求一个指标的最优将会忽略另一个指标的实现效果。因此要想实现个人听音区域声效的较好实现,需要综合考虑声音对比度和声压重建误差这两个指标。
上述内容仅用于辅助理解本发明的技术方案,并不代表承认上述内容是现有技术。
发明内容
本发明的主要目的在于提供了一种听音区域声效恢复方法及装置,旨在解决现有技术中听音区域声效恢复准确度低的技术问题。
为实现上述目的,本发明提供了一种听音区域声效恢复方法,所述方法包括以下步骤:
获取原始声源在预设坐标系中的声源位置坐标,以及所述原始声源对应的频域信号;
获取扬声器阵列中各扬声器在所述预设坐标系中的扬声器位置坐标;
根据预设中心坐标确定所述扬声器阵列对应的声区,并获取所述声区内各声音采样点的采样点坐标;
根据所述各声音采样点的采样点坐标、所述声源位置坐标和所述频域信号确定所述原始声源在所述声区内产生的径向质点速度;
根据所述各扬声器的扬声器位置坐标和所述各声音采样点的采样点坐标确定所述扬声器阵列在所述声区内产生的径向质点速度;
根据所述原始声源在所述声区内产生的径向质点速度和所述扬声器阵列在所述声区内产生的径向质点速度确定所述扬声器阵列的控制信号;
根据所述控制信号控制所述扬声器阵列完成听音区域声效恢复。
可选地,所述预设中心坐标包括第一预设中心坐标和第二预设中心坐标,所述声区包括明声区和暗声区,所述根据预设中心坐标确定所述扬声器阵列对应的声区,并获取所述声区内各声音采样点的采样点坐标,包括:
根据第一预设中心坐标和第一预设半径确定所述扬声器阵列对应的明声区;
根据第二预设中心坐标和第二预设半径确定所述扬声器阵列对应的暗声区;
在所述明声区内选取预设数量的明声区声音采样点,在所述暗声区内选取预设数量的暗声区声音采样点;
获取各明声区声音采样点的坐标和各暗声区声音采样点的坐标。
可选地,所述根据所述各声音采样点的采样点坐标、所述声源位置坐标和所述频域信号确定所述原始声源在所述声区内产生的径向质点速度,包括:
根据各声音采样点的采样点坐标、声源位置坐标、频域信号、声源频率和单位径向向内向量,通过第一预设公式确定所述原始声源在所述声区内产生的径向质点速度;
其中,所述第一预设公式为:
Figure BDA0003083070540000021
Figure BDA0003083070540000022
式中,
Figure BDA0003083070540000023
为采样点坐标,i为虚数单位,k=2πf/c,f表示声源频率,c表示声音传播速度,
Figure BDA0003083070540000024
表示声源位置坐标,
Figure BDA0003083070540000025
表示单位径向向内向量,s(ω)表示频域信号。
可选地,所述根据所述各扬声器的扬声器位置坐标和所述各声音采样点的采样点坐标确定所述扬声器阵列在所述声区内产生的径向质点速度,包括:
根据所述各扬声器的扬声器位置坐标、所述各采样点坐标、声源频率、单位径向向内向量,通过第二预设公式确定所述扬声器阵列在所述声区内产生的径向质点速度;
其中所述第二预设公式为:
u2=Eq
Figure BDA0003083070540000031
Figure BDA0003083070540000032
式中,
Figure BDA0003083070540000033
为采样点坐标,i为虚数单位,k=2πf/c,f表示声源频率,c表示声音传播速度,
Figure BDA0003083070540000034
为扬声器位置坐标,
Figure BDA0003083070540000035
表示单位径向向内向量,q为扬声器阵列的控制信号。
可选地,所述原始声源在所述声区内产生的径向质点速度包括原始声源在明声区内产生的径向质点速度和原始声源在暗声区内产生的径向质点速度,所述扬声器阵列在所述声区内产生的径向质点速度包括扬声器阵列在明声区内产生的径向质点速度和扬声器阵列在暗声区内产生的径向质点速度,所述根据所述原始声源在所述声区内产生的径向质点速度和所述扬声器阵列在所述声区内产生的径向质点速度确定所述扬声器阵列的控制信号,包括:
根据原始声源在明声区内产生的径向质点速度、扬声器阵列在明声区内产生的径向质点速度、扬声器阵列在暗声区内产生的径向质点速度和权重因子,通过预设代价函数确定所述扬声器阵列的控制信号;
其中,所述预设代价函数为:
F=σud1 Hud1+(1-σ)(ub1-ub0)H(ub1-ub0)
式中,σ为权重因子,ud1为扬声器阵列在暗声区内产生的径向质点速度,ub1为扬声器阵列在明声区内产生的径向质点速度,ub0为原始声源在明声区内产生的径向质点速度。
此外,为实现上述目的,本发明还提出一种听音区域声效恢复装置,所述装置包括:第一获取模块、第二获取模块、声区确定模块、第一径向质点速度确定模块、第二径向质点速度确定模块、控制信号确定模块和听音区域声效恢复模块:
所述第一获取模块,用于获取原始声源在预设坐标系中的声源位置坐标,以及所述原始声源对应的频域信号;
所述第二获取模块,用于获取扬声器阵列中各扬声器在所述预设坐标系中的扬声器位置坐标;
所述声区确定模块,用于根据预设中心坐标确定所述扬声器阵列对应的声区,并获取所述声区内各声音采样点的采样点坐标;
所述第一径向质点速度确定模块,用于根据所述各声音采样点的采样点坐标、所述声源位置坐标和所述频域信号确定所述原始声源在所述声区内产生的径向质点速度;
所述第二径向质点速度确定模块,用于根据所述各扬声器的扬声器位置坐标和所述各声音采样点的采样点坐标确定所述扬声器阵列在所述声区内产生的径向质点速度;
所述控制信号确定模块,用于根据所述原始声源在所述声区内产生的径向质点速度和所述扬声器阵列在所述声区内产生的径向质点速度确定所述扬声器阵列的控制信号;
所述听音区域声效恢复模块,用于根据所述控制信号控制所述扬声器阵列完成听音区域声效恢复。
可选地,所述预设中心坐标包括第一预设中心坐标和第二预设中心坐标,所述声区包括明声区和暗声区,所述声区确定模块还用于:
根据第一预设中心坐标和第一预设半径确定所述扬声器阵列对应的明声区;
根据第二预设中心坐标和第二预设半径确定所述扬声器阵列对应的暗声区;
在所述明声区内选取预设数量的明声区声音采样点,在所述暗声区内选取预设数量的暗声区声音采样点;
获取各明声区声音采样点的坐标和各暗声区声音采样点的坐标。
可选地,所述第一径向质点速度确定模块还用于:
根据各声音采样点的采样点坐标、声源位置坐标、频域信号、声源频率和单位径向向内向量,通过第一预设公式确定所述原始声源在所述声区内产生的径向质点速度;
其中,所述第一预设公式为:
Figure BDA0003083070540000051
Figure BDA0003083070540000052
式中,
Figure BDA0003083070540000053
为采样点坐标,i为虚数单位,k=2πf/c,f表示声源频率,c表示声音传播速度,
Figure BDA0003083070540000054
表示声源位置坐标,
Figure BDA0003083070540000055
表示单位径向向内向量,s(ω)表示频域信号。
可选地,所述第二径向质点速度确定模块还用于:
根据所述各扬声器的扬声器位置坐标、所述各采样点坐标、声源频率、单位径向向内向量,通过第二预设公式确定所述扬声器阵列在所述声区内产生的径向质点速度;
其中所述第二预设公式为:
u2=Eq
Figure BDA0003083070540000056
Figure BDA0003083070540000057
式中,
Figure BDA0003083070540000058
为采样点坐标,i为虚数单位,k=2πf/c,f表示声源频率,c表示声音传播速度,
Figure BDA0003083070540000059
为扬声器位置坐标,
Figure BDA00030830705400000510
表示单位径向向内向量,q为扬声器阵列的控制信号。
可选地,所述原始声源在所述声区内产生的径向质点速度包括原始声源在明声区内产生的径向质点速度和原始声源在暗声区内产生的径向质点速度,所述扬声器阵列在所述声区内产生的径向质点速度包括扬声器阵列在明声区内产生的径向质点速度和扬声器阵列在暗声区内产生的径向质点速度,所述控制信号确定模块还用于:
根据原始声源在明声区内产生的径向质点速度、扬声器阵列在明声区内产生的径向质点速度、扬声器阵列在暗声区内产生的径向质点速度和权重因子,通过预设代价函数确定所述扬声器阵列的控制信号;
其中,所述预设代价函数为:
F=σud1 Hud1+(1-σ)(ub1-ub0)H(ub1-ub0)
式中,σ为权重因子,ud1为扬声器阵列在暗声区内产生的径向质点速度,ub1为扬声器阵列在明声区内产生的径向质点速度,ub0为原始声源在明声区内产生的径向质点速度。
本发明通过获取原始声源在预设坐标系中的声源位置坐标,以及所述原始声源对应的频域信号;获取扬声器阵列中各扬声器在所述预设坐标系中的扬声器位置坐标;根据预设中心坐标确定所述扬声器阵列对应的声区,并获取所述声区内各声音采样点的采样点坐标;根据所述各声音采样点的采样点坐标、所述声源位置坐标和所述频域信号确定所述原始声源在所述声区内产生的径向质点速度;根据所述各扬声器的扬声器位置坐标和所述各声音采样点的采样点坐标确定所述扬声器阵列在所述声区内产生的径向质点速度;根据所述原始声源在所述声区内产生的径向质点速度和所述扬声器阵列在所述声区内产生的径向质点速度确定扬声器阵列的控制信号;根据控制信号控制扬声器阵列完成原始声源的听音区域声效恢复,提高了听音区域声效恢复的准确度。
附图说明
图1是本发明实施例方案涉及的硬件运行环境的听音区域声效恢复设备的结构示意图;
图2为本发明听音区域声效恢复方法第一实施例的流程示意图;
图3为本发明听音区域声效恢复方法一实施例平面直角坐标系示意图;
图4为本发明听音区域声效恢复方法第二实施例的流程示意图;
图5为本发明听音区域声效恢复装置第一实施例的结构框图。
本发明目的的实现、功能特点及优点将结合实施例,参照附图做进一步说明。
具体实施方式
应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
参照图1,图1为本发明实施例方案涉及的硬件运行环境的听音区域声效恢复设备结构示意图。
如图1所示,该听音区域声效恢复设备可以包括:处理器1001,例如中央处理器(Central Processing Unit,CPU),通信总线1002、用户接口1003,网络接口1004,存储器1005。其中,通信总线1002用于实现这些组件之间的连接通信。用户接口1003可以包括显示屏(Display)、输入单元比如键盘(Keyboard),可选用户接口1003还可以包括标准的有线接口、无线接口。网络接口1004可选的可以包括标准的有线接口、无线接口(如无线保真(WIreless-FIdelity,WI-FI)接口)。存储器1005可以是高速的随机存取存储器(RandomAccess Memory,RAM)存储器,也可以是稳定的非易失性存储器(Non-Volatile Memory,NVM),例如磁盘存储器。存储器1005可选的还可以是独立于前述处理器1001的存储装置。
本领域技术人员可以理解,图1中示出的结构并不构成对听音区域声效恢复设备的限定,可以包括比图示更多或更少的部件,或者组合某些部件,或者不同的部件布置。
如图1所示,作为一种存储介质的存储器1005中可以包括操作系统、数据存储模块、网络通信模块、用户接口模块以及听音区域声效恢复程序。
在图1所示的听音区域声效恢复设备中,网络接口1004主要用于与网络服务器进行数据通信;用户接口1003主要用于与用户进行数据交互;本发明听音区域声效恢复设备中的处理器1001、存储器1005可以设置在听音区域声效恢复设备中,所述听音区域声效恢复设备通过处理器1001调用存储器1005中存储的听音区域声效恢复程序,并执行本发明实施例提供的听音区域声效恢复方法。
本发明实施例提供了一种听音区域声效恢复方法,参照图2,图2为本发明听音区域声效恢复方法第一实施例的流程示意图。
本实施例中,所述听音区域声效恢复方法包括以下步骤:
步骤S10:获取原始声源在预设坐标系中的声源位置坐标,以及所述原始声源对应的频域信号。
需要说明的是,本实施例的执行主体可以是一种具有数据处理、网络通信以及程序运行功能的计算服务设备,例如平板电脑、个人电脑、手机等,或者是一种能够实现上述功能的电子设备、听音区域声效恢复设备等。以下以听音区域声效恢复设备为例,对本实施例及下述各实施例进行说明。
可以理解的是,预设坐标系可以是平面直角坐标系、空间直角坐标系,也可以是其他坐标系,可以根据实际应用场景确定具体的坐标系,建立平面直角坐标系,坐标系的原点为O,原点坐标为(0,0),若建立空间直角坐标系,坐标系原点为O,则原点坐标为(0,0,0),本实施例以预设坐标系为平面直角坐标系为例进行说明,并不构成本实施例的限制。
应该理解的是,原始声源是需要进行听音区域声效恢复的声源,原始声源可以是一个,也可以是多个,若有一个原始声源时,则获取该原始声源的位置坐标和该原始声源的频域信号;若有多个原始声源时,则获取各原始声源的位置坐标以及各原始声源对应的频域信号,本实施例以一个原始声源为例进行说明,并不构成对本实施例的限制。
步骤S20:获取扬声器阵列中各扬声器在所述预设坐标系中的扬声器位置坐标。
可以理解的是,扬声器阵列是进行听音区域声效恢复的扬声器组成的阵列,在扬声器阵列中有多个扬声器,可以根据实际使用场景确定扬声器的个数,扬声器阵列中的各扬声器可以被设置在一个圆形平面上,也可以设置在一个球面上,在扬声器阵列中的各扬声器设置在一个圆形平面上时,可以使用平面直角坐标系,在扬声器阵列中的各扬声器设置在一个球面上时,可以使用空间直角坐标系,在具体的使用中,可以存在多种使用情况,根据具体的使用场景选取合适的坐标系,本实施例对此不作限制。
步骤S30:根据预设中心坐标确定所述扬声器阵列对应的声区,并获取所述声区内各声音采样点的采样点坐标。
可以理解的是,在对原始声源进行听音区域声效恢复的时候,需要将目标声场聚焦于特定区域内,并且降低其他区域的声能量。
应该理解的是,预设中心坐标是提前设置的需要将目标声场聚焦或降低的特定区域的中心坐标,可以预设中心坐标为圆心,以设定的长度为半径确定声区,可将需要声场聚焦的区域设定为明声区,需要将声音能量降低的区域设定为暗声区。
可以理解的是,可以在声区内进行随机采样,获得若干个采样点,采样点的个数可根据实际使用场景对听音区域声效恢复精度的要求进行确定,在得到若干个采样点后,获取各采样点在预设坐标系中的采样点坐标,其中在明声区中采样点的坐标为明声区声音采样点坐标,在暗声区中采样点的坐标为暗声区声音采样点坐标。
步骤S40:根据所述各声音采样点的采样点坐标、所述声源位置坐标和所述频域信号确定所述原始声源在所述声区内产生的径向质点速度。
可以理解的是,根据明声区声音采样点坐标、声源位置坐标和原始声源对应的频域信号可以确定原始声源在明声区内产生的径向质点速度;根据暗声区声音采样点坐标、原始声源的位置坐标和原始声源对应的频域信号后,可以确定原始声源在暗声区内的径向质点速度。
步骤S50:根据所述各扬声器的扬声器位置坐标和所述各声音采样点的采样点坐标确定所述扬声器阵列在所述声区内产生的径向质点速度。
可以理解的是,在获取到扬声器位置坐标和明声区声音采样点坐标后,可以根据获取到的数据确定扬声器阵列在明声区内产生的径向质点速度;在获取到扬声器位置坐标和暗声区声音采样点坐标后,可以根据获取到的数据确定扬声器阵列在暗声区内产生的径向质点速度。
步骤S60:根据所述原始声源在所述声区内产生的径向质点速度和所述扬声器阵列在所述声区内产生的径向质点速度确定所述扬声器阵列的控制信号。
可以理解的是,在得到原始声源在明声区内产生的径向质点速度、原始声源在暗声区内产生的径向质点速度、扬声器阵列在明声区内产生的径向质点速度和扬声器阵列在暗声区内产生的径向质点速度后,使得扬声器阵列的控制信号对应的代价函数F最小,此时根据代价函数的最小值可求得对应的扬声器阵列的控制信号。
步骤S70:根据所述控制信号控制所述扬声器阵列完成听音区域声效恢复。
可以理解的是,在扬声器阵列的控制信号确定后,可根据所述控制信号完成原始声源的听音区域声效恢复。
进一步地,为了确定原始声源在声区内产生的径向质点速度,所述步骤S40包括:根据各声音采样点的采样点坐标、声源位置坐标、频域信号、声源频率和单位径向向内向量,通过第一预设公式确定所述原始声源在所述声区内产生的径向质点速度;其中,所述第一预设公式为:
Figure BDA0003083070540000101
Figure BDA0003083070540000102
式中,
Figure BDA0003083070540000103
为采样点坐标,i为虚数单位,k=2πf/c,f表示声源频率,c表示声音传播速度,
Figure BDA0003083070540000104
表示声源位置坐标,
Figure BDA0003083070540000105
表示单位径向向内向量,s(ω)表示频域信号。
可以理解的是,参照图3,图3为在平面直角坐标系下的示意图,也可以是在空间直角坐标系或其他坐标系,本实施例以平面直角坐标系为例进行说明,并不构成对本实施例的限制,假设
Figure BDA0003083070540000106
为明声区的中心点坐标,
Figure BDA0003083070540000107
为暗声区的中心点坐标,在明声区内随机选取t个声音采样点,可分别用
Figure BDA0003083070540000108
表示,对应的平面直角坐标为(h1x,h1y),(h2x,h2y),...(htx,hty),每一个采样点对应一个单位径向向量,可分别用
Figure BDA0003083070540000109
表示,将明声区内的采样点坐标,声源频率,声源位置坐标,单位径向向量和声源的频域信号代入如上公式,即可得到原始声源在明声区内产生的径向质点速度;在暗声区内随机选取l个声音采样点,可分别用
Figure BDA00030830705400001010
表示,对应的平面直角坐标为(k1x,k1y),(k2x,k2y),...(klx,kly),每一个采样点对应的单位径向向量可分别用
Figure BDA00030830705400001011
表示,将暗声区内的采样点坐标,声源频率,声源位置坐标,单位径向向量和声源频率代入如上公式,即可得到原始声源在暗声区内产生的径向质点速度。
进一步地,为了确定扬声器阵列在声区内产生的径向质点速度,所述步骤S50包括:根据所述各扬声器的扬声器位置坐标、所述各采样点坐标、声源频率、单位径向向内向量,通过第二预设公式确定所述扬声器阵列在所述声区内产生的径向质点速度;
其中所述第二预设公式为:
u2=Eq
Figure BDA0003083070540000111
Figure BDA0003083070540000112
式中,
Figure BDA0003083070540000113
为采样点坐标,i为虚数单位,k=2πf/c,f表示声源频率,c表示声音传播速度,
Figure BDA0003083070540000114
为扬声器位置坐标,
Figure BDA0003083070540000115
表示单位径向向内向量,q为扬声器阵列的控制信号。
可以理解的是,参照图3,假设扬声器整列中有n个扬声器,将n个扬声器的坐标,明声区的声音采样点坐标,声源频率,明声区的声音采样点坐标对应的单位径向向量代入上式,可得到扬声器阵列在明声区内产生的径向质点速度ub1;将n个扬声器的坐标,暗声区的声音采样点坐标,声源频率,暗声区的声音采样点坐标对应的单位径向向量代入上式,可得到扬声器阵列在暗声区内产生的径向质点速度ud1
进一步地,为了提高听音区域声效恢复的准确度,所述步骤S60包括:根据原始声源在明声区内产生的径向质点速度、扬声器阵列在明声区内产生的径向质点速度、扬声器阵列在暗声区内产生的径向质点速度和权重因子,通过预设代价函数确定所述扬声器阵列的控制信号;其中,所述预设代价函数为:
F=σud1 Hud1+(1-σ)(ub1-ub0)H(ub1-ub0)
式中,σ为权重因子,ud1为扬声器阵列在暗声区内产生的径向质点速度,ub1为扬声器阵列在明声区内产生的径向质点速度,ub0为原始声源在明声区内产生的径向质点速度。
可以理解的是,权重因子σ所属的区间为[0,1],可根据实际使用场景确定权重因子的取值,将权重因子,扬声器阵列在暗声区内产生的径向质点速度,扬声器阵列在明声区内产生的径向质点速度和原始声源在明声区内产生的径向质点速度代入预设代价函数,在代价函数F取最小值时,对应的q即为扬声器阵列的控制信号。
本实施例通过获取原始声源在预设坐标系中的声源位置坐标,以及所述原始声源对应的频域信号;获取扬声器阵列中各扬声器在所述预设坐标系中的扬声器位置坐标;根据预设中心坐标确定所述扬声器阵列对应的声区,并获取所述声区内各声音采样点的采样点坐标;根据所述各声音采样点的采样点坐标、所述声源位置坐标和所述频域信号确定所述原始声源在所述声区内产生的径向质点速度;根据所述各扬声器的扬声器位置坐标和所述各声音采样点的采样点坐标确定所述扬声器阵列在所述声区内产生的径向质点速度;根据所述原始声源在所述声区内产生的径向质点速度和所述扬声器阵列在所述声区内产生的径向质点速度确定所述扬声器阵列的控制信号;根据所述控制信号控制所述扬声器阵列完成听音区域声效恢复,由于本实施例是通过原始声源在声区内产生的径向质点速度和扬声器阵列在声区内产生的径向质点速度确定扬声器阵列的控制信号,根据控制信号控制扬声器阵列完成原始声源的听音区域声效恢复,提高了听音区域声效恢复的准确度。
参考图4,图4为本发明听音区域声效恢复方法第二实施例的流程示意图。
基于上述第一实施例,在本实施例中,所述步骤S30包括:
步骤S301:根据第一预设中心坐标和第一预设半径确定所述扬声器阵列对应的明声区。
可以理解的是,可以预先设定的第一预设中心坐标为圆心,以第一预设半径为圆的半径确定一个圆形区域,该圆形区域即为明声区。
应该理解的是,扬声器阵列对应的声区可以有多个明声区,若有多个明声区,可以根据预先设定的中心坐标和半径,确定其他明声区的区域范围,明声区个数的多少可根据实际使用场景设定,本实施例对此不作限制。
步骤S302:根据第二预设中心坐标和第二预设半径确定所述扬声器阵列对应的暗声区。
可以理解的是,可以预先设定的第二预设中心坐标为圆心,以第二预设半径为圆的半径确定一个圆形区域,该圆形区域即为暗声区。
应该理解的是,扬声器阵列对应的声区可以有多个暗声区,若有多个暗声区,可以根据预先设定的中心坐标和半径,确定其他暗声区的区域范围,暗声区个数的多少可根据实际使用场景设定,本实施例对此不作限制。
步骤S303:在所述明声区内选取预设数量的明声区声音采样点,在所述暗声区内选取预设数量的暗声区声音采样点。
可以理解的是,在明声区区域内,可以随机采集预设数量的明声区声音采样点,在暗声区内,可以随机采集预设数量的暗声区声音采样点。
步骤S304:获取各明声区声音采样点的坐标和各暗声区声音采样点的坐标。
可以理解的是,在明声区声音采样点和暗声区声音采样点确定后,可以获取各采样点在预设坐标系中的坐标。
本实施例根据第一预设中心坐标和第一预设半径确定所述扬声器阵列对应的明声区;根据第二预设中心坐标和第二预设半径确定所述扬声器阵列对应的暗声区;在所述明声区内选取预设数量的明声区声音采样点,在所述暗声区内选取预设数量的暗声区声音采样点;获取各明声区声音采样点的坐标和各暗声区声音采样点的坐标,能够满足不同的使用场景。
参照图5,图5为本发明听音区域声效恢复装置第一实施例的结构框图。
如图5所示,本发明实施例提出的听音区域声效恢复装置包括:第一获取模块10、第二获取模块20、声区确定模块30、第一径向质点速度确定模块40、第二径向质点速度确定模块50、控制信号确定模块60和听音区域声效恢复模块70;
所述第一获取模块10,用于获取原始声源在预设坐标系中的声源位置坐标,以及所述原始声源对应的频域信号;
所述第二获取模块20,用于获取扬声器阵列中各扬声器在所述预设坐标系中的扬声器位置坐标;
所述声区确定模块30,用于根据预设中心坐标确定所述扬声器阵列对应的声区,并获取所述声区内各声音采样点的采样点坐标;
所述第一径向质点速度确定模块40,用于根据所述各声音采样点的采样点坐标、所述声源位置坐标和所述频域信号确定所述原始声源在所述声区内产生的径向质点速度;
所述第二径向质点速度确定模块50,用于根据所述各扬声器的扬声器位置坐标和所述各声音采样点的采样点坐标确定所述扬声器阵列在所述声区内产生的径向质点速度;
所述控制信号确定模块60,用于根据所述原始声源在所述声区内产生的径向质点速度和所述扬声器阵列在所述声区内产生的径向质点速度确定所述扬声器阵列的控制信号;
所述听音区域声效恢复模块70,用于根据所述控制信号控制所述扬声器阵列完成听音区域声效恢复。
基于本发明上述听音区域声效恢复装置第一实施例,提出本发明听音区域声效恢复装置的第二实施例。
在本实施例中,所述声区确定模块30还用于根据第一预设中心坐标和第一预设半径确定所述扬声器阵列对应的明声区;根据第二预设中心坐标和第二预设半径确定所述扬声器阵列对应的暗声区;在所述明声区内选取预设数量的明声区声音采样点,在所述暗声区内选取预设数量的暗声区声音采样点;获取各明声区声音采样点的坐标和各暗声区声音采样点的坐标。
所述第一径向质点速度确定模块40还用于根据各声音采样点的采样点坐标、声源位置坐标、频域信号、声源频率和单位径向向内向量,通过第一预设公式确定所述原始声源在所述声区内产生的径向质点速度;其中,所述第一预设公式为:
Figure BDA0003083070540000141
Figure BDA0003083070540000142
式中,
Figure BDA0003083070540000143
为采样点坐标,i为虚数单位,k=2πf/c,f表示声源频率,c表示声音传播速度,
Figure BDA0003083070540000144
表示声源位置坐标,
Figure BDA0003083070540000145
表示单位径向向内向量,s(ω)表示频域信号。
所述第二径向质点速度确定模块50还用于根据所述各扬声器的扬声器位置坐标、所述各采样点坐标、声源频率、单位径向向内向量,通过第二预设公式确定所扬声器阵列在所述声区内产生的径向质点速度;其中所述第二预设公式为:
u2=Eq
Figure BDA0003083070540000151
Figure BDA0003083070540000152
式中,
Figure BDA0003083070540000153
为采样点坐标,i为虚数单位,k=2πf/c,f表示声源频率,c表示声音传播速度,
Figure BDA0003083070540000154
为扬声器位置坐标,
Figure BDA0003083070540000155
表示单位径向向内向量,q为扬声器阵列的控制信号。
所述控制信号确定模块60还用于根据原始声源在明声区内产生的径向质点速度、扬声器阵列在明声区内产生的径向质点速度、扬声器阵列在暗声区内产生的径向质点速度和权重因子,通过预设代价函数确定所述扬声器阵列的控制信号;其中,所述预设代价函数为:
F=σud1 Hud1+(1-σ)(ub1-ub0)H(ub1-ub0)
式中,σ为权重因子,ud1为扬声器阵列在暗声区内产生的径向质点速度,ub1为扬声器阵列在明声区内产生的径向质点速度,ub0为原始声源在明声区内产生的径向质点速度。
本发明听音区域声效恢复装置的其他实施例或具体实现方式可参照上述各方法实施例,此处不再赘述。
需要说明的是,在本文中,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者系统不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者系统所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括该要素的过程、方法、物品或者系统中还存在另外的相同要素。
上述本发明实施例序号仅仅为了描述,不代表实施例的优劣。
通过以上的实施方式的描述,本领域的技术人员可以清楚地了解到上述实施例方法可借助软件加必需的通用硬件平台的方式来实现,当然也可以通过硬件,但很多情况下前者是更佳的实施方式。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质(如只读存储器/随机存取存储器、磁碟、光盘)中,包括若干指令用以使得一台终端设备(可以是手机,计算机,服务器,空调器,或者网络设备等)执行本发明各个实施例所述的方法。
以上仅为本发明的优选实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换,或直接或间接运用在其他相关的技术领域,均同理包括在本发明的专利保护范围内。

Claims (10)

1.一种听音区域声效恢复方法,其特征在于,所述方法包括:
获取原始声源在预设坐标系中的声源位置坐标,以及所述原始声源对应的频域信号;
获取扬声器阵列中各扬声器在所述预设坐标系中的扬声器位置坐标;
根据预设中心坐标确定所述扬声器阵列对应的声区,并获取所述声区内各声音采样点的采样点坐标;
根据所述各声音采样点的采样点坐标、所述声源位置坐标和所述频域信号确定所述原始声源在所述声区内产生的径向质点速度;
根据所述各扬声器的扬声器位置坐标和所述各声音采样点的采样点坐标确定所述扬声器阵列在所述声区内产生的径向质点速度;
根据所述原始声源在所述声区内产生的径向质点速度和所述扬声器阵列在所述声区内产生的径向质点速度确定所述扬声器阵列的控制信号;
根据所述控制信号控制所述扬声器阵列完成听音区域声效恢复。
2.如权利要求1所述的方法,其特征在于,所述预设中心坐标包括第一预设中心坐标和第二预设中心坐标,所述声区包括明声区和暗声区,所述根据预设中心坐标确定所述扬声器阵列对应的声区,并获取所述声区内各声音采样点的采样点坐标,包括:
根据第一预设中心坐标和第一预设半径确定所述扬声器阵列对应的明声区;
根据第二预设中心坐标和第二预设半径确定所述扬声器阵列对应的暗声区;
在所述明声区内选取预设数量的明声区声音采样点,在所述暗声区内选取预设数量的暗声区声音采样点;
获取各明声区声音采样点的坐标和各暗声区声音采样点的坐标。
3.如权利要求1所述的方法,其特征在于,所述根据所述各声音采样点的采样点坐标、所述声源位置坐标和所述频域信号确定所述原始声源在所述声区内产生的径向质点速度,包括:
根据各声音采样点的采样点坐标、声源位置坐标、频域信号、声源频率和单位径向向内向量,通过第一预设公式确定所述原始声源在所述声区内产生的径向质点速度;
其中,所述第一预设公式为:
Figure FDA0003083070530000021
Figure FDA0003083070530000022
式中,
Figure FDA0003083070530000023
为采样点坐标,i为虚数单位,k=2πf/c,f表示声源频率,c表示声音传播速度,
Figure FDA0003083070530000024
表示声源位置坐标,
Figure FDA0003083070530000025
表示单位径向向内向量,s(ω)表示频域信号。
4.如权利要求3所述的方法,其特征在于,所述根据所述各扬声器的扬声器位置坐标和所述各声音采样点的采样点坐标确定所述扬声器阵列在所述声区内产生的径向质点速度,包括:
根据所述各扬声器的扬声器位置坐标、所述各采样点坐标、声源频率、单位径向向内向量,通过第二预设公式确定所述扬声器阵列在所述声区内产生的径向质点速度;
其中所述第二预设公式为:
u2=Eq
Figure FDA0003083070530000026
Figure FDA0003083070530000027
式中,
Figure FDA0003083070530000028
为采样点坐标,i为虚数单位,k=2πf/c,f表示声源频率,c表示声音传播速度,
Figure FDA0003083070530000029
为扬声器位置坐标,
Figure FDA00030830705300000210
表示单位径向向内向量,q为扬声器阵列的控制信号。
5.如权利要求4所述的方法,其特征在于,所述原始声源在所述声区内产生的径向质点速度包括原始声源在明声区内产生的径向质点速度和原始声源在暗声区内产生的径向质点速度,所述扬声器阵列在所述声区内产生的径向质点速度包括扬声器阵列在明声区内产生的径向质点速度和扬声器阵列在暗声区内产生的径向质点速度,所述根据所述原始声源在所述声区内产生的径向质点速度和所述扬声器阵列在所述声区内产生的径向质点速度确定所述扬声器阵列的控制信号,包括:
根据原始声源在明声区内产生的径向质点速度、扬声器阵列在明声区内产生的径向质点速度、扬声器阵列在暗声区内产生的径向质点速度和权重因子,通过预设代价函数确定所述扬声器阵列的控制信号;
其中,所述预设代价函数为:
F=σud1 Hud1+(1-σ)(ub1-ub0)H(ub1-ub0)
式中,σ为权重因子,ud1为扬声器阵列在暗声区内产生径向质点速度,ub1为扬声器阵列在明声区内产生的径向质点速度,ub0为原始声源在明声区内产生的径向质点速度。
6.一种听音区域声效恢复装置,其特征在于,所述装置包括:第一获取模块、第二获取模块、声区确定模块、第一径向质点速度确定模块、第二径向质点速度确定模块、控制信号确定模块和听音区域声效恢复模块;
所述第一获取模块,用于获取原始声源在预设坐标系中的声源位置坐标,以及所述原始声源对应的频域信号;
所述第二获取模块,用于获取扬声器阵列中各扬声器在所述预设坐标系中的扬声器位置坐标;
所述声区确定模块,用于根据预设中心坐标确定所述扬声器阵列对应的声区,并获取所述声区内各声音采样点的采样点坐标;
所述第一径向质点速度确定模块,用于根据所述各声音采样点的采样点坐标、所述声源位置坐标和所述频域信号确定所述原始声源在所述声区内产生的径向质点速度;
所述第二径向质点速度确定模块,用于根据所述各扬声器的扬声器位置坐标和所述各声音采样点的采样点坐标确定所述扬声器阵列在所述声区内产生的径向质点速度;
所述控制信号确定模块,用于根据所述原始声源在所述声区内产生的径向质点速度和所述扬声器阵列在所述声区内产生的径向质点速度确定所述扬声器阵列的控制信号;
所述听音区域声效恢复模块,用于根据所述控制信号控制所述扬声器阵列完成听音区域声效恢复。
7.如权利要求6所述的装置,其特征在于,所述预设中心坐标包括第一预设中心坐标和第二预设中心坐标,所述声区包括明声区和暗声区,所述声区确定模块还用于:
根据第一预设中心坐标和第一预设半径确定所述扬声器阵列对应的明声区;
根据第二预设中心坐标和第二预设半径确定所述扬声器阵列对应的暗声区;
在所述明声区内选取预设数量的明声区声音采样点,在所述暗声区内选取预设数量的暗声区声音采样点;
获取各明声区声音采样点的坐标和各暗声区声音采样点的坐标。
8.如权利要求6所述的装置,其特征在于,所述第一径向质点速度确定模块还用于:
根据各声音采样点的采样点坐标、声源位置坐标、频域信号、声源频率和单位径向向内向量,通过第一预设公式确定所述原始声源在所述声区内产生的径向质点速度;
其中,所述第一预设公式为:
Figure FDA0003083070530000041
Figure FDA0003083070530000042
式中,
Figure FDA0003083070530000043
为采样点坐标,i为虚数单位,k=2πf/c,f表示声源频率,c表示声音传播速度,
Figure FDA0003083070530000051
表示声源位置坐标,
Figure FDA0003083070530000052
表示单位径向向内向量,s(ω)表示频域信号。
9.如权利要求8所述装置,其特征在于,所述第二径向质点速度确定模块还用于:
根据所述各扬声器的扬声器位置坐标、所述各采样点坐标、声源频率、单位径向向内向量,通过第二预设公式确定所述扬声器阵列在所述声区内产生的径向质点速度;
其中所述第二预设公式为:
u2=Eq
Figure FDA0003083070530000053
Figure FDA0003083070530000054
式中,
Figure FDA0003083070530000055
为采样点坐标,i为虚数单位,k=2πf/c,f表示声源频率,c表示声音传播速度,
Figure FDA0003083070530000056
为扬声器位置坐标,
Figure FDA0003083070530000057
表示单位径向向内向量,q为扬声器阵列的控制信号。
10.如权利要求9所述的装置,其特征在于,所述原始声源在所述声区内产生的径向质点速度包括原始声源在明声区内产生的径向质点速度和原始声源在暗声区内产生的径向质点速度,所述扬声器阵列在所述声区内产生的径向质点速度包括扬声器阵列在明声区内产生的径向质点速度和扬声器阵列在暗声区内产生的径向质点速度,所述控制信号确定模块还用于:
根据原始声源在明声区内产生的径向质点速度、扬声器阵列在明声区内产生的径向质点速度、扬声器阵列在暗声区内产生的径向质点速度和权重因子,通过预设代价函数确定所述扬声器阵列的控制信号;
其中,所述预设代价函数为:
F=σud1 Hud1+(1-σ)(ub1-ub0)H(ub1-ub0)
式中,σ为权重因子,ud1为扬声器阵列在暗声区内产生的径向质点速度,ub1为扬声器阵列在明声区内产生的径向质点速度,ub0为原始声源在明声区内产生的径向质点速度。
CN202110575028.8A 2021-05-25 2021-05-25 一种听音区域声效恢复方法及装置 Active CN113949983B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110575028.8A CN113949983B (zh) 2021-05-25 2021-05-25 一种听音区域声效恢复方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110575028.8A CN113949983B (zh) 2021-05-25 2021-05-25 一种听音区域声效恢复方法及装置

Publications (2)

Publication Number Publication Date
CN113949983A true CN113949983A (zh) 2022-01-18
CN113949983B CN113949983B (zh) 2023-09-22

Family

ID=79327369

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110575028.8A Active CN113949983B (zh) 2021-05-25 2021-05-25 一种听音区域声效恢复方法及装置

Country Status (1)

Country Link
CN (1) CN113949983B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20160057539A1 (en) * 2014-08-20 2016-02-25 National Tsing Hua University Method for recording and reconstructing three-dimensional sound field
CN106559730A (zh) * 2016-11-25 2017-04-05 武汉大学 一种多声道系统声场重建方法及装置
CN107251579A (zh) * 2015-04-08 2017-10-13 华为技术有限公司 驱动扬声器阵列的装置和方法
CN108966114A (zh) * 2018-07-13 2018-12-07 武汉轻工大学 声场重建方法、音频设备、存储介质及装置
CN111464932A (zh) * 2020-04-07 2020-07-28 武汉轻工大学 基于多个听音点的声场重建方法、装置、设备及存储介质

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20160057539A1 (en) * 2014-08-20 2016-02-25 National Tsing Hua University Method for recording and reconstructing three-dimensional sound field
CN107251579A (zh) * 2015-04-08 2017-10-13 华为技术有限公司 驱动扬声器阵列的装置和方法
CN106559730A (zh) * 2016-11-25 2017-04-05 武汉大学 一种多声道系统声场重建方法及装置
CN108966114A (zh) * 2018-07-13 2018-12-07 武汉轻工大学 声场重建方法、音频设备、存储介质及装置
CN111464932A (zh) * 2020-04-07 2020-07-28 武汉轻工大学 基于多个听音点的声场重建方法、装置、设备及存储介质

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
王松,张聪: "一种非中心点处声场重建方法", 计算机应用与软件, vol. 37, no. 12, pages 101 - 107 *

Also Published As

Publication number Publication date
CN113949983B (zh) 2023-09-22

Similar Documents

Publication Publication Date Title
CN108733342B (zh) 音量调节方法、移动终端及计算机可读存储介质
CN110809214B (zh) 音频播放方法、音频播放装置及终端设备
CN110414232B (zh) 恶意程序预警方法、装置、计算机设备及存储介质
US11222646B2 (en) Apparatus and method for generating audio signal with noise attenuated based on phase change rate
CN110827843A (zh) 音频处理方法、装置、存储介质及电子设备
US11189298B2 (en) Acoustic zooming
US11250870B2 (en) Electronic device for supporting audio enhancement and method for the same
CN112233689B (zh) 音频降噪方法、装置、设备及介质
US10812031B2 (en) Electronic device and method for adjusting gain of digital audio signal based on hearing recognition characteristics
US11238563B2 (en) Noise processing method and apparatus
CN109961802B (zh) 音质比较方法、装置、电子设备及存储介质
CN112233688B (zh) 音频降噪方法、装置、设备及介质
US20220212108A1 (en) Audio frequency signal processing method and apparatus, terminal and storage medium
CN113949983A (zh) 一种听音区域声效恢复方法及装置
CN111464932A (zh) 基于多个听音点的声场重建方法、装置、设备及存储介质
CN110827850A (zh) 音频分离方法、装置、设备及计算机可读存储介质
CN115835119A (zh) 声效播放方法、装置、设备及存储介质
JP2020522189A (ja) インコヒーレント冪等アンビソニックスレンダリング
CN111916105B (zh) 语音信号处理方法、装置、电子设备及存储介质
CN112633218B (zh) 人脸检测方法、装置、终端设备及计算机可读存储介质
CN114264365A (zh) 风噪检测方法、装置、终端设备及存储介质
CN113920979A (zh) 语音数据的获取方法、装置、设备及计算机可读存储介质
CN109218920B (zh) 一种信号处理方法、装置及终端
CN116935883B (zh) 声源定位方法、装置、存储介质及电子设备
NL2025740B1 (en) Increasing low frequency extension for microspeakers using a volume dependent linkwitz transform and multiband compressor

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