System and method for providing noise suppression utilizing null processing noise subtraction

Information

  • Patent Grant
  • 9185487
  • Patent Number
    9,185,487
  • Date Filed
    Monday, June 30, 2008
    16 years ago
  • Date Issued
    Tuesday, November 10, 2015
    9 years ago
Abstract
Systems and methods for noise suppression using noise subtraction processing are provided. The noise subtraction processing comprises receiving at least a primary and a secondary acoustic signal. A desired signal component may be calculated and subtracted from the secondary acoustic signal to obtain a noise component signal. A determination may be made of a reference energy ratio and a prediction energy ratio. A determination may be made as to whether to adjust the noise component signal based partially on the reference energy ratio and partially on the prediction energy ratio. The noise component signal may be adjusted or frozen based on the determination. The noise component signal may then be removed from the primary acoustic signal to generate a noise subtracted signal which may be outputted.
Description
CROSS-REFERENCE TO RELATED APPLICATION

The present application is related to U.S. patent application Ser. No. 11/825,563, filed Jul. 6, 2007 and entitled “System and Method for Adaptive Intelligent Noise Suppression,” (now U.S. Pat. No. 8,774,844), and U.S. patent application Ser. No. 12/080,115, filed Mar. 31, 2008 and entitled “System and Method for Providing Close Microphone Adaptive Array Processing,” (now U.S. Pat. No. 8,204,252), both of which are herein incorporated by reference.


The present application is also related to U.S. patent application Ser. No. 11/343,524, filed Jan. 30, 2006 and entitled “System and Method for Utilizing Inter-Microphone Level Differences for Speech Enhancement,” (now U.S. Pat. No. 8,345,890), and U.S. patent application Ser. No. 11/699,732, filed Jan. 29, 2007 and entitled “System and Method for Utilizing Omni-Directional Microphones for Speech Enhancement,” (now U.S. Pat. No. 8,194,880), both of which are herein incorporated by reference.


BACKGROUND OF THE INVENTION

1. Field of Invention


The present invention relates generally to audio processing and more particularly to adaptive noise suppression of an audio signal.


2. Description of Related Art


Currently, there are many methods for reducing background noise in an adverse audio environment. One such method is to use a stationary noise suppression system. The stationary noise suppression system will always provide an output noise that is a fixed amount lower than the input noise. Typically, the stationary noise suppression is in the range of 12-13 decibels (dB). The noise suppression is fixed to this conservative level in order to avoid producing speech distortion, which will be apparent with higher noise suppression.


In order to provide higher noise suppression, dynamic noise suppression systems based on signal-to-noise ratios (SNR) have been utilized. This SNR may then be used to determine a suppression value. Unfortunately, SNR, by itself, is not a very good predictor of speech distortion due to existence of different noise types in the audio environment. SNR is a ratio of how much louder speech is than noise. However, speech may be a non-stationary signal which may constantly change and contain pauses. Typically, speech energy, over a period of time, will comprise a word, a pause, a word, a pause, and so forth. Additionally, stationary and dynamic noises may be present in the audio environment. The SNR averages all of these stationary and non-stationary speech and noise. There is no consideration as to the statistics of the noise signal; only what the overall level of noise is.


In some prior art systems, an enhancement filter may be derived based on an estimate of a noise spectrum. One common enhancement filter is the Wiener filter. Disadvantageously, the enhancement filter is typically configured to minimize certain mathematical error quantities, without taking into account a user's perception. As a result, a certain amount of speech degradation is introduced as a side effect of the noise suppression. This speech degradation will become more severe as the noise level rises and more noise suppression is applied. That is, as the SNR gets lower, lower gain is applied resulting in more noise suppression. This introduces more speech loss distortion and speech degradation.


Some prior art systems invoke a generalized side-lobe canceller. The generalized side-lobe canceller is used to identify desired signals and interfering signals comprised by a received signal. The desired signals propagate from a desired location and the interfering signals propagate from other locations. The interfering signals are subtracted from the received signal with the intention of cancelling interference.


Many noise suppression processes calculate a masking gain and apply this masking gain to an input signal. Thus, if an audio signal is mostly noise, a masking gain that is a low value may be applied (i.e., multiplied to) the audio signal. Conversely, if the audio signal is mostly desired sound, such as speech, a high value gain mask may be applied to the audio signal. This process is commonly referred to as multiplicative noise suppression.


SUMMARY OF THE INVENTION

Embodiments of the present invention overcome or substantially alleviate prior problems associated with noise suppression and speech enhancement. In exemplary embodiments, at least a primary and a secondary acoustic signal are received by a microphone array. The microphone array may comprise a close microphone array or a spread microphone array.


A noise component signal may be determined in each sub-band of signals received by the microphone by subtracting the primary acoustic signal weighted by a complex-valued coefficient σ from the secondary acoustic signal. The noise component signal, weighted by another complex-valued coefficient α, may then be subtracted from the primary acoustic signal resulting in an estimate of a target signal (i.e., a noise subtracted signal).


A determination may be made as to whether to adjust α. In exemplary embodiments, the determination may be based on a reference energy ratio (g1) and a prediction energy ratio (g2). The complex-valued coefficient α may be adapted when the prediction energy ratio is greater than the reference energy ratio to adjust the noise component signal. Conversely, the adaptation coefficient may be frozen when the prediction energy ratio is less than the reference energy ratio. The noise component signal may then be removed from the primary acoustic signal to generate a noise subtracted signal which may be outputted.





BRIEF DESCRIPTION OF THE DRAWINGS


FIG. 1 is an environment in which embodiments of the present invention may be practiced.



FIG. 2 is a block diagram of an exemplary audio device implementing embodiments of the present invention.



FIG. 3 is a block diagram of an exemplary audio processing system utilizing a spread microphone array.



FIG. 4 is a block diagram of an exemplary noise suppression system of the audio processing system of FIG. 3.



FIG. 5 is a block diagram of an exemplary audio processing system utilizing a close microphone array.



FIG. 6 is a block diagram of an exemplary noise suppression system of the audio processing system of FIG. 5.



FIG. 7
a is a block diagram of an exemplary noise subtraction engine.



FIG. 7
b is a schematic illustrating the operations of the noise subtraction engine.



FIG. 8 is a flowchart of an exemplary method for suppressing noise in an audio device.



FIG. 9 is a flowchart of an exemplary method for performing noise subtraction processing.





DESCRIPTION OF EXEMPLARY EMBODIMENTS

The present invention provides exemplary systems and methods for adaptive suppression of noise in an audio signal. Embodiments attempt to balance noise suppression with minimal or no speech degradation (i.e., speech loss distortion). In exemplary embodiments, noise suppression is based on an audio source location and applies a subtractive noise suppression process as opposed to a purely multiplicative noise suppression process.


Embodiments of the present invention may be practiced on any audio device that is configured to receive sound such as, but not limited to, cellular phones, phone handsets, headsets, and conferencing systems. Advantageously, exemplary embodiments are configured to provide improved noise suppression while minimizing speech distortion. While some embodiments of the present invention will be described in reference to operation on a cellular phone, the present invention may be practiced on any audio device.


Referring to FIG. 1, an environment in which embodiments of the present invention may be practiced is shown. A user acts as a speech (audio) source 102 to an audio device 104. The exemplary audio device 104 may include a microphone array. The microphone array may comprise a close microphone array or a spread microphone array.


In exemplary embodiments, the microphone array may comprise a primary microphone 106 relative to the audio source 102 and a secondary microphone 108 located a distance away from the primary microphone 106. While embodiments of the present invention will be discussed with regards to having two microphones 106 and 108, alternative embodiments may contemplate any number of microphones or acoustic sensors within the microphone array. In some embodiments, the microphones 106 and 108 may comprise omni-directional microphones.


While the microphones 106 and 108 receive sound (i.e., acoustic signals) from the audio source 102, the microphones 106 and 108 also pick up noise 110. Although the noise 110 is shown coming from a single location in FIG. 1, the noise 110 may comprise any sounds from one or more locations different than the audio source 102, and may include reverberations and echoes. The noise 110 may be stationary, non-stationary, or a combination of both stationary and non-stationary noise.


Referring now to FIG. 2, the exemplary audio device 104 is shown in more detail. In exemplary embodiments, the audio device 104 is an audio receiving device that comprises a processor 202, the primary microphone 106, the secondary microphone 108, an audio processing system 204, and an output device 206. The audio device 104 may comprise further components (not shown) necessary for audio device 104 operations. The audio processing system 204 will be discussed in more details in connection with FIG. 3.


In exemplary embodiments, the primary and secondary microphones 106 and 108 are spaced a distance apart in order to allow for an energy level difference between them. Upon reception by the microphones 106 and 108, the acoustic signals may be converted into electric signals (i.e., a primary electric signal and a secondary electric signal). The electric signals may, themselves, be converted by an analog-to-digital converter (not shown) into digital signals for processing in accordance with some embodiments. In order to differentiate the acoustic signals, the acoustic signal received by the primary microphone 106 is herein referred to as the primary acoustic signal, while the acoustic signal received by the secondary microphone 108 is herein referred to as the secondary acoustic signal.


