CN102195616A - 高效子带自适应有限脉冲响应滤波 - Google Patents
高效子带自适应有限脉冲响应滤波 Download PDFInfo
- Publication number
- CN102195616A CN102195616A CN2011100500367A CN201110050036A CN102195616A CN 102195616 A CN102195616 A CN 102195616A CN 2011100500367 A CN2011100500367 A CN 2011100500367A CN 201110050036 A CN201110050036 A CN 201110050036A CN 102195616 A CN102195616 A CN 102195616A
- Authority
- CN
- China
- Prior art keywords
- subband
- signal
- filter
- fir filter
- filter coefficient
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000001914 filtration Methods 0.000 title claims abstract description 11
- 230000003044 adaptive effect Effects 0.000 title description 13
- 238000000034 method Methods 0.000 claims abstract description 36
- 238000001228 spectrum Methods 0.000 claims abstract description 8
- 230000004044 response Effects 0.000 claims description 33
- 238000005070 sampling Methods 0.000 claims description 9
- 238000013461 design Methods 0.000 claims description 8
- 230000008859 change Effects 0.000 claims description 6
- 238000011144 upstream manufacturing Methods 0.000 claims description 6
- 230000001419 dependent effect Effects 0.000 abstract 2
- 230000006870 function Effects 0.000 description 36
- 238000012546 transfer Methods 0.000 description 23
- 230000006978 adaptation Effects 0.000 description 22
- 238000012545 processing Methods 0.000 description 12
- 238000009826 distribution Methods 0.000 description 11
- 238000004458 analytical method Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 5
- 239000002131 composite material Substances 0.000 description 5
- 230000005236 sound signal Effects 0.000 description 5
- 238000011156 evaluation Methods 0.000 description 4
- 230000001965 increasing effect Effects 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000001105 regulatory effect Effects 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 230000006399 behavior Effects 0.000 description 2
- 230000001364 causal effect Effects 0.000 description 2
- 230000001934 delay Effects 0.000 description 2
- 238000012423 maintenance Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000013519 translation Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000021615 conjugation Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H21/00—Adaptive networks
- H03H21/0012—Digital adaptive filters
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
- G10L19/00—Speech 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/02—Speech 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 spectral analysis, e.g. transform vocoders or subband vocoders
- G10L19/0204—Speech 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 spectral analysis, e.g. transform vocoders or subband vocoders using subband decomposition
Landscapes
- Filters That Use Time-Delay Elements (AREA)
Abstract
本发明公开了一种用于设计一套子带FIR滤波器的方法。其中,每个FIR滤波器具有许多的滤波器系数,并被连接到可调节延迟线的上游或下游。该方法包括:将输入信号划分成许多子带信号,这些子带信号的频谱组成输入信号的频谱;对于每个子带,提供各自的取决于期望信号的期望子带信号;使用相应的FIR滤波器和延迟线对每个子带信号进行滤波和延迟,从而生成滤波后的信号;对于每个子带,提供取决于滤波信号和相应的期望信号的误差信号;对每个子带FIR滤波器的滤波器系数进行调整,使得各自的滤波信号与相应的期望子带信号近似匹配;以及,对于每个子带,改变各自的延迟线的延迟,以便最小化或最大化第一质量准则。
Description
技术领域
本发明涉及在音频应用中的子带信号处理,尤其涉及一种用于设计一套子带FIR滤波器的方法,该方法提供了期望目标传递函数的计算机执行的高效的实施方式。
背景技术
专利公开EP 0 746 133 A2描述了包括时间变化补偿的多级回声消除器。该回声消除器包括两个分别包括自适应滤波器的级联回声消除器级(stage),从而在第二级中,可调节延迟元件被连接到自适应滤波器的上游。
与IIR滤波器(无限脉冲响应滤波器)相比,FIR滤波器(有限脉冲响应滤波器)提供了实现具有可任意定义的期望传递函数(即,幅度和相位响应)的数字滤波器的可能性。因此,该期望的传递函数不仅可被设计成具有最小相位,而且也可以具有最大相位或混合相位。进一步地,也可实现线性相位传递函数,这通常是在音频信号处理中所期望的。
然而,尤其在音频信号处理中,当使用FIR滤波器时,所需滤波器长度通常相当高,这增加了操作过程中的计算工作量以及存储需求。这主要归因于两个物理事实,第一,在大多数音频应用中存在的室内脉冲响应的衰减时间(decay time)相当长,这导致了在使用FIR滤波器时相应的滤波器长度很长。第二,人类的听觉***(human auditory system)适于在不同频带上提供不等的频率分辨率。人类的听觉***非常善于分辨低频率,即,在低的绝对频率处可完全识别频率差异,而对高频率则不是那么容易辨别。例如,100Hz的音调和200Hz的音调可轻易地被完全区分开,而人耳要将5000Hz的音调与5100Hz的音调区分开则是很困难的,尽管这两种情况下的频率差均为100Hz。也就是说,人类听觉***的频率分辨率是随着频率的增加而降低的。这种现象是众所周知的,并且形成了适用于人类听觉***的心理声学频率标度(scale)的基础,例如,Bark标度、Mel标度,以及ERB(等效矩形频宽)标度。
研究表明,取决于内部(地毯、布置的家具等),收听室的室内脉冲响应相当长,尤其是在低频率处,因为能量的递降很慢。这种效果甚至被增强,因为音频复制***生成的声压处于低音频率(即,低于200Hz)范围的最大处,然而人类听觉***对于低频率音频信号较为不敏感。
所有这些因素的综合观点指示,人类听觉***的特性以及典型的收听室的特性要求FIR滤波器的长度不能小于特定的最小长度,以便在音频信号处理***中提供足够好的音频质量。例如,为在低音频率范围内提供约10Hz的必需的频率分辨率,在44100HZ的采样频率下,每个音频信道需要具有4410个滤波器系数的FIR滤波器,由此现代音频***具有多达16个信道。
只要FIR滤波器在操作过程中带来相当大的计算工作量和/或相当高的存储需求,就存在着对这样一种用于在音频应用中高效实施FIR滤波器的需要,以便能够使用标准数字音频信号处理硬件,并因此造成了对适当的滤波器设计方法的需要。
发明内容
本文公开了一种用于设计一套子带FIR滤波器的方法。其中,每个FIR滤波器具有许多滤波器系数,并且每个FIR滤波器被连接到可调节延迟线上游或下游。该方法包括:将输入信号划分成许多子带信号,这些子带信号的频谱组成输入信号的频谱;对于每个子带,提供各自的取决于期望信号的期望子带信号;将每个子带信号进行滤波和延迟,从而使用相应的FIR滤波器和延迟线生成滤波信号;对于每个子带,提供取决于滤波信号的误差信号和相应的期望信号;对每个子带FIR滤波器的滤波器系数进行调整,使得各自的滤波信号与相应的期望子带信号近似匹配;以及,对于每个子带改变各自延迟线的延迟,以便最小化或最大化第一质量准则。
附图说明
参考后面的附图和描述能够更好地理解本发明。在附图中,相同的参考数字指示相应的部件。在附图中:
图1示出了FIR滤波器的滤波器系数的自适应计算的基本信号处理结构,该FIR滤波器代表与预定义目标函数P(z)近似匹配的传递函数G(z);
图2示出了图1中的结构的调整,由此图1的全波段FIR滤波器G(z)可以用一套子带FIR滤波器Gm(z)代替;
图3示出了包括可替换全波段FIR滤波器G(z)的一套子带FIR滤波器Gm(z)的信号处理结构;
图4示出了包括许多子带的自适应FIR滤波器设计的信号处理结构,每个子带包括子带FIR滤波器和延迟线;
图5示出了当调整子带FIR滤波器的长度时,基于强调低频率子带的Bark标度的加权因子;以及
图6(a)和6(b)示意性地示出了偶数堆栈滤波器组和奇数堆栈滤波器组的幅度响应。
具体实施方式
图1示出了用于自适应计算滤波器20的滤波器系数gk(k=0,1,...,K-1)的基本信号处理结构,其中下标k表示滤波器系数的索引,而K表示滤波器长度。FIR滤波器20具有(离散)传递函数G(z),在成功调整滤波器系数gk之后,传递函数G(z)与参考***10的预定义的目标函数P(z)近似匹配。为了执行自适应滤波器的设计程序,为参考滤波器10和FIR滤波器20提供测试信号(输入信号x[n]),该测试信号为,例如,白噪声或具有包括目标传递函数P(z)的通带的带宽的任意其它信号。从参考***10的输出信号,即,从期望信号d[n],中减去(减法器30)FIR滤波器20的输出信号y[n]。差值d[n]-y[n]被用作误差信号e[n]并被提供给调整单元21。调整单元21被配置成在每个样本时间段中根据误差信号和输入信号x[n](在该情况下还被表示为参考信号)来计算一组更新的FIR滤波器系数gk。对于滤波器系数的调整可采用最小均方(LMS)算法或归一化最小均方(NLMS)算法。然而,为了达到该目的也可利用不同的调整算法。在调整算法收敛之后,FIR滤波器系数gk代表传递函数G(z),传递函数G(z)是目标传递函数P(z)的最佳近似。
当使用FIR滤波器时,减少计算工作量的一种选择是将要滤波的信号的频谱划分为许多窄带信号(子带信号),并对每个窄带信号进行独立的滤波。将全波段信号划分成若干子带信号可通过所谓的分析滤波器组(AFB)来实现。类似地,子带信号可被(再次)组合成具有相应的合成滤波器组(SFB)的单个全波段信号。在下面,全波段信号被表示为没有下标,例如,期望信号d[n],其中n是时间索引。进一步地,具有下标的信号,例如dm[n],表示的是一套子带信号,这些子带信号是相应的全波段信号d[n]的分解。因此,下标m表示子带序号(m=1,2,...,M)。相似地,离散全波段传递函数G(z)可被分解成许多个子带传递函数Gm(z)。
图2示出了类似于图1的基本信号处理结构,由此自适应FIR滤波器20被一套子带FIR滤波器20’代替。为了该目的,通过使用AFB 22,将全波段输入信号x[n]划分成M个子带输入信号xm[n](m=1,2,...,M)。相似地,使用AFB11,将全波段期望信号d[n]分成M个子带信号dm[n](同样,m=1,2,...,M)。实现窄带传递函数Gm(z)的每个子带FIR滤波器,下标m同样表示子带的序号。每个子带滤波器Gm(z)还可以用其滤波器系数gmk来代表,由此,同样,k表示范围从k=0到k=Km-1(Km是滤波器Gm(z)在第m个子带中的滤波器长度)的滤波器系数的索引。每个FIR滤波器Gm(z)与调整单元(该套调整单元在图2中用附图标记21’表示)相关联,该调整单元21’接收相应的误差信号em[n]=dm[n]-ym[n],并且为各自的子带m计算各自的一套更新的滤波器系数gmk(k=1,2,...,km-1)。
类似于图1中的例子,对于M个子带FIR滤波器Gm(z)中的每一个的滤波器系数gmk进行调整,使得,在调整算法收敛之后,由所有子带传递函数Gm(z)组合获得的总传递特性与预定义的目标函数P(z)匹配。
在计算适合的滤波器系数gmk之后,该套FIR滤波器Gm(z)20’可以在分析滤波器组(AFB 22)和相应的合成滤波器组(SFB 22’)之间工作,如图3所示,用于对音频信号进行滤波。在这种情况下,AFB 22、FIR滤波器组20’和SFB 22’一起实现与目标函数P(z)近似匹配的传递函数G(z),该传递函数可代表,例如,在音频***中的均衡滤波器。在现代音频***中,不仅幅度,而且相位也是均衡的对象,以便给听众生成期望的声音印象。因此,目标函数P(z)通常代表具有非线性相位特性的非最小相位滤波器。
图3中的包括分析滤波器组、一套子带滤波器和合成滤波器组的信号处理结构是否比具有相同传递函数的“普通”FIR滤波器更加高效(在计算工作量和存储器要求方面)呢?除其它因素之外,该问题的答案取决于分析和合成滤波器组的高效实施的有效性。为了解决人类听觉***的频率分辨率不均匀的问题,需要使用这样的滤波器组,其中,具有低中心频率的子带的带宽比具有较高中心频率的子带的带宽要窄。现有若干种方法用于实现这种心理声学推动的将全波段信号划分成一套子带信号,这些子带信号的带宽取决于相应的子带在可听得见的频率范围内的位置。然而,并没有已知的高效滤波器组允许将输入频谱不均匀地划分成一套不同带宽的子带。然而,本领域公知的其它方法允许将全波段信号划分成一套相等带宽的子带信号。一个例子是如S.Weiss等人描述的过取样(oversampled)GDFT滤波器组的快速实现(见S.Weiss,R.W.Stewart的“Fast Implementation of OversampledModulated Filter Banks”,IEE Electronics Letters,vol.36,pp.1502-1503,2000),其中使用了单个的原型滤波器(prototype filter)和在几乎每种信号处理环境下都可使用的FFT算法。
综上所述,滤波器组将被用于操作相等带宽的子带,因为没有用于处理非均匀带宽的子带的高效实现。然而,为了减轻相等宽度子带的不足,可以选择分配给各自的子带的不同的FIR滤波器的滤波器长度。即,在需要低频率分辨率的子带中的FIR滤波器与需要高频率分辨率的子带中的相比,包括较少的滤波器系数。后者的子带通常是那些位于可听得见频率范围的较低部分的子带。这样,可通过使用具有相等宽度的子带的高效滤波器组,来实现与人类听觉***的频率分辨率相对应的频率分辨率。
如上所述,目标函数P(z)一般为在频率上具有非线性组延迟特性的非最小相位滤波器。为了补偿由于在不同子带中的不同组延迟引起的不同的信号传播延迟,延迟线可被连接到每个子带FIR滤波器的上游或下游。因此,该延迟均衡不需要使用额外的FIR滤波器系数来获得,避免了计算效率低下。由于延迟的值以及滤波器系数的数量取决于要实现的目标传递函数P(z)(即,幅度和相位相应),因此对于每个子带可自适应地确定滤波器系数的数量(即,滤波器“接头(taps)”的数量)以及延迟值,这在下文中被描述为“自适应接头分配”算法和“自适应延迟分配”算法。因此,当设计M个子带FIR滤波器Gm(z)时,不仅自适应地确定滤波器系数(见图1中的系数gmk),而且还有系数KM的数量和额外延迟Δm。图4示出了子带FIR滤波器设计使用的信号处理结构。图4的例子是图2的结构的增强版本,在每个子带信号路径中具有额外延迟和全局“自适应接头分配和延迟分配单元”。
如在图2的例子中,将全波段输入信号x[n](例如,有限波段的白噪声)提供给具有目标传递函数P(z)的***10,从而生成期望信号d[n]。该期望信号d[n],以及输入信号x[n],被分别划分成M个子带信号dm[n]和xm[n]。图4的例子只示出了与第一个和最后一个(m=1和m=M)子带相关联的分量和信号。子带输入信号xm[n]被提供给具有传递函数Gm(z)的自适应FIR滤波器,从而生成滤波后的子带信号ym[n]。从对应的期望信号dm[n]中减去每个滤波后的子带信号ym[n],得到每个子带的误差信号em[n]。调整单元被分配给每个FIR滤波器Gm(z),以优化滤波器系数gmk(即,滤波器的脉冲相应(gm0,gm1,...,gm(k-1))),其中通过最佳的一套滤波器系数gmk将各自的误差信号em[n]的范数(norm)(例如,功率)最小化。
提供延迟Δm的延迟线被连接到每个子带FIR滤波器Gm(z)的上游或下游。进一步地,提供了“自适应接头分配和自适应延迟分配单元”40,其被配置成依据在下面更详细地描述的自适应接头分配和自适应延迟分配算法,动态地调整FIR滤波器Gm(k)的滤波器长度Km以及延迟线Δm的相应的延迟值。
对于以上提及的自适应接头分配(即,FIR滤波器长度的自适应)可以考虑不同的方法。一种方法是改变子带FIR滤波器Gm(z)的滤波器长度Km,直到总误差信号e[n](其中e[n]=e1[n]+e2[n]+...+eM[n])达到最小值为止。这样的方法在实践中产生了良好的结果,但是相当耗时,因为每次滤波器系数的数量改变之后,自适应滤波器都需要时间来再次收敛。产生与前面的方法差不多好的结果但是时间消耗少得多的另一种方法考虑最后S个滤波器系数gm(km-1),gm(km-2),...,gm(km-S)的能量。改变子带滤波器Gm(z)的滤波器长度Km,直到所提及的最后S个滤波器系数的能量近似相等为止。这种方法需要子带滤波器的脉冲响应随时间成指数衰减,在现实界***中应该一直是这样的情况。对比每个子带滤波器的最后S个滤波器接头的能量,以评定子带滤波器Gm(z)有多接近目标函数P(z),并提供在子带滤波器之间重新分配滤波器系数的准则,以便获得信号衰减行为与目标函数P(z)的脉冲响应的信号衰减行为相似的子带滤波器脉冲响应,其可以被看作最佳的最小化误差。
下面,更加详细地描述了自适应接头分配算法的例子。应该注意的是,对于实数值的全波段输入信号x[n](见图2),只有M/2子带需要被处理,因为另外的M/2子带信号是在第一个M/2子带上的信号的共轭复数复制。因此,各自的子带FIR滤波器传递函数遵守以下关系:
Gm(z)=GM-m+1(z)*,m=1,...,M/2, (1)
其中星号表示复共轭算子。
据此,子带滤波器Gm(z)的滤波器长度Km(其中m=1,2,...,M/2)随着Q采样(即,在子带***中的采样)的周期而被调整。但是,所有子带滤波器Gm(z)的滤波器系数gmk[n]的总数保持恒定。也就是说,如果一个或多个子带滤波器的滤波器长度增加,则另一滤波器的滤波器长度必须减少,以保持滤波器系数的总数恒定。据此,随着Q采样的周期,每个M/2子带FIR滤波器的长度减少ΔK个系数。因此,依据在以下进一步表述的一定的准则,有ΔK·M/2个“自由”系数将在M个子带滤波器之间被再次分配。
上述的“再次分配”可表达为以下等式:
其中m=1,2,...,M/2表示子带的数量。表达式Cm[n/Q]代表上述滤波器接头(即,滤波器系数)的分配准则。如上所述,一个有用的准则是子带误差信号em[n]的能量。在这种情况下Cm[n/Q]可被表示为:
其中m=1,2,...,M/2,并且R为采样数量,误差信号在该采样数量上被平均。但是,在对等式(3)求值之前,自适应子带FIR滤波器必须收敛,并且因此,R需要远小于Q(即,R<<Q)。另一个准则考虑了各自的子带FIR滤波器的最后S个滤波器系数的能量。在这种情况下,Cm[n/Q]可表示为:
其中m=1,2,...,M/2,并且Km为在各自子带中的当前滤波器长度Km[n]。可替代地,子带输入信号xm[n]的能量与最后S个滤波器系数可以一起被考虑(见等式5)。在这种情况下Cm[n/Q]可被表示成:
如上所述,依据等式(3)的准则,产生了最好的结果,但是求值是很耗时的。
当目标***是时变的并且输入信号是任意颜色时,例如在AEC(声音回波消除)***的情况下,应使用依据等式(5)的准则。在该情况下,即,对于子带FIR滤波器设计,输入信号x[n]可被设计者选择为白噪声,等式(4)得到高质量的结果并同时允许进行快速的自适应。
为了解决心理声学方面的问题,等式(3)、(4)或(5)中定义的表达式cm可以用相应的加权因子wm加权,即,在等式(3)、(4)或(5)中的表达式cm[n/Q]被替换成wm·cm[n/Q]。该加权因子wm的选择应该考虑人类听觉***的频率分辨率。使用Bark标度,因子wm可被计算如下:
其中fc,m表示第m个子带的中心频率(单位为Hz),其可以被计算成fc,m=(2m-1)·fs/(2·M),fs是单位为Hz的采样频率。在图5中示出了依据等式(6)计算的一套归一化的加权因子wm的例子,其中考虑了16个子带(M/2=16)。
如上所述,依据等式(2)的系数的再次分配可对子带误差信号em[n]的能量(对应等式(3))或者子带滤波器系数gmk(其中k={Km-S,...,Km-1},对应等式(4))的一部分(例如,最后S个滤波器系数)的能量,或者子带滤波器系数gmk的能量和子带输入信号xm[n]的能量的组合(对应等式(5))进行“均衡”。可以采用基于不同能量的准则。因此,该算法的“目标”是在各子带之间均等分配信号能量。由于该目标是优化处理的期望结果,所以上述算法也可被看作是最小化(或最大化)任务。
在上述准则的理想情况下,cm[n/Q](见等式(3)至(5))在每个子带m中具有相等的值,其中
也就是说,在M/2个被考虑的子带的每一个中,准则cm[n]是比在所有子带上求和的cm[n]总值小2/M的因子。返回上述的最小化任务,被最小化的质量函数可写成:
如上文所提及的,FIR滤波器Gm(z)在不同的子带内通常具有不同组延迟,导致不同的信号传播延迟。为了补偿这样的延迟差异,每个FIR滤波器被连接到每个子带内的可调节延迟线的上游或下游。再次,在该子带内的延迟值可如下面所描述的“自适应延迟分配”那样以迭代方式被优化。
如上面的“自适应接头分配”中所解释的,在子带中被自适应地分配给FIR滤波器的FIR滤波器系数可被认为是无限脉冲响应,该无限脉冲响应通过乘以矩形窗口函数而被截短,得到有限脉冲响应。对于每个子带FIR滤波器,这种矩形窗口函数可(沿着时间轴)被平移,其中各自的时间位移(timeshift)代表连接到相应的FIR滤波器的延迟线的有效延迟。
在此描述的自适应延迟分配算法的目标是,对于每个子带,找到导致各自的FIR滤波器的最大能量有限脉冲响应的延迟值(即,以上提到的矩形窗口的时间平移)。也就是说,当将延迟值分配到各自的延迟线时,只考虑能量或相应的FIR滤波器系数。但是,这并不必然是唯一的成功的方法。可替代地,可选择延迟值,使得总误差信号e[n]=e1[n]+e2[n]+...+eM/2[n]的范数被最小化。如关于自适应接头分配算法已经提到的那样,使用这样的准则将导致收敛很慢,同时会得到足够好的结果。可替代地,可选择延迟值,使得各自的FIR滤波器系数的总能量在最大值处。
但是,由于子带滤波器的FIR脉冲响应(即,FIR滤波器系数)通常包括与转换到子带有因果性和非因果性的分量,所以用于寻找最佳延迟值的其它准则的目标是均衡各自的脉冲响应的因果性和非因果性分量,即,使得各自的子带FIR脉冲响应的第一半的能量和第二半的能量之间的差异最小化。这样的优化策略可被简短地总结如下:
对于每个子带m,定义各自的子带FIR脉冲响应gmk的第一部分和第二部分,并计算该FIR脉冲响应的该第一部分和第二部分的能量。以V个采样为周期重复以下步骤:
(A)如果FIR脉冲响应gmk的第一部分的能量大于该脉冲响应第二部分的能量,则将各自延迟线的延迟Δm增加P个接头(例如,P=1)。
(B)如果FIR脉冲响应gmk的第一部分的能量低于该脉冲响应的第二部分的能量,则将各自延迟线的延迟Δm减少P个接头(例如,P=1)。但最小延迟为零。
(C)如果两部分能量相等,则保持延迟不变。
人们可能会试图将FIR脉冲响应划分成相等长度的两部分。但这只会对脉冲响应的最大值近似位于脉冲响应的中间的脉冲响应产生良好的结果。考虑到最大值在左侧的因果性脉冲响应,以上的优化策略将使上述的矩形窗口在时间上向右平移,直到能量被均等分配为止,这样,脉冲响应的左侧被***零,这显然是不期望的。因此,如果上述的FIR脉冲响应的第一和第二部分未覆盖各自脉冲响应的全部长度,而只是覆盖了由变量U定义的上部和下部边缘,此处0<U<1,则实际上获得了更好的结果。作为一个例子,如果U=0.6,则计算能量时,只考虑滤波器系数的60%,其中下部的30%形成以上提到的第一部分,而上部的30%对应于第二部分。这样,在第m个子带中,第一部分和第二部分的能量EmL和EmU分别为:
这样,如果EmU<EmL,则延迟值Δm被增加,如果EmU>EmL,则延迟值Δm被减少,而如果EmU=EmL则延迟值Δm保持不改变,或者表达成数学形式:
自适应延迟分配算法通常也可被看作最小化任务,其中,在每个子带中,改变各自的延迟值Δm,直到质量准则C′m[n/V]=|EmU-EmL|达到最小值为止,由此索引n/V表示该准则只在每V个采样时间上被求值。
应该注意的是,当在一个子带中将延迟线的延迟值Δm调节(例如,增加)一定数量的接头时,相应的滤波器脉冲响应(FIR滤波器系数)被向左平移,从而保持总延迟(包括延迟线的延迟和FIR滤波器的延迟)恒定。这避免了由于延迟值Δm的变化而对FIR滤波器系数的漫长的再次自适应的需要。然而,执行自适应延迟分配和自适应接头分配的自适应周期应当比LMS/NLMS算法的自适应步长要长得多,以使得当自适应接头分配和自适应延迟分配开始时,FIR滤波器系数gmk的自适应处于稳定状态。
虽然在每个子带中的滤波器长度的自适应调节和额外的延迟的自适应调节可被有益地组合,但这两种调整方法(自适应接头分配和自适应延迟分配)可彼此独立地被单独使用。例如,对于每个子带(例如,依照Bark标度)可预定义滤波器系数的数量,并且延迟Δm可以如上面解释的那样被自适应地调节。进一步,滤波器长度的自适应确定可依照上面解释的自适应接头分配方法来执行,而没有单独的延迟单元或带有固定的预定义延迟。这尤其可适于期望目标函数P(z)为线性相位或最小相位传递函数的情况。
正如在上面已经提到的,分析滤波器组以及相应的合成滤波器组(对应于图3中的滤波器组22和22’)将实数值的全波段输入信号x[n]划分成一套复数值的子带信号xm[n],并且将复数值的滤波后的子带信号ym[n]再次组合成实数值的输出信号y[n]。然而,也可以使用只操作实数值的子带信号的其它滤波器组。在这种情况下,子带FIR滤波器可更高效地被实施,然而,实数值的滤波器组需要相应复数值滤波器组的两倍的MAC(乘-累加操作)。进一步地,可在使用复数子带信号时通过因子M对子带进行下采样(under-sampled),相反,当使用实数值的滤波器组时,最大下采样因子仅为M/2。最后,所需要的子带FIR滤波器的滤波器长度仅为使用实数值的子带信号时所需的长度的一半。这表示复数值AFB和SFB的使用比实数值的滤波器组需要的计算工作量更少。但是,两种类型的滤波器组都可被应用于本发明的使用中。
当使用复数值滤波器组时,可在偶数堆栈或奇数堆栈滤波器组之间进行选择。这两种类型都可以通过通用离散傅立叶变换算法(GDFT算法)来高效地实施。图6示出了偶数堆栈(图6a)以及奇数堆栈(图6b)滤波器组的子带。乍看来,使用奇数堆栈GDFT滤波器组似乎更加优选,因为M个子带的一半是该M个子带的另一半的复共轭复制,并且因此仅需要处理数量为M/2个复数值的子带信号。当使用偶数堆栈的GDFT滤波器组时,最下部的和最上部的子带为实数值的,但是只具有其它子带一半的带宽。因此,必须处理M/2+1个子带信号,其中最下部子带中的信号为实数值的。结果,与复数子带相比,在最下部子带中仅需要一半的存储器以及四分之一的MAC。考虑到该事实,最下部子带的FIR滤波器(G1(z))具有最大的滤波器系数g1k(当考虑例如依据Bark或任意其它心理声学调整的频率标度的频率分辨率时),在存储器需求方面的总体效率提高了25%,在MAC需求方面提高了50%,尽管要处理的子带的总数量从M/2增加到M/2+1。
虽然已经详细地描述了本发明及其优点,但应该理解,可以做出各种变化、置换和替代,而不会背离由所附的权利要求书所限定的本发明的精神和范围。
而且,本申请的保护范围并不是要限于在本说明书中描述的处理、机械、制造、物质合成、装置、方法和步骤的特定实施例。本领域普通技术人员从本发明的揭示内容将容易想到,根据本发明,可以采用当前存在的或稍后将被开发的,足以执行与在此描述的相应的实施例相同的功能,或足以获得与在此描述的相应的实施例相同的结果的处理、机械、生产、物质合成、装置、方法或步骤。因此,所附的权利要求旨在将这样的处理、机械、生产、物质合成、装置、方法或步骤包括在其保护范围以内。
Claims (12)
1.一种用于设计一套子带FIR滤波器(Gm(z))的方法,每个FIR滤波器(Gm(z))具有许多(Km)个滤波器系数(gmk),并被连接到可调节延迟线的上游或下游;所述方法包括:
将输入信号(x[n])划分成许多(M)个子带信号(xm[n]),所述子带信号的频谱组成所述输入信号的频谱;
对于每个子带,提供各自的期望子带信号(dm[Z]),该期望子带信号(dm[Z])取决于期望信号(d[n]);
分别使用相应的FIR滤波器(Gm(z))和延迟线对每个子带信号(xi[n])进行滤波和延迟,从而生成滤波后的信号(ym[n]);
对于每个子带,提供取决于所述滤波后的信号(ym[n])和相应的期望信号(dm[n])的误差信号(em[n]);
调整每个子带FIR滤波器(Gm(z))的滤波器系数(gmk),使得各自的滤波后的信号(ym[n])与相应的期望子带信号(dm[n])近似匹配;以及
对于每个子带,改变各自的延迟线的延迟(Δm),以便最小化或最大化第一质量准则(cm’[n/V])。
2.如权利要求1所述的方法,其中,依照心理声学推动的频率标度,尤其依照Bark标度、Mel标度或ERB标度,所述子带FIR滤波器(Gm(z))的脉冲响应彼此不同。
3.如权利要求1或2所述的方法,进一步包括:
对于每个子带,改变各自的FIR滤波器(Gm(z))的滤波器系数(gmk)的数量(Km),以便最小化或最大化第二质量准则(cm(n/Q))。
4.如权利要求1至3中的任一项所述的方法,其中,对于每个子带,所述第一质量准则(cm’[n/V])取决于各自的子带FIR滤波器Gm(z)的滤波器系数的能量。
5.如权利要求4所述的方法,进一步包括:
对于每个子带,定义各自的子带FIR滤波器的滤波器系数的第一子集;
对于每个子带,定义各自的子带FIR滤波器的滤波器系数的第二子集;
如果,在各自的子带中,所述第一子集的滤波器系数的能量与所述第二子集的滤波器系数的能量相等,则认为所述第一质量准则(cm’[n/V])最小。
6.如权利要求1至3中的任一项所述的方法,进一步包括:
计算每个子带的误差信号的总和的范数(‖e[n]‖=‖e1[n]+e2[n]+...+eM/2[n]‖),
如果所述范数变为最小,则认为所述第一质量准则(cm’[n/V])最小。
7.如权利要求1至6中的任一项所述的方法,其中,所述将输入信号划分成许多个子带信号的步骤包括,使用偶数堆栈概括的DFT滤波器组对所述输入信号进行滤波。
8.如权利要求7所述的方法,其中所述子带输入信号(xm[n])和相应的FIR滤波器系数在最低和最高子带(M=1)为实数值,并且在其它子带(M=2,...,M/2+1)为复数值。
9.如权利要求3所述的方法,其中所述第二质量准则(cm(n/Q))取决于所述子带误差信号的能量或所述子带输入信号的能量或所述子带FIR滤波器系数的能量或以上的组合。
10.如权利要求9所述的方法,其中所述第二质量准则(cm(n/Q))取决于包括最后一个滤波器系数的连续的滤波器系数的子集的能量。
11.如权利要求9所述的方法,其中所述第二质量准则(cm(n/Q))取决于子带误差信号的能量(em[n],em[n-1],......,em[n-R]),只考虑包括最后一个采样的各自的误差信号的给定数量的连续采样。
12.如权利要求1至11中的任一项所述的方法,其中所述第二质量准则(cm(n/Q))在每个子带中被加权,其中所述加权因子代表人类听觉***的频率分辨率。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP10155192.7 | 2010-03-02 | ||
EP10155192.7A EP2365630B1 (en) | 2010-03-02 | 2010-03-02 | Efficient sub-band adaptive fir-filtering |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102195616A true CN102195616A (zh) | 2011-09-21 |
CN102195616B CN102195616B (zh) | 2014-12-17 |
Family
ID=42705160
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110050036.7A Active CN102195616B (zh) | 2010-03-02 | 2011-03-02 | 高效子带自适应有限脉冲响应滤波 |
Country Status (3)
Country | Link |
---|---|
US (1) | US8788554B2 (zh) |
EP (1) | EP2365630B1 (zh) |
CN (1) | CN102195616B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104579243A (zh) * | 2015-02-11 | 2015-04-29 | 桂林电子科技大学 | 一种dft调制滤波器组的快速设计方法 |
CN105425582A (zh) * | 2015-11-04 | 2016-03-23 | 北京航空航天大学 | 一种基于卡尔曼滤波的Stewart机构在线标定方法 |
CN106170049A (zh) * | 2016-05-12 | 2016-11-30 | 西南交通大学 | 一种有偏补偿的归一化子带自适应回声消除方法 |
CN107277669A (zh) * | 2017-07-31 | 2017-10-20 | 歌尔科技有限公司 | 耳机的数字降噪滤波器生成方法及装置 |
CN107483029A (zh) * | 2017-07-28 | 2017-12-15 | 广州多益网络股份有限公司 | 一种自适应滤波器的长度调节方法及装置 |
CN110266287A (zh) * | 2019-05-05 | 2019-09-20 | 深圳信息职业技术学院 | 电子耳蜗分数延迟滤波器构造方法、存储介质及电子耳蜗 |
CN113676156A (zh) * | 2021-08-09 | 2021-11-19 | 成都玖锦科技有限公司 | 一种基于lms的任意幅频响应fir滤波器设计方法 |
CN113711304A (zh) * | 2019-03-29 | 2021-11-26 | 伯斯有限公司 | 用于具有部分非因果传递函数的***的子带自适应滤波器 |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2488599B (en) * | 2011-03-04 | 2017-11-29 | Snell Advanced Media Ltd | Adaptive signal processing |
WO2015041477A1 (ko) * | 2013-09-17 | 2015-03-26 | 주식회사 윌러스표준기술연구소 | 오디오 신호 처리 방법 및 장치 |
US10204630B2 (en) | 2013-10-22 | 2019-02-12 | Electronics And Telecommunications Research Instit Ute | Method for generating filter for audio signal and parameterizing device therefor |
BR112016014892B1 (pt) | 2013-12-23 | 2022-05-03 | Gcoa Co., Ltd. | Método e aparelho para processamento de sinal de áudio |
EP4294055A1 (en) * | 2014-03-19 | 2023-12-20 | Wilus Institute of Standards and Technology Inc. | Audio signal processing method and apparatus |
CN108966111B (zh) | 2014-04-02 | 2021-10-26 | 韦勒斯标准与技术协会公司 | 音频信号处理方法和装置 |
US9704471B1 (en) * | 2016-03-30 | 2017-07-11 | Bose Corporation | Adaptive modeling of secondary path in an active noise control system |
GB201709849D0 (en) * | 2017-06-20 | 2017-08-02 | Nokia Technologies Oy | Processing audio signals |
CN111277245B (zh) * | 2020-03-24 | 2023-03-10 | 西安电子科技大学 | 一种滤波正交频分复用***的低阶子带滤波器设计方法 |
US11610597B2 (en) * | 2020-05-29 | 2023-03-21 | Shure Acquisition Holdings, Inc. | Anti-causal filter for audio signal processing |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5272656A (en) * | 1990-09-21 | 1993-12-21 | Cambridge Signal Technologies, Inc. | System and method of producing adaptive FIR digital filter with non-linear frequency resolution |
EP0746133A2 (en) * | 1995-05-31 | 1996-12-04 | AT&T Corp. | Multistage echo canceller including time variation compensation |
CN1790902A (zh) * | 2004-12-13 | 2006-06-21 | 上海无线通信研究中心 | 一种自适应滤波方法及装置 |
CN101512899A (zh) * | 2006-07-04 | 2009-08-19 | 杜比瑞典公司 | 滤波器压缩器以及用于产生压缩子带滤波器冲激响应的方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5553014A (en) * | 1994-10-31 | 1996-09-03 | Lucent Technologies Inc. | Adaptive finite impulse response filtering method and apparatus |
US5956674A (en) * | 1995-12-01 | 1999-09-21 | Digital Theater Systems, Inc. | Multi-channel predictive subband audio coder using psychoacoustic adaptive bit allocation in frequency, time and over the multiple channels |
US6915264B2 (en) * | 2001-02-22 | 2005-07-05 | Lucent Technologies Inc. | Cochlear filter bank structure for determining masked thresholds for use in perceptual audio coding |
EP2200180B1 (en) * | 2008-12-08 | 2015-09-23 | Harman Becker Automotive Systems GmbH | Subband signal processing |
-
2010
- 2010-03-02 EP EP10155192.7A patent/EP2365630B1/en active Active
-
2011
- 2011-03-02 US US13/038,964 patent/US8788554B2/en active Active
- 2011-03-02 CN CN201110050036.7A patent/CN102195616B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5272656A (en) * | 1990-09-21 | 1993-12-21 | Cambridge Signal Technologies, Inc. | System and method of producing adaptive FIR digital filter with non-linear frequency resolution |
EP0746133A2 (en) * | 1995-05-31 | 1996-12-04 | AT&T Corp. | Multistage echo canceller including time variation compensation |
CN1790902A (zh) * | 2004-12-13 | 2006-06-21 | 上海无线通信研究中心 | 一种自适应滤波方法及装置 |
CN101512899A (zh) * | 2006-07-04 | 2009-08-19 | 杜比瑞典公司 | 滤波器压缩器以及用于产生压缩子带滤波器冲激响应的方法 |
Non-Patent Citations (1)
Title |
---|
S WEISS, RW STEWART: "Fast implementation of oversampled modulated filter banks", 《IEEE ELECTRONICS LETTERS》 * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104579243A (zh) * | 2015-02-11 | 2015-04-29 | 桂林电子科技大学 | 一种dft调制滤波器组的快速设计方法 |
CN104579243B (zh) * | 2015-02-11 | 2017-10-27 | 桂林电子科技大学 | 一种dft调制滤波器组的快速设计方法 |
CN105425582A (zh) * | 2015-11-04 | 2016-03-23 | 北京航空航天大学 | 一种基于卡尔曼滤波的Stewart机构在线标定方法 |
CN105425582B (zh) * | 2015-11-04 | 2018-03-13 | 北京航空航天大学 | 一种基于卡尔曼滤波的Stewart机构在线标定方法 |
CN106170049A (zh) * | 2016-05-12 | 2016-11-30 | 西南交通大学 | 一种有偏补偿的归一化子带自适应回声消除方法 |
CN106170049B (zh) * | 2016-05-12 | 2019-05-17 | 西南交通大学 | 一种有偏补偿的归一化子带自适应回声消除方法 |
CN107483029B (zh) * | 2017-07-28 | 2021-12-07 | 广州多益网络股份有限公司 | 一种voip通讯中的自适应滤波器的长度调节方法及装置 |
CN107483029A (zh) * | 2017-07-28 | 2017-12-15 | 广州多益网络股份有限公司 | 一种自适应滤波器的长度调节方法及装置 |
CN107277669A (zh) * | 2017-07-31 | 2017-10-20 | 歌尔科技有限公司 | 耳机的数字降噪滤波器生成方法及装置 |
CN113711304A (zh) * | 2019-03-29 | 2021-11-26 | 伯斯有限公司 | 用于具有部分非因果传递函数的***的子带自适应滤波器 |
US11770114B2 (en) | 2019-03-29 | 2023-09-26 | Bose Corporation | Subband adaptive filter for systems with partially acausal transfer functions |
CN113711304B (zh) * | 2019-03-29 | 2023-12-08 | 伯斯有限公司 | 用于具有部分非因果传递函数的***的子带自适应滤波器 |
CN110266287A (zh) * | 2019-05-05 | 2019-09-20 | 深圳信息职业技术学院 | 电子耳蜗分数延迟滤波器构造方法、存储介质及电子耳蜗 |
CN113676156A (zh) * | 2021-08-09 | 2021-11-19 | 成都玖锦科技有限公司 | 一种基于lms的任意幅频响应fir滤波器设计方法 |
CN113676156B (zh) * | 2021-08-09 | 2024-01-26 | 成都玖锦科技有限公司 | 一种基于lms的任意幅频响应fir滤波器设计方法 |
Also Published As
Publication number | Publication date |
---|---|
EP2365630B1 (en) | 2016-06-08 |
US20110216807A1 (en) | 2011-09-08 |
CN102195616B (zh) | 2014-12-17 |
EP2365630A1 (en) | 2011-09-14 |
US8788554B2 (en) | 2014-07-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102195616B (zh) | 高效子带自适应有限脉冲响应滤波 | |
CN102811037B (zh) | 高效子频带自适应有限脉冲响应滤波 | |
US7818079B2 (en) | Equalization based on digital signal processing in downsampled domains | |
CN100477704C (zh) | 用于与自适应波束形成组合的回声抵消的方法和设备 | |
KR100312636B1 (ko) | 보상필터 | |
EP0615350B1 (en) | Method of identifying an unknown system with a band-splitting adaptive filter and a device thereof | |
KR100584491B1 (ko) | 다수의 소스들을 갖는 오디오 처리 장치 | |
EP1495463B1 (en) | Active noise control system in unrestricted space | |
US20070253577A1 (en) | Equalizer bank with interference reduction | |
Ramos et al. | Filter design method for loudspeaker equalization based on IIR parametric filters | |
EP1919251B1 (en) | Beamforming weights conditioning for efficient implementations of broadband beamformers | |
US20080228470A1 (en) | Signal separating device, signal separating method, and computer program | |
US7203324B2 (en) | Audio feedback processing system | |
EP2392149A2 (en) | Method for determining inverse filter from critically banded impulse response data | |
JPH0879137A (ja) | 帯域分割適応フィルタによる未知システム同定の方法及び装置 | |
US9997171B2 (en) | Multi-band signal processor for digital audio signals | |
CN116709117A (zh) | 用于扬声器群时延补偿的全通滤波器及其设计方法 | |
CN101401454A (zh) | 立体声成像 | |
Sebastian et al. | A low complex 10-band non-uniform FIR digital filter bank using frequency response masking technique for hearing aid | |
Vary | An adaptive filter-bank equalizer for speech enhancement | |
Sebastian et al. | Digital filter bank for hearing aid application using FRM technique | |
Xiuli et al. | Design of digital parametric equalizer based on second-order function | |
EP2941020B1 (en) | A multi-band signal processor for digital audio signals | |
Malathi et al. | FPGA Implementation of Adaptive NMLS Algorithm: Timbre Based Filtering from Multiple Harmonics using FIR Filters | |
Liu et al. | A variable step-size pre-filter-bank adaptive algorithm |
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 |