CN103929706B - 用于测量声源产生的直达声的装置及方法 - Google Patents
用于测量声源产生的直达声的装置及方法 Download PDFInfo
- Publication number
- CN103929706B CN103929706B CN201410079512.1A CN201410079512A CN103929706B CN 103929706 B CN103929706 B CN 103929706B CN 201410079512 A CN201410079512 A CN 201410079512A CN 103929706 B CN103929706 B CN 103929706B
- Authority
- CN
- China
- Prior art keywords
- wave
- produce
- coefficient
- rad
- direct
- 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
- H04R29/00—Monitoring arrangements; Testing arrangements
- H04R29/001—Monitoring arrangements; Testing arrangements for loudspeakers
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H3/00—Measuring characteristics of vibrations by using a detector in a fluid
- G01H3/10—Amplitude; Power
- G01H3/12—Amplitude; Power by electric means
- G01H3/125—Amplitude; Power by electric means for representing acoustic field distribution
-
- 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
-
- G—PHYSICS
- G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
- G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
- G03H3/00—Holographic processes or apparatus using ultrasonic, sonic or infrasonic waves for obtaining holograms; Processes or apparatus for obtaining an optical image from them
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Otolaryngology (AREA)
- Signal Processing (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明提供了—种用于测量受声学环境噪声源Q1和声学边界(例如,房间墙壁)反射的影响的受验声源(例如,扬声器)产生的直达声wrad的装置及方法。采集装置14通过传感器1测量扫描范围Gm内多个测量点rm处声场的状态变量pt(rm),并且产生扫描数据集根据该数据集,分析器55确定与波方程的解的展开函数相关的系数识别器16利用扫描数据集为分析器55产生参数信息P,这是分离直达声wrad与房间反射wref以及受验声源表面散射的其它波wsec的基出。外推器(11)利用波展开的系数预测扫描范围Gm外的任意点直达声wrad的状态变量
Description
技术领域
本发明涉及—种用于确定受声学环境噪声源和声学边界(例如,房间墙壁)反射的影响的受验声源(例如,扬声器)产生的直达声的装置及方法。
因此,对直达声的评估需要不受限制的声场条件,这在消声室(anechoic room)中能够实现。对用在专业应用(例如,线阵列)中的大型扬声器系统,其远声场响应的测量需要具有被特殊处理过的墙壁的大型消声室,以确保低频声音被充分吸收。这里需要在非消声环境中进行准确的测量。
背景技术
下面将根据多个作者的近期研究活动来描述本发明及本发明的潜在问题,但是,本发明并不限定于该类型的应用。
M.Melon等人在第55卷第12期音频工程师协会杂志的第1077-1091页(2007年)中的“四种低音炮测量技术的比较”的论文中,探讨了多种可替换的方法。可通过对后到达房间的反射wref的脉冲响应进行加窗来分离受验声源Q0产生的高频直达声wrad。这种方法不能用于低频直达声,因为低频直达声的相应波长没有小到可比于声源Q0和周围边界面之间的最小距离。
G.Weinreich和E.Amold在第68卷第2期的美国声学学会杂志的第404-411页(1980年)发表的“用于测量声辐射场方法”的论文中,提出了一种根据在两个同心球面中测量的声压的阐述,测量受验声源Q0直达声的全息输入输出场分离方法(IOFS)。
Melon在2012年4月23-27日法国南斯(Nantes,France)举办的2012声学学报年会会刊中发表的“采用场分离方法的低音炮测量:p-p和p-v公式的比较”论文中,将这种IOFS方法应用到扬声器中。如果表面S0没有小到可比于扫描面S1的大小,受验声源Q0表面S0反射和散射的声波会引起测量误差。
为了克服这个问题,C.-X.Bi在2012年4月23-27日法国南斯(Nantes,France)举办的2012声学学报年会会刊中发表的“利用球面波叠加方法进行自由场恢复”一文中,提出了—种全息原始声场分离方法(PSFS),采用声阻抗和表面S0的形状作为附加的输入信息。由于扬声器设计中采用的形状很复杂并且材料很多样,因此,不容易获取这些信息,并且因此实际上该方法通常不实用。
C.Langrenne在1997年Universite du Maine他的论文“Methodes deregularisation du probleme inverse acoustique pour l’indentification desources en milieu confine et pertube”中,提出了—种用于将直达声与表面S0的反射和散射波进行分离的扰动方法。这种方法需要额外测量多个测量点rm的声压pt(rm),其中m=1,...M,与此同时关闭受验声源Q0并且激活外扫描表面S2与房间边界SB之间的空间内的多个位置re处的扰动声源Q2(re),其中e=1,...,E。由于需要大量的测量值(E和M的乘积)来确定传递矩阵H,因此,这种扰动方法非常耗时。
已知的采用球面谐波的场分离方法,将频率限定在截止频率fG≈55N/r1以下,其中展开的最大阶数(order)N足够大,以足以接近半径为r1的球面扫描面S1围绕的场。
Melon通过将受验声源设置在声学硬地板上,扫描半空间中的球面上的声压。扫描面的中心对应于球面坐标系统的原点,并且作为根据球面谐波进行波展开的展开点r0。因此,展开点r0与受验声源的声学中心并不相同。利用安装在大型封闭空间内的多个换能器对扬声器系统进行全息测量,需要与大量测量点相关联的波展开具有大的阶数N,从而以足够的准确度描述声源的近场声压。
这就需要找到更简单/更可靠地可以测量受房间反射和环境噪声影响的受验声源Q0直达声的方法,以克服已知方法的缺点。
发明内容
本发明提供一种装置及方法,用于测量围绕受验声源Q0的多个测量点rm的声压pt(rm)、声强或者声场的任意其它状态变量,测量点rm通过围绕受验声源Q0的扫描范围Gm的扫描矢量R=[r1,r2,...rm,...rM]来定义。扫描范围Gm的内边界面Sin可以是紧密围绕受验声源而不切割(cutting)其表面S0的球面、圆柱面或者其它任意合适的形状。面Sin与S0之间应当存在最小空气体积(air volume),用于评估受验声源Q0的近场以及实现最高截止频率fG,这里适用IOFS方法。近场测量的进一步优点在于直达声wrad的幅值大于声学边界SB与扫描范围Gm外表面Sout之间的外部声源Q1产生的房间反射wref和环境噪声的幅值。
测量的状态变量汇总为扫描数据集并且通过输出声波wout与输入声波win的波展开来描述
采用系数
以及展开函数ΨGm,out和ΨGm,in。
本发明公开了—种新的PSFS方法,用于将直达声wrad与表面S0散射的第二声波wsec进行分离。PSFS方法利用通过声学扫描处理所提供的信息。该方法不能确定根据声阻抗和S0的几何形状传播到表面S0的速度根据本发明,测量描述输入声波win波系数Cin与第二声波wsec波系数Csec之间的关系的传递函数f
Csec=f(Cin)。 (3)
当直达声被削弱并且仅仅剩下房间反射波、驻波以及其它混响声波分量时,利用脉冲响应的后面部分(1ate part),从扫描数据集推导出传递函数f的参数P。任意测量点rm的混响声波部分的声压
可以通过线性滤波器将激励u(t)与窗口脉冲响应
进行卷积来确定。
使用展开函数,利用第二IO场分离方法确定输入和输出混响声波的波系数
这些波系数是确定第二声波的波系数的基础
该第二声波的透射参数为
可供造择地,第二波系数
可利用透过系数(transparency factor)
来确定。
透射参数kx和透过系数kt都取决于受验声源Q0相对于边界面SB的特殊的位置和方向,这与在扰动方法中传递矩阵H不依赖于边界面SB正好相反。
受验声源Q0发射的直达声wrad的波系数
可利用IOFS方法和新的PSFS方法的结果来确定。
本发明进一步的特点在于通过IOFS误差矢量
来监测IOFS方法的使用可靠性,该IOFS误差矢量描述了扫描范围Gm内的测量点rM的测量声压与扩展声压之间的不匹配度。
通过输入声波与混响分量之间的相关系数
可对新的PSFS方法的有效性进行和评估。
如果IO误差矢量的范数(norm)||eio||超过预先没定的阈值或者相关系数cx很小,那么通过场分离方法确定的波系数就非常不准确。这发生在高于由最大阶数N和扫描范围Gm的大小共同决定的截止频率fG时。
根据本发明,确定在频率方面直达声的波系数高于交叉频率fC,通过可替换的方法,产生包含声压
的滤波数据集声压作为激励u(t)和由时间窗Wdir选通的脉冲响应的卷积。这里给出了直达声的波系数
可替换的方法的有效性可通过窗口误差矢量
来评估,窗口误差矢量描述了扫描范围Gm内的外表面Sout的测量声压与扩展声压之间的偏离。当窗口Wdir不能将直达声wrod与房间反射wref分离开时,窗口误差范数||ew||就上升到较低频率。
通过在可替换的方法中选择产生最小误差矢量范数||ew||和||eio||的—种方法,可以确定可供选择方法中的最优交叉频率fC。
相关联的误差系数展开的波系数
能够被用于外推出相对误差
以及扫描范围Gm外的任意观测点r的所预测的直达声的有效性。
如果采用多个相同性能的麦克风,通过同步多通道采样,在所有点rm,其中,m=1,...,M,同时测量声压pt(rm),那么,场分离方法仅能够分离原始声波wrad与环境噪声源Q1产生的非相干噪声wn。ise。
—种成本有效的替换是利用单个麦克风进行连续扫描处理,因为传感器的幅值和相位响应所引起的任意线性失真都不会影响全息处理。此外,与角分辨率相关的测量点最大数值和展开阶数仅仅通过扫描处理时间和对扫描数据集的数值处理来限定。
连续扫描处理需要—种新技术,用于对被与激励u(t)不相干的环境噪声wnoise损坏的无效测量值进行检测。根据本发明,在扫描范围Gm内,在比测量声压pt(rm)的第一麦克风距离受验声源Q0更远的距离处,第二麦克风测量声压pt(ra)。如果信噪比
低于预先设定的阈值,那么测量声压pt(rm)就是无效的。通过重复测量相同激励u(t)下的同一测量点rm并且对测量的声压pt(rm)取平均,可以提高SNR(信噪比)。通过取平均,排除—部分被脉冲扰动损坏并且具有非常小的SNR的声压信号pt(rm),以及可以仅仅将多个测量值中的有效部分整合为有效的响应,存储到所述扫描数据集中。
由于连续扫描非常耗时,因此这里要求需要的测量点rm的数量M最小。但是,测量点rm的数量M,其中,m=1,...,M,限定了展开阶数N和测量指向性的角分辨率。例如,根据球面参考函数的展开需要至少M=2(N+1)2个测量点,用于分离输入和输出的声场并且用于确定与阶数N相关的直达声的波系数。
与现有技术相反,全息测量被实现为自适应技术,研究扫描分析与全息分析之间的互相关性。初步全息分析利用在中采集的第—扫描数据,用于确定受验声源Q0的几何性能和声学性能,并且适应性地确定冗余测量点的最优位置。
根据受验声源Q0表面S0的形状分别推导出扫描范围Gm的内表面Sin和外表面Sout。例如,小口径多路扬声器适合的扫描面是圆柱面。
在η(ω,n)≤η0≤η(ω,n+1)的范围内,寻找阶数n来确定,作为频率ω的函数的展开阶数N(ω)=n(ω)的最优值,其中系数cn,m对总的声波能量的相对贡献
刚好在预先设定的界限η0之下。
扫描范围Gm内测量点rm的密度和分布以及它们与原点的距离rm取决于受验声源Q0的声学性能。例如,具有高指向性的专业扬声器需要在扬声器的主辐射方向比扬声器后侧具有更多的测量点。对于这些装置,有利地是,可将扫描范围Gm划分为多个子区间G1,G2,...,这些子区间的测量点密度各不相同,对应于要求的角分辨率。通过插入低阶的波场展开产生的点,可以扩大包含低密度测量点的扫描数据集。所有子区间的测量数据和实际数据者将结合到随后的高阶波展开中并受其支配,其中,系数的总数量(N+1)2/2可以大于测量点的总数量M。
展开点r0(ω)的位置对于以最低阶数N产生可接受拟合误差的波扩展的快速收敛至关重要。受验声源Q0的声学中心是一个好的选择,其在大多数采用多个换能器的扬声器系统中是频率ω的函数。通过由脉冲响应推导出的群延迟、扫描范围Gm内的最大声压、波展开的最小拟合误差、或者通过对波展开系数的解读,可以检测出声学中心。例如,通过利用
球面波展开的零阶和一阶系数就给出了对与最优展开点对应的声学中心的很好的估值。与声学中心相关的信息还可用于在连续扫描处理中,优化扫描范围Gm的形状以及优化测量点在扫描矢量R[l+1]中的位置。
根据本发明,通过利用指向性图案的轴对称和镜像对称信息以及扫描数据集的全息分析提供的对称轴取向信息A[l],可以进一步减少测量点的数量。例如,具有大致轴向对称指向性的声源,可以通过球面展开函数的子集、以及由显著较少数量M的测量点就能确定的系数cn,0,d,其中n=0,1,...,N,进行建模,其中M>N。对称信息A[l]已经能够根据在连续扫描处理过程中采集的第一测量数据推导出来,并且可用于优化扫描处理和全息分析。
如果受验声源不对称,那么通过去掉(remove)对直达声wrda辐射声功率贡献很小的标识符(identification)中的展开函数和相关的波系数,就可以减少测量点的总数量M。通过贡献率
可以评估每个系数Cj,rad在其它波系数中的贡献。
根据本发明,利用展开函数的稀疏集,通过迭代处理对扫描数据集进行全息分析,这里在每个步骤l中将具有低贡献率rj[j]的系数Cj,rad[l]替换为高阶展开函数的系数,这可以提高展开的准确度和角分辨率。
所附的权利要求书中公开了(setforthin)本发明的具体实施例。参照下文中描述的实施例,本发明的这些方面和其它方面将会是显而易见的并且将得到阐述。
附图说明
参照附图,仅仅通过示例,将描述本发明进一步的细节、方面和实施例。为了简明和清楚地描述了附图中的元件(elements),并且这些元件没有必要按比例画出(drawn toscale)。
图1示意性地示出了采用IOFS方法的直达声全息测量。
图2示出了基于表面声阻抗的PSFS方法,用于分离散射波。
图3示出了基于外部声源Q2扰动的PSFS方法。
图4示意性地示出了基于改进PSFS方法的全息测量。
图5示出了基于混响声产生传递参数P的识别器(identifier)的实施例。
图6示出了使用IOFS方法和改进PSFS方法的分析器的实施例。
图7示出了确定最优展开点的测量系统的实施例。
图8示出了采用具有测量点密度的两个扫描范围G1和G2的测量系统的实施例。
图9示出了使用所选择的展开函数的分析器的实施例。
在所有的附图中,相同的或者至少具有相同功能的图形元件,特征及信号都具有相同的附图标记,另有明确说明的除外。
具体实施方式
图1示出了—种IOFS方法,用于测量处于试验2下的声源Q0在受到发生器8产生的激励信号u(t)的激励时在观测点r处产生的直达声采集装置4利用传感器1测量在靠近处于试验2下的声源在两个球面扫描面S1和S2上的声压pt(rm)或其它声学状态变量,并且分别利用存储装置7和9收集两个扫描数据集和中的声压信号。
直达声wrad在边界SB的房间反射产生的混响声wref以及外部声源Q1产生的噪声wnoise,进入到输入波win中的扫描面,并且通过反射、衍射和散射在处于试验2下的装置的表面S0上产生第二声波wsec。第二声波wsec与直达声wrad的叠加在两个扫描面形成输出波wout。
在频率为ω的声场中,点r处的总声压
被展开为分别对应于输出波毗ut和输入波w加的声压分量和两个分量被展开在包含展开函数Ψj,out和Ψj,in的序列中,展开函数Ψj,out和Ψj,in的权重系数为Cj,out和Cj,in。
波方程在笛卡尔、圆柱面以及球面坐标中的解,最佳候选是展开函数。例如,展开函数描述了在球面坐标中,球面谐波的角度相关性以及在第一类和第二类汉克尔函数和在展开点的径向相关性。展开点r0是坐标系的原点。
输出波wout的声压
包括直达声的声压和第二波wsec的声压
可替换的展开描述了总声压
是由系数表示的驻波wsw和系数表示的冗余波wex的叠加。
冗余声压
包括受验声源Q0产生的直达声wrad和在表面S0散射和反射的声波wscat。
测量点rm,其中,m=1,...,M,在扫描面S∈{S1,S2,...}上的总声压可被汇总为扫描数据集
其中展开利用声压矢量
波系数
以及展开矩阵
具有表示声波分量的指标
d∈{t,in,out,ex,scat,sw,sec,rad} (31)
以及表示声源的指标
Q∈{Q0,Q1,Q2,...}。 (32)
方程式(27)中的展开对应于波系数之间的以下关系:
在IOFS方法13中波系数和通过
或者通过在球面扫描面S1和S2进行积分而得到,由E.Williams描述在学术出版社1999年出版的傅立叶声学一书的第7.4章中。
外推器11利用分析器55中的冗余声波系数确定扫描面外任意观测点r处直达声的声压
图2示出了结合IOFS方法和PSFS方法的全息测量。通过IOFS方法13估算系数和将该系数提供给接下来的PSFS方法10,在直达声压的展开中产生系数
估算与散射声波wscat相关的系数
利用在处于试验2下的声源的表面S0上的速度
驻波的展开矩阵
输出波的展开矩阵
以及表面S0的声学导纳Y。
图3示出了—种扰动方法,利用放置在扫描范围Gm外表面Sout与房间边界SB之间的空间内的多个点re处的外部声源Q2,其中e=1,...,E。开关27将发生器8产生的激励提供给受验声源Q1或外部声源Q2。开关29将声压信号pt(rm)提供给存储装置7和9产生扫描数据集或者提供给存储装置15和16产生取决于位置re的扫描数据集
IOFS模块19具有产生波系数
波系数被提供给子系统21,以确定传递矩阵
传递矩阵具有
以及
IOFS模块18根据存储装置7和9提供的扫描数据集和产生波系数
接下来的PSFS模块23利用IOFS方法和传递矩阵H的结果产生直达声的系数
图4示意性地示出了根据本发明的全息测量方法的一个实施例。采集装置14产生激励u(t),使用声源Q0激励处于试验2下的装置,并且将传感器1放置在位于扫描范围Gm内的测量点rm上。与现有技术相反,扫描范围Gm取决于处于试验2下的装置的表面S0的形状。向识别器16提供扫描数据集识别器16将产生描述处于试验2下的装置的声学性能的参数P。根据该参数P以及扫描数据集分析器55产生表示直达声wrad的波系数和表示直达声的预测误差的误差系数Erad。根据该信息,外推器11根据方程式(35)产生直达声的声压并且根据方程式(18)产生观测点r处的相对误差erad(r)。
图5示出了利用扫描数据集的识别器16的实施例,这与现有技术仅仅利用输入信息相反。线性滤波器39根据方程式(4)产生经过滤波的扫描数据集接下来的IOFS模块19根据方程式(6)产生分别表示混响声波的输入和输出场的波系数和将波系数和作为参数P提供给分析器55中的PSFS模块24,分析器55利用通过IOFS模块18从扫描数据集产生的波系数和根据方程式(7)、(9)、(11)产生原始声波的波系数
根据扫描数据集线性滤波器111根据方程式(14)产生包括具有小部分时间延迟的部件的滤波扫描数据集自由场展开器113根据方程式(15)产生波系数和误差矢量ew,作为参数P提供给分析器55。分析器55包括估算器117,通过方程式(17)根据来自自由场展开器113的误差矢量ew和来自IOFS模块18的误差矢量eio产生交叉频率fC和误差系数Erad。交叉混合器115接收交叉频率fC和波系数并且根据方程式(15)产生波系数
图6示出了根据本发明的PSFS模块24的实施例。相关器97相据方程式(8)基于波系数和产生透射参数kx,其中透射参数kx描述了总的输入波win的指向性与输入波win混向部分的指向性的一致性。
比较器99根据方程式(10)基于波系数和产生透过参数kt,其描述了被扫描范围(Gm)封闭的空间的声学透过性。合成器93根据方程式(9)基于总的输入波系数透过参数kt以及透射参数kx产生与第二声波wsec的展开相关的第二波系数加法器95根据方程式(15)基于总的输出波系数和第二波系数产生直达波系数
图7示出了应用到多通路扬声器系统2中的测量系统的实施例,多通道扬声器系统2包括低音炮41、中端换能器43和高音炮45。扫描范围Gm的内、外表面Sin和Sout分别都是圆柱面,以适应扬声器的小口径形状。通过定位装置49,利用圆柱面坐标r,z,φ中的三个执行器11、47、51将传感器1放置在测量点rm处。第二传感器76放置在扫描范围Gm外部,其与受验声源Q0的距离比与传感器1的距离远,用于监测表示由噪声源Q1产生的声学干扰wnoise的坏境声学信号pt(ra)。噪声识别器77将扫描信号pt(rm)与环境声学信号pt(ra)进行比较,并且检测出被噪声源Q1破坏的无效测量值。在将扫描信号pt(rm)的有效部分存储在扫描数据集之后,噪声识别器77在当前传感器位置rm处进行重复测量,直到完成有效部分的收集。
在第一步l=1中,分析器55嵌入到迭代处理中,根据初始扫描数据集开始确定波系数展开点识别器59检测处于试验2下的装置的声学中心z0(ω),声学中心是频率ω的函数,并且确定展开点r0[l+1]的坐标。对称性识别器60产生对称性参数A[l+1],,表示声场的轴对称和反射对称以及与主辐射方向相关的处于试验2下的装置的指向性。根据被确定的展开点r0[l+1]和波系数的对称性参数A[l+1],通过引入设置在扫描范围Gm内的最优位置处的附加测量点,产生器53扩展扫描矢量R[l+1]。考虑到处于试验2下的装置被确定的性能,根据扩展的扫描矢量R[l+1],定位装置49进行自适应扫描处理。因此,当利用的测量点数量最少时,扫描数据集可为波展开提供足够的信息。
图8示出了采用具有不同测量点密度的两个扫描范围G1和G2的分析器55的实施例。采集装置14收集第一扫描数据集和第二扫描数据集中的测量传感器输出pt(rm)。第一子分析器63通过方程式(15)根据第二扫描数据集产生第二波系数其中,第二扫描范围G2中测量点的数量M2限制了第二波系数的最大阶数N2。内插器(65)根据外推法,通过方程式(35)根据第二波系数产生插入扫描数据集插入扫描数据集中元素的数量M2s大于第二扫描范围G2中测量点的数量M2。加法器67通过将第—扫描数据集与插入扫描数据集合并,产生统一的数据集其中第二扫描范围G2中的采样点密度等于第—扫描范围G1中的测量点密度。第二分析器69通过方程式(15)根据统一的数据集产生直达声的波系数其中展开的阶数Nrad大于第二波系数的阶数N2。但是,插入并不能增加第二扫描范围G2,中被确定的指向性图案的分辨率。
图9示出了分析器55可替换的实施例,用于根据测量点的最小数量M得到最大分辨率的指向性图案。估算器81利用在第l步中确定的展开矩阵Ψd,rea[l],通过方程式(15)产生缩减波系数
根据缩减波系数反向转换器85利用预先设定的选择矩阵S[l]产生迭代波系数
其中,所述迭代波系数中元素的数量大于所述缩减波系数中元素的数量。
根据迭代波系数估算器87根据方程式(22)产生包含贡献率rj[l]的贡献矢量r[l]。选择器89利用贡献矢量r[l]产生选择矩阵S[l],并且产生与展开函数相关的系数的最高阶数N[l+1],其中系数Cj,rad被设定为0,如果其贡献率rj[l]低于临界阈值r0。
产生器84产生完全展开矩阵Ψd[l+1],要考虑到
-展开函数的最大阶数N[l+1]和/或
-所述扫描范围(Gm)内的测量点(rm)的数量和位置R[l]和/或
-表示受验声源Q0声学中心的位置的展开点r0[l]和/或
-表示受验声源Q0辐射的主方向的方向信息A[l]和/或
-受验声源Q0产生的声场的旋转和/或反射对称性A[l]。
根据完全展开矩阵Ψd[l+1]以及选择矩阵S[l+1],转换器83产生缩减展开矩阵的更新值
Ψd,red[l+1]=S[l+1]Ψd[l+1] (48)
其中,缩减展开矩阵Ψd,red[l+1]中元素的数量小于完全展开矩阵Ψd[l+1]中元素的数量。缩减展开矩阵Ψd,red[l+1]是包含缩减数量的系数的稀疏波展开的基础,系数可通过限定数量的测量点进行评估。
如果在迭代中两个迭代步骤中迭代波系数的值之间的差异低于预先设定的阈值,那么,根据迭代波系数控制器91产生直达波系数
本发明的优点
本发明测量在环境噪声的影响下,由受验装置Q0在有回响的声学环境中产生的直达声。由于免去了消声室,因此减少了成本并且在扬声器及其它声学装置的发展中给出了更好的灵活性。新的测量方法提供了描述扫描范围Gm之外的任意点产生的声场的综合数据集近场信息对于评估移动电话、笔记本电脑以及其它个人声学装置很重要。相比于已知的现有技术的方法,本发明的其它益处包括测量的简化,稳健以及测量速度的提高。新的PSFS方法免去了耗时的扰动,并且无需受验装置Q0的表面的形状和声学性能的信息。该测量技术利用在扫描数据集中提供的冗余信息,并且通过相对误差测量对通过任意观测点测量结果的一致性和准确性进行了描述。本发明利用与短的测量时间相关的最小数量的测量点,进行了场分离和波展开。
在前面的说明书中,已经参照本发明特定实施例示例,对本发明进行了描述。但是,显而易见的是,在不脱离本发明更广泛的精神和范围的情形下,这里做出了各种修改和变型,如在所附的权利要求中所提出的。例如,连接可以是适合于从各自节点、单元或装置传输信号,或者传输信号到各自节点、单元或装置的任—种连接类型,例如,通过中间装置。因此,除非是不言而喻的或者另外声明的,连接可以是例如直接连接或者间接连接。
由于执行本发明的装置,其中大部分,者是由本领域技术人员已知的电子元器件和电路组成的,因此,为了理解和了解本发明的基本概念,而不是用于模糊或分散本发明的教导,电路及其元器件的细节将不会比前述必须探讨的有更深的程度。
上述一些实施例,如果适用的话,可以利用一系列不同的电路元器件执行。例如,附图中提供的示例性拓扑结构以及对其的探讨仅仅是为讨论本发明的不同方面提供的有用参考。当然,为了进行讨论,已经简化了对拓扑结构的说明,并且这仅仅是可用于本发明的许多不同类型的合适拓扑结构中的—种。本领域技术人员需要认识到,逻辑模块之间的边界仅仅是说明性的,并且可替换的实施例可以是逻辑模块或电路元件的合并,或者对各种逻辑模块或电路元件进行可替换的功能分解。
因此,需要理解的是,这里描述的结构仅仅是示例性的,并且事实上可以采用完成相同功能的许多其他总体结构。简言之,但毫无疑问,完成相同功能的元器件的任意配置必须是有效“连接的”,从而使得能够完成要求的功能。因此,这里完成特定功能的任意两个元器件组合可被认为是相互“连接”,使得可以完成要求的功能,不管是总体结构还是中间元器件。同样,任意两个元器件的如此连接也可被认为是相互之间“可操作地连接”,或者“可操作地耦合”,以完成要求的功能。
而且,本发明并不限定于在不可编程硬件中执行的物理装置或单元,而是也可应用到可编程装置或单元中,以能够根据合适的程序代码执行要求的装置功能。此外,装置可以是物理上有多个装置的分布,而功能上仅仅作为单个装置操作。功能上形成分离的装置可以被集成到单个物理装置中。
在权利要求中,在括号中放置的附图标记不能构成对权利要求的限定。词语“包含,,并没有排除权利要求中列出那些之外的其它元件或步骤。此外,这里使用的术语“一”或“一个”,限定为一个或多于一个。而且,权利要求中引入惯用语,例如“至少一个”以及“一个或多个”,的使用,并不应该构成暗示,在其它权利要求元件中引入非限定性冠词“一”或“一个”,限定了对仅仅包含一个这种元件的发明引入包含这种权利要求元件的任意特定权利要求,即使当同一权利要求中引入了引入性短语“一个或多个”或“至少一个,以及非限定性冠词例如“一”或“一个”。同样地,这在对限定性冠词的使用中也是有效的。除非另外声明,术语例如“第一”以及“第二”都用于任意区分这些术语描述的元件之间区别。因此,这些术语并不一定意在表示这些元件的临时的或其它的优先级。在相互不同的权利要求中仅仅通过事实列举的某些测量并不表示不可以采用这些测量的组合的优点。权利要求中的方法步骤的页序并不是实际执行步骤的顺序的偏见,除非在权利要求特别地列举出来。
本领域技术人员将要意识到的是,附图中说明的元件是为了简化和清楚,并且没有必要按比例绘制。例如,仅仅使用选择的元件用于改进对功能以及本发明不同实施例中这些元件的结构的理解。而且,为了便于对本发明这些不同实施例有一个较为简要的理解,在商业上可行的实施例中有用的或必须的这些常见但很好理解的元件大部分没有描述。需要进一步意识到,在描述的方法中的某些行为和/或步骤可以特定发生顺序进行描述或描绘,而本领域技术人员需要理解这些关于顺序的说明并不是实际需要的。还需要理解的是,除非这里另外提出了其他特定的意思,本发明说明书中使用的术语和表述具有通常的含义,因为它们依照了这些术语和表述在其相应的各自领域的调查和研究。
Claims (30)
1.一种用于确定由受验声源(Q0)产生的直达声(wrad)的装置,包括:
采集装置(1,14,77),通过测量多个测量点(rm)处围绕受验声源(Q0)的声场的状态变量(pt(rm))产生扫描数据集其中测量点(rm)设置在扫描范围(Gm)内,并且扫描数据集描述直达声(wrad)与至少一个以下声波分量的叠加:
-由环境噪声源(Q1)或由在外部边界(SB)上反射的直达声(wrad)所产生的输入声波(win),
-由所述输入声波(win)与受验声源(Q0)的表面(S0)的相互作用所产生的第二声波(wsec);
识别器(16),产生与扫描数据集的展开相关的后输入声波系数和后输出声波系数所述系数确定了各展开函数(Ψ)的权重,所述展开函数是波方程的解,所述识别器(16)还包括了:
-滤波器(39),根据所述扫描数据集产生滤波扫描数据集其中保留扫描数据集中的混响声波部分,并且削弱扫描数据集中的直达声波部分(wrad);
-IO场分离模块(19),根据滤波扫描数据集产生与输入声波(win)的混响部分的展开相关的所述后输入波系数并且在直达声(wrad)被削弱情况下产生与第二声波(wsec)的后部分的展开相关的后输出波系数
分析器(55),根据扫描数据集产生与直达声(wrad)的展开相关的直达波系数其中分析器(55)包括:
-IO场分离模块(18),根据扫描数据集产生与总输入声(win)的展开相关的总输入波系数并且产生与总输出声的展开相关的总输出波系数所述输出波系数包含第二声波(wsec)和直达声波(wrad);以及
-PS场分离模块(24),根据总输入波系数和总输出波系数以及后输入波系数和后输出波系数产生与直达声的展开相关的直达波系数其中所述直达波系数确定了各展开函数(Ψ)的权重,所述展开函数(Ψ)是波方程的解;
外推器(11),根据直达波系数产生预测状态变量(prad(r)),所述预测 状态变量(prad(r))描述扫描范围(Gm)的内边界(Sin)外侧预先设定的观测点(r)处的直达声场(wrad)。
2.根据权利要求1的装置,其中,所述PS场分离模块(24)包括:
比较器(99),根据后输入波系数和后输出波系数产生透过参数(kt),所述透过参数(kt)描述被扫描范围(Gm)封闭的空间的声学透过性;
合成器(93),根据总输入波系数和透过参数(kt),产生与受验声源(Q0)产生的第二声波(wsec)的展开相关的第二波系数以及
加法器(95),根据总输出波系数和第二波系数产生所述直达波系数
3.根据权利要求1的装置,其中,所述PS场分离模块(24)包括
相关器(97),根据总输入波系数和后输入波系数产生透射参数(kx),所述透射参数(kx)描述总输入声波(win)的指向性与输入声波(win)混响部分的指向性的一致性;
合成器(93),根据总输入波系数和透过参数(kt),产生与受验声源(Q0)产生的第二声波(wsec)的展开相关的第二波系数以及
加法器(95),根据总输出波系数和第二波系数产生所述直达波系数
4.一种用于确定由受验声源(Q0)产生的直达声(wrad)的装置,包括:
采集装置(1,14,77),通过测量多个测量点(rm)处围绕受验声源(Q0)的声场的状态变量(pt(rm))产生扫描数据集其中测量点设置在扫描范围(Gm)内,并且扫描数据集描述直达声(wrad)与至少一个以下声波分量的叠加:
-由直达声(wrad)在外部边界(SB)上反射产生的输入声波(win),
-由所述输入声波(win)和受验声源(Q0)的表面(S0)相互作用产生的第二声波(wsec);
识别器(17),包含:
-滤波器(111),根据扫描数据集产生滤波扫描数据集其中所述滤波器保留具有小于预先设定值的时间延迟的直达声(wrad)并削弱具有大于预先设定值的时间延迟的信号分量;以及
-自由场展开器(113),产生与滤波数据集的展开相关的滤波系数 其中所述滤波系数确定了各展开函数(Ψ)的权重,展开函数(Ψ)是波方程的解,并且所述展开方式是一种对信号频率高于预先设定的截止频率(fg)的直达声(wrad)的有效表述;
分析器(55),根据扫描数据集产生与直达声(wrad)的展开相关的直达波系数所述分析器(55)包括至少一个以下部件:
-场分离模块(18;24)在低于指定截止频率(fg)的频段分离扫描数据集中的直达波(wrad),并且产生与直达波(wrad)的展开相关的分离直达波系数其中所述分离直达波系数确定了各展开函数(Ψ)的权重,展开函数(Ψ)是波方程的解;
-交叉混合器(115),为信号频率低于交叉频率(fC)的直达波系数分配分离直达波系数并且为信号频率高于或等于所述交叉频率(fC)的直达波系数分配滤波系数
以及外推器(11),根据直达波系数产生预测状态变量(prad(r)),所述预测状态变量(prad(r))描述扫描范围(Gm)的内边界(Sin)外侧预先设定的观测点(r)处的直达声(wrad)。
5.根据权利要求4的装置,其中
所述分析器(55)包括估算器(117),根据场分离模块(18;24)或识别器(16)提供的估算信息(eio;ew),产生所述截止频率(fg),所述截止频率(fg)在所述交叉混合器(115)中用作交叉频率(fC);
-自由场展开器(113),根据滤波扫描数据集产生所述估算信息(ew);其中所述估算信息(ew)描述输入声波(win)的幅值。
6.根据权利要求4的装置,
所述分析器(55),根据扫描数据集产生与直达声(wrad)的展开相关的直达波系数其中,分析器(55)接收描述由受验声源(Q0)所产生的声场的性能的场信息(A[l+1];r0[l+1]);
识别器(59,60),根据直达波系数或扫描数据集产生所述场信息(A[l+1];r0[l+1]),所述场信息(A[l+1];r0[l+1])包含其中至少一个:
-与受验声源(Q0)的声学中心的位置相关的展开点(r0[l+1]),
-与受验声源(Q0)的主辐射方向相关的方向信息(A[l+1]),
-受验声源(Q0)所产生的声场的旋转对称信息(A[l+1]),
-受验声源(Q0)所产生的声场的反射对称信息(A[l+1])。
7.根据权利要求6的装置,其中所述采集装置包括以下之一:
至少一个扫描传感器(1),每个扫描传感器测量扫描传感器(1)的当前位置(rm)处的声场的状态变量(pt(rm));
扫描产生器(53),根据场信息(A[l+1];r0[l+1])产生扫描矢量(R[l+1]),其中扫描矢量包括至少一个附加测量点(rm)的位置,其给出了有关受验声源(Q0)产生的直达声的特定信息;
定位装置(49;47;11;51),将扫描传感器(1)移动到由扫描矢量(R[l+1])限定的附加测量点(rm)。
8.根据权利要求6的装置,其中所述分析器(55)包括以下中的至少一个:
产生器(84),根据场信息(A[l+1];r0[l+1])产生展开矩阵(Ψd[l+1]);其中,展开矩阵(Ψd[l+1])利用与受验声源(Q0)的位置和方向对齐的坐标系或者利用声场的对称性;以及
估算器(83;81;85;91),通过利用包含与波方程的解相关的展开函数的所述展开矩阵(Ψd[l+1]),根据扫描数据集产生所述直达波系数
9.根据权利要求4的装置,
所述采集装置(1;14;77),测量在两个不重叠扫描范围(G1,G2)内分离的多个测量点(rm)处围绕所述受验声源(Q0)的声场的状态变量(pt(rm));其中,产生第一扫描数据集的第一扫描范围(G1)的测量点(rm)密度高于产生第二扫描数据集的第二扫描范围(G2)的测量点(rm)密度;
分析器(55),根据第一扫描数据集和第二扫描数据集产生与直达声(wrad)的展开相关的直达波系数
10.根据权利要求9的装置,其中所述分析器(55)包括:
第一子分析器(63),根据第二扫描数据集产生第二波系数其中第二扫描范围(G2)中的测量点的数量(M2)限制了第二波系数的最大阶数(N2);
内插器(65);根据第二波系数产生内插扫描数据集所述内插扫描数据集包含多个(M2S)元素,所述元素的数量(M2S)大于第二扫描范围(G2)内测量点的数量(M2);
加法器(67),根据第一扫描数据集和内插扫描数据集产生统一的数据 集
第二子分析器(69),根据统一的数据集产生所述直达波系数所述第二子分析器(69)以最大阶数(Nrad)提供至少一个波系数,所述波系数的最大阶数(Nrad)大于第二波系数的最大阶数(N2)。
11.根据权利要求4的装置,其中
所述直达波系数包含至少一个具有最大阶数(N)的系数(Cj,rad),其中
-所述直达波系数中的至少一个系数(Ci,rad=0)接近等于0,其中所述系数(Ci,rad=0)的阶数(Ni)小于所述最大阶数(N),或者
-扫描数据集中测量点(rm)的数量(M)小于所述直达波系数中系数(Cj,rad)的最大数量((N+1)2)。
12.根据权利要求11的装置,其中所述分析器(55)包括以下中的至少一个:
估算器(81),利用缩减展开矩阵(Ψd,red[l])根据扫描数据集产生缩减波系数
反向转换器(85),利用预先设定的选择矩阵(S[l]),根据缩减波系数产生迭代波系数其中所述迭代波系数中的元素的数量大于所述缩减波系数中的元素的数量;
估算器(87),根据迭代波系数产生贡献矢量(γ[l]),所述贡献矢量(γ[l])描述迭代波系数中的元素(Cj,rad)对受验声源(Q0)产生的总声波能量的贡献率(γj[l]);
选择器(89),根据贡献矢量(γ[l])产生所述预先设定的选择矩阵(S[l]),并且产生展开函数的最大阶数(N[l+1]);其中,如果其贡献率(γj[l])低于临界阈值(γ0),设置至少一个元素(Cj,rad)为0;
产生器(84),利用以下中的至少一个产生完全展开矩阵(Ψd[l+1]):
-展开函数的最大阶数N[l+1],
-所述扫描范围(Gm)内的测量点(rm)的数量和位置(R[l]),
-表示受验声源(Q0)的声学中心的位置的展开点(r0[l]),
-表示受验声源(Q0)的辐射的主方向的方向信息(A[l]),
-受验声源(Q0)产生的声场的旋转对称性(A[l]),
-受验声源(Q0)产生的声场的反射对称性(A[l]);
转换器(83),根据完全展开矩阵(Ψd[l+1])产生所述缩减展开矩阵(Ψd,red[l+1])的更新值,其中缩减展开矩阵(Ψd,red[l+1])中元素的数量小于完全展开矩阵(Ψd[l+1])中元素的数量;以及
控制器(91),如果在迭代中两个迭代步骤之间的迭代波系数的变化低于预先设定的阈值,其根据迭代波系数产生所述直达波系数
13.根据权利要求4的装置,
其中所述采集装置,包括
-至少一个扫描传感器(1),产生扫描输出(pt(rm)),所述扫描输出(pt(rm))表示扫描传感器(1)当前位置(rm)处被测量的声场的状态变量,
-定位装置(49;47;11;51),将扫描传感器(1)移动到扫描范围(Gm)内的至少一个测量点(rm),
-环境噪声传感器(76),产生表示扫描范围(Gm)外的环境位置(ra)处的声场的状态变量的环境噪声输出(pt(ra)),其中环境噪声传感器(76)与受验声源(Q0)的距离比扫描传感器(1)与受验声源(Q0)的距离远;
噪声识别器(77),根据扫描输出(pt(rm))和环境噪声输出(pt(ra)),产生有效的扫描数据集其中所述噪声识别器(77)通过分析环境噪声输出(pt(ra))检测扫描输出(pt(rm))中被所述环境噪声源(Q1)破坏的无效部分,并且从扫描数据集中排除检测到的无效部分。
14.根据权利要求1的装置,其中
如果扫描输出(pt(rm))包含被所述环境噪声源(Q1)破坏的无效部分,则所述噪声识别器(77)重复测量在测量点(rm)处声场的状态变量,或者
所述噪声识别器(77)包含存储装置,存储多次测量后的扫描输出(pt(rm))的有效部分,并且产生分配给测量点(rm)处包含有效信息和完全信息的扫描数据集的合并输出。
15.根据权利要求4的装置,其中
所述分析器(55),根据扫描数据集产生表示扫描数据集展开中的误差(ew)的声波误差系数(Erad);以及
外推器(11),根据声波误差系数(Erad)产生描述观测点(r)处的状态变量(prad(r)) 的误差(erad(r))的信息。
16.一种用于确定由受验声源(Q0)产生的直达声(wrad)的方法,包括:
测量围绕受验声源(Q0)的声场的状态变量(pt(rm));
通过收集在多个测量点(rm)处测量到的状态变量(pt(rm)),产生扫描数据集其中所述测量点设置在扫描范围(Gm)内,并且所述扫描数据集描述了直达声(wrad)和至少一个以下声波分量叠加:
-由环境噪声源(Q1)或由在外部边界(SB)上反射的直达声(wrad)所产生的输入声波(win),
-由所述输入声波(win)与受验声源(Q0)表面(S0)的相互作用所产生的第二声波(wsec);
通过对所述扫描数据集的滤波,产生滤波扫描数据集其中保留混响部分并削弱直达声部分;
产生与所述滤波扫描数据集中的输入波(win)的展开相关的后输入波系数其中后输入波系数确定了各展开函数(Ψ)的权重,展开函数(Ψ)是波方程的解;
产生与滤波扫描数据集中的第二声波(wsec)的展开相关的后输出波系数其中后输出波系数确定了各展开函数(Ψ)的权重,展开函数是波方程的解;
产生与所述扫描数据集中的总输入波(win)的展开相关的总输入波系数其中所述总输入波系数确定了各展开函数(Ψ)的权重,展开函数是波方程的解;
根据扫描数据集产生总输出波系数其中总输出波系数 与由第二声波(wsec)和直达声(wrad)构成的总输出波的展开相关,其中所述输出波系数确定了各展开函数(Ψ)的权重,展开函数是波方程的解;
根据总输入波系数和总输出波系数以及后输入波系数和后输出波系数产生直达波系数其中所述直达波系数与直达声(wrad)的展开相关并且确定了各展开函数(Ψ)的权重,展开函数是波方程的 解;
根据直达波系数产生预测状态变量(prad(r)),所述预测状态变量描述了扫描范围(Gm)内边界(Sin)外侧的预先设定的观测点(r)处的直达声场(wrad)。
17.根据权利要求16的方法,进一步包括以下中的至少一个:
根据后输入波系数和后输出波系数产生透过参数(kt),其中透过参数(kt)描述被扫描范围(Gm)封闭的空间的声学透过性;
根据总输入波系数和透过参数(kt)产生第二波系数其中第二波系数与受验声源(Q0)产生的第二声波(wsec)的展开相关;
根据总输出波系数和第二波系数产生所述直达波系数
18.根据权利要求16的方法,进一步包括以下中的至少一个:
通过使总输入波系数和后输入波系数互相关产生透射参数(kx);其中透射参数(kx)描述总输入声波(win)的指向性与输入声波(win)混响部分的指向性的一致性;
根据总输入波系数和透射参数(kx),合成第二波系数其中第二波系数与受验声源(Q0)产生的第二声波(wsec)的展开相关;
根据总输出波系数和第二波系数产生所述直达波系数
19.一种用于确定由受验声源(Q0)产生的直达声(wrad)的方法,包括:
测量围绕受验声源(Q0)的声场的状态变量(pt(rm));其中所述状态变量描述直达声(wrad)和至少一个以下声波分量的叠加:
-由环境噪声源(Q1)或在外部边界(SB)上反射的直达声(wrad)所产生的输入声波(win),
-由所述输入声波(win)与受验声源(Q0)表面(S0)的相互作用所产生的第二声波(wsec);
通过收集设置在扫描范围(Gm)内的多个测量点(rm)处测量到的状态变量(pt(rm)),产生扫描数据集
通过对扫描数据集进行滤波产生滤波扫描数据集其中所述滤波保存直达声(wrad)并且削弱其它声波分量(win,wsec),其它声波分量(win,wsec)与直达声(wrad)相比拥有更大的时间延迟;
产生与滤波扫描数据集的展开相关的滤波系数其中滤波系数确定了各展开函数(Ψ)的权重,展开函数是波方程的解;所述展开方式是对信号频率高于设定的截止频率(fg)的直达声(wrad)的有效表示;
利用表示输入声波(win)和输出声波(wout)的分离展开函数,将信号频率低于截止频率(fg)的扫描数据集进行展开;
产生与输出声波(wout)的展开相关的分离直达波系数其中分离直达波系数确定了各展开函数(Ψ)的权重,展开函数是波方程的解;所述展开方式是对信号频率低于设定的截止频率(fg)的直达声(wrad)的有效表示;
通过以下方式,产生与直达声(wrad)的展开相关的直达波系数
-为信号频率低于设定交叉频率(fC)的直达波系数分配分离直达波系数
-为信号频率高于或等于所述交叉频率(fC)的直达波系数分配滤波系数
以及根据直达波系数产生预测状态变量(prad(r)),所述预测状态变量(prad(r))描述扫描范围(Gm)内边界(Sin)外侧的预先设定的观测点(r)处产生的直达声(wrad)。
20.根据权利要求19的方法,进一步包括:
根据与所述展开函数相关的所述分离直达波系数通过评估扫描数据集 与展开数据集之间的不匹配,产生第一误差(eio);
通过用于输入的滤波声波(win,w)与输出的滤波声波(wout,w)的分离展开函数,对滤波扫描数据集进行展开;
通过比较输入的滤波声波(win,w)与输出的滤波声波(wout,w)之间的能量,产生滤波系数的第二误差(ew);
将第一误差(eio)和第二误差(ew)汇总为总误差;
产生给出最小总误差的所述截止频率(fg)的最优值;以及
校正交叉频率(fC)到截止频率(fg)。
21.根据权利要求19的方法,还包括:
根据扫描数据集产生场信息(A[l+1];r0[l+1]),所述场信息包含以下中的至少一个:
-与受验声源(Q0)的声学中心的位置相关的展开点(r0[l+1]),
-与受验声源(Q0)的主辐射方向相关的方向信息(A[l+1]),
-由受验声源(Q0)产生的声场的旋转对称信息(A[l+1]),
-由受验声源(Q0)产生的声场的反射对称信息(A[l+1]);
根据描述受验声源(Q0)产生的声场的性能的场信息(A[l+1];r0[l+1]),产生与直达声(wrad)的展开相关的所述直达波系数
22.根据权利要求21的方法,其中所述的产生扫描数据集包含以下中的至少一个:
测量扫描传感器(1)的当前位置(rm)处的声场的状态变量(pt(rm));
根据场信息(A[l+1];r0[l+1])产生扫描矢量(R[l+1]);其中扫描矢量(R[l+1])包括附加测量点(rm)的位置,其给出了有关受验声源(Q0)产生的直达声的特定信息;
将扫描传感器(1)移动到由扫描矢量(R[l+1])限定的附加测量点(rm);
测量附加测量点(rm)处的声场的状态变量(pt(rm))。
23.根据权利要求21的方法,其中产生直达波系数包括
根据场信息(A[l+1];r0[l+1])产生展开矩阵(Ψd[l+1]);其中,展开矩阵(Ψd[l+1])通过在展开中利用与受验声源(Q0)的位置和方向对齐的坐标系或者利用声场的对称性而包括最少的元素;以及
通过利用所述展开矩阵(Ψd[l+1]),根据扫描数据集产生所述直达波系数
24.根据权利要求19的方法,还包括:
测量分离的两个不重叠扫描范围(G1,G2)内的多个测量点(rm)处围绕所述受验声源(Q0)的声场的状态变量(pt(rm));其中,第一扫描范围的测量点(rm)密度高于第二扫描范围(G2)的测量点(rm)密度;
产生表示第一扫描范围(G1)的第一扫描数据集
产生表示第二扫描范围(G2)的第二扫描数据集
根据第一扫描数据集和第二扫描数据集产生与直达声(wrad)的展开相关的直达波系数
25.根据权利要求24的方法,其中产生直达波系数包括:
根据第二扫描数据集产生第二波系数其中第二扫描范围(G2)测量点的数量(M2)限制了第二波系数的系数的最大阶数(N2);
根据第二波系数产生内插扫描数据集其中内插扫描数据集 包含多个(M2S)元素,所述元素的数量(M2S)大于第二扫描范围(G2)内测量点的数量(M2);
根据第一扫描数据集和内插扫描数据集产生统一的数据集以及
根据统一的数据集产生所述直达波系数其中所述直达波系数包含至少一个具有最大阶数(Nrad)的系数,所述系数的最大阶数(Nrad)大于第二波系数 中的最大阶数(N2)。
26.根据权利要求19的方法,还包括:
根据扫描数据集产生与直达声(wrad)的展开相关的直达波系数所述直达波系数包含最大阶数(N),其中;
-所述直达波系数中的至少一个系数(Ci,rad=0)接近等于0,其中所述系数(Ci,rad=0)的阶数(Ni)小于所述最大阶数(N),或者
-扫描数据集中测量点(rm)的数量(M)小于所述直达波系数中系数(Cj,rad)的最大数量((N+1)2)。
27.根据权利要求26的方法,其中所述的产生直达波系数包括以下中的至少一个:
通过利用缩减展开矩阵(Ψd,red[l]),根据扫描数据集产生缩减波系数
通过利用预先设定的选择矩阵(S[l]),根据缩减波系数产生迭代波系数其中所述迭代波系数中的元素的数量大于所述缩减波系数中的元素的数量;
根据迭代波系数产生贡献矢量(γ[l]),所述贡献矢量(γ[l])描述系数(Cj,rad)对受验声源(Q0)产生的总声波能量的贡献率(γj[l]);
根据贡献矢量(γ[l])产生所述选择矩阵(S[l]);其中,如果其贡献率(γj[l])低于临界阈值(γ0),选择矩阵(S[l])设置系数(Cj,rad)为0;
根据贡献矢量(γ[l])产生与展开相关的系数的最大阶数(N[l+1];
利用以下中的至少一个产生完全展开矩阵(Ψd[l+1]):
-展开的最大阶数N[l+1],
-所述扫描范围(Gm)内的测量点(rm)的数量和位置(R[l]),
-表示受验声源(Q0)的声学中心的位置的展开点(r0[l]),
-表示受验声源(Q0)的辐射的主方向的方向信息(A[l]),
-由受验声源(Q0)产生的声场的旋转对称性信息(A[l]),
-由受验声源(Q0)产生的声场的反射对称性信息(A[l]);
根据完全展开矩阵(Ψd[l+1])产生所述缩减展开矩阵(Ψd,red[l+1])的更新值,其中缩减展开矩阵(Ψd,red[l+1])中元素的数量小于完全展开矩阵(Ψd[l+1])中元素的数量;
如果在迭代中两个迭代步骤的迭代波系数值之间的差低于预先设定的阈值,根据迭代波系数产生所述直达波系数以及
根据直达波系数产生预测状态变量(prad(r)),所述预测状态变量描述扫描范围(Gm)的内边界(Sin)外侧预先设定的观测点(r)处产生的直达声(wrad)。
28.根据权利要求19的方法,还包括:
将扫描传感器(1)移动到位于扫描范围(Gm)内的至少一个测量点(rm);
产生表示扫描传感器(1)的当前位置(rm)处被测量的声场的状态变量的扫描输出(pt(rm));
产生表示扫描范围(Gm)外的环境位置(ra)处的声场的状态变量的环境噪声输出(pt(ra)),其中环境噪声传感器(76)与受验声源(Q0)之间的距离比扫描传感器(1)与受验声源(Q0)之间的距离远;
分析环境噪声输出(pt(ra))和扫描输出(pt(rm));
将由所述环境噪声源(Q1)损坏的扫描输出(pt(rm))的无效部分与有效部分分离;
收集扫描数据集中扫描输出(pt(rm))的有效部分;
根据扫描数据集产生与直达声(wrad)的展开相关的直达波系数
29.根据权利要求28的方法,其中所述的收集扫描输出(pt(rm))包括:
如果扫描输出(pt(rm))包含被所述环境噪声源(Q1)破坏的至少一个无效部分,重复测量在测量点(rm)处声场的状态变量;
存储扫描输出(pt(rm))的有效部分;
产生合并多个测量值有效部分的合并输出;以及
如果合并输出包括测量点(rm)处的有效信息和完全信息,分配该合并输出到扫描数据集
30.根据权利要求19的方法,还包括:
产生表示扫描数据集展开中的误差(ew)的声波误差系数(Erad);
根据直达波系数产生预测状态变量(prad(r)),所述预测状态变量描述扫描范围(Gm)的内边界(Sin)外侧预先设定的观测点(r)处产生的直达声(wrad);以及
根据声波误差系数(Erad)产生描述观测点(r)处的预测数量(prad(r))的误差(erad(r))的信息。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DE201310000684 DE102013000684B3 (de) | 2013-01-11 | 2013-01-11 | Anordnung und Verfahren zur holografischen Bestimmung des Direktschalles akustischer Quellen |
DE102013000684.8 | 2013-01-11 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103929706A CN103929706A (zh) | 2014-07-16 |
CN103929706B true CN103929706B (zh) | 2017-05-31 |
Family
ID=49880092
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410079512.1A Active CN103929706B (zh) | 2013-01-11 | 2014-01-10 | 用于测量声源产生的直达声的装置及方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US9584939B2 (zh) |
CN (1) | CN103929706B (zh) |
DE (1) | DE102013000684B3 (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104469647B (zh) * | 2014-12-18 | 2017-10-20 | 中国电子科技集团公司第三研究所 | 一种宽带阵列流型的测算系统及其测算方法 |
CN104568137B (zh) * | 2015-01-12 | 2018-04-27 | 国家电网公司 | 噪声测试中特定噪声源贡献率的判定方法 |
JP6631713B2 (ja) * | 2016-07-22 | 2020-01-15 | ヤマハ株式会社 | タイミング予想方法、タイミング予想装置、及び、プログラム |
CN107785025B (zh) * | 2016-08-25 | 2021-06-22 | 上海英波声学工程技术股份有限公司 | 基于房间脉冲响应重复测量的噪声去除方法及装置 |
CN109375197B (zh) * | 2018-10-23 | 2023-02-10 | 哈尔滨工程大学 | 一种小尺寸矢量阵低频散射校正方法 |
CN109644456B (zh) * | 2018-11-29 | 2023-10-10 | 北京小米移动软件有限公司 | 波束扫描范围确定方法、装置、设备及存储介质 |
CN111157952B (zh) * | 2020-01-06 | 2022-04-01 | 北京大学 | 一种基于移动麦克风阵列的房间边界估计方法 |
CN111561991B (zh) * | 2020-06-08 | 2021-03-30 | 中国海洋大学 | 基于边缘填充和傅里叶变换的近场声全息方法 |
CN114034379B (zh) * | 2021-11-08 | 2023-04-14 | 北京理工大学 | 一种基于直达声场的封闭空腔噪声试验平台搭建方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6738479B1 (en) * | 2000-11-13 | 2004-05-18 | Creative Technology Ltd. | Method of audio signal processing for a loudspeaker located close to an ear |
JP2007047539A (ja) * | 2005-08-11 | 2007-02-22 | Sony Corp | 音場補正装置及び音場補正方法 |
CN101137250A (zh) * | 2007-10-08 | 2008-03-05 | 卢灿光 | 一种车载音响自适应音场调节方法及装置 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7599248B2 (en) * | 2006-12-18 | 2009-10-06 | The United States Of America As Represented By The Secretary Of The Navy | Method and apparatus for determining vector acoustic intensity |
KR100838239B1 (ko) * | 2007-04-17 | 2008-06-17 | (주)에스엠인스트루먼트 | 음질 표시 장치, 음질 표시 방법, 음질 표시 프로그램을 기록한 컴퓨터로 읽을 수 있는 매체 |
JP2010011433A (ja) * | 2008-05-30 | 2010-01-14 | Nittobo Acoustic Engineering Co Ltd | 音源分離及び表示方法並びにシステム |
US8077540B2 (en) * | 2008-06-13 | 2011-12-13 | The United States Of America As Represented By The Secretary Of The Navy | System and method for determining vector acoustic intensity external to a spherical array of transducers and an acoustically reflective spherical surface |
WO2010003837A1 (en) * | 2008-07-08 | 2010-01-14 | Brüel & Kjær Sound & Vibration Measurement A/S | Reconstructing an acoustic field |
-
2013
- 2013-01-11 DE DE201310000684 patent/DE102013000684B3/de active Active
-
2014
- 2014-01-10 US US14/152,556 patent/US9584939B2/en active Active
- 2014-01-10 CN CN201410079512.1A patent/CN103929706B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6738479B1 (en) * | 2000-11-13 | 2004-05-18 | Creative Technology Ltd. | Method of audio signal processing for a loudspeaker located close to an ear |
JP2007047539A (ja) * | 2005-08-11 | 2007-02-22 | Sony Corp | 音場補正装置及び音場補正方法 |
CN101137250A (zh) * | 2007-10-08 | 2008-03-05 | 卢灿光 | 一种车载音响自适应音场调节方法及装置 |
Non-Patent Citations (1)
Title |
---|
"Recovery of the free field using the spherical wave superposition method";C-X.Bi et al;《ACOUSTICS 2012》;20120427;全文 * |
Also Published As
Publication number | Publication date |
---|---|
DE102013000684B3 (de) | 2014-01-23 |
US20140198921A1 (en) | 2014-07-17 |
CN103929706A (zh) | 2014-07-16 |
US9584939B2 (en) | 2017-02-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103929706B (zh) | 用于测量声源产生的直达声的装置及方法 | |
CN102089633B (zh) | 用于重建声学场的方法 | |
Sarradj | Three‐Dimensional Acoustic Source Mapping with Different Beamforming Steering Vector Formulations | |
JP4891910B2 (ja) | ノイズ源の遠距離音場分析 | |
Landschoot et al. | Model-based Bayesian direction of arrival analysis for sound sources using a spherical microphone array | |
Fernandez-Grande et al. | Generative models for sound field reconstruction | |
Bucaro et al. | Acoustic identification of buried underwater unexploded ordnance using a numerically trained classifier (L) | |
Müller-Trapet et al. | On the in situ impedance measurement with pu-probes—Simulation of the measurement setup | |
Park et al. | Noise localization method for model tests in a large cavitation tunnel using a hydrophone array | |
Javier et al. | Analysis of the underwater radiated noise generated by hull vibrations of the ships | |
Lim et al. | Potential-based methodology for active sound control in three dimensional settings | |
Du et al. | Modeling of multiply connected sound spaces by the surface coupling approach | |
WO2005071576A1 (en) | Reconstruction of transient acoustic radiation from a finite object subject to arbitrarily time-dependent excitation | |
Caiti et al. | Mapping ocean sediments by RBF networks | |
Liu et al. | Sound field reconstruction using multipole equivalent source model with un-fixed source locations | |
Williams et al. | Experimental estimation of in vacuo structural admittance using random sources in a non-anechoic room | |
Sanalatii et al. | Estimation of loudspeaker frequency response and directivity using the radiation-mode method | |
Zhang et al. | Separation of nonstationary sound fields based on the time-domain equivalent source method with single layer pressure-velocity measurements | |
Bastine et al. | Time-frequency-dependent directional analysis of room reflections using eigenbeam processing and von Mises–Fisher clustering | |
Sternini et al. | Acoustic digital twin for passive structural health monitoring | |
Flynn et al. | High-resolution acoustic localization of changes in spatially-distributed coherent sources for structural health monitoring | |
Feng et al. | A study on the asymmetric cylinder wall thickness difference discrimination by dolphins | |
Pezzoli | Space-time Parametric approach to Extended Audio Reality (SP-EAR) | |
Cheng et al. | Sound field separation technique for nonplanar sources based on field prediction by using statistically optimized near-field acoustic holography | |
CN108303699A (zh) | 超声相控阵远场超分辨率成像方法、装置、存储介质和系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |