CN102104816B - 基于群时延的低音管理 - Google Patents

基于群时延的低音管理 Download PDF

Info

Publication number
CN102104816B
CN102104816B CN201010601232.4A CN201010601232A CN102104816B CN 102104816 B CN102104816 B CN 102104816B CN 201010601232 A CN201010601232 A CN 201010601232A CN 102104816 B CN102104816 B CN 102104816B
Authority
CN
China
Prior art keywords
group delay
response
loud speaker
filter
frequency
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
Application number
CN201010601232.4A
Other languages
English (en)
Other versions
CN102104816A (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.)
Harman Becker Automotive Systems GmbH
Original Assignee
Harman Becker Automotive Systems GmbH
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 Harman Becker Automotive Systems GmbH filed Critical Harman Becker Automotive Systems GmbH
Publication of CN102104816A publication Critical patent/CN102104816A/zh
Application granted granted Critical
Publication of CN102104816B publication Critical patent/CN102104816B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • H04S7/301Automatic calibration of stereophonic sound system, e.g. with test microphone
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04RLOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
    • H04R2499/00Aspects covered by H04R or H04S not otherwise provided for in their subgroups
    • H04R2499/10General applications
    • H04R2499/13Acoustic transducers and sound field adaptation in vehicles
    • 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/307Frequency adjustment, e.g. tone control

Landscapes

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

Abstract

基于群时延的低音管理。本发明提供一种用于改善收听空间中低音频率范围内的音频再现的全通滤波器设计方法。该收听空间包含至少一个扬声器和至少一个收听位置。该方法包含:针对每个扬声器,提供需要均衡且与收听空间中一个预定位置相关的群时延响应;针对每一个布置在一个相应扬声器上游的全通滤波器计算滤波器系数;该全通滤波器具有使得相应群时延响应匹配预定目标群时延响应的传递特性。

Description

基于群时延的低音管理
技术领域
本发明涉及一种用于自动均衡音频系统产生的低音频率(低音)范围中的群时延的方法和系统,本发明还涉及“低音管理”方法或系统。
背景技术
直到现在,例如在机动车辆中手动地听觉优化专用系统是常见的实践。尽管已经付出很大的努力来使得该手动过程自动化,这些方法和系统在实践中显示出不足或者是极其复杂和昂贵的。在诸如汽车内部这样的小的高反射区域,仅仅实现了听觉中的微小改进。在一些情况中,结果甚至更糟。
尤其在约100至150赫兹以下的频率范围中,小的高反射空间内部的驻波可能在诸如机动车辆中的两个前排座和两个后排乘客座位这样的各个收听位置导致极不相同的声压级(SPL)。这些不同的声压级使得个人的音频感觉与他/她的收听位置相关。然而,专业声学家的工作已经证明这一事实:即使使用简单的方式也能实现好的听觉结果。
已知一种允许几乎在任意空间建模任意音响系统的方法。然而,这种所谓的波前合成需要诸如计算功率、存储器、扬声器、放大器通道等大量资源。由于成本和可用性等原因,该技术因而不适于很多应用,尤其不适于用在汽车产业中。
而且,自动低音管理系统是已知的,这种自动低音管理系统的目标在于均衡且同时最大化收听空间内听者位置处的低音频率范围中的声压级。然而,在听力测试中,结果评估为不足,指示执行SPL均衡可能仅仅是改进低音频率级中的声音再现的质量的一个步骤。
需要一种自动低音管理系统,其能够适当地代替原先使用的通过有经验的声学家实施的手工均衡的复杂过程,且可靠地改善低音频率范围中的声印象。
发明内容
本发明公开了一种用于改善收听空间中低音频率范围内音频再现的方法。收听空间包括至少一个扬声器和至少一个收听位置。该方法包括:为每个扬声器提供与收听空间内的一个预定位置相关的需要均衡的群时延响应;为每个布置在一个相应扬声器上游的(多个)全通滤波器计算滤波器系数,该(多个)全通滤波器具有使得(多个)相应群时延响应匹配预定目标群时延响应的传递特性。
附图说明
参考下面的附图和说明能够更好地理解本发明。在附图中,相似的参考标号表示相应的部件。附图中:
图1是说明在汽车车厢中4个不同收听位置测量的频率上以分贝表示的声压级的视图,该汽车具有供应到扬声器的未经修改的音频信号。
图2是说明汽车车厢中的驻声波的示意性侧视图,该驻声波是造成收听位置之间的声压级(SPL)中大差异的原因。
图3是说明机动车辆车厢中的座位布置以及收听位置以及扬声器布置的示意性顶视图;
图4说明限定目标全通滤波器的群时延的频率相关限制的频率上的群时延约束函数的示例;以及
图5是说明扬声器上游的音频通道中的群时延均衡滤波器的布置的示意性顶视图。
具体实施方式
在借助汽车中的扬声器或一组扬声器再现音频信号时,即使扬声器贯穿汽车对称布置,对于不同收听位置存在的声压级(SPL),汽车车厢中的测量得出明显不同的结果。图1的视图说明了这种结果。在视图中,示出了4条曲线,每条曲线说明在向扬声器提供音频信号的同时,车厢中的4个不同收听位置(即,两个前排和两个后排座位的头枕附近)测量的频率上以分贝(dB)表示的声压级。可以看出,依赖于应用的频率,在空间前排的收听位置测量的声压级和在后排的收听位置测量的声压级的差异高达15dB。然而,典型地可以在作为低音频率范围的一部分的约40至90赫兹的频率范围内观察到SPL曲线之间的最大差距。
“低音频率范围”并不是很好定义的术语,但是广泛地用在例如从0至80赫兹、0至100赫兹或甚至0至150赫兹的范围中的低频的音响系统中。尤其是,当使用具有布置在后窗架中或后备箱中的低音炮的汽车音响系统时,可以观察到收听空间内的声压级的不当分配。60至70赫兹(参考图1)之间的SPL最大值可能倾向于认为是嗡嗡声且令后排乘客不悦。
可以观察的汽车的前排和后排的收听位置中的不同位置中声压级的大差异的频率范围与收听空间的尺寸相关。原因可以参考图2解释,图2是汽车的示意性侧视图。半波长(表示为λ/2)与车厢的纵向一致。当假设c=340m/s的声速时,λ/2=2.5m的典型长度产生f=c/λ=68Hz的频率。从图1可以看出,大约在该频率,在后排收听位置可以观察到最大SPL。这指示在汽车(收听室)内部的纵向和横向中若干驻波的交叠可能是造成收听空间中不均匀SPL分布的原因。
例如从公报EP2051543A1和2043384A1获知自动低音管理系统。这种系统的目标是均衡和(可选地)同时最大化在收听空间的听者位置处低音频率范围的声压级。然而,在听力测试中,所得的低音再现被评估为不足(即,淘汰或疲软),这指示执行SPL均衡可能只是改善低音频率级中声音再现质量的一个步骤。此处描述的新颖的低音管理系统考虑低音频率范围中再现的音频信号的群时延。
图3说明贯穿小且回响的收听空间(诸如机动车辆车厢),收听位置FR、FL、RR、RL和扬声器的简单布置。然而,本发明不应限制于汽车应用且可应用于任意收听空间。而且,本领域技术人员将理解,本实例可以容易地调适为考虑多于或少于4个的收听位置。
图3中示出的4个收听位置FL、FR、RL、RR代表机动车辆车厢中的左前、右前、左后、右后收听位置。在本示例中,贯穿车厢布置5个扬声器LS1至LS5,诸如左前扬声器LS1、右前扬声器LS2、左后扬声器LS3、右后扬声器LS4以及中后扬声器LS5(例如低音炮)。当向扬声器LS1至LS5应用不同频率的测试信号(或宽带测试信号)时,在每个收听位置可以观察到所得的脉冲响应h[k]、频率响应H(ω)(即,幅度|H(ω)|和相位的传递函数)以及群时延(τG(ω))响应。这种“系统识别”方法在声学领域是公知的。频率响应是脉冲响
H(ω)=PFT{h[k]},(1)
应的傅立叶变换,且通过快速傅立叶变换(FFT)粗略估计:
而且,群时延限定为:
在每个收听位置FL、FR、RL、RR观察的频率响应Hx(ω)(X∈{FL、FR、RL、RR})是从每单个扬声器LS1至LS5得到的频率响应的叠加,即:
HX(ω)=Sum{HX-LSi(ω)},fori=1,...,5,(3)
其中HX-LSi(ω)是描述在收听位置X观察的声学信号与供应到扬声器LSi且从扬声器LSi(参见图3)辐射的相应音频信号之间的关系的系统传递函数。类似地,在收听位置X观察的群时延响应τGX(ω)被认为是成分τGX-LSi(ω)的叠加,其中i=1,...,5且X∈{FL、FR、RL、RR},在本示例中:
τGX(ω)=Sum{τGX-LSi(ω)},fori=1,...,5,(4)
从心理声学研究(例如见J.Blauert,P.Laws:Perceptibilityofgroupdelaydistortions,in:J.Acoust.Soc.Am.,Vol.63,No.5,1978)可以知道:超过给定频率相关阈值的群时延失真可以被人类听者感觉。因而,通过减小群时延失真,即,通过均衡低音频率范围中的群时延响应,可以改善高保真度的音频再现的质量。
可以采用供应扬声器LS1、LS2、...LS5的音频通道中的相位滤波器(全通滤波器HAP1、HAP2、...HAP5,见图5)来均衡在收听空间所需位置的群时延响应。这种所需位置可以是收听位置,或者,为了计算多于一个收听位置,可以是两个或更多收听位置之间的位置。类似地,如果需要改善多于一个收听位置的声音印象,平均群时延响应可以被均衡,该平均群时延响应可以通过在4个收听位置FL、FR、RL、RR观察的4个群时延响应的平均所表达。
为了进一步讨论,经历均衡的群时延响应一般表示为τG(ω),相应的传递函数(频率响应)是H(ω)。如上所述,群时延响应τG(ω)可以是在收听空间中给定位置观察的群时延响应,或者是从各个(已知)收听位置观察的两个或更多群时延响应计算的平均群时延响应。
如等式4所示,考虑的群时延响应τG(ω)可以分解成很多被加数:
τG(ω)=τG1(ω)+τG2(ω)+...+τGN(ω),(5)
其中被加数的数目等于收听室中布置的扬声器的数目N,每个被加数τGi(ω)对应于定义的扬声器LSi。对于相应的相位,可以完成相同的分解:
其中相位响应是复数传递函数H(ω)的相位,即应当注意,相位被加数以及群时延被加数τGi(ω)可以容易地从测量的脉冲响应得出,该脉冲响应限定从每个扬声器到每个考虑的收听位置的传递特性。仅仅是给出示例,经历均衡的群时延τG(ω)可以是在收听位置FL、FR、RL、RR其中每一个观察的群时延(τGFL(ω)、τGFR(ω)、τGRL(ω)、τGRR(ω))的平均;这些群时延τGX(ω)(X∈{FL、FR、RL、RR})其中每一个是与单个扬声器LS1、LS2、...LS5相关的群时延的和τGX-L51(ω)+τGX-L52(ω)+τGX-L53(ω)+τGX-L54(ω)+τGX-L56(ω)。类似地,等式6中的相位响应可以是在各个收听位置FL、FR、RL、RR且涉及扬声器LSi观察的相位响应 的平均值。
对于群时延均衡,布置在供应扬声器LSi的每个音频通道中的全通滤波器被设计为具有这种相位响应使得等式5中的每个所得的群时延响应τGi(ω)(i=1,2,...)匹配预定目标(即所需)群时延响应τTARGET(ω)。因而,具有相位响应的全通滤波器HAPi(ω)可以被认为是群时延均衡滤波器。目标群时延响应τTARGET(ω)直接与目标相位响应相关,且相应地,布置在扬声器LSi上游的音频通道中的全通滤波器的寻找的相位响应是:
其中N是扬声器的数目(在图3的示例中,N=5)。全通滤波器的幅度响应|HAPi(ω)|当然是|HAPi(ω)|=1。对于本领域技术人员而言,存在很多可能性来从等式7的相位响应计算相应的全通脉冲响应(即FIR滤波器系数)hAPi[k]。下面给出一个示例。
复数全通传递函数的实部和虚部设置如下:
复数全通传递函数HAPi(ω)因而可以写成:
其中j是-1的方根。用于低音频率范围以上频率的相位值被设置为0,以避免低音频率范围之外的宽带相位失真,即:
等式10的传递函数HAPi(ω)可以借助反FFT变换到(离散)时域。在变换到时域之前,必须确保是对称的,即:
real{HAPi(ω)}=real{HAPi(-ω)}and(12)
imag{HAPi(ω)}=-imag{HAPi(-ω)},(13)
以获得实值脉冲响应hAPi[k]。一般而言,所得的全通滤波器脉冲响应hAPi[k]将是非因果联系的。为了获得具有有限脉冲响应的因果滤波器,脉冲响应hAPi[k]在时域中设计时必须时移且被截取。备选地,传递函数HAPi(ω)可以与窗口函数相乘以本质上实现相同的结果(还参见Oppenheim,Schafer:“DesignofFIRFiltersbyWindowing”,in:Discrete-TimeSignalProcessing.2ndEd.,Section7.2,PrenticeHall,1999)。
然而,声音测试得出:使用如上所述经典FIR滤波器设计方案设计的全通滤波器(即,相位均衡滤波器)并没有带来音频质量的所需改善。不希望的音频失真恶化了高保真声音再现。这种失真是当使用标准设计方案设计时全通滤波器可以呈现的显著预振铃(pre-ring)的结果。已经发现,新颖的FIR全通滤波器设计方法可以解决上述问题,且显著增强音频再现的质量,特别是在低音频率范围。
根据本发明的一个示例,全通滤波器并不使用上述经典方法设计,而是使用如下所述的迭代优化方法设计。如果全通滤波器设计为使得根据限定(与频率相关)区间的群时延约束函数限制所得的群时延响应,则是有利的。即,所得的全通滤波器(与每个扬声器LSi相关的一个全通滤波器HAPi)的群时延响应保留在以CL(ω)和CU(ω)表示的约束函数所限定的范围内。
所需的相位响应由等式(7)给出且指定为在迭代滤波器设计过程的开始,相应的全通滤波器HAPi(ω)被初始化,例如,初始化为HAPi(ω)=exp(0)=1。而且,解决了下面的最小化任务(用于最小化误差函数E):
考虑附加条件:
d(arg(HApi(jω)))/dω<CU(ω)foranyω,and(14a)
d(arg(HApi(jω)))/dω>CL(ω)foranyω.(14b)
可以使用任意常见最小搜索方法。在测试中,使用Matlab(tm)函数“fminsearch”提供的Nelder-Mead单纯形法,以用于发现最佳全通滤波器系数HAPiOPY(ω)。
应当注意,可以选择等式(14)中使用的用于计算最小化误差的范式||·||,从而得出二次误差,即:
||x(ω)||=x(ω1)2+x(ω2)2+...+x(ωk)2(15)
其中K是离散频率值ωk的数目且因而是FIR全通滤波器的长度,例如K=4096。
约束函数CU和CL的一个示例在图4中示出。一般而言,约束函数的形状(例如,对于群时延上限,图4中的虚线)可以被描述为指数衰减曲线,诸如:
cu(ω)=a·exp(ω/p)+b(16)
其中a、p和b是常量参数,参数b限定渐近线。图4中说明的FIR滤波器“体时延(bulkdelay)”对应于全通FIR滤波器的一半长度。在本示例中,全通滤波器长度K是4096个阶,且相应地,对于44.1kHz的采样频率,体时延是对应于46.44ms的2048个阶。在图4的示例中,限定群时延上限的约束函数c(ω)为:
c(ω)=3.39ms·exp(ω/(2π·820Hz))+46.44ms,(17)
应当注意,限定下限的约束函数CL(ω)与函数CU(ω)相对于代表体时延的水平线是对称的。
在图4中示出了整体系统结构。全通滤波器分别布置在扬声器LS1、LS2、LS3、LS4、LS5其中每一个的上游的每个音频通道(HAP1、HAP2、HAP3、HAP4和HAP5)中。为了简单目的,在说明中省略功率放大器,由此全通传递函数HAP1、HAP2、HAP3、HAP4和HAP5如上面所解释设计为均衡与一个或更多收听位置相关的给定群时延响应以匹配预定目标群时延响应(例如,恒定群时延)。附加的线性(或恒定)相位滤波器可以布置在每个音频通道中以用于全局均衡,以实现所需的声音印象。当然,这些滤波器可以与音频通道中已经存在的其他滤波器组合(即卷积)以用于其他目的。
下面总结了图5中示出的系统的一些重要方面以及相应的均衡方法。如上面所讨论,采用图4中示出的系统以用于改善收听室内低音频率范围的音频再现。收听室包含至少一个扬声器和至少一个收听位置。在本示例中,存在在机动车辆车厢中提供的4个收听位置FL、FR、RL、RR和5个扬声器LSi(i∈{1,2,3,4,5})。相对于收听室中的预定位置需要均衡的群时延τG1(ω)、τG2(ω)、τG3(ω)、τG4(ω)、τG5(ω)与每个扬声器LS1、LS2、LS3、LS4和LS5相关。该预定收听位置可以是收听室中的任意位置,例如,4个收听位置之间的中间位置(距离每个收听位置FL、FR、RL、RR等间距)。然而,预定收听位置还可以是“虚拟”收听位置,其中需要均衡的相关群时延(用于每个扬声器)是与实际收听位置FL、FR、RL、RR相关的群时延响应的平均。例如,对于扬声器LSi,需要均衡的群时延响应限定为:
其中τGX-LSi(ω)(X∈{FL、FR、RL、RR})代表与收听位置X和扬声器LSi相关的群时延响应。如上面所讨论,需要均衡的每个群时延响应τGi(ω)可以转换成相应的相位响应
在每个扬声器上游的音频通道中布置一个群时延均衡滤波器。每个滤波器是全通滤波器,其传递特性由其滤波器系数限定。每个滤波器的滤波器系数设置为使得所得的群时延响应τGi(ω)匹配预定目标群时延响应τGTarget(ω)。实际中,可以通过设置滤波器系数使得相位响应(对应于群时延响应τGi(ω))匹配代表上述目标群时延响应τGTarget(ω)的目标相位响应来执行这种均衡。
用于改善收听室中低音频率范围内的音频再现的方法包括针对每个扬声器LSi提供需要均衡的群时延响应τGi(ω)的步骤,其中每个群时延响应τGi(ω)与收听室内一个预定位置相关。如上面所解释,这种预定位置可以是收听室中的任意实际位置,并且当(多个)平均群时延τGi(ω)需要均衡时可以是“虚拟”收听位置。该方法还包括计算用于全通滤波器的滤波器系数HAPi(ω)的步骤。在每个扬声器LSi的上游的相应音频通道中布置一个滤波器。全通滤波器HAPi(ω)各自具有传递特性,以使得所得的(多个)群时延响应τGi(ω)匹配预定目标群时延响应τTarget(ω)。
如上所述,可以通过设置滤波器的相位响应使得所得的相位响应(对应于群时延响应τGi(ω))匹配预定目标相位响应(对应于目标群时延响应τGTarget(ω)),来执行这种均衡。提供需要均衡的群时延响应τGi(ω)的步骤还包括针对每一对收听位置和扬声器X-LSi(X∈{FL、FR、RL、RR},i∈{1,2,3,4,5})提供相位响应该相位响应代表从扬声器LSi到相应收听位置X的音频信号的相位传递特性。由此,每个相位响应代表相应群时延τGX-LSi(ω)。然后,依赖于(多个)群时延响应τGX-LSi(ω),可以为每个扬声器LSi提供均衡的群时延响应τGi(ω)。这可以包括如上所述的加权平均。
最后,计算滤波器系数的上述步骤可以包括提供代表目标群时延响应τGTarget(ω)的目标相位响应以及针对每个扬声器计算代表要被均衡的群时延响应的相位响应与目标相位响应之间的与频率相关的相位差 且最后,针对每个扬声器,使用作为滤波器设计中所需的滤波器相位响应的计算的(多个)相位差来计算全通滤波器系数。
所得的群时延均衡滤波器可以与预定全局均衡滤波器卷积以调节整体声音印象。预定全局均衡滤波器可以具有任意希望的幅度响应和恒定或线性相位响应。
尽管详细描述了本发明及其优点,应当理解,可以做出各种改变、替换和变更而不偏离所附权利要求限定的本发明的精神和范围。
然而,本申请的范围并不旨在限制为说明书中描述的处理、加工、制造、物质组成、装置、方法以及步骤。本领域技术人员将容易从本发明的公开意识到,根据本发明,可以使用当前存在或以后发展的执行与此处描述的相应实施例具有基本相同功能或实现基本相同结果的处理、加工、制造、物质组成、装置、方法或步骤。因此,所附权利要求旨在在其范围内包括这种处理、加工、制造、物质组成、装置、方法或步骤。

Claims (15)

1.一种用于改善收听空间中低音频率范围内的音频再现的全通滤波器设计方法,该收听空间包含至少一个扬声器(LSi)和至少一个收听位置(FL),所述方法包括:
针对每个扬声器(LSi),提供需要均衡且与收听空间中一个预定位置相关的群时延响应(τGi(ω));
针对每一个被布置在一个相应扬声器(LSi)上游的全通滤波器(HAPi(ω))计算滤波器系数,全通滤波器(HAPi(ω))具有使得相应群时延响应(τGi(ω))匹配预定目标群时延响应(τGTarget(ω))的传递特性,其中,计算滤波器系数的步骤包括针对与一个相应扬声器(LSi)相关的每个全通滤波器(HAPi(ω)):
提供与频率有关的群时延约束,该与频率有关的群时延约束限定约束全通滤波器(HAPi(ω))的群时延响应(τGi(ω))的有限范围;
迭代地计算更新的滤波器系数,使得在与群时延约束相符时误差范数最小,该误差范数代表各个全通滤波器(HAPi(ω))的群时延响应(τGi(ω))从相应目标群时延响应(τGTarget(ω))的偏离。
2.根据权利要求1所述的方法,其中所述与频率有关的群时延约束限定随着频率的增加而指数衰减的与频率有关的区间。
3.根据权利要求2所述的方法,其中所述区间关于对应于所述全通滤波器的半个滤波器长度的全通体时延对称布置,所述全通滤波器是FIR滤波器。
4.根据权利要求2或3所述的方法,其中所述区间随着频率的增加渐进地接近恒定区间。
5.根据权利要求4所述的方法,其中所述区间由上限cu(ω)=a·exp(ω/p)+b和下限cL(ω)=-a·exp(ω/p)+b限定,
其中ω是以rad/s表示的频率,b是代表全通体时延的常量参数,且a和p是描述区间的指数缩窄的常量参数。
6.根据权利要求1所述的方法,其中提供需要均衡的群时延响应(τGi(ω))的步骤包括:
针对每一对收听位置(FL)和扬声器(LSi),提供代表从扬声器(LSi)到相应收听位置(FL)的音频信号的相位传递特性的相位响应每个相位响应代表相应的群时延响应(τGFL-LSi(ω));
基于群时延响应(τGFL-LSi(ω)),针对每个扬声器(LSi)提供需要均衡的群时延响应(τGi(ω))。
7.根据权利要求1所述的方法,其中针对每个扬声器提供需要均衡的群时延响应的步骤还包括:
针对每个扬声器(LSi),计算所有考虑的收听位置上与考虑的扬声器(LSi)相关的相位响应的加权平均,所得的平均相位响应代表需要均衡的群时延响应。
8.根据权利要求1所述的方法,其中计算滤波器系数的步骤包括:
提供代表目标群时延响应(τGTarget(ω))的目标相位响应
针对每个扬声器,计算代表需要均衡的群时延响应的相位响应与目标相位响应之间的与频率有关的相位差
针对每个扬声器,使用计算的相位差作为所需的滤波器相位响应来计算全通滤波器系数。
9.根据权利要求1所述的方法,还包括:
将全通滤波器系数的每个计算序列与预定全局均衡滤波器的滤波器系数的序列进行卷积。
10.根据权利要求9所述的方法,其中所述预定全局均衡滤波器是具有预定幅度响应的线性相位或恒定相位滤波器。
11.一种用于改善包含至少一个扬声器(LSi)和至少一个收听位置(FL)的收听空间中的低音频率范围内的音频再现的系统,关于收听空间内预定位置需要均衡的群时延响应与每个扬声器关联,该系统包括:
群时延均衡滤波器,布置在每个扬声器上游,每个滤波器是其传递特性由其滤波器系数限定的全通滤波器,
其中每个全通滤波器的滤波器系数设置为使得所得的群时延响应(τGi(ω))匹配预定目标群时延响应(τGTarget(ω)),
其中,提供与频率有关的群时延约束,该与频率有关的群时延约束限定约束全通滤波器的群时延响应的有限范围;
其中,迭代地计算更新的滤波器系数,使得在与群时延约束相符时误差范数最小,该误差范数代表各个全通滤波器的群时延响应从相应目标群时延响应的偏离。
12.根据权利要求11所述的系统,其中滤波器系数具有由与频率有关的群时延约束限定的群时延响应,所述与频率有关的群时延约束限定随着频率的增加而指数衰减的与频率有关的区间。
13.根据权利要求11或12所述的系统,其中,针对每个扬声器,需要均衡的群时延响应对应于相应的相位响应,该相位响应是根据与每对收听位置(FL)和扬声器(LSi)相关的相位特性计算出来的。
14.根据权利要求13所述的系统,其中,针对每个扬声器,需要均衡的群时延响应对应于相应的相位响应,该相位响应是与每对收听位置(FL)和扬声器(LSi)相关的相位响应的加权平均。
15.一种用于改善收听空间中的低音频率范围内的音频再现的系统,该收听空间包括至少一个扬声器和连接在每个扬声器上游的至少一个全通FIR滤波器,该FIR滤波器是根据权利要求1至9的方法之一设计的。
CN201010601232.4A 2009-12-22 2010-12-22 基于群时延的低音管理 Active CN102104816B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
EP09180411.2 2009-12-22
EP09180411A EP2357846A1 (en) 2009-12-22 2009-12-22 Group-delay based bass management

