The present invention relates generally to communication, and more specifically to techniques for suppressing noise and interference in communication and voice recognition systems using an array microphone.
Communication and voice recognition systems are commonly used for many applications, such as hands-free car kit, cellular phone, hands-free voice control devices, telematics, teleconferencing system, and so on. These systems may be operated in noisy environments, such as in a vehicle or a restaurant. For each of these systems, one or multiple microphones in the system pick up the desired voice signal as well as noise and interference. The noise typically refers to local ambient noise. The interference may be from acoustic echo, reverberation, unwanted voice, and other artifacts.
Noise suppression is often required in many communication and voice recognition systems to suppress ambient noise and remove unwanted interference. For a communication or voice recognition system operating in a noisy environment, the microphone(s) in the system pick up the desired voice as well as noise. The noise is more severe for a hands-free system whereby the loudspeaker and microphone may be located some distance away from a talking user. The noise degrades communication quality and speech recognition rate if it is not dealt with in an appropriate manner.
For a system with a single microphone, noise suppression is conventionally achieved using a spectral subtract technique. For this technique, which performs signal processing in the frequency domain, the noise power spectrum of a noisy voice signal is estimated and subtracted from the power spectrum of the noisy voice signal to obtain an enhanced voice signal. The phase of the enhanced voice signal is set equal to the phase of the noisy voice signal. This technique is somewhat effective for stationary noise or slow-varying non-stationary (such as air-conditioner noise or fan noise, which does not change over time) but may not be effective for fast-varying non-stationary noise. Moreover, even for stationary noise, this technique can cause voice distortion if the noisy voice signal has a low signal-to-noise ratio (SNR). Conventional noise suppression for stationary noise is described in various literatures including U.S. Pat. Nos. 4,185,168 and 5,768,473.
For a system with multiple microphones, an array microphone is formed by placing these microphones at different positions sufficiently far apart. The array microphone forms a signal beam that is used to suppress noise and interference outside of the beam. Conventionally, the spacing between the microphones needs to be greater than a certain minimum distance D in order to form the desired beam. This spacing requirement prevents the array microphone from being used in many applications where space is limited. Moreover, conventional beam-forming with the array microphone is typically not effective at suppressing noise in an environment with diffused noise. Conventional systems with array microphone are described in various literatures including U.S. Pat. Nos. 5,371,789, 5,383,164, 5,465,302 and 6,002,776.
As can be seen, techniques that can effectively suppress noise and interference in communication and voice recognition systems are highly desirable.
Techniques are provided herein to suppress both stationary and non-stationary noise and interference using an array microphone and a combination of time-domain and frequency-domain signal processing. These techniques are also effective at suppressing diffuse noise, which cannot be handled by a single microphone system and a conventional array microphone system. The inventive techniques can provide good noise and interference suppression, high voice quality, and faster voice recognition rate, all of which are highly desirable for hands-free full-duplex applications in communication or voice recognition systems.
The array microphone is composed of a combination of omni-directional microphones and uni-directional microphones. The microphones may be placed close to each other (i.e., closer than the minimum distance required by a conventional array microphone). This allows the array microphone to be used in various applications. The array microphone forms a signal beam at a desired direction. This beam is then used to suppress stationary and non-stationary noise and interference.
A specific embodiment of the invention provides a noise suppression system that includes an array microphone, at least one voice activity detector (VAD), a reference generator, a beam-former, and a multi-channel noise suppressor. The array microphone is composed of multiple microphones, which include at least one omni-directional microphone and at least one uni-directional microphone. Each microphone provides a respective received signal. One of the received signals is designated as the main signal, and the remaining received signal(s) are designated as secondary signal(s). The VAD(s) provide at least one voice detection signal, which is used to control the operation of the reference generator, the beam-former, and the multi-channel noise suppressor. The reference generator provides a reference signal based on the main signal, a first set of at least one secondary signal, and an intermediate signal from the beam-former. The beam-former provides the intermediate signal and a beam-formed signal based on the main signal, a second set of at least one secondary signal, and the reference signal. Depending on the number of microphones used for the array microphone, the first and second sets may include the same or different secondary signals. The reference signal has the desired voice signal suppressed, and the beam-formed signal has the noise and interference suppressed. The multi-channel noise suppressor further suppresses noise and interference in the beam-formed signal to provide an output signal having much of the noise and interference suppressed.
In one embodiment, the array microphone is composed of three microphones—one omni-directional microphone and two uni-directional microphones (which may be placed close to each other). The omni-directional microphone is referred to as the main microphone/channel and its received signal is the main signal a(n). One of the uni-directional microphones faces toward a desired talker and is referred to as a first secondary microphone/channel. Its received signal is the first secondary signal s1(n). The other uni-directional microphone faces away from the desired talker and is referred to as a second secondary microphone/channel. Its received signal is the second secondary signal s2(n).
In another embodiment, the array microphone is composed of two microphones—one omni-directional microphone and one uni-directional microphone (which again may be placed close to each other). The uni-directional microphone faces toward the desired talker and its received signal is the main signal a(n). The omni-directional microphone is the secondary microphone/channel and its received signal is the secondary signal s(n).
Various other aspects, embodiments, and features of the invention are also provided, as described in further detail below.
For clarity, various signals and controls described herein are labeled with lower case and upper case symbols. Time-variant signals and controls are labeled with “(n)” and “(m)”, where n denotes sample time and m denotes frame index. A frame is composed of L samples. Frequency-variant signals and controls are labeled with “(k,m)”, where k denotes frequency bin. Lower case symbols (e.g., s(n) and d(m)) are used to denote time-domain signals, and upper case symbols (e.g., B(k,m)) are used to denote frequency-domain signals.
The N received signals, provided by N microphones 112a through 112n placed at different positions, carry information for the differences in the microphone positions. The N digitized signals s1(n) through SN(n) are provided to a beam-former 118 and used to form a signal beam. This beam is used to suppress noise and interference outside of the beam and to enhance the desired voice within the beam. Beam-former 118 may be a fixed beam-former (e.g., a delay-and-sum beam-former) or an adaptive beam-former (e.g., an adaptive sidelobe cancellation beam-former). These various types of beam-former are well known in the art. Conventional array microphone system 100 is associated with several limitations that curtail its use and/or effectiveness, including (1) requirement of a minimum distance of D for the spacing between microphones and (2) marginal effectiveness for diffused noise.
In the embodiment shown in
Microphones 212a, 212b, and 212c provide three received signals, which are amplified by amplifiers 214a, 214b, and 214c, respectively. An ADC 216a receives and digitizes the amplified signal from amplifier 214a and provides a first secondary signal s1(n). An ADC 216b receives and digitizes the amplified signal from amplifier 214b and provides a main signal a(n). An ADC 216c receives and digitizes the amplified signal from amplifier 214c and provides a second secondary signal s2(n).
A first voice activity detector (VAD1) 220 receives the main signal a(n) and the first secondary signal s1(n). VAD 1220 detects for the presence of near-end voice based on a metric of total power over noise power, as described below. VAD1220 provides a first voice detection signal d1(n), which indicates whether or not near-end voice is detected.
A second voice activity detector (VAD2) 230 receives the main signal a(n) and the second secondary signal s2(n). VAD2230 detects for the absence of near-end voice based on a metric of the cross-correlation between the main signal and the desired voice signal over the total power, as described below. VAD2230 provides a second voice detection signal d2(n), which also indicates whether or not near-end voice is absent.
A reference generator 240 receives the main signal a(n), the first secondary signal s1(n), the first voice detection signal d1(n), and a first beam-formed signal b1(n). Reference generator 240 updates its coefficients based on the first voice detection signal d1(n), detects for the desired voice signal in the first secondary signal s1(n) and the first beam-formed signal b2(n), cancels the desired voice signal from the main signal a(n), and provides two reference signals r1(n) and r2(n). The reference signals r1(n) and r2(n) both contain mostly noise and interference. However, the reference signal r2(n) is more accurate than r1(n) in order to estimate the presence of noise and interference.
A beam-former 250 receives the main signal a(n), the second secondary signal s2(n), the second reference signal r2(n), and the second voice detection signal d2(n). Beam-former 250 updates its coefficients based on the second voice detection signal d2(n), detects for the noise and interference in the second secondary signal s2(n) and the second reference signal r2(n), cancels the noise and interference from the main signal a(n), and provides the two beam-formed signals b1(n) and b2(n). The beam-formed signal b2(n) is more accurate than b1(n) to represent the desired signal.
A delay unit 242 delays the second reference signal r2(n) by a delay of Ta and provides a third reference signal r3(n), which is r3(n)=r2(n−Ta). The delay Ta synchronizes (i.e., time-aligns) the third reference signal r3(n) with the second beam-formed signal b2(n).
A third voice activity detector (VAD3) 260 receives the third reference signal r3(n) and the second beam-formed signal b2(n). VAD3260 detects for the presence of near-end voice based on a metric of desired voice power over noise power, as described below. VAD3260 provides a third voice detection signal d3(m) to dual-channel noise suppressor 280, which also indicates whether or not near-end voice is detected. The third voice detection signal d3(m) is a function of frame index m instead of sample index n.
A dual-channel FFT unit 270 receives the second beam-formed signal b2(n) and the third reference signal r3(n). FFT unit 270 transforms the signal b2(n) from the time domain to the frequency domain using an L-point FFT and provides a corresponding frequency-domain beam-formed signal B(k,m). FFT unit 270 also transforms the signal r3(n) from the time domain to the frequency domain using the L-point FFT and provides a corresponding frequency-domain reference signal R(k,m).
A dual-channel noise suppressor 280 receives the frequency-domain signals B(k,m) and R(k,m) and the third voice detection signal d3(m). Noise suppressor 280 further suppresses noise and interference in the signal B(k,m) and provides a frequency-domain output signal Bo(k,m) having much of the noise and interference suppressed.
An inverse FFT unit 290 receives the frequency-domain output signal Bo(k,m), transforms it from the frequency domain to the time domain using an L-point inverse FFT, and provides a corresponding time-domain output signal bo(n). The output signal bo(n) may be converted to an analog signal, amplified, filtered, and so on, and provided to a speaker.
Within VAD 220x, a subtraction unit 310 subtracts the first secondary signal s1(n) from the main signal a(n) and provides a first difference signal e1(n), which is e1(n)=a(n)−s1(n). The first difference signal e1(n) contains mostly noise and interference. High-pass filters 312 and 314 respectively receive the signals a(n) and e1(n), filter these signals with the same set of filter coefficients to remove low frequency components, and provide filtered signals ã1(n) and {tilde over (e)}1(n), respectively. Power calculation units 316 and 318 then respectively receive the filtered signals ã1(n) and {tilde over (e)}1(n), compute the powers of the filtered signals, and provide computed powers pa1(n) and pe1(n), respectively. Power calculation units 316 and 318 may further average the computed powers. In this case, the averaged computed powers may be expressed as:
pa1(n)=a1·pa1(n−1)+(1−a1)·ã1(n)·ã1(n), and Eq (1a)
pe1(n)=a1·pe1(n−1)+(1−a1)·{tilde over (e)}1(n)·{tilde over (e)}1(n), Eq(1b)
where α1 is a constant that determines the amount of averaging and is selected such that 1>α1>0. A large value for α1 corresponds to more averaging and smoothing. The term pa1(n) includes the total power from the desired voice signal as well as noise and interference. The term pe1(n) includes mostly noise and interference power.
A divider unit 320 then receives the averaged powers pa1(n) and pe1(n) and calculates a ratio h1(n) of these two powers. The ratio h1(n) may be expressed as:
The ratio h1(n) indicates the amount of total power relative to the noise power. A large value for h1(n) indicates that the total power is large relative to the noise power, which may be the case if near-end voice is present. A larger value for h1(n) corresponds to higher confidence that near-end voice is present.
A smoothing filter 322 receives and filters or smoothes the ratio h1(n) and provides a smoothed ratio hs1(n). The smoothing may be expressed as:
hs1(n)=αh1·hs1(n−1)+(1−αh1)·h1(n), Eq (3)
where αh1 is a constant that determines the amount of smoothing and is selected as 1>αh1>0.
A threshold calculation unit 324 receives the instantaneous ratio h1(n) and the smoothed ratio hs1(n) and determines a threshold q1(n). To obtain q1(n), an initial threshold q1′(n) is first computed as:
where β1 is a constant that is selected such that β1>0. In equation (4), if the instantaneous ratio h1(n) is greater than β1hs1(n), then the initial threshold q1′(n) is computed based on the instantaneous ratio h1(n) in the same manner as the smoothed ratio hs1(n). Otherwise, the initial threshold for the prior sample period is retained (i.e., q1′(n)=q1′(n−1)) and the initial threshold q1′(n) is not updated with h1(n). This prevents the threshold from being updated under abnormal condition for small values of h1(n).
The initial threshold q1′(n) is further constrained to be within a range of values defined by Qmax1 and Qmin1. The threshold q1(n) is then set equal to the constrained initial threshold q1′(n), which may be expressed as:
where Qmax1 and Qmin1 are constants selected such that Qmax1>Qmin1.
The threshold q1(n) is thus computed based on a running average of the ratio h1(n), where small values of h1(n) are excluded from the averaging. Moreover, the threshold q1(n) is further constrained to be within the range of values defined by Qmax1 and Qmin1. The threshold q1(n) is thus adaptively computed based on the operating environment.
A comparator 326 receives the ratio h1(n) and the threshold q1(n), compares the two quantities h1(n) and q1(n), and provides the first voice detection signal d1(n) based on the comparison results. The comparison may be expressed as:
The voice detection signal d1(n) is set to 1 to indicate that near-end voice is detected and set to 0 to indicate that near-end voice is not detected.
Within VAD 230x, a subtraction unit 410 subtracts the main signal a(n) from the second secondary signal s2(n) and provides a second difference signal e2(n), which is e2(n)=s2(n)−a(n). High-pass filters 412 and 414 respectively receive the signals a(n) and e2(n), filter these signals with the same set of filter coefficients to remove low frequency components, and provide filtered signals ã2(n) and {tilde over (e)}2(n), respectively. The filter coefficients used for high-pass filters 412 and 414 may be the same or different from the filter coefficients used for high-pass filters 312 and 314.
A power calculation unit 416 receives the filtered signal ã2(n), computes the power of this filtered signal, and provides the computed power pa2(n). A correlation calculation unit 418 receives the filtered signals ã2(n) and {tilde over (e)}2(n), computes their cross correlation, and provides the correlation pae(n). Units 416 and 418 may further average their computed results. In this case, the averaged computed power from unit 416 and the averaged correlation from unit 418 may be expressed as:
pa2(n)=α2·pa2(n−1)+(1−α2)·ã2(n)·ã2(n), and Eq (7a)
pae(n)=α2·pae(n−1)+(1−α2)·ã2(n)·{tilde over (e)}2(n), Eq (7b)
where α2 is a constant that is selected such that 1>α2>0. The constant α2 for VAD2230x may be the same or different from the constant α1 for VAD1220x. The term pa2(n) includes the total power for the desired voice signal as well as noise and interference. The term pae(n) includes the correlation between a(n) and e2(n), which is typically negative if near-end voice is present.
A divider unit 420 then receives pa2(n) and pae(n) and calculates a ratio h2(n) of these two quantities, as follows:
A smoothing filter 422 receives and filters the ratio h2(n) to provide a smoothed ratio hs2(n), which may be expressed as:
hs2(n)=αh2·hs2(n−1)+(1−αh2)·h2(n), Eq(9)
where αh2 is a constant that is selected such that 1>αh2>0. The constant αh2 for VAD2230x may be the same or different from the constant αh1 for VAD1220x.
A threshold calculation unit 424 receives the instantaneous ratio h2(n) and the smoothed ratio hs2(n) and determines a threshold q2(n). To obtain q2(n), an initial threshold q2′(n) is first computed as:
where β2 is a constant that is selected such that β2>0. The constant β2 for VAD2230x may be the same or different from the constant β1 for VAD 1220x. In equation (10), if the instantaneous ratio h2(n) is greater than β2hs2(n), then the initial threshold q2′(n) is computed based on the instantaneous ratio h2(n) in the same manner as the smoothed ratio hs2(n). Otherwise, the initial threshold for the prior sample period is retained.
The initial threshold q2′(n) is further constrained to be within a range of values defined by Qmax2 and Qmin2. The threshold q2(n) is then set equal to the constrained initial threshold q2′(n), which may be expressed as:
where Qmax2 and Qmin2 are constants selected such that Qmax2>Qmin2.
A comparator 426 receives the ratio h2(n) and the threshold q2(n), compares the two quantities h2(n) and q2(n), and provides the second voice detection signal d2(n) based on the comparison results. The comparison may be expressed as:
The voice detection signal d2(n) is set to 1 to indicate that near-end voice is absent and set to 0 to indicate that near-end voice is present.
Within reference generator 240x, a delay unit 512 receives and delays the main signal a(n) by a delay of T1 and provides a delayed signal a(n−T1). The delay T1 accounts for the processing delays of an adaptive filter 520. For linear FIR-type adaptive filter, T1 is set to equal to half the filter length. Adaptive filter 520 receives the delayed signal a(n−T1) at its xin input, the first secondary signal s1(n) at its xref input, and the first voice detection signal d1(n) at its control input. Adaptive filter 520 updates its coefficients only when the first voice detection signal d1(n) is 1. These coefficients are then used to isolate the desired voice component in the first secondary signal s1(n). Adaptive filter 520 then cancels the desired voice component from the delayed signal a(n−T1) and provides the first reference signal r1(n) at its xout output. The first reference signal r1(n) contains mostly noise and interference. An exemplary design for adaptive filter 520 is described below.
A delay unit 522 receives and delays the first reference signal r1(n) by a delay of T2 and provides a delayed signal r1(n−T2). The delay T2 accounts for the difference in the processing delays of adaptive filters 520 and 540 and the processing delay of an adaptive filter 530. Adaptive filter 530 receives the first beam-formed signal b1(n) at its xref input, the delayed signal r1(n−T2) at its xin input, and the first voice detection signal d1(n) at its control input. Adaptive filter 530 updates its coefficients only when the first voice detection signal d1(n) is 1. These coefficients are then used to isolate the desired voice component in the first beam-formed signal b1(n). Adaptive filter 530 then further cancels the desired voice component from the delayed signal r1(n−T2) and provides the second reference signal r2(n) at its xout output. The second reference signal r2(n) contains mostly noise and interference. The use of two adaptive filters 520 and 530 to generate the reference signals can provide improved performance.
Within beam-former 250x, a delay unit 532 receives and delays the main signal a(n) by a delay of T3 and provides a delayed signal a(n−T3). The delay T3 accounts for the processing delays of adaptive filter 540. For linear FIR-type adaptive filter, T3 is set to equal to half the filter length. Adaptive filter 540 receives the delayed signal a(n−T3) at its xin input, the second secondary signal s2(n) at its xref input, and the second voice detection signal d2(n) at its control input. Adaptive filter 540 updates its coefficients only when the second voice detection signal d2(n) is 1. These coefficients are then used to isolate the noise and interference component in the second secondary signal s2(n). Adaptive filter 540 then cancels the noise and interference component from the delayed signal a(n−T3) and provides the first beam-formed signal b1(n) at its xout output. The first beam-formed signal b1(n) contains mostly the desired voice signal.
A delay unit 542 receives and delays the first beam-formed signal b1(n) by a delay of T4 and provides a delayed signal b1(n−T4). The delay T4 accounts for the total processing delays of adaptive filters 530 and 550. Adaptive filter 550 receives the delayed signal b1(n−T4) at its xin input, the second reference signal r2(n) at its xref input, and the second voice detection signal d2(n) at its control input. Adaptive filter 550 updates its coefficients only when the second voice detection signal d2(n) is 1. These coefficients are then used to isolate the noise and interference component in the second reference signal r2(n). Adaptive filter 550 then cancels the noise and interference component from the delayed signal b1(n−T4) and provides the second beam-formed signal b2(n) at its xout output. The second beam-formed signal b2(n) contains mostly the desired voice signal.
Within VAD 260x, high-pass filters 612 and 614 respectively receive the second beam-formed signal b2(n) from beam-former 250 and the third reference signal r3(n) from delay unit 242, filter these signals with the same set of filter coefficients to remove low frequency components, and provide filtered signals {tilde over (b)}2(n) and {tilde over (r)}3(n), respectively. Power calculation units 616 and 618 then respectively receive the filtered signals {tilde over (b)}2(n) and {tilde over (r)}3(n), compute the powers of the filtered signals, and provide computed powers pb2(n) and pr3(n), respectively. Power calculation units 616 and 618 may further average the computed powers. In this case, the averaged computed powers may be expressed as:
pb2(n)=α3·pb2(n−1)+(1−α3)·{tilde over (b)}2(n)·{tilde over (b)}2(n), and Eq(13a)
pr3(n)=α3·pr3(n−1)+(1−α3)·{tilde over (r)}3(n)·{tilde over (r)}3(n), Eq(13b)
where α3 is a constant that is selected such that 1>α3>0. The constant α3 for VAD3260x may be the same or different from the constant α2 for VAD2230x and the constant α1 for VAD1220x.
A divider unit 620 then receives the averaged powers pb2(n) and pr3(n) and calculates a ratio h3(n) of these two powers, as follows:
The ratio h3(n) indicates the amount of desired voice power relative to the noise power.
A smoothing filter 622 receives and filters the ratio h3(n) to provide a smoothed ratio hs3(n), which may be expressed as:
hs3(n)=αh3·hs3(n−1)+(1−αh3)·h3(n), Eq (15)
where αh3 is a constant that is selected such that 1>αh3>0. The constant αh3 for VAD3260x may be the same or different from the constant αh2 for VAD2230x and the constant αh1 for VAD1220x.
A threshold calculation unit 624 receives the instantaneous ratio h3(n) and the smoothed ratio hs3(n) and determines a threshold q3(n). To obtain q3(n), an initial threshold q3′(n) is first computed as:
where β3 is a constant that is selected such that β3>0. In equation (16), if the instantaneous ratio h3(n) is greater than β3hs3(n), then the initial threshold q3′(n) is computed based on the instantaneous ratio h3(n) in the same manner as the smoothed ratio hs3(n). Otherwise, the initial threshold for the prior sample period is retained.
The initial threshold q3(n) is further constrained to be within a range of values defined by Qmax3 and Qmin3. The threshold q3(n) is then set equal to the constrained initial threshold q3′(n), which may be expressed as:
where Qmax3 and Qmin3 are constants selected such that Qmax3>Qmin3.
A comparator 626 receives the ratio h3(n) and the threshold q3(n) and averages these quantities over each frame m. For each frame, the ratio h3(m) is obtained by accumulating L values for h3(n) for that frame and dividing by L. The threshold q3(m) is obtained in similar manner. Comparator 626 then compares the two averaged quantities h3(m) and q3(m) for each frame m and provides the third voice detection signal d3(m) based on the comparison result. The comparison may be expressed as:
The third voice detection signal d3(m) is set to 1 to indicate that near-end voice is detected and set to 0 to indicate that near-end voice is not detected. However, the metric used by VAD3 is different from the metrics used by VAD1 and VAD2.
Within noise suppressor 280x, a noise estimator 710 receives the frequency-domain beam-formed signal B(k,m) from FFT unit 270, estimates the magnitude of the noise in the signal B(k,m), and provides a frequency-domain noise signal N1(k,m). The noise estimation may be performed using a minimum statistics based method or some other method, as is known in the art. The minimum statistics based method is described by R. Martin, in a paper entitled “Spectral subtraction based on minimum statistics,” EUSIPCO'94, pp. 1182–1185, September 1994. A noise estimator 720 receives the noise signal N1(k,m), the frequency-domain reference signal R(k,m), and the third voice detection signal d3(m). Noise estimator 720 determines a final estimate of the noise in the signal B(k,m) and provides a final noise estimate N2(k,m), which may be expressed as:
where γa1, γa2, γb1, and γb2 are constants and are selected such that γa1>γb1>0 and γb2>γa2>0. As shown in equation (19), the final noise estimate N2(k,m) is set equal to the sum of a first scaled noise estimate, γx1·N1(k,m), and a second scaled noise estimate, γx2·|R(k,m)|, where γx can be equal to γa or γb. The constants γa1, γa2, γb1, and γb2 are selected such that the final noise estimate N2(k,m) includes more of the noise estimate N1(k,m) and less of the reference signal magnitude |R(k,m)| when d3(m)=1, indicating that near-end voice is detected. Conversely, the final noise estimate N2(k,m) includes less of the noise estimate N1(k,m) and more of the reference signal magnitude |R(k,m)| when d3(m)=0, indicating that near-end voice is not detected.
A noise suppression gain computation unit 730 receives the frequency-domain beam-formed signal B(k,m), the final noise estimate N2(k,m), and the frequency-domain output signal Bo(k, m−1) for a prior frame from a delay unit 734. Computation unit 730 computes a noise suppression gain G(k,m) that is used to suppress additional noise and interference in the signal B(k,m).
To obtain the gain G(k,m), an SNR estimate G′SNR,B(k,m) for the beam-formed signal B(k,m) is first computed as follows:
The SNR estimate G′SNR,B(k,m) is then constrained to be a positive value or zero, as follows:
A final SNR estimate GSNR(k,m) is then computed as follows:
where λ is a positive constant that is selected such that 1>λ>0. As shown in equation (22), the final SNR estimate GSNR(k,m) includes two components. The first component is a scaled version of an SNR estimate for the output signal in the prior frame, i.e., λ·|Bo(k, m−1)|/N2(k,m). The second component is a scaled version of the constrained SNR estimate for the beam-formed signal, i.e., (1−λ)·GSNR,B(k,m). The constant λ determines the weighting for the two components that make up the final SNR estimate GSNR(k,m).
The gain G(k,m) is then computed as:
The gain G(k,m) is a real value and its magnitude is indicative of the amount of noise suppression to be performed. In particular, G(k,m) is a small value for more noise suppression and a large value for less noise suppression.
A multiplier 732 then multiples the frequency-domain beam-formed signal B(k,m) with the gain G(k,m) to provide the frequency-domain output signal Bo(k,m), which may be expressed as:
Bo(k,m)=B(k,m)·G(k,m) Eq (24)
Within FIR filter 810, the digital samples for the reference signal xref(n) are provided to M−1 series-coupled delay elements 812b through 812m, where M is the number of taps of the FIR filter. Each delay element provides one sample period of delay. The reference signal xref(n) and the outputs of delay elements 812b through 812m are provided to multipliers 814a through 814m, respectively. Each multiplier 814 also receives a respective filter coefficient hi(n) from coefficient calculation unit 820, multiplies its received samples with its filter coefficient hi(n), and provides output samples to a summer 816. For each sample period n, summer 816 sums the M output samples from multipliers 814a through 814m and provides a filtered sample for that sample period. The filtered sample xfir(n) for sample period n may be computed as:
where the symbol “*” denotes a complex conjugate. Summer 818 receives and subtracts the FIR signal xfir(n) from the input signal xin(n) and provides the output signal xout(n).
Coefficient calculation unit 820 provides the set of M coefficients for FIR filter 810, which is denoted as H*(n)=[h0*(n), h1*(n), . . . hM−1*(n)]. Unit 820 further updates these coefficients based on a particular adaptive algorithm, which may be a least mean square (LMS) algorithm, a normalized least mean square (NLMS) algorithm, a recursive least square (RLS) algorithm, a direct matrix inversion (DMI) algorithm, or some other algorithm. The NLMS and other algorithms are described by B. Widrow and S. D. Sterns in a book entitled “Adaptive Signal Processing,” Prentice-Hall Inc., Englewood Cliffs, N.J., 1986. The LMS, NLMS, RLS, DMI, and other adaptive algorithms are described by Simon Haykin in a book entitled “Adaptive Filter Theory”, 3rd edition, Prentice Hall, 1996. Coefficient update unit 820 also receives the control signal d(n) from VAD1 or VAD2, which controls the manner in which the filter coefficients are updated. For example, the filter coefficients may be updated only when voice activity is detected (i.e., when d(n)=1) and may be maintained when voice activity is not detected (i.e., when d(n)=0).
For clarity, a specific design for the small array microphone system has been described above, as shown in
The array microphone and noise suppression techniques described herein may be implemented by various means. For example, these techniques may be implemented in hardware, software, or a combination thereof. For a hardware implementation, the processing units used to implement the array microphone and noise suppression may be implemented within one or more application specific integrated circuits (ASICs), digital signal processors (DSPs), digital signal processing devices (DSPDs), programmable logic devices (PLDs), field programmable gate arrays (FPGAs), processors, controllers, micro-controllers, microprocessors, other electronic units designed to perform the functions described herein, or a combination thereof.
For a software implementation, the array microphone and noise suppression techniques may be implemented with modules (e.g., procedures, functions, and so on) that perform the functions described herein. The software codes may be stored in a memory unit (e.g., memory unit 1032 in
The previous description of the disclosed embodiments is provided to enable any person skilled in the art to make or use the present invention. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the generic principles defined herein may be applied to other embodiments without departing from the spirit or scope of the invention. Thus, the present invention is not intended to be limited to the embodiments shown herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.
This application claims the benefit of provisional U.S. Application Ser. No. 60/426,715, entitled “Small Array Microphone for Beam-forming,” filed Nov. 15, 2002, which is incorporated herein by reference in its entirety for all purposes. This application is further related to U.S. application Ser. No. 10/076,201, entitled “Noise Suppression for a Wireless Communication Device,” filed on Feb. 12, 2002, U.S. application Ser. No. 10/076,120, entitled “Noise Suppression for Speech Signal in an Automobile”, filed on Feb. 12, 2002, and U.S. patent application Ser. No. 10/371,150, entitled “Small Array Microphone for Acoustic Echo Cancellation and Noise Suppression,” filed Feb. 21, 2003, all of which are assigned to the assignee of the present application and incorporated herein by reference in their entirety for all purposes.
Number | Name | Date | Kind |
---|---|---|---|
6339758 | Kanazawa et al. | Jan 2002 | B1 |
6937980 | Krasny et al. | Aug 2005 | B2 |
20030027600 | Krasny et al. | Feb 2003 | A1 |
20030063759 | Brennan et al. | Apr 2003 | A1 |
Number | Date | Country | |
---|---|---|---|
60426715 | Nov 2002 | US |