CN102820036A - 一种自适应消除噪声的方法和装置 - Google Patents
一种自适应消除噪声的方法和装置 Download PDFInfo
- Publication number
- CN102820036A CN102820036A CN2012103304758A CN201210330475A CN102820036A CN 102820036 A CN102820036 A CN 102820036A CN 2012103304758 A CN2012103304758 A CN 2012103304758A CN 201210330475 A CN201210330475 A CN 201210330475A CN 102820036 A CN102820036 A CN 102820036A
- Authority
- CN
- China
- Prior art keywords
- wave filter
- noise
- microphone
- signal
- 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
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000001914 filtration Methods 0.000 claims abstract description 66
- 238000012546 transfer Methods 0.000 claims description 50
- 230000009467 reduction Effects 0.000 claims description 34
- 238000004422 calculation algorithm Methods 0.000 claims description 32
- 238000013459 approach Methods 0.000 claims description 18
- 230000008676 import Effects 0.000 claims description 6
- 230000000694 effects Effects 0.000 abstract description 4
- 230000004044 response Effects 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 15
- 238000004458 analytical method Methods 0.000 description 8
- 230000008569 process Effects 0.000 description 7
- 238000012545 processing Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000021615 conjugation Effects 0.000 description 2
- 230000002950 deficient Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000010276 construction Methods 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
- 230000006872 improvement Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000017105 transposition 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
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10K—SOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
- G10K11/00—Methods or devices for transmitting, conducting or directing sound in general; Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
- G10K11/16—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
- G10K11/175—Methods or devices for protecting against, or for damping, noise or other acoustic waves in general using interference effects; Masking sound
-
- 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
-
- 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
- G10L21/0216—Noise filtering characterised by the method used for estimating noise
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
- H04R3/00—Circuits for transducers, loudspeakers or microphones
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
- H04R3/00—Circuits for transducers, loudspeakers or microphones
- H04R3/005—Circuits for transducers, loudspeakers or microphones for combining the signals of two or more microphones
-
- 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
- G10L21/0216—Noise filtering characterised by the method used for estimating noise
- G10L2021/02161—Number of inputs available containing the signal or the noise to be suppressed
- G10L2021/02165—Two microphones, one receiving mainly the noise signal and the other one mainly the speech signal
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04R—LOUDSPEAKERS, MICROPHONES, GRAMOPHONE PICK-UPS OR LIKE ACOUSTIC ELECTROMECHANICAL TRANSDUCERS; DEAF-AID SETS; PUBLIC ADDRESS SYSTEMS
- H04R2410/00—Microphones
- H04R2410/05—Noise reduction with a separate noise microphone
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Multimedia (AREA)
- Signal Processing (AREA)
- Health & Medical Sciences (AREA)
- Quality & Reliability (AREA)
- Computational Linguistics (AREA)
- Audiology, Speech & Language Pathology (AREA)
- Human Computer Interaction (AREA)
- General Health & Medical Sciences (AREA)
- Otolaryngology (AREA)
- Circuit For Audible Band Transducer (AREA)
- Soundproofing, Sound Blocking, And Sound Damping (AREA)
Abstract
本发明公开一种自适应消除噪声方法和装置,所述方法包括:使用第一滤波器对第一麦克风接收的信号进行滤波,使用第二滤波器对第二麦克风接收的信号进行滤波,将滤波后的信号相减得出降噪后信号;其中,在噪声段,使用降噪后信号分别对第一滤波器和第二滤波器的系数进行更新,使得第一滤波器滤波后的信号中包含的噪声分量和第二滤波器滤波后的信号中包含的噪声分量趋于相同;并且,在带噪语音段,分别保持第一滤波器和第二滤波器的系数不变,第一滤波器和第二滤波器分别使用上次噪声段更新后的系数对第一麦克风和第二麦克风接收的信号进行滤波。本发明能够解决现有技术由于使用一个FIR滤波器无法逼近消除噪声的最优解而导致消除噪声效果较差的问题。
Description
技术领域
本发明涉及信号处理领域,特别涉及一种自适应消除噪声的方法和装置。
背景技术
现有技术的LMS(Least Mean Square,最小均方)算法采用单滤波器结构,如图1所示。其原理如图2所示,对其中一路麦克风接收到的信号进行滤波,滤波后的信号与另一路麦克风接收到的信号进行相减得到降噪后的语音。该单滤波器结构的滤波器仅在噪声段进行更新,在带噪语音段,滤波器保持不变。
标准的时域LMS算法,对于卷积的非加性干扰噪声,其计算复杂度比较大,为了减小计算复杂度,Ferrara提出了FBLMS(Fast Block LMS,快速块最小均方)算法,该算法采用的是时频域结合的方法,即将原来在时域进行的卷积运算变换成频域的乘积运算,从而大大减小了计算的复杂度。
以下对现有技术中的单滤波器结构LMS算法存在的缺陷进行说明。
通过分析单滤波器结构中滤波器的理论最优解来阐述单滤波器结构存在的缺陷。由于在频域可以清晰地分析出滤波器的最优解,因此对滤波器理论最优解的分析计算在频域进行。
如图3所示,为单滤波器结构中滤波器频域最优解的分析示意图。图3中S1代表信号源,S2代表噪声源。由于FIR(Finite Impulse Response,有限脉冲响应)滤波器可以较为准确的表征信源到麦克风的传递函数,因此在分析中,采用FIR滤波器分别模拟信号源与第一麦克风间的信道传递函数H11,噪声源与第一麦克风间的信道传递函数H12,信号源与第二麦克风间的信道传递函数H21,噪声源与第二麦克风间的信道传递函数H22。第一麦克风接收到的信号是X1,第二麦克风接收到的信号是X2,W为滤波器,Y1是降噪后的信号。
可以获得如下公式:
X1=S1×H11+S2×H12 公式1
X2=S1×H21+S2×H22 公式2
Y1=X1-X2×W=(S1×H11+S2×H12)-(S1×H21+S2×H22)×W
公式3
=S1×(H11-H21×W)+S2×(H12-H22×W)
由于当W取最优解时,噪声源S2将被完全消除,从而可以推出W的最优解如公式4所示。
Y1=S1×(H11-H21×W)=S1×(H11-H21×H12/H22) 公式5
由公式5可知Y1是S1经过某种滤波后的形式,并不含S2的任何分量。
从上面得到的W=H12/H22最优解的形式可以得知,W的最优解并不是FIR滤波器,但在实际应用中为了保证滤波器的稳定且易实现性,通常采用FIR滤波器,这就会引入较大的误差,其原因是使用一个FIR滤波器不能很好的去逼近一个非FIR滤波器。
标准单滤波器结构LMS算法中滤波器的最优解是非FIR滤波器,而该结构中的滤波器在实际应用中通常采用FIR滤波器去逼近这个最优解,从而会引入较大的误差,造成噪声消除效果较差。
发明内容
本发明提供了一种自适应消除噪声的方法和装置,以解决现有技术由于使用一个FIR滤波器无法逼近消除噪声的最优解而导致消除噪声效果较差的问题。
本发明公开了一种自适应消除噪声方法,所述方法包括:
使用第一滤波器对第一麦克风接收的信号进行滤波,使用第二滤波器对第二麦克风接收的信号进行滤波,将滤波后的信号相减得出降噪后信号;
其中,在噪声段,使用降噪后信号分别对第一滤波器的系数和第二滤波器的系数进行更新,使得第一滤波器滤波后的信号中包含的噪声分量和第二滤波器滤波后的信号中包含的噪声分量趋于相同;并且,
在带噪语音段,分别保持第一滤波器的系数和第二滤波器的系数不变,第一滤波器使用上次噪声段更新后的系数对第一麦克风接收的信号进行滤波,第二滤波器使用上次噪声段更新后的系数对第二麦克风接收的信号进行滤波。
其中,所述在噪声段,使用降噪后信号分别对第一滤波器的系数和第二滤波器的系数进行更新,使得第一滤波器滤波后的信号中包含的噪声分量和第二滤波器滤波后的信号中包含的噪声分量趋于相同具体包括:
在噪声段,按如下方式更新第一滤波器的系数和第二滤波器的系数,第一滤波器的传递函数与第二滤波器的传递函数的比值,逼近噪声源和第二麦克风间的信道传递函数与噪声源和第一麦克风间的信道传递函数的比值。
其中,所述第一滤波器的传递函数与第二滤波器的传递函数的比值,逼近噪声源和第二麦克风间的信道传递函数与噪声源和第一麦克风间的信道传递函数的比值具体包括:
第一滤波器的传递函数逼近噪声源与第二麦克风间的信道传递函数,第二滤波器的传递函数逼近噪声源与第一麦克风间的信道传递函数。
其中,所述第一滤波器的传递函数与第二滤波器的传递函数的比值,逼近噪声源和第二麦克风间的信道传递函数与噪声源和第一麦克风间的信道传递函数的比值具体包括:
第一滤波器的传递函数逼近噪声源与第二麦克风间的信道传递函数与常量的乘积,第二滤波器的传递函数逼近噪声源与第一麦克风间的信道传递函数与所述常量的乘积。
其中,所述使用降噪后信号分别对第一滤波器的系数和第二滤波器的系数进行更新具体包括:
通过最小均方算法或快速块最小均方算法,使用降噪后信号分别更新第一滤波器的系数和第二滤波器的系数。
本发明还公开了一种自适应消除噪声装置,所述装置包括:第一麦克风、第二麦克风、第一滤波器、第二滤波器、和减法器,
第一麦克风将接收信号输入第一滤波器,第一滤波器将滤波后信号输入减法器;
第二麦克风将接收信号输入第二滤波器,第二滤波器将滤波后信号输入减法器;
减法器将第一滤波器和第二滤波器滤波后的信号相减得出降噪后信号;
其中,在噪声段,第一滤波器的系数和第二滤波器系数分别根据降噪后信号进行更新,使得第一滤波器滤波后的信号中包含的噪声分量和第二滤波器滤波后的信号中包含的噪声分量趋于相同;并且,
在带噪语音段,第一滤波器的系数和第二滤波器的系数分别保持不变,第一滤波器对第一麦克风接收的信号进行滤波时使用的系数为上次噪声段更新后的系数,第二滤波器对第二麦克风接收的信号进行滤波时使用的系数为上次噪声段更新后的系数。
其中,第一滤波器的传递函数与第二滤波器的传递函数的比值,逼近噪声源和第二麦克风间的信道传递函数与噪声源和第一麦克风间的信道传递函数的比值。
其中,第一滤波器的传递函数逼近噪声源与第二麦克风间的信道传递函数,第二滤波器的传递函数逼近噪声源与第一麦克风间的信道传递函数。
其中,第一滤波器的传递函数逼近噪声源与第二麦克风间的信道传递函数与常量的乘积,第二滤波器的传递函数逼近噪声源与第一麦克风间的信道传递函数与所述常量的乘积。
其中,第一滤波器的系数具体通过最小均方算法或快速块最小均方算法,根据降噪后信号进行更新;
第二滤波器的系数具体通过最小均方算法或快速块最小均方算法,根据降噪后信号进行更新。
本发明的有益效果是:在噪声段,使用降噪后信号分别对第一滤波器和第二滤波器的系数进行更新,使得第一滤波器滤波后的信号中包含的噪声分量和第二滤波器滤波后的信号中包含的噪声分量趋于相同;并且,在带噪语音段,分别保持第一滤波器的系数和第二滤波器的系数不变,第一滤波器和第二滤波器分别使用上次噪声段更新后的系数对第一麦克风和第二麦克风接收的信号进行滤波;进而在将两个滤波器滤波后的信号相减时,使信号中的噪音分量基本被相互抵消掉,从而消除噪声效果被增强。
附图说明
图1为现有技术的LMS采用单滤波器消除噪声方法的示意图。
图2为现有技术的LMS采用单滤波器消除噪声方法的原理图。
图3为现有技术的LMS采用单滤波器消除噪声的频域最优解的原理分析示意图。
图4为本发明实施例的自适应消除噪声方法的流程图。
图5为本发明实施例的自适应消除噪声方法的原理图。
图6为本发明实施例的自适应消除噪声方法的原理分析示意图。
图7为本发明实施例的自适应消除噪声方法的时域处理流程图。
图8为本发明实施例的自适应消除噪声方法的示意图。
图9为本发明实施例的自适应消除噪声方法的频域处理流程图。
图10为本发明实施例的自适应消除噪声装置的结构图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明实施方式作进一步地详细描述。
实施例一
参见图4,为本发明实施例的自适应消除噪声方法的流程图。所述方法包括如下步骤。
步骤S100,第一麦克风接收信号,第二麦克风接收信号。
步骤S200,在噪声段,使用降噪后信号分别对第一滤波器的系数和第二滤波器的系数进行更新,使得第一滤波器滤波后的信号中包含的噪声分量和第二滤波器滤波后的信号中包含的噪声分量趋于相同,并使用第一滤波器对第一麦克风接收的信号进行滤波,使用第二滤波器对第二麦克风接收的信号进行滤波,将滤波后的信号相减得出降噪后信号。
步骤S300,在带噪语音段,分别保持第一滤波器的系数和第二滤波器的系数不变,第一滤波器使用上次噪声段更新后的系数对第一麦克风接收的信号进行滤波,第二滤波器使用上次噪声段更新后的系数对第二麦克风接收的信号进行滤波。
实施例二
在实施例二中,滤波器的更新过程如下所述。
所述在噪声段,对第一滤波器的系数和第二滤波器的系数进行更新具体包括:在噪声段,按如下方式更新第一滤波器的系数和第二滤波器的系数,
第一滤波器的传递函数与第二滤波器的传递函数的比值,逼近噪声源和第二麦克风间的信道传递函数与噪声源和第一麦克风间的信道传递函数的比值。
下面对本实施例中,自适应消除噪声的方法的原理进行说明。图5为本发明实施例的自适应消除噪声方法的原理图。图6为本发明实施例的自适应消除噪声方法的原理分析示意图。
参照图6,S1代表信号源,S2代表噪声源,由X1是第一麦克风上接收到的信号的频域值,X2第二麦克风上接收到的信号的频域值,W1,W2分别为第一滤波器和第二滤波器的传递函数,Y1是降噪后的信号的频域值。
可以获得如下公式。
X1=S1×H11+S2×H12 公式6
X2=S1×H21+S2×H22 公式7
Y1=X1×W1-X2×W2=(S1×H11+S2×H12)×W1-(S1×H21+S2×H22)×W2
=S1×(H11×W1-H21×W2)+S2×(H12×W1-H22×W2) 公式8
由于当W取最优解时,噪声源S2将被完全消除,因此存在公式9所示的两个滤波器W1和W2间关系。
当两个滤波器的传递函数关系满足公式9时,降噪后的信号为
Y1是S1经过某种滤波后的形式,由上述分析可知Y1不含S2的任何分量。
在本实施例中,可以通过多种方式使得第一滤波器与第二滤波器的传递函数比值逼近噪声源与第二麦克风和噪声源与第一麦克风间的信道传递函数的比值。
例如,第一滤波器的传递函数逼近噪声源与第二麦克风间的信道传递函数,第二滤波器的传递函数逼近噪声源与第一麦克风间的信道传递函数。
如图6所示,为该举例中自适应消除噪声的方法的原理分析示意图。
第一滤波器的传递函数为W1,W1=H22,第二滤波器的传递函数为W2,W2=H12。此时,两个滤波器滤波后信号中的噪声分量相同。因此,在本举例中,使W1逼近H22,使W2逼近H12,能够保证第一滤波器和第二滤波器滤波后的信号中包含的噪声分量尽可能相同,从而有效消除噪声。
又例如,第一滤波器的传递函数逼近噪声源与第二麦克风间的信道传递函数与常量的乘积,第二滤波器的传递函数逼近噪声源与第一麦克风间的信道传递函数与所述常量的乘积。该常量可以为常数或某个传递函数。即W1=H22·H,W2=H12·H,H为某一传递函数或常数。
在该举例中,同样保证第一滤波器和第二滤波器滤波后的信号中包含的噪声分量尽可能相同,从而有效消除噪声。
其中,使用最小均方算法或快速块最小均方算法更新滤波器(第一滤波器或第二滤波器)系数,使得滤波器逼近对应的传递函数。
由于当两个滤波器的传递函数关系满足公式9时,能够将信号中噪声进行消除,因此使用两个FIR滤波器使其相互关系逼近公式9,这样引入的误差就明显减少,从而大大提升了降噪效果。
在该方式中,每次使用上次噪声段最新更新的滤波器系数进行滤波,两个滤波器滤波后的信号中噪声分量趋于相同,两者相互抵消,使得降噪后信号中噪声分量不断减少,输出语音质量不断被提高。
实施例三
本实施例中采用时域LMS算法对滤波器中系数进行更新,本实施例的自适应消除噪声方法的时域处理流程如图7所示,本实施例自适应消除噪声方法的示意图如图8所示,其中采用双滤波器消除噪声。
步骤S701,第一麦克风和第二麦克风接收信号。
步骤S702,判断信号是否为噪声段,如果是,执行步骤S703,否则,执行步骤S704。
如果信号为带噪语音段的信号,则不进行滤波器系数更新,滤波器使用上一次在噪声段更新过的系数。
步骤S703,更新第一和第二滤波器系数。
步骤S704,使用滤波器将信号在时域进行滤波。
步骤S705,将两个滤波器滤波后信号相减,输出消除噪声后的信号。
下面根据图8所示的示意图具体说明步骤S703中更新第一和第二滤波器系数的过程。
采用时域LMS算法更新双滤波器结构中滤波器系数。第一滤波器滤波后信号为y(n),如公式11所示,其为输入信号经第一滤波器后的带噪信号。第二滤波器滤波后信号为d(n),如公式12所示,其为输入信号经第二滤波器后的带噪信号。两个滤波器信号相减后输出信号为e(n),如公式13所示。
e(n)=d(n)-y(n) 公式13
采用LMS算法对滤波器的传递函数进行更新,第一滤波器的传递函数按公式14进行更新,第二滤波器的传递函数按公式15进行更新。
其中W1(n),W2(n),X1(n),X2(n)均表示列向量,上标T表示转置,且
X1(n)=[x1(n)x1(n-1)……x1(n-N+1)]T
X2(n)=[x2(n)x2(n-1)……x2(n-N+1)]T
其中e(n)表示降噪后的信号,d(n)表示第一滤波器滤波后信号,y(n)表示第二滤波器滤波后信号,W1(n)表示第一滤波器的传递函数,W2(n)表示第二滤波器的传递函数,μ表示步长因子,X1(n)表示第一麦克风接收到的信号向量,X2(n)表示第二麦克风接收到的信号向量,N表示滤波器的阶数。
实施例四
本实施例中采用时频域结合的FBLMS算法对滤波器中系数进行更新,本实施例中自适应消除噪声方法的频域处理流程如图9所示。
步骤S901,第一麦克风和第二麦克风接收信号。
步骤S902,将第一麦克风和第二麦克风接收信号分块,并变换到频域。
步骤S903,判断信号是否为噪声段,如果是,执行步骤S904,否则,执行步骤S905。
如果信号为带噪语音段的信号,则不进行滤波器参数更新,滤波器使用上一次在噪声段更新过的系数。
步骤S904,在频域更新第一和第二滤波器系数。
步骤S905,在频域滤波,将滤波后信号变换到时域。
步骤S906,将两个滤波器滤波后信号相减,输出消除噪声后的信号。
参考图5所示的原理图,具体说明步骤S904中更新第一和第二滤波器系数的过程。
下面给出采用双滤波器结构的FBLMS算法的滤波器更新公式,其中“*”表示卷积。
其中,第一滤波器滤波后信号为y(n),如公式16所示,其为输入信号经第一滤波器后的带噪信号。第二滤波器滤波后信号为d(n),如公式17所示,其为输入信号经第二滤波器后的带噪信号。两个滤波器信号相减后输出信号为e(n),如公式18所示。
y(n)=w1(n)*x1(n) 公式16
d(n)=w2(n)*x2(n) 公式17
e(n)=d(n)-y(n) 公式18
将公式18做FFT(Fast Fourier Transform,快速傅里叶变换)变换到频域如公式19所示
E(k)=D(k)-Y(k)=W2(k)·X2(k)-W1(k)·X1(k) 公式19
采用FBLMS算法的原理的如下公式。
其中e(n)表示降噪后的信号,E(k)为e(n)的频域表示,d(n)表示第一滤波器滤波后信号,D(k)是d(n)的频域表示,y(n)表示第二滤波器滤波后信号,Y(k)是y(n)的频域表示,X1(k)为第一麦克风接收到的信号的频域表示,X2(k)为第二麦克风接收的信号的频域表示,W1、W2表示自适应滤波器的传递函数的频域表示,μ表示步长因子,表示X1(k)的共轭,表示X2(k)的共轭。
基于公式22和公式23,采用FBLMS算法进行滤波器系数更新。
1、滤波
设长度为N的两个频域滤波器为wF1(k)、wF2(k),将第一麦克风和第二麦克风接收到的信号前后均填充N个零,而后进行分块得到长度为L+N-1的块信号块间重叠N个数据。
2、误差估计
e(m)=d(N:L+N-1)-y(N:L+N-1) 公式28
其中,m=1:L(表示1至L);d(N:L+N-1)是公式27中d(k)的最后L个元素,与图5中的d(n)对应;y(N:L+N-1)是公式26中y(k)的最后L个元素,与图5中的y(n)对应。e(m)是降噪后的信号。
3、滤波器更新
4、滤波器约束
公式32
公式30与公式31中的滤波器传递函数中包含多余错误数据,通过公式32与公式33,在传递函数中抛去多余错误数据后填充零。
参见图10,为本发明实施例的自适应消除噪声的装置的结构图。
所述装置包括:第一麦克风110、第二麦克风120、第一滤波器210、第二滤波器220、和减法器300。
第一麦克风110将接收信号输入第一滤波器210,第一滤波器210将滤波后信号输入减法器300;
第二麦克风120将接收信号输入第二滤波器220,第二滤波器220将滤波后信号输入减法器300;
减法器300将第一滤波器210和第二滤波器220滤波后的信号相减得出降噪后信号。
其中,在噪声段,第一滤波器210的系数和第二滤波器220的系数分别根据降噪后信号进行更新,使得第一滤波器210滤波后的信号中包含的噪声分量和第二滤波器220滤波后的信号中包含的噪声分量趋于相同;
并且,在带噪语音段,第一滤波器210的系数和第二滤波器220的系数分别保持不变,第一滤波器210对第一麦克风110接收的信号进行滤波时使用的系数为上次噪声段更新后的系数,第二滤波器220对第二麦克风120接收的信号进行滤波时使用的系数为上次噪声段更新后的系数。
进一步地,第一滤波器210的传递函数与第二滤波器220的传递函数的比值,逼近噪声源和第二麦克风120间的信道传递函数与噪声源和第一麦克风110间的信道传递函数的比值。
进一步地,第一滤波器210的传递函数逼近噪声源与第二麦克风120间的信道传递函数,第二滤波器220的传递函数逼近噪声源与第一麦克风110间的信道传递函数。
进一步地,第一滤波器210的传递函数逼近噪声源与第二麦克风120间的信道传递函数与常量的乘积,第二滤波器220的传递函数逼近噪声源与第一麦克风110间的信道传递函数与所述常量的乘积。
进一步地,第一滤波器210的系数具体通过最小均方算法或快速块最小均方算法,根据降噪后信号进行更新;
第二滤波器220的系数具体通过最小均方算法或快速块最小均方算法,根据降噪后信号进行更新。
以上所述仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内所作的任何修改、等同替换、改进等,均包含在本发明的保护范围内。
Claims (10)
1.一种自适应消除噪声方法,其特征在于,所述方法包括:
使用第一滤波器对第一麦克风接收的信号进行滤波,使用第二滤波器对第二麦克风接收的信号进行滤波,将滤波后的信号相减得出降噪后信号;
其中,在噪声段,使用降噪后信号分别对第一滤波器的系数和第二滤波器的系数进行更新,使得第一滤波器滤波后的信号中包含的噪声分量和第二滤波器滤波后的信号中包含的噪声分量趋于相同;并且,
在带噪语音段,分别保持第一滤波器的系数和第二滤波器的系数不变,第一滤波器使用上次噪声段更新后的系数对第一麦克风接收的信号进行滤波,第二滤波器使用上次噪声段更新后的系数对第二麦克风接收的信号进行滤波。
2.根据权利要求1所述的方法,其特征在于,
所述在噪声段,使用降噪后信号分别对第一滤波器的系数和第二滤波器的系数进行更新,使得第一滤波器滤波后的信号中包含的噪声分量和第二滤波器滤波后的信号中包含的噪声分量趋于相同具体包括:
在噪声段,按如下方式更新第一滤波器的系数和第二滤波器的系数,第一滤波器的传递函数与第二滤波器的传递函数的比值,逼近噪声源和第二麦克风间的信道传递函数与噪声源和第一麦克风间的信道传递函数的比值。
3.根据权利要求2所述的方法,其特征在于,
所述第一滤波器的传递函数与第二滤波器的传递函数的比值,逼近噪声源和第二麦克风间的信道传递函数与噪声源和第一麦克风间的信道传递函数的比值具体包括:
第一滤波器的传递函数逼近噪声源与第二麦克风间的信道传递函数,第二滤波器的传递函数逼近噪声源与第一麦克风间的信道传递函数。
4.根据权利要求2所述的方法,其特征在于,
所述第一滤波器的传递函数与第二滤波器的传递函数的比值,逼近噪声源和第二麦克风间的信道传递函数与噪声源和第一麦克风间的信道传递函数的比值具体包括:
第一滤波器的传递函数逼近噪声源与第二麦克风间的信道传递函数与常量的乘积,第二滤波器的传递函数逼近噪声源与第一麦克风间的信道传递函数与所述常量的乘积。
5.根据权利要求1所述的方法,其特征在于,
所述使用降噪后信号分别对第一滤波器的系数和第二滤波器的系数进行更新具体包括:
通过最小均方算法或快速块最小均方算法,使用降噪后信号分别更新第一滤波器的系数和第二滤波器的系数。
6.一种自适应消除噪声装置,其特征在于,所述装置包括:第一麦克风、第二麦克风、第一滤波器、第二滤波器、和减法器,
第一麦克风将接收信号输入第一滤波器,第一滤波器将滤波后信号输入减法器;
第二麦克风将接收信号输入第二滤波器,第二滤波器将滤波后信号输入减法器;
减法器将第一滤波器和第二滤波器滤波后的信号相减得出降噪后信号;
其中,在噪声段,第一滤波器的系数和第二滤波器系数分别根据降噪后信号进行更新,使得第一滤波器滤波后的信号中包含的噪声分量和第二滤波器滤波后的信号中包含的噪声分量趋于相同;并且,
在带噪语音段,第一滤波器的系数和第二滤波器的系数分别保持不变,第一滤波器对第一麦克风接收的信号进行滤波时使用的系数为上次噪声段更新后的系数,第二滤波器对第二麦克风接收的信号进行滤波时使用的系数为上次噪声段更新后的系数。
7.根据权利要求6所述的装置,其特征在于,
第一滤波器的传递函数与第二滤波器的传递函数的比值,逼近噪声源和第二麦克风间的信道传递函数与噪声源和第一麦克风间的信道传递函数的比值。
8.根据权利要求7所述的装置,其特征在于,
第一滤波器的传递函数逼近噪声源与第二麦克风间的信道传递函数,第二滤波器的传递函数逼近噪声源与第一麦克风间的信道传递函数。
9.根据权利要求7所述的装置,其特征在于,
第一滤波器的传递函数逼近噪声源与第二麦克风间的信道传递函数与常量的乘积,第二滤波器的传递函数逼近噪声源与第一麦克风间的信道传递函数与所述常量的乘积。
10.根据权利要求6所述的装置,其特征在于,
第一滤波器的系数具体通过最小均方算法或快速块最小均方算法,根据降噪后信号进行更新;
第二滤波器的系数具体通过最小均方算法或快速块最小均方算法,根据降噪后信号进行更新。
Priority Applications (7)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210330475.8A CN102820036B (zh) | 2012-09-07 | 2012-09-07 | 一种自适应消除噪声的方法和装置 |
PCT/CN2013/082791 WO2014036918A1 (zh) | 2012-09-07 | 2013-09-02 | 一种自适应消除噪声的方法和装置 |
EP13835187.9A EP2814030B1 (en) | 2012-09-07 | 2013-09-02 | Method and device for self-adaptive noise reduction |
KR1020147026030A KR101538282B1 (ko) | 2012-09-07 | 2013-09-02 | 자기 적응적으로 잡음을 제거하기 위한 방법 및 장치 |
US14/412,080 US9570062B2 (en) | 2012-09-07 | 2013-09-02 | Method and device for self-adaptively eliminating noises |
JP2014555935A JP5762650B2 (ja) | 2012-09-07 | 2013-09-02 | 適応型ノイズ除去方法及び装置 |
DK13835187.9T DK2814030T3 (en) | 2012-09-07 | 2013-09-02 | Method and device for self-adjusting noise reduction |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210330475.8A CN102820036B (zh) | 2012-09-07 | 2012-09-07 | 一种自适应消除噪声的方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102820036A true CN102820036A (zh) | 2012-12-12 |
CN102820036B CN102820036B (zh) | 2014-04-16 |
Family
ID=47304121
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210330475.8A Active CN102820036B (zh) | 2012-09-07 | 2012-09-07 | 一种自适应消除噪声的方法和装置 |
Country Status (7)
Country | Link |
---|---|
US (1) | US9570062B2 (zh) |
EP (1) | EP2814030B1 (zh) |
JP (1) | JP5762650B2 (zh) |
KR (1) | KR101538282B1 (zh) |
CN (1) | CN102820036B (zh) |
DK (1) | DK2814030T3 (zh) |
WO (1) | WO2014036918A1 (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014036918A1 (zh) * | 2012-09-07 | 2014-03-13 | 歌尔声学股份有限公司 | 一种自适应消除噪声的方法和装置 |
CN105139853A (zh) * | 2015-08-13 | 2015-12-09 | 深圳市双平泰科技有限公司 | 一种体征检测装置的控制方法及装置 |
CN106373586A (zh) * | 2015-07-24 | 2017-02-01 | 南宁富桂精密工业有限公司 | 噪声滤除电路 |
CN110970052A (zh) * | 2019-12-31 | 2020-04-07 | 歌尔科技有限公司 | 降噪方法、装置、头戴显示设备和可读存储介质 |
CN112151002A (zh) * | 2019-06-27 | 2020-12-29 | 株洲中车时代电气股份有限公司 | 一种降噪座椅和降噪方法 |
CN112397080A (zh) * | 2020-10-30 | 2021-02-23 | 浙江大华技术股份有限公司 | 回声消除方法及装置、语音设备及计算机可读存储介质 |
CN115188358A (zh) * | 2022-07-21 | 2022-10-14 | 广东浦尔顿科技有限公司 | 一种直流充电模块的噪音处理方法 |
CN115841820A (zh) * | 2023-02-23 | 2023-03-24 | 中国电子科技集团公司第十研究所 | 一种机载话音处理平台 |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104581534B (zh) * | 2015-01-22 | 2017-12-12 | 四川省阳子森环保设备有限公司 | 一种高效噪音控制装置 |
US10586552B2 (en) * | 2016-02-25 | 2020-03-10 | Dolby Laboratories Licensing Corporation | Capture and extraction of own voice signal |
CN108010543B (zh) * | 2018-01-04 | 2024-02-23 | 澜邦科技(成都)集团有限公司 | 一种电压跟随低失真音频播放系统 |
JP6635396B1 (ja) * | 2019-04-08 | 2020-01-22 | パナソニックIpマネジメント株式会社 | 音響ノイズ抑圧装置及び音響ノイズ抑圧方法 |
US11398241B1 (en) * | 2021-03-31 | 2022-07-26 | Amazon Technologies, Inc. | Microphone noise suppression with beamforming |
CN114007167B (zh) * | 2021-10-29 | 2023-08-25 | 中电科航空电子有限公司 | 一种模拟音频双向通信系统及通信方法 |
US11741934B1 (en) | 2021-11-29 | 2023-08-29 | Amazon Technologies, Inc. | Reference free acoustic echo cancellation |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040071284A1 (en) * | 2002-08-16 | 2004-04-15 | Abutalebi Hamid Reza | Method and system for processing subband signals using adaptive filters |
CN101192411A (zh) * | 2007-12-27 | 2008-06-04 | 北京中星微电子有限公司 | 大距离麦克风阵列噪声消除的方法和噪声消除系统 |
US7660425B1 (en) * | 1999-05-25 | 2010-02-09 | British Telecommunications Plc | Acoustic echo cancellation |
CN102543060A (zh) * | 2011-12-27 | 2012-07-04 | 瑞声声学科技(深圳)有限公司 | 有源噪声控制系统及其设计方法 |
CN202838949U (zh) * | 2012-09-07 | 2013-03-27 | 歌尔声学股份有限公司 | 一种自适应消除噪声装置 |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2950260B2 (ja) * | 1996-11-22 | 1999-09-20 | 日本電気株式会社 | 雑音抑圧送話装置 |
JP2000312395A (ja) * | 1999-04-28 | 2000-11-07 | Alpine Electronics Inc | マイクロホンシステム |
US8503691B2 (en) * | 2007-06-13 | 2013-08-06 | Aliphcom | Virtual microphone arrays using dual omnidirectional microphone array (DOMA) |
US7206418B2 (en) * | 2001-02-12 | 2007-04-17 | Fortemedia, Inc. | Noise suppression for a wireless communication device |
DE10118653C2 (de) * | 2001-04-14 | 2003-03-27 | Daimler Chrysler Ag | Verfahren zur Geräuschreduktion |
KR100480789B1 (ko) * | 2003-01-17 | 2005-04-06 | 삼성전자주식회사 | 피드백 구조를 이용한 적응적 빔 형성방법 및 장치 |
US20080152157A1 (en) * | 2006-12-21 | 2008-06-26 | Vimicro Corporation | Method and system for eliminating noises in voice signals |
EP2196988B1 (en) * | 2008-12-12 | 2012-09-05 | Nuance Communications, Inc. | Determination of the coherence of audio signals |
CN103229237B (zh) * | 2010-10-12 | 2016-05-18 | 日本电气株式会社 | 信号处理设备、信号处理方法 |
JP5561195B2 (ja) | 2011-02-07 | 2014-07-30 | 株式会社Jvcケンウッド | ノイズ除去装置およびノイズ除去方法 |
CN102820036B (zh) | 2012-09-07 | 2014-04-16 | 歌尔声学股份有限公司 | 一种自适应消除噪声的方法和装置 |
-
2012
- 2012-09-07 CN CN201210330475.8A patent/CN102820036B/zh active Active
-
2013
- 2013-09-02 US US14/412,080 patent/US9570062B2/en active Active
- 2013-09-02 KR KR1020147026030A patent/KR101538282B1/ko active IP Right Grant
- 2013-09-02 EP EP13835187.9A patent/EP2814030B1/en active Active
- 2013-09-02 WO PCT/CN2013/082791 patent/WO2014036918A1/zh active Application Filing
- 2013-09-02 DK DK13835187.9T patent/DK2814030T3/en active
- 2013-09-02 JP JP2014555935A patent/JP5762650B2/ja active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7660425B1 (en) * | 1999-05-25 | 2010-02-09 | British Telecommunications Plc | Acoustic echo cancellation |
US20040071284A1 (en) * | 2002-08-16 | 2004-04-15 | Abutalebi Hamid Reza | Method and system for processing subband signals using adaptive filters |
CN101192411A (zh) * | 2007-12-27 | 2008-06-04 | 北京中星微电子有限公司 | 大距离麦克风阵列噪声消除的方法和噪声消除系统 |
CN102543060A (zh) * | 2011-12-27 | 2012-07-04 | 瑞声声学科技(深圳)有限公司 | 有源噪声控制系统及其设计方法 |
CN202838949U (zh) * | 2012-09-07 | 2013-03-27 | 歌尔声学股份有限公司 | 一种自适应消除噪声装置 |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014036918A1 (zh) * | 2012-09-07 | 2014-03-13 | 歌尔声学股份有限公司 | 一种自适应消除噪声的方法和装置 |
US9570062B2 (en) | 2012-09-07 | 2017-02-14 | Goertek Inc | Method and device for self-adaptively eliminating noises |
CN106373586A (zh) * | 2015-07-24 | 2017-02-01 | 南宁富桂精密工业有限公司 | 噪声滤除电路 |
TWI622284B (zh) * | 2015-07-24 | 2018-04-21 | 新加坡商雲網科技新加坡有限公司 | 噪音濾除電路 |
CN106373586B (zh) * | 2015-07-24 | 2020-03-17 | 南宁富桂精密工业有限公司 | 噪声滤除电路 |
CN105139853A (zh) * | 2015-08-13 | 2015-12-09 | 深圳市双平泰科技有限公司 | 一种体征检测装置的控制方法及装置 |
CN112151002A (zh) * | 2019-06-27 | 2020-12-29 | 株洲中车时代电气股份有限公司 | 一种降噪座椅和降噪方法 |
CN112151002B (zh) * | 2019-06-27 | 2024-02-23 | 株洲中车时代电气股份有限公司 | 一种降噪座椅和降噪方法 |
CN110970052A (zh) * | 2019-12-31 | 2020-04-07 | 歌尔科技有限公司 | 降噪方法、装置、头戴显示设备和可读存储介质 |
CN112397080A (zh) * | 2020-10-30 | 2021-02-23 | 浙江大华技术股份有限公司 | 回声消除方法及装置、语音设备及计算机可读存储介质 |
CN115188358A (zh) * | 2022-07-21 | 2022-10-14 | 广东浦尔顿科技有限公司 | 一种直流充电模块的噪音处理方法 |
CN115841820A (zh) * | 2023-02-23 | 2023-03-24 | 中国电子科技集团公司第十研究所 | 一种机载话音处理平台 |
Also Published As
Publication number | Publication date |
---|---|
JP5762650B2 (ja) | 2015-08-12 |
KR20140127331A (ko) | 2014-11-03 |
US9570062B2 (en) | 2017-02-14 |
EP2814030A1 (en) | 2014-12-17 |
DK2814030T3 (en) | 2017-02-13 |
EP2814030B1 (en) | 2016-11-09 |
CN102820036B (zh) | 2014-04-16 |
WO2014036918A1 (zh) | 2014-03-13 |
JP2015511330A (ja) | 2015-04-16 |
US20150179160A1 (en) | 2015-06-25 |
KR101538282B1 (ko) | 2015-07-20 |
EP2814030A4 (en) | 2015-09-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102820036B (zh) | 一种自适应消除噪声的方法和装置 | |
CN102903368B (zh) | 用于卷积盲源分离的方法和设备 | |
Kuo et al. | Frequency-domain delayless active sound quality control algorithm | |
EP3175452B1 (en) | Method and apparatus for processing an audio signal, audio decoder, and audio encoder | |
CN101083640A (zh) | 低复杂度的降噪方法 | |
US8868416B2 (en) | Apparatus and method for cancelling echo in joint time domain and frequency domain | |
CN104954595A (zh) | 残留回声消除方法和装置 | |
CN202838949U (zh) | 一种自适应消除噪声装置 | |
CN106601266A (zh) | 回声消除方法、装置及系统 | |
CN102377454B (zh) | 一种回声消除的方法和装置 | |
US20120099637A1 (en) | Systems and methods for performing dynamic channel estimation | |
CN105324762A (zh) | 滤波器系数群计算设备以及滤波器系数群计算方法 | |
Wang et al. | Hybrid time–frequency algorithm for active sound quality control of vehicle interior noise based on stationary discrete wavelet transform | |
CN113889068A (zh) | 一种主动噪声消除方法、装置、电子设备和存储介质 | |
EP2730026B1 (en) | Low-delay filtering | |
EP2811758B1 (en) | Audio signal mixing | |
CN115801145A (zh) | 一种混合信号的时延估计方法、装置和电子设备 | |
KR100848789B1 (ko) | 크로스토크를 제거하기 위한 후처리 방법 | |
CN111402905A (zh) | 音频数据恢复方法、装置及蓝牙设备 | |
Duan et al. | A computational-efficient active sound tuning system for steady-state and transient vehicle powertrain response | |
JP4948019B2 (ja) | 適応信号処理装置およびその適応信号処理方法 | |
CN103270772A (zh) | 信号处理设备、信号处理方法和信号处理程序 | |
CN105957534A (zh) | 自适应滤波方法和自适应滤波器 | |
US20150358042A1 (en) | Nonlinear Filtering Using Polyphase Filter Banks | |
EP3770907B1 (en) | Systems and methods for estimating noise |
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 | ||
C56 | Change in the name or address of the patentee | ||
CP01 | Change in the name or title of a patent holder |
Address after: 261031 Dongfang Road, Weifang high tech Industrial Development Zone, Shandong, China, No. 268 Patentee after: Goertek Inc. Address before: 261031 Dongfang Road, Weifang high tech Industrial Development Zone, Shandong, China, No. 268 Patentee before: Goertek Inc. |