Publications (2)

Publication Number Publication Date
CN102104816A CN102104816A (zh) 2011-06-22
CN102104816B true CN102104816B (zh) 2016-01-13

Family

ID=42133680

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201010601232.4A Active CN102104816B (zh) 2009-12-22 2010-12-22 基于群时延的低音管理

Country Status (3)

Country Link
US (1) US9191766B2 (zh)
EP (2) EP2357846A1 (zh)
CN (1) CN102104816B (zh)

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9661190B2 (en) * 2012-05-31 2017-05-23 Dolby Laboratories Licensing Corporation Low latency and low complexity phase shift network
US9467111B2 (en) * 2012-06-29 2016-10-11 Audyssey Laboratories Operator adjustable full-bandwidth audio spectral shifting control with a simple listener interface
DE102013105375A1 (de) * 2013-05-24 2014-11-27 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. Tonsignalerzeuger, Verfahren und Computerprogramm zum Bereitstellen eines Tonsignals
GB2541639B (en) 2015-06-15 2019-06-12 Meridian Audio Ltd Asymmetric stereophonic bass compensation
CN105262503B (zh) * 2015-07-16 2018-04-24 中国电子科技集团公司第四十一研究所 一种基于群时延校准的多径时延产生装置及方法
BR112018008699B1 (pt) * 2015-10-30 2022-03-03 Dirac Research Ab Redução da diferença de fase entre canais de áudio em posições espaciais múltiplas
US10075789B2 (en) * 2016-10-11 2018-09-11 Dts, Inc. Gain phase equalization (GPEQ) filter and tuning methods for asymmetric transaural audio reproduction
EP3509320A1 (en) * 2018-01-04 2019-07-10 Harman Becker Automotive Systems GmbH Low frequency sound field in a listening environment
WO2020052756A1 (de) * 2018-09-12 2020-03-19 Ask Industries Gmbh Verfahren zum betrieb einer kraftfahrzeugseitigen audioausgabeeinrichtung
CN109089203B (zh) * 2018-09-17 2020-10-02 中科上声(苏州)电子有限公司 汽车音响系统的多声道信号转换方法及汽车音响系统
BR112021011597A2 (pt) * 2018-12-21 2021-08-31 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. Sistema de reprodução/simulação de som, aparelho para determinar um ou mais parâmetros de processamento e métodos
CN111526455A (zh) * 2020-05-21 2020-08-11 菁音电子科技(上海)有限公司 车载音响的校正增强方法及系统
CH719150A1 (de) 2021-11-17 2023-05-31 Rocket Science Ag Verfahren zum Eliminieren von Raummoden und digitaler Signalprozessor sowie Lautsprecher dafür.
EP4322554A1 (en) * 2022-08-11 2024-02-14 Bang & Olufsen A/S Method and system for managing the low frequency content in a loudspeaker system
CN117676418B (zh) * 2023-12-06 2024-05-24 广州番禺职业技术学院 一种用于混合相位系统中的声场均衡方法及系统

