Ambient noise root mean square (RMS) detector

Information

  • Patent Grant
  • 9107010
  • Patent Number
    9,107,010
  • Date Filed
    Friday, February 8, 2013
    11 years ago
  • Date Issued
    Tuesday, August 11, 2015
    8 years ago
Abstract
An RMS detector uses the concept of the k-NN (classifying using nearest neighbors)—algorithm in order to obtain RMS values. A rms detector using first-order regressor with a variable smoothing factor is modified to penalize samples from center of data in order to obtain RMS values. Samples which vary greatly from the background noise levels, such as speech, scratch, wind and other noise spikes, are dampened in the RMS calculation. When background noise changes, the system will track the changes in background noise and include the changes in the calculation of the corrected RMS value. A minimum tracker runs more often (e.g. two or three times) than the rate as in prior art detectors and methods, tracks the minimum rms value, which is to compute a normalized distance value, which in turn is used to normalize the smoothing factor. From this data, a corrected or revised RMS value is determined as the function of the previous RMS value multiplied by one minus the smoothing factor plus the smooth factor times the minimum rms value to output the corrected RMS for the present invention. The rms value is used to generate a reset signal for the minimum tracker and is used to avoid deadlock in the tracker, for example, when the background signal increases/decreases over time.
Description
FIELD OF THF INVENTION

The present invention relates to an ambient noise Root Mean Square (RMS) level detector. In particular, the present invention is directed toward an improved noise RMS detector that is robust to speech presence, wind noise, and other sudden variations in noise levels.


BACKGROUND OF THF INVENTION

A personal audio device, such as a wireless telephone, includes an adaptive noise canceling (ANC) circuit that adaptively generates an anti-noise signal from a reference microphone signal and injects the anti-noise signal into the speaker or other transducer output to cause cancellation of ambient audio sounds. An error microphone is also provided proximate the speaker to measure the ambient sounds and transducer output near the transducer, thus providing an indication of the effectiveness of the noise canceling. A processing circuit uses the reference and/or error microphone, optionally along with a microphone provided for capturing near-end speech, to determine whether the ANC circuit is incorrectly adapting or may incorrectly adapt to the instant acoustic environment and/or whether the anti-noise signal may be incorrect and/or disruptive and then takes actions in the processing circuit to prevent or remedy such conditions.


Examples of such adaptive noise cancellation systems are disclosed in published U.S. Patent Application 2012/0140943, published on Jun. 7, 2012, and in published U.S. Patent Application 2012/0207317, published on Aug. 16, 2012, both of which are incorporated herein by reference. Both of these references are assigned to the same assignee as the present application and name at least one inventor in common and thus are not prior art to the present application, but are provided to facilitate the understating of ANC circuits as applied in the field of use.


Referring now to FIG. 1, a wireless telephone 10 is illustrated in accordance with an embodiment of the present invention is shown in proximity to a human ear 5. Wireless telephone 10 includes a transducer, such as speaker SPKR that reproduces distant speech received by wireless telephone 10, along with other local audio events such as ring tones, stored audio program material, injection of near-end speech (i.e., the speech of the user of wireless telephone 10) to provide a balanced conversational perception, and other audio that requires reproduction by wireless telephone 10, such as sources from web-pages or other network communications received by wireless telephone 10 and audio indications such as battery low and other system event notifications. A near-speech microphone NS is provided to capture near-end speech, which is transmitted from wireless telephone 10 to the other conversation participant(s).


Wireless telephone 10 includes adaptive noise canceling (ANC) circuits and features that inject an anti-noise signal into speaker SPKR to improve intelligibility of the distant speech and other audio reproduced by speaker SPKR. A reference microphone R is provided for measuring the ambient acoustic environment and is positioned away from the typical position of a user's/talker's mouth, so that the near-end speech is minimized in the signal produced by reference microphone R. A third microphone, error microphone E, is provided in order to further improve the ANC operation by providing a measure of the ambient audio combined with the audio reproduced by speaker SPKR close to ear 5, when wireless telephone 10 is in close proximity to ear 5. Exemplary circuit 14 within wireless telephone 10 includes an audio CODEC integrated circuit 20 that receives the signals from reference microphone R, near speech microphone NS, and error microphone E and interfaces with other integrated circuits such as an RF integrated circuit 12 containing the wireless telephone transceiver.


In general, the ANC techniques measure ambient acoustic events (as opposed to the output of speaker SPKR and/or the near-end speech) impinging on reference microphone R, and by also measuring the same ambient acoustic events impinging on error microphone E, the ANC processing circuits of illustrated wireless telephone 10 adapt an anti-noise signal generated from the output of reference microphone R to have a characteristic that minimizes the amplitude of the ambient acoustic events at error microphone E. Since acoustic path P(z) (also referred to as the Passive Forward Path) extends from reference microphone R to error microphone E, the ANC circuits are essentially estimating acoustic path P(z) combined with removing effects of an electro-acoustic path S(z) (also referred to as Secondary Path) that represents the response of the audio output circuits of CODEC IC 20 and the acoustic/electric transfer function of speaker SPKR including the coupling between speaker SPKR and error microphone E in the particular acoustic environment, which is affected by the proximity and structure of ear 5 and other physical objects and human head structures that may be in proximity to wireless telephone 10, when wireless telephone is not firmly pressed to ear 5.


Such adaptive noise cancellation (ANC) systems may employ a Root Mean Square (rms) detector to detect average background noise levels. Such an RMS detector needs to track background noise levels slowly but not so slowly as to become insensitive to environmental variations. An ideal RMS detector should be robust to speech presence, robust to scratching (contact) on the microphone, robust to wind noise, and a have a low computational complexity. For the purposes of describing the present ambient noise RMS detector, the lower case rms variable is utilized to refer to the prior art techniques and the upper case RMS to represent the corrected signal of the present ambient noise RMS detector, as set forth below. The present ambient noise RMS detector may utilize the prior art rms value in generating the RMS signal.


