CN103931211B - 处理刚性球上的球面麦克风阵列的信号的方法及装置 - Google Patents
处理刚性球上的球面麦克风阵列的信号的方法及装置 Download PDFInfo
- Publication number
- CN103931211B CN103931211B CN201280055175.1A CN201280055175A CN103931211B CN 103931211 B CN103931211 B CN 103931211B CN 201280055175 A CN201280055175 A CN 201280055175A CN 103931211 B CN103931211 B CN 103931211B
- Authority
- CN
- China
- Prior art keywords
- noise
- signal
- transmission function
- ambisonics
- microphone array
- 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
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/027—Spatial or constructional arrangements of microphones, e.g. in dummy heads
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
- H04R1/00—Details of transducers, loudspeakers or microphones
- H04R1/20—Arrangements for obtaining desired frequency or directional characteristics
- H04R1/32—Arrangements for obtaining desired frequency or directional characteristics for obtaining desired directional characteristic only
- H04R1/326—Arrangements for obtaining desired frequency or directional characteristics for obtaining desired directional characteristic only for microphones
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
- H04R3/00—Circuits for transducers, loudspeakers or microphones
- H04R3/005—Circuits for transducers, loudspeakers or microphones for combining the signals of two or more microphones
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
- H04R1/00—Details of transducers, loudspeakers or microphones
- H04R1/20—Arrangements for obtaining desired frequency or directional characteristics
- H04R1/32—Arrangements for obtaining desired frequency or directional characteristics for obtaining desired directional characteristic only
- H04R1/40—Arrangements for obtaining desired frequency or directional characteristics for obtaining desired directional characteristic only by combining a number of identical transducers
- H04R1/406—Arrangements for obtaining desired frequency or directional characteristics for obtaining desired directional characteristic only by combining a number of identical transducers microphones
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
- H04R2201/00—Details of transducers, loudspeakers or microphones covered by H04R1/00 but not provided for in any of its subgroups
- H04R2201/40—Details of arrangements for obtaining desired directional characteristic by combining a number of identical transducers covered by H04R1/40 but not provided for in any of its subgroups
- H04R2201/401—2D or 3D arrays of transducers
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
- H04R29/00—Monitoring arrangements; Testing arrangements
- H04R29/004—Monitoring arrangements; Testing arrangements for microphones
- H04R29/005—Microphone arrays
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04S—STEREOPHONIC SYSTEMS
- H04S2400/00—Details of stereophonic systems covered by H04S but not provided for in its groups
- H04S2400/15—Aspects of sound capture and related signal processing for recording or reproduction
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Signal Processing (AREA)
- Health & Medical Sciences (AREA)
- Otolaryngology (AREA)
- General Health & Medical Sciences (AREA)
- Circuit For Audible Band Transducer (AREA)
- Obtaining Desirable Characteristics In Audible-Bandwidth Transducers (AREA)
- Stereophonic System (AREA)
- Soundproofing, Sound Blocking, And Sound Damping (AREA)
Abstract
球面麦克风阵列捕捉三维声场{P(Ωc,t))以产生高保真度立体声响复制表示{Amn(t)),其中,通过阵列的炭精盒对在球的表面上的压力分布进行采样。使用逆麦克风传递函数移除麦克风对所捕捉的声场的影响。麦克风阵列的传递函数的等化是大问题,因为传递函数的倒数对在传递函数中的较小的值造成高增益并且这些较小的值受到换能器噪声影响。本发明通过在频域中使用维纳滤波器处理(34)使该噪声最小化,其中,根据麦克风阵列的信噪比针对每个波数自动地控制(33)该处理。
Description
技术领域
本发明涉及处理用于产生声场的高保真度立体声响复制(Ambisonics)表示的刚性球上的球面麦克风阵列的信号的方法及装置,其中,对逆麦克风阵列响应施加校正滤波器。
背景技术
球面麦克风阵列提供了捕捉三维声场的能力。一种存储并处理声场的方式是Ambisonics表示。Ambisonics使用正交球面函数描述在原点(也被称为“最佳听音位置(sweet spot)”)周围的区域中的声场。描述的准确度取决于Ambisonics阶N,其中,有限数量的Ambisonics系数描述声场。球面阵列的最高Ambisonics阶受限于麦克风炭精盒的数量,其中,该数量必须等于或大于Ambisonics系数的数量O=(N+1)2。
Ambisonics表示的一个有利之处在于,声场的再现可以单独地适用于任何指定的扬声器布置。而且,该表示使得能够在后期制作时使用音束成形技术对不同的麦克风特征进行仿真。
B格式是Ambisonics的一个已知示例。B格式麦克风需要在四面体上有四个炭精盒以捕捉Ambisonics阶为1的声场。
阶大于1的Ambisonics被称为高阶Ambisonics(HOA),典型地,HOA麦克风是在刚性球上的球面麦克风阵列,例如mhAcoustics的Eigenmike。为了进行Ambisonics处理,阵列的炭精盒对在球的表面上的压力分布进行采样。然后,将采样压力转换为Ambisonics表示。这种Ambisonics表示描述了声场,但是包含了麦克风阵列的影响。使用将平面波的声场变换为在麦克风炭精盒测量的压力的逆麦克风阵列响应移除麦克风对所捕捉的声场的影响。对炭精盒的方向性和麦克风阵列对声场的干扰进行仿真。
麦克风阵列的传递函数的等化对于HOA记录是大问题。如果阵列响应的Ambisonics表示是已知的,则可以通过Ambisonics表示与逆阵列响应的相乘来移除影响。然而,使用传递函数的倒数可能在传递函数中对较小的值和零造成高增益。因此,应当考虑健壮的逆传递函数来设计麦克风阵列。例如,B格式麦克风使用心形(cardioid)炭精盒来克服在全向炭精盒的传递函数中的零。
发明内容
本发明涉及在刚性球上的球面麦克风阵列。刚性球的遮挡效应使具有相对于阵列直径的较小的波长的频率能够具有良好的方向性。另一方面,这些麦克风阵列的滤波器响应对于低频和高Ambisonics阶(亦即大于1)具有非常小的值。因此,捕捉的压力的Ambisonics表示具有较小的更高阶系数,其表示与阵列大小相比较长的波长的在炭精盒上的小压差。压差受换能器噪声影响,并且因此更高阶系数也受换能器噪声影响。因此,对于低频,逆滤波器响应主要是放大噪声而不是放大更高阶Ambisonics系数。
用于克服该问题的已知技术是使低频的高阶逐渐消失(或者高通滤波器)(亦即,在该处限制滤波器增益),其一方面减少了低频的空间分辨率,但另一方面移除了(高度失真的)HOA系数,从而破坏了完整的Ambisonics表示。在Sébastien Moreau、Daniel、Stéphanie Bertet的《3D Sound field Recording with Higher Order Ambisonics-Objective Measurements and Validation of 4th Order Spherical Microphone》(Audio Engineering Society会议论文,2006年5月20日至23日第120届会议,法国巴黎)的第4节中描述了一种尝试使用Tikhonov正则化滤波器解决该问题的对应的补偿滤波器设计。Tikhonov正则化滤波器使由Ambisonics阶的限制导致的均方误差最小化。然而,Tikhonov滤波器需要必须通过“试错”手动地调整以适应所记录的信号的特征的正则化参数,并且没有定义该参数的解析式。根据Boaz Rafaely的《Analysis and Design ofSpherical Microphone Arrays》(IEEE Transactions on Speech and AudioProcessing,第13卷,第1号,第135至143页,2005年)的关于球面麦克风阵列的分析,本发明示出如何根据麦克风信号的信号统计自动地得到正则化参数。
本发明要解决的问题是使在布置在刚性球上的球面麦克风阵列的信号的Ambisonics表示中的噪声(特别是低频噪声)最小化。该问题通过在权利要求1中公开的方法得以解决。在权利要求2中公开了一种利用该方法的装置。
本发明的处理用于依赖平均声场功率与麦克风炭精盒的噪声功率的信噪比计算正则化Tikhonov参数,亦即,根据记录的麦克风阵列信号的信噪比计算该最优化参数。对最优化或正则化的参数的计算包含以下步骤:
-将表示在所述麦克风阵列的表面上的压力的麦克风炭精盒信号P(Ωc,t)转换为球面谐波(或者相等的Ambisonics)表示
-针对每个波数k,计算麦克风炭精盒信号P(Ωc,t)的时变信噪比SNR(k)的估计,其中使用从麦克风阵列记录的平面波的平均源功率|P0(k)|2以及表示由在麦克风阵列中的模拟处理产生的空间不相关的噪声的对应的噪声功率|Pnoise(k)|2,亦即,包含通过分别计算参考信号和噪声信号来计算平均空间功率,其中,参考信号是能够通过所使用的麦克风阵列创建的声场的表示,并且噪声信号是由麦克风阵列的模拟处理产生的空间不相关的噪声;
-通过使用对每个阶n对离散的有限的波数k根据信噪比估计SNR(k)设计的时变维纳(Wiener)滤波器,将维纳滤波器的传递函数乘以麦克风阵列的逆传递函数以便得到适应传递函数Fn,array(k);
-使用线性滤波器处理将对球面谐波表示施加该适应传递函数
Fn,array(k),得到适应方向系数
滤波器设计需要对声场的平均功率的估计以便得到记录的SNR。根据在球面谐波表示中的对在阵列的炭精盒上的平均信号功率的仿真推出该估计。该估计包含在球面谐波表示中的对炭精盒信号的空间一致性的计算。已知根据平面波的连续的表示计算空间一致性,但是根据本发明对在刚性球上的球面阵列计算空间一致性,因为不能以连续的表示计算在刚性球上的平面波的声场。亦即,根据本发明,根据炭精盒信号估计SNR。
本发明包含以下优点:
-Ambisonics表示的阶最优地适合于每个频率子带的记录的SNR。这降低了在Ambbisonics表示的再现时的可听噪声。
-滤波器设计需要对SNR的估计。其可以通过使用查找表以较低的计算复杂度实现。这有助于以可管理的计算工作量进行时变自适应滤波器设计。
-通过降噪,为低频部分地恢复方向信息。
基本上,本发明的方法适于处理在刚性球上的球面麦克风阵列的麦克风炭精盒信号,所述方法包含以下步骤:
-将表示在所述麦克风阵列的表面上的压力的所述麦克风炭精盒信号P(Ωc,t)转换为球面谐波或者Ambisonics表示
-针对每个波数k,使用从所述麦克风阵列记录的平面波的平均源功率|P0(k)|2以及表示由在所述麦克风阵列中的模拟处理产生的空间不相关的噪声的对应的噪声功率|Pnoise(k)|2,计算所述麦克风炭精盒信号P(Ωc,t)的时变信噪比SNR(k)的估计;
-通过使用对每个阶n对离散的有限的波数k根据所述信噪比估计SNR(k)设计的时变维纳滤波器,将所述维纳滤波器的传递函数乘以所述麦克风阵列的逆传递函数以便得到适应传递函数Fn,array(k);
-使用线性滤波器处理对所述球面谐波表示施加所述适应传递函数Fn,array(k),得到适应方向系数
基本上,本发明的装置适于处理在刚性球上的球面麦克风阵列的麦克风炭精盒信号,所述装置包含:
-适用于将表示在所述麦克风阵列的表面上的压力的所述麦克风炭精盒信号P(Ωc,t)转换为球面谐波或者Ambisonics表示约部件;
-适用于针对每个波数k,使用从所述麦克风阵列记录的平面波的平均源功率|P0(k)|2以及表示由在所述麦克风阵列中的模拟处理产生的空间不相关的噪声的对应的噪声功率|Pnoise(k)|2,计算所述麦克风炭精盒信号P(Ωc,t)的时变信噪比SNR(k)的估计的部件;
-适用于通过使用对每个阶n对离散的有限的波数k根据所述信噪比估计SNR(k)设计的时变维纳滤波器,将所述维纳滤波器的传递函数乘以所述麦克风阵列的逆传递函数以便得到适应传递函数Fn,array(k)的部件;
-适用于使用线性滤波器处理对所述球面谐波表示施加所述适应传递函数Fn,array(k),从而得到适应方向系数的部件。
在各个从属权利要求中公开了本发明的有利的另外的实施例。
附图说明
参考附图说明本发明的示例性的实施例,其中:
图1示出在刚性球上具有32个炭精盒的麦克风阵列的得到的扬声器权数中的参考、混叠及噪声分量的功率;
图2示出对于SNR(k)=20dB的降噪滤波器;
图3示出基于框的自适应Ambisonics处理的框图;
图4示出在图2的最优化滤波器之后的权数分量的平均功率。
具体实施方式
在以下部分中说明球面麦克风阵列处理。
Ambisonics理论
通过假设辐射平面波的声场的扬声器定义Ambisonics解码,参见M.A Poletti的《Three-Dimensional Surround Sound Systems Based on Spherical Harmonics》(Journal Audio Engineering Society,第53卷,第11号,第1004至1025页,2005年):
L个扬声器的布置重构了存储在Ambisonics系数中的三维声场。分别对每个波数
执行该处理,其中,f是频率,csound是声音的速度。索引n从0到有限阶N,而对于每个索引n,索引m从-n到n。因此,系数的总数为O=(N+1)2。在球面坐标中通过方向向量Ωl=[Θl,Φl]T定义扬声器位置,并且[·]T表示向量的转置形式。
等式(1)定义了Ambisonics系数到扬声器权数w(Ωl,k)的转换。这些权数是扬声器的驱动函数。所有扬声器权数的叠加重构了声场。
解码系数描述一般的Ambisonics解码处理。这包含在Morag Agmon、BoazRafaely的《Beamforming for a Spherical-Aperture Microphone》(IEEEI,第227至230页,2008年)的第3节中示出的音束方向图的共轭复系数以及在上述的M.A Poletti的论文的第3.2节中给出的模式匹配解码矩阵的行。在Johann-Markus Batke、FlorianKeiler的《Using VBAP-Derived Panning Functions for 3D Ambisonics Decoding》(Proc.Of the 2nd International Symposimum on Ambisonics and SphericalAcoustics,2010年5月6日至7日,法国巴黎)的第4节中描述的不同的处理方式使用基于向量的幅度平移计算用于任意三维扬声器布置的解码矩阵。还通过系数描述这些矩阵的行元素。
如在Boaz Rafaely的《Plane-wave decomposition of the sound field on asphere by spherical convolution》(J.Acoustical Society of America,第116卷,第4号,第2149至2157页,2004年)的第3节中所述,可以总是将Ambisonics系数分解为平面波的叠加。因此,该分析可能限制为从方向Ωs入射的平面波的系数:
定义平面波的系数,以采用辐射平面波的声场的扬声器。通过波数k的P0(k)定义在原点处的压力。共轭复球面谐波表示平面波的方向系数。使用在上述的M.A.Poletti的论文中给出的球面谐波的定义。
该球面谐波是Ambisonics表示的正交基函数,并且满足
其中,
是三角形脉冲。
球面麦克风阵列对在球的表面上的压力进行采样,其中,采样点的数量必须等于或者大于N阶Ambisonics的Ambisonics系数的数量O=(N+1)2。而且,采样点必须均匀分布在球的表面上,其中,仅确切知道关于阶N=1的O个点的最优分布。对于更高阶,存在球的采样的优良的近似值,参见:mh acoustics主页http://www.mhacoustics.com(访问日:2007年2月1日)以及F.Zotter的《Sampling Strategies for Acoustic Holography/Holophonyon the Sphere》(Proceedings of the NAG-DAGA,2009年3月23日至26日,鹿特丹)。
对于最优的采样点Ωc,根据等式(4)的积分等于根据等式(6)的离散和:
其中,对于C≥(N+1)2,n′≤N并且n≤N,C是炭精盒的总数。
为了获得非最优采样点的稳定结果,可以用从L×O的球面谐波矩阵得到的伪逆矩阵的列替换共轭复球面谐波,其中,球面谐波的O个系数是的行元素,参见上述的Moreau/Daniel/Bertet的论文中的第3.2.2节:
以下,定义用表示的列元素,使得根据等式(6)的正交条件还满足
其中,对于C≥(N+1)2,n′≤N并且n≤N。
如果假设球面麦克风阵列具有在球的表面上接近均匀分布的炭精盒并且炭精盒的数量大于O,则
成为有效表达式。将(9)代入(8),得到将要在下面被考虑的正交条件
其中,对于C≥(N+1)2,n′≤N并且n≤N。
处理的仿真
在刚性(硬性、固定)球上的球面麦克风阵列的完整的HOA处理链包含:估计炭精盒的压力,计算HOA系数,以及对扬声器权数进行解码。其基于:对平面波从麦克风阵列重构的权数w(k)必须等于从在等式(3)中给出的平面波系数重构的参考权数wref(k)。
以下部分介绍将w(k)分解为参考权数wref(k)、空间混叠权数walias(k)和噪声权数wnoise(k)。混叠是由对有限的阶N的连续的声场进行的采样造成的,并且噪声对为每个炭精盒引入的空间不相关的信号部分进行仿真。对于指定的麦克风阵列,无法移除空间混叠。
炭精盒信号的仿真
在上述的M.A.Poletti的论文的第2.2节的等式(19)中定义了在刚性球的表面上的麦克风阵列的入射平面波的传递函数:
其中,是第一类Hankel函数,半径r等于球的半径R。根据使得压力分散在刚性球上的物理原理推出该传递函数,这表示辐射速率在刚性球的表面上消失。换言之,进来的和分散的声场的辐射导出量(radial derivation)的叠加为零,参见《FourierAcoustics》一书的第6.10.3节。
因此,对于从Ωs入射的平面波,在位置Ω处的在球的表面上的压力在Moreau/Daniel/Bertet的论文的第3.2.1节的等式(21)中给出:
添加各向同性的噪声信号Pnoise(Ωc,k)以对换能器噪声进行仿真,其中,“各向同性”表示炭精盒的噪声信号是空间不相关的,这不包含在时域中的相关性。压力可以被分为对麦克风阵列的最大阶N计算的压力Pref(Ωc,kR)和来自剩余的阶的压力,参见在上述的Rafaely的论文《Analysis and design……》中的第7节的等式(24)。来自剩余的阶的压力Palias(Ωc,kR)被称为空间混叠压力,因为麦克风阵列的阶不足以重构这些信号分量。因此,在炭精盒c记录的全部压力被定义为:
Ambisonics编码
通过在等式(14a)中给出的对等式(12)进行的求逆从在炭精盒的压力得到Ambisonics系数参见上述的Moreau/Daniel/Bertet的论文的第3.2.2节的等式(26)。使用等式(8)通过对球面谐波求逆,并且通过它的逆,使传递函数bn(kR)等化:
如等式(14b)和(14c)所示,使用等式(14a)和(13a)可以将Ambisonics系数分为参考系数混叠系数以及噪声系数
Ambisonics解码
最优化使用在原点处的得到的扬声器权数w(k)。假设所有的扬声器到原点距离均相同,从而所有扬声器权数的和得到w(k)。等式(15)根据等式(1)和(14b)提供w(k),其中,L是扬声器的数量:
等式(15b)示出w(k)还可以被分为三个权数wref(k)、walias(k)和wnoise(k)。为了简化,此处不考虑在上述的Rafaely的论文《Analysis and design……》的第7节的等式(24)中给出的定位误差。
在解码中,参考系数是阶n的综合产生的平面波将创建的权数。在下面的等式(16a)中,将来自等式(13b)的参考压力Pref(Ωc,kR)代入等式(15a),由此压力信号Palias(Ωc,kR)和Pnoise(Ωc,k)被忽略(亦即,被设为零):
可以使用等式(8)消去c、n′和m′的和,使得等式(16a)可以被简化为根据等式(3)的Ambisonics表示中的平面波的权数的和。因此,如果忽略混叠和噪声信号,则可以根据麦克风阵列记录完美地重构阶N的平面波的理论系数。
根据等式(15a)并仅使用来自等式(13b)的Pnoise(Ωc,k),得到的噪声信号wnoise(k)的权数通过下式给出:
将来自等式(13b)的Palias(Ωc,kR)项代入到等式(15a)并忽略其他压力信号,得到:
因为索引n′大于N,所以不能通过来自等式(8)的正交化条件简化该得到的混叠权数walias(k)。
对混叠权数的仿真需要以足够的准确度表示炭精盒信号的Ambisonics阶。在上述的Moreau/Daniel/Bertet的论文的第2.2.2节的等式(14)中给出了关于Ambisonics声场重构的截断误差的分析。对于
可以得到该声场的合理的准确度,其中,表示向上舍入到最近的整数。该准确度用于仿真的频率上限fmax。因此,Ambisonics阶
被用于对每个波数的混叠压力进行的仿真。这得到在频率上限的可接受的准确度,甚至对于低频,准确度也提高了。
对扬场器权数的分析
图1示出关于在刚性球上具有32个炭精盒的麦克风阵列(根据上述Agmon/Rafaely的论文的Eigenmike被用于仿真)的来自方向Ωs=[0,0]T的平面波的得到的扬声器权数中的扬声器权数分量分量a)wref(k)、b)wnoise(k)和c)walias(k)的功率。麦克风炭精盒在R=4.2cm的球的表面上均匀地分布,使得正交条件得以满足。该阵列支持的最大的Ambisonics阶N为4。根据Fliege、Ulrike Maier的《A Two-Stage Approach for ComputingCubature Formulae for the Sphere》(技术报告,1996年,Fachbereich Mathematik,Dortmund,德国),在上述的M.A Poletti的论文中描述的模式匹配处理被用于得到关于25个均匀分布的扬声器位置的解码系数http://www.mathematik.uni-dortmund.de/lsx/research/projects/fliege/nodes/nodes.html示出了节点数量。
参考功率wref(k)在整个频率范围内是不变的。得到的噪声权数wnoise(k)在低频上显示出高功率,并在更高的频率上减少。通过方差为20dB(亦即,低于平面波功率的20dB)的正态分布的无偏伪随机噪声对噪声信号或功率进行仿真。混叠噪声walias(k)在低频上可以被忽略,但是随着不断上升的频率而增加,并且在10kHz以上时超过参考功率。混叠功率曲线的斜率取决于平面波的方向。然而,对于所有的方向,平均趋势是一致的。
两个误差信号wnoise(k)和walias(k)使在不同的频率范围内的参考权数失真。另外,误差信号彼此独立。因此,提出不考虑混叠信号使噪声信号最小化。
对所有的进来的平面波的方向,使在参考权数和失真的参考权数之间的均方误差最小化。忽略混叠信号walias(k)中的权数,因为在由Ambisonics表示的阶进行了空间限带之后无法校正walias(k)。这等同于其中无法从被采样并经空间限带的时间信号中移除混叠的时域混叠。
最优化-降噪
降噪使由噪声信号引入的均方误差最小化。在频域中使用维纳滤波器处理计算每个阶n的补偿滤波器的频率响应。对于每个波数k,根据参考权数wref(k)和经滤波器并失真的权数wref(k)+wnoise(k)得到误差信号。如前所述,此处忽略混叠误差walias(k)。使用最优化传递函数F(k)对失真的权数进行滤波,其中,在频域中通过失真信号和传递函数F(k)的相乘实施该处理。通过使参考权数与经滤波并失真的权数之间的平方误差的期望值最小化推出零相位传递函数F(k):
E{|wref(k)-F(k)(wref(k)+wnoise(k))|2}= (21a)
=E{|wref(k)|2}-2F(k)E{|wref(k)|2}
+F(k)2(E{|wref(k)|2}+E{|wnoise(k)|2}) (21b)
然后,通过下式给出作为维纳滤波器而众所周知的该解决方案:
平方绝对权数的期望值E表示权数的平均信号功率。因此,wnoise(k)和wref(k)的功率的分式表示关于每个波数k的经重构的权数的信噪比的倒数。在以下部分解释wnoise(k)和wref(k)的功率的计算。
根据上述的Rafaely的论文《Analysis and design……》的附录部分的等式(34)从等式(16)得到参考权数wref(k)的功率:
等式(24c)示出功率等于对所有的扬声器加在一起的平方绝对HOA系数的和。假设|P0(k)|2是平均声场能量并且P0(k)对于所有的Ωs均不变。这表示wref(k)的功率可以被分为每个阶n的功率的和。如果这对于wnoise(k)的期望值也是真的,则可以对每个阶n分别根据等式(21)使误差信号最小化,以便得到全局的最小值。
在上述的Rafaely的论文《Analysis and design……》的第7节的等式(28)中给出wnoise(k)的功率的推导。因为噪声信号是空间不相关的,所以可以为每个炭精盒独立地计算期望值。根据等式(17)通过下式推出噪声权数的期望的功率:
为了达到从每个阶n的功率的和中分离噪声功率权数,制定一些限制。如果扬声器C的和可以被简化为等式(10),则可以得到该分离。
因此,炭精盒位置必须接近均等地分布在球的表面上,使得等式(9)中的条件得以满足。而且,对于所有炭精盒,噪声压力的功率必须是不变的。然后,噪声功率独立于Ωc并且可以从C的和中排除。这样,对于所有炭精盒,通过下式定义不变的噪声功率:
施加这些限制,等式(25b)简化为:
因为阵列应当均匀地对在球上的压力进行采样,所以对于球面麦克风阵列,一般均满足对炭精盒位置的限制。对于由模拟处理(例传感器噪声或放大)以及对每个麦克风信号的模数转换产生的噪声,可以一直假设不变的噪声功率。因此,所述限制对于一般的球面麦克风阵列是有效的。
根据等式(21b)的期望值是参考功率与噪声功率的线性叠加。每个权数的功率均可以被分为每个阶n的功率的和。因此,根据等式(21b)的期望值也可以被分为每个阶n的叠加。这表示可以根据每个阶n的最小值推出全局的最小值,使得可以为每个阶n定义一个最优化传递函数Fn(k):
通过结合等式(23)、(24)和(25)从传递函数F(k)得到传递函数Fn(k)。N+1个最优化传递函数被定义为:
传递函数Fn(k)取决于炭精盒的数量以及波数k的信噪比:
另一方面,传递函数独立于Ambisonics解码器,这表示其对于三维Ambisonics解码和定向音束成形是有效的。因此,还可以不考虑解码系数的和,根据Ambisionics系数的均方误差推出传递函数。因为功率|P0(k)|2随时间变化,可以根据记录的信号的当前SNR(k)设计出自适应的传递函数。在“最优的Ambisonics处理”一节中进一步地说明该传递函数设计。
传递函数Fn(k)与来自在上述的Moreau/Daniel/Bertet论文中的第4节的等式(32)的Tikhonov正则化传递函数的比较表明可以从等式(29c)推出正则化参数λ。Tikhonov正则化的对应的参数
关于指定的SNR(k),使Ambisonics记录的平均重构误差最小化。
在图2中,传递函数Fn(k)被分别示为Ambisonics阶0到4的函数“a”到“e”,其中,传递函数具有对于每个阶n,对更高阶增大截止频率这样的高通特征。20dB的不变的SNR(k)被用于传递函数设计。截止频率以在上述的Moreau/Daniel/Bertet的论文中的第4.1.2节中所述的正则化参数λ衰减。因此,对于低频,需要高SNR(k)以得到更高阶Ambisonics系数。
根据下式计算最优化的权数w′(k):
最优的Ambisonics处理
在对Ambisonics麦克风阵列处理的实际的实现中,根据下式得到最优化的Ambisonics系数
其中,包含炭精盒c的和以及每个阶n的自适应传递函数以及波数k。该和将在球的表面上所采样的压力分布转换为Ambisonics表示,对于宽带信号,可以在时域中实施。该处理步骤将时域压力信号P(Ωc,t)转换为第一Ambisonics表示
在第二处理步骤中,最优化传递函数
根据第一Ambisonics表示重构方向信息项。传递函数bn(kR)的倒数将转换为方向系数其中,假设所采样的声场是通过分散在球的表面上的平面波的叠加而创建的。系数表示在上述的Rafaely的论文《Plane-wavedecomposition……》的第3节的等式(14)中所述的声场的平面波分解,这种表示主要用于Ambisonics信号的传送。取决于SNR(k),最优化传递函数Fn(k)降低更高阶系数的影响,以便移除被噪声掩盖的HOA系数。
可以将对系数的处理看作是线性滤波处理,其中,用Fn,array(k)确定滤波器的传递函数。这可以在频域以及时域中实施。FFT可以被用于将系数变换到频域,用于连续乘以传递函数Fn,array(k)。乘积的逆FFT得到时域系数该传递函数处理也被称为使用重叠相加或重叠保留方法的快速卷积。可替代地,可以通过FIR滤波器逼近线性滤波器,其系数可以通过以下方式根据传递函数Fn,array(k)计算出:使用逆FFT将其变换到时域,实施环形移位,并且对得到的滤波器脉冲响应施加锥形窗口以平滑对应的传递函数。然后,在时域中通过传递函数Fn,array(k)的时域系数与n和m的每个组合的系数的卷积实施线性滤波处理。
在图3中图示了本发明的自适应的基于框的Ambisonics处理。在上部的信号路径中,在步骤或阶段31中使用等式(14a)将麦克风炭精盒信号的时域压力信号P(Ωc,t)转换为Ambisonics表示由此,由麦克风传递函数bn(kR)进行的划分不执行(从而,计算而不是),而是在步骤/阶段32中执行。然后,步骤/阶段32在时域或频域中实施所述的线性滤波操作以便得到系数第二处理路径被用于传递函数Fn,array(k)的自动的自适应滤波器设计。步骤/阶段33实施对所考虑的时间段(亦即,一组采样)的信噪比SNR(k)的估计。对于有限数量的离散波数k在频域中实施该估计。因此,必须例如使用FFT将所涉及的压力信号P(Ωc,t)变换到频域。SNR(k)值是由两个功率信号|Pnoise(k)|2和|P0(k)|2指定的。噪声信号的功率|Pnoise(k)|2对于指定的阵列是不变的,其表示由炭精盒产生的噪声。必须根据压力信号P(Ωc,t)估计出平面波的功率|P0(k)|2。在“SNR估计”一节中进一步地说明该估计。根据所估计的SNR(k),在步骤/阶段34中设计具有n≤N的传递函数Fn,array(k)。滤波器设计包含在等式(29c)中给出的维纳滤波器以及逆阵列响应或逆传递函数1/bn(kR)的设计。有利地,维纳滤波器限制了逆阵列响应的传递函数的高度放大。这得到传递函数Fn,array(k)的可管理的放大。然后,该滤波器实现适用于在步骤/阶段32的时域或频域中的对应的线性滤波器处理。
SNR估计
将根据所记录的炭精盒信号估计SNR(k)值:其取决于平面波的平均功率|P0(k)|2以及噪声功率|Pnoise(k)|2。
在没有任何声源使得可以假设|P0(k)|2=0的安静的环境中根据等式(26)得到噪声功率。对于可调节的麦克风放大器,应当为若干放大器增益测量噪声功率。然后,该噪声功率可以适用于对若干记录所使用的放大器增益。
根据在炭精盒测量的压力Pmic(Ωc,k)估计平均源功率|P0(k)|2。这通过比较来自等式(13)的在炭精盒的压力的期望值与由下式定义的在炭精盒的所测量的平均信号功率来实施:
必须从所测量的功率中减于噪声功率|Pnoise(k)|2以得到期望值Psig(k)。
还可以通过下式为来自等式(13)的在炭精盒的压力的Ambisonics表示估计期望值Psig(k):
在等式(36b)中,可以对绝对量的扩展施加来自等式(4)的正交条件以推出等式(36c)。由此,根据球面谐波的互相关性估计出平均信号功率。结合传递函数bn(kR),这表示在炭精盒位置处的压力场的相干性。
等式(35)和(36)的等化根据所记录的压力信号Pmic(Ωc,k)和所估计的噪声功率|Pnoise(k)|2得到对|P0(k)|2的估计,其在等式(37)中呈现:
等式(37)中的分母对于指定的麦克风阵列的每个波数k是不变的。因此,可以只针对要被存储在查找表中的或为每个波数k存储的Ambisonics阶Nmax计算一次。
最后,通过下式根据炭精盒信号P(Ωc,kR)得到SNR(k)值:
还可以根据线性麦克风阵列处理得知根据指定炭精盒信号的平均源功率的估计。炭精盒信号的互相关性被称为声场的空间相干性。对于线性阵列处理,根据平面波的连续表示确定该空间相干性。仅以Ambisonics表示的形式得知对在刚性球上的分散的声场的描述。因此,对SNR(k)的所呈现的估计基于确定在刚性球的表面上的空间相干性的新处理。
因此,对于模式匹配的Ambisonics解码器,在图4中示出从图2的最优化滤波器得到的平均功率分量w′(k)。噪声功率降低至-35dB,直至频率为1kHz时。在高于1kHz时,噪声功率线性地增加到-10dB。得到的噪声功率小于Pnoise(Ωc,k)=-20dB,直至频率约为8kHz时。在高于10kHz时,全部功率上升10dB,这是由混叠功率造成的。在高于10kHz时,对于半径等于R的球,麦克风阵列的HOA阶不能充分地描述在表面上的压力分布。因此,由所得到的Ambisonics系数造成的平均功率大于参考功率。
Claims (8)
1.一种处理在刚性球上的球面麦克风阵列的麦克风炭精盒信号(P(Ωc,t))的方法,所述方法包含以下步骤:
将表示在所述麦克风阵列的表面上的压力的所述麦克风炭精盒信号(P(Ωc,t))转换(31)为球面谐波或者高保真度立体声响复制表示
针对每个波数k,使用从所述麦克风阵列记录的平面波的平均源功率|P0(k)|2以及表示由在所述麦克风阵列中的模拟处理产生的空间不相关的噪声的对应的噪声功率|Pnoise(k)|2,计算(33)所述麦克风炭精盒信号(P(Ωc,t))的时变信噪比的估计SNR(k);
通过使用(34)对每个阶n对离散的有限的波数k根据所述信噪比的估计SNR(k)设计的时变维纳滤波器,将所述维纳滤波器的传递函数乘以(34)所述麦克风阵列的逆传递函数以便得到适应传递函数Fn,arrary(k);
使用线性滤波器处理对所述球面谐波或者高保真度立体声响复制表示施加(32)所述适应传递函数Fn,arrary(k),得到适应方向系数
其中,索引n从0到有限阶N,而对于每个索引n,索引m从-n到n。
2.根据权利要求1所述的方法,其中,在没有任何声源的安静的环境中得到所述噪声功率|Pnoise(k)|2使得|P0(k)|2=0。
3.根据权利要求1或2所述的方法,其中,根据在麦克风炭精盒测量的压力Pmic(Ωc,k)通过比较在麦克风炭精盒上的压力的期望值与所测量的在麦克风炭精盒上的平均信号功率估计所述平均源功率|P0(k)|2。
4.根据权利要求1或2所述的方法,其中,在频域中确定阵列的所述传递函数Fn,arrary(k),包含:
使用FFT将系数变换到频域,然后乘以所述传递函数Fn,arrary(k);
实施乘积的逆FFT以得到时域系数
或者,在时域中通过FIR滤波器实施逼近,包含:
实施逆FFT;
实施循环移位;
对得到的滤波器脉冲响应施加锥形窗口以便平滑对应的传递函数;
对于n与m的每个组合实施得到的滤波器系数与系数的卷积。
5.一种用于处理在刚性球上的球面麦克风阵列的麦克风炭精盒信号(P(Ωc,t))的装置,所述装置包含:
部件(31),适用于将表示在所述麦克风阵列的表面上的压力的所述麦克风炭精盒信号(P(Ωc,t))转换为球面谐波或者高保真度立体声响复制表示
部件(33),适用于针对每个波数k,使用从所述麦克风阵列记录的平面波的平均源功率|P0(k)|2以及表示由在所述麦克风阵列中的模拟处理产生的空间不相关的噪声的对应的噪声功率|Pnoise(k)|2,计算所述麦克风炭精盒信号(P(Ωc,t))的时变信噪比的估计SNR(k);
部件(34),适用于通过使用对每个阶n对离散的有限的波数k根据所述信噪比的估计SNR(k)设计的时变维纳滤波器,将所述维纳滤波器的传递函数乘以所述麦克风阵列的逆传递函数以便得到适应传递函数Fn,arrary(k);
部件(32),适用于使用线性滤波器处理对所述球面谐波或者高保真度立体声响复制表示施加所述适应传递函数Fn,arrary(k),得到适应方向系数
其中,索引n从0到有限阶N,而对于每个索引n,索引m从-n到n。
6.根据权利要求5所述的装置,其中,在没有任何声源的安静的环境中得到所述噪声功率|Pnoise(k)|2使得|P0(k)|2=0。
7.根据权利要求5或6所述的装置,其中,根据在麦克风炭精盒测量的压力Pmic(Ωc,k)通过比较在麦克风炭精盒上的压力的期望值与所测量的在麦克风炭精盒上的平均信号功率估计所述平均源功率|P0(k)|2。
8.根据权利要求5或6所述的装置,其中,在频域中确定阵列的所述传递函数Fn,arrary(k),包含:
使用FFT将系数变换到频域,然后乘以所述传递函数Fn,arrary(k);
实施乘积的逆FFT以得到时域系数
或者,在时域中通过FIR滤波器实施逼近,包含:
实施逆FFT;
实施循环移位;
对得到的滤波器脉冲响应施加锥形窗口以便平滑对应的传递函数;
对于n与m的每个组合实施得到的滤波器系数与系数的卷积。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP11306471.1A EP2592845A1 (en) | 2011-11-11 | 2011-11-11 | 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 |
EP11306471.1 | 2011-11-11 | ||
PCT/EP2012/071535 WO2013068283A1 (en) | 2011-11-11 | 2012-10-31 | 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 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103931211A CN103931211A (zh) | 2014-07-16 |
CN103931211B true CN103931211B (zh) | 2017-02-15 |
Family
ID=47143887
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201280055175.1A Active CN103931211B (zh) | 2011-11-11 | 2012-10-31 | 处理刚性球上的球面麦克风阵列的信号的方法及装置 |
Country Status (6)
Country | Link |
---|---|
US (1) | US9503818B2 (zh) |
EP (2) | EP2592845A1 (zh) |
JP (1) | JP6030660B2 (zh) |
KR (1) | KR101938925B1 (zh) |
CN (1) | CN103931211B (zh) |
WO (1) | WO2013068283A1 (zh) |
Families Citing this family (35)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10021508B2 (en) * | 2011-11-11 | 2018-07-10 | Dolby Laboratories Licensing Corporation | 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 |
EP2592846A1 (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 |
US9466305B2 (en) | 2013-05-29 | 2016-10-11 | Qualcomm Incorporated | Performing positional analysis to code spherical harmonic coefficients |
US9769586B2 (en) | 2013-05-29 | 2017-09-19 | Qualcomm Incorporated | Performing order reduction with respect to higher order ambisonic coefficients |
US20150127354A1 (en) * | 2013-10-03 | 2015-05-07 | Qualcomm Incorporated | Near field compensation for decomposed representations of a sound field |
EP2866475A1 (en) | 2013-10-23 | 2015-04-29 | Thomson Licensing | Method for and apparatus for decoding an audio soundfield representation for audio playback using 2D setups |
DE102013223201B3 (de) * | 2013-11-14 | 2015-05-13 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. | Verfahren und Vorrichtung zum Komprimieren und Dekomprimieren von Schallfelddaten eines Gebietes |
US9922656B2 (en) | 2014-01-30 | 2018-03-20 | Qualcomm Incorporated | Transitioning of ambient higher-order ambisonic coefficients |
US9489955B2 (en) | 2014-01-30 | 2016-11-08 | Qualcomm Incorporated | Indicating frame parameter reusability for coding vectors |
JP6243060B2 (ja) * | 2014-03-21 | 2017-12-06 | ドルビー・インターナショナル・アーベー | 高次アンビソニックス(hoa)信号を圧縮する方法、圧縮されたhoa信号を圧縮解除する方法、hoa信号を圧縮する装置および圧縮されたhoa信号を圧縮解除する装置 |
EP2922057A1 (en) | 2014-03-21 | 2015-09-23 | Thomson Licensing | Method for compressing a Higher Order Ambisonics (HOA) signal, method for decompressing a compressed HOA signal, apparatus for compressing a HOA signal, and apparatus for decompressing a compressed HOA signal |
KR102429841B1 (ko) | 2014-03-21 | 2022-08-05 | 돌비 인터네셔널 에이비 | 고차 앰비소닉스(hoa) 신호를 압축하는 방법, 압축된 hoa 신호를 압축 해제하는 방법, hoa 신호를 압축하기 위한 장치, 및 압축된 hoa 신호를 압축 해제하기 위한 장치 |
US9620137B2 (en) | 2014-05-16 | 2017-04-11 | Qualcomm Incorporated | Determining between scalar and vector quantization in higher order ambisonic coefficients |
US20150332682A1 (en) * | 2014-05-16 | 2015-11-19 | Qualcomm Incorporated | Spatial relation coding for higher order ambisonic coefficients |
US10770087B2 (en) | 2014-05-16 | 2020-09-08 | Qualcomm Incorporated | Selecting codebooks for coding vectors decomposed from higher-order ambisonic audio signals |
US9852737B2 (en) * | 2014-05-16 | 2017-12-26 | Qualcomm Incorporated | Coding vectors decomposed from higher-order ambisonics audio signals |
EP3172541A4 (en) * | 2014-07-23 | 2018-03-28 | The Australian National University | Planar sensor array |
TWI584657B (zh) * | 2014-08-20 | 2017-05-21 | 國立清華大學 | 一種立體聲場錄音以及重建的方法 |
KR101586364B1 (ko) * | 2014-09-05 | 2016-01-18 | 한양대학교 산학협력단 | 공간 음향 분할을 통한 동적 방향성 임펄스 응답을 생성하기 위한 방법, 장치 및 컴퓨터 판독 가능한 기록 매체 |
US9747910B2 (en) | 2014-09-26 | 2017-08-29 | Qualcomm Incorporated | Switching between predictive and non-predictive quantization techniques in a higher order ambisonics (HOA) framework |
US9560441B1 (en) * | 2014-12-24 | 2017-01-31 | Amazon Technologies, Inc. | Determining speaker direction using a spherical microphone array |
EP3073488A1 (en) | 2015-03-24 | 2016-09-28 | Thomson Licensing | Method and apparatus for embedding and regaining watermarks in an ambisonics representation of a sound field |
RU2687882C1 (ru) * | 2016-03-15 | 2019-05-16 | Фраунхофер-Гезеллшафт Цур Фёрдерунг Дер Ангевандтен Форшунг Е.В. | Устройство, способ формирования характеристики звукового поля и машиночитаемый носитель информации |
US10492000B2 (en) | 2016-04-08 | 2019-11-26 | Google Llc | Cylindrical microphone array for efficient recording of 3D sound fields |
US11218807B2 (en) * | 2016-09-13 | 2022-01-04 | VisiSonics Corporation | Audio signal processor and generator |
US10516962B2 (en) * | 2017-07-06 | 2019-12-24 | Huddly As | Multi-channel binaural recording and dynamic playback |
CN109963249B (zh) * | 2017-12-25 | 2021-12-14 | 北京京东尚科信息技术有限公司 | 数据处理方法及其***、计算机***和计算机可读介质 |
CN112292870A (zh) | 2018-08-14 | 2021-01-29 | 阿里巴巴集团控股有限公司 | 音频信号处理装置及方法 |
JP6969793B2 (ja) | 2018-10-04 | 2021-11-24 | 株式会社ズーム | アンビソニックスのためのa/bフォーマット変換装置、a/bフォーマット変換ソフトウェア、レコーダー、再生ソフトウェア |
CN110133579B (zh) * | 2019-04-11 | 2021-02-05 | 南京航空航天大学 | 适用于球面麦克风阵列声源定向的球谐波阶数自适应选择方法 |
KR102154553B1 (ko) * | 2019-09-18 | 2020-09-10 | 한국표준과학연구원 | 지향성이 향상된 마이크로폰 어레이 및 이를 이용한 음장 취득 방법 |
CN112530445A (zh) * | 2020-11-23 | 2021-03-19 | 雷欧尼斯(北京)信息技术有限公司 | 高阶Ambisonic音频的编解码方法及芯片 |
CN113395638B (zh) * | 2021-05-25 | 2022-07-26 | 西北工业大学 | 一种基于等效源法的室内声场扬声器重放方法 |
CN113281900B (zh) * | 2021-05-26 | 2022-03-18 | 复旦大学 | 一种基于汉克尔变换与波束传播法的光学建模与计算方法 |
US11349206B1 (en) | 2021-07-28 | 2022-05-31 | King Abdulaziz University | Robust linearly constrained minimum power (LCMP) beamformer with limited snapshots |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101185368A (zh) * | 2005-05-27 | 2008-05-21 | 马丁坎托拉咨询有限公司 | 用于声换能器的组件、***以及方法 |
CN101627641A (zh) * | 2007-03-05 | 2010-01-13 | 格特朗尼克斯公司 | 具有信号处理功能的小器件封装麦克风模块 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7123727B2 (en) * | 2001-07-18 | 2006-10-17 | Agere Systems Inc. | Adaptive close-talking differential microphone array |
US20030147539A1 (en) * | 2002-01-11 | 2003-08-07 | Mh Acoustics, Llc, A Delaware Corporation | Audio system based on at least second-order eigenbeams |
US7558393B2 (en) * | 2003-03-18 | 2009-07-07 | Miller Iii Robert E | System and method for compatible 2D/3D (full sphere with height) surround sound reproduction |
EP1737271A1 (en) * | 2005-06-23 | 2006-12-27 | AKG Acoustics GmbH | Array microphone |
JP4671303B2 (ja) * | 2005-09-02 | 2011-04-13 | 国立大学法人北陸先端科学技術大学院大学 | マイクロホンアレイ用ポストフィルタ |
GB0906269D0 (en) * | 2009-04-09 | 2009-05-20 | Ntnu Technology Transfer As | Optimal modal beamformer for sensor arrays |
EP2592846A1 (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 |
US9197962B2 (en) * | 2013-03-15 | 2015-11-24 | Mh Acoustics Llc | Polyhedral audio system based on at least second-order eigenbeams |
-
2011
- 2011-11-11 EP EP11306471.1A patent/EP2592845A1/en not_active Withdrawn
-
2012
- 2012-10-31 WO PCT/EP2012/071535 patent/WO2013068283A1/en active Application Filing
- 2012-10-31 KR KR1020147015362A patent/KR101938925B1/ko active IP Right Grant
- 2012-10-31 US US14/356,185 patent/US9503818B2/en active Active
- 2012-10-31 CN CN201280055175.1A patent/CN103931211B/zh active Active
- 2012-10-31 JP JP2014540395A patent/JP6030660B2/ja active Active
- 2012-10-31 EP EP12783190.7A patent/EP2777297B1/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101185368A (zh) * | 2005-05-27 | 2008-05-21 | 马丁坎托拉咨询有限公司 | 用于声换能器的组件、***以及方法 |
CN101627641A (zh) * | 2007-03-05 | 2010-01-13 | 格特朗尼克斯公司 | 具有信号处理功能的小器件封装麦克风模块 |
Non-Patent Citations (1)
Title |
---|
"3D Sound field recording with higher order ambisonics--objective measurements and validation of a 4th order spherical microphone";Jerome daniel;《Audio engineering society》;20060523;全文 * |
Also Published As
Publication number | Publication date |
---|---|
KR101938925B1 (ko) | 2019-04-10 |
EP2777297A1 (en) | 2014-09-17 |
KR20140091578A (ko) | 2014-07-21 |
JP2014535231A (ja) | 2014-12-25 |
WO2013068283A1 (en) | 2013-05-16 |
US9503818B2 (en) | 2016-11-22 |
EP2777297B1 (en) | 2016-06-08 |
JP6030660B2 (ja) | 2016-11-24 |
EP2592845A1 (en) | 2013-05-15 |
CN103931211A (zh) | 2014-07-16 |
US20140286493A1 (en) | 2014-09-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103931211B (zh) | 处理刚性球上的球面麦克风阵列的信号的方法及装置 | |
CN104041074B (zh) | 处理用于产生声场的高保真度立体声响复制表示的刚性球上的球形麦克风阵列的信号的方法和装置 | |
CN103856866B (zh) | 低噪微分麦克风阵列 | |
Teutsch et al. | Acoustic source detection and localization based on wavefield decomposition using circular microphone arrays | |
CN103888862B (zh) | 声音捕获*** | |
JP6389259B2 (ja) | マイクロホンアレイを使用した残響音の抽出 | |
CN103339961B (zh) | 用于通过声波三角测量进行空间性选择声音获取的装置及方法 | |
US10659873B2 (en) | Spatial encoding directional microphone array | |
KR100856246B1 (ko) | 실제 잡음 환경의 특성을 반영한 빔포밍 장치 및 방법 | |
WO2017218399A1 (en) | Spatial encoding directional microphone array | |
Sakamoto et al. | Sound-space recording and binaural presentation system based on a 252-channel microphone array | |
US10021508B2 (en) | 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 | |
CN110907892A (zh) | 一种球麦克风阵列语音信号到达角估计方法 | |
EP2757811B1 (en) | Modal beamforming | |
Hoffmann et al. | Theoretical study of acoustic circular arrays with tangential pressure gradient sensors | |
JP2014045440A (ja) | 音源別音声強調装置、方法、プログラム | |
JP7182168B2 (ja) | 音情報処理装置及びプログラム | |
Kordon et al. | Optimization of spherical microphone array recordings | |
Liu et al. | Virtual audition technology for head-mounted electronic sentinel system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C41 | Transfer of patent application or patent right or utility model | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20160714 Address after: Amsterdam Applicant after: Dolby International AB Address before: I Si Eli Murli Nor, France Applicant before: Thomson Licensing SA |
|
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |