CN114051636A - 具有经改进频率分辨率的低延迟音频滤波器组 - Google Patents

具有经改进频率分辨率的低延迟音频滤波器组 Download PDF

Info

Publication number
CN114051636A
CN114051636A CN202080045434.7A CN202080045434A CN114051636A CN 114051636 A CN114051636 A CN 114051636A CN 202080045434 A CN202080045434 A CN 202080045434A CN 114051636 A CN114051636 A CN 114051636A
Authority
CN
China
Prior art keywords
response
filter
ripple
generating
generate
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.)
Pending
Application number
CN202080045434.7A
Other languages
English (en)
Inventor
D·S·麦格拉思
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.)
Dolby Laboratories Licensing Corp
Original Assignee
Dolby Laboratories Licensing Corp
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 Dolby Laboratories Licensing Corp filed Critical Dolby Laboratories Licensing Corp
Publication of CN114051636A publication Critical patent/CN114051636A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10LSPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L19/00Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis
    • G10L19/04Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis using predictive techniques
    • G10L19/26Pre-filtering or post-filtering
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00Networks using digital techniques
    • H03H17/02Frequency selective networks
    • H03H17/0248Filters characterised by a particular frequency response or filtering method
    • H03H17/0264Filter sets with mutual related characteristics
    • H03H17/0266Filter banks
    • H03H17/0267Filter banks comprising non-recursive filters
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10LSPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L19/00Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis
    • G10L19/008Multichannel audio signal coding or decoding using interchannel correlation to reduce redundancy, e.g. joint-stereo, intensity-coding or matrixing

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Computer Hardware Design (AREA)
  • Health & Medical Sciences (AREA)
  • Signal Processing (AREA)
  • Computational Linguistics (AREA)
  • Audiology, Speech & Language Pathology (AREA)
  • Human Computer Interaction (AREA)
  • Acoustics & Sound (AREA)
  • Multimedia (AREA)
  • Stereophonic System (AREA)
  • Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
  • Circuit For Audible Band Transducer (AREA)
  • Tone Control, Compression And Expansion, Limiting Amplitude (AREA)

Abstract

一种适于在每一频带中以动态增益修改音频信号的滤波器组经建构使得所感知的延迟较小,而在低频下施加较大群延时,以在较低频带中实现较高频率分辨率。可通过将全通滤波器插入到经重建滤波器响应中来实现在低频下的所述较高群延时。

Description

