CN1197243C - 一种利用频谱减法抑制噪声的数字滤波器设计方法和装置 - Google Patents
一种利用频谱减法抑制噪声的数字滤波器设计方法和装置 Download PDFInfo
- Publication number
- CN1197243C CN1197243C CN00812582.1A CN00812582A CN1197243C CN 1197243 C CN1197243 C CN 1197243C CN 00812582 A CN00812582 A CN 00812582A CN 1197243 C CN1197243 C CN 1197243C
- Authority
- CN
- China
- Prior art keywords
- estimated value
- density estimated
- high frequency
- power spectrum
- spectrum density
- 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
- 230000003595 spectral effect Effects 0.000 title claims abstract description 60
- 238000013461 design Methods 0.000 title claims abstract description 19
- 238000000034 method Methods 0.000 title claims description 24
- 230000001629 suppression Effects 0.000 title abstract description 3
- 238000001228 spectrum Methods 0.000 claims abstract description 84
- 238000012546 transfer Methods 0.000 claims abstract description 15
- 238000002592 echocardiography Methods 0.000 claims description 4
- 230000011218 segmentation Effects 0.000 abstract description 9
- 238000012935 Averaging Methods 0.000 abstract 2
- 230000008901 benefit Effects 0.000 description 3
- 101100382854 Arabidopsis thaliana CCD7 gene Proteins 0.000 description 2
- 101100129499 Arabidopsis thaliana MAX2 gene Proteins 0.000 description 2
- 102100022002 CD59 glycoprotein Human genes 0.000 description 2
- 101000897400 Homo sapiens CD59 glycoprotein Proteins 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000008030 elimination Effects 0.000 description 2
- 238000003379 elimination reaction Methods 0.000 description 2
- 230000008676 import Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 101100129496 Arabidopsis thaliana CYP711A1 gene Proteins 0.000 description 1
- 241000581364 Clinitrachus argentatus Species 0.000 description 1
- 101100083446 Danio rerio plekhh1 gene Proteins 0.000 description 1
- 241000406668 Loxodonta cyclotis Species 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Images
Classifications
-
- 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
- G10L21/00—Speech or voice signal processing techniques to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
- G10L21/02—Speech enhancement, e.g. noise reduction or echo cancellation
- G10L21/0208—Noise filtering
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H17/00—Networks using digital techniques
- H03H17/02—Frequency selective networks
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Audiology, Speech & Language Pathology (AREA)
- Human Computer Interaction (AREA)
- Computational Linguistics (AREA)
- Quality & Reliability (AREA)
- Signal Processing (AREA)
- Health & Medical Sciences (AREA)
- Computer Hardware Design (AREA)
- Mathematical Physics (AREA)
- Acoustics & Sound (AREA)
- Multimedia (AREA)
- Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
- Noise Elimination (AREA)
- Picture Signal Circuits (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Image Processing (AREA)
Abstract
一种用于利用频谱减法抑制噪声的数字滤波器设计装置,包括一个第一频谱估算器(12),用于从带有噪声的语音信号块确定高频率分辨率的带有噪声的语音功率谱密度估算值。一个第二频谱估算器(24)从背景噪声信号块确定高频率分辨率的背景噪声功率谱密度估算值。求平均单元(20,26)形成分段式恒定的带有噪声的语音功率谱密度估算值和分段式恒定的背景噪声功率谱密度估算值。这些求平均的单元由装置(14,16,18)控制,使单个段的长度与高频率分辨率的有噪声的语音功率谱密度估算值的形状相适应,并对两个分段式恒定的估算值使用相同的段。装置(28)利用基于分段式恒定的带有噪声的语音功率谱密度估算值和分段式恒定的背景噪声功率谱密度估算值的频谱减法来确定分段式恒定数字滤波器传递函数。
Description
技术领域
本发明涉及一种数字滤波器设计,尤其是频域内滤波器的设计。
背景技术
在几例应用中,频域内数字滤波器H[k]是在“匆忙中”设计的。一个例子是利用频谱减法的噪声抑制。另一个例子是用于回波消除的频率选择非线性处理器的设计。这些应用的特有特征是滤波器设计方法非常复杂。由于滤波器频繁更新,给执行这些设计算法的硬件/软件带来了沉重的负担。
参考文献[1]描述了一种方法,该方法把频域k分成等长或不等长的几段,并且对每一段内的H[K]和噪声或回波污染的语音信号的基本功率谱密度估算值Φx[k]作用一个常量值。这就减少了复杂性,因为只需对频段而不是对每一个频率元(frequency bin)k确定滤波器H[k]。但是,这一方法也有缺点,即可能把H[k]的一个峰拆分为两个不同的段。这可能会引起波动峰,而产生恼人的“音乐噪声”。也会降低频谱的锐度,进一步降低语音质量。
发明内容
本发明的一个目的是减少或消除现有技术的这些缺点。
根据所附的权利要求,可实现这一目的。
简要地,本发明使段的长度和位置动态地与语音信号功率谱的当前形状相适配。该谱的峰和谷被确定,该方法保证在频域内分配各段时峰不被拆分在不同的段当中。优选地,每一个峰由其中心对在峰上的段所覆盖。段长优选地由人类听觉系统的频率特性来控制。
这一方法的优点是,既降低滤波器计算的复杂性,又不牺牲在重要谱峰的精确性。而且,该方法也减小了帧与帧之间谱的变化,从而提高语音质量。
根据本发明的一方面,提供了一种用于利用频谱减法抑制噪声的数字滤波器设计方法,包括步骤:根据带有噪声的语音信号块来确定高频率分辨率的带有噪声的语音功率谱密度估算值;通过在所述带有噪声的语音高频率分辨率的功率谱密度估算值的频率元段内对功率密度求平均来形成分段式常数的带有噪声的语音功率谱密度估算值;根据背景噪声信号块来确定高频率分辨率的背景噪声功率谱密度估算值;通过在所述高频率分辨率的背景噪声功率谱密度估算值的频率元段内对功率密度求平均来形成分段式常数的背景噪声功率谱密度估算值;利用基于所述分段式常数的带有噪声的语音功率谱密度估算值和所述分段式常数的背景噪声功率谱密度估算值的频谱减法来确定分段式常数的数字滤波器传递函数,其特征在于:将各个的段的长度以及位置与所述高频率分辨率的带有噪声的语音功率谱密度估算值的形状动态地相适配,使在频率域上分配各段时峰不被拆分在不同的段当中;并且对所述高频率分辨率的带有噪声的语音功率谱密度估算值和所述高频率分辨率的背景噪声功率谱密度估算值使用相同的段。
根据本发明的另一方面,提供了一种用于非线性回波消除的数字滤波器设计方法,其步骤包括:根据包含语音的残余回波信号块来确定高频率分辨率的语音功率谱密度估算值;通过在所述高频率分辨率的语音功率谱密度估算值的频率元段内对功率密度求平均来形成分段式常数的语音功率谱密度估算值;根据回波信号块来确定高频率分辨率的残余回波功率谱密度估算值;通过在所述高频率分辨率的残余回波功率谱密度估算值的频率元段内对功率密度求平均来形成分段式常数的残余回波功率谱密度估算值;利用所述分段式常数的语音功率谱密度估算值和所述分段式常数的残余回波功率谱密度估算值来确定分段式常数数字滤波器传递函数,其特征在于:将各个段的长度以及位置与所述高频率分辨率的语音功率谱密度估算值的形状动态地相适配,使在频率域上分配各段时峰不被拆分在不同的段当中;并且对所述高频率分辨率的语音功率谱密度估算值和所述高频率分辨率的残余回波功率谱密度估算值使用相同的段。
根据本发明的另一方面,提供了一种利用频谱减法抑制噪声的数字滤波器装置,包括:用于根据带有噪声的语音信号块来确定高频率分辨率的带有噪声的语音功率谱密度估算值的装置;用于通过在所述带有噪声的语音的高频率分辨率的功率谱密度估算值的频率元段内对功率密度求平均来形成分段式常数的带有噪声的语音功率谱密度估算值的装置;用于根据背景噪声信号块来确定高频率分辨率的背景噪声功率谱密度估算值的装置;用于通过在所述高频率分辨率的背景噪声功率谱密度估算值的频率元段内对功率密度求平均来形成分段式常数的背景噪声功率谱密度估算值的装置;用于利用基于所述分段式常数的带有噪声的语音功率谱密度估算值和所述分段式常数的背景噪声功率谱密度估算值的频谱减法来确定分段式常数的数字滤波器传递函数的装置,其特征在于:用于将各个段的长度以及位置与所述高频率分辨率的带有噪声的语音功率谱密度估算值的形状动态地相适配使在频率域内分配各段时峰不被拆分在不同的段当中的装置(14,18);和用于对所述高频率分辨率的带有噪声的语音功率谱密度估算值和所述高频率分辨率的背景噪声功率谱密度估算值使用相同的段的装置(18)。
根据本发明的另一方面,提供了一种用于非线性回波消除的数字滤波器装置,包括:用于根据包含语音的残余回波信号块来确定高频率分辨率语音功率谱密度估算值的装置;用于通过在所述高频率分辨率语音功率谱密度估算值的频率元段内对功率密度求平均来形成分段式常数语音功率谱密度估算值的装置;用于根据回波信号块来确定高频率分辨率的残余回波功率谱密度估算值的装置;用于通过在所述高频率分辨率的残余回波功率谱密度估算值的频率元段内对功率密度求平均来形成分段式常数的残余回波功率谱密度估算值的装置;用于利用所述分段式常数语音功率谱密度估算值和所述分段式常数的残余回波功率谱密度估算值来确定分段式常数的数字滤波器传递函数的装置,其特征在于:用于将各个段的长度以及位置与所述高频率分辨率语音功率谱密度估算值的形状动态地相适配使在频率域内分配各段时峰不被拆分在不同的段当中的装置(14,18);和用于对所述高频率分辨率的语音功率谱密度估算值和所述高频率分辨率的残余回波功率谱密度估算值使用相同的段的装置(18)。
附图说明
参考下列附图和其说明,本发明及其其它的目的和优点可得到很好的理解,其中:
图1示出有噪声的语音信号的功率谱密度估算值图;
图2示出根据现有技术分段和求平均后,图1中功率谱密度估算值图;
图3是与图1相似的图;
图4示出根据本发明分段和求平均后,图1中功率谱密度估算值图;
图5是根据本发明的方法的典型实施方案的流程图;和
图6是根据本发明的滤波器设计和滤波设备的典型实施方案的方框图。
具体实施方式
在一些应用中,滤波器在频域中确定。例如,在电话应用中,经常使用基于频谱相减的噪声抑制(见[2,3])。这种情况下滤波器作为频率的函数H(ω)被确定:
x[n]=y[n]+v[n]其中v[n]是噪声信号,x[n]是有噪声的语音信号,y[n]是欲得信号。通过把H(ω)代表的滤波器应用于有噪声的信号x[n]就得到了欲得信号y[n]的估算值。
另一个滤波器在频域中被确定的应用实例是用于回波消除的频率选择非线性处理器。在这种情况中,滤波器由下述函数定义:
其中
代表被残余回波污染的信号x[n]的功率谱密度的估算值,
代表残余回波信号e[n]的功率谱密度的估算值。一个恰当函数f的例子是上述的用于噪声抑制函数,
在此代表被残余回波污染的信号的功率谱密度估算值,
则被残余回波功率谱密度的估算值
所代替。滤波器H(ω)基于以下模型:
x[n]=y[n]+e[n]其中y[n]是欲得信号。通过把H(ω)代表的滤波器应用于被残余回波污染的信号x[n],就得到欲得信号y[n]的估算值。
在上述的实例中,滤波器由一个实值连续频率传递函数H(ω)描述。这一函数被采样得到离散频率传递函数H[k]。典型地当估算值基于参数估算法时,执行该步骤。但是,也可能直接得到离散频率传递函数H[k],例如通过使用基于估算法的周期图。参数估算法的一个优点是,其估算值在帧与帧之间的变化一般低于基于方法的周期图的估算值。
下列描述将限于噪声抑制,但是可以理解的是,同样的原理也可用于其他应用,如回波消除。
离散频率功率谱估算值
首先公知地具有高频率分辨率,典型地为128或256频率元。图1示出有噪声的语音信号的功率谱密度估算值
在这里具有256频率元。该频谱由参数估算法(典型地为一种自动回归模型,例如指令10)得到。如果滤波器传递函数H[k]欲以同样的分辨率来确定,这将给执行上述噪声抑制算法的硬件/软件带来巨大的计算量。根据[1],频率范围因此被分成常数长度的段,并且在每段中形成
的平均值,如图2所示。在计算H[k]时,该平均值
代替
相似的分段的和求平均的估算值
代替
这样,由
H[Segment]代表,并由公式
定义的H[k]的单个值就可以用于整个段而不仅仅是一个k值。这一方法的缺点是由于段的常数长度和位置,使得频谱的重要峰可能被拆分在几段内,如图2所示的MAX2和MAX3。这会引起这些峰的低分辨率。而且,由于从语音帧到语音帧,这些峰的位置可能移动,则它们有时被分开而有时不被分开,这将导致非常恼人的“音乐噪声”。
如下列典范的算法所示,本发明通过使段的长度和位置动态与当前估算值
的形状相适配来解决这一问题。简要地,算法由找到的局部极大(峰位)和局部极小(谷位)而开始。然后使段对中在每个极大上,并在峰之间分配段以覆盖谷。当定位极大和极小时,基于自动回归的参数频谱估算尤其吸引人,因为这样的频谱确保具有至多M/2个峰,其中M是模型指令。
优选地,使各个段的长度适配于人类听觉系统的特性,这在[4]中已经研究过。基于这一信息,可得到下列段中心fc和段长度的关系(其中频率范围为256元,采样频率为8000Hz,这样获得的频率分辨率为31.25Hz/bin)。
转化为离散频率域k,则给出:
利用这一关系,下列在图5中也示出的算法可利用动态确定的段长和位置来确一个分段的和求平均的滤波器传递函数H[k]:
S1:得到下一个x[n]的信号块
S4:设第一极大的k值为kmax
S5:设第一极小的k值为kmin
S6:设kc=kmax
用求平均的
和
确定
H[segment]
S8:如果kc-segment(kc)/2>kmin,则执行S9-S10
S9:设kc=kc-segment(kc)(在segment(kc)中用旧值kc)
S10:在以kc为中心的Segment[kc]内确定
的平均值。
确定同一段中
的平均值
用求平均的
和
确定
H[segment]
跳转到S8
S11:设kc=kmax
S12:设下一个极小的k值为kmin
S13:如果kc+segment(kc)/2<kmin,则执行S14-S15
S14:设kc=kc-segment(kc)(在segment(kc)中用旧值kc)
S15:在以kc为中心的Segment(kc)内确定
的平均值。
确定同一段中
的平均值
跳转到S13
S16:如果kmax是最后一个极大值,跳转到S1
S17:设下一个极大的k值为kmax,并跳转到S6
把这一算法用于图1的频谱,产生图4中分段的和求平均的频谱。如其所示,图4中局部极大位于:
MAX1:k=20
MAX2:k=41
MAX3:k=73
局部极小位于:
MIN1:k=0
MIN2:k=31
MIN3:k=61
MIN4:k=128
例如,把以上算法应用于第二个极大值k=41处,给出以k=41为中心的一个5元段,再给出该极大值左侧的两个5元段,以及该极大值的右侧的两个7元段。当把图2与图1比较时,会发现覆盖峰的段总是以峰为中心。还会发现,越低的频率导致越短的段。如上所见,如果从帧到帧,峰改变了位置,该算法将保证段仍以峰为中心,并且段宽适合于峰的位置。
图6示出的是根据本发明的滤波器设计装置的典范的实施方案方框图,这里用于利用频率谱减法的噪声抑制。有噪声的语音信号采样值x[n]流被传入收集信号瞬时值块或帧的缓存器10。频谱估算器12找到该块的自动回归(AR)参数并利用这些参数来确定有噪声的语音信号x[n]当前块的功率谱密度估算值
典型地,这一估算值具有128或256个信号采样值。极大-极小检测器14在估算值中搜索局部极大值和极小值。局部极大值和极小值的位置被传入根据参考图5所述方法分配各段的段分配器18。段的位置和长度被转入平均器20。平均器20接收估算值
的采样值,并在每一个特定的段形成平均值。在一个没有语音的块期间,背景噪声信号v[n]的一个块被收集于缓存器22中。频谱估算器24找到该块的自动回归参数并利用这些参数来确定该背景噪声信号v[n]块的功率谱密度估算值
这一估算值具有与估算值
同样数量的采样值。来自段分配器18的段的位置和长度也被传入另一个平均器26。平均器26接收估算值
的采样值,并在每个特定的段形成平均值。平均器20、26都把每一段的平均值传入确定每一段滤波器传递函数值的滤波器计算器28。这将产生由块30代表的分段的滤波器
H[k]。该滤波器被传入乘法器32的一个输入端。缓存器10中的信号块也被传入到把块变换到频域的快速傅立叶变换(FFT)块34。变换后的块的长度与分段的滤波器的长度相同。变换后的信号传入乘法器32的另一个输入端在乘法器32中它与分段后的滤波器相乘。最后,滤波后的信号在快速傅立叶反变换(IFFT)块36中被变换回到时间域。
典型地,图6中不同的块由一个或几个微缩处理器或,微/信号处理器的组合来实现。但是,它们也可由一个或几个ASIC(应用专门集成电路)来实现。
相似的结构可用于回波消除中的非线性滤波。这种情况下,x[n]代表残余回波污染的信号,v[n]则由残余回波e[n]的估算值所代替。另一个差别是,这种情况下,估算值
和
出自同一个语音帧(在利用频谱减法的噪声抑制中,噪声频谱被认为是静态的且在语音停顿期间进行估算)。
那些本领域的技术人员应当理解的是,不偏离所附权利要求其所限定的范围,可对本发明作出各种各样的修改和改变。
参考文献
[1]美国专利第5,839,101号(A.Vhitalo et al)。
[2]J.S.Lim和V.V.Oppenheim,“含噪声的语音增强和带宽压缩”,IEEE学报,卷67,第12,1979期,1586-1604页。
[3]S.F.Boll,“利用频谱减法抑制语音中的噪声”,IEEE期刊,并于声学,语音和信号处理,卷ASSP27,第21979号,113-120页。
[4]U.Zlser,“数字音频信号处理”,John Wiley & Sons,Chichester,英国,1997,252-253页。
Claims (12)
1.一种用于利用频谱减法抑制噪声的数字滤波器设计方法,包括步骤:
根据带有噪声的语音信号块来确定高频率分辨率的带有噪声的语音功率谱密度估算值,
通过在所述带有噪声的语音高频率分辨率的功率谱密度估算值的频率元段内对功率密度求平均来形成分段式常数的带有噪声的语音功率谱密度估算值,
根据背景噪声信号块来确定高频率分辨率的背景噪声功率谱密度估算值,
通过在所述高频率分辨率的背景噪声功率谱密度估算值的频率元段内对功率密度求平均来形成分段式常数的背景噪声功率谱密度估算值,
利用基于所述分段式常数的带有噪声的语音功率谱密度估算值和所述分段式常数的背景噪声功率谱密度估算值的频谱减法来确定分段式常数的数字滤波器传递函数,其特征在于
将各个的段的长度以及位置与所述高频率分辨率的带有噪声的语音功率谱密度估算值的形状动态地相适配,使在频率域上分配各段时峰不被拆分在不同的段当中;并且
对所述高频率分辨率的带有噪声的语音功率谱密度估算值和所述高频率分辨率的背景噪声功率谱密度估算值使用相同的段。
2.权利要求1的数字滤波器设计方法,其特征在于使段的中心位于所述高频率分辨率的带有噪声的语音功率谱密度估算值的局部极大处。
3.权利要求1或2的数字滤波器设计方法,其特征在于对于高频率根据人类听觉系统增加所述段的长度。
4.一种用于非线性回波消除的数字滤波器设计方法,其步骤包括:
根据包含语音的残余回波信号块来确定高频率分辨率的语音功率谱密度估算值,
通过在所述高频率分辨率的语音功率谱密度估算值的频率元段内对功率密度求平均来形成分段式常数的语音功率谱密度估算值,
根据回波信号块来确定高频率分辨率的残余回波功率谱密度估算值,
通过在所述高频率分辨率的残余回波功率谱密度估算值的频率元段内对功率密度求平均来形成分段式常数的残余回波功率谱密度估算值,
利用所述分段式常数的语音功率谱密度估算值和所述分段式常数的残余回波功率谱密度估算值来确定分段式常数数字滤波器传递函数,其特征在于
将各个段的长度以及位置与所述高频率分辨率的语音功率谱密度估算值的形状动态地相适配,使在频率域上分配各段时峰不被拆分在不同的段当中;并且
对所述高频率分辨率的语音功率谱密度估算值和所述高频率分辨率的残余回波功率谱密度估算值使用相同的段。
5.权利要求4的数字滤波器设计方法,其特征在于使段的中心位于所述高频率分辨率语音功率谱密度估算值的局部极大处。
6.权利要求4或5的数字滤波器设计方法,其特征在于对于高频率根据人类听觉系统增加所述段的长度。
7.一种利用频谱减法抑制噪声的数字滤波器装置,包括:
用于根据带有噪声的语音信号块来确定高频率分辨率的带有噪声的语音功率谱密度估算值的装置,
用于通过在所述带有噪声的语音的高频率分辨率的功率谱密度估算值的频率元段内对功率密度求平均来形成分段式常数的带有噪声的语音功率谱密度估算值的装置,
用于根据背景噪声信号块来确定高频率分辨率的背景噪声功率谱密度估算值的装置,
用于通过在所述高频率分辨率的背景噪声功率谱密度估算值的频率元段内对功率密度求平均来形成分段式常数的背景噪声功率谱密度估算值的装置,
用于利用基于所述分段式常数的带有噪声的语音功率谱密度估算值和所述分段式常数的背景噪声功率谱密度估算值的频谱减法来确定分段式常数的数字滤波器传递函数的装置,其特征在于
用于将各个段的长度以及位置与所述高频率分辨率的带有噪声的语音功率谱密度估算值的形状动态地相适配使在频率域内分配各段时峰不被拆分在不同的段当中的装置(14,18);和
用于对所述高频率分辨率的带有噪声的语音功率谱密度估算值和所述高频率分辨率的背景噪声功率谱密度估算值使用相同的段的装置(18)。
8.权利要求7的数字滤波器装置,其特征在于使段的中心位于所述高频率分辨率的带有噪声的语音功率谱密度估算值的局部极大处的装置(18)。
9.权利要求7或8的数字滤波器装置,其特征在于用于对于高频率根据人类听觉系统增加所述段的长度的装置(18)。
10.一种用于非线性回波消除的数字滤波器装置,包括:
用于根据包含语音的残余回波信号块来确定高频率分辨率语音功率谱密度估算值的装置,
用于通过在所述高频率分辨率语音功率谱密度估算值的频率元段内对功率密度求平均来形成分段式常数语音功率谱密度估算值的装置,
用于根据回波信号块来确定高频率分辨率的残余回波功率谱密度估算值的装置,
用于通过在所述高频率分辨率的残余回波功率谱密度估算值的频率元段内对功率密度求平均来形成分段式常数的残余回波功率谱密度估算值的装置,
用于利用所述分段式常数语音功率谱密度估算值和所述分段式常数的残余回波功率谱密度估算值来确定分段式常数的数字滤波器传递函数的装置,其特征在于
用于将各个段的长度以及位置与所述高频率分辨率语音功率谱密度估算值的形状动态地相适配使在频率域内分配各段时峰不被拆分在不同的段当中的装置(14,18);和
用于对所述高频率分辨率的语音功率谱密度估算值和所述高频率分辨率的残余回波功率谱密度估算值使用相同的段的装置(18)。
11.权利要求10的数字滤波器装置,其特征在于使段的中心位于所述高频率分辨率语音功率谱密度估算值的局部极大处的装置(18)。
12.权利要求10或11的数字滤波器装置,其特征在于用来对于高频率根据人类听觉系统增加所述段的长度的装置(20)。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
SE99031601 | 1999-09-07 | ||
SE9903160A SE514875C2 (sv) | 1999-09-07 | 1999-09-07 | Förfarande och anordning för konstruktion av digitala filter |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1373930A CN1373930A (zh) | 2002-10-09 |
CN1197243C true CN1197243C (zh) | 2005-04-13 |
Family
ID=20416885
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN00812582.1A Expired - Fee Related CN1197243C (zh) | 1999-09-07 | 2000-08-23 | 一种利用频谱减法抑制噪声的数字滤波器设计方法和装置 |
Country Status (8)
Country | Link |
---|---|
US (1) | US6564184B1 (zh) |
EP (1) | EP1221197B1 (zh) |
CN (1) | CN1197243C (zh) |
AT (1) | ATE445260T1 (zh) |
AU (1) | AU6884800A (zh) |
DE (1) | DE60043115D1 (zh) |
SE (1) | SE514875C2 (zh) |
WO (1) | WO2001018961A1 (zh) |
Families Citing this family (28)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6859531B1 (en) * | 2000-09-15 | 2005-02-22 | Intel Corporation | Residual echo estimation for echo cancellation |
EP1278185A3 (en) * | 2001-07-13 | 2005-02-09 | Alcatel | Method for improving noise reduction in speech transmission |
US7187719B2 (en) * | 2003-06-18 | 2007-03-06 | Mindspeed Technologies, Inc. | Method and system for data rate optimization in a digital communication system |
CN100336307C (zh) * | 2005-04-28 | 2007-09-05 | 北京航空航天大学 | 接收机射频系统电路内部噪声的分配方法 |
CN100370519C (zh) * | 2005-10-17 | 2008-02-20 | 西安交通大学 | 一种增强电子喉语音的方法及其系统 |
US20070237341A1 (en) * | 2006-04-05 | 2007-10-11 | Creative Technology Ltd | Frequency domain noise attenuation utilizing two transducers |
US8761387B2 (en) | 2006-05-04 | 2014-06-24 | Mindspeed Technologies, Inc. | Analog transmit crosstalk canceller |
US8300849B2 (en) * | 2007-11-06 | 2012-10-30 | Microsoft Corporation | Perceptually weighted digital audio level compression |
WO2009082299A1 (en) * | 2007-12-20 | 2009-07-02 | Telefonaktiebolaget L M Ericsson (Publ) | Noise suppression method and apparatus |
KR20090110244A (ko) * | 2008-04-17 | 2009-10-21 | 삼성전자주식회사 | 오디오 시맨틱 정보를 이용한 오디오 신호의 부호화/복호화 방법 및 그 장치 |
CN101615538B (zh) * | 2008-06-27 | 2013-01-30 | 上海亿盟电气自动化技术有限公司 | 一种脱扣器滤波方法 |
CN101320566B (zh) * | 2008-06-30 | 2010-10-20 | 中国人民解放军第四军医大学 | 基于多带谱减法的非空气传导语音增强方法 |
CN101778183B (zh) * | 2009-01-13 | 2013-06-05 | 华为终端有限公司 | 一种残留回声抑制方法及设备 |
US8229126B2 (en) * | 2009-03-13 | 2012-07-24 | Harris Corporation | Noise error amplitude reduction |
DE112010004988B4 (de) * | 2009-12-25 | 2023-03-30 | Mitsubishi Electric Corporation | Störungsbeseitungsvorrichtung und Störungsbeseitigungsprogramm |
US8600312B2 (en) * | 2010-01-25 | 2013-12-03 | Qualcomm Incorporated | Method and apparatus for spectral sensing |
CN102158197A (zh) * | 2010-02-11 | 2011-08-17 | 北京国能子金电气技术有限公司 | 变频率信号的自适应低通数字滤波器和方法 |
US9648421B2 (en) | 2011-12-14 | 2017-05-09 | Harris Corporation | Systems and methods for matching gain levels of transducers |
CN103021420B (zh) * | 2012-12-04 | 2015-02-25 | 中国科学院自动化研究所 | 一种基于相位调整和幅值补偿的多子带谱减法的语音增强方法 |
CN103632677B (zh) | 2013-11-27 | 2016-09-28 | 腾讯科技(成都)有限公司 | 带噪语音信号处理方法、装置及服务器 |
CN104143337B (zh) * | 2014-01-08 | 2015-12-09 | 腾讯科技(深圳)有限公司 | 一种提高音频信号音质的方法和装置 |
EP3125429A1 (en) * | 2015-07-28 | 2017-02-01 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. | Partitioned block frequency domain adaptive filter device comprising adaptation modules and correction modules |
CN105812995B (zh) * | 2016-03-09 | 2018-11-09 | 西安航空学院 | 一种噪音滤除方法和装置 |
CN105877706A (zh) * | 2016-03-31 | 2016-08-24 | 济南大学 | 一种基于改进谱减法的心音增强方法 |
CN105827218B (zh) * | 2016-05-17 | 2018-04-06 | 中国电子科技集团公司第四十一研究所 | 一种对fft数据实施的噪声带限滤波方法 |
US10481831B2 (en) * | 2017-10-02 | 2019-11-19 | Nuance Communications, Inc. | System and method for combined non-linear and late echo suppression |
EP3667420A1 (en) * | 2018-12-13 | 2020-06-17 | ASML Netherlands B.V. | Method of measuring a parameter of a lithographic process, metrology apparatus |
CN111239839B (zh) * | 2020-02-10 | 2021-04-09 | 山东大学 | 一种针对频率域电磁法的频谱密度计算方法及装置 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4658426A (en) * | 1985-10-10 | 1987-04-14 | Harold Antin | Adaptive noise suppressor |
IT1184023B (it) * | 1985-12-17 | 1987-10-22 | Cselt Centro Studi Lab Telecom | Procedimento e dispositivo per la codifica e decodifica del segnale vocale mediante analisi a sottobande e quantizzazione vettorariale con allocazione dinamica dei bit di codifica |
CN1062963C (zh) * | 1990-04-12 | 2001-03-07 | 多尔拜实验特许公司 | 用于产生高质量声音信号的解码器和编码器 |
DK173156B1 (da) | 1994-10-05 | 2000-02-14 | Damixa As | Etgrebsblandingsarmatur |
SE505156C2 (sv) * | 1995-01-30 | 1997-07-07 | Ericsson Telefon Ab L M | Förfarande för bullerundertryckning genom spektral subtraktion |
FI100840B (fi) | 1995-12-12 | 1998-02-27 | Nokia Mobile Phones Ltd | Kohinanvaimennin ja menetelmä taustakohinan vaimentamiseksi kohinaises ta puheesta sekä matkaviestin |
US5729694A (en) * | 1996-02-06 | 1998-03-17 | The Regents Of The University Of California | Speech coding, reconstruction and recognition using acoustics and electromagnetic waves |
FR2748342B1 (fr) * | 1996-05-06 | 1998-07-17 | France Telecom | Procede et dispositif de filtrage par egalisation d'un signal de parole, mettant en oeuvre un modele statistique de ce signal |
-
1999
- 1999-09-07 SE SE9903160A patent/SE514875C2/sv not_active IP Right Cessation
-
2000
- 2000-08-23 CN CN00812582.1A patent/CN1197243C/zh not_active Expired - Fee Related
- 2000-08-23 WO PCT/SE2000/001609 patent/WO2001018961A1/en active Application Filing
- 2000-08-23 DE DE60043115T patent/DE60043115D1/de not_active Expired - Lifetime
- 2000-08-23 EP EP00957195A patent/EP1221197B1/en not_active Expired - Lifetime
- 2000-08-23 AT AT00957195T patent/ATE445260T1/de not_active IP Right Cessation
- 2000-08-23 AU AU68848/00A patent/AU6884800A/en not_active Abandoned
- 2000-09-06 US US09/656,295 patent/US6564184B1/en not_active Expired - Lifetime
Also Published As
Publication number | Publication date |
---|---|
US6564184B1 (en) | 2003-05-13 |
ATE445260T1 (de) | 2009-10-15 |
WO2001018961A1 (en) | 2001-03-15 |
DE60043115D1 (de) | 2009-11-19 |
EP1221197A1 (en) | 2002-07-10 |
AU6884800A (en) | 2001-04-10 |
CN1373930A (zh) | 2002-10-09 |
SE9903160D0 (sv) | 1999-09-07 |
SE9903160L (sv) | 2001-03-08 |
SE514875C2 (sv) | 2001-05-07 |
EP1221197B1 (en) | 2009-10-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1197243C (zh) | 一种利用频谱减法抑制噪声的数字滤波器设计方法和装置 | |
Wang et al. | The unimportance of phase in speech enhancement | |
CN1175709C (zh) | 利用频谱减除的双拾音器信号降噪系统和方法 | |
CN1083183C (zh) | 用来降低语音信号中噪声的方法和装置 | |
JP3484801B2 (ja) | 音声信号の雑音低減方法及び装置 | |
CN1145931C (zh) | 降低语音信号中的噪声的方法和采用该方法的系统和电话机 | |
CN101080766A (zh) | 使用bark频带weiner滤波器和线性衰减的噪声降低和舒适噪声增益控制 | |
CN1277500A (zh) | 噪声抑制装置 | |
CN1286788A (zh) | 关于低比特率语音编码器的噪声抑制 | |
CA2436318A1 (fr) | Procede et dispositif de reduction de bruit | |
CN1746974A (zh) | 增强语音质量的方法及其装置 | |
CN101038744A (zh) | 音频信号噪声降低设备和方法 | |
US20040148166A1 (en) | Noise-stripping device | |
EP0998166A1 (fr) | Dispositif de traitement audio récepteur et procédé pour filtrer un signal utile et le restituer en présence de bruit ambiant | |
CN1649374A (zh) | 具有任意重放采样率的回波消除系统和方法 | |
CA2413867A1 (en) | Noise reduction method and apparatus | |
Butucea et al. | Sharp optimality in density deconvolution with dominating bias. I | |
JP2001265367A (ja) | 音声区間判定装置 | |
Sun et al. | Average sampling in shift invariant subspaces with symmetric averaging functions | |
JP2004341339A (ja) | 雑音抑圧装置 | |
Nuha et al. | Noise reduction and speech enhancement using wiener filter | |
EP0534837B1 (fr) | Procédé de traitement de la parole en présence de bruits acoustiques utilisant la sous traction spectrale non-linéaire et les modèles de Markov cachés | |
US6804646B1 (en) | Method and apparatus for processing a sound signal | |
JP3279254B2 (ja) | スペクトル雑音除去装置 | |
US7177805B1 (en) | Simplified noise suppression circuit |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
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: 20050413 Termination date: 20150823 |
|
EXPY | Termination of patent right or utility model |