The output device 206 is any device which provides an audio output to the user. For example, the output device 206 may comprise an earpiece of a headset or handset, or a speaker on a conferencing device.



FIG. 3 is a detailed block diagram of the exemplary audio processing system 204a according to one embodiment of the present invention. In exemplary embodiments, the audio processing system 204a is embodied within a memory device. The audio processing system 204a of FIG. 3 may be utilized in embodiments comprising a spread microphone array.


In operation, the acoustic signals received from the primary and secondary microphones 106 and 108 are converted to electric signals and processed through a frequency analysis module 302. In one embodiment, the frequency analysis module 302 takes the acoustic signals and mimics the frequency analysis of the cochlea (i.e., cochlear domain) simulated by a filter bank. In one example, the frequency analysis module 302 separates the acoustic signals into frequency sub-bands. A sub-band is the result of a filtering operation on an input signal where the bandwidth of the filter is narrower than the bandwidth of the signal received by the frequency analysis module 302. Alternatively, other filters such as short-time Fourier transform (STFT), sub-band filter banks, modulated complex lapped transforms, cochlear models, wavelets, etc., can be used for the frequency analysis and synthesis. Because most sounds (e.g., acoustic signals) are complex and comprise more than one frequency, a sub-band analysis on the acoustic signal determines what individual frequencies are present in the complex acoustic signal during a frame (e.g., a predetermined period of time). According to one embodiment, the frame is 8 ms long. Alternative embodiments may utilize other frame lengths or no frame at all. The results may comprise sub-band signals in a fast cochlea transform (FCT) domain.


Once the sub-band signals are determined, the sub-band signals are forwarded to a noise subtraction engine 304. The exemplary noise subtraction engine 304 is configured to adaptively subtract out a noise component from the primary acoustic signal for each sub-band. As such, output of the noise subtraction engine 304 is a noise subtracted signal comprised of noise subtracted sub-band signals. The noise subtraction engine 304 will be discussed in more detail in connection with FIG. 7a and FIG. 7b. It should be noted that the noise subtracted sub-band signals may comprise desired audio that is speech or non-speech (e.g., music). The results of the noise subtraction engine 304 may be output to the user or processed through a further noise suppression system (e.g., the noise suppression engine 306). For purposes of illustration, embodiments of the present invention will discuss embodiments whereby the output of the noise subtraction engine 304 is processed through a further noise suppression system.


The noise subtracted sub-band signals along with the sub-band signals of the secondary acoustic signal are then provided to the noise suppression engine 306a. According to exemplary embodiments, the noise suppression engine 306a generates a gain mask to be applied to the noise subtracted sub-band signals in order to further reduce noise components that remain in the noise subtracted speech signal. The noise suppression engine 306a will be discussed in more detail in connection with FIG. 4 below.


The gain mask determined by the noise suppression engine 306a may then be applied to the noise subtracted signal in a masking module 308. Accordingly, each gain mask may be applied to an associated noise subtracted frequency sub-band to generate masked frequency sub-bands. As depicted in FIG. 3, a multiplicative noise suppression system 312a comprises the noise suppression engine 306a and the masking module 308.


Next, the masked frequency sub-bands are converted back into time domain from the cochlea domain. The conversion may comprise taking the masked frequency sub-bands and adding together phase shifted signals of the cochlea channels in a frequency synthesis module 310. Alternatively, the conversion may comprise taking the masked frequency sub-bands and multiplying these with an inverse frequency of the cochlea channels in the frequency synthesis module 310. Once conversion is completed, the synthesized acoustic signal may be output to the user.


Referring now to FIG. 4, the noise suppression engine 306a of FIG. 3 is illustrated. The exemplary noise suppression engine 306a comprises an energy module 402, an inter-microphone level difference (ILD) module 404, an adaptive classifier 406, a noise estimate module 408, and an adaptive intelligent suppression (AIS) generator 410. It should be noted that the noise suppression engine 306a is exemplary and may comprise other combinations of modules such as that shown and described in U.S. patent application Ser. No. 11/343,524, which is incorporated by reference.


According to an exemplary embodiment of the present invention, the AIS generator 410 derives time and frequency varying gains or gain masks used by the masking module 308 to suppress noise and enhance speech in the noise subtracted signal. In order to derive the gain masks, however, specific inputs are needed for the AIS generator 410. These inputs comprise a power spectral density of noise (i.e., noise spectrum), a power spectral density of the noise subtracted signal (herein referred to as the primary spectrum), and an inter-microphone level difference (ILD).


According to exemplary embodiment, the noise subtracted signal (c′(k)) resulting from the noise subtraction engine 304 and the secondary acoustic signal (f′(k)) are forwarded to the energy module 402 which computes energy/power estimates during an interval of time for each frequency band (i.e., power estimates) of an acoustic signal. As can be seen in FIG. 7b, f′(k) may optionally be equal to f(k). As a result, the primary spectrum (i.e., the power spectral density of the noise subtracted signal) across all frequency bands may be determined by the energy module 402. This primary spectrum may be supplied to the AIS generator 410 and the ILD module 404 (discussed further herein). Similarly, the energy module 402 determines a secondary spectrum (i.e., the power spectral density of the secondary acoustic signal) across all frequency bands which is also supplied to the ILD module 404. More details regarding the calculation of power estimates and power spectrums can be found in co-pending U.S. patent application Ser. No. 11/343,524 and co-pending U.S. patent application Ser. No. 11/699,732, which are incorporated by reference.


In two microphone embodiments, the power spectrums are used by an inter-microphone level difference (ILD) module 404 to determine an energy ratio between the primary and secondary microphones 106 and 108. In exemplary embodiments, the ILD may be a time and frequency varying ILD. Because the primary and secondary microphones 106 and 108 may be oriented in a particular way, certain level differences may occur when speech is active and other level differences may occur when noise is active. The ILD is then forwarded to the adaptive classifier 406 and the AIS generator 410. More details regarding one embodiment for calculating ILD may be can be found in co-pending U.S. patent application Ser. No. 11/343,524 and co-pending U.S. patent application Ser. No. 11/699,732. In other embodiments, other forms of ILD or energy differences between the primary and secondary microphones 106 and 108 may be utilized. For example, a ratio of the energy of the primary and secondary microphones 106 and 108 may be used. It should also be noted that alternative embodiments may use cues other then ILD for adaptive classification and noise suppression (i.e., gain mask calculation). For example, noise floor thresholds may be used. As such, references to the use of ILD may be construed to be applicable to other cues.


The exemplary adaptive classifier 406 is configured to differentiate noise and distractors (e.g., sources with a negative ILD) from speech in the acoustic signal(s) for each frequency band in each frame. The adaptive classifier 406 is considered adaptive because features (e.g., speech, noise, and distractors) change and are dependent on acoustic conditions in the environment. For example, an ILD that indicates speech in one situation may indicate noise in another situation. Therefore, the adaptive classifier 406 may adjust classification boundaries based on the ILD.


According to exemplary embodiments, the adaptive classifier 406 differentiates noise and distractors from speech and provides the results to the noise estimate module 408 which derives the noise estimate. Initially, the adaptive classifier 406 may determine a maximum energy between channels at each frequency. Local ILDs for each frequency are also determined. A global ILD may be calculated by applying the energy to the local ILDs. Based on the newly calculated global ILD, a running average global ILD and/or a running mean and variance (i.e., global cluster) for ILD observations may be updated. Frame types may then be classified based on a position of the global ILD with respect to the global cluster. The frame types may comprise source, background, and distractors.


Once the frame types are determined, the adaptive classifier 406 may update the global average running mean and variance (i.e., cluster) for the source, background, and distractors. In one example, if the frame is classified as source, background, or distracter, the corresponding global cluster is considered active and is moved toward the global ILD. The global source, background, and distractor global clusters that do not match the frame type are considered inactive. Source and distractor global clusters that remain inactive for a predetermined period of time may move toward the background global cluster. If the background global cluster remains inactive for a predetermined period of time, the background global cluster moves to the global average.


Once the frame types are determined, the adaptive classifier 406 may also update the local average running mean and variance (i.e., cluster) for the source, background, and distractors. The process of updating the local active and inactive clusters is similar to the process of updating the global active and inactive clusters.


Based on the position of the source and background clusters, points in the energy spectrum are classified as source or noise; this result is passed to the noise estimate module 408.


In an alternative embodiment, an example of an adaptive classifier 406 comprises one that tracks a minimum ILD in each frequency band using a minimum statistics estimator. The classification thresholds may be placed a fixed distance (e.g., 3 dB) above the minimum ILD in each band. Alternatively, the thresholds may be placed a variable distance above the minimum ILD in each band, depending on the recently observed range of ILD values observed in each band. For example, if the observed range of ILDs is beyond 6 dB, a threshold may be place such that it is midway between the minimum and maximum ILDs observed in each band over a certain specified period of time (e.g., 2 seconds). The adaptive classifier is further discussed in the U.S. nonprovisional application entitled “System and Method for Adaptive Intelligent Noise Suppression,” Ser. No. 11/825,563, filed Jul. 6, 2007, which is incorporated by reference.


In exemplary embodiments, the noise estimate is based on the acoustic signal from the primary microphone 106 and the results from the adaptive classifier 406. The exemplary noise estimate module 408 generates a noise estimate which is a component that can be approximated mathematically by

N(t,ω)=λ1(t,ω)E1(t,ω)+(1−λ1(t,ω))min[N(t−1,ω),E1(t,ω)]

according to one embodiment of the present invention. As shown, the noise estimate in this embodiment is based on minimum statistics of a current energy estimate of the primary acoustic signal, E1(t,ω) and a noise estimate of a previous time frame, N(t−1, ω). As a result, the noise estimation is performed efficiently and with low latency.


λ1(t,ω) in the above equation may be derived from the ILD approximated by the ILD module 404, as








λ
I



(

t
,
ω

)


=

{





0





if






ILD


(

t
,
ω

)



<
threshold







1





if






ILD


(

t
,
ω

)



>
threshold










That is, when the primary microphone 106 is smaller than a threshold value (e.g., threshold=0.5) above which speech is expected to be, λ1 is small, and thus the noise estimate module 408 follows the noise closely. When ILD starts to rise (e.g., because speech is present within the large ILD region), λ1 increases. As a result, the noise estimate module 408 slows down the noise estimation process and the speech energy does not contribute significantly to the final noise estimate. Alternative embodiments, may contemplate other methods for determining the noise estimate or noise spectrum. The noise spectrum (i.e., noise estimates for all frequency bands of an acoustic signal) may then be forwarded to the AIS generator 410.


The AIS generator 410 receives speech energy of the primary spectrum from the energy module 402. This primary spectrum may also comprise some residual noise after processing by the noise subtraction engine 304. The AIS generator 410 may also receive the noise spectrum from the noise estimate module 408. Based on these inputs and an optional ILD from the ILD module 404, a speech spectrum may be inferred. In one embodiment, the speech spectrum is inferred by subtracting the noise estimates of the noise spectrum from the power estimates of the primary spectrum. Subsequently, the AIS generator 410 may determine gain masks to apply to the primary acoustic signal. More detailed discussion of the AIS generator 410 may be found in U.S. patent application Ser. No. 11/825,563 entitled “System and Method for Adaptive Intelligent Noise Suppression,” which is incorporated by reference. In exemplary embodiments, the gain mask output from the AIS generator 410, which is time and frequency dependent, will maximize noise suppression while constraining speech loss distortion.


It should be noted that the system architecture of the noise suppression engine 306a is exemplary. Alternative embodiments may comprise more components, less components, or equivalent components and still be within the scope of embodiments of the present invention. Various modules of the noise suppression engine 306a may be combined into a single module. For example, the functionalities of the ILD module 404 may be combined with the functions of the energy module 402.


Referring now to FIG. 5, a detailed block diagram of an alternative audio processing system 204b is shown. In contrast to the audio processing system 204a of FIG. 3, the audio processing system 204b of FIG. 5 may be utilized in embodiments comprising a close microphone array. The functions of the frequency analysis module 302, masking module 308, and frequency synthesis module 310 are identical to those described with respect to the audio processing system 204a of FIG. 3 and will not be discussed in detail.


The sub-band signals determined by the frequency analysis module 302 may be forwarded to the noise subtraction engine 304 and an array processing engine 502. The exemplary noise subtraction engine 304 is configured to adaptively subtract out a noise component from the primary acoustic signal for each sub-band. As such, output of the noise subtraction engine 304 is a noise subtracted signal comprised of noise subtracted sub-band signals. In the present embodiment, the noise subtraction engine 304 also provides a null processing (NP) gain to the noise suppression engine 306a. The NP gain comprises an energy ratio indicating how much of the primary signal has been cancelled out of the noise subtracted signal. If the primary signal is dominated by noise, then NP gain will be large. In contrast, if the primary signal is dominated by speech, NP gain will be close to zero. The noise subtraction engine 304 will be discussed in more detail in connection with FIG. 7a and FIG. 7b below.


In exemplary embodiments, the array processing engine 502 is configured to adaptively process the sub-band signals of the primary and secondary signals to create directional patterns (i.e., synthetic directional microphone responses) for the close microphone array (e.g., the primary and secondary microphones 106 and 108). The directional patterns may comprise a forward-facing cardioid pattern based on the primary acoustic (sub-band) signals and a backward-facing cardioid pattern based on the secondary (sub-band) acoustic signal. In one embodiment, the sub-band signals may be adapted such that a null of the backward-facing cardioid pattern is directed towards the audio source 102. More details regarding the implementation and functions of the array processing engine 502 may be found (referred to as the adaptive array processing engine) in U.S. patent application Ser. No. 12/080,115 entitled “System and Method for Providing Close Microphone Array Noise Reduction,” which is incorporated by reference. The cardioid signals (i.e., a signal implementing the forward-facing cardioid pattern and a signal implementing the backward-facing cardioid pattern) are then provided to the noise suppression engine 306b by the array processing engine 502.


The noise suppression engine 306b receives the NP gain along with the cardioid signals. According to exemplary embodiments, the noise suppression engine 306b generates a gain mask to be applied to the noise subtracted sub-band signals from the noise subtraction engine 304 in order to further reduce any noise components that may remain in the noise subtracted speech signal. The noise suppression engine 306b will be discussed in more detail in connection with FIG. 6 below.


The gain mask determined by the noise suppression engine 306b may then be applied to the noise subtracted signal in the masking module 308. Accordingly, each gain mask may be applied to an associated noise subtracted frequency sub-band to generate masked frequency sub-bands. Subsequently, the masked frequency sub-bands are converted back into time domain from the cochlea domain by the frequency synthesis module 310. Once conversion is completed, the synthesized acoustic signal may be output to the user. As depicted in FIG. 5, a multiplicative noise suppression system 312b comprises the array processing engine 502, the noise suppression engine 306b, and the masking module 308.


Referring now to FIG. 6, the exemplary noise suppression engine 306b is shown in more detail. The exemplary noise suppression engine 306b comprises the energy module 402, the inter-microphone level difference (ILD) module 404, the adaptive classifier 406, the noise estimate module 408, and the adaptive intelligent suppression (AIS) generator 410. It should be noted that the various modules of the noise suppression engine 306b functions similar to the modules in the noise suppression engine 306a.


In the present embodiment, the primary acoustic signal (c″(k)) and the secondary acoustic signal (f″(k)) are received by the energy module 402 which computes energy/power estimates during an interval of time for each frequency band (i.e., power estimates) of an acoustic signal. As a result, the primary spectrum (i.e., the power spectral density of the primary sub-band signals) across all frequency bands may be determined by the energy module 402. This primary spectrum may be supplied to the AIS generator 410 and the ILD module 404. Similarly, the energy module 402 determines a secondary spectrum (i.e., the power spectral density of the secondary sub-band signal) across all frequency bands which is also supplied to the ILD module 404. More details regarding the calculation of power estimates and power spectrums can be found in co-pending U.S. patent application Ser. No. 11/343,524 and co-pending U.S. patent application Ser. No. 11/699,732, which are incorporated by reference.


As previously discussed, the power spectrums may be used by the ILD module 404 to determine an energy difference between the primary and secondary microphones 106 and 108. The ILD may then be forwarded to the adaptive classifier 406 and the AIS generator 410. In alternative embodiments, other forms of ILD or energy differences between the primary and secondary microphones 106 and 108 may be utilized. For example, a ratio of the energy of the primary and secondary microphones 106 and 108 may be used. It should also be noted that alternative embodiments may use cues other then ILD for adaptive classification and noise suppression (i.e., gain mask calculation). For example, noise floor thresholds may be used. As such, references to the use of ILD may be construed to be applicable to other cues.


The exemplary adaptive classifier 406 and noise estimate module 408 perform the same functions as that described in accordance with FIG. 4. That is, the adaptive classifier differentiates noise and distractors from speech and provides the results to the noise estimate module 408 which derives the noise estimate.


The AIS generator 410 receives speech energy of the primary spectrum from the energy module 402. The AIS generator 410 may also receive the noise spectrum from the noise estimate module 408. Based on these inputs and an optional ILD from the ILD module 404, a speech spectrum may be inferred. In one embodiment, the speech spectrum is inferred by subtracting the noise estimates of the noise spectrum from the power estimates of the primary spectrum. Additionally, the AIS generator 410 uses the NP gain, which indicates how much noise has already been cancelled by the time the signal reaches the noise suppression engine 306b (i.e., the multiplicative mask) to determine gain masks to apply to the primary acoustic signal. In one example, as the NP gain increases, the estimated SNR for the inputs decreases. In exemplary embodiments, the gain mask output from the AIS generator 410, which is time and frequency dependent, may maximize noise suppression while constraining speech loss distortion.


It should be noted that the system architecture of the noise suppression engine 306b is exemplary. Alternative embodiments may comprise more components, less components, or equivalent components and still be within the scope of embodiments of the present invention.



FIG. 7
a is a block diagram of an exemplary noise subtraction engine 304. The exemplary noise subtraction engine 304 is configured to suppress noise using a subtractive process. The noise subtraction engine 304 may determine a noise subtracted signal by initially subtracting out a desired component (e.g., the desired speech component) from the primary signal in a first branch, thus resulting in a noise component. Adaptation may then be performed in a second branch to cancel out the noise component from the primary signal. In exemplary embodiments, the noise subtraction engine 304 comprises a gain module 702, an analysis module 704, an adaptation module 706, and at least one summing module 708 configured to perform signal subtraction. The functions of the various modules 702-708 will be discussed in connection with FIG. 7a and further illustrated in operation in connection with FIG. 7b.


Referring to FIG. 7a, the exemplary gain module 702 is configured to determine various gains used by the noise subtraction engine 304. For purposes of the present embodiment, these gains represent energy ratios. In the first branch, a reference energy ratio (g1) of how much of the desired component is removed from the primary signal may be determined. In the second branch, a prediction energy ratio (g2) of how much the energy has been reduced at the output of the noise subtraction engine 304 from the result of the first branch may be determined. Additionally, an energy ratio (i.e., NP gain) may be determined that represents the energy ratio indicating how much noise has been canceled from the primary signal by the noise subtraction engine 304. As previously discussed, NP gain may be used by the AIS generator 410 in the close microphone embodiment to adjust the gain mask.


The exemplary analysis module 704 is configured to perform the analysis in the first branch of the noise subtraction engine 304, while the exemplary adaptation module 706 is configured to perform the adaptation in the second branch of the noise subtraction engine 304.


Referring to FIG. 7b, a schematic illustrating the operations of the noise subtraction engine 304 is shown. Sub-band signals of the primary microphone signal c(k) and secondary microphone signal f(k) are received by the noise subtraction engine 304 where k represents a discrete time or sample index. c(k) represents a superposition of a speech signal s(k) and a noise signal n(k). f(k) is modeled as a superposition of the speech signal s(k), scaled by a complex-valued coefficient σ, and the noise signal n(k), scaled by a complex-valued coefficient ν. ν represents how much of the noise in the primary signal is in the secondary signal. In exemplary embodiments, ν is unknown since a source of the noise may be dynamic.


In exemplary embodiments, σ is a fixed coefficient that represents a location of the speech (e.g., an audio source location). In accordance with exemplary embodiments, σ may be determined through calibration. Tolerances may be included in the calibration by calibrating based on more than one position. For a close microphone, a magnitude of a may be close to one. For spread microphones, the magnitude of σ may be dependent on where the audio device 102 is positioned relative to the speaker's mouth. The magnitude and phase of the σ may represent an inter-channel cross-spectrum for a speaker's mouth position at a frequency represented by the respective sub-band (e.g., Cochlea tap). Because the noise subtraction engine 304 may have knowledge of what σ is, the analysis module 704 may apply σ to the primary signal (i.e., σ(s(k)+n(k)) and subtract the result from the secondary signal (i.e., σs(k)+ν(k)) in order to cancel out the speech component σ s(k) (i.e., the desired component) from the secondary signal resulting in a noise component out of the summing module 708. In an embodiment where there is not speech, α is approximately 1/(ν−σ), and the adaptation module 706 may freely adapt.


If the speaker's mouth position is adequately represented by σ, then f(k)−σc(k)=(ν−σ)n(k). This equation indicates that signal at the output of the summing module 708 being fed into the adaptation module 706 (which, in turn, applies an adaptation coefficient α(k)) may be devoid of a signal originating from a position represented by σ (e.g., the desired speech signal). In exemplary embodiments, the analysis module 704 applies σ to the secondary signal f(k) and subtracts the result from c(k). Remaining signal (referred to herein as “noise component signal”) from the summing module 708 may be canceled out in the second branch.


The adaptation module 706 may adapt when the primary signal is dominated by audio sources 102 not in the speech location (represented by σ). If the primary signal is dominated by a signal originating from the speech location as represented by σ, adaptation may be frozen. In exemplary embodiments, the adaptation module 706 may adapt using one of a common least-squares method in order to cancel the noise component n(k) from the signal c(k). The coefficient may be update at a frame rate according to on embodiment.


In an embodiment where n(k) is white and a cross-correlation between s(k) and n(k) is zero within a frame, adaptation may happen every frame with the noise n(k) being perfectly cancelled and the speech s(k) being perfectly unaffected. However, it is unlikely that these conditions may be met in reality, especially if the frame size is short. As such, it is desirable to apply constraints on adaptation. In exemplary embodiments, the adaptation coefficient α(k) may be updated on a per-tap/per-frame basis when the reference energy ratio g1 and the prediction energy ratio g2 satisfy the follow condition:

g2·γ>g1

where γ>0. Assuming, for example, that {circumflex over (σ)}(k)=σ, α(k)=1/(ν−σ), and s(k) and n(k) are uncorrelated, the following may be obtained:










g
1

=



E


{


(


s


(
k
)


+

n


(
k
)



)

2

}








v
-
σ



2

·
E



{


n
2



(
k
)


}



=





S
+
N






v
-
σ



2

·
N







and










g
2

=








v
-
σ



2

·
E



{


n
2



(
k
)


}



E


{


s
2



(
k
)


}



=







v
-
σ



2

·

N
S




,








where E{ . . . } is an expected value, S is a signal energy, and N is a noise energy. From the previous three equations, the following may be obtained:

SNR2+SNR<γ2|ν−σ|4,

where SNR=S/N. If the noise is in the same location as the target speech (i.e., σ=ν), this condition may not be met, so regardless of the SNR, adaptation may never happen. The further away from the target location the source is, the greater |ν−σ|4 and the larger the SNR is allowed to be while there is still adaptation attempting to cancel the noise.


In exemplary embodiments, adaptation may occur in frames where more signal is canceled in the second branch as opposed to the first branch. Thus, energies may be calculated after the first branch by the gain module 702 and g1 determined. An energy calculation may also be performed in order to determine g2 which may indicate if α is allowed to adapt. If γ2|ν−σ|4>SNR2+SNR4 is true, then adaptation of a may be performed. However, if this equation is not true, then α is not adapted.


The coefficient γ may be chosen to define a boundary between adaptation and non-adaptation of α. In an embodiment where a far-field source at 90 degree angle relative to a straight line between the microphones 106 and 108. In this embodiment, the signal may have equal power and zero phase shift between both microphones 106 and 108 (e.g., ν=1). If the SNR=1, then γ2|ν−σ|4=2, which is equivalent to γ=sqrt(2)/|1−σ|4.


Lowering γ relative to this value may improve protection of the near-end source from cancellation at the expense of increased noise leakage; raising γ has an opposite effect. It should be noted that in the microphones 106 and 108, ν=1 may not be a good enough approximation of the far-field/90 degrees situation and may have to substituted by a value obtained from calibration measurements.



FIG. 8 is a flowchart 800 of an exemplary method for suppressing noise in an audio device. In step 802, audio signals are received by the audio device 102. In exemplary embodiments, a plurality of microphones (e.g., primary and secondary microphones 106 and 108) receive the audio signals. The plurality of microphones may comprise a close microphone array or a spread microphone array.


In step 804, the frequency analysis on the primary and secondary acoustic signals may be performed. In one embodiment, the frequency analysis module 302 utilizes a filter bank to determine frequency sub-bands for the primary and secondary acoustic signals.


Noise subtraction processing is performed in step 806. Step 806 will be discussed in more detail in connection with FIG. 9 below.


Noise suppression processing may then be performed in step 808. In one embodiment, the noise suppression processing may first compute an energy spectrum for the primary or noise subtracted signal and the secondary signal. An energy difference between the two signals may then be determined. Subsequently, the speech and noise components may be adaptively classified according to one embodiment. A noise spectrum may then be determined. In one embodiment, the noise estimate may be based on the noise component. Based on the noise estimate, a gain mask may be adaptively determined.


The gain mask may then be applied in step 810. In one embodiment, the gain mask may be applied by the masking module 308 on a per sub-band signal basis. In some embodiments, the gain mask may be applied to the noise subtracted signal. The sub-bands signals may then be synthesized in step 812 to generate the output. In one embodiment, the sub-band signals may be converted back to the time domain from the frequency domain. Once converted, the audio signal may be output to the user in step 814. The output may be via a speaker, earpiece, or other similar devices.


Referring now to FIG. 9, a flowchart of an exemplary method for performing noise subtraction processing (step 806) is shown. In step 902, the frequency analyzed signals (e.g., frequency sub-band signals or primary signal) are received by the noise subtraction engine 304. The primary acoustic signal may be represented as c(k)=s(k)+n(k) where s(k) represents the desired signal (e.g., speech signal) and n(k) represents the noise signal. The secondary frequency analyzed signal (e.g., secondary signal) may be represented as f(k)=σs(k)+νn(k).


In step 904, σ may be applied to the primary signal by the analysis module 704. The result of the application of σ to the primary signal may then be subtracted from the secondary signal in step 906 by the summing module 708. The result comprises a noise component signal.


In step 908, the gains may be calculated by the gain module 702. These gains represent energy ratios of the various signals. In the first branch, a reference energy ratio (g1) of how much of the desired component is removed from the primary signal may be determined. In the second branch, a prediction energy ratio (g2) of how much the energy has been reduce at the output of the noise subtraction engine 304 from the result of the first branch may be determined.


In step 910, a determination is made as to whether α should be adapted. In accordance with one embodiment if SNR2+SNR<γ2|ν−σ|4 is true, then adaptation of α may be performed in step 912. However, if this equation is not true, then α is not adapted but frozen in step 914.


The noise component signal, whether adapted or not, is subtracted from the primary signal in step 916 by the summing module 708. The result is a noise subtracted signal. In some embodiments, the noise subtracted signal may be provided to the noise suppression engine 306 for further noise suppression processing via a multiplicative noise suppression process. In other embodiments, the noise subtracted signal may be output to the user without further noise suppression processing. It should be noted that more than one summing module 708 may be provided (e.g., one for each branch of the noise subtraction engine 304).


In step 918, the NP gain may be calculated. The NP gain comprises an energy ratio indicating how much of the primary signal has been cancelled out of the noise subtracted signal. It should be noted that step 918 may be optional (e.g., in close microphone systems).


The above-described modules may be comprised of instructions that are stored in storage media such as a machine readable medium (e.g., a computer readable medium). The instructions may be retrieved and executed by the processor 202. Some examples of instructions include software, program code, and firmware. Some examples of storage media comprise memory devices and integrated circuits. The instructions are operational when executed by the processor 202 to direct the processor 202 to operate in accordance with embodiments of the present invention. Those skilled in the art are familiar with instructions, processors, and storage media.


The present invention is described above with reference to exemplary embodiments. It will be apparent to those skilled in the art that various modifications may be made and other embodiments may be used without departing from the broader scope of the present invention. For example, the microphone array discussed herein comprises a primary and secondary microphone 106 and 108. However, alternative embodiments may contemplate utilizing more microphones in the microphone array. Therefore, there and other variations upon the exemplary embodiments are intended to be covered by the present invention.

Claims
  • 1. A method for suppressing noise, comprising: receiving at least a primary acoustic signal from a primary microphone and a secondary acoustic signal from a different, secondary microphone;applying a coefficient to the primary acoustic signal to generate a desired signal component, the coefficient representing a source location, the desired signal component not being a function of the secondary acoustic signal;subtracting the desired signal component from the secondary acoustic signal to obtain a noise component signal;performing a first determination of at least one energy ratio related to the desired signal component and the noise component signal;performing a second determination of whether to adjust the noise component signal based on the at least one energy ratio;adjusting the noise component signal based on the second determination;subtracting the adjusted noise component signal from the primary acoustic signal to generate a noise subtracted signal; andoutputting the noise subtracted signal.
  • 2. The method of claim 1 wherein the at least one energy ratio comprises a reference energy ratio and a prediction energy ratio.
  • 3. The method of claim 2 further comprising adapting an adaptation coefficient applied to the noise component signal when the prediction energy ratio is greater than the reference energy ratio.
  • 4. The method of claim 2 further comprising freezing an adaptation coefficient applied to the noise component signal when the prediction energy ratio is less than the reference energy ratio.
  • 5. The method of claim 1 further comprising determining a NP gain based on the at least one energy ratio, the NP gain indicating how much of the primary acoustic signal has been cancelled out of the noise subtracted signal.
  • 6. The method of claim 5 further comprising providing the NP gain to a multiplicative noise suppression system.
  • 7. The method of claim 1 wherein the primary and secondary acoustic signals are separated into sub-band signals.
  • 8. The method of claim 1 wherein outputting the noise subtracted signal comprises outputting the noise subtracted signal to a multiplicative noise suppression system.
  • 9. The method of claim 8 wherein the multiplicative noise suppression system comprises generating a gain mask based at least on the noise subtracted signal.
  • 10. The method of claim 9 further comprising applying the gain mask to the noise subtracted signal to generate an audio output signal.
  • 11. A system for suppressing noise, comprising: a microphone array configured to receive at least a primary acoustic signal from a primary microphone and a secondary acoustic signal from a different, secondary microphone;an analysis module configured to generate a desired signal component which may be subtracted from the secondary acoustic signal to obtain a noise component signal, the analysis module being further configured to apply a coefficient to the primary acoustic signal to generate the desired signal component, the coefficient representing a source location, the desired signal component not being a function of the secondary acoustic signal;a gain module configured to perform a first determination of at least one energy ratio related to the desired signal component and the noise component signal;an adaptation module configured to perform a second determination of whether to adjust the noise component signal based on the at least one energy ratio, the adaption module further configured to adjust the noise component signal based on the second determination; andat least one summing module configured to subtract the desired signal component from the adjusted secondary acoustic signal and to subtract the noise component signal from the primary acoustic signal to generate a noise subtracted signal.
  • 12. The system of claim 11 wherein the at least one energy ratio comprises a reference energy ratio and a prediction energy ratio.
  • 13. The system of claim 12 wherein the adaptation module is configured to adapt an adaptation coefficient applied to the noise component signal when the prediction energy ratio is greater than the reference energy ratio.
  • 14. The system of claim 12 wherein the adaptation module is configured to freeze an adaptation coefficient applied to the noise component signal when the prediction energy ratio is less than the reference energy ratio.
  • 15. The system of claim 11 wherein further comprising a gain module configured to determine a NP gain based on the at least one energy ratio, the NP gain indicating how much of the primary acoustic signal has been cancelled out of the noise subtracted signal.
  • 16. A non-transitory machine readable storage medium having embodied thereon a program, the program providing instructions executable by a processor for suppressing noise using noise subtraction processing method, the method comprising: receiving at least a primary acoustic signal from a primary microphone and a secondary acoustic signal from a different, secondary microphone;applying a coefficient to the primary acoustic signal to generate a desired signal component, the coefficient representing a source location, the desired signal component not being a function of the secondary acoustic signal;subtracting the desired signal component from the secondary acoustic signal to obtain a noise component signal;performing a first determination of at least one energy ratio related to the desired signal component and the noise component signal;performing a second determination of whether to adjust the noise component signal based on the at least one energy ratio;adjusting the noise component signal based on the second determination;subtracting the adjusted noise component signal from the primary acoustic signal to generate a noise subtracted signal; andoutputting the noise subtracted signal.
  • 17. The non-transitory machine readable storage medium of claim 16 wherein the at least one energy ratio comprises a reference energy ratio and a prediction energy ratio.
  • 18. The non-transitory machine readable storage medium of claim 17 wherein the method further comprises adapting an adaptation coefficient applied to the noise component signal when the prediction energy ratio is greater than the reference energy ratio.
  • 19. The non-transitory machine readable storage medium of claim 17 wherein the method further comprises freezing an adaptation coefficient applied to the noise component signal when the prediction energy ratio is less than the reference energy ratio.
  • 20. A method for suppressing noise, comprising: receiving at least a primary acoustic signal from a primary microphone and a secondary acoustic signal from a different, secondary microphone;applying a coefficient to the primary acoustic signal to generate a desired signal component, the coefficient representing a source location, the desired signal component not being a function of the secondary acoustic signal;subtracting the desired signal component from the secondary acoustic signal to obtain a noise component signal;performing a first determination of at least one energy ratio related to the desired signal component and the noise component signal, wherein the at least one energy ratio comprises a reference energy ratio and a prediction energy ratio;performing a second determination of whether to adjust the noise component signal based on the at least one energy ratio;adjusting the noise component signal based on the second determination; andsubtracting adjusted the noise component signal from the primary acoustic signal to generate a noise subtracted signal.
US Referenced Citations (268)
Number Name Date Kind
3976863 Engel Aug 1976 A
3978287 Fletcher et al. Aug 1976 A
4137510 Iwahara Jan 1979 A
4433604 Ott Feb 1984 A
4516259 Yato et al. May 1985 A
4535473 Sakata Aug 1985 A
4536844 Lyon Aug 1985 A
4581758 Coker et al. Apr 1986 A
4628529 Borth et al. Dec 1986 A
4630304 Borth et al. Dec 1986 A
4649505 Zinser, Jr. et al. Mar 1987 A
4658426 Chabries et al. Apr 1987 A
4674125 Carlson et al. Jun 1987 A
4718104 Anderson Jan 1988 A
4811404 Vilmur et al. Mar 1989 A
4812996 Stubbs Mar 1989 A
4864620 Bialick Sep 1989 A
4920508 Yassaie et al. Apr 1990 A
5027410 Williamson et al. Jun 1991 A
5054085 Meisel et al. Oct 1991 A
5058419 Nordstrom et al. Oct 1991 A
5099738 Hotz Mar 1992 A
5119711 Bell et al. Jun 1992 A
5142961 Paroutaud Sep 1992 A
5150413 Nakatani et al. Sep 1992 A
5175769 Hejna, Jr. et al. Dec 1992 A
5187776 Yanker Feb 1993 A
5208864 Kaneda May 1993 A
5210366 Sykes, Jr. May 1993 A
5224170 Waite, Jr. Jun 1993 A
5230022 Sakata Jul 1993 A
5319736 Hunt Jun 1994 A
5323459 Hirano Jun 1994 A
5341432 Suzuki et al. Aug 1994 A
5371800 Komatsu et al. Dec 1994 A
5381473 Andrea et al. Jan 1995 A
5381512 Holton et al. Jan 1995 A
5400409 Linhard Mar 1995 A
5402493 Goldstein Mar 1995 A
5402496 Soli et al. Mar 1995 A
5471195 Rickman Nov 1995 A
5473702 Yoshida et al. Dec 1995 A
5473759 Slaney et al. Dec 1995 A
5479564 Vogten et al. Dec 1995 A
5502663 Lyon Mar 1996 A
5544250 Urbanski Aug 1996 A
5574824 Slyh et al. Nov 1996 A
5583784 Kapust et al. Dec 1996 A
5587998 Velardo, Jr. et al. Dec 1996 A
5590241 Park et al. Dec 1996 A
5602962 Kellermann Feb 1997 A
5675778 Jones Oct 1997 A
5682463 Allen et al. Oct 1997 A
5694474 Ngo et al. Dec 1997 A
5706395 Arslan et al. Jan 1998 A
5717829 Takagi Feb 1998 A
5729612 Abel et al. Mar 1998 A
5732189 Johnston et al. Mar 1998 A
5749064 Pawate et al. May 1998 A
5757937 Itoh et al. May 1998 A
5774837 Yeldener et al. Jun 1998 A
5792971 Timis et al. Aug 1998 A
5796819 Romesburg Aug 1998 A
5806025 Vis et al. Sep 1998 A
5809463 Gupta et al. Sep 1998 A
5819215 Dobson et al. Oct 1998 A
5825320 Miyamori et al. Oct 1998 A
5839101 Vahatalo et al. Nov 1998 A
5920840 Satyamurti et al. Jul 1999 A
5933495 Oh Aug 1999 A
5943429 Handel Aug 1999 A
5956674 Smyth et al. Sep 1999 A
5974380 Smyth et al. Oct 1999 A
5978824 Ikeda Nov 1999 A
5983139 Zierhofer Nov 1999 A
5990405 Auten et al. Nov 1999 A
6002776 Bhadkamkar et al. Dec 1999 A
6061456 Andrea et al. May 2000 A
6072881 Linder Jun 2000 A
6097820 Turner Aug 2000 A
6108626 Cellario et al. Aug 2000 A
6122610 Isabelle Sep 2000 A
6134524 Peters et al. Oct 2000 A
6137349 Menkhoff et al. Oct 2000 A
6140809 Doi Oct 2000 A
6173255 Wilson et al. Jan 2001 B1
6180273 Okamoto Jan 2001 B1
6205421 Morii Mar 2001 B1
6216103 Wu et al. Apr 2001 B1
6222927 Feng et al. Apr 2001 B1
6223090 Brungart Apr 2001 B1
6226616 You et al. May 2001 B1
6263307 Arslan et al. Jul 2001 B1
6266633 Higgins et al. Jul 2001 B1
6317501 Matsuo Nov 2001 B1
6339758 Kanazawa et al. Jan 2002 B1
6355869 Mitton Mar 2002 B1
6363345 Marash et al. Mar 2002 B1
6381570 Li et al. Apr 2002 B2
6430295 Handel et al. Aug 2002 B1
6434417 Lovett Aug 2002 B1
6449586 Hoshuyama Sep 2002 B1
6469732 Chang et al. Oct 2002 B1
6487257 Gustafsson et al. Nov 2002 B1
6496795 Malvar Dec 2002 B1
6513004 Rigazio et al. Jan 2003 B1
6516066 Hayashi Feb 2003 B2
6529606 Jackson, Jr. II et al. Mar 2003 B1
6549630 Bobisuthi Apr 2003 B1
6584203 Elko et al. Jun 2003 B2
6622030 Romesburg et al. Sep 2003 B1
6717991 Gustafsson et al. Apr 2004 B1
6718309 Selly Apr 2004 B1
6738482 Jaber May 2004 B1
6760450 Matsuo Jul 2004 B2
6785381 Gartner et al. Aug 2004 B2
6792118 Watts Sep 2004 B2
6795558 Matsuo Sep 2004 B2
6798886 Smith et al. Sep 2004 B1
6810273 Mattila et al. Oct 2004 B1
6882736 Dickel et al. Apr 2005 B2
6915264 Baumgarte Jul 2005 B2
6917688 Yu et al. Jul 2005 B2
6944510 Ballesty et al. Sep 2005 B1
6978159 Feng et al. Dec 2005 B2
6982377 Sakurai et al. Jan 2006 B2
6999582 Popovic et al. Feb 2006 B1
7016507 Brennan Mar 2006 B1
7020605 Gao Mar 2006 B2
RE39080 Johnston Apr 2006 E
7031478 Belt et al. Apr 2006 B2
7054452 Ukita May 2006 B2
7058572 Nemer Jun 2006 B1
7065485 Chong-White et al. Jun 2006 B1
7065486 Thyssen Jun 2006 B1
7076315 Watts Jul 2006 B1
7092529 Yu et al. Aug 2006 B2
7092882 Arrowood et al. Aug 2006 B2
7099821 Visser et al. Aug 2006 B2
7142677 Gonopolskiy et al. Nov 2006 B2
7146013 Saito et al. Dec 2006 B1
7146316 Alves Dec 2006 B2
7155019 Hou Dec 2006 B2
7164620 Hoshuyama Jan 2007 B2
7171008 Elko Jan 2007 B2
7171246 Mattila et al. Jan 2007 B2
7174022 Zhang et al. Feb 2007 B1
7206418 Yang et al. Apr 2007 B2
7209567 Kozel et al. Apr 2007 B1
7225001 Eriksson et al. May 2007 B1
7242762 He et al. Jul 2007 B2
7246058 Burnett Jul 2007 B2
7254242 Ise et al. Aug 2007 B2
7254535 Kushner et al. Aug 2007 B2
7359520 Brennan et al. Apr 2008 B2
7412379 Taori et al. Aug 2008 B2
7433907 Nagai et al. Oct 2008 B2
7516067 Seltzer et al. Apr 2009 B2
7555434 Nomura et al. Jun 2009 B2
7574352 Quatieri, Jr. Aug 2009 B2
7925502 Droppo et al. Apr 2011 B2
7949522 Hetherington et al. May 2011 B2
8175291 Chan et al. May 2012 B2
8213597 Hjelm Jul 2012 B2
8705759 Wolff et al. Apr 2014 B2
8718290 Murgia et al. May 2014 B2
8744844 Klein Jun 2014 B2
8774423 Solbach Jul 2014 B1
20010016020 Gustafsson et al. Aug 2001 A1
20010031053 Feng et al. Oct 2001 A1
20020002455 Accardi et al. Jan 2002 A1
20020009203 Erten Jan 2002 A1
20020041693 Matsuo Apr 2002 A1
20020080980 Matsuo Jun 2002 A1
20020106092 Matsuo Aug 2002 A1
20020116187 Erten Aug 2002 A1
20020133334 Coorman et al. Sep 2002 A1
20020147595 Baumgarte Oct 2002 A1
20020184013 Walker Dec 2002 A1
20030014248 Vetter Jan 2003 A1
20030026437 Janse et al. Feb 2003 A1
20030033140 Taori et al. Feb 2003 A1
20030039369 Bullen Feb 2003 A1
20030040908 Yang et al. Feb 2003 A1
20030061032 Gonopolskiy Mar 2003 A1
20030063759 Brennan et al. Apr 2003 A1
20030072382 Raleigh et al. Apr 2003 A1
20030072460 Gonopolskiy et al. Apr 2003 A1
20030095667 Watts May 2003 A1
20030099345 Gartner et al. May 2003 A1
20030101048 Liu May 2003 A1
20030103632 Goubran et al. Jun 2003 A1
20030128851 Furuta Jul 2003 A1
20030138116 Jones et al. Jul 2003 A1
20030147538 Elko Aug 2003 A1
20030169891 Ryan et al. Sep 2003 A1
20030228023 Burnett et al. Dec 2003 A1
20040013276 Ellis et al. Jan 2004 A1
20040047464 Yu et al. Mar 2004 A1
20040057574 Faller Mar 2004 A1
20040078199 Kremer et al. Apr 2004 A1
20040102967 Furuta et al. May 2004 A1
20040131178 Shahaf et al. Jul 2004 A1
20040133421 Burnett et al. Jul 2004 A1
20040165736 Hetherington et al. Aug 2004 A1
20040196989 Friedman et al. Oct 2004 A1
20040263636 Cutler et al. Dec 2004 A1
20050025263 Wu Feb 2005 A1
20050027520 Mattila et al. Feb 2005 A1
20050049864 Kaltenmeier et al. Mar 2005 A1
20050060142 Visser et al. Mar 2005 A1
20050114123 Lukac et al. May 2005 A1
20050152559 Gierl et al. Jul 2005 A1
20050152563 Amada et al. Jul 2005 A1
20050185813 Sinclair et al. Aug 2005 A1
20050213778 Buck et al. Sep 2005 A1
20050216259 Watts Sep 2005 A1
20050228518 Watts Oct 2005 A1
20050240399 Makinen Oct 2005 A1
20050276423 Aubauer et al. Dec 2005 A1
20050278171 Suppappola et al. Dec 2005 A1
20050288923 Kok Dec 2005 A1
20060072768 Schwartz et al. Apr 2006 A1
20060074646 Alves et al. Apr 2006 A1
20060098809 Nongpiur et al. May 2006 A1
20060120537 Burnett et al. Jun 2006 A1
20060133621 Chen et al. Jun 2006 A1
20060149535 Choi et al. Jul 2006 A1
20060184363 McCree et al. Aug 2006 A1
20060198542 Benjelloun Touimi et al. Sep 2006 A1
20060222184 Buck et al. Oct 2006 A1
20070021958 Visser et al. Jan 2007 A1
20070027685 Arakawa et al. Feb 2007 A1
20070033020 (Kelleher) Francois et al. Feb 2007 A1
20070067166 Pan et al. Mar 2007 A1
20070078649 Hetherington et al. Apr 2007 A1
20070094031 Chen Apr 2007 A1
20070100612 Ekstrand et al. May 2007 A1
20070116300 Chen May 2007 A1
20070150268 Acero et al. Jun 2007 A1
20070154031 Avendano et al. Jul 2007 A1
20070165879 Deng et al. Jul 2007 A1
20070195968 Jaber Aug 2007 A1
20070230712 Belt et al. Oct 2007 A1
20070276656 Solbach et al. Nov 2007 A1
20080019548 Avendano Jan 2008 A1
20080033723 Jang et al. Feb 2008 A1
20080140391 Yen et al. Jun 2008 A1
20080201138 Visser et al. Aug 2008 A1
20080228474 Huang et al. Sep 2008 A1
20080228478 Hetherington et al. Sep 2008 A1
20080260175 Elko Oct 2008 A1
20090012783 Klein Jan 2009 A1
20090012786 Zhang et al. Jan 2009 A1
20090089054 Wang et al. Apr 2009 A1
20090129610 Kim et al. May 2009 A1
20090220107 Every et al. Sep 2009 A1
20090238373 Klein Sep 2009 A1
20090253418 Makinen Oct 2009 A1
20090271187 Yen et al. Oct 2009 A1
20100036659 Haulick et al. Feb 2010 A1
20100094622 Cardillo et al. Apr 2010 A1
20100094643 Avendano et al. Apr 2010 A1
20100278352 Petit et al. Nov 2010 A1
20110178800 Watts Jul 2011 A1
20110286605 Furuta et al. Nov 2011 A1
20110305345 Bouchard et al. Dec 2011 A1
20130034243 Yermeche et al. Feb 2013 A1
Foreign Referenced Citations (28)
Number Date Country
62110349 May 1987 JP
04184400 Jul 1992 JP
5053587 Mar 1993 JP
05-172865 Jul 1993 JP
06269083 Sep 1994 JP
H07248793 Sep 1995 JP
10-313497 Nov 1998 JP
11-249693 Sep 1999 JP
2004053895 Feb 2004 JP
2004531767 Oct 2004 JP
2004533155 Oct 2004 JP
2005110127 Apr 2005 JP
2005148274 Jun 2005 JP
2005518118 Jun 2005 JP
2005195955 Jul 2005 JP
2007006525 Jan 2007 JP
526468 Apr 2003 TW
I279776 Apr 2007 TW
0174118 Oct 2001 WO
02080362 Oct 2002 WO
02103676 Dec 2002 WO
03043374 May 2003 WO
03069499 Aug 2003 WO
03069499 Aug 2003 WO
2004010415 Jan 2004 WO
2007081916 Jul 2007 WO
2007140003 Dec 2007 WO
2010005493 Jan 2010 WO
Non-Patent Literature Citations (73)
Entry
Boll, Steven F. “Suppression of Acoustic Noise in Speech using Spectral Subtraction”, IEEE Transactions on Acoustics, Speech and Signal Processing, vol. ASSP-27, No. 2, Apr. 1979, pp. 113-120.
Dahl, Mattias et al., “Simultaneous Echo Cancellation and Car Noise Suppression Employing a Microphone Array”, 1997 IEEE International Conference on Acoustics, Speech, and Signal Processing, Apr. 21-24, pp. 239-242.
“Ent 172.” Instructional Module. Prince George's Community College Department of Engineering Technology. Accessed: Oct. 15, 2011. Subsection: “Polar and Rectangular Notation”. <http://academic.ppgcc.edu/ent/ent172—instr—mod.html>.
Fulghum, D. P. et al., “LPC Voice Digitizer with Background Noise Suppression”, 1979 IEEE International Conference on Acoustics, Speech, and Signal Processing, pp. 220-223.
Graupe, Daniel et al., “Blind Adaptive Filtering of Speech from Noise of Unknown Spectrum Using a Virtual Feedback Configuration”, IEEE Transactions on Speech and Audio Processing, Mar. 2000, vol. 8, No. 2, pp. 146-158.
Haykin, Simon et al. “Appendix A.2 Complex Numbers.” Signals and Systems. 2nd Ed. 2003. p. 764.
Hermansky, Hynek “Should Recognizers Have Ears?”, In Proc. ESCA Tutorial and Research Workshop on Robust Speech Recognition for Unknown Communication Channels, pp. 1-10, France 1997.
Hohmann, V. “Frequency Analysis and Synthesis Using a Gammatone Filterbank”, ACTA Acustica United with Acustica, 2002, vol. 88, pp. 433-442.
Jeffress, Lloyd A. et al. “A Place Theory of Sound Localization,” Journal of Comparative and Physiological Psychology, 1948, vol. 41, p. 35-39.
Jeong, Hyuk et al., “Implementation of a New Algorithm Using the STFT with Variable Frequency Resolution for the Time-Frequency Auditory Model”, J. Audio Eng. Soc., Apr. 1999, vol. 47, No. 4., pp. 240-251.
Kates, James M. “A Time-Domain Digital Cochlear Model”, IEEE Transactions on Signal Processing, Dec. 1991, vol. 39, No. 12, pp. 2573-2592.
Lazzaro, John et al., “A Silicon Model of Auditory Localization,” Neural Computation Spring 1989, vol. 1, pp. 47-57, Massachusetts Institute of Technology.
Lippmann, Richard P. “Speech Recognition by Machines and Humans”, Speech Communication, Jul. 1997, vol. 22, No. 1, pp. 1-15.
Martin, Rainer “Spectral Subtraction Based on Minimum Statistics”, in Proceedings Europe. Signal Processing Conf., 1994, pp. 1182-1185.
Mitra, Sanjit K. Digital Signal Processing: a Computer-based Approach. 2nd Ed. 2001. pp. 131-133.
Watts, Lloyd Narrative of Prior Disclosure of Audio Display on Feb. 15, 2000 and May 31, 2000.
Cosi, Piero et al. (1996), “Lyon's Auditory Model Inversion: a Tool for Sound Separation and Speech Enhancement,” Proceedings of ESCA Workshop on ‘The Auditory Basis of Speech Perception,’ Keele University, Keele (UK), Jul. 15-19, 1996, pp. 194-197.
Rabiner, Lawrence R. et al. “Digital Processing of Speech Signals”, (Prentice-Hall Series in Signal Processing). Upper Saddle River, NJ: Prentice Hall, 1978.
Weiss, Ron et al., “Estimating Single-Channel Source Separation Masks: Revelance Vector Machine Classifiers vs. Pitch-Based Masking”, Workshop on Statistical and Perceptual Audio Processing, 2006.
Schimmel, Steven et al., “Coherent Envelope Detection for Modulation Filtering of Speech,” 2005 IEEE International Conference on Acoustics, Speech, and Signal Processing, vol. 1, No. 7, pp. 221-224.
Slaney, Malcom, “Lyon's Cochlear Model”, Advanced Technology Group, Apple Technical Report #13, Apple Computer, Inc., 1988, pp. 1-79.
Slaney, Malcom, et al. “Auditory Model Inversion for Sound Separation,” 1994 IEEE International Conference on Acoustics, Speech and Signal Processing, Apr. 19-22, vol. 2, pp. 77-80.
Slaney, Malcom. “An Introduction to Auditory Model Inversion”, Interval Technical Report IRC 1994-014, http://coweb.ecn.purdue.edu/˜maclom/interval/1994-014/, Sep. 1994, accessed on Jul. 6, 2010.
Solbach, Ludger “An Architecture for Robust Partial Tracking and Onset Localization in Single Channel Audio Signal Mixes”, Technical University Hamburg-Harburg, 1998.
Syntrillium Software Corporation, “Cool Edit User's Manual”, 1996, pp. 1-74.
Tchorz, Jurgen et al., “SNR Estimation Based on Amplitude Modulation Analysis with Applications to Noise Suppression”, IEEE Transactions on Speech and Audio Processing, vol. 11, No. 3, May 2003, pp. 184-192.
Watts, Lloyd, “Robust Hearing Systems for Intelligent Machines,” Applied Neurosystems Corporation, 2001, pp. 1-5.
Yoo, Heejong et al., “Continuous-Time Audio Noise Suppression and Real-Time Implementation”, 2002 IEEE International Conference on Acoustics, Speech, and Signal Processing, May 13-17, pp. IV3980-IV3983.
International Search Report dated Jun. 8, 2001 in Application No. PCT/US01/08372.
International Search Report dated Apr. 3, 2003 in Application No. PCT/US02/36946.
International Search Report dated May 29, 2003 in Application No. PCT/US03/04124.
International Search Report and Written Opinion dated Sep. 16, 2008 in Application No. PCT/US07/12628.
International Search Report and Written Opinion dated May 11, 2009 in Application No. PCT/US09/01667.
International Search Report and Written Opinion dated May 20, 2010 in Application No. PCT/US09/06754.
Fast Cochlea Transform, US Trademark Reg. No. 2,875,755 (Aug. 17, 2004).
Dahl, Mattias et al., “Acoustic Echo and Noise Cancelling Using Microphone Arrays”, International Symposium on Signal Processing and its Applications, ISSPA, Gold coast, Australia, Aug. 25-30, 1996, pp. 379-382.
Demol, M. et al. “Efficient Non-Uniform Time-Scaling of Speech With WSOLA for Call Applications”, Proceedings of InSTIL/ICALL2004—NLP and Speech Technologies in Advanced Language Learning Systems—Venice Jun. 17-19, 2004.
Laroche, Jean. “Time and Pitch Scale Modification of Audio Signals”, in “Applications of Digital Signal Processing to Audio and Acoustics”, The Kluwer International Series in Engineering and Computer Science, vol. 437, pp. 279-309, 2002.
Moulines, Eric et al., “Non-Parametric Techniques for Pitch-Scale and Time-Scale Modification of Speech”, Speech Communication, vol. 16, pp. 175-205, 1995.
Verhelst, Werner, “Overlap-Add Methods for Time-Scaling of Speech”, Speech Communication vol. 30, pp. 207-221, 2000.
Avendano, C., “Frequency-Domain Techniques for Source Identification and Manipulation in Stereo Mixes for Enhancement, Suppression and Re-Panning Applications,” in Proc. IEEE Workshop on Application of Signal Processing to Audio and Acoustics, Waspaa, 03, New Paltz, NY, 2003.
Elko, Gary W., “Differential Microphone Arrays,”Audio Signal Processing for Next-Generation Multimedia Communication Systems, 2004, pp. 12-65, Kluwer Academic Publishers, Norwell, Massachusetts, USA.
B. Widrow et al., “Adaptive Antenna Systems,” Proceedings IEEE, vol. 55, No. 12, pp. 2143-2159, Dec. 1967.
Allen, Jont B. “Short Term Spectral Analysis, and Modification by Discrete Fourier Transform”, IEEE Transactions on Acoustics, Speech, and Signal Processing. vol. ASSP-25, 3. Jun. 1977. pp. 235-238.
Allen, Jont B. et al. “A Unified Approach to Short-Time Fourier Analysis and Synthesis”, Proceedings of the IEEE. vol. 65, 11, Nov. 1977. pp. 1558-1564.
Boll, Steven F. “Suppression of Acoustic Noise in Speech Using Spectral Subtraction”, Dept. of Computer Science, University of Utah Salt Lake City, Utah, Apr. 1979, pp. 18-19.
Boll, Steven et al. “Suppression of Acoustic Noise in Speech Using Two Microphone Adaptive Noise Cancellation”, source(s): IEEE Transactions on Acoustic, Speech, and Signal Processing. vol. v ASSSP-28, n 6, Dec. 1980, pp. 752-753.
Chen, Jingdong et al. “New Insights into the Noise Reduction Wierner Filter”, source(s): IEEE Transactions on Audio, Speech, and Language Processing. vol. 14, 4, Jul. 2006, pp. 1218-1234.
Cohen, Isreal, “Mutichannel Post-Filtering in Nonstationary Noise Environment”, source(s): IEEE Transactions on Signal Processing. vol. 52, 5, May 2004, pp. 1149-1160.
Cohen et al. “Microphone Array Post-Filtering for Non-Stationary Noise”, source(s): IEEE, May 2002.
Fuchs, Martin et al. “Noise Suppression for Automotive Applications Based on Directional Information”, source(s): 2004 IEEE. pp. 237-240.
Goubran, R.A. . “Acoustic Noise Suppression Using Regression Adaptive Filtering”, source(s): 1990 IEEE. pp. 48-53.
Liu, Chen et al. “A two-microphone dual delay-line approach for extraction of a speech sound in the pressence of multiple interferers”, source(s): Acoustical Society of America. vol. 110, 6, Dec. 2001, pp. 3218-3231.
Martin, Rainer et al. “Combined Acoustic Echo Cancellation, Derverberation and Noise Reduction: A two Microphone Approach”, source(s): Annles des Telecommunications of Telecommunications. vol. 29, 7-8, Jul.-Aug. 1994, pp. 429-438.
Mizumachi, Mitsunori et al. “Noise Reduction by Paired-Microphones Using Spectral Subtraction”, source(s): 1998 IEEE. pp. 1001-1004.
Moonen, Marc et at. “Multi-Microphone Signal Enhancement Techniques for Noise Suppression and Dereverbration,” source(s): http://www.esat.kuleuven.ac.be/sista/yearreport97/node37.html.
Parra, Lucas et al. “Convolutive blind Separation of Non-Stationary”, source(s): IEEE Transactions on Speech and Audio Processing. vol. 8, 3, May 2008, pp. 320-327.
Tashev, Ivan et al. “Microphone Array of Headset with Spatial Noise Suppressor”, source(s): http://research.microsoft.com/users/ivantash/Documents/Tashev—MAforHeadset—HSCMA—05.pdf. (4 pages).
Valin, Jean-Marc et al. “Enhanced Robot Audition Based on Micophone Array Source Separation with Post-Filter”, source(s): Proceedings of 2004 IEEE/RSJ International Conference on Intelligent Robots and Systems, Sep. 28-Oct. 2, 2004, Sendai, Japan. pp. 2123-2128.
Stahl, V.; Fischer, A.; Bippus, R.; “Quantile based noise estimation for spectral subtraction and Wiener filtering,” Acoustics, Speech, and Signal Processing, 2000. ICASSP '00. Proceedings. 2000 IEEE International Conference on, vol. 3, no., pp. 1875-1878 vol. 3, 2000.
International Search Report and Written Opinion dated Aug. 27, 2009 in Application No. PCT/US09/03813.
International Search Report and Written Opinion dated Oct. 19, 2007 in Application No. PCT/US07/00463.
International Search Report and Written Opinion dated Oct. 1, 2008 in Application No. PCT/US08/08249.
International Search Report and Written Opinion dated Apr. 9, 2008 in Application No. PCT/US07/21654.
Notice of Allowance, Jun. 5, 2014, U.S. Appl. No. 12/228,034, filed Aug. 8, 2008.
Office Action, May 13, 2014, U.S. Appl. No. 12/962,519, filed Dec. 7, 2010.
Office Action, Jul. 15, 2014, U.S. Appl. No. 13/432,490, filed Mar. 28, 2012.
Notice of Allowance, Jul. 16, 2014, U.S. Appl. No. 13/426,436, filed Mar. 21, 2012.
Notice of Allowance, Jun. 19, 2014, U.S. Appl. No. 13/705,132, filed Dec. 4, 2012.
Allowance mailed May 21, 2014 in Finnish Patent Application 20100001, filed Jan. 4, 2010.
Office Action mailed May 2, 2014 in Taiwanese Patent Application 098121933, filed Jun. 29, 2009.
Office Action mailed Jun. 27, 2014 in Korean Patent Application No. 10-2010-7000194, filed Jan. 6, 2010.
Office Action mailed Jun. 18, 2014 in Finnish Patent Application No. 20080428, filed Jul. 4, 2008.
Related Publications (1)
Number Date Country
20090323982 A1 Dec 2009 US