具有经改进频率分辨率的低延迟音频滤波器组
相关申请案的交叉引用
本申请案主张2020年5月22日申请的第63/028,966号美国临时专利申请案及2019年6月26日申请的第62/866,823号美国临时专利申请案的优先权的权益,所述申请案中的每一者的全部内容特此以引用的方式并入本文中。
技术领域
本发明涉及音频处理,且特定来说,涉及动态音频滤波器。
背景技术
除非本文中另有指示,否则本章节中所描述的方法并非本申请案中的权利要求书的现有技术且不因包含在本章节中而被承认是现有技术。
当处理输入音频信号时,通常期望使用滤波器组来更改音频信号使得通过一组动态系数指定在各种频率下的增益。可通过利用若干预先计算的带通滤波器来实施动态音频滤波器,其中音频滤波器响应由带通滤波器响应的加权和形成。权重可动态地变化。
一般来说,频带不具有相同带宽,而是较低频率滤波器比较高频率滤波器窄(且更紧密间隔)。因此,较高频率滤波器的脉冲响应将比较低频率滤波器的脉冲响应(其能量散布在较大量样本内)更紧凑(其能量散布在较少量样本内)。
发明内容
现有滤波器组的一个问题是延迟。针对较低频率滤波器,为了准确地表示脉冲响应,必须使用定义数目个样本,这与滤波器的延迟对应。较高频率滤波器不需要使用与较低频率滤波器一样多的样本,且因此具有较低延迟。然而,在含有低频滤波器及高频滤波器两者的滤波器组中,为了重组经滤波频带,需要使高频滤波器延时以匹配低频滤波器的延迟。因此,低频滤波器的延迟约束滤波器组的总延迟。
鉴于上文,需要在保留滤波器的脉冲响应的同时减少在低频下的延迟。本文中描述与低延迟滤波器设计相关的技术。
根据一实施例,一种音频处理方法包含从多个理想脉冲响应生成多个经修改脉冲响应,其中所述多个理想脉冲响应分别与多个频率对应,其中生成所述多个经修改脉冲响应包含:对所述多个理想脉冲响应中的至少一者执行衰落操作(fade operation)及逆时操作。所述方法进一步包含运用所述多个经修改脉冲响应对输入信号进行滤波以生成输出信号。
生成所述多个经修改脉冲响应可包含:基于第一理想脉冲响应生成前涟波响应;基于所述前涟波响应生成后涟波响应;及加上所述第一理想脉冲响应,减去所述前涟波响应,且加上所述后涟波响应以生成第一经修改脉冲响应。
生成所述多个经修改脉冲响应可包含:基于第一理想脉冲响应生成第一滤波器响应;基于所述第一滤波器响应生成前涟波响应;基于所述前涟波响应生成中间响应;基于所述中间响应生成第二滤波器响应;基于所述第二滤波器响应生成后涟波响应;及加上所述第一滤波器响应,加上所述后涟波响应,且减去所述前涟波响应以生成第一经修改脉冲响应。
根据另一实施例,一种设备包含处理器及存储器。所述处理器经配置以控制所述设备以从多个理想脉冲响应生成多个经修改脉冲响应,其中所述多个理想脉冲响应分别与多个频率对应,其中生成所述多个经修改脉冲响应包含:对所述多个理想脉冲响应中的至少一者执行衰落操作及逆时操作。所述处理器进一步经配置以控制所述设备以运用所述多个经修改脉冲响应对输入信号进行滤波而生成输出信号。所述设备可额外地包含与本文中所描述的方法的一或多者的细节类似的细节。
根据另一实施例,一种非暂时性计算机可读媒体存储计算机程序,所述计算机程序在由处理器执行时控制设备以执行包含本文中所描述的方法的一或多者的处理。
以下具体实施方式及附图提供对各种实施方案的性质及优点的进一步理解。
附图说明
图1是滤波器组100的框图。
图2是展示一组滤波器组频率响应的曲线图200。
图3是展示示范性组合滤波器响应302的曲线图300。
图4是滤波器组400的框图。
图5是滤波器组500的框图。
图6是脉冲响应600的曲线图。
图7是滤波器阵列700的框图。
图8是展示一组示范性时域滤波器组脉冲响应的曲线图800。
图9是展示修改脉冲响应的改进方法的信号的曲线图900。
图10是生成滤波器响应的方法1000的框图。
图11是展示与图10相关的各种信号的曲线图1100。
图12是展示相位失真1202及群延时1204的一组曲线图。
图13是展示各种脉冲响应的曲线图1300。
图14是生成滤波器响应的方法1400的框图。
图15是展示与图14相关的各种信号的曲线图1500。
图16是音频处理的方法1600的流程图。
图17是音频处理的方法1700的流程图。
图18是音频处理的方法1800的流程图。
具体实施方式
本文中描述与音频滤波器相关的技术。在以下描述中,出于说明的目的,阐述许多实例及特定细节以提供对本发明的透彻理解。然而,对于所属领域的技术人员来说将显而易见的是,如由权利要求书定义的本发明可单独地或与下文描述的其它特征组合地包含这些实例中的一些或所有特征,且可进一步包含本文中所描述的特征及概念的修改及等效物。
在以下描述中,详述各种方法、处理过程及过程。尽管可以特定顺序描述特定步骤,但这顺序主要是为了方便及清楚。特定步骤可重复一次以上,可在其它步骤之前或之后发生(即使原本以另一顺序描述所述步骤),且可与其它步骤并行发生。仅在必须在第二步骤开始之前完成第一步骤时,才要求第二步骤跟随第一步骤。当从上下文不清楚时,将明确指出此情形。
在本文献中,使用术语“及”、“或”及“及/或”。此类术语应被理解为具有包含性含义。例如,“A及B”可至少意味着以下每一者:“A及B两者”、“至少A及B两者”。作为另一实例,“A或B”可至少意味着以下每一者:“至少A”、“至少B”、“A及B两者”、“至少A及B两者”。作为另一实例,“A及/或B”可至少意味着以下每一者:“A及B”、“A或B”。当旨在排他性或时,这将被明确指出(例如,“A或B中的任一者”、“至多A及B中的一者”)。
本文献描述与例如块、元件、组件、电路等的结构相关联的各种处理功能。一般来说,这些结构可由受一或多个计算机程序控制的处理器实施。
图1是滤波器组100的框图。滤波器组100包含未个别地展示的若干滤波器。滤波器组100经配置具有若干权重102(也称为加权系数),其中权重102中的每一者与特定频带的增益对应。滤波器组100接收输入信号104,将权重102应用到输入信号104,且生成输出信号106。权重102可随时间变化。
输入信号104(及输出信号106)可为单通道信号,在所述情况中,滤波器组100包含若干滤波器,其中每一滤波器与权重102中的一者及特定频带对应。输入信号104(及输出信号106)可为多通道信号,在所述情况中,滤波器组100包含若干滤波器阵列,其中每一滤波器阵列与所述通道中的一者对应,且给定阵列中的每一滤波器与权重102中的一者及特定频带对应。输入信号104及输出信号106可为具有不同数目个通道的多通道信号,其中每一滤波器阵列与输入通道中的一者及输出通道中的一者对应。
一般来说,本文中所描述的滤波技术可应用到滤波器阵列100中的每一个别滤波器。
图2是展示一组滤波器组频率响应的曲线图200。在曲线图200中,x轴为频率且y轴为增益。例如,可预先计算一组B个带通脉冲响应h1(n)、h2(n)、...、hB(n),其具有带通中心频率fc1、fc2、...、fcB。理想地,滤波器hb(n)将具有频率响应Hb(f),其(近似)满足:
等式(1):
Figure BDA0003423219840000041
滤波器H2(f)的频率响应202被展示为在f=fc2下具有增益1,且在所有其它值f=fcb,b≠2下具有增益0。其它滤波器Hb(f)在其它频率fcb下具有增益1。频率fs/2为采样率(奈奎斯速率(Nyquist rate))的一半的频率。
接着,可使用一组加权系数w1、w2、...、wB来形成组合滤波器响应。等式的时域及频域版本如下:
等式(2):
时域:
Figure BDA0003423219840000042
等式(3):
频域:
Figure BDA0003423219840000051
图3是展示示范性组合滤波器响应302的曲线图300。x轴为频率(展示五个频率值)且y轴为增益(展示五个对应权重)。作为实例,在f=fc2下的增益304等于w2。其它频率fcb具有与其相应权重wb对应的增益。在实际应用中,此组合滤波器响应可以多种方式实施,如图4到5中展示。
图4是滤波器组400的框图。滤波器组400包含若干卷积块402、若干乘法块404及加法块406。滤波器组400接收输入信号410(被展示为X且也称为x(n),其参考给定样本n的输入信号),且生成输出信号412(被展示为Y且也称为y(n))。
卷积块402中的每一者大体上与与脉冲响应420(被展示为hb)相关联的给定频率、频带或频率范围对应。给定卷积块402将输入信号410(例如,x(n))与脉冲响应420(例如,hb(n))进行卷积以生成经滤波信号422(例如,xb(n))。
乘法块404中的每一者与权重430(被展示为wb)相关联且与卷积块402中的一者相关联(因此同样与频带相关联)。给定乘法块404将经滤波信号422(例如,xb(n))与权重430(例如,wb)相乘以生成经加权信号432。
加法块406将来自乘法块404中的每一者的经加权信号432相加以生成输出信号412。
表示这些运算的等式如下:
等式(4):
经滤波子带信号:
Figure BDA0003423219840000052
等式(5):
Figure BDA0003423219840000053
等式(6):
经滤波输出:
Figure BDA0003423219840000054
等式(4)中所使用的
Figure BDA0003423219840000055
运算符指示h(n)与x(n)的卷积。
等式(5)展示卷积运算导致生成子带信号x1(n)、x2(n)、...、xB(n)的细节。接着,将输出信号412(例如,y(n))形成为这些子带信号乘以其相应权重的总和,如等式(6)中所展示。
根据此方法,权重可为时变的,因此在等式(6)中,wb也可由wb(n)取代。
图5是滤波器组500的框图。滤波器组500包含若干乘法块502、加法块504,及卷积块506。滤波器组500接收输入信号510(被展示为X且也称为x(n),其参考给定样本n的输入信号),且生成输出信号512(被展示为Y且也称为y(n))。
乘法块502中的每一者大体上对应于与脉冲响应520(被展示为hb)相关联的给定频率、频带或频率范围,且也与权重530(被展示为wb)相关联。给定乘法块502将脉冲响应520(例如,hb(n))及加权系数530(例如,wb(n))相乘以生成经加权响应532中的一者。
加法块504将来自乘法块502中的每一者的经加权响应532相加以生成组合滤波器响应540。
卷积块506将输入信号510(例如,x(n))与组合滤波器响应540进行卷积以生成输出信号512(例如,y(n))。
一般来说,滤波器组500实施动态滤波器,所述动态滤波器执行滤波器响应的周期性计算。例如,可在(重叠)音频块中处理输入音频,且在块#k的时间处,可根据等式(3)来运算(频域中的)滤波器响应。
在一个实施例中,我们可使用在音频块之间具有平滑交叉衰落的重叠音频块来在频域中实施实际卷积运算,这可允许逐块地改变滤波器响应。
在一个常见应用中,滤波器组响应是因果滤波器,其中滤波器输出仅取决干过去的及现在的输入(使得针对每一b∈{1,2,...,B},脉冲响应hb(n)=0,
Figure BDA0003423219840000061
)。如果所有滤波器组响应皆为因果关系,那么组合响应(h(n),根据等式(2))也将为因果关系。
图6是脉冲响应600的曲线图。脉冲响应600起因于应用滤波器使得针对每一b∈{1,2,...,B},脉冲响应hb(n)=0,
Figure BDA0003423219840000062
其中L为正整数。实际上,可通过在滤波器实施方案中添加L个样本的额外延时而实施具有此性质的滤波器(使得将脉冲响应在时间上偏移达L个样本以允许将其实施为因果滤波器)。
图7是滤波器阵列700的框图。滤波器阵列700包含若干滤波器702(被展示为Filterny,nx)及若干加法块704。滤波器阵列700接收若干输入信号710(被展示为X1到Xnx)且生成若干输出信号712(被展示为Y1到Yny)。一般来说,滤波器阵列700说明多通道输入、多通道输出系统的操作,其中滤波器将每一输入连接到每一输出。滤波器700还可被概念化为多通道滤波器组,其中针对每一频带,存在定义输入的混合以形成输出的[ny乘nx]二维阵列。
滤波器702布置成组,其中每一滤波器组与输入信号710中的一者(例如,通道)相关联。滤波器702中的每一者一般与频带、滤波器响应及权重相关联,这与本文中所论述的其它滤波器类似。例如,滤波器702中的每一者可为例如滤波器组400(参见图4)、滤波器组500(参见图5)等的滤波器。
就频域滤波器响应来说,我们可依据以下等式描述滤波器阵列700的处理:
等式(7):
Figure BDA0003423219840000071
本文中所描述的实施例涉及用于实施滤波器702中的一或多者的方法。将明白,实施例同样适用于如图7中所展示的滤波器阵列。
返回参考图3,展示示范性滤波器的所要频率响应,其中根据一组预定义预定义频率fc1、fc2、…、及对应增益值w1、w2、…定义依据频率而变化的滤波器响应302的增益。
在图3的实例中,滤波器在频率fc2下的增益由w2设置,如被展示为增益204。
在实施例中,通过对若干预定义滤波器组响应进行加权求和而实现图3的频率响应。示范性响应在图2中展示,其中展示频带2的频率响应202。
我们将这些预定义滤波器组响应中的每一者的响应称为Hb(f)(在频域中)或替代地称为hb(n)(在时域中)。
所要滤波器响应(例如图3中的302)可由预定义滤波器组响应的加权和形成。这可被表达为时域或频域求和:
等式(8):
时域:
Figure BDA0003423219840000081
等式(9):
频域:
Figure BDA0003423219840000082
我们可选择坚持认为预定义滤波器组响应中的每一者应表示因果滤波器(使得h(n)=0,
Figure BDA0003423219840000083
)。然而,为了方便起见,我们代替性地坚持认为h(n)=0,
Figure BDA0003423219840000084
且因此我们将被要求将L个延迟样本添加到我们的滤波器的实施方案,使得最终(可实现)滤波器输出y′(n)将根据以下计算为理想响应y(n)的延时版本:
等式(10):
Figure BDA0003423219840000085
等式(11):
y′(n)=y(n-L)
等式(12):
Figure BDA0003423219840000086
应注意,理想信号y(n)的输出样本n的计算使用至多x(n+L)的输入样本(包含L个未来输入样本)来运算。相比之下,经延时信号y′(n)的输出样本n的计算仅使用至多x(n)的输入样本来运算。
在图6中展示具有性质h(n)=0,
Figure BDA0003423219840000087
的示范性脉冲响应。
在实施例中,滤波器组的频率分辨率(相邻滤波器组中心频率之间的间隔,例如fcb+1-fcb)较小,而延迟L也较小。
在一个特定实施例中,低延迟滤波器组用于提供低延迟L及高分辨率。在另一实施例中,全通低延迟滤波器组用于提供低感知延迟L及高分辨率。
理想滤波器组响应
可通过多种手段来建构滤波器组,包含其中如图2中所展示般生成理想频率响应且接着将其转换成线性相位的理想时域脉冲响应的方法。
图8是展示一组示范性时域滤波器组脉冲响应的曲线图800。在曲线图800中,x轴为时间(以样本为单位)且y轴为频率。在曲线图800中展示许多滤波器响应,其中每一滤波器响应与如通过名称F1、F2、...、FB指示的给定频率对应,其中滤波器组中的滤波器的数目在此实例中为B=10。最高数目(例如,10)与最高频率对应且最低数目(例如,1)与最低频率对应。滤波器Fb的理想脉冲响应为h′b(n)(针对b∈{1,2,...,B})。为了进行下文的进一步论述而指示两个滤波器响应802及804;滤波器响应802与第三滤波器(h′3(n))对应且滤波器响应804与最低频率滤波器响应(h′1(n))对应。
根据常见实践,理想滤波器组响应展现完美求和性质,如下:
等式(13):
Figure BDA0003423219840000091
等式(14):
其中
Figure BDA0003423219840000092
优选地,我们可利用其中较低频率滤波器比较高频率滤波器窄(且更紧密间隔)的滤波器组,如图2中所说明,例如,其中fc1与fc2之间的间隔小于fC4与fc5之间的间隔。因此,较高频率滤波器的脉冲响应将比较低频率滤波器的脉冲响应更紧凑(其能量散布在较少量样本内),如图8中可观察到。
理论上,理想滤波器组的脉冲响应的长度可为无限的,且特定来说在图8中可见,滤波器F3的脉冲响应802不满足所要低延迟性质:h′3(n)=0,
Figure BDA0003423219840000093
因此,图8中所展示的滤波器不适于实施低延迟滤波器组。
用于低延迟滤波器响应的技术
根据一种方法,可通过简单地截断理想脉冲响应而生成低延迟滤波器组(具有延迟L),如下:
等式(15):
Figure BDA0003423219840000101
然而,此截断将影响在低频下的所得响应,且在真实世界系统中使用的情况下可导致听觉伪差(auditory artifact)。
图9是展示修改脉冲响应的改进方法的信号的曲线图900。在曲线图900中,x轴为时间(以样本为单位)。y轴为量值;应注意,每一信号在y轴上彼此独立(例如,y轴独立地展示每一信号的量值范围,而非信号之间的任何种类的相对量值或比较)。信号802与滤波器F3的理想脉冲响应(例如,h′3(n))对应,也如图8中所展示。信号904与衰落函数对应。信号906与信号802乘以信号904对应以生成低延迟滤波器响应,如下:
等式(16):
hb(n)=h′b(n)q(n)
等式(17):
其中:
Figure BDA0003423219840000102
所属领域的技术人员将明白,根据两种技术(等式(15)及等式(16)),所得脉冲响应hb(n)具有期望性质hb(0)=h′b(0),且这保证将满足完美重建准则(等式(13))。
然而,这些滤波器的频率响应(如由等式(15)或等式(16)产生)将不一定匹配原始(理想)响应h′b(n)且可运用改进技术来实现,如下文所论述。
用于低延迟滤波器响应的改进方法
图10是生成滤波器响应的方法1000的框图。方法1000可用于生成对本文中所论述的滤波器组(例如,图4的滤波器组400、图5的滤波器组500、图7的滤波器阵列700中的滤波器702等)中的任一者中的滤波器中的一或多者的滤波器响应。方法1000包含衰落块1002、逆时块1004及求和块1006。方法1000的块可以各种方式来实施,例如由电路元件、由执行一或多个计算机程序的处理器等。
衰落块1002接收理想脉冲响应1010,应用衰落函数,且生成前涟波响应1012。理想脉冲响应1010可与上文所论述的理想脉冲响应中的一或多者(例如,如图8或9中所论述的h′b(n)对应,例如信号802)。在图11中展示衰落函数的实例。
逆时块1004对前涟波响应1012执行逆时操作以生成后涟波响应1014。逆时操作大体上与围绕信号的零样本镜像前涟波响应1012对应。在图11中展示逆时操作的实例。
求和块1006加上理想脉冲响应1010,减去前涟波响应1012,且加上后涟波响应1014,以生成滤波器响应1016。
图11是展示与图10相关的各种信号的曲线图1100。在曲线图1100中,x轴为时间(以样本为单位)且y轴为量值;正如图9,应注意,每一信号的量值与每一其它信号无关。信号802与滤波器F3的理想脉冲响应(例如,h′3(n))对应,也如图8到9中所展示。也参见图10中的理想脉冲响应1010。
信号1102与衰落函数对应。也参见图9中的信号904。信号1104与信号802乘以信号1102对应以生成前涟波响应。也参见图10中的前涟波响应1012。
信号1106与信号1104的逆时对应,称为后涟波响应。也参见图10中的后涟波响应1014。
信号1108与滤波器组响应对应,且是通过加上信号802(理想脉冲响应1010),减去信号1104(前涟波响应1012)且加上信号1106(后涟波响应1014)而生成。也参见图10中的滤波器响应1016。
以下等式展示此处理过程:
等式(18):
前涟波:
Figure BDA0003423219840000111
等式(19):
衰落函数:
Figure BDA0003423219840000112
等式(20):
后涟波:
Figure BDA0003423219840000113
等式(21):
低延迟响应:
Figure BDA0003423219840000114
应注意,根据等式(20),
Figure BDA0003423219840000115
且因此,低延迟脉冲响应hb(n)(根据等式(21)将具有期望性质hb(0)=h′b(0),且这保证将满足完美重建准则(等式(13))。
低延迟全通滤波器响应
可在收听者未感知到滤波的群延时的情况下使用各种类型的音频滤波来修改音频信号。作为实例,典型高通滤波器(例如用于从录制声音移除不想要的低频噪声的滤波器)将在低频下引入将不被收听者认为有害的相位失真。收听者将无法感知的另一常见类型的滤波为全通滤波器。
图12是展示相位失真1202及群延时1204的一组曲线图。相位失真1202及群延时1204起因于使用如等式(22)中所展示的传递函数处理以48kHz取样的音频信号。
等式(22):
Figure BDA0003423219840000121
针对以每秒48000个样本的采样率操作的滤波器,等式(22)中的常数0.994986与大约37Hz的过渡频率的典型值(其中一阶全通滤波器的相移为90度)对应。可选择此常数的替代值,且针对以每秒48000个样本操作的滤波器,常数可在0.9490与0.9993之间变化,所述值分别与400Hz及5Hz的过渡频率对应。
假定此相位失真1202在输出音频信号中产生收听者可接受的伪差,那么因此可接受将相同相位失真应用到理想脉冲响应以产生一组新的全通滤波器响应,如图13中所展示。
图13是展示各种脉冲响应的曲线图1300。在曲线图1300中,x轴为时间(以样本为单位)且y轴在各种频率Fn下布置脉冲响应。滤波器F3的脉冲响应1302看起来与图8中的滤波器F3的原始(理想)响应802类似。相比之下,滤波器F1的脉冲响应1304看起来与图8中的滤波器F1的原始(理想)响应804不同。这差异归因于全通滤波器的效应,其影响较低频率滤波器而对较高频率滤波器的影响可被忽略不计。下文更详细地论述全通滤波器。
可见,与区域1308(在时间零之后的100个样本以上)相比,滤波器F1的脉冲响应1304在区域1306(在时间零之前的100个样本以上)中具有更低振幅。脉冲响应的能量向脉冲响应中的较晚时间偏移为全通滤波器的理想属性,因为其产生更紧密地接近延迟约束(其中在图13的实例中,L=100)的滤波器组。
图14是生成滤波器响应的方法1400的框图。方法1400可用于生成对本文中所论述的滤波器组(例如,图4的滤波器组400、图5的滤波器组500、图7的滤波器阵列700中的滤波器702等)中的任一者中的滤波器中的一或多者的滤波器响应。方法1400包含全通滤波器1402、衰落块1404、逆时块1406、全通滤波器1408、全通滤波器1410及求和块1412。方法1400的块可以各种方式来实施,例如由电路元件、由执行一或多个计算机程序的处理器等。
全通滤波器1402接收理想脉冲响应1420,执行全通滤波,且生成滤波器响应1422。理想脉冲响应1420可与上文所论述的理想脉冲响应中的一或多者(例如,如图8或9中所论述的h′b(n),例如信号802)对应。在下文等式中更详细地论述全通滤波器的脉冲响应。
衰落块1404将衰落函数应用到滤波器响应1422以生成前涟波响应1424。在图15中展示衰落函数的实例。
逆时块1406对前涟波响应1424执行逆时操作以生成中间响应1426。在图15中展示逆时操作的实例。
全通滤波器1408对中间响应1426执行全通滤波以生成滤波器响应1428。
全通滤波器1410对滤波器响应1428执行全通滤波以生成后涟波响应1430。
求和块1412加上滤波器响应1422,加上后涟波响应1430,且减去前涟波响应1424,以生成滤波器响应1432。
在实施例中,为了正确地满足延迟约束,我们在调适计算的情况下应用等式(21)的处理过程,以允许引入全通滤波器。在以下等式中,标记法
Figure BDA0003423219840000131
用于指示脉冲响应a(n)及b(n)的卷积。低延迟全通滤波器响应的生成如下实行:
等式(23):
全通:
Figure BDA0003423219840000132
等式(24):
前涟波:
Figure BDA0003423219840000133
等式(25):
衰减函数:
Figure BDA0003423219840000134
等式(26):
后涟波:
Figure BDA0003423219840000141
等式(27):
低延迟响应:
Figure BDA0003423219840000142
在上文等式中,g(n)为全通滤波器的脉冲响应(也参见等式(22))。
在实施例中,全通滤波器是在10Hz与200Hz之间的频率下具有单个实极点的一阶滤波器。针对对以每秒Fs个样本的速率取样的音频操作的具有延迟L的滤波器组,全通极点的频率可被近似设置为:
等式(28):
全通极点频率:
Figure BDA0003423219840000143
如等式(28)中所展示的极点频率的计算为实例。将明白,其它极点频率也将在滤波器组的设计中提供改进延迟/准确度折衷的优点。
当将低延迟全通滤波器组中的B个滤波器的组的脉冲响应加总在一起时,最终结果(重建脉冲响应)近似等于全通响应:
等式(29):
Figure BDA0003423219840000144
等式(30):
≈g(n)
图14的过程可应用到滤波器组中的每一滤波器,因此例如针对b∈{1,2,...,B},可根据图14处理理想脉冲响应h′b(n)以产生低延迟全通滤波器hb(n)。
图15是展示与图14相关的各种信号的曲线图1500。在曲线图1500中,x轴为时间(以样本为单位)且y轴为量值;正如图9及11,应注意,每一信号的量值与每一其它信号无关。信号1304与滤波器F1的全通滤波脉冲响应(例如,h′1(n))对应,也如图13中所展示。也参见图14中的滤波器响应1422及等式(23)。
信号1502与衰落函数对应。也参见图9中的信号904及图11中的信号1102及等式(25)。信号1504与信号1304乘以信号1502以生成前涟波响应对应。也参见图14中的前涟波响应1424及等式(24)。
信号1506与信号1504的逆时对应,称为中间响应。也参见图14中的中间响应1426。信号1508与将全通滤波器应用到信号1506两次对应,称为后涟波响应。也参见图14的后涟波响应1430及等式(26)。
信号1510与滤波器组响应对应,且通过加上信号1304(滤波器响应1422),加上信号1508(后涟波响应1430),且减去信号1504(前涟波响应1424)来生成。也参见等式(27)。
图10的方法1000与图14的方法1400之间的区别在于可如下那样概念化。首先,应注意,方法1000对滤波器响应h′b(n)进行操作,且方法1400对全通滤波版本
Figure BDA0003423219840000151
进行操作。方法1000的处理过程依赖于以下事实:前涟波分量及后涟波分量是通过围绕t=0进行镜像且结合我们已知理想脉冲响应1010关于时间零对称而得到。
将此知识扩展到方法1400,想象我们可得到新滤波器分量:
Figure BDA0003423219840000152
Figure BDA0003423219840000153
这意味着着我们采用原始pre_ripple(例如,中间响应1426或信号1506),且应用全通函数g的逆函数。
接着,我们得到另一新滤波器分量:postA(t)=preA(-t)。这简单为preA的镜像版本。
现在,我们可得到另一滤波器分量:finalA=h′b(n)+postA-preA。此步骤非常像图10中所使用的处理过程,且最重要的是,此步骤遵循我们的口诀(mantra):前涟波分量及后涟波分量是通过围绕t=0进行镜像且结合我们已知原始滤波器关于时间零对称而得到。
接着,我们可建立最终滤波器:
Figure BDA0003423219840000154
Figure BDA0003423219840000155
现在,我们知道:
Figure BDA0003423219840000156
Figure BDA0003423219840000157
(这是因为全通滤波器的时间镜像与全通滤波器的逆滤波器相同而成立。)
因此,我们得出:
Figure BDA0003423219840000158
Figure BDA0003423219840000159
图16是音频处理的方法1600的流程图。方法1600可由例如根据一或多个计算机程序执行指令的处理器来实施。
在1602,从理想脉冲响应生成经修改脉冲响应。理想脉冲响应分别与若干频率对应。生成经修改脉冲响应包含对理想脉冲响应中的至少一者执行衰落操作及逆时操作。例如,方法1000(参见图10)可对理想脉冲响应1010使用衰落块1002执行衰落操作且使用逆时块1004执行逆时操作。作为另一实例,方法1400(参见图14)可对理想脉冲响应1420使用衰落块1404执行衰落操作且使用逆时块1406执行逆时操作。
在1604,运用经修改脉冲响应(参见1602)对输入信号进行滤波以生成输出信号。例如,滤波器组100(参见图1)可包含具有在1602生成的经修改脉冲响应的若干滤波器,且可从输入信号104生成输出信号106。
在图17到18中进一步详述方法1600的一般步骤。
图17是音频处理的方法1700的流程图。方法1700可由例如根据一或多个计算机程序执行指令的处理器来实施。方法1700与方法1600(参见图16)类似,其具有特定于方法1000(参见图10)的操作的更多细节。
在1702,从理想脉冲响应生成经修改脉冲响应。这包含子步骤1704到1708。
在1704,基于第一理想脉冲响应生成前涟波响应。例如,衰落块1002(参见图10)可对理想脉冲响应1010执行衰落操作以生成前涟波响应1012。也参见等式(18)。
在1706,基于前涟波响应生成后涟波响应。例如,逆时块1004(参见图10)可对前涟波响应1012执行逆时操作以生成后涟波响应1014。也参见等式(20)。
在1708,通过加上第一理想脉冲响应,减去前涟波响应且加上后涟波响应而生成第一经修改脉冲响应。例如,求和块1006(参见图10)可加上理想脉冲响应1010,减去前涟波响应1012,且加上后涟波响应1014而生成滤波器响应1016。也参见等式(21)。
在1710,运用经修改脉冲响应(参见1702到1708)对输入信号进行滤波以生成输出信号。例如,滤波器组100(参见图1)可包含具有在1702生成的经修改脉冲响应的若干滤波器,且可从输入信号104生成输出信号106。
图18是音频处理的方法1800的流程图。方法1800可通过例如根据一或多个计算机程序执行指令的处理器来实施。方法1800与方法1600(参见图16)类似,其具有特定于方法1400(参见图14)的操作的更多细节。
在1802,从理想脉冲响应生成经修改脉冲响应。这包含子步骤1804到1814。
在1804,基于第一理想脉冲响应生成第一滤波器响应。例如,全通滤波器1402(参见图14)可对理想脉冲响应1420执行全通滤波操作以生成滤波器响应1422。也参见等式(23)。
在1806,基于第一滤波器响应生成前涟波响应。例如,衰落块1404(参见图14)可对滤波器响应1422执行衰落操作以生成前涟波响应1424。也参见等式(24)。
在1808,基于前涟波响应生成中间响应。例如,逆时块1406(参见图14)可对前涟波响应1424执行逆时操作以生成中间响应1426。也参见等式(26)。
在1810,基于中间响应生成第二滤波器响应。例如,全通滤波器1408(参见图14)可对中间响应1426执行全通滤波操作以生成滤波器响应1428。也参见等式(27)。
在1812,基于第二滤波器响应生成后涟波响应。例如,全通滤波器1410(参见图14)可对滤波器响应1428执行全通滤波操作以生成后涟波响应1430。也参见等式(27)。
在1814,通过加上第一滤波器响应,加上后涟波响应且减去前涟波响应而生成第一经修改脉冲响应。例如,求和块1412(参见图14)可加上滤波器响应1422,加上后涟波响应1430,且减去前涟波响应1424以生成滤波器响应1432。也参见等式(27)。
在1816,运用经修改脉冲响应(参见1802到1814)对输入信号进行滤波以生成输出信号。例如,滤波器组100(参见图1)可包含具有在1802生成的经修改脉冲响应的若干滤波器,且可从输入信号104生成输出信号106。
在实施例中,可处理输入音频信号以产生输出音频信号,其中在若干预定频率下处理的增益近似等于一组对应增益系数。处理可根据滤波器响应实施,所述滤波器响应是对应数目个低延迟全通滤波器响应的加权和,其中权重是通过对应增益系数确定。低延迟全通滤波器响应已根据图18的方法确定。
在替代实施例中,处理输入音频信号以产生输出音频信号,其中所述输出音频信号是通过以总脉冲响应对所述输入音频信号进行滤波而形成,所述总脉冲响应是一组滤波器组脉冲响应的加权和。所述滤波器组脉冲响应的总和是全通滤波器。全通滤波器可在近似500Hz的频率以上具有近似恒定的相位响应,及在较低频率下上升的群延时。
在替代实施例中,处理输入音频信号以产生输出音频信号,其中所述输出音频信号是通过以总脉冲响应对所述输入音频信号进行滤波而形成,所述总脉冲响应是一组滤波器组脉冲响应的加权和。重建脉冲响应可由滤波器组脉冲响应的总和形成,使得重建脉冲响应的群延时在高于500Hz的频率范围内近似恒定且在低于500Hz时增加。例如,群延时在高于500Hz时可在感知上恒定(变化不超过0.1ms),且在较低频率(例如,小于500Hz)下可变化不超过0.5/f秒。来自所述组滤波器组脉冲响应的一或多个脉冲响应的带宽可根据重建脉冲响应的增加的群延时而更窄。
实施方案细节
实施例可以硬件、存储在计算机可读媒体上的可执行模块或两者的组合(例如,可编程逻辑阵列)来实施。除非另有指定,否则由实施例执行的步骤不需要固有地与任何特定计算机或其它设备相关,尽管它们可在特定实施例中。特定来说,各种通用机器可与根据本文中的教示写入的程序一起使用,或可方便地建构更专用设备(例如,集成电路)以执行所需方法步骤。因此,实施例可以在相应包括至少一个处理器、至少一个数据存储系统(包含易失性及非易失性存储器及/或存储元件)、至少一个输入设备或端口及至少一个输出装置或端口的一或多个可编程计算机系统上执行的一或多个计算机程序来实施。将程序代码应用到输入数据以执行本文中所描述的功能且生成输出信息。将输出信息以已知方式应用到一或多个输出装置。
每一此计算机程序优选地存储在通用或专用可编程计算机可读的存储媒体或装置(例如,固态存储器或媒体,或磁性或光学媒体)上或下载到所述存储媒体或装置,以在计算机系统读取所述存储媒体或装置时配置及操作计算机以执行本文中所描述的程序。本发明系统也可被视为实施为配置有计算机程序的计算机可读存储媒体,其中如此配置的存储媒体引起计算机系统以特定且预定的方式操作以执行本文中所描述的功能。(软件本身及无形或暂时性信号在其为不能取得专利的标的物的情况下被排除在外。)
上文描述说明本发明的各个实施例以及可如何实施本发明的方面的实例。上文实例及实施例不应被视为仅有实施例,而是为了说明如以下权利要求书定义的本发明的灵活性及优点而呈现。基于上文公开内容及以下权利要求书,其它布置、实施例、实施方案及等效物对于所属领域的技术人员来说将为显而易见的且可在不脱离如由权利要求书定义的本发明的精神及范围的情况下采用。

Claims (20)

1.一种音频处理方法,所述方法包括:
从多个理想脉冲响应生成多个经修改脉冲响应,其中所述多个理想脉冲响应分别与多个频率对应,其中生成所述多个经修改脉冲响应包含:对所述多个理想脉冲响应中的至少一者执行衰落操作及逆时操作;及
运用所述多个经修改脉冲响应对输入信号进行滤波以生成输出信号。
2.根据权利要求1所述的方法,其中生成所述多个经修改脉冲响应包含:
基于第一理想脉冲响应来生成前涟波响应;
基于所述前涟波响应来生成后涟波响应;及
加上所述第一理想脉冲响应,减去所述前涟波响应,且加上所述后涟波响应以生成第一经修改脉冲响应。
3.根据权利要求2所述的方法,其中生成所述前涟波响应包括:
对所述第一理想脉冲响应执行所述衰落操作以生成所述前涟波响应。
4.根据权利要求2到3中任一权利要求所述的方法,其中生成所述后涟波响应包括:
对所述前涟波响应执行所述逆时操作以生成所述后涟波响应。
5.根据权利要求1所述的方法,其中生成所述多个经修改脉冲响应包含:
基于第一理想脉冲响应来生成第一滤波器响应;
基于所述第一滤波器响应来生成前涟波响应;
基于所述前涟波响应来生成中间响应;
基于所述中间响应来生成第二滤波器响应;
基于所述第二滤波器响应来生成后涟波响应;及
加上所述第一滤波器响应,加上所述后涟波响应,且减去所述前涟波响应以生成第一经修改脉冲响应。
6.根据权利要求5所述的方法,其中生成所述第一滤波器响应包括:
对所述第一理想脉冲响应执行第一全通滤波操作以生成所述第一滤波器响应。
7.根据权利要求6所述的方法,其中所述第一全通滤波操作将传递函数Gap应用到所述第一理想脉冲响应,其中
Figure FDA0003423219830000021
8.根据权利要求6到7中任一权利要求所述的方法,其中所述第一全通滤波操作是具有在10Hz与200Hz之间的单个实极点的一阶滤波器。
9.根据权利要求6到8中任一权利要求所述的方法,其中所述第一全通滤波操作具有以每秒Fs个样本的速率取样的音频信号的L个样本的延迟,且其中根据等式来计算所述全通滤波操作的极点的频率f
Figure FDA0003423219830000022
10.根据权利要求5到9中任一权利要求所述的方法,其中生成所述前涟波响应包括:
对所述第一滤波器响应执行所述衰落操作以生成所述前涟波响应。
11.根据权利要求5到10中任一权利要求所述的方法,其中生成所述中间响应包括:
对所述前涟波响应执行所述逆时操作以生成所述中间响应。
12.根据权利要求5到11中任一权利要求所述的方法,其中生成所述第二滤波器响应包括:
对所述中间响应执行第二全通滤波操作以生成所述第二滤波器响应。
13.根据权利要求5到12中任一权利要求所述的方法,其中生成所述后涟波响应包括:
对所述第二滤波器响应执行第三全通滤波操作以生成所述后涟波响应。
14.根据权利要求1到13中任一权利要求所述的方法,其进一步包括:
通过将多个权重应用到所述多个经修改脉冲响应来生成多个经加权的经修改脉冲响应,
其中对所述输入信号进行滤波包括:运用所述多个经加权的经修改脉冲响应对所述输入信号进行滤波以生成所述输出信号。
15.根据权利要求14所述的方法,其中重建脉冲响应与所述经加权的经修改脉冲响应的总和对应,所述重建脉冲响应在高于500Hz的频率范围内近似恒定,且在低于500Hz时增加。
16.根据权利要求14到15中任一权利要求所述的方法,其中所述多个经加权的经修改脉冲响应中的至少一者的带宽根据所述重建脉冲响应的增加的群延时为更窄。
17.根据权利要求14到16中任一权利要求所述的方法,其中所述多个权重是时变的。
18.根据权利要求1到17中任一权利要求所述的方法,其中所述输入信号是多个输入信号中的一者,其中所述多个经修改脉冲响应通过多个滤波器组施加,其中所述输出信号是多个输出信号中的一者,且其中给定滤波器组对给定输入信号进行滤波以生成给定输出信号。
19.一种非暂时性计算机可读媒体,其存储计算机程序,所述计算机程序在由处理器执行时控制设备以执行包含根据权利要求1到18中任一权利要求所述的方法的处理。
20.一种用于音频处理的设备,所述设备包括:
处理器;及
存储器,
其中所述处理器经配置以控制所述设备以从多个理想脉冲响应生成多个经修改脉冲响应,其中所述多个理想脉冲响应分别与多个频率对应,其中生成所述多个经修改脉冲响应包含:对所述多个理想脉冲响应中的至少一者执行衰落操作及逆时操作;及
其中所述处理器经配置以控制所述设备以运用所述多个经修改脉冲响应对输入信号进行滤波以生成输出信号。
CN202080045434.7A 2019-06-26 2020-06-25 具有经改进频率分辨率的低延迟音频滤波器组 Pending CN114051636A (zh)

Applications Claiming Priority (5)

Application Number Priority Date Filing Date Title
US201962866823P 2019-06-26 2019-06-26
US62/866,823 2019-06-26
US202063028966P 2020-05-22 2020-05-22
US63/028,966 2020-05-22
PCT/US2020/039472 WO2020264064A1 (en) 2019-06-26 2020-06-25 Low latency audio filterbank having improved frequency resolution

Publications (1)

Publication Number Publication Date
CN114051636A true CN114051636A (zh) 2022-02-15

Family

ID=71662320

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202080045434.7A Pending CN114051636A (zh) 2019-06-26 2020-06-25 具有经改进频率分辨率的低延迟音频滤波器组

Country Status (13)

Country Link
US (1) US20220231669A1 (zh)
EP (1) EP3991294A1 (zh)
JP (1) JP2022537654A (zh)
KR (1) KR20220024761A (zh)
CN (1) CN114051636A (zh)
AU (1) AU2020301244A1 (zh)
BR (1) BR112021023987A2 (zh)
CA (1) CA3142163A1 (zh)
IL (1) IL288362A (zh)
MX (1) MX2021015712A (zh)
SG (1) SG11202113159VA (zh)
TW (1) TW202105908A (zh)
WO (1) WO2020264064A1 (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA3211264A1 (en) * 2021-03-11 2022-09-15 Sascha Disch Audio decorrelator, processing system and method for decorrelating an audio signal

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101169934A (zh) * 2006-10-24 2008-04-30 华为技术有限公司 时域听觉阈值加权滤波器的构造方法和设备、编解码器
CN102349235A (zh) * 2009-04-10 2012-02-08 杜比实验室特许公司 利用前向-后向滤波获得期望的非零相移
TW201234753A (en) * 2011-02-14 2012-08-16 Richtek Technology Corp Frequency control circuit and method for a ripple regulator
CN103155591A (zh) * 2010-10-14 2013-06-12 杜比实验室特许公司 使用自适应频域滤波和动态快速卷积的自动均衡
US20130287226A1 (en) * 2012-04-30 2013-10-31 Conexant System, Inc. Reduced-delay subband signal processing system and method
US9026451B1 (en) * 2012-05-09 2015-05-05 Google Inc. Pitch post-filter
GB201522560D0 (en) * 2015-12-21 2016-02-03 Craven Peter G And Law Malcolm Lossless bandsplitting and bandjoining using allpass filters
CN105393553A (zh) * 2013-03-26 2016-03-09 拉克伦·保罗·巴拉特 虚拟采样率增加的音频滤波
CN109525218A (zh) * 2009-02-18 2019-03-26 杜比国际公司 低延迟调制滤波器组

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3236626B1 (en) * 2016-04-21 2020-09-23 Institut Mines Telecom / Telecom Bretagne Filter for linear modulation based communication systems

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101169934A (zh) * 2006-10-24 2008-04-30 华为技术有限公司 时域听觉阈值加权滤波器的构造方法和设备、编解码器
CN109525218A (zh) * 2009-02-18 2019-03-26 杜比国际公司 低延迟调制滤波器组
CN102349235A (zh) * 2009-04-10 2012-02-08 杜比实验室特许公司 利用前向-后向滤波获得期望的非零相移
CN103155591A (zh) * 2010-10-14 2013-06-12 杜比实验室特许公司 使用自适应频域滤波和动态快速卷积的自动均衡
TW201234753A (en) * 2011-02-14 2012-08-16 Richtek Technology Corp Frequency control circuit and method for a ripple regulator
US20130287226A1 (en) * 2012-04-30 2013-10-31 Conexant System, Inc. Reduced-delay subband signal processing system and method
US9026451B1 (en) * 2012-05-09 2015-05-05 Google Inc. Pitch post-filter
CN105393553A (zh) * 2013-03-26 2016-03-09 拉克伦·保罗·巴拉特 虚拟采样率增加的音频滤波
GB201522560D0 (en) * 2015-12-21 2016-02-03 Craven Peter G And Law Malcolm Lossless bandsplitting and bandjoining using allpass filters

Also Published As

Publication number Publication date
BR112021023987A2 (pt) 2022-01-25
EP3991294A1 (en) 2022-05-04
AU2020301244A1 (en) 2022-01-06
TW202105908A (zh) 2021-02-01
MX2021015712A (es) 2022-02-03
KR20220024761A (ko) 2022-03-03
SG11202113159VA (en) 2021-12-30
CA3142163A1 (en) 2020-12-30
US20220231669A1 (en) 2022-07-21
WO2020264064A1 (en) 2020-12-30
IL288362A (en) 2022-01-01
JP2022537654A (ja) 2022-08-29

Similar Documents

Publication Publication Date Title
JP5718251B2 (ja) 分解されたオーディオ信号の再構成のためのシステムおよび方法
JP6487383B2 (ja) 反響装置およびオーディオ信号を反響させる方法
US8150065B2 (en) System and method for processing an audio signal
EP2530840B1 (en) Efficient sub-band adaptive FIR-filtering
WO2010019192A1 (en) Sound field widening and phase decorrelation system and method
TW202116077A (zh) 使用頻譜正交音訊分量之副頻帶空間及串擾處理
KR20220044566A (ko) 심리음향적 주파수 범위 확장을 위한 비선형 적응성 필터뱅크
CN114051636A (zh) 具有经改进频率分辨率的低延迟音频滤波器组
US9036752B2 (en) Low-delay filtering
TWI421858B (zh) 用於處理音頻訊號的系統及方法
RU2807607C2 (ru) Банк аудиофильтров с малой задержкой и повышенной разрешающей способностью по частоте
US11488574B2 (en) Method and system for implementing a modal processor
WO2003047097A1 (fr) Procede de conception d'un filtre numerique, appareil de projet, programme de projet pour la realisation d'un tel filtre, et filtre numerique y relatif
Isen DSP for MATLABTM and LabVIEWTM I: Fundamentals of Discrete Signal Processing
Deppisch Plug-In for Frequency Dependent Control of Microphone Polar Patterns
KR20240011251A (ko) 심리음향 주파수 범위 확장을 위해 스케일 의존적 비선형성을 사용하는 적응형 필터뱅크
TW202334938A (zh) 正交鏡像濾波器域中之沉浸式音訊及視訊服務空間重建濾波器庫
Korohoda et al. A Practical Multirate FIR-Based Equalizer-A Case Study
JP6289041B2 (ja) イコライザ
Shah An Efficient and Reconfigurable Hardware Architecture for Narrow Bandwidth Low-pass and Polyphase Recursive Filters
Lin et al. High-Definition Audio Spectrum Separation via Frequency Response Masking Filter Banks
US20110007917A1 (en) Two Chip Solution Band Filtering
Akbaş A Zero-Assignment Approach to Two-Channel Filter Banks and Wavelets
Chaguaro Aldaz Digital Signal Processing Filtering Algorithm: Audio Equalization Using Matlab
Lin Neuro-fuzzy systems with nonuniform spectral structure for multichannel audio synthesis and delivery

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
REG Reference to a national code

Ref country code: HK

Ref legal event code: DE

Ref document number: 40065139

Country of ref document: HK