Perhaps the most well-known background noise estimation method, based on minimum statistics, was the rms detector introduced by Ranier Martin. See, Martin, Ranier, Noise Power Spectral Density Estimation Based on Optimal Smoothing and Minimum Statistics, IEEE Transactions on Speech and Audio Processing, Col. 9, No, 5, July 2001, incorporated herein by reference, as well as Martin, Ranier, Spectral Subtraction Based on Minimum Statistics, in Proc. 7th EUSIPCO '94, Edinburgh, U.K., Sep. 13-16, 1994, pp/. 1182-1195, also incorporated herein by reference. Israel Cohen has made another RMS detector based on the Martin design. See, Cohen, Israel, Noise Spectrum Estimation in Adverse Environments: Improved Minima Controlled Recursive Averaging, IEEE Transactions on Speech and Audio Processing, Vol. 11, Issue 5, September 2003, incorporated herein by reference as well as Cohen, Israel, Noise Estimation by Minima Controlled Recursive Averaging for Robust Speech Enhancement, IEEE Signal Processing Letters, Vol. 9, No. 1, January 2002, also incorporated herein by reference. Both the Martin and Cohen methods and designs employ a method to track the minimum RMS value. Both methods also use a first-order regressor with a variable smoothing factor.


The Cohen design may be less complex compared and provides better performance compared to the Martin design. The Cohen design depends on a couple of thresholds and parameters that should be adjusted for different applications. The Cohen design also uses less memory than the Martin design in that previous values of rms are kept to find the minimum value. The problem with the Cohen design is that it is susceptible to non-stationary noise such as spike noise. For example, when used in an adaptive noise cancellation system (ANC) on a cellular phone or the like, spike noise such as wind noise or scratching (user's/talker's hand scratching or rubbing the case) may create spikes to which the Cohen design would over-react. As a result, the performance of an ANC system, for example, in a cellular telephone or the like, may be degraded, as the rms detector over-reacts to these spike noises.


A simple rms detector based on a first order regression may produce an output illustrated in FIG. 2. This first order regression may be calculated as shown in equation (1):











rms


(
n
)


=



(

1
-
α

)

·

rms


(

n
-
1

)



+

α
·



input


(
n
)













α
=