Family Cites Families (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4771466A (en) * 1983-10-07 1988-09-13 Modafferi Acoustical Systems, Ltd. Multidriver loudspeaker apparatus with improved crossover filter circuits
JPH03216097A (ja) * 1990-01-22 1991-09-24 Sony Corp スピーカシステム
GB9026906D0 (en) * 1990-12-11 1991-01-30 B & W Loudspeakers Compensating filters
US5325436A (en) * 1993-06-30 1994-06-28 House Ear Institute Method of signal processing for maintaining directional hearing with hearing aids
US6760451B1 (en) * 1993-08-03 2004-07-06 Peter Graham Craven Compensating filters
US7277554B2 (en) * 2001-08-08 2007-10-02 Gn Resound North America Corporation Dynamic range compression using digital frequency warping
US7769183B2 (en) * 2002-06-21 2010-08-03 University Of Southern California System and method for automatic room acoustic correction in multi-channel audio environments
US7567675B2 (en) * 2002-06-21 2009-07-28 Audyssey Laboratories, Inc. System and method for automatic multiple listener room acoustic correction with low filter orders
US8280076B2 (en) * 2003-08-04 2012-10-02 Harman International Industries, Incorporated System and method for audio system configuration
US7720237B2 (en) * 2004-09-07 2010-05-18 Audyssey Laboratories, Inc. Phase equalization for multi-channel loudspeaker-room responses
US9008331B2 (en) * 2004-12-30 2015-04-14 Harman International Industries, Incorporated Equalization system to improve the quality of bass sounds within a listening area
US8355510B2 (en) * 2004-12-30 2013-01-15 Harman International Industries, Incorporated Reduced latency low frequency equalization system
JP2009530915A (ja) * 2006-03-15 2009-08-27 ドルビー・ラボラトリーズ・ライセンシング・コーポレーション 立体音響サウンドイメージ化
EP2051543B1 (en) * 2007-09-27 2011-07-27 Harman Becker Automotive Systems GmbH Automatic bass management
US8194885B2 (en) * 2008-03-20 2012-06-05 Dirac Research Ab Spatially robust audio precompensation
CA2735244C (en) * 2009-05-18 2015-10-27 Ryan J. Mihelich Efficiency optimized audio system

Also Published As

Publication number Publication date
EP2357846A1 (en) 2011-08-17
EP2357847B1 (en) 2016-08-10
US9191766B2 (en) 2015-11-17
EP2357847A2 (en) 2011-08-17
US20110150241A1 (en) 2011-06-23
EP2357847A3 (en) 2011-12-28
CN102104816A (zh) 2011-06-22

Similar Documents

Publication Publication Date Title
CN102104816B (zh) 基于群时延的低音管理
EP1843635B1 (en) Method for automatically equalizing a sound system
US9049533B2 (en) Audio system phase equalization
CN104185134B (zh) 听音室内个别声音区的产生
EP2806664B1 (en) Sound system for establishing a sound zone
US20090086995A1 (en) Automatic bass management
JPH10509565A (ja) 録音及び再生システム
JP2006191562A (ja) リスニング領域内での、バスサウンドの音質を改善するためのイコライゼーションシステム
US20180146290A1 (en) Individual delay compensation for personal sound zones
CN107079229A (zh) 个人多声道音频预补偿控制器设计
EP3446493A1 (en) An active monitoring headphone and a method for calibrating the same
AU2016240348B2 (en) Apparatus and method for processing stereo signals for reproduction in cars to achieve individual three-dimensional sound by frontal loudspeakers
EP1843636B1 (en) Method for automatically equalizing a sound system
CN116367076A (zh) 车辆内音频处理方法、设备及存储介质
CN109923877B (zh) 对立体声音频信号进行加权的装置和方法
CN117676418B (zh) 一种用于混合相位系统中的声场均衡方法及系统
Johansson et al. Sound field control using a limited number of loudspeakers
EP4207802B1 (en) Sound collection loudspeaker apparatus, method and program for the same
Lin et al. Partial update adaptive filtering based on head-related model for in-vehicle audio enhancement
Zhang et al. Hybrid multi-channel system in automotive audio
Mannerheim Visually adaptive virtual sound imaging using loudspeakers
Lakhdhar et al. Iterative equalization of room transfer function using biquadratic filters
Bai et al. Signal Processing Implementation and Comparison of Automotive Spatial Sound Rendering Strategies
Wallace et al. Developing an in-car 3D audio system using the latest Virtual Audio Methods

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant