CN101695148B - 一种用于数字助听器的多通道宽动态范围压缩系统 - Google Patents
一种用于数字助听器的多通道宽动态范围压缩系统 Download PDFInfo
- Publication number
- CN101695148B CN101695148B CN200910236221A CN200910236221A CN101695148B CN 101695148 B CN101695148 B CN 101695148B CN 200910236221 A CN200910236221 A CN 200910236221A CN 200910236221 A CN200910236221 A CN 200910236221A CN 101695148 B CN101695148 B CN 101695148B
- Authority
- CN
- China
- Prior art keywords
- pass
- mrow
- channel
- transformation
- filter bank
- 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.)
- Expired - Fee Related
Links
- 230000009466 transformation Effects 0.000 claims abstract description 58
- 230000006835 compression Effects 0.000 claims abstract description 43
- 238000007906 compression Methods 0.000 claims abstract description 43
- 238000004458 analytical method Methods 0.000 claims abstract description 37
- 230000008447 perception Effects 0.000 claims abstract description 34
- 230000003321 amplification Effects 0.000 claims abstract description 19
- 238000003199 nucleic acid amplification method Methods 0.000 claims abstract description 19
- 238000000034 method Methods 0.000 claims description 28
- 238000003786 synthesis reaction Methods 0.000 claims description 16
- 238000004364 calculation method Methods 0.000 claims description 15
- 238000006243 chemical reaction Methods 0.000 claims description 14
- 230000008569 process Effects 0.000 claims description 14
- 230000015572 biosynthetic process Effects 0.000 claims description 13
- 238000012545 processing Methods 0.000 claims description 11
- 238000005070 sampling Methods 0.000 claims description 11
- 238000001514 detection method Methods 0.000 claims description 10
- 230000004044 response Effects 0.000 claims description 6
- 210000005069 ears Anatomy 0.000 claims description 4
- 239000012634 fragment Substances 0.000 claims description 4
- 230000001502 supplementing effect Effects 0.000 claims description 3
- 108010076504 Protein Sorting Signals Proteins 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 7
- 238000004088 simulation Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 208000032041 Hearing impaired Diseases 0.000 description 3
- 230000005236 sound signal Effects 0.000 description 3
- 206010011878 Deafness Diseases 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 231100000888 hearing loss Toxicity 0.000 description 2
- 230000010370 hearing loss Effects 0.000 description 2
- 208000016354 hearing loss disease Diseases 0.000 description 2
- 230000002194 synthesizing effect Effects 0.000 description 2
- 206010011891 Deafness neurosensory Diseases 0.000 description 1
- 208000009966 Sensorineural Hearing Loss Diseases 0.000 description 1
- 230000001364 causal effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 231100000879 sensorineural hearing loss Toxicity 0.000 description 1
- 208000023573 sensorineural hearing loss disease Diseases 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Landscapes
- Compression, Expansion, Code Conversion, And Decoders (AREA)
Abstract
本发明涉及一种基于听觉感知模型的多通道宽动态范围压缩系统,包括:模拟听觉感知模型的分析滤波器组、声压级检测模块、压缩放大增益计算模块、乘法器和模拟听觉感知模型的综合滤波器组;音频数字信号x(n)经分析滤波器组后分成K个通道,声压级检测模块检测各个通道的声压级,压缩放大增益计算模块计算得到各个通道的具体增益值,乘法器将通道增益值与对应子带信号相乘,所得结果通过综合滤波器组综合成一路输出信号y(n);该滤波器组通过加权叠接相加结构与全通变换相结合的方式,在加权叠接相加结构的分析滤波器组和综合滤波器组中分别进行全通变换和全通反变换,实现在通道数目较少的情况下模拟人耳听觉分辨率。
Description
技术领域
本发明涉及数字助听器中的数字信号处理技术,特别地,涉及一种用于数字助听器的多通道宽动态范围压缩系统,用于压缩放大数字助听器的音频信号。
背景技术
听力损失,尤其是感音神经性听力损失除了表现为听阈上升外,大部分患者的不舒适阈不变甚至下降,造成听觉宽动态范围减小。针对这一特点,数字助听器通常采用宽动态范围压缩(WDRC),将输入信号宽动态范围按一定的比例压缩到听障者残余的听力动态范围内。另外,听障者可能对某一频率范围的声音感知和普通人相同,但却无法和普通人一样听到其它频率的声压级较低的声音,解决这个问题的方法是进行多通道宽动态范围压缩。典型的多通道WDRC系统是先用滤波器组将信号划分为多个通道,对每个通道的信号单独进行压缩处理,这样可以根据该频率对应的听力损失情况设计合适的压缩比和压缩阈值,更灵活地进行听力补偿。最后,将处理后的各通道信号综合成一路信号。
因此,数字滤波器组的选择对于数字助听器的多通道宽动态范围压缩系统的信号质量、计算复杂度和信号延迟有重要影响。数字滤波器组按照带宽情况可以分为均匀滤波器组和非均匀滤波器组两种,在均匀滤波器组中,离散傅立叶变换(DFT)分析-综合滤波器组应用比较广泛,这种复调制滤波器组可以用多相(polyphase)结构或者加权叠接相加(weighted overlap-add,WOLA)结构高效实现,其中,WOLA结构通过改变DFT滤波器组升/降采样的位置,将卷积运算转化为加权运算且多通道共享加权过程,和用快速傅立叶变换(FFT)实现DFT,使得DFT滤波器组的实现效率大大提高。
然而,DFT滤波器组的频率分辨率不适合人类的听觉系统,因为人耳的频率分辨率是非均匀的,随着频率增大而降低,其实际感知频率尺度(听觉频率尺度)与普通频率尺度之间呈非线性关系,可以用Bark频率尺度表示。若采用均匀滤波器组,要使频率分辨率达到Bark频率尺度,要求通道划分很细,而这会造成高频部分的通道划分过细,浪费资源。因此,许多研究者提出在语音处理系统中用非均匀分析-综合滤波器组来处理音频信号,这种方法可行的原因有二:一是近似Bark频率尺度的非均匀滤波器组结合了听觉感知模型;二是,一般情况下,语音信号的能量主要集中在低频部分。
用非均匀滤波器组获得近似Bark频率尺度的方法之一是采用离散小波(包)变换,实现这一变换可以采用树形结构;另一种方法是采用非线性频率变换的分析-综合滤波器组,这种方法的信号延时和计算复杂度往往比树形结构滤波器组低。
设计非线性频率变换的分析-综合滤波器组通常采用一阶全通变换(allpasstransformation)法,通过调节一阶全通滤波器 的极点b的取值,可以很好地模拟Bark频率尺度变换,这使得全通变换分析-综合滤波器组达到人耳频率分辨率所需的通道数目明显小于均匀滤波器组,所以更适合用于语音信号处理系统或类似的其它应用中。
目前,非线性频率变换的分析-综合滤波器组主要有两类实现方式,一是将整个信号进行全通变换,再用均匀滤波器组处理,这样,在变换后的频率尺度上做均匀通道划分,就相当于在原频率尺度上做非均匀通道划分,此类方法需要对整个信号做时间翻转运算,因此,不能用于实时信号处理;二是用全通滤波器A(z)代替多相结构或WOLA结构中的延迟单元z-1,这样做的缺点是全通滤波器会造成相位失真,补偿这些失真使得综合滤波器组的设计更加复杂,同时会增加计算复杂度和信号延时。
多通道WDRC系统设计中的一个关键问题是如何使系统的频率分辨率与人类听觉系统的频率分辨率相匹配,又可以尽量减少通道数目。针对以离散傅立叶变换(DFT)为基础的多通道WDRC系统只能提供均匀通道划分的问题,专利号为US7277554的发明提出了用一阶全通滤波器A(z)代替离散傅立叶变换中的延迟单元z-1,从而实现在非线性频率变换域上计算压缩放大增益的宽动态范围压缩系统(Warped WDRC),适当选择全通滤波器A(z)的极点位置,该系统可以较好地模拟听觉系统。但是,由于全通滤波器的群延迟与频率有关,该系统会造成信号的不同频率成分通过该系统的时间不同,导致信号失真,这种由于频率相关的群延迟产生的信号失真在通道数目较多时尤为严重。
发明内容
本发明的目的是,针对上述几种方法存在的问题,提供一种用于数字助听器的多通道宽动态范围压缩系统,其中使用的模拟听觉感知模型的滤波器组(以下简称Warped WOLA滤波器组)将均匀DFT滤波器组的高效实现方法加权叠接相加(weighted overlap-add,WOLA)结构与全通变换相结合,以用较少的通道数目的滤波器组对音频信号实现近似人耳听觉分辨率的非均匀通道划分,同时,保证信号重建效果,并使其具有较低的计算复杂度。模拟听觉感知模型的滤波器组只对有限长信号片段做全通变换,所以可以应用于实时信号处理。而且,模拟听觉感知模型的滤波器组因为在综合过程中加入了全通反变换操作,所以基本不会产生由于全通变换引起的相位失真。另外,模拟听觉感知模型的滤波器组可以用较少的通道数目模拟人耳听觉分辨率,无需均匀DFT滤波器组中先做较多数目的分通道,再将高频通道进行合并来模拟人耳听觉分辨率的过程。
为了实现上述目的,本发明所述的模拟听觉感知模型的滤波器组实现信号分通道和信号综合,该滤波器组采用一种用加权叠接相加结构与全通变换相结合的方式,在加权叠接相加结构的分析和综合过程中适当的位置分别加入全通变换和全通反变换,从而实现在通道数目较少的情况下模拟人耳听觉分辨率,同时保证信号重建效果,另外,采用加权叠接相加结构使本发明的一种模拟听觉感知模型的滤波器组具有较低的计算复杂度。
本发明涉及一种基于听觉感知模型的多通道宽动态范围压缩系统,包括:模拟听觉感知模型的分析滤波器组2、声压级检测模块3、压缩放大增益计算模块4和模拟听觉感知模型的综合滤波器组6;音频数字信号x(n)1经分析滤波器组2后分成K个通道,经声压级检测模块3检测各个通道的声压级、压缩放大增益计算模块4计算得到各个通道的具体增益值后,乘法器5将通道增益值与对应子带信号相乘,所得结果通过综合滤波器组6综合成一路输出信号y(n)7;该滤波器组通过加权叠接相加结构与全通变换相结合的方式,在加权叠接相加结构的分析滤波器组2和综合滤波器组6中分别进行全通变换和全通反变换,实现在通道数目较少的情况下模拟人耳听觉分辨率。所述的模拟听觉感知模型的分析滤波器组2通过全通变换,以及调整全通变换参数,得到非均匀通道划分,实现步骤包括:
1)截取信号的步骤,截取有限长信号片段,长度P的取值为不大于(1+|b|)L/(1-|b|)的最大整数,截取步长为D,其中L为分析原型滤波器的长度,b为全通变换参数,D为降采样率;
2)全通变换的步骤,对有限长信号片段进行全通变换;所述的全通变换为L-1阶、全通变换参数为b的全通变换,L-1为全通滤波器的个数;
3)加权处理的步骤,用分析原型滤波器的时间翻转h(-r)对全通变换的结果做加权处理;
4)时间叠加的步骤,将加权处理的结果平均分成L/K个长度为K的小片段,并将这些小片段相加,L为分析原型滤波器的长度;
5)离散傅立叶变换的步骤,进行K点离散傅立叶变换,K为通道数目;
上述技术方案中,所述的步骤2)的全通变换,步骤包括:
21)对于有限长信号s(n)13,先做时间翻转14得到s(N-n);
在单位圆上,有A(ejω)=ejθ(ω),其中
式中ω=2πf/fs,fs为采样率。
所述的离散傅立叶变换可以用快速傅立叶变换实现。
所述的截取信号片段的长度P与分析原型滤波器的长度L之间满足:P是不大于(1+|b|)L/(1-|b|)的最大整数。
所述的模拟听觉感知模型的综合滤波器组6通过全通反变换,将全通变换后的频率尺度再变换成原来的频率尺度,实现步骤包括:
2)离散傅立叶反变换的步骤,进行K点离散傅立叶反变换(IDFT);
3)将离散傅立叶反变换的结果复制L/K次,形成一个长度为L的序列;
4)用综合原型滤波器f(r)对上述长度为L的序列进行加权处理;
5)全通反变换的步骤,对加权处理的结果进行P-1阶参数为b的全通反变换;
6)叠加的步骤,将全通反变换的结果叠加至长度为P的输出缓存,输出缓存的左侧移出D个点作为输出,同时,输出缓存的右侧补充D个零,另外,将输出结果通过频率响应为一阶全通滤波器A(z)的群延迟的倒数的滤波器,得到最终的输出结果。
所述的离散傅立叶反变换可以用快速傅立叶反变换(IFFT)实现。
所述的步骤5)的P-1阶参数为b的全通反变换,可以通过P-1阶参数为-b的全通变换实现。
对于输入信号x(n)1为实信号的情况,所述的声压级检测模块3和压缩放大增益计算模块4仅需对第0至K/2通道信号进行处理。
本发明在分析过程中,先将输入信号段进行全通变换,再对变换后的信号用分析原型滤波器的时间翻转加权,时间叠加,离散傅立叶变换(可以采用快速傅立叶变换),最后乘以中心频率调整系数,得到各个通道信号。在分析过程中引入了全通变换,实现输入信号的频率的非线性变换,这样在变换后的频率尺度上的均匀通道划分相当于在原频率尺度上的非均匀通道划分。适当选择全通变换系数,可以模拟Bark频率尺度,O.Smith等在1999年给出了根据采样率选择适当全通变换参数的具体计算公式。
根据声压级检测模块的结果和各压缩放大增益计算参数,计算各个通道信号的压缩放大增益值,并将其与对应通道信号相乘。
在综合过程中,先将各通道信号乘以另一组中心频率调整系数,再做离散傅立叶反变换(可以用快速傅立叶反变换实现),将上述结果复制若干次后,用综合原型滤波器对其进行加权处理,再对加权结果进行全通反变换,将结果加在输出缓存中,输出缓存每次左侧移出D个点,作为输出。最后,将输出结果通过频率响应为一阶全通滤波器A(z)的群延迟的倒数的滤波器,得到最终的输出结果。在综合过程中,加入了全通反变换,将变换后的频率尺度,变换成原来的频率尺度,避免由分析过程中频率变换引起的综合信号的相位失真,从而保证好的信号重建效果。
本发明的优点在于,本发明中采用的模拟听觉感知模型的滤波器组结合了DFT滤波器组的高效实现方法加权叠接相加结构和全通变换,既具有加权叠接相加结构的高效率,同时克服了目前频率变换滤波器组方法中存在的无法实时实现和相位失真的问题,可以用较少的通道数目实现人耳的频率分辨率,通道数目的减少不仅可以降低计算复杂度,同时可以降低系统的群延迟。该系统的一个应用是数字助听器。
本发明设计的基于听觉感知模型的多通道宽动态范围压缩系统可以根据实际需求调节通道数目、降采样率、分析原型/综合原型滤波器的截止频率和长度、全通变换参数、各个通道压缩放大增益计算参数等参数。通过调节全通变换参数,可以实现在通道数目较少的情况下模拟听觉感知模型的信号非均匀分通道。
附图说明
图1是本发明的基于听觉感知模型的多通道宽动态范围压缩系统原理图;
图2a是全通变换参数b=0时16通道划分情况图;
图2b是全通变换参数b=0.4时16通道划分情况图;
图3是模拟听觉感知模型的分析滤波器组的流程图;
图4是全通变换原理图;
图5a是全通滤波器A(z)的相位响应图;
图5b是全通滤波器A(z)的群延迟图;
图6是模拟听觉感知模型的综合滤波器组的流程图。
图7是算法仿真所用听力图
图8a是输入信号波形图;
图8b是WOLA-64WDRC输出波形图;
图8c是Warped WOLA-16 WDRC输出波形图;
图9a是输入信号语谱图;
图9b是WOLA-64 WDRC输出语谱图;
图9c是Warped WOLA-16 WDRC输出语谱图。
具体实施方式
如图1所示,本发明公开的基于听觉感知模型的多通道宽动态范围压缩系统包括模拟听觉感知模型的分析滤波器组2,声压级检测模块3,压缩放大增益计算模块4,乘法器5和模拟听觉感知模型的综合滤波器组6五个主要部分。输入数字信号x(n)1,经模拟听觉感知模型的分析滤波器组2后分成K个通道,声压级检测模块3检测各个通道的声压级,压缩放大增益计算模块4计算得到各个通道的具体增益值,乘法器5将通道增益值与对应子带信号相乘,所得结果通过模拟听觉感知模型的综合滤波器组6综合成一路输出信号y(n)7。
如图3所示,为了提高效率,模拟听觉感知模型的分析滤波器组2采用了加权叠接相加结构中加入了全通变换操作,调整全通变换参数,可以得到非均匀通道划分。例如,对于采样率为8kHZ,通道数目为K=16,分析原型滤波器h(n)的截止频率为π/K,图2a为全通变换参数b=0时的通道划分情况,此时为均匀滤波器组。图2b为全通变换参数b=0.4时的通道划分情况,为非均匀滤波器组。这样,可以直接对各个通道信号计算压缩放大增益。而对于基于WOLA滤波器组的宽动态范围压缩系统来说,需要先将信号分成较多数目的通道,然后将高频通道合并,以模拟听觉感知模型,之后,才能对合并后的通道进行压缩放大增益计算。
模拟听觉感知模型的分析滤波器组2具体实现流程如图3所示,实现步骤如下:
(1)长度为P的输入缓存每次右侧移入x(n)的D个点,左侧移出D个点;
(2)计算输入缓存中信号的L-1阶参数为b的全通变换;
(3)用分析原型滤波器的时间翻转h(-r)对(2)的结果加权;
(4)将(3)的结果平均分成L/K个长度为K的小片段,并将这些小片段相加;
(5)计算(4)的结果的K点快速傅立叶变换;
其中,全通变换的实现流程如图4所示。对于有限长信号s(n)13,先做时间翻转14得s(N-n),然后经过全通滤波器A(z)15链,当n=N时取值,得到变换结果16。A(z)15的个数称为全通变换的阶数。其中,
为实的、稳定的、因果的一阶全通滤波器,其极点b称为全通变换参数。在单位圆上,有A(ejω)=ejθ(ω),其中
式中ω=2πf/fs,fs为采样率。
图4的过程实现了θ:ω→θ(ω)的非线性频率变换。由A(z)的性质可知,实现其反变换θ-1:θ(ω)→ω只要将A(z)的极点改为-b即可。
图5a显示了极点b取不同值时θ(ω)与ω之间的关系,可见当b≠0时,二者是非线性关系,对θ(ω)进行均匀通道划分,相当于对ω进行非均匀通道划分。图5b显示了极点b取不同值时A(z)的群延迟τb(ω)=(1-b2)/(1-2bcosω+b2)随ω和b变化的情况。通过对A(z)的群延迟的分析可知,当b>0时,频率越高群延迟越小,要使中包含s(k)中的所有ω=π的信息,要求序列s(k)的长度P与A(z)链长度L之间满足:L≥(1+|b|)P/(1-|b|)。而如果只要求中包含s(k)中的所有ω=0的信息,则只需满足P≥(1+|b|)L/(1-|b|)。
对于实信号,正、负频率成分是对称的,所以此时声压级检测模块(3)和压缩放大增益计算模块(4)实际上只需对第0至K/2通道信号进行处理。
如图6所示,模拟听觉感知模型的综合滤波器组6的具体步骤为:
(2)计算(1)的结果的K点快速傅立叶反变换;
(3)将(2)的结果复制L/K次,形成一个长度为L的序列;
(4)用综合原型滤波器f(r)对(3)的结果加权;
(5)计算(4)的结果的P-1阶参数为-b的全通变换(即P-1阶参数为b的全通反变换);
(6)将(5)的结果叠加至长度为P的输出缓存序列,然后将输出缓存的左侧移出D个点,作为该步骤的输出,最后将输出缓存的右侧补充D个零。
(7)将(6)的输出通过频率响应为1/τb(ω)的滤波器,得到最后的输出y(n)。
为了与基于WOLA滤波器组的宽动态范围压缩系统进行对比,我们对下面的例子进行了仿真实验。
设一个听障者的听力图如图7所示。
用64通道的WOLA滤波器组(记作WOLA-64)和16通道的Warped WOLA滤波器组(记作Warped WOLA-16)对信号进行分通道,具体参数如表1所示。
表1WOLA-64和Warped WOLA-16的参数情况
滤波器组参数 | WOLA-64 | Warped WOLA-16 |
通道数目K | 64 | 16 |
降采样率D | 16 | 4 |
原型分析/综合滤波器长度L | 256 | 64 |
原型分析滤波器截止频率 | π/64 | π/16 |
原型综合滤波器截止频率 | π/16 | π/11 |
全通变换参数b | 无 | 0.5657 |
对WOLA-64高频通道进行适当合并,得到这两个滤波器组对应于正频率的9个通道的中心频率和带宽情况如表2所示。
表2WOLA-64和Warped WOLA-16的中心频率和带情况
根据听阈图(图7)设计相应通道的WDRC参数如表3所示。
表3各通道的WDRC参数
通道序号 | 1.2 | 3 | 4 | 5,6 | 7,8 | 9 |
启动时间(ms) | 5 | 5 | 5 | 5 | 5 | 5 |
恢复时间(ms) | 70 | 70 | 70 | 70 | 70 | 70 |
下拐点(dB) | 50 | 50 | 50 | 50 | 50 | 50 |
上拐点(dB) | 100 | 100 | 100 | 100 | 100 | 100 |
50dB增益(dB) | 20 | 25 | 30 | 35 | 40 | 45 |
80dB增益(dB) | 10 | 10 | 15 | 17 | 22 | 25 |
压缩比 | 1.5∶1 | 2∶1 | 2∶1 | 2.5∶1 | 2.5∶1 | 3∶1 |
测试信号为女声“school”,采样率为16kHz,其波形如图8a所示,语谱图如图9a所示。Warped WOLA-16和WOLA-64 WDRC系统的输出波形分别为图8b和图8c。图9b和图9c分别为WOLA-64 WDRC系统输出和Warped WOLA-16 WDRC系统输出的语谱图。从上述图像的比较可以看出,Warped WOLA-16 WDRC系统和WOLA-64 WDRC系统具有相似的压缩放大效果。
考察数字助听器的信号处理算法的重要指标有频率分辨率、计算复杂度和系统群延迟。从上面的分析可以看出,WOLA-64 WDRC和Warped WOLA-16 WDRC系统在获得相似频率分辨率的前提下,后者具有通道数目更少的优势。另外,WOLA-64延迟为272个采样点,而Warped WOLA-16的延迟为234个采样点。这是因为通道数目的减少可以降低对分析/综合原型滤波器长度的要求,从而降低滤波器组的延迟。滤波器组的延迟是WDRC系统群延迟的重要组成部分,因此Warped WOLA-16WDRC系统具有更小的群延迟。
最后所应说明的是,以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,对本发明的技术方案进行修改或者等同替换,都不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (7)
1.一种用于数字助听器的多通道宽动态范围压缩系统,该多通道宽动态范围压缩系统基于一种模拟听觉感知模型的滤波器组,包括:模拟听觉感知模型的分析滤波器组(2)、声压级检测模块(3)、压缩放大增益计算模块(4)、乘法器(5)和模拟听觉感知模型的综合滤波器组(6);音频数字信号x(n)(1)经分析滤波器组(2)后分成K个通道,经声压级检测模块(3)检测各个通道的声压级、压缩放大增益计算模块(4)计算得到各个通道的具体增益值后,乘法器(5)将通道增益值与对应子带信号相乘,所得结果通过综合滤波器组(6)综合成一路输出信号y(n)(7);其特征在于,该模拟听觉感知模型的滤波器组通过加权叠接相加结构与全通变换相结合的方式,在加权叠接相加结构的分析滤波器组(2)和综合滤波器组(6)中分别进行全通变换和全通反变换,实现在通道数目较少的情况下模拟人耳听觉分辨率;
所述的模拟听觉感知模型的分析滤波器组(2)通过全通变换,以及调整全通变换参数,得到非均匀通道划分,实现步骤包括:
1)截取信号的步骤,截取有限长信号片段,长度P的取值为不大于(1+|b|)L/(1-|b|)的最大整数,截取步长为D,其中L为分析原型滤波器的长度,b为全通变换参数,D为降采样率;
2)全通变换的步骤,对有限长信号片段进行全通变换;所述的全通变换为L-1阶、全通变换参数为b的全通变换,L-1为全通滤波器的个数;
3)加权处理的步骤,用分析原型滤波器的时间翻转h(-r)对全通变换的结果做加权处理;
4)时间叠加的步骤,将加权处理的结果平均分成L/K个长度为K的小片段,并将这些小片段相加,L为分析原型滤波器的长度;
5)离散傅立叶变换的步骤,进行K点离散傅立叶变换,K为通道数目;
6)用中心频率调整系数exp(jmDθ-1(ωk))乘以离散傅立叶变换的结果,得到各个通道信号的序列其为第k个分量对应第k个通道m时刻的输出,其中ωk=2kπ/K,k=0,...,K-1;
所述的模拟听觉感知模型的综合滤波器组(6)通过全通反变换,将全通变换后的频率尺度再变换成原来的频率尺度,实现步骤包括:
8)离散傅立叶反变换的步骤,进行K点离散傅立叶反变换;
9)将离散傅立叶反变换的结果复制L/K次,形成一个长度为L的序列;
10)加权处理的步骤,用综合原型滤波器f(r)对3)中的结果进行加权处理;
11)全通反变换的步骤,对加权处理的结果进行P-1阶参数为b的全通反变换;
12)叠加的步骤,将全通反变换的结果叠加至长度为P的输出缓存,并将输出缓存的左侧移出D个点作为输出,同时,将输出缓存右侧补充D个零,另外,将输出结果通过频率响应为一阶全通滤波器A(z)的群延迟的倒数的滤波器,得到最终的输出结果。
3.根据权利要求1所述的多通道宽动态范围压缩系统,其特征在于,所述的离散傅立叶变换用快速傅立叶变换实现。
4.根据权利要求1所述的多通道宽动态范围压缩系统,其特征在于,所述的截取信号片段的长度P与分析原型滤波器的长度L之间满足:P是不大于(1+|b|)L/(1-|b|)的最大整数。
5.根据权利要求1所述的多通道宽动态范围压缩系统,其特征在于,所述的离散傅立叶反变换用快速傅立叶反变换实现。
6.根据权利要求1所述的多通道宽动态范围压缩系统,其特征在于,所述的步骤5)的P-1阶参数为b的全通反变换,可以通过P-1阶参数为-b的全通变换实现。
7.根据权利要求1所述的多通道宽动态范围压缩系统,其特征在于,对于输入信号x(n)(1)为实信号的情况,所述的声压级检测模块(3)和压缩放大增益计算模块(4)仅需对第0至K/2通道信号进行处理。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200910236221A CN101695148B (zh) | 2009-10-23 | 2009-10-23 | 一种用于数字助听器的多通道宽动态范围压缩系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN200910236221A CN101695148B (zh) | 2009-10-23 | 2009-10-23 | 一种用于数字助听器的多通道宽动态范围压缩系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101695148A CN101695148A (zh) | 2010-04-14 |
CN101695148B true CN101695148B (zh) | 2012-09-05 |
Family
ID=42094078
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200910236221A Expired - Fee Related CN101695148B (zh) | 2009-10-23 | 2009-10-23 | 一种用于数字助听器的多通道宽动态范围压缩系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101695148B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TWI594233B (zh) * | 2017-01-05 | 2017-08-01 | 元鼎音訊股份有限公司 | 具聽力保護功能之聲音調整裝置及其聲音調整之方法 |
TWI609367B (zh) * | 2016-10-20 | 2017-12-21 | 宏碁股份有限公司 | 電子裝置及利用窗化濾波器差異之特定頻段補償增益方法 |
TWI609365B (zh) * | 2016-10-20 | 2017-12-21 | 宏碁股份有限公司 | 助聽器及其寬動態範圍壓縮的恢復時間動態調整方法 |
TWI609366B (zh) * | 2016-10-21 | 2017-12-21 | 宏碁股份有限公司 | 電子裝置及其分頻濾波增益優化方法 |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102984634A (zh) * | 2011-11-22 | 2013-03-20 | 南京工程学院 | 一种数字助听器非等宽子带自动增益控制方法 |
CN104183242A (zh) * | 2014-08-07 | 2014-12-03 | 深圳市微纳集成电路与系统应用研究院 | 基于模块复用的语音处理方法的硬件设计方法及系统 |
CN104202098B (zh) * | 2014-08-08 | 2016-06-15 | 中国科学院上海微系统与信息技术研究所 | 一种基于多通道压缩采样的宽带功率谱估计方法 |
CN104902420B (zh) * | 2014-12-26 | 2019-04-09 | 深圳市微纳集成电路与系统应用研究院 | 基于曲线压缩实现的宽动态范围压缩方法和系统 |
CN104661168B (zh) * | 2015-01-13 | 2017-02-08 | 欧仕达听力科技(厦门)有限公司 | 一种基于子带宽动态范围压缩的最佳验配方法与系统 |
US10491179B2 (en) * | 2017-09-25 | 2019-11-26 | Nuvoton Technology Corporation | Asymmetric multi-channel audio dynamic range processing |
CN109951786A (zh) * | 2019-03-27 | 2019-06-28 | 钰太芯微电子科技(上海)有限公司 | 一种纯数字架构的助听器系统 |
CN113473287A (zh) * | 2021-06-23 | 2021-10-01 | 天津大学 | 一种新型蓝牙耳机输出声压级控制方法 |
CN113993053B (zh) * | 2021-10-27 | 2022-06-24 | 南京工程学院 | 一种通道自适应的数字助听器宽动态范围压缩方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1640190A (zh) * | 2001-08-08 | 2005-07-13 | Gn瑞声达公司 | 使用数字频率扭曲的动态范围压缩 |
CN1976226A (zh) * | 2006-12-20 | 2007-06-06 | 北京中星微电子有限公司 | 正交滤波器组设计方法及装置 |
-
2009
- 2009-10-23 CN CN200910236221A patent/CN101695148B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1640190A (zh) * | 2001-08-08 | 2005-07-13 | Gn瑞声达公司 | 使用数字频率扭曲的动态范围压缩 |
CN1976226A (zh) * | 2006-12-20 | 2007-06-06 | 北京中星微电子有限公司 | 正交滤波器组设计方法及装置 |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TWI609367B (zh) * | 2016-10-20 | 2017-12-21 | 宏碁股份有限公司 | 電子裝置及利用窗化濾波器差異之特定頻段補償增益方法 |
TWI609365B (zh) * | 2016-10-20 | 2017-12-21 | 宏碁股份有限公司 | 助聽器及其寬動態範圍壓縮的恢復時間動態調整方法 |
TWI609366B (zh) * | 2016-10-21 | 2017-12-21 | 宏碁股份有限公司 | 電子裝置及其分頻濾波增益優化方法 |
TWI594233B (zh) * | 2017-01-05 | 2017-08-01 | 元鼎音訊股份有限公司 | 具聽力保護功能之聲音調整裝置及其聲音調整之方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101695148A (zh) | 2010-04-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101695148B (zh) | 一种用于数字助听器的多通道宽动态范围压缩系统 | |
KR101294634B1 (ko) | 오디오 신호 처리 시스템 및 방법 | |
EP1433359B1 (en) | Dynamic range compression using digital frequency warping | |
JP5718251B2 (ja) | 分解されたオーディオ信号の再構成のためのシステムおよび方法 | |
CN101646123B (zh) | 一种模拟听觉感知模型的滤波器组 | |
JP5341128B2 (ja) | 補聴器における安定性の改善 | |
CN102638755A (zh) | 基于频率压缩搬移的数字助听器响度补偿方法 | |
US8948424B2 (en) | Hearing device and method for operating a hearing device with two-stage transformation | |
CN103731780A (zh) | 数字助听器信号处理方法 | |
EP2720477B1 (en) | Virtual bass synthesis using harmonic transposition | |
US11516581B2 (en) | Information processing device, mixing device using the same, and latency reduction method | |
Subbulakshmi et al. | A survey of filter bank algorithms for biomedical applications | |
CN103714825A (zh) | 基于听觉感知模型的多通道语音增强方法 | |
TWI421858B (zh) | 用於處理音頻訊號的系統及方法 | |
CN115299075B (zh) | 扬声器的低音增强 | |
Søndergaard et al. | On the relationship between multi-channel envelope and temporal fine structure | |
US20130216062A1 (en) | Hearing apparatus having an adaptive filter and method for filtering an audio signal | |
Parfieniuk et al. | Tunable non-uniform filter bank mixing cosine modulation with perceptual frequency warping by allpass transformation | |
Sokolova et al. | Multirate audiometric filter bank for hearing aid devices | |
RU2819779C1 (ru) | Усиление низких частот для громкоговорителей | |
Swamy et al. | Real-time Implementation of Delay Efficient DCT Based Hearing Aid Algorithm Using TMS320C5505 DSP Processor | |
Shrivastav et al. | An optimized frequency response masking reconfigurable filter to enhance the performance of the hearing aid system | |
BR112022018207B1 (pt) | Método de processamento de áudio implementado por computador, meio legível por computador não transitório e aparelho para processamento de áudio | |
Ying | Design of Computationally Efficient Digital FIR Filters and Filter Banks |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120905 Termination date: 20151023 |
|
EXPY | Termination of patent right or utility model |