{




α
att






input


>

rms


(

n
-
1

)








α
dec



else









(
1
)








where α represents a smoothing factor, rms(n) represents the rms value for the sample n and input(n) represents the input signal for sample n, and n is a sample integer number. Thus, the rms value in equation (1) is calculated by multiplying a smoothing factor (subtracted from one) times the previous rms value and then adding the absolute value of the input value times this same smoothing factor. The smoothing factor α may be selected from one of two values, αatt or αdec depending on whether the absolute value of the input signal is greater or less than the previous rms value.


The problem with such a simple rms detector is that it not only tracks background noise, but also speech, scratch, and wind noise. As illustrated in FIG. 2, the outer darker line 210 represents a speech signal, with occasional spike noise 220 as shown. The lighter line 230 represents the rms signal, calculated with a slow attack and fast decay, as shown in Equation (1). As can be seen in FIG. 2, the rms value 230 calculated using Equation (1) ends up tracking these spike signals 220, which maybe undesirable for an adaptive noise cancellation (ANC) circuit. By tracking the spike signals 220, the ANC circuit may end up generating inappropriate anti-noise, and as a result, create artifacts in the reproduced audio signal for the user.


SUMMARY OF THE INVENTION

The present ambient noise RMS detector represents an improvement over the prior art rms detector from a adaptive or machine learning perspective. The present ambient noise RMS detector uses the concept of a k-NN (classifying using nearest neighbors) algorithm in order to obtain RMS values. The k-nearest neighbor algorithm (k-NN) is a method for classifying objects based on closest training examples in the feature space. k-NN is a type of instance-based learning, or lazy learning where the function is only approximated locally and all computation is deferred until classification. An object is classified by a majority vote of its neighbors, with the object being assigned to the class most common amongst its k nearest neighbors (k is a positive integer, typically small). If k=1, then the object is simply assigned to the class of its nearest neighbor.


The same method can be used for regression, by simply assigning the property value for the object to be the average of the values of its k nearest neighbors. It can be useful to weight the contributions of the neighbors, so that the nearer neighbors contribute more to the average than the more distant ones. (A common weighting scheme is to give each neighbor a weight of 1/d, where d is the distance to the neighbor. This scheme is a generalization of linear interpolation.)


The present invention incorporates a prior art rms detector using first-order regressor with a variable smoothing factor but adds additional features to penalize samples from center of data in order to obtain RMS values. Thus, samples which vary greatly from the background noise levels, such as speech, scratch, and other noise spikes, are dampened in the RMS calculation. However, when background noise increases/decreases (changes in general), the system will track this change in background noise and include that in the calculation of the corrected RMS value.


Output from a prior art rms detector using a first-order regressor with a variable smoothing factor is fed to a minimum tracker, which is also known in the art. The minimum tracker tracks the minimum rms value, Rmin over time. This revised minimum value is used to compute a normalized distance value d, which represents the ratio expressed as the absolute value of the difference between the previously calculated rms value, and the RMS value calculated in the present ambient noise RMS detector divided by the RMS value calculated by the present ambient noise RMS detector. This value d in turn is used to normalize the smoothing factor α by dividing the smoothing factor by the maximum of d or 1.


Once these values are calculated, a corrected or revised RMS value can be determined as the function of the previous RMS value multiplied by one minus the smoothing factor plus the smooth factor times the minimum rms value to output the corrected RMS for the present ambient noise RMS detector. The rms value may be used to generate a reset signal for the minimum tracker. This reset signal may be operated on an order of 0.1 to 1 seconds and is used to avoid deadlock in the tracker, for example, when the background signal increases over time.


The effect of the present ambient noise RMS detector, as demonstrated in the Figures attached herewith, is to provide a background RMS value which is largely immune from sudden spikes in value, such as due to speech, “scratching” (when a person physically touches the microphone, for example), or wind noise, particularly when compared to the prior art techniques.


While discussed herein in the context of cellular telephones and adaptive noise cancellation circuits used therein, the present ambient noise RMS detector has applications for a number of audio devices and the like. For example, the RMS detector of the present invention may be applied to audio and audio-visual recording equipment, computing devices equipped with microphones, speech recognition systems, speech activated systems (e.g., in automobiles), and even event detectors, such as alarm systems, where it may be desirable to filter background sounds from sudden noises, such as glass-break or speech by intruders. While disclosed in the context of cellular phones and adaptive noise cancellation circuits, the present ambient noise RMS detector should in no way be construed as being limited to that particular application.





BRIEF DESCRIPTION OF THE DRAWINGS


FIG. 1 is a diagram illustrating how dual microphones may be used in an adaptive noise cancellation circuit in a cellular telephone.



FIG. 2 is a graph illustrating a voice signal with spike components and the resulting rms signal calculation using the techniques of the prior art.



FIG. 3 is a block diagram of an embodiment of the present ambient noise RMS detector.



FIG. 4 is a graph illustrating how the minimum RMS value is tracked.



FIG. 5A is a graph illustrating instantaneous RMS and ambient RMS for a sample input signal comprising background noise with speech.



FIG. 5B is a graph illustrating the value a calculated from the instantaneous RMS according to equation (7) and block 160 in FIG. 3.



FIG. 5C is a graph illustrating the calculation of distance value d according to equation (6) and block 150 of FIG. 3.



FIG. 5D is a graph illustrating the value of resulting Rmin as determined from equation (2) below and block 140 of FIG. 3.



FIG. 6 is a graph comparing a signal containing background noise with speech, showing a comparison between the old method of the prior art and the technique and apparatus of the present ambient noise RMS detector.



FIG. 7 is a graph comparing a signal containing background noise with a “scratch” signal in the background noise, showing a comparison between the old method of the prior art and the technique and apparatus of the present ambient noise RMS detector.





DETAILED DESCRIPTION OF THE INVENTION

The present ambient noise RMS detector improves upon the techniques of prior art rms detectors such as taught by Martin and Cohen by using an improved algorithm in the RMS detector. FIG. 3 is a block diagram of the present ambient noise RMS detector. Referring to FIG. 3, a raw rms value is calculated from the input signal using known prior art techniques. Blocks 110, 120, and 130 are elements of a first-order regressor with a variable smoothing factor. The input signal, which in this instance may be a background noise signal with speech, is fed to block 110 where the absolute value of the signal is taken. This absolute value signal in turn is fed to low-pass-filter 120 and then to downsampler 130. The net effect is to output a raw rms value such as described above in connection with Equation (1). As these first three elements of the block diagram are known in the art, they will not be described in further detail.


Both the Martin and Cohen methods and designs discussed above also employ a method to track the minimum rms value, Rmin, and tracking the minimum rms value is one function of the present ambient noise RMS detector. Speech, scratching (physical contact) on the microphone, wind noise, and any spike noise are all unlikely background noise in that they are not always present but appear as noise spikes in the ambient noise signal. This fact can be leveraged by comparing a short-term minimum RMS value with a long-term one to determine whether such a spike has occurred. FIG. 4 is a graph illustrating how the minimum RMS value is tracked. For every instantaneous transition, short-term rms values Rmin and Rtmp may be calculated as:









{






R
min



(
l
)


=

min


{



R
min



(

l
-
1

)


,

rms


(
l
)



}










R
tmp



(
l
)


=

min


{



R
tmp



(

l
-
1

)


,

rms


(
l
)



}










(
2
)








where Rmin is the minimum rms value over time, and Rtmp is a temporary minimum rms value to track background noise changes.


The reset mechanism for the ambient noise detector is then calculated simultaneously with equation (2). This reset mechanism calculates a long-term rms value every 0.1 to 1 seconds for values Rmin and Rtmp as:









{






R
min



(
l
)


=

min


{



R
tmp



(

l
-
1

)


,

rms


(
l
)



}










R
tmp



(
l
)


=

rms


(
l
)










(
3
)







As illustrated in FIG. 4, this approach has the effect of delaying the change in minimum RMS value Rmin in response to changes in the base rms calculation of background noise rms value BK rms. As the background rms signal increases from level A to level B, the temporary minimum value Rtmp, calculated according to Equations (2) and (3) above, rises from level A to level B, delayed over time, as illustrated in FIG. 4. The value of minimum RMS value Rmin rises from level A to level B delayed even further (the same is true for decreasing from level B to level A), as illustrated in FIG. 4. Although FIG. 4 only shows the case where level A is less than level B, the same effect occurs when level A is greater than level B as well.


In Cohen's method from this minimum RMS value Rmin calculation, it may be possible to calculate RMS using a first approach based on the probability of the presence of disturbance in the background noise signal:

RMS(l)=αd(l)·RMS(l−1)+(1−αd(l))·|input|
αd(l)=β+(1−β)*p(l)
p(l)→1→αd(l)→1  (4)


Here, p(1) is the probability of the presence of any disturbance (e.g., speech presence), and as this probability approaches one, the smoothing factor value approaches one. This probability value may be calculated as follows:











p


(
l
)


=



α
p

·

p


(

l
-
1

)



+


(

1
-

α
p


)

·

I


(
l
)












I


(
l
)


=

{



1





rms


(
l
)




R
min



(
l
)



>
δ





0


else









(
5
)








where αp represents a smoothing factor, and δ is the threshold which determines the level of any disturbance compared to Rmin(1).


One problem with this RMS tracking technique is that there are too many parameters to adjust. In addition, its reaction time is slow and is not robust. Speech rms can leak to the background RMS value. While the prior art Cohen design has additional components to make the system more robust, the system still suffers from these same operational problems. Thus, the present ambient noise RMS detector improves on the algorithms of equations (4) and (5) to provide an improved minimum RMS value Rmin tracking technique and RMS calculation.


Referring back to FIG. 3, in the present ambient noise RMS detector, the output raw rms value is then fed to a minimum tracker 140. In block 150, the normalized distance d between the current RMS and the instantaneous rms value is computed as:









d
=





rms


(
l
)


-

RMS


(

l
-
1

)






RMS


(

l
-
1

)







(
6
)








where rms(l) is a raw rms value for sample l and RMS(l−1) is a previous corrected RMS value.


In block 160, the smooth factor is normalized with this distance d:











α
d



(
l
)


=


α
0


max


(

d
,
1

)







(
7
)








where αd(l) represents the normalized smoothing factor for sample l and α0 represents a standard smoothing factor, and max(d,1) is the maxima of the normalized distance and 1. The normalized smoothing factor is then fed to block 170:

RMS(l)=(1−αd(l))·RMS(l−1)+αd(lRmin(l)  (8)

where RMS(l) is the corrected RMS value, and RMS(l−1) is a previous corrected RMS value, αd(l) represents the normalized smoothing factor for sample l as calculated in equation (7) and minimum RMS value Rmin is the minimum rms value calculated in equation (3).


The raw rms value is also fed to block 190, which then generates a reset signal Reset. The reset signal Reset is triggered in order to reset the system to avoid any deadlock, for example, when the background noise signal rises gradually. The reset mechanism is shown in equation (3) as discussed previously.



FIGS. 4-6 are graphs illustrating the operation of the present ambient noise RMS detector. In FIG. 5A, the instantaneous RMS and ambient RMS are shown for a sample input signal comprising background noise with speech. In FIG. 5A, the background noise appears as the baseline signal 510 and the speech portion appears in the center as the elevated portion 520. The instantaneous rms appears as the thick line (510, 520), while the final calculated ambient RMS appears as the thin line 530 below the thick line. In FIG. 5B, the value a is shown calculated from the instantaneous rms according to equation (7) above and block 160 in FIG. 3. FIG. 5C shows the calculation of d according to equation (6) above and block 150 of FIG. 3. FIG. 5D shows the resulting minimum RMS value Rmin as determined from equation (8) above and block 170 of FIG. 3.



FIG. 6 is a graph comparing a signal containing background noise with speech, showing a comparison between the old method of the prior art and the technique and apparatus of the present invention. The rms(l) signal is shown as the wide dark signal 610 in FIG. 6 with the speech disturbance 620 in the central portion. The rms calculation using the prior art method is shown as the wavy light line 630 in the center of that signal. As shown in FIG. 6, spikes occur in this signal in relationship to the source signal. As illustrated in FIG. 6, the prior art technique is sensitive to speech in the background noise signal. The bottom line 640 represents the RMS value calculated using the technique of the present ambient noise RMS detector. As illustrated in FIG. 6, the technique of the present ambient noise RMS detector is far less responsive to transient spikes than the prior art technique.



FIG. 7 is a graph comparing a signal containing background noise 710 with a scratch signal 720 in the background noise, and showing a comparison between the old method of the prior art and the technique and apparatus of the present ambient noise RMS detector. The scratch signals 720 are more pronounced than the speech signals 620 of FIG. 6. The rms(l) signal is shown as the wide dark signal 710 in FIG. 7. The rms calculation using the prior art method is shown as the wavy light line 730 in the center of that signal. As shown in FIG. 7, spikes 720 occur in this signal in relationship to the source signal 710. The bottom line 740 represents the RMS value calculated using the technique for the present ambient noise RMS detector. As illustrated in FIG. 7, the technique for the present ambient noise RMS detector is far less responsive to transient spikes than the prior art technique.


The present ambient noise RMS detector has thus been proven to more accurately calculate RMS values from an input signal, while being relatively immune to speech, wind noise, scratch, and other signal spikes. This improved RMS value calculation provides a better input value for an adaptive noise cancellation (ANC) circuit for use, for example, in a cellular telephone or the like. This improved value in turn allows for better operation of the ANC circuit, creating fewer artifacts or dropped out audio (e.g., due to the ANC circuit overcompensating and muting desired audio signals) in the audio output to the user.


While embodiments of the present ambient noise RMS detector have been disclosed and described in detail herein, it may be apparent to those skilled in the art that various changes in form and detail may be made therein without departing from the spirit and scope thereof.

Claims
  • 1. A root mean square (RMS) detector detecting an RMS level of a background noise input signal while being substantially immune to voice, wind, scratch sounds, and any spike noise, the RMS detector comprising: a raw rms detector receiving a background noise input signal and outputting a raw rms value;a minimum rms tracker receiving the raw rms value and tracking a minimum rms value of the raw rms value;a normalized distance tracker receiving the minimum rms value and calculating a distance value between the minimum rms value and a previous corrected RMS value;a normalized smoothing factor calculator normalizing a smoothing factor by dividing the smoothing factor by a maximum of the distance value or 1; andan RMS value calculator determining a corrected RMS value from the minimum rms value, a previous corrected RMS value, and the normalized smoothing factor, and outputting a corrected RMS value.
  • 2. The RMS detector of claim 1, further comprising a reset generator receiving the raw rms value and generating a reset signal to the minimum rms tracker to reset the minimum rms tracker when the raw rms value changes in value over time to prevent the minimum rms tracker from locking up.
  • 3. The RMS detector of claim 2, wherein the raw rms detector determines raw rms by adding a previous raw rms value to an input signal value.
  • 4. The RMS detector of claim 3, wherein the absolute value of the input signal value is multiplied by a smoothing factor prior to being added to the previous raw rms value.
  • 5. The RMS detector of claim 4, wherein the previous rms value is multiplied by one minus the smoothing factor prior to being added to the input signal value.
  • 6. The RMS detector of claim 5 wherein the smoothing factor is selected from one of two predetermined values depending on whether the absolute value of the input signal is greater or less than the previous raw rms value.
  • 7. The RMS detector of claim 2, where in the raw rms detector determines raw rms by:
  • 8. The RMS detector of claim 2, wherein the minimum tracker determines a short-term minimum rms value by taking the minimum of the previous minimum rms value and the current raw rms value, and for every 0.1 to 1 seconds, calculating a long-term minimum rms value as the minimum of a previous temporary minimum rms value and the present raw rms value to reset the detector, where the temporary rms value tracks background noise changes.
  • 9. The RMS detector of claim 8, wherein the minimum tracker sets the temporary rms value to a current raw rms value and the minimum rms value to a minimum of a previous temporary rms value and the current raw rms value at every 0.1 to 1 seconds to more closely track the minimum rms value.
  • 10. The RMS detector of claim 9, wherein the normalized distance is calculated by dividing the difference between the current raw rms value and the previous corrected RMS value by the previous corrected RMS value.
  • 11. The RMS detector of claim 10, wherein the normalized smoothing factor is calculated by dividing a standard predetermined smoothing factor by the maxima of the normalized distance and one.
  • 12. The RMS detector of claim 11, wherein the corrected RMS value output by the RMS detector is calculated by the sum of the normalized smoothing factor times the minimum rms value determined by the minimum rms value tracker and the product of the previous corrected RMS value times one minus the normalized smoothing factor.
  • 13. The RMS detector of claim 2, wherein the minimum tracker determines the minimum rms value by taking the minimum of the previous minimum rms value and the current raw rms value
  • 14. The RMS detector of claim 13, wherein the normalized distance d is calculated by:
  • 15. The RMS detector of claim 14, wherein the normalized smoothing factor is calculated by:
  • 16. The RMS detector of claim 15, wherein the corrected RMS value output by the RMS detector is calculated by: RMS(l)=(1−αd(l))·RMS(l−1)+αd(l)·Rmin(l)where RMS(l) is the corrected RMS value, and RMS(l−1) is a previous corrected RMS value, αd(l) represents the normalized smoothing factor for sample l, determined by the normalized smoother factor calculator, and Rmin is the minimum rms value determined by the minimum rms value tracker.
  • 17. In an RMS detector, a method of detecting RMS level of a background noise input signal while being substantially immune to voice, scratch, wind sounds, and any spike noise, the method comprising: generating in an initial RMS detector receiving a background noise input signal, a raw rms value;tracking in a minimum rms tracker receiving the raw rms value, a minimum rms value of the raw rms value;calculating in a normalized distance tracker receiving the minimum rms value, a distance value between the minimum rms value and a previous corrected RMS value;normalizing, in a normalized smoothing factor calculator, a smoothing factor by dividing the smoothing factor by a maximum of the distance value or 1; andcalculating in an RMS value calculator, a corrected RMS value by determining a corrected RMS value from the minimum rms value, a previous corrected RMS value, and the normalized smoothing factor.
  • 18. The method of claim 17, further comprising: generating in a reset generator receiving the raw rms value, a reset signal to the minimum rms tracker to reset the minimum rms tracker when the raw rms value changes in value over time to prevent the minimum rms tracker from locking up.
  • 19. The method of claim 18, wherein the raw rms detector determines raw rms by adding a previous raw rms value to an input signal value.
  • 20. The method of claim 19, wherein the absolute value of the input signal value is multiplied by a smoothing factor prior to being added to the previous raw rms value.
  • 21. The method of claim 20, wherein the previous raw rms value is multiplied by one minus the smoothing factor prior to being added to the input signal value.
  • 22. The method of claim 21, wherein the smoothing factor is selected from one of two predetermined values depending on whether the absolute value of the input signal is greater or less than the previous raw rms value.
  • 23. The method of claim 18, where in the raw rms detector determines raw rms by:
  • 24. The method of claim 18, wherein the minimum tracker determines a short-term minimum rms value by taking the minimum of the previous minimum rms value and the current raw rms value, and for every 0.1 to 1 seconds, calculating a long-term minimum rms value as the minimum of a previous temporary minimum rms value and the present raw rms value to reset the detector, where the temporary rms value tracks background noise changes.
  • 25. The method of claim 24, wherein the minimum tracker sets the temporary rms value to a current raw rms value and the minimum rms value to a minimum of a previous temporary rms value and the current raw rms value at every 0.1 to 1 seconds to more closely track the minimum rms value.
  • 26. The method of claim 25, wherein the normalized distance is calculated by dividing the difference between the current raw rms value and the previous corrected RMS value by the previous corrected RMS value.
  • 27. The method of claim 26, wherein the normalized smoothing factor is calculated by dividing a standard predetermined smoothing factor by the maxima of the normalized distance and one.
  • 28. The method of claim 27, wherein the corrected RMS value output by the RMS detector is calculated by the sum of the normalized smoothing factor times the minimum rms value determined by the minimum rms value tracker, and the product of the previous corrected RMS value times one minus the normalized smoothing factor.
  • 29. The method of claim 18, wherein the minimum tracker determines the minimum rms value by taking the minimum of the previous minimum rms value and the current raw rms value
  • 30. The method of claim 29, wherein the normalized distance d is calculated by:
  • 31. The RMS detector of claim 30, wherein the normalized smoothing factor is calculated by:
  • 32. The RMS detector of claim 31, wherein the corrected RMS value output by the RMS detector is calculated by: RMS(l)=(1−αd(l))·RMS(l−1)+αd(l)·Rmin(l)where RMS(l) is the corrected RMS value, and RMS(l−1) is a previous corrected RMS value, αd(l) represents the normalized smoothing factor for sample l, determined by the normalized smoother factor calculator, and Rmin is the minimum rms value determined by the minimum rms value tracker.
US Referenced Citations (181)
Number Name Date Kind
5251263 Andrea et al. Oct 1993 A
5278913 Delfosse et al. Jan 1994 A
5321759 Yuan Jun 1994 A
5337365 Hamabe et al. Aug 1994 A
5359662 Yuan et al. Oct 1994 A
5410605 Sawada et al. Apr 1995 A
5425105 Lo et al. Jun 1995 A
5445517 Kondou et al. Aug 1995 A
5465413 Enge et al. Nov 1995 A
5548681 Gleaves et al. Aug 1996 A
5586190 Trantow et al. Dec 1996 A
5640450 Watanabe Jun 1997 A
5699437 Finn Dec 1997 A
5706344 Finn Jan 1998 A
5740256 Castello Da Costa et al. Apr 1998 A
5768124 Stothers et al. Jun 1998 A
5815582 Claybaugh et al. Sep 1998 A
5832095 Daniels Nov 1998 A
5946391 Dragwidge et al. Aug 1999 A
5991418 Kuo Nov 1999 A
6041126 Terai et al. Mar 2000 A
6118878 Jones Sep 2000 A
6219427 Kates et al. Apr 2001 B1
6278786 McIntosh Aug 2001 B1
6282176 Hemkumar Aug 2001 B1
6418228 Terai et al. Jul 2002 B1
6434246 Kates et al. Aug 2002 B1
6434247 Kates et al. Aug 2002 B1
6522746 Marchok et al. Feb 2003 B1
6683960 Fujii et al. Jan 2004 B1
6766292 Chandran Jul 2004 B1
6768795 Feltstrom et al. Jul 2004 B2
6850617 Weigand Feb 2005 B1
6940982 Watkins Sep 2005 B1
7103188 Jones Sep 2006 B1
7181030 Rasmussen et al. Feb 2007 B2
7330739 Somayajula Feb 2008 B2
7365669 Melanson Apr 2008 B1
7680456 Muhammad et al. Mar 2010 B2
7742790 Konchitsky et al. Jun 2010 B2
7817808 Konchitsky et al. Oct 2010 B2
8019050 Mactavish et al. Sep 2011 B2
8249262 Chua et al. Aug 2012 B2
8290537 Lee et al. Oct 2012 B2
8325934 Kuo Dec 2012 B2
8379884 Horibe et al. Feb 2013 B2
8401200 Tiscareno et al. Mar 2013 B2
8442251 Jensen et al. May 2013 B2
20010053228 Jones Dec 2001 A1
20020003887 Zhang et al. Jan 2002 A1
20030063759 Brennan et al. Apr 2003 A1
20030185403 Sibbald Oct 2003 A1
20040047464 Yu et al. Mar 2004 A1
20040165736 Hetherington et al. Aug 2004 A1
20040167777 Hetherington et al. Aug 2004 A1
20040202333 Csermak et al. Oct 2004 A1
20040264706 Ray et al. Dec 2004 A1
20050004796 Trump et al. Jan 2005 A1
20050018862 Fisher Jan 2005 A1
20050117754 Sakawaki Jun 2005 A1
20050207585 Christoph Sep 2005 A1
20050240401 Ebenezer Oct 2005 A1
20060035593 Leeds Feb 2006 A1
20060069556 Nadjar et al. Mar 2006 A1
20060153400 Fujita et al. Jul 2006 A1
20070030989 Kates Feb 2007 A1
20070033029 Sakawaki Feb 2007 A1
20070038441 Inoue et al. Feb 2007 A1
20070047742 Taenzer et al. Mar 2007 A1
20070053524 Haulick et al. Mar 2007 A1
20070076896 Hosaka et al. Apr 2007 A1
20070154031 Avendano et al. Jul 2007 A1
20070258597 Rasmussen et al. Nov 2007 A1
20070297620 Choy Dec 2007 A1
20080019548 Avendano Jan 2008 A1
20080101589 Horowitz et al. May 2008 A1
20080107281 Togami et al. May 2008 A1
20080144853 Sommerfeldt et al. Jun 2008 A1
20080177532 Greiss et al. Jul 2008 A1
20080181422 Christoph Jul 2008 A1
20080226098 Haulick et al. Sep 2008 A1
20080240455 Inoue et al. Oct 2008 A1
20080240457 Inoue et al. Oct 2008 A1
20090012783 Klein Jan 2009 A1
20090034748 Sibbald Feb 2009 A1
20090041260 Jorgensen et al. Feb 2009 A1
20090046867 Clemow Feb 2009 A1
20090060222 Jeong et al. Mar 2009 A1
20090080670 Solbeck et al. Mar 2009 A1
20090086990 Christoph Apr 2009 A1
20090175466 Elko et al. Jul 2009 A1
20090196429 Ramakrishnan et al. Aug 2009 A1
20090220107 Every et al. Sep 2009 A1
20090238369 Ramakrishnan et al. Sep 2009 A1
20090245529 Asada et al. Oct 2009 A1
20090254340 Sun et al. Oct 2009 A1
20090290718 Kahn et al. Nov 2009 A1
20090296965 Kojima Dec 2009 A1
20090304200 Kim et al. Dec 2009 A1
20090311979 Husted et al. Dec 2009 A1
20100014683 Maeda et al. Jan 2010 A1
20100014685 Wurm Jan 2010 A1
20100061564 Clemow et al. Mar 2010 A1
20100069114 Lee et al. Mar 2010 A1
20100082339 Konchitsky et al. Apr 2010 A1
20100098263 Pan et al. Apr 2010 A1
20100124335 Wessling et al. May 2010 A1
20100124336 Shridhar et al. May 2010 A1
20100124337 Wertz et al. May 2010 A1
20100150367 Mizuno Jun 2010 A1
20100158330 Guissin et al. Jun 2010 A1
20100166203 Peissig et al. Jul 2010 A1
20100195838 Bright Aug 2010 A1
20100195844 Christoph et al. Aug 2010 A1
20100246855 Chen Sep 2010 A1
20100266137 Sibbald et al. Oct 2010 A1
20100272276 Carreras et al. Oct 2010 A1
20100272283 Carreras et al. Oct 2010 A1
20100274564 Bakalos et al. Oct 2010 A1
20100284546 De Brunner et al. Nov 2010 A1
20100291891 Ridgers et al. Nov 2010 A1
20100296666 Lin Nov 2010 A1
20100296668 Lee et al. Nov 2010 A1
20100310086 Magrath et al. Dec 2010 A1
20100322430 Isberg Dec 2010 A1
20110007907 Park et al. Jan 2011 A1
20110106533 Yu May 2011 A1
20110129098 Delano et al. Jun 2011 A1
20110130176 Magrath et al. Jun 2011 A1
20110142247 Fellers et al. Jun 2011 A1
20110144984 Konchitsky Jun 2011 A1
20110158419 Theverapperuma et al. Jun 2011 A1
20110206214 Christoph et al. Aug 2011 A1
20110222698 Asao et al. Sep 2011 A1
20110249826 Van Leest Oct 2011 A1
20110288860 Schevciw et al. Nov 2011 A1
20110293103 Park et al. Dec 2011 A1
20110299695 Nicholson Dec 2011 A1
20110305347 Wurm Dec 2011 A1
20110317848 Ivanov et al. Dec 2011 A1
20120135787 Kusunoki et al. May 2012 A1
20120140917 Nicholson et al. Jun 2012 A1
20120140942 Loeda Jun 2012 A1
20120140943 Hendrix et al. Jun 2012 A1
20120148062 Scarlett et al. Jun 2012 A1
20120155666 Nair Jun 2012 A1
20120170766 Alves et al. Jul 2012 A1
20120207317 Abdollahzadeh Milani et al. Aug 2012 A1
20120215519 Park et al. Aug 2012 A1
20120250873 Bakalos et al. Oct 2012 A1
20120259626 Li et al. Oct 2012 A1
20120263317 Shin et al. Oct 2012 A1
20120300958 Klemmensen Nov 2012 A1
20120300960 Mackay et al. Nov 2012 A1
20120308021 Kwatra et al. Dec 2012 A1
20120308024 Alderson et al. Dec 2012 A1
20120308025 Hendrix et al. Dec 2012 A1
20120308026 Kamath et al. Dec 2012 A1
20120308027 Kwatra Dec 2012 A1
20120308028 Kwatra et al. Dec 2012 A1
20120310640 Kwatra et al. Dec 2012 A1
20130010982 Elko et al. Jan 2013 A1
20130083939 Fellers et al. Apr 2013 A1
20130243198 Van Rumpt Sep 2013 A1
20130243225 Yokota Sep 2013 A1
20130272539 Kim et al. Oct 2013 A1
20130287218 Alderson et al. Oct 2013 A1
20130287219 Hendrix et al. Oct 2013 A1
20130301842 Hendrix et al. Nov 2013 A1
20130301846 Alderson et al. Nov 2013 A1
20130301847 Alderson et al. Nov 2013 A1
20130301848 Zhou et al. Nov 2013 A1
20130301849 Alderson et al. Nov 2013 A1
20130343556 Bright Dec 2013 A1
20130343571 Rayala et al. Dec 2013 A1
20140044275 Goldstein et al. Feb 2014 A1
20140050332 Nielsen et al. Feb 2014 A1
20140086425 Jensen et al. Mar 2014 A1
20140177851 Kitazawa et al. Jun 2014 A1
20150010403 Kwatra et al. Apr 2015 A1
20150092953 Abdollahzadeh Milani et al. Apr 2015 A1
Foreign Referenced Citations (21)
Number Date Country
102011013343 Sep 2012 DE
1880699 Jan 2008 EP
1947642 Jul 2008 EP
2133866 Dec 2009 EP
2216774 Aug 2011 EP
2395500 Dec 2011 EP
2395501 Dec 2011 EP
2401744 Nov 2004 GB
2455821 Jun 2009 GB
2455824 Jun 2009 GB
2455828 Jun 2009 GB
2484722 Apr 2012 GB
H06-186985 Jul 1994 JP
03015074 Feb 2003 WO
WO03015275 Feb 2003 WO
WO2004009007 Jan 2004 WO
WO2004017303 Feb 2004 WO
2007007916 Jan 2007 WO
2007113487 Nov 2007 WO
2010117714 Oct 2010 WO
2012134874 Oct 2012 WO
Non-Patent Literature Citations (53)
Entry
Kates, James M., “Principles of Digital Dynamic Range Compression,” Trends in Amplification, Spring 2005, pp. 45-76, vol. 9, No. 2, Sage Publications.
James G. Ryan, Rafik A. Goubran, “Optimum near-field performance of microphone arrays subject to a far-field beampattern constraint”, 2248 J. Acoust. Soc. Am. 108, Nov. 2000.
I. Cohen, B. Berdugo, “Noise Estimation by Minima Controlled Recursive Averaging for Robust Speech Enhancement”, IEEE Signal Processing Letters, vol. 9, No. 1, Jan. 2002.
R. Martin, “Noise Power Spectral Density Estimation Based on Optimal Smoothing and Minimum Statistics”, IEEE Trans. on Speech and Audio Processing, Col. 9, No. 5, Jul. 2001.
R. Martin, “Spectral Subtraction Based on Minimum Statistics”, Proc. 7th EUSIPCO '94, Edinburgh, U.K., Sep. 13-16, 1994, pp. 1182-1195.
I. Cohen, “Noise Spectrum Estimation in Adverse Environments: Improved Minima Controlled Recursive Averaging”, IEEE Trans. on Speech & Audio Proc., vol. 11, Issue 5, Sep. 2003.
A.A. Miliani, G. Kannan, and I.M.S. Panahi, “On maximum achievable noise reduction in ANC systems”, in Proc. ICASSP, 2010, pp. 349-352, Mar. 2010.
Jin, et al. “A simultaneous equation method-based online secondary path modeling algorithm for active noise control”, Journal of Sound and Vibration, Apr. 25, 2007, pp. 455-474, vol. 303, No. 3-5, London, GB.
International Search Report and Written Opinion of the International Searching Authority, International Patent Application No. PCT/US2014/017343, mailed Aug. 8, 2014, 22 pages.
International Search Report and Written Opinion of the International Searching Authority, International Patent Application No. PCT/US2014/018027, mailed Sep. 4, 2014, 14 pages.
International Search Report and Written Opinion of the International Searching Authority, International Patent Application No. PCT/US2014/017374, mailed Sep. 8, 2014, 13 pages.
International Search Report and Written Opinion of the International Searching Authority, International Patent Application No. PCT/US2014/019395, mailed Sep. 9, 2014, 14 pages.
International Search Report and Written Opinion of the International Searching Authority, International Patent Application No. PCT/US2014/019469, mailed Sep. 12, 2014, 13 pages.
Feng, Jinwei et al., “A broadband self-tuning active noise equaliser”, Signal Processing, Elsevier Science Publishers B.V. Amsterdam, NL, vol. 62, No. 2, Oct. 1, 1997, pp. 251-256.
Zhang, Ming et al., “A Robust Online Secondary Path Modeling Method with Auxiliary Noise Power Scheduling Strategy and Norm Constraint Manipulation”, IEEE Transactions on Speech and Audio Processing, IEEE Service Center, New York, NY, vol. 11, No. 1, Jan. 1, 2003.
Lopez-Gaudana, Edgar et al., “A hybrid active noise cancelling with secondary path modeling”, 51st Midwest Symposium on Circuits and Systems, 2008, MWSCAS 2008, Aug. 10, 2008, pp. 277-280.
Widrow, B. et al., Adaptive Noise Cancelling; Principles and Applications, Proceedings of the IEEE, IEEE, New York, NY, U.S. vol. 63, No. 13, Dec. 1975, pp. 1692-1716.
Morgan, Dennis R. et al., A Delayless Subband Adaptive Filter Architecture, IEEE Transactions on Signal Processing, IEEE Service Center, New York, New York. US, vol. 43, No. 8, Aug. 1995, pp. 1819-1829.
International Search Report and Written Opinion of the International Searching Authority, International Patent Application No. PCT/US2014/040999, mailed Oct. 18, 2014, 12 pages.
International Search Report and Written Opinion of the International Searching Authority, International Patent Application No. PCT/US2034/049407, mailed Jun. 18, 2914, 13 pages.
Campbell, Mikey, “Apple looking into self-adjusting earbud headphones with noise cancellation tech”, Apple Insider, Jul. 4, 2013, pp. 1-10 (10 pages in pdf), downloaded on May 14, 2014 from http://appleinsider.com/articles/13/07/04/apple-looking-into-self-adjusting-earbud-headphones-with-noise-cancellation-tech.
Erkelens et al., “Tracking of Nonstationary Noise Based on Data-Driven Recursive Noise Power Estimation”, IEEE Transactions on Audio Speech, and Language Processing, vol. 16, No. 6, Aug. 2008.
Rao et al., “A Novel Two Stage Single Channle Speech Enhancement Technique”, India Conference (INDICON) 2011 Annual IEEE, IEEE, Dec. 15, 2011.
Rangachari et al., “A noise-estimation algorithm for highly non-stationary environments” Speech Communication, Elsevier Science Publishers, vol. 48, No. 2, Feb. 1, 2006.
Hurstm et al., “An improved double sampling scheme for switched-capacitator delta-sigma modulators”, 1992 IEEE Int. Symp. on Circuits and Systems, May 10-13, 1992, vol. 3, pp. 1179-1182, San Diego, CA.
Senderowicz, et al., “Low-Voltage Double-Sampled Delta-Sigma Converters”, IEEE Journal on Solid-State Circuits, Dec. 1997, pp. 1907-1919, vol. 32, No. 12, Piscataway, NJ.
Lopez-Caudana, Edgar Omar, “Active Noise Cancellation: The Unwanted Signal and the Hybrid Solution”, Adaptive Filtering Applications, Dr. Lino Garcia (Ed.), Jul. 2011, pp. 49-84, ISBN: 978-953-307-306-4, InTech.
Kuo, et al., “Residual noise shaping technique for active noise control systems”, J. Acoust. Soc. Am. 95 (3), Mar. 1994, pp. 1665-1668.
Booij, et al., “Virtual sensors for local, three dimensional, broadband multiple-channel active noise control and the effects on the quiet zones”, Proceedings of the International Conference on Noise and Vibration Engineering, ISMA 2010, Sep. 20-22, 2010, pp. 151-166, Leuven.
Black, John W., “An Application of Side-Tone in Subjective Tests of Microphones and Headsets”, Project Report No. NM 001 064.01.20, Research Report of the U.S. Naval School of Aviation Medicine, Feb. 1, 1954, 12 pages. (pp. 1-12 in pdf), Pensacola, FL, US.
Peters, Robert W., “The Effect of High-Pass and Low-Pass Filtering of Side-Tone Upon Speaker Intelligibility”, Project Report No. NM 001 064.01.25, Research Report of the U.S. Naval School of Aviation Medicine, Aug. 16, 1954, 13 pages. (pp. 1-13 in pdf), Pensacola, FL, US.
Lane, et al., “Voice Level: Autophonic Scale, Perceived Loudness, and the Effects of Sidetone”, The Journal of the Acoustical Society of America, Feb. 1961, pp. 160-167, vol. 33, No. 2., Cambridge, MA, US.
Liu, et al., “Compensatory Responses to Loudness-shifted Voice Feedback During Production of Mandarin Speech”, Journal of the Acoustical Society of America, Oct. 2007, pp. 2405-2412, vol. 122, No. 4.
Paepcke, et al., “Yelling in the Hall: Using Sidetone to Address a Problem with Mobile Remote Presence Systems”, Symposium on User Interface Software and Technology, Oct. 16-19, 2011, 10 pages. (pp. 1-10 in pdf), Santa Barbara, CA, US.
Therrien, et al., “Sensory Attenuation of Self-Produced Feedback: The Lombard Effect Revisited”, PLOS ONE, Nov. 2012, pp. 1-7, vol. 7, Issue 11, e49370, Ontario, Canada.
Toochinda, et al. “A Single-Input Two-Output Feedback Formulation for ANC Problems,” Proceedings of the 2001 American Control Conference, Jun. 2001, pp. 923-928, vol. 2, Arlington, VA.
Kuo, et al., “Active Noise Control: A Tutorial Review,” Proceedings of the IEEE, Jun. 1999, pp. 943-973, vol. 87, No. 6, IEEE Press, Piscataway, NJ.
Johns, et al., “Continuous-Time LMS Adaptive Recursive Filters,” IEEE Transactions on Circuits and Systems, Jul. 1991, pp. 769-778, vol. 38, No. 7, IEEE Press, Piscataway, NJ.
Shoval, et al., “Comparison of DC Offset Effects in Four LMS Adaptive Algorithms,” IEEE Transactions on Circuits and Systems II: Analog and Digital Processing, Mar. 1995, pp. 176-185, vol. 42, Issue 3, IEEE Press, Piscataway, NJ.
Mali, Dilip, “Comparison of DC Offset Effects on LMS Algorithm and its Derivatives,” International Journal of Recent Trends in Engineering, May 2009, pp. 323-328, vol. 1, No. 1, Academy Publisher.
Gao, et al., “Adaptive Linearization of a Loudspeaker,” IEEE International Conference on Acoustics, Speech, and Signal Processing, Apr. 14-17, 1991, pp. 3589-3592, Toronto, Ontario, CA.
Silva, et al., “Convex Combination of Adaptive Filters With Different Tracking Capabilities,” IEEE International Conference on Acoustics, Speech, and Signal Processing, Apr. 15-20, 2007, pp. III 925-928, vol. 3, Honolulu, HI, USA.
Akhtar, et al., “A Method for Online Secondary Path Modeling in Active Noise Control Systems,” IEEE International Symposium on Circuits and Systems, May 23-26, 2005, pp. 264-267, vol. 1, Kobe, Japan.
Davari, et al., “A New Online Secondary Path Modeling Method for Feedforward Active Noise Control Systems,” IEEE International Conference on Industrial Technology, Apr. 21-24, 2008, pp. 1-6, Chengdu, China.
Lan, et al., “An Active Noise Control System Using Online Secondary Path Modeling With Reduced Auxiliary Noise,” IEEE Signal Processing Letters, Jan. 2002, pp. 16-18, vol. 9, Issue 1, IEEE Press, Piscataway, NJ.
Liu, et al., “Analysis of Online Secondary Path Modeling With Auxiliary Noise Scaled by Residual Noise Signal,” IEEE Transactions on Audio, Speech and Language Processing, Nov. 2010, pp. 1978-1993, vol. 18, Issue 8, IEEE Press, Piscataway, NJ.
Pfann, et al., “LMS Adaptive Filtering with Delta-Sigma Modulated Input Signals”, IEEE Signal Processing Letters, vol. 5, No. 4, Apr. 1998.
Milani, et al., “On Maximum Achievable Noise Reduction in ANC Systems”, Proceedings of the IEEE International Conference on Acoustics, Speech, and Signal Processing, ICASSP 2010, Mar. 14-19, 2010 pp. 349-352.
Ryan, et al., “Optimum near-field performance of microphone arrays subject to a far-field beampattern constraint”, 2248 J. Acoust. Soc. Am. 108, Nov. 2000.
Cohen, et al., “Noise Estimation by Minima Controlled Recursive Averaging for Robust Speech Enhancement”, IEEE Signal Processing Letters, vol. 9, No. 1, Jan. 2002.
Martin, “Noise Power Spectral Density Estimation Based on Optimal Smoothing and Minimum Statistics”, IEEE Trans. on Speech and Audio Processing, col. 9, No. 5, Jul. 2001.
Martin, “Spectral Subtraction Based on Minimum Statistics”, Proc. 7th EUSIPCO '94, Edinburgh, U.K., Sep. 13-16, 1994, pp. 1182-1195.
Cohen, “Noise Spectrum Estimation in Adverse Environments: Improved Minima Controlled Recursive Averaging”, IEEE Trans. on Speech & Audio Proc., vol. 11, Issue 5, Sep. 2003.
Related Publications (1)
Number Date Country
20140226827 A1 Aug 2014 US