CN1296258A - 噪声抑制装置 - Google Patents
噪声抑制装置 Download PDFInfo
- Publication number
- CN1296258A CN1296258A CN00122136.1A CN00122136A CN1296258A CN 1296258 A CN1296258 A CN 1296258A CN 00122136 A CN00122136 A CN 00122136A CN 1296258 A CN1296258 A CN 1296258A
- Authority
- CN
- China
- Prior art keywords
- frequency spectrum
- noise
- amplitude
- spectrum
- amplitude frequency
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000001228 spectrum Methods 0.000 claims abstract description 560
- 230000001629 suppression Effects 0.000 claims abstract description 48
- 238000012937 correction Methods 0.000 claims abstract description 27
- 238000009499 grossing Methods 0.000 claims description 109
- 230000003595 spectral effect Effects 0.000 claims description 93
- 230000008859 change Effects 0.000 claims description 13
- 230000000452 restraining effect Effects 0.000 claims description 7
- 230000009466 transformation Effects 0.000 claims description 4
- 230000006866 deterioration Effects 0.000 abstract 1
- 230000002349 favourable effect Effects 0.000 abstract 1
- 230000000694 effects Effects 0.000 description 15
- 238000012545 processing Methods 0.000 description 15
- 238000000034 method Methods 0.000 description 13
- 238000001914 filtration Methods 0.000 description 8
- 230000008569 process Effects 0.000 description 8
- 230000005236 sound signal Effects 0.000 description 8
- 230000035807 sensation Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 239000000203 mixture Substances 0.000 description 6
- 239000000284 extract Substances 0.000 description 5
- 230000005764 inhibitory process Effects 0.000 description 5
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 4
- 230000002159 abnormal effect Effects 0.000 description 4
- 238000005070 sampling Methods 0.000 description 4
- 238000000605 extraction Methods 0.000 description 3
- 230000008676 import Effects 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 230000008034 disappearance Effects 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 230000011218 segmentation Effects 0.000 description 2
- 238000005728 strengthening Methods 0.000 description 2
- 101000822695 Clostridium perfringens (strain 13 / Type A) Small, acid-soluble spore protein C1 Proteins 0.000 description 1
- 101000655262 Clostridium perfringens (strain 13 / Type A) Small, acid-soluble spore protein C2 Proteins 0.000 description 1
- 101000655256 Paraclostridium bifermentans Small, acid-soluble spore protein alpha Proteins 0.000 description 1
- 101000655264 Paraclostridium bifermentans Small, acid-soluble spore protein beta Proteins 0.000 description 1
- 206010038743 Restlessness Diseases 0.000 description 1
- 230000005534 acoustic noise Effects 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000012790 confirmation Methods 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000003467 diminishing effect Effects 0.000 description 1
- 230000009467 reduction 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
Landscapes
- Engineering & Computer Science (AREA)
- Computational Linguistics (AREA)
- Quality & Reliability (AREA)
- Signal Processing (AREA)
- Health & Medical Sciences (AREA)
- Audiology, Speech & Language Pathology (AREA)
- Human Computer Interaction (AREA)
- Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Multimedia (AREA)
- Noise Elimination (AREA)
- Compression, Expansion, Code Conversion, And Decoders (AREA)
Abstract
提供听觉上良好可抑制噪声即使在高噪声情况下品质劣化也小的噪声抑制装置。该噪声抑制装置频率分析输入信号,用噪声振幅频谱计算装置算出用于该帧的输入振幅频谱的噪声振幅频谱,用频谱校正增益计算装置算出噪声振幅频谱校正增益和噪声除去频谱校正增益,用频谱减法装置输出第一噪声除去频谱,频谱振幅抑制装置将第二噪声除去频谱变换成时间轴信号。
Description
本发明涉及在各种噪声环境下使用的声音通信***和声音确认***等中抑制目标信号以外噪声的噪声抑制装置。
作为抑制在声音信号上重叠的噪声等目标外信号的噪声抑制装置,例如披露于特开平8-221093号中。它着眼于文献1(Steven F.Boll,“Suppression of Acoustic noise in speech using spectralsubtraction”,IEE ETrans.ASSP,Vol.ASSP-27,No.2,April 1997)中所示的振幅频谱,是基于所谓的频谱减法(SpectralSubtraction:SS)方法的装置。
使用图13来说明在特开平8-221093号中披露的以往的噪声抑制装置的结构。在图13中,101是帧化处理部分,102是开窗口处理部分,103是快速傅立叶变换部分,104是频带分割部分,105是噪声估计部分,106是NR值计算部分,107是Hn值计算部分,108是滤波器处理部分,109是频带变换部分,110是频谱修正部分,111是逆快速傅立叶变换部分,112是重叠加法部分,113是声音信号输入端子,114是声音信号输出端子,114是运算装置,115是输出信号端。此外,作为噪声估计部分105的内部结构,121是RMS计算部分,122是相对能量计算部分,123是最大RMS计算部分,124是估计噪声电平计算部分,125是最大SNR计算部分,126是噪声频谱估计部分。
以下说明以往的噪声抑制装置的工作原理。
在声音信号输入端子113上,输入包括声音和噪声成分的输入声音信号y[t]。该输入信号y[t]例如是抽样频率为FS的数字信号,向帧化处理部分101传送,帧长度被分割成每个FL抽样的帧,进行对以下各帧的处理。
在开窗口处理部分102中,在快速傅立叶变换处理部分103的计算之前,对于由帧化处理部分101传送的各帧化信号yframe[j,k]进行开窗口处理。但是,j是抽样号码,k是帧号码。
在快速傅立叶变换部分103中,例如实施256点的快速傅立叶变换,将获得的频率频谱振幅值用频带分割部分104例如分割成18个频带。该频带分割的输入信号频谱Y[w,k]被传送至频谱修正部分110、噪声估计部分105内的噪声频谱估计部分126和Hn值计算部分107。W是频带号码。
接着,在噪声估计部分105中,从帧化信号yframe[j,k]中区分为噪声和声音,检测噪声和估计的帧,并将估计噪声电平值和最大SN比传送至NR计算部分106。
在RMS计算部分121中,进行每帧的各信号成分的平方平均值的平方根(RMS:均方根)的计算,作为RMS值RMS[k]来输出。
在相对能量计算部分122中,计算、输出与前帧的衰减能量有关的第k帧的相对能量。
在最大RMS计算部分123中,求出后述的估计噪声电平值与信号电平和估计噪声电平之间比的最大值,即求出估算所谓的最大SN比所需要的最大RMS值,作为最大RMS值MaxRMS[k]来输出。
在估计噪声电平计算部分124中,在评价基底噪声或背景噪声电平上最合适的最小RMS值从当前帧的过去5个局部极小值(局部最小值)内的最小值来选择,作为估计噪声电平值MinRMS[k]输出。
在最大SN比计算部分125中,使用最大RMS值MaxRMS[k]和估计噪声值MinRMS[k],算出最大SN比MaxSNR[k]。
在噪声频谱估计部分126中,使用RMS值RMS[k]、相对能量、估计噪声电平MinRMS[k]、最大RMS值MaxRMS[k],算出并输出背景噪声频谱的时间平均估计值N[w,k]。
在NR值计算部分106中,算出为避免滤波器响应急剧变化使用的值NR[w,k]。
在Hn值计算部分107中,使用频带分割后的输入信号频谱Y[w,k]、噪声频谱的时间平均估计值N[w,k]和NR值计算部分106输出的NR[w,k],生成用于从输入信号中除去噪声信号的滤波器Hn[w,k]。其中,生成的滤波器Hn[w,k],如果噪声成分比声音成分大,那么表示增强抑制的响应,相反地,如果声音成分比噪声成分大,那么表示减弱抑制的响应。
在滤波器处理部分108中,将滤波器Hn[w,k]值在频率轴方向和时间轴方向上进行平滑。对频率轴方向的平滑进行中值滤波处理,对于时间轴方向的平滑来说,仅在声音区间或噪声区间进行AR平滑处理,而对于过度信号来说不进行处理。
在频带变换部分109中,对来自滤波处理部分108的频带分割的滤波值进行可适合反快速傅立叶变换处理部分的输入的内插处理。频谱修正部分110将快速傅立叶变换处理部分103的输出信号与上述内插的滤波值相乘,进行频谱修正即降低噪声成分的处理,输出噪声降低信号。
在反快速傅立叶变换处理部分111中,使用由频谱处理部分110获得的噪声降低信号,进行反快速傅立叶变换处理,输出得到的IFFT信号。在重叠加法部分112中,进行对于各帧的IFFT信号的帧边界部分的重合,将获得的输出声音信号从声音信号输出端子114输出。
以上说明的噪声抑制装置根据估计噪声信号与输入声音信号的比(估计SN比)和噪声信号电平,按照帧内的声音信号和噪声信号的分布量来控制用于除去来自输入频谱的噪声频谱的滤波特性,使用该滤波特性进行频谱抑制处理,将目的信号的失真抑制到最小限度,并且可以出色地确保大的噪声抑制量,但却存在以下问题。
由于进行依赖于估计噪声信号电平和估计SN比的控制,所以在估计噪声信号电平的估计不同的情况下,存在不能进行适当的噪声抑制,而会进行过度抑制的问题。
此外,使用估计噪声信号的抑制量控制由于从判定为过去的噪声的帧的平均频谱中来生成估计噪声信号,例如在句首等急剧变化输入声音信号电平的部分上由于在滤波控制中产生时滞,所以会进行过度的抑制,其结果,存在产生句首消失感等异常噪声的问题。
本发明用于解决这种课题,目的在于提供听感好可抑制噪声并且即使在高噪声下品质劣化小的噪声抑制装置。
本发明的噪声抑制装置从不需要的噪声信号混入目标信号的输入信号中除去噪声信号,输出目标信号,该噪声抑制装置包括:
时间/频率变换装置,在每一帧中频率分析输入信号,变换为振幅频谱和相位频谱;
类噪声分析装置,判定输入信号帧的类噪声;
噪声振幅频谱计算装置,根据所述类噪声分析装置输出的判定结果,使用该帧的输入振幅频谱算出噪声振幅频谱;
频谱校正增益计算装置,使用所述输入振幅频谱和所述噪声振幅频谱及第一预定系数算出噪声振幅频谱校正增益,同时使用相同的所述输入振幅频谱和所述噪声振幅频谱及第二预定系数算出噪声除去频谱校正增益;
频谱减法装置,从所述输入振幅频谱中减去所述噪声振幅频谱与所述校正增益计算装置输出的噪声振幅频谱校正增益的乘积,输出第一噪声除去频谱;
频谱振幅抑制装置,将所述第一噪声除去频谱乘以所述校正增益计算装置输出的除去噪声频谱校正增益,输出第二噪声除去频谱;和
频率/时间变换装置,将所述第二噪声除去频谱变换为时间轴信号。
此外,本发明的噪声抑制装置包括:频谱校正增益限制值计算装置,所述频谱校正增益计算装置使用输入振幅频谱和噪声振幅频谱,算出限制噪声振幅频谱和噪声除去频谱的校正增益的频谱校正增益限制值;和
校正增益计算装置,使用输入振幅频谱和噪声振幅频谱及频谱校正增益限制值,算出噪声振幅频谱的各振幅中进行各频率成分校正的噪声振幅频谱校正增益和噪声除去频谱的各振幅值中进行各频率成分校正的噪声除去频谱校正增益。
此外,本发明的噪声抑制装置配有频谱频带分割装置,将来自所述时间/频率变换装置的输入振幅频谱分割成多个频带,算出每个频带的平均频谱,同时将来自噪声振幅频谱计算装置的噪声振幅频谱分割成多个频带,算出每个频带的平均频谱;
频谱校正增益限制值计算装置和校正增益计算装置使用频谱频带分割装置输出的输入振幅频谱和噪声振幅频谱的每个频带的平均频谱代替输入振幅频谱和噪声振幅频谱,算出噪声振幅频谱校正增益、噪声除去频谱校正增益。
此外,本发明的噪声抑制装置包括:频谱平滑系数计算装置,根据输入信号的状态,算出输入振幅频谱和噪声振幅频谱的平滑系数;
频谱平滑装置,使用所述频谱平滑系数,在时间/频率方向平滑输入振幅频谱和噪声振幅频谱,输出平滑化的输入振幅频谱和平滑化的噪声振幅频谱;和
校正增益计算装置,使用所述平滑输入振幅频谱和平滑噪声振幅频谱算出噪声振幅频谱的各振幅值中进行各频率成分校正的噪声振幅频谱校正增益和噪声除去频谱的各振幅值中进行各频率成分校正的噪声除去频谱校正增益。
此外,本发明的噪声抑制装置配有频谱频带分割装置,将来自所述时间/频率变换装置的输入振幅频谱分割成多个频带,算出各频带的平均频谱,同时将来自噪声振幅频谱计算装置的噪声振幅频谱分割成多个频带,算出各频带的平均频谱;
频谱平滑系数计算装置使用来自频谱频带分割装置的各频带的输入振幅平均频谱和各频带的噪声振幅平均频谱算出输入振幅频谱和噪声振幅频谱的平滑系数;
频谱平滑装置使用来自频谱频带分割装置的各频带的输入振幅平均频谱和各频带的噪声振幅平均频谱算出平滑输入振幅频谱和平滑噪声振幅频谱。
此外,本发明的噪声抑制装置配有:频谱平滑系数计算装置,根据输入信号的状态,算出输入振幅频谱和噪声振幅频谱的平滑系数;和
频谱平滑装置,使用所述频谱平滑系数,在时间/频率方向上平滑输入振幅频谱和噪声振幅频谱,输出平滑输入振幅频谱和平滑噪声振幅频谱;
所述校正增益计算装置使用平滑输入振幅频谱和平滑噪声振幅频谱及频谱校正增益限制值代替输入振幅频谱和噪声振幅频谱,算出噪声振幅频谱校正增益和噪声除去频谱校正增益。
此外,本发明的噪声抑制装置配有频谱频带分割装置,将来自所述时间/频率变换装置的输入振幅频谱分割成多个频带,算出各频带的平均频谱,同时将来自噪声振幅频谱计算装置的噪声振幅频谱分割成多个频带,算出各频带的平均频谱;所述频谱平滑系数计算装置、频谱平滑装置、频谱校正增益限制值计算装置和校正增益计算装置使用来自所述频谱频带分割装置的输出代替输入振幅频谱和噪声振幅频谱。
此外,本发明的噪声抑制装置的所述频谱平滑系数计算装置根据类噪声分析装置输出的判定结果算出输入振幅频谱和噪声振幅频谱的平滑系数。
图1是表示本发明实施例1结构的方框图。
图2是表示本发明实施例2结构的方框图。
图3是表示本发明实施例3结构的方框图。
图4是表示本发明实施例4结构的方框图。
图5是表示本发明实施例5结构的方框图。
图6是表示本发明实施例6结构的方框图。
图7是噪声振幅频谱校正增益限制值的对全频带的SN比变化图。
图8是噪声除去频谱校正增益限制值对输入信号功率的变化图。
图9是噪声振幅频谱校正增益的变化图。
图10是噪声除去频谱校正增益的变化图。
图11是与噪声振幅频谱校正增益有关的听觉加权Wα的示意图。
图12是与噪声除去频谱校正增益有关的听觉加权Wβ的示意图。
图13是表示以往的噪声抑制装置结构的方框图。
实施例1
以下,参照附图说明本发明噪声抑制装置的实施例1。
图1是表示本发明噪声抑制装置的实施例1结构的方框图。1是输入信号端子,2是时间/频率变换装置,3是类噪声分析装置,4是噪声振幅频谱计算装置,5是频谱校正增益限制值计算装置,6是校正增益计算装置,7是频谱减法装置,8是频谱抑制装置,9是频率/时间变换装置,10是输出信号端子。
再有,在本实施例1中,由频谱校正增益限制值计算装置5和校正增益计算装置6构成频谱校正增益计算装置。
以下,根据图1说明本发明的噪声抑制装置的工作原理。
按预定的抽样频率(例如8kHz)来抽样,从输入信号端子1输入分割成预定帧单位(例如20ms)的输入信号s[t]。该输入信号有只有背景噪声的情况,也有背景噪声混入声音信号的情况。
时间/频率变换装置2使用例如256点FFT(快速傅立叶变换),将输入信号s[t]变换为振幅频谱S[f]和相位频谱P[f]。再有,由于FFT是公知的方法,所以省略说明。
类噪声分析装置3由线性预测分析装置14、低通滤波器11、反滤波器12、自相关分析装置13、更新速度系数决定装置15构成。
首先,用低通滤波器11获得进行输入信号滤波处理的低通滤波信号。该滤波器的截止频率例如为2kHz。通过进行低通滤波处理,可以除去高频带噪声影响并能进行稳定分析。
接着,用线性预测分析装置14进行低通滤波信号的线性预测分析,获得线性预测系数(例如10次的α参数)。反滤波器12使用线性预测系数,将低通滤波信号进行反滤波处理,输出低通线性预测残留误差信号(以下略记为低通残留误差信号)。接着,用自相关分析装置13进行低通残留误差信号的自相关分析,求出正的峰值,以该值作为RACmax。
更新速度系数决定装置15例如使用前面的RACmax、当前帧的低通残留误差信号的功率Rpow、来自输入信号端子1的当前帧的全频带功率Fpow,算出类噪声电平Nlevel,根据它算出噪声振幅频谱更新速度系数r。
类噪声电平Nlevel使用RACmax、Rpow和Fpow按照以下规则来决定。但是,RACth、Rpowth、Fpowth分别为自相关最大值的阈值、低通残留误差功率阈值、帧全频带功率阈值,各自为预定的常数。
start:Nlevel=0;;;将类噪声电平值清零if(RACmax>RACth) Nlevel=Nlevel+2if(Rpow>Rpowth) Nlevel=Nlevel+1if(Fpow>Fpowth) Nlevel=Nlevel+1output Nlevel;;;输出类噪声电平值end:
噪声振幅频谱更新速度系数r与各个类噪声电平Nlevel对应,例如象表1那样来提供。再有,r越大,在作为保持的过去噪声振幅频谱的平均的噪声振幅频谱N[f]中,当前帧的输入振幅频谱成分就被越大地反映。
【表1】
类噪声电平 | 噪声性 | 更新速度系数r |
01234 | 噪声性大...噪声性小 | 0.50.60.80.950.999 |
噪声振幅频谱计算装置4根据类噪声分析装置3输出的噪声振幅频谱更新系数r和时间/频率变换装置2输出的输入振幅频谱S[f]进行式(1)那样的噪声振幅频谱N[f]更新。Nold[f]是更新前的噪声振幅频谱,而Nnew[f]是更新后的噪声振幅频谱。以下,噪声振幅频谱N[f]指更新后的噪声振幅频谱Nnew[f]。
Nnew[f]=(1-r)·Nold[f]+r·S[f] …式1
再有,噪声振幅频谱N[f]的初期值以噪声振幅频谱更新系数r为1.0按式(1)求出。
频谱校正增益限制值计算装置5使用时间/频率变换装置2输出的输入振幅频谱S[f]和噪声振幅频谱计算装置4输出的噪声振幅频谱N[f],算出噪声振幅频谱校正增益的限制值Lα和噪声除去频谱校正增益的限制值Lβ。
首先,按照式(2)求出输入振幅频谱S[f]的功率Ps(分贝值)。
Ps(dB)=10log10(∑(S[f]·S[f])) …式2
接着,按式(3)求出噪声振幅频谱N[f]的功率Pn(分贝值)。再有,将Pn限制在PnMIN≤Pn≤0的范围内。PnMIN是表示噪声信号功率的最小值(分贝值)的预定常数。MIN(a,b)是选择两个自变量a、b中大的一个的函数。
Pn(dB)=MAX(-10log10(∑(N[f]·N[f]),PnMIN) …式3
接着,使用前面的Ps和Pn按式(4)求出当前帧全频带的输入信号与噪声信号的SN比snrall。
Snrall(dB)=Ps+Pn …式4
接着,使用按式(4)求出的全频带SN比snrall,根据式(5)决定并输出噪声振幅频谱校正增益限制值Lα。式(5)中的αMIN是噪声振幅频谱校正增益的最大值(dB),αMIN是噪声振幅频谱校正增益的最小值(dB),分别为预定的常数。此外,SNR1、SNRn是与全频带SN比有关的阈值,都为预定的常数。再有,Lα是后述的频谱减法装置7中从输入振幅频谱中减去噪声振幅频谱的抽取率的最大值极限。图7表示式(5)的Lα对snrall的变化的大致情况。
接着,使用式(6),算出输入信号功率Ps和阈值Psth之间的差dPs。再有,Psth是作为输入信号功率阈值的预定常数。
dPs(dB)=Ps-Psth …式6
在求出输入信号功率和阈值之间的差dPs后,使用式7,决定并输出噪声除去频谱校正增益β[f]的限制值Lβ。再有,Lβ是后述的频谱抑制装置中的振幅抑制量的最大值的极限。图8表示Ps与式(7)的Lβ的变化的大致情况。
校正增益计算装置6使用输入振幅频谱S[f]、噪声振幅频谱N[f]、噪声振幅频谱校正增益限制值Lα和噪声除去频谱校正增益限制值Lβ算出用于进行对噪声振幅频谱N[f]中每个频率成分的校正的噪声频谱校正增益α[f]、以及用于进行对后述的第一噪声除去频谱Ss[t]中每个频率成分校正的噪声除去频谱校正增益α[f]。
首先,根据式(8)求出输入振幅频谱和噪声振幅频谱的各频率成分的SN比snrsp[f],fn是乃奎斯特频率。
根据式(9),使用按式(8)求出的每个频率的SN比snrsp[f]、噪声功率的最小值PnMIN、噪声振幅频谱校正增益限制值Lα、听觉加权Wα算出噪声振幅频谱校正增益α[f]。式(9)中的噪声功率的最小值PnMIN是预定的常数。此外,MIN(a,b)是两个自变量a、b中返回小的值的函数。
gainα=MIN(snrsp[f]·Wα[f]+Ph,0)
α[f]=Lα·{(PnMIN+gainα)/PnMIN} …式9
通过使用式(9),随着snrsp[f]变大即每个频率的SN比变高,gainα变大,其结果,由于噪声振幅频谱校正增益α[f]变大,所以在后述的频谱减法装置7中,SN比高的频谱成分使从输入信号频谱中抽取噪声频谱的比率变大,而SN比低的频谱成分使抽取率变小。图9表示α[f]对snrsp[f]的变化的大概情况。
再有,听觉加权Wα[f]是以频率f为参数的预定加权系数,频率越高,其值越大。通过进行这种加权,因在高频区域中α[f]的值变小而防止高频区域中的过度抑制,防止发生异常声音。图11表示Wα[f]的大概情况。
接着,根据式(10),使用输入振幅频谱S[f]、噪声振幅频谱N[f]、听觉加权Wβ[f]和噪声除去频谱校正增益限制值Lβ算出用于校正第二噪声除去频谱Sr[f]的各振幅值的噪声除去频谱校正增益β[f]。
gainβ=MIN(snrsp[f]·Wβ[f]+Lβ,0)
通过使用式(10),随着snrsp[f]变大即SN比变高,gainβ变小,其结果,由于噪声除去频谱校正增益β[f]变大,所以在后述的频谱抑制装置8中,SN比高的频谱成分使输出的噪声除去频谱振幅变大,而SN比低的频谱成分使输出振幅变小。图10表示β[f]对snrsp[f]变化的大概情况。
再有,听觉加权Wβ与前面的听觉加权Wα[f]一样,是以频率f为参数的预定加权系数,频率越高,其值越大。通过进行这种加权,因在高频区域中β[f]的值变大而防止高频区域中的过度抑制,防止发生异常声音。图12表示Wβ[f]的大概情况。
频谱减法装置7从式(11)那样的输入振幅频谱S[f]中减去噪声振幅频谱N[f]与校正增益计算装置6输出的噪声振幅频谱校正增益α[f]的乘积,输出第一噪声除去频谱Ss[f]。此外,在第一噪声除去频谱Ss[f]变负的情况下,进行***0或预定的低电平噪声n[f]回填处理。通过将校正增益α[f]与噪声频谱相乘,在SN比低的情况下可以抑制噪声频谱成分的抽取,而在SN比高的情况下,通过加强抽取噪声频谱成分,可以抑制低SN时的过度频谱抽取。
频谱抑制装置8使用式(12),将所述第一噪声除去频谱Ss[f]乘以校正增益计算装置6输出的噪声除去频谱校正增益β[f]后,输出第二噪声除去频谱Sr[f]。通过将校正增益β[f]与第一噪声除去频谱相乘,具有抑制上述频谱减法装置7中频谱抽取残留的残留噪声和频谱减法的结果产生的人工噪声(音乐噪声)的效果,此外,由于在低SN时可减弱振幅抑制,而在高SN时可增强振幅抑制,所以可以防止低SN时的过度振幅抑制。
Sr[f]=β[f]·Ss[f] …式12
频率/时间变换装置9采用时间/频率变换装置2的相反程序,例如进行逆FFT,使用第二噪声除去频谱Sr[f]和相位频谱P[f]变换成时间信号sr[t],进行与前帧的边界部分的信号重合,由输出信号端子10输出噪声抑制信号。
通过将噪声振幅频谱校正增益α[f]与噪声振幅频谱相乘,在SN比低的情况下抑制噪声振幅频谱成分的消失,而在SN比高的情况下通过加强减去噪声频谱成分,可以抑制低SN时过度的频谱抽取,此外,通过将噪声除去频谱校正增益β[f]与第一噪声除去频谱相乘,具有抑制频谱减法中频谱抽取残留的残留噪声和频谱减法的结果产生的人工噪声(音乐噪声)的效果。
此外,由于在低SN时可减弱振幅抑制,而在高SN时可增强振幅抑制,所以可以防止低SN时的过度振幅抑制,而且,即使在句首等输入声音信号电平急剧变化的部分中,由于进行不仅与噪声信号电平而且与输入信号电平对应的频谱减法和频谱振幅抑制处理,所以同时可以防止因过度的频谱减法和抑制处理造成的句首消失感和频谱变形感等的发生,极大地原样保持噪声区间的噪声抑制量,防止声音区间中过度的频谱抑制,可以进行良好的噪声抑制。
实施例2
下面,参照附图说明本发明噪声抑制装置的实施例2。
图2是表示实施例2结构的方框图。与图1相比,在结构上除去了频谱校正增益限制值计算装置5,而新设置了频谱平滑系数计算装置21和频谱平滑装置22。由于对于其它结构来说与实施例1相同,所以省略说明。以下,按照图2说明实施例2的工作原理。
例如,频谱平滑系数计算装置21根据类噪声判定装置3输出的输入信号的类噪声判定结果的电平算出对频谱的时间轴方向的平滑使用的频谱时间方向平滑系数γt、对频谱的频率方向的平滑使用的频谱频率方向平滑系数γf。
作为按照类噪声算出平滑系数的方法,例如,可以使用表参照表2所示的与类噪声对应的平滑系数。由此,在噪声性大的情况下,即在噪声区间中选择可增强平滑的平滑系数γt、γf,在噪声性小时,即在声音部分可以选择可减弱平滑的平滑系数γt、γf并进行输出。
表2
类噪声电平 | 类噪声 | 平滑系数γt | 平滑系数γf |
01234 | 噪声性大...噪声性小 | 0.50.60.70.80.9 | 0.70.80.850.90.95 |
频谱平滑装置22按照式(13)、式(14),使用时间方向平滑系数γt、频率方向平滑系数γf在时间方向和频率方向上平滑输入振幅频谱S[f]和噪声振幅频谱N[f],算出平滑输入振幅频谱Ssm[f]和平滑噪声振幅频谱Nsm[f]。
首先,使用式(13)进行输入振幅频谱Ssm[f]和噪声振幅频谱Nsm[f]的时间方向的平滑,算出时间方向平滑输入振幅频谱St[f]、时间方向平滑噪声振幅频谱Nt[f]。式(13)中的Spre[f]、Npre[f]分别是前帧的输入振幅频谱、噪声振幅频谱。Fn是乃奎斯特(Nyquist)频率。
St[f]=γt·S[f]+(1-γt)·Spre[f],f=0,…,fn
Nt[f]=γt·N[f]+(1-γt)·Npre[f],f=0,…,fn
…式(13)
接着,使用式(14)进行由式(13)得到的时间方向平滑输入振幅频谱St[f]、时间方向平滑噪声振幅频谱Nt[f]的频率方向的平滑,算出作为频谱平滑装置输出结果的平滑输入振幅频谱Ssm[f]、平滑噪声振幅频谱Nsm[f]。
Ssm[f]=γf·St[f]+(1-γf)·St[f-1],f=1,…,fn
Nsm[f]=γf·Nt[f]+(1-γf)·Nt[f-1],f=1,…,fn
…式(14)
校正增益计算装置6使用平滑输入振幅频谱Ssm[f]和平滑噪声振幅频谱Nsm[f]代替输入振幅频谱S[f]和噪声振幅频谱N[f],算出噪声振幅频谱校正增益α[f]和噪声除去频谱校正增益β[f]。
接着,使用平滑SN比snrsp-sm[f],根据式(16)和式(17),算出平滑噪声振幅频谱校正增益αsm[f]和平滑噪声除去频谱校正增益βsm[f]。
gainα=MIN(snrsp-sm[f]·Wα[f]+Pn,0)
αsm[f]=αMAX·{(PnMIN+gainα)/PnMIN) …式16
gainβ=MIN(snrsp-sm[f]·Wβ[f]+Pn(=βMIN),0)
βsm[f]=10(gain,/20)
…式17
通过使用平滑SN比snrsm[f]求出校正增益,输入声音的电平比即SN比对噪声信号小的噪声区间中,极大地抑制频谱校正增益的变动,而在声音区间那样的SN比高的区间中,成为不大抑制校正增益变动的动作。
再有,在式(16)和式(17)中,与上述实施例1中式(9)和式(10)不同的部分在于未使用噪声振幅频谱校正增益限制值Lα和噪声除去频谱校正增益限制值Lβ。此外,式中的αMAX为噪声振幅频谱校正增益最大值,βMIN为噪声除去频谱校正增益最小值(βMIN=Pn),分别为预定的常数。
通过进行与类噪声电平对应的频谱平滑系数的控制,由于在噪声性大的情况下增强平滑,在噪声性小即声音部分中选择可减弱平滑的平滑系数,在噪声性大即在噪声区间中选择可增强平滑的平滑系数,所以可以进一步进行频谱校正增益的适当控制,可以进行良好的噪声控制。
上述噪声除去频谱的不连续感缓和效果在高噪声电平等频谱校正增益的精度变低的低SN时特别有效。
实施例3
作为实施例1的另一形态,在实施例1的结构中可以导入用实施例2说明的频谱平滑处理。图3是表示本实施例3的结构的方框图。
频谱校正增益限制值计算装置5按照用实施例2说明的步骤,使用频谱平滑装置22算出的平滑输入振幅频谱Ssm[f]和平滑噪声振幅频谱Nsm[f],按与实施例1相同的步骤,算出噪声振幅频谱校正增益的限制值Lα和噪声除去频谱校正增益的限制值Lβ。
校正增益计算装置6使用来自频谱平滑装置22的平滑输入振幅频谱Ssm[f]和平滑噪声振幅频谱Nsm[f],以及来自频谱校正增益限制值计算装置5的噪声振幅频谱校正增益限制值Lα和噪声除去频谱校正增益限制值Lβ,按照与上述实施例1相同的式(9)和式(10),算出噪声振幅频谱校正增益α[f]和噪声除去频谱校正增益β[f]。
由于其它结构与在实施例1和实施例2中说明的结构相同,所以省略说明。
通过获得本实施例的结构,除了实施例1的效果以外,利用实施例2的相乘效果,可以进行更合适的噪声抑制。
实施例4
作为根据输入声音的状态算出频谱平滑系数的方法,例如,也可以使用当前帧的SN比。图4是表示本实施例4结构的方框图。
接着,按照式19,使用帧的SN比SNRfr,求出对频谱的时间方向平滑使用的频谱时间方向平滑系数的伪系数γt’、对频谱的频率方向平滑使用的频谱频率方向平滑系数的伪系数γf’。
接着,按照式(20),使用前帧的平滑系数γ(old)t、γ(old)f将伪平滑系数γt’、γf’进行AR平滑,输出频谱时间方向平滑系数γt和频谱频率方向平滑系数γf。
γt=0.8·γt’+0.2·γ(old)t
γf=0.8·γf’+0.2·γ(old)f …式20
如本实施例所示,通过使用与输入信号的SN比对应的频谱平滑系数来平滑输入振幅频谱和噪声振幅频谱,采用使用它们算出的频谱校正增益来进行噪声抑制处理,可以根据输入信号的SN比来控制频谱校正增益的变动,例如,在噪声区间等的低SN时,由于可以缓和噪声除去频谱的时间方向和频率方向的不连续感,所以可以抑制输出声音发生异常声音,可以进行稳定的噪声抑制。
实施例5
作为发明的实施例1的另一形态,将输入振幅频谱不仅分割成其每个频率成分,还分割成多个频带,使用每个频带的平均频谱,算出噪声振幅频谱校正增益和噪声除去频谱校正增益,利用这些增益可以进行各自的频谱校正。
在本实施例5中,在频谱校正增益限制值计算装置5的前级中,配有频谱频带分割装置,将来自所述时间/频率变换装置的输入振幅频谱分割成多个频带,算出各频带的平均频谱,同时将来自噪声振幅频谱计算装置的噪声振幅频谱分割成多个频带,算出各频带的平均频谱。
频谱频带分割装置将输入振幅频谱分割成例如16个频带,如式(21)所示,对于各个频带(称为信道ch)求出输入信号的平均频谱Save[ch]和噪声信号的平均频谱Nave[ch]。nch是信道ch中的频谱成分数。
接着,频谱校正增益限制值计算装置5使用按式(21)求出的不同信道的平均频谱Save[ch]、Nave[ch],按照式(22)求出输入信号功率Psave和噪声信号功率Pnave,求出全频带SN比snrall-ave。PnMIN是最低噪声功率的预定常数。
Psave(dB)=10log10(∑Save[ch]·Save[ch])
Pnave(dB)=MAX(-10log10(∑Nave[ch]·Nave[ch],PnMIN)
snrall-ave=Psave+Pnave …式22
接着,代替上述发明的实施例1中的Ps、Pn,使用上述得到的输入信号功率Psave、噪声信号功率Pnave,算出噪声振幅频谱校正增益限制值Lα、噪声除去频谱校正增益限制值Lβ。
在校正增益计算装置6中,按式(23)求出各信道的SN比snrsp[ch],利用它算出各信道中噪声振幅频谱校正增益α[ch]、噪声除去频谱校正增益β[ch]。Nch是信道总数。
在频谱减法装置7、频谱抑制装置8中,根据各自输入的校正增益,展开与频谱各成分对应的值,进行频谱减法和频谱振幅抑制。
通过获得本实施例的结构,除了发明的实施例1具有的效果以外,还具有可以消减进行频谱校正增益计算的运算量,以及存储频谱校正增益的存储量的效果。
实施例6
作为实施例4的另一形态,不仅将输入振幅频谱分割成其每个频率成分,还分割成多个频带,使用每个频带的平均频谱,也可以算出频谱平滑系数。图5是本实施例6的结构图。
在图5中,23是频谱频带分割装置,将来自所述时间/频率变换装置的输入振幅频谱分割成多个频带,算出各频带的平均频谱,同时将来自噪声振幅频谱计算装置的噪声振幅频谱分割成多个频带,算出各频带的平均频谱。
频谱频带分割装置23与上述式(21)一样将输入振幅频谱分割成例如16个频带,对于各个频带(称为信道ch)求出输入信号的平均频谱Save[ch]和噪声信号的平均频谱Nave[ch]。
接着,使用上述发明的实施例2中的式(14)和式(15),代替上述帧的SN比SNRfr,采用使用平均频谱算出的帧的SN比SNRfr-ave,输出频谱时间方向平滑系数γt和频谱频率方向平滑系数γf。
频谱平滑装置22按照式(25)、式(26),使用由上述平均频谱得到的时间方向平滑系数γt、频率方向平滑系数γf在时间方向和频率方向上平滑输入信号的平均频谱Save[ch]和噪声信号的平均频谱Nave[ch],算出平滑输入平均频谱Ssm-ave[ch]和平滑噪声平均频谱Nsm- ave[ch]。
首先,使用式(25),进行输入信号的平均频谱Save[ch]和噪声信号的平均频谱Nave[ch]的时间方向的平滑,算出时间方向平滑输入信号的平均频谱St-ave[ch]、时间方向平滑噪声信号的平均频谱Nt- ave[ch]。式(25)中的Spre-ave[ch]、Npre-ave[ch]分别是前帧的输入信号平均频谱、噪声信号平均频谱。Nch是最大信道数。
St-ave[ch]=γt·Save[ch]+(1-γt)·Spre-ave[ch],ch=0,…Nch
Nt-ave[ch]=γt·Nave[ch]+(1-γt)·Npre-ave[ch],ch=0,…Nch
…式25
接着,使用式(26),进行由式(25)得到的时间方向平滑输入信号的平均频谱St-ave[ch]、时间方向平滑噪声信号的平均频谱Nt- ave[ch]的频率方向的平滑,算出作为频谱平滑装置的输出结果的平滑输入振幅频谱Ssm-ave[ch]、平滑噪声振幅频谱Nsm-ave[ch]。
Ssm-ave[ch]=γf·St-ave[ch]+(1-γf)·St-ave[ch-1],ch=0,…,Nch
Nsm-ave[ch]=γf·Nt-ave[ch]+(1-γf)·Nt-ave[ch-1],ch=0,…,Nch
…式26
校正增益计算装置6使用平滑输入信号的平均频谱Ssm-ave[ch]和平滑噪声振幅频谱Nsm-ave[ch]代替平滑输入振幅频谱Ssm[f]和平滑噪声振幅频谱Nsm[f],算出每个信道的噪声振幅频谱校正增益α[ch]和噪声除去频谱校正增益β[ch]。
首先,使用式(27),采用平滑输入信号的平均频谱Ssm-ave[ch]和平滑噪声振幅频谱Nsm-ave[ch],求出各信道的平滑SN比snrsm- ave[f]。
接着,使用平滑SN比snrch-sm[ch],根据式(28)和式(29),算出平滑噪声振幅频谱校正增益αsm[ch]和平滑噪声除去频谱校正增益βsm[ch]。
gainα=MIN(snrch-sm[ch]·Wα[ch]+Pn,0)
αsm[ch]=αMAX·{(PnMIN+gainα)/PnMIN} …式28
gainβ=MIN(snrch-sm[ch]·Wβ[ch]+Pn(=βMIN,0)
βsrn[ch]=10(gainβ/20)
…式29
使用以上得到的平滑噪声振幅频谱校正增益αsm[ch]和平滑噪声除去频谱校正增益βsm[ch],进行频谱的减法和频谱抑制处理。
通过形成本实施例的结构,除了发明的实施例2具有的效果以外,还具有可以消减进行频谱平滑系数的计算和频谱平滑的运算量,以及存储频谱平滑系数的存储量的效果。
实施例7
作为实施例3的另一形态,也可以有组合上述实施例5和实施例6结构的结构。图6是本实施例7的结构图。
频谱频带分割装置23与实施例6一样将输入振幅频谱分割成多个频带,算出各频带的平均频谱,同时将来自噪声振幅频谱计算装置的噪声振幅频谱分割成多个频带,算出各频带的平均频谱。
频谱平滑装置22使用从频谱平滑系数计算装置21得到的时间方向平滑系数γt、频率方向平滑系数γf在时间方向和频率方向上平滑输入信号的各频带的平均频谱Save[ch]和噪声信号的各频带的平均频谱Nave[ch],算出平滑输入平均频谱Ssm-ave[ch]和平滑噪声平均频谱Nsm- ave[ch]。
接着,频谱校正增益限制值计算装置5使用平滑输入平均频谱Ssm-ave[ch]和平滑噪声平均频谱Nsm-ave[ch],根据式(22)求出输入信号功率Psave和噪声信号功率Pnave,求出全频带SN比snrall-ave。PnMIN是最低噪声功率的预定常数。
接着,代替上述发明的实施例1中的Ps、Pn,使用上述获得的输入信号功率Psave、噪声信号功率Pnave,算出噪声振幅频谱校正增益限制值Lα、噪声除去频谱校正增益限制值Lβ。
在校正增益计算装置6中,按式(23)求出各信道的SN比snrsp[ch],用它算出各信道中噪声振幅频谱校正增益α[ch]、噪声除去频谱校正增益β[ch]。Nch是信道总数。
由于其它结构与实施例5和实施例6中说明的结构相同,所以省略说明。
通过形成本实施例的结构,除了实施例3具有的效果以外,还具有可以消减算出频谱校正增益、频谱平滑系数和用于进行频谱平滑的运算量,以及存储频谱校正增益、频谱平滑系数的存储量的效果。
按照如上所述的本发明的噪声抑制装置,由于根据输入信号帧的类噪声的判定结果,使用该帧的输入振幅频谱算出的噪声振幅频谱和使用输入振幅频谱和各个预定系数算出噪声振幅频谱校正增益及噪声除去频谱校正增益,从所述输入振幅频谱中减去所述噪声振幅频谱与所述校正增益计算装置输出的噪声振幅频谱校正增益的乘积,将输出的第一噪声除去频谱与所述校正增益计算装置输出的噪声除去频谱校正增益相乘,输出第二噪声除去频谱,将所述第二噪声除去频谱变换成时间轴信号,所以即使在句首等输入声音信号电平急剧变化的部分中,由于进行不仅与噪声信号电平而且与输入信号电平对应的频谱减法和频谱振幅抑制处理,所以可以防止因过度的频谱减法和抑制处理造成的句首消失感和频谱变形感等的发生,极大地原样保持噪声区间的噪声抑制量,防止声音区间中过度的频谱抑制,可以进行良好的噪声抑制。
此外,通过使噪声除去频谱校正增益与第一噪声除去频谱相乘,具有抑制频谱减法中抽取残留频谱的残留噪声和频谱减法的结果产生的人工噪声(音乐噪声)的效果。
此外,按照本发明的噪声抑制装置,使用与输入信号状态对应的输入振幅频谱和噪声振幅频谱的平滑系数,在时间/频率方向上平滑输入振幅频谱和噪声振幅频谱,算出平滑输入振幅频谱和平滑噪声振幅频谱,算出使用该平滑输入振幅频谱和平滑噪声振幅频谱的噪声振幅频谱校正增益和噪声除去频谱校正增益,通过进行与类噪声电平对应的频谱平滑系数的控制,由于在噪声性大的情况下增强平滑,在噪声性小即在声音部分中选择使平滑减弱的平滑系数,在噪声性大即在噪声区间选择可增强平滑的平滑系数,所以可以进一步进行频谱校正增益的适当控制,可以进行良好的噪声抑制。
此外,按照本发明的噪声抑制装置,配有频谱频带分割装置,将输入振幅频谱分割成多个频带,算出各频带的平均频谱,同时将噪声振幅频谱分割成多个频带,算出各频带的平均频谱,通过将各频带的平均频谱使用于频谱的平滑系数和平滑频谱的计算出,可以防止因过度的频谱减法和抑制处理造成的句首消失感和频谱变形感等的发生,极大地原样保持噪声区间的噪声抑制量,防止声音区间中过度的频谱抑制,可以进行良好的噪声抑制。
此外,由于通过进行与类噪声电平对应的频谱平滑系数的控制,可以在噪声性大的情况下增强平滑,在噪声性小即在声音部分中选择使平滑减弱的平滑系数,在噪声性大即在噪声区间选择可增强平滑的平滑系数,所以可以进一步进行频谱校正增益的适当控制,可以进行良好的噪声抑制。
此外,按照本发明的噪声抑制装置,通过使用与输入信号状态对应的频谱平滑系数,将输入振幅频谱和噪声振幅频谱平滑,使用用这些频谱算出的频谱校正增益来进行噪声抑制处理,可以根据输入信号的状态来控制频谱校正增益的变动,例如,在噪声区间等低SN时,由于可以缓和噪声除去频谱的时间方向和频率方向的不连续感,所以可以抑制输出声音的异常声音的发生,可以进行稳定的噪声抑制。
此外,按照本发明的噪声抑制装置,由于根据输入信号的状态使用输入振幅频谱和噪声振幅频谱的平滑系数,在时间/频率方向上平滑输入振幅频谱和噪声振幅频谱,算出平滑输入振幅频谱和平滑噪声振幅频谱,使用该平滑输入振幅频谱和平滑噪声振幅频谱以及频谱校正增益限制值,算出噪声振幅频谱校正增益和噪声除去频谱校正增益,所以可以防止因过度的频谱减法和抑制处理造成的句首消失感和频谱变形感等的发生,极大地原样保持噪声区间的噪声抑制量,防止声音区间中过度的频谱抑制,可以达到良好的噪声抑制,此外,具有可以消减进行频谱校正增益计算的运算量,以及可以消减存储频谱校正增益的存储量的效果。
此外,按照本发明的噪声抑制装置,由于将输入振幅频谱分割成多个频带,算出各频带的平均频谱,同时将噪声振幅频谱分割成多个频带,使用各频带的平均频谱算出输入振幅频谱和噪声振幅频谱的平滑系数,使用该各频带的输入振幅平均频谱和各频带的噪声振幅平均频谱,算出平滑输入振幅频谱和平滑噪声振幅频谱,由于可以进行与类噪声电平对应的频谱平滑系数的控制,可以在噪声性大的情况下增强平滑,在噪声性小即在声音部分中选择使平滑减弱的平滑系数,在噪声性大即在噪声区间选择增强平滑的平滑系数,所以可以进一步进行频谱校正增益的适当控制,可以进行良好的噪声抑制。
此外,具有可以消减进行频谱平滑系数的计算和频谱平滑的运算量,以及可以消减存储频谱平滑系数存储量的效果。
此外,按照本发明的噪声抑制装置,所述频谱平滑系数计算装置、频谱平滑装置、频谱校正增益限制值计算装置和校正增益计算装置通过将输入振幅频谱分割成多个频带,使用各频带的平均频谱,和将噪声振幅频谱分割成多个频带,使用各频带的平均频谱,代替输入振幅频谱和噪声振幅频谱,可以防止因过度的频谱减法和抑制处理造成的句首消失感和频谱变形感等的发生,极大地原样保持噪声区间的噪声抑制量,防止声音区间中过度的频谱抑制,可以进行良好的噪声抑制。
此外,由于通过进行与类噪声电平对应的频谱平滑系数的控制,可以在噪声性大的情况下增强平滑,在噪声性小即在声音部分中选择使平滑减弱的平滑系数,在噪声性大即在噪声区间选择增强平滑的平滑系数,所以可以进一步进行频谱校正增益的适当控制,可以进行良好的噪声抑制。
此外,具有可以消减算出频谱校正增益、频谱平滑系数的运算量和进行频谱平滑的运算量,以及可以消减存储频谱校正增益、频谱平滑系数的存储量的效果。
Claims (8)
1.一种噪声抑制装置,从不需要的噪声信号混入目标信号的输入信号中除去噪声信号,输出目标信号,其特征在于,该噪声抑制装置包括:
时间/频率变换装置,在每一帧中频率分析输入信号,变换成振幅频谱和相位频谱;
类噪声分析装置,判定输入信号帧的类噪声;
噪声振幅频谱计算装置,根据所述类噪声分析装置输出的判定结果,使用该帧的输入振幅频谱算出噪声振幅频谱;
频谱校正增益计算装置,使用所述输入振幅频谱和所述噪声振幅频谱及第一预定系数算出噪声振幅频谱校正增益,同时使用相同的所述输入振幅频谱和所述噪声振幅频谱及第二预定系数算出噪声除去频谱校正增益;
频谱减法装置,从所述输入振幅频谱中减去所述噪声振幅频谱与所述校正增益计算装置输出的噪声振幅频谱校正增益的乘积,输出第一噪声除去频谱;
频谱振幅抑制装置,将所述第一噪声除去频谱乘以所述校正增益计算装置输出的噪声除去频谱校正增益,输出第二噪声除去频谱;和
频率/时间变换装置,将所述第二噪声除去频谱变换为时间轴信号。
2.如权利要求1所述的噪声抑制装置,其特征在于,包括:
频谱校正增益限制值计算装置,所述频谱校正增益计算装置使用输入振幅频谱和噪声振幅频谱,算出限制噪声振幅频谱和噪声除去频谱的校正增益的频谱校正增益限制值;和
校正增益计算装置,使用输入振幅频谱和噪声振幅频谱及频谱校正增益限制值,算出在噪声振幅频谱的各振幅中进行各频率成分校正的噪声振幅频谱校正增益和在噪声除去频谱的各振幅值中进行各频率成分校正的噪声除去频谱校正增益。
3.如权利要求2所述的噪声抑制装置,其特征在于包括:
频谱频带分割装置,将来自所述时间/频率变换装置的输入振幅频谱分割成多个频带,算出每个频带的平均频谱,同时将来自噪声振幅频谱计算装置的噪声振幅频谱分割成多个频带,算出每个频带的平均频谱;
形成频谱校正增益计算装置的频谱校正增益限制值计算装置和校正增益计算装置使用频谱频带分割装置输出的输入振幅频谱和噪声振幅频谱的每个频带的平均频谱代替输入振幅频谱和噪声振幅频谱,算出频谱振幅限制值、噪声振幅频谱校正增益、噪声除去频谱校正增益。
4.如权利要求1所述的噪声抑制装置,其特征在于包括:
频谱平滑系数计算装置,根据输入信号的状态,算出输入振幅频谱和噪声振幅频谱的平滑系数;
频谱平滑装置,使用所述频谱平滑系数,在时间/频率方向平滑输入振幅频谱和噪声振幅频谱,输出平滑输入振幅频谱和平滑噪声振幅频谱;
所述频谱校正增益计算装置配有校正增益计算装置,使用所述平滑输入振幅频谱和平滑噪声振幅频谱算出在噪声振幅频谱的各振幅值中进行各频率成分校正的噪声振幅频谱校正增益和在噪声除去频谱的各振幅值中进行各频率成分的校正的噪声除去频谱校正增益。
5.如权利要求4所述的噪声抑制装置,其特征在于配有频谱频带分割装置,将来自所述时间/频率变换装置的输入振幅频谱分割成多个频带,算出各频带的平均频谱,同时将来自噪声振幅频谱计算装置的噪声振幅频谱分割成多个频带,算出各频带的平均频谱;
频谱平滑系数计算装置使用来自频谱频带分割装置的各频带的输入振幅平均频谱和各频带的噪声振幅平均频谱算出输入振幅频谱和噪声振幅频谱的平滑系数;
频谱平滑装置使用来自频谱频带分割装置的各频带的输入振幅平均频谱和各频带的噪声振幅平均频谱算出平滑输入振幅频谱和平滑噪声振幅频谱。
6.如权利要求2所述的噪声抑制装置,其特征在于配有:
频谱平滑系数计算装置,根据输入信号的状态,算出输入振幅频谱和噪声振幅频谱的平滑系数;和
频谱平滑装置,使用所述频谱平滑系数,在时间/频率方向上平滑输入振幅频谱和噪声振幅频谱,输出平滑输入振幅频谱和平滑噪声振幅频谱;
所述校正增益计算装置使用平滑输入振幅频谱和平滑噪声振幅频谱及频谱校正增益限制值代替输入振幅频谱和噪声振幅频谱,算出噪声振幅频谱校正增益和噪声除去频谱校正增益。
7.如权利要求6所述的噪声抑制装置,其特征在于配有频谱频带分割装置,将来自所述时间/频率变换装置的输入振幅频谱分割成多个频带,算出各频带的平均频谱,同时将来自噪声振幅频谱计算装置的噪声振幅频谱分割成多个频带,算出各频带的平均频谱;
所述频谱平滑系数计算装置、频谱平滑装置、频谱校正增益限制值计算装置和校正增益计算装置使用来自所述频谱频带分割装置的输出来代替输入振幅频谱和噪声振幅频谱。
8.如权利要求4至权利要求7中任何一项所述的噪声抑制装置,其特征在于,所述频谱平滑系数计算装置根据类噪声分析装置输出的判定结果算出输入振幅频谱和噪声振幅频谱的平滑系数。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP31912699A JP3454206B2 (ja) | 1999-11-10 | 1999-11-10 | 雑音抑圧装置及び雑音抑圧方法 |
JP319126/1999 | 1999-11-10 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1296258A true CN1296258A (zh) | 2001-05-23 |
CN1192360C CN1192360C (zh) | 2005-03-09 |
Family
ID=18106764
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB001221361A Expired - Fee Related CN1192360C (zh) | 1999-11-10 | 2000-07-31 | 噪声抑制装置 |
Country Status (6)
Country | Link |
---|---|
US (1) | US7158932B1 (zh) |
EP (1) | EP1100077B1 (zh) |
JP (1) | JP3454206B2 (zh) |
CN (1) | CN1192360C (zh) |
DE (1) | DE60040895D1 (zh) |
HK (1) | HK1037052A1 (zh) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1302462C (zh) * | 2003-12-03 | 2007-02-28 | 富士通株式会社 | 降噪装置和降噪方法 |
CN1308914C (zh) * | 2001-06-06 | 2007-04-04 | 三菱电机株式会社 | 噪声抑制装置和方法 |
WO2008067735A1 (fr) * | 2006-12-05 | 2008-06-12 | Huawei Technologies Co., Ltd. | Procédé et dispositif de classement pour un signal sonore |
CN1892822B (zh) * | 2005-05-31 | 2010-06-09 | 日本电气株式会社 | 噪声抑制方法和设备 |
CN102270454A (zh) * | 2010-06-07 | 2011-12-07 | 宇达电脑(上海)有限公司 | 改善音频输出的方法及装置 |
CN101237814B (zh) * | 2005-06-10 | 2012-11-14 | 纽诺莫尼可斯有限公司 | 用于频谱地调整数字音频信号的数字回放设备和方法及装置 |
CN102906813A (zh) * | 2010-05-24 | 2013-01-30 | 日本电气株式会社 | 信号处理方法、信息处理装置和信号处理程序 |
CN103718241A (zh) * | 2011-11-02 | 2014-04-09 | 三菱电机株式会社 | 噪音抑制装置 |
CN105103230A (zh) * | 2013-04-11 | 2015-11-25 | 日本电气株式会社 | 信号处理装置、信号处理方法、信号处理程序 |
CN105445549A (zh) * | 2015-12-11 | 2016-03-30 | 中国电子科技集团公司第四十一研究所 | 一种微弱电磁信号频谱测量方法 |
CN111147983A (zh) * | 2018-11-06 | 2020-05-12 | 展讯通信(上海)有限公司 | 扬声器控制方法及装置、可读存储介质 |
CN111312277A (zh) * | 2014-03-03 | 2020-06-19 | 三星电子株式会社 | 用于带宽扩展的高频解码的方法及设备 |
US11688406B2 (en) | 2014-03-24 | 2023-06-27 | Samsung Electronics Co., Ltd. | High-band encoding method and device, and high-band decoding method and device |
Families Citing this family (55)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FI19992453A (fi) | 1999-11-15 | 2001-05-16 | Nokia Mobile Phones Ltd | Kohinanvaimennus |
JP4282227B2 (ja) * | 2000-12-28 | 2009-06-17 | 日本電気株式会社 | ノイズ除去の方法及び装置 |
US7103539B2 (en) * | 2001-11-08 | 2006-09-05 | Global Ip Sound Europe Ab | Enhanced coded speech |
US8271279B2 (en) | 2003-02-21 | 2012-09-18 | Qnx Software Systems Limited | Signature noise removal |
US8326621B2 (en) | 2003-02-21 | 2012-12-04 | Qnx Software Systems Limited | Repetitive transient noise removal |
US7725315B2 (en) * | 2003-02-21 | 2010-05-25 | Qnx Software Systems (Wavemakers), Inc. | Minimization of transient noises in a voice signal |
US7885420B2 (en) * | 2003-02-21 | 2011-02-08 | Qnx Software Systems Co. | Wind noise suppression system |
US7895036B2 (en) * | 2003-02-21 | 2011-02-22 | Qnx Software Systems Co. | System for suppressing wind noise |
US7949522B2 (en) | 2003-02-21 | 2011-05-24 | Qnx Software Systems Co. | System for suppressing rain noise |
US8073689B2 (en) * | 2003-02-21 | 2011-12-06 | Qnx Software Systems Co. | Repetitive transient noise removal |
JP2004341339A (ja) * | 2003-05-16 | 2004-12-02 | Mitsubishi Electric Corp | 雑音抑圧装置 |
JP4608650B2 (ja) * | 2003-05-30 | 2011-01-12 | 独立行政法人産業技術総合研究所 | 既知音響信号除去方法及び装置 |
WO2005070011A2 (en) * | 2004-01-23 | 2005-08-04 | Advanced Cell Technology, Inc. | Improved modalities for the treatment of degenerative diseases of the retina |
US7725314B2 (en) * | 2004-02-16 | 2010-05-25 | Microsoft Corporation | Method and apparatus for constructing a speech filter using estimates of clean speech and noise |
JP4542399B2 (ja) * | 2004-09-15 | 2010-09-15 | 日本放送協会 | 音声スペクトル推定装置および音声スペクトル推定プログラム |
EP1806739B1 (en) * | 2004-10-28 | 2012-08-15 | Fujitsu Ltd. | Noise suppressor |
US7957964B2 (en) | 2004-12-28 | 2011-06-07 | Pioneer Corporation | Apparatus and methods for noise suppression in sound signals |
KR100657948B1 (ko) * | 2005-02-03 | 2006-12-14 | 삼성전자주식회사 | 음성향상장치 및 방법 |
US7596231B2 (en) * | 2005-05-23 | 2009-09-29 | Hewlett-Packard Development Company, L.P. | Reducing noise in an audio signal |
JP4520970B2 (ja) * | 2005-11-21 | 2010-08-11 | 日本電信電話株式会社 | 類似度計算装置、エコー消去装置、類似度計算方法、エコー消去方法 |
JP4863713B2 (ja) * | 2005-12-29 | 2012-01-25 | 富士通株式会社 | 雑音抑制装置、雑音抑制方法、及びコンピュータプログラム |
JP2007199017A (ja) * | 2006-01-30 | 2007-08-09 | Omron Corp | 異音明瞭化方法、異音明瞭化装置および異音検査装置 |
JP4965891B2 (ja) * | 2006-04-25 | 2012-07-04 | キヤノン株式会社 | 信号処理装置およびその方法 |
JP4836720B2 (ja) | 2006-09-07 | 2011-12-14 | 株式会社東芝 | ノイズサプレス装置 |
JP4753821B2 (ja) | 2006-09-25 | 2011-08-24 | 富士通株式会社 | 音信号補正方法、音信号補正装置及びコンピュータプログラム |
JP2008148179A (ja) * | 2006-12-13 | 2008-06-26 | Fujitsu Ltd | 音声信号処理装置および自動利得制御装置における雑音抑圧処理方法 |
PL2118889T3 (pl) | 2007-03-05 | 2013-03-29 | Ericsson Telefon Ab L M | Sposób i sterownik do wygładzania stacjonarnego szumu tła |
JP5156260B2 (ja) * | 2007-04-27 | 2013-03-06 | ニュアンス コミュニケーションズ,インコーポレイテッド | 雑音を除去して目的音を抽出する方法、前処理部、音声認識システムおよびプログラム |
JP2008309955A (ja) * | 2007-06-13 | 2008-12-25 | Toshiba Corp | ノイズサプレス装置 |
ATE501506T1 (de) * | 2007-09-12 | 2011-03-15 | Dolby Lab Licensing Corp | Spracherweiterung mit anpassung von geräuschpegelschätzungen |
WO2009119460A1 (ja) * | 2008-03-24 | 2009-10-01 | 日本ビクター株式会社 | オーディオ信号処理装置及びオーディオ信号処理方法 |
US9159335B2 (en) | 2008-10-10 | 2015-10-13 | Samsung Electronics Co., Ltd. | Apparatus and method for noise estimation, and noise reduction apparatus employing the same |
JP5245714B2 (ja) * | 2008-10-24 | 2013-07-24 | ヤマハ株式会社 | 雑音抑圧装置及び雑音抑圧方法 |
CN102150206B (zh) * | 2008-10-24 | 2013-06-05 | 三菱电机株式会社 | 噪音抑制装置以及声音解码装置 |
JP5131149B2 (ja) * | 2008-10-24 | 2013-01-30 | ヤマハ株式会社 | 雑音抑圧装置及び雑音抑圧方法 |
JP5252563B2 (ja) * | 2009-01-20 | 2013-07-31 | Necシステムテクノロジー株式会社 | 電子機器、および電子機器における機構部品の診断方法、ならびに診断プログラム |
CN101986386B (zh) * | 2009-07-29 | 2012-09-26 | 比亚迪股份有限公司 | 一种语音背景噪声的消除方法和装置 |
JP5754899B2 (ja) | 2009-10-07 | 2015-07-29 | ソニー株式会社 | 復号装置および方法、並びにプログラム |
US8666092B2 (en) * | 2010-03-30 | 2014-03-04 | Cambridge Silicon Radio Limited | Noise estimation |
JP5850216B2 (ja) | 2010-04-13 | 2016-02-03 | ソニー株式会社 | 信号処理装置および方法、符号化装置および方法、復号装置および方法、並びにプログラム |
JP6075743B2 (ja) | 2010-08-03 | 2017-02-08 | ソニー株式会社 | 信号処理装置および方法、並びにプログラム |
JP2012058358A (ja) * | 2010-09-07 | 2012-03-22 | Sony Corp | 雑音抑圧装置、雑音抑圧方法およびプログラム |
JP5707842B2 (ja) | 2010-10-15 | 2015-04-30 | ソニー株式会社 | 符号化装置および方法、復号装置および方法、並びにプログラム |
JP2013148724A (ja) * | 2012-01-19 | 2013-08-01 | Sony Corp | 雑音抑圧装置、雑音抑圧方法およびプログラム |
US20150287406A1 (en) * | 2012-03-23 | 2015-10-08 | Google Inc. | Estimating Speech in the Presence of Noise |
JP6140685B2 (ja) | 2012-03-23 | 2017-05-31 | パナソニック インテレクチュアル プロパティ コーポレーション オブ アメリカPanasonic Intellectual Property Corporation of America | 帯域パワー算出装置及び帯域パワー算出方法 |
JP2014123011A (ja) * | 2012-12-21 | 2014-07-03 | Sony Corp | 雑音検出装置および方法、並びに、プログラム |
JP6263890B2 (ja) * | 2013-07-25 | 2018-01-24 | 沖電気工業株式会社 | 音声信号処理装置及びプログラム |
CN105531762B (zh) | 2013-09-19 | 2019-10-01 | 索尼公司 | 编码装置和方法、解码装置和方法以及程序 |
CA3162763A1 (en) | 2013-12-27 | 2015-07-02 | Sony Corporation | Decoding apparatus and method, and program |
JP6379839B2 (ja) * | 2014-08-11 | 2018-08-29 | 沖電気工業株式会社 | 雑音抑圧装置、方法及びプログラム |
RU2573274C1 (ru) * | 2014-11-05 | 2016-01-20 | Общество с ограниченной ответственностью Научно-производственное предприятие "ЭКРА" | Способ удаления импульсного шума (выброса) из электрической величины |
JP6447357B2 (ja) | 2015-05-18 | 2019-01-09 | 株式会社Jvcケンウッド | オーディオ信号処理装置、オーディオ信号処理方法及びオーディオ信号処理プログラム |
JP6668995B2 (ja) | 2016-07-27 | 2020-03-18 | 富士通株式会社 | 雑音抑圧装置、雑音抑圧方法及び雑音抑圧用コンピュータプログラム |
CN111144347B (zh) * | 2019-12-30 | 2023-06-27 | 腾讯科技(深圳)有限公司 | 一种数据处理方法、装置、平台及存储介质 |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4811404A (en) * | 1987-10-01 | 1989-03-07 | Motorola, Inc. | Noise suppression system |
SG49709A1 (en) | 1993-02-12 | 1998-06-15 | British Telecomm | Noise reduction |
JPH07193548A (ja) * | 1993-12-25 | 1995-07-28 | Sony Corp | 雑音低減処理方法 |
JP3484757B2 (ja) * | 1994-05-13 | 2004-01-06 | ソニー株式会社 | 音声信号の雑音低減方法及び雑音区間検出方法 |
JP3484801B2 (ja) | 1995-02-17 | 2004-01-06 | ソニー株式会社 | 音声信号の雑音低減方法及び装置 |
JP3591068B2 (ja) * | 1995-06-30 | 2004-11-17 | ソニー株式会社 | 音声信号の雑音低減方法 |
JPH09212196A (ja) | 1996-01-31 | 1997-08-15 | Nippon Telegr & Teleph Corp <Ntt> | 雑音抑圧装置 |
FR2768547B1 (fr) * | 1997-09-18 | 1999-11-19 | Matra Communication | Procede de debruitage d'un signal de parole numerique |
JPH11162240A (ja) | 1997-11-28 | 1999-06-18 | Matsushita Electric Works Ltd | 照明器具 |
JP3454190B2 (ja) * | 1999-06-09 | 2003-10-06 | 三菱電機株式会社 | 雑音抑圧装置および方法 |
US6671667B1 (en) * | 2000-03-28 | 2003-12-30 | Tellabs Operations, Inc. | Speech presence measurement detection techniques |
-
1999
- 1999-11-10 JP JP31912699A patent/JP3454206B2/ja not_active Expired - Lifetime
-
2000
- 2000-06-21 US US09/599,367 patent/US7158932B1/en not_active Expired - Fee Related
- 2000-07-13 EP EP00115227A patent/EP1100077B1/en not_active Expired - Lifetime
- 2000-07-13 DE DE60040895T patent/DE60040895D1/de not_active Expired - Lifetime
- 2000-07-31 CN CNB001221361A patent/CN1192360C/zh not_active Expired - Fee Related
-
2001
- 2001-10-22 HK HK01107378A patent/HK1037052A1/xx not_active IP Right Cessation
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1308914C (zh) * | 2001-06-06 | 2007-04-04 | 三菱电机株式会社 | 噪声抑制装置和方法 |
CN1302462C (zh) * | 2003-12-03 | 2007-02-28 | 富士通株式会社 | 降噪装置和降噪方法 |
CN1892822B (zh) * | 2005-05-31 | 2010-06-09 | 日本电气株式会社 | 噪声抑制方法和设备 |
CN101237814B (zh) * | 2005-06-10 | 2012-11-14 | 纽诺莫尼可斯有限公司 | 用于频谱地调整数字音频信号的数字回放设备和方法及装置 |
WO2008067735A1 (fr) * | 2006-12-05 | 2008-06-12 | Huawei Technologies Co., Ltd. | Procédé et dispositif de classement pour un signal sonore |
CN102906813A (zh) * | 2010-05-24 | 2013-01-30 | 日本电气株式会社 | 信号处理方法、信息处理装置和信号处理程序 |
CN102270454A (zh) * | 2010-06-07 | 2011-12-07 | 宇达电脑(上海)有限公司 | 改善音频输出的方法及装置 |
CN103718241B (zh) * | 2011-11-02 | 2016-05-04 | 三菱电机株式会社 | 噪音抑制装置 |
CN103718241A (zh) * | 2011-11-02 | 2014-04-09 | 三菱电机株式会社 | 噪音抑制装置 |
CN105103230A (zh) * | 2013-04-11 | 2015-11-25 | 日本电气株式会社 | 信号处理装置、信号处理方法、信号处理程序 |
US10431243B2 (en) | 2013-04-11 | 2019-10-01 | Nec Corporation | Signal processing apparatus, signal processing method, signal processing program |
CN111312277A (zh) * | 2014-03-03 | 2020-06-19 | 三星电子株式会社 | 用于带宽扩展的高频解码的方法及设备 |
US11676614B2 (en) | 2014-03-03 | 2023-06-13 | Samsung Electronics Co., Ltd. | Method and apparatus for high frequency decoding for bandwidth extension |
CN111312277B (zh) * | 2014-03-03 | 2023-08-15 | 三星电子株式会社 | 用于带宽扩展的高频解码的方法及设备 |
US11688406B2 (en) | 2014-03-24 | 2023-06-27 | Samsung Electronics Co., Ltd. | High-band encoding method and device, and high-band decoding method and device |
CN105445549A (zh) * | 2015-12-11 | 2016-03-30 | 中国电子科技集团公司第四十一研究所 | 一种微弱电磁信号频谱测量方法 |
CN111147983A (zh) * | 2018-11-06 | 2020-05-12 | 展讯通信(上海)有限公司 | 扬声器控制方法及装置、可读存储介质 |
Also Published As
Publication number | Publication date |
---|---|
EP1100077A2 (en) | 2001-05-16 |
DE60040895D1 (de) | 2009-01-08 |
EP1100077A3 (en) | 2002-07-10 |
EP1100077B1 (en) | 2008-11-26 |
US7158932B1 (en) | 2007-01-02 |
CN1192360C (zh) | 2005-03-09 |
HK1037052A1 (en) | 2002-01-25 |
JP2001134287A (ja) | 2001-05-18 |
JP3454206B2 (ja) | 2003-10-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1192360C (zh) | 噪声抑制装置 | |
CN1282155C (zh) | 噪声抑制装置和方法 | |
CN1308914C (zh) | 噪声抑制装置和方法 | |
CN1146155C (zh) | 噪声抑制装置 | |
CN1164036C (zh) | 声音回波和噪声的消除 | |
CN1113335A (zh) | 降低语音信号中噪声的方法和检测噪声域的方法 | |
CN1193644C (zh) | 用频谱减少来降低双送话器信号噪声的***和方法 | |
CN1192358C (zh) | 声音信号加工方法和声音信号加工装置 | |
CN1905006A (zh) | 噪声抑制***与方法及程序 | |
CN1918461A (zh) | 在存在背景噪声时用于语音增强的方法和设备 | |
CN1302462C (zh) | 降噪装置和降噪方法 | |
CN1169788A (zh) | 谱削减噪声抑制方法 | |
CN1669074A (zh) | 话音增强装置 | |
CN1409577A (zh) | 台词分量强调装置 | |
CN1141548A (zh) | 用来降低语音信号中噪声的方法和装置 | |
CN1951147A (zh) | 啸声检测设备和方法 | |
CN1606687A (zh) | 音频解码设备和方法 | |
CN1159703C (zh) | 声音识别*** | |
CN1969320A (zh) | 噪声抑制装置及噪声抑制方法 | |
CN1223991C (zh) | 声频信号处理设备及方法 | |
CN101031963A (zh) | 处理有噪声的声音信号的方法以及实现该方法的装置 | |
CN101079266A (zh) | 基于多统计模型和最小均方误差实现背景噪声抑制的方法 | |
CN1624762A (zh) | 主动噪音消除头盔、利用该头盔的车辆***、以及头盔内噪音消除方法 | |
CN1310431C (zh) | 用于编码音频信号的设备和方法 | |
CN1493073A (zh) | 噪声抑制方法和设备 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C06 | Publication | ||
PB01 | Publication | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20050309 Termination date: 20170731 |