CORE ESTIMATOR AND ADAPTIVE GAINS FROM SIGNAL TO NOISE RATIO IN A HYBRID SPEECH ENHANCEMENT SYSTEM

Abstract

A speech enhancement system receives noisy speech and produces enhanced speech. The noisy speech is characterized by a spectral amplitude spanning a plurality of frequency bins. The speech enhancement system modifies the spectral amplitude of the noisy speech without affecting the phase of the noisy speech. The speech enhancement system includes a core estimator that applies to the noisy speech one of a first set of gains for each frequency bin. A noise adaptation module segments the noisy speech into noise-only and signal-containing frames, maintains a current estimate of the noise spectrum and an estimate of the probability of signal absence in each frequency bin. A signal-to-noise ratio estimator measures an a-posteriori signal-to-noise ratio and estimates an a-priori signal-to-noise ratio based on the noise estimate. Each one of the first set of gains is based on the a-priori signal-to-noise ratio, as well as the probability of signal absence in each bin and a level of aggression of the speech enhancement. A soft decision module computes a second set of gains that is based on the a-posteriori signal-to-noise ratio and the a-priori signal-to-noise ratio, and the probability of signal absence in each frequency bin.

Description

BACKGROUND OF THE INVENTION

[0002] There are many environments where noisy conditions interfere with speech, such as the inside of a car, a street, or a busy office. The severity of background noise varies from the gentle hum of a fan inside a computer to a cacophonous babble in a crowded cafe. This background noise not only directly interferes with a listener's ability to understand a speaker's speech, but can cause further unwanted distortions if the speech is encoded or otherwise processed. Speech enhancement is an effort to process the noisy speech for the benefit of the intended listener, be it a human, speech recognition module, or anything else. For a human listener, it is desirable to increase the perceptual quality and intelligibility of the perceived speech, so that the listener understands the communication with minimal effort and fatigue.

[0003] It is usually the case that for a given speech enhancement scheme, a trade-off must be made between the amount of noise removed and the distortion introduced as a side effect. If too much noise is removed, the resulting distortion can result in listeners preferring the original noise scenario to the enhanced speech. Preferences are based on more than just the energy of the noise and distortion: unnatural sounding distortions become annoying to humans when just audible, while a certain elevated level of “natural sounding” background noise is well tolerated. Residual background noise also serves to perceptually mask slight distortions, making its removal even more troublesome.

[0004] Speech enhancement can be broadly defined as the removal of additive noise from a corrupted speech signal in an attempt to increase the intelligibility or quality of speech. In most speech enhancement techniques, the noise and speech are generally assumed to be uncorrelated. Single channel speech enhancement is the simplest scenario, where only one version of the noisy speech is available, which is typically the result of recording someone speaking in a noisy environment with a single microphone.

[0005] FIG. 1 illustrates a speech enhancement setup for N noise sources for a single-channel system. For the single channel case illustrated in FIG. 1, exact reconstruction of the clean speech signal is usually impossible in practice. So speech enhancement algorithms must strike a balance between the amount of noise they attempt to remove and the degree of distortion that is introduced as a side effect. Since any noise component at the microphone cannot in general be distinguished as coming from a specific noise source, the sum of the responses at the microphone from each noise source is denoted as a single additive noise term.

[0006] Speech enhancement has a number of potential applications. In some cases, a human listener observes the output of the speech enhancement directly, while in others speech enhancement is merely the first stage in a communications channel and might be used as a preprocessor for a speech coder or speech recognition module. Such a variety of different application scenarios places very different demands on the performance of the speech enhancement module, so any speech enhancement scheme ought to be developed with the intended application in mind. Additionally, many well-known speech enhancement processes perform very differently with different speakers and noise conditions, making robustness in design a primary concern. Implementation issues such as delay and computational complexity are also considered.

I. Modified MMSE-LSA Approach

[0007] The modified Minimum Mean-Square Error Log-Spectral Amplitude (modified MMSE-LSA) estimator for speech enhancement was designed by David Malah and draws upon three main ideas: the Minimum Mean Square Error Log-Spectral Amplitude (MMSE-LSA) estimator (Y. Ephraim and D. Malah, “Speech Enhancement Using a Minimum Mean-Square Error Log-Spectral Amplitude Estimator,” IEEE Transactions on Acoustics, Speech, and Signal Processing, vol. ASSP-33, pp. 443-445, 1985); the soft decision approach (R. J. McAulay and M. L. Malpass, “Speech Enhancement Using a Soft-Decision Noise Suppression Filter,” IEEE Transactions on Acoustics, Speech, and Signal Processing, vol. ASSP-28, pp. 137-145, 1980); and a novel noise adaptation scheme. The modified MMSE-LSA speech enhancement system is a member of the class of STSA enhancement techniques and is schematically depicted in FIG. 2.

[0008] With reference to FIG. 2, the MMSE-LSA estimator 10 operates in the frequency domain and applies a gain to each DFT coefficient of the noisy speech that is computed from signal-to-noise ratio (SNR) estimates 12. A soft decision module 14 applies an additional gain in the frequency domain that accounts for signal presence uncertainty. A noise adaptation scheme 16 supplies estimates of current noise characteristics for use in the SNR calculations.

I.A. The MMSE-LSA Estimator

[0009] We begin by assuming additive independent noise and that the DFT coefficients of both the clean speech and the noise are zero-mean, statistically independent, Gaussian random variables. We formulate the speech enhancement problem as

y[n]=x[n]+w[n]   (1)

[0010] Taking the DFT of (1), we obtain

Y k =X k +W k   (2)

[0011] We express the complex clean and noisy speech DFT coefficients in exponential form as

Xk=AkeJφk  (3)

Yk=RkeJθk  (4)

[0012] Now the MMSE-LSA estimate of Ak is the amplitude that minimizes the difference between log Ak and the logarithm of that amplitude in a MMSE sense: 1$\begin{array}{cc}{\hat{A}}_k=\underset{B}{\arg \min }E\left[{\left(\log A_k-\log B\right)}^2\right]& \left(5\right)\end{array}$

[0013] The solution to (5) is the exponential of the conditional expectation (A. Papoulis, Probability, Random Variables, and Stochastic Processes, 3 ed. New York: McGraw-Hill, Inc., 1991):

 k =exp (E[log Ak|Yk])  (6)

[0014] Therefore, to implement the MMSE-LSA estimator 10, we must scale the noisy speech DFT coefficients Yk so that they have the estimated amplitude Âk. Our estimate of the clean speech in the frequency domain is now 2$\begin{array}{cc}{\hat{X}}_k={\hat{A}}_k\frac{Y_k}{\&LeftBracketingBar;Y_k\&RightBracketingBar;}& \left(7\right)\end{array}$

[0015] We are using the “noisy phase” in (7), since the phase of the DFT coefficients of the noisy speech is used in our estimate of the clean speech. The MMSE complex exponential estimator does not have a modulus of 1. (Y. Ephraim and D. Malah, “Speech Enhancement Using a Minimum Mean-Square Error Short-Time Spectral Amplitude Estimator,” IEEE Transactions on Acoustics, Speech, and Signal Processing, vol. ASSP-32, pp. 1109-1121, 1984). So when an optimal complex exponential estimator is combined with an optimal amplitude estimator, the resulting amplitude estimate is no longer optimal. When the estimate's modulus is constrained to be unity, however, the MMSE complex exponential estimator is the exponent of the noisy phase. In addition, the optimal estimator of the principal value of the phase is the noisy phase itself. This provides justification for using the MMSE-LSA estimator 10 to estimate Ak and to leave the noisy phase untouched, as indicated in (7).

[0016] The computation of the expectation in (6) is non-trivial and presented in the article by Y. Ephraim and D. Malah, “Speech Enhancement Using a Minimum Mean-Square Error Log-Spectral Amplitude Estimator,” IEEE Transactions on Acoustics, Speech, and Signal Processing, vol. ASSP-33, pp. 443-445, 1985, where Âk is shown to be:

 k =G (ξk,γk)·Rk  (8)

[0017] where 3$\begin{array}{cc}G\left({\xi }_k,{\gamma }_k\right)=\frac{{\xi }_k}{1+{\xi }_k}\exp \left(\frac{1}{2}{\int }_{{\upsilon }_k}^{\infty }\frac{e^{-t}}{t}dt\right)& \left(9\right)\\ {\upsilon }_k=\frac{{\xi }_k}{1+{\xi }_k}{\gamma }_k& \left(10\right)\end{array}$

 ξk=λx(k)/λw(k)  (11)

γk=Rk2/λw(k)  (12)

λx(k)=E[|Xk|2]=E[|Ak|2]  (13)

λw(k)=E[|Wk|2]  (14)

[0018] Here λx(k) and λw(k) defined in (13) and (14) are the energy spectral coefficients of the clean speech and the noise, respectively. As defined in (11) and (12), the quantities εk and γk can be interpreted as signal-to-noise ratios. We will denote εk as the a-priori SNR, as it is the ratio of the energy spectrum of speech to that of the noise prior to the contamination of the speech by the noise. Similarly, we will call γk the a-posteriori SNR, as it is the ratio of the energy of the current frame of noisy speech to the energy spectrum of the noise, after the speech has been contaminated.

[0019] In order to compute G(εk,γk) as given in (9), we must first estimate these SNR's εk and γk. Malah's noise adaptation scheme 16 provides an estimate of λw(k), so the a-posteriori SNR γk is straightforward to estimate since Rk is readily computed from the noisy speech. However, the a-priori SNR εk is somewhat more difficult to estimate. It turns out that the Maximum Likelihood (ML) estimate of εk does not work very well. In the article by Y. Ephraim and D. Malah, “Speech Enhancement Using a Minimum Mean-Square Error Short-Time Spectral Amplitude Estimator,” IEEE Transactions on Acoustics, Speech, and Signal Processing, vol. ASSP-32, pp. 1109-1121, 1984, the shortcomings of the ML estimate are discussed and a “decision directed” estimation approach is considered. The key idea is that under our assumption of Gaussian DFT coefficients, the a-priori SNR can be expressed in terms of the a-posteriori SNR as

γk=E[γk−1]  (15)

[0020] For each frame of noisy speech, we can then take a convex combination of our two expressions (11) and (15) for εk, after dropping the expectation operators, to obtain an estimate for the a-priori SNR using previous values of Âk and {circumflex over (λ)}k. For the nth frame we have 4$\begin{array}{cc}{\hat{\xi }}_k\left(n\right)=\alpha \frac{{\hat{A}}_k^2\left(n-1\right)}{{\hat{\lambda }}_w\left(k,n-1\right)}+\left(1-\alpha \right)P\left[{\hat{\gamma }}_k\left(n\right)-1\right]\mathrm{where}P\left[x\right]=\{\begin{array}{cc}x& \mathrm{if}x\ge 0\\ 0& \mathrm{otherwise}\end{array}& \left(16\right)\end{array}$

[0021] The P[x] function is used to clip the a-posteriori SNR γk to 1 if a smaller value is calculated, and 0≦α≦1.

[0022] This “decision directed” estimate is mainly responsible for the elimination of musical noise artifacts that plague earlier speech enhancement algorithms. (0. Cappé, “Elimination of the Musical Noise Phenomenon with the Ephraim and Malah Noise Suppressor,” IEEE Transactions on Speech and Audio Processing, vol. 2, pp. 345-349, 1994). The intuition behind this mechanism is that for large a-posteriori SNRs, the a-priori SNR follows the a-posteriori SNR with a single frame delay. This allows the enhancement scheme to adapt quickly to any sudden changes in the noise characteristics that the noise adaptation scheme perceives. However, for small a-posteriori SNRs, the a-priori SNR is a highly smoothed version of the a-posteriori SNR. Since the a-priori SNR has a major impact in determining the gain as seen in (9), there are no sudden fluctuations in gain at any fixed frequency from frame to frame when there is a good deal of noise present. This greatly reduces the musical noise phenomenon.

[0023] We can choose α to trade-off between the degree of noise reduction and the overall distortion. α must be close to 1 (>0.98) in order to achieve the greatest musical noise reduction effect. (O. Cappé, “Elimination of the Musical Noise Phenomenon with the Ephraim and Malah Noise Suppressor,” IEEE Transactions on Speech and Audio Processing, vol. 2, pp. 345-349, 1994). The higher a, however, the more aggressive the algorithm is in removing the residual noise, which causes additional speech distortion. In fact, the easiest way to trade-off between aggression and distortion is through changing a, which has the awkward side effect of disturbing the smoothing properties discussed above.

I.B. Signal Presence Uncertainty

[0024] The above analysis assumes that there is speech present in every frequency bin of every frame of the noisy speech. This is generally not the case, and there are two well-established ways of taking advantage of this situation.

[0025] The first, called “hard decision”, treats the presence of speech in some frequency bin as a time-varying deterministic condition that can be determined using classical detection theory. The second, “soft decision”, treats the presence of speech as a stochastic process with a changing binary probability distribution. (R. J. McAulay and M. L. Malpass, “Speech Enhancement Using a Soft-Decision Noise Suppression Filter,” IEEE Transactions on Acoustics, Speech, and Signal Processing, vol. ASSP-28, pp. 137-145, 1980). The soft decision approach has been found to be more successful in speech enhancement. (Y. Ephraim and D. Malah, “Speech Enhancement Using a Minimum Mean-Square Error Short-Time Spectral Amplitude Estimator,” IEEE Transactions on Acoustics, Speech, and Signal Processing, vol. ASSP-32, pp. 1109-1121, 1984). A hard decision approach can in fact lead to musical noise. When the decision oscillates between signal presence and absence in time for some frequency bin, an enhancement scheme that greedily eliminates frequency components containing only noise would produce tonal artifacts at that frequency. Following this outline, we define two states for each frequency bin k. H0k denotes the state where the speech signal is absent in the kth bin, while H1k is the state where the signal is present in the kth bin. Now our estimate of log Ak is given by

E[log A k |Y k ,H 1 k ]Pr (H1k|Yk)+E[log Ak|Yk,H0k]Pr(H0k|Yk)  (17)

[0026] Since E[log Ak|Yk,H0k]=0, soft decision entails weighting our previous estimate of log Ak by Pr(H1k|Yk). To compute this weighting factor, we first expand Pr(H1k,Yk) in two different ways:

Pr (H1k|Yk)·Pr(Yk)=Pr(Yk|H1k)·Pr(H1k)  (18)

[0027] Also,

Pr (Yk)=Pr(Yk|H1k)·Pr(H1k)+Pr(Yk|H0k)·Pr(H0k)  (19)

[0028] From (18) and (19) we can solve for Pr(H1k|Yk) and express it in terms of the likelihood function Λ(k): 5$\begin{array}{cc}\mathrm{Pr}\left(H_1^k|Y_k\right)=\frac{\Lambda \left(k\right)}{1+\Lambda \left(k\right)}\mathrm{where}& \left(20\right)\\ \Lambda \left(k\right)={\mu }_k\frac{\mathrm{Pr}\left(Y_k|H_1^k\right)}{\mathrm{Pr}\left(Y_k|H_0^k\right)}& \left(21\right)\\ {\mu }_k=\mathrm{Pr}\left(H_1^k\right)/\mathrm{Pr}\left(H_0^k\right)=\frac{1-q_k}{q_k}& \left(22\right)\end{array}$

[0029] Here qk is the a-priori probability of signal absence in the kth bin, and A(k) is clearly the likelihood function from classical detection theory. (A. Papoulis, Probability, Random Variables, and Stochastic Processes, 3 ed. New York: McGraw-Hill, Inc., 1991). With our Gaussian distribution assumptions on Yk, it is straightforward to calculate Λ(k): 6$\begin{array}{cc}\Lambda \left(k\right)=\frac{1-q_k}{q_k}·\frac{1}{1+{\eta }_k}\exp {\left(\frac{{\eta }_k}{1+{\eta }_k}{\gamma }_k\right)}_{{\eta }_k=\frac{{\xi }_k}{1-q_k}}& \left(23\right)\end{array}$

[0030] where the SNR's γk and εk can be estimated in the same manner as described in Section I.A.

I.C. Noise Adaptation

[0031] An important development for the modified MMSE-LSA speech enhancement technique is the noise adaptation scheme 16, which allows the speech enhancement technique to handle non-stationary noise. The adaptation proceeds in two steps; the first identifies all the spectral coefficients in the current frame that are reasonably good representations of the noise, and the second adapts the current noise estimate to this new information.

[0032] Direct spectral information about the noise can become available when a frame of the noisy speech is a “noise-only” frame, meaning that the speech contribution during that time period is negligible. In this case, the entire noise spectrum estimate can be updated. Additionally, even if a frame contains both speech and noise, there may still be some “noise-only” frequency bins so that the speech contribution within certain frequency ranges is negligible during the current frame. Here we can update the corresponding spectral components of our noise estimate accurately.

[0033] The process of deciding whether a given frame is a noise-only frame is dubbed “segmentation”, and the decision is based on the a-posteriori SNR estimates γk. Under our Gaussian distribution assumptions on Yk, we can compute the probability density function ƒ(γk) for γk, which turns out to be an exponential distribution with mean and standard deviation 1+εk given by 7$\begin{array}{cc}f\left({\gamma }_k\right)=\frac{1}{1+{\xi }_k}\exp \left(\frac{-{\gamma }_k}{1+{\xi }_k}\right)& \left(24\right)\end{array}$

[0034] We declare a frame of speech to be noise-only if both our average (over k) estimate of the a-posteriori SNRs is low and the average of our estimate of the variance of the a-posteriori SNR estimator is low. That is, a frame is noise-only when

{overscore (γ)}≦{overscore (γ)}Threshold and {overscore (ξ)}≦σThreshold−1  (25)

[0035] When a noise-only frame is discovered, we update all the spectral components of our noise estimate by averaging our estimates for the previous frame with our new estimates. So our noise spectral estimate for the kth frequency bin and the nth frame is given by:

{circumflex over (λ)}w(k,n)=αw{circumflex over (λ)}w(k,n−1)+(1−αw)Rw2  (26)

[0036] where αw is the forgetting factor of the update equation, which is dynamically updated based on the average estimate of γk. In this manner, the forgetting factor is directly related to the current value of {circumflex over (γ)} so that the lower {circumflex over (γ)} is, the better our estimate of the noise spectrum, and therefore we discard our previous noise spectral estimates more quickly.

[0037] The situation for dealing with noise-only frequency bins for frames with signal present is quite similar, except the individual SNR estimates for each frequency bin are used instead of their averages. There is one main difference; since we have an estimate of the probability that each bin contains no signal present (qk from our soft decision discussion in Section I.B.), we can use this to refine our update of the forgetting factor for each frequency bin.

[0038] The impact of this noise adaptation scheme 16 is dramatic. The complete modified MMSE-LSA enhancement technique is capable of adapting to great changes in noise volume in only a few frames of speech, and has demonstrated promising performance in dealing with highly non-stationary noise, such as music.

II. Signal Subspace Approach

[0039] Yariv Ephraim and Harry L. Van Trees developed a signal subspace approach (Y. Ephraim and H. L. V. Trees, “A Signal Subspace Approach for Speech Enhancement,” IEEE Transactions on Speech and Audio Processing, vol. 3, pp. 251-266, 1995) that provides a theoretical framework for understanding a number of classical speech enhancement techniques, and allows for the application of external criteria to control enhancement performance. The basic idea is that the vector space of the noisy speech can be decomposed into a signal-plus-noise subspace and a noise-only subspace. Once identified, the noise-only subspace can be eliminated and then the speech estimated from the remaining signal-plus-noise subspace. We assume that the full space has dimension K and the signal-plus-noise subspace has dimension M<K.

[0040] Say we have clean speech x[n] that is corrupted by independent additive noise w[n] to produce a noisy speech signal y[n]. We constrain ourselves to estimating x[n] using a linear filter H, and will initially consider w[n] to be a white noise process with variance σw2. In vector notation, we have

y=x+w   (27)

[0041] {circumflex over (x)}=Hy  (28)

[0042] We can decompose the residual error into a term solely dependent on the clean speech, called the signal distortion rx, and a term solely dependent on the noise, called the residual noise rw: 8$\begin{array}{cc}\begin{array}{c}r=\hat{x}-x\\ =\left(H-I\right)x+\mathrm{Hw}\\ =r_x+r_w\end{array}& \left(29\right)\end{array}$

[0043] In (29) we have explicitly identified the trade-off between residual noise and speech distortion. Since different applications could require different trade-offs between these two factors, it is desirable to perform a constrained minimization using functions of the distortion and residual noise vectors. Then the constraints can be selected to meet the application requirements.

II.A. Time Domain Constrained Estimator

[0044] Two different frameworks for performing a constrained minimization using functions of the residual noise and signal distortion are presented in the article by Y. Ephraim and H. L. V. Trees, “A Signal Subspace Approach for Speech Enhancement,” IEEE Transactions on Speech and Audio Processing, vol. 3, pp. 251-266, 1995. The first examines the energy in these vectors and results in a time domain constrained estimator. We define

{overscore (ε)}x2=trE[rxrx#]=tr{(H−I)Ry(H−I)#}  (30)

[0045] to be the energy of the signal distortion vector rx, and similarly define

{overscore (ε)}w2=trE[rwrw#]=σw2tr{HH#}  (31)

[0046] to be the energy of the residual noise vector rw.

[0047] We desire to minimize the energy of the signal distortion while constraining the energy of the residual noise to be less than some fraction Kα of the noise variance σw2: 9$\begin{array}{cc}\underset{H}{\min }{\stackrel{\_}{\epsilon }}_x^2\mathrm{subject}\mathrm{to}{\stackrel{\_}{\epsilon }}_w^2/K\le {\mathrm{\alpha \sigma }}_w^2& \left(32\right)\end{array}$

[0048] The solution to the constrained minimization problem in (32) involves first the projection of the noisy speech signal onto the signal-plus-noise subspace, followed by a gain applied to each eigenvalue, and finally the reconstruction of the signal from the signal-plus-noise subspace. The gain for the mth eigenvalue is a function of the Lagrange multiplier μ, and is given by 10$\begin{array}{cc}{g}_{\mu }\left(m\right)=\frac{{\lambda }_x\left(m\right)}{{\lambda }_x\left(m\right)+{\mathrm{\mu \sigma }}_w^2}& \left(33\right)\end{array}$

[0049] where λx(m) is the mth eigenvalue of the clean speech.

[0050] Thus, the enhancement system, which is schematically illustrated in FIG. 3, can be implemented as a Karhunen-Loève Transform (KLT) 24 which receives a noisy signal, followed by a set of gains (G1, . . . , GN) 26, and ending with an inverse KLT 28 which outputs an enhanced signal.

[0051] Ephraim shows that μ is uniquely determined by our choice of the constraint α, and demonstrates how the generalized Wiener filter in (33) can implement linear MMSE estimation and spectral subtraction for specific values of μ and certain approximations to the KLT.

II.B. Spectral Domain Constrained Estimator

[0052] To provide a tighter means of control over the trade-off between residual noise and signal distortion, Ephraim derives a spectral domain constrained estimator (Y. Ephraim and H. L. V. Trees, “A Signal Subspace Approach for Speech Enhancement,” IEEE Transactions on Speech and Audio Processing, vol. 3, pp. 251-266, 1995) which minimizes the energy of the signal distortion while constraining each of the eigenvalues of the residual noise by a different constant proportion of the noise variance: 11$\begin{array}{cc}\underset{H}{\min }{\stackrel{\_}{\epsilon }}_x^2\mathrm{subject}\mathrm{to}E\left[{\&LeftBracketingBar;u_k^{\#}{r}_w\&RightBracketingBar;}^2\right]\le {\alpha }_k{\sigma }_w^2& \left(34\right)\end{array}$

[0053] Here uk is the kth eigenvector of the noisy speech, and the constraint is applied for each k in the signal-plus-noise subspace. The form of the solution to this constrained minimization is very similar to the time domain constrained estimator illustrated in FIG. 3; the only difference is that the eigenvalue gains are given by

g (m)={square root}{square root over (αk)}  (35)

[0054] instead of the result in (33).

[0055] Now with such freedom over the constraints αk, the difficulty arises as to how to optimally choose these constants to obtain a reasonable speech enhancement system. One choice Ephraim investigated is

αk=exp{−νσw2/λx(k)}  (36)

[0056] where ν is a constant that determines the level of noise suppression, or the aggression level of the enhancement algorithm. The constraints in (36) effectively shape the noise so it resembles the clean speech, which takes advantage of the masking properties of the human auditory system. This choice of functional form for αk is an aggressive one.

[0057] There is no treatment of noise distortion in this signal subspace approach, and it turns out that the residual noise in the enhanced signal can contain artifacts so annoying that the result is less desirable than the original noisy speech. Therefore, when using this signal subspace framework it is desirable to aggressively reduce the residual noise at the possibly severe cost of increased signal distortion.

II.C. Reverse Spectral Domain Constrained Estimator

[0058] The spectral domain constrained estimator can be placed in a framework that will substantially reduce the noise distortion. In such scenarios, it might be advantageous to use a variant of Ephraim's spectral domain constrained estimator. Here we minimize the residual noise with the signal distortion constrained: 12$\begin{array}{cc}\underset{H}{\min }{\stackrel{\_}{\epsilon }}_w^2\mathrm{such}\mathrm{that}E\left[{\&LeftBracketingBar;u_k^{\#}{r}_y\&RightBracketingBar;}^2\right]\le {\alpha }_k{\lambda }_{y,k}& \left(37\right)\end{array}$

[0059] Since H could have complex entries, we set the Jacobians of both the real and imaginary parts of the Lagrangian from (37) to zero in order to obtain the first order conditions, expressed in matrix form as

HR w +UΛ μ U # (H−I)Ry=0  (38)

[0060] where Λμ=diag(μ1, . . . , μK) is a diagonal matrix of Lagrange multipliers. Applying the eigendecomposition of Ry and using the assumption that the noise is white, we obtain:

σw2Q+ΛμQΛy=ΛμΛy  (39)

[0061] where

Q=U#HU  (40)

[0062] We note that a possible solution to the constrained minimization is obtained when Q is diagonal with elements given by 13$\begin{array}{cc}{q}_{\mathrm{kk}}=\{\begin{array}{cc}\frac{{\mu }_k{\lambda }_{y,k}}{{\sigma }_w^2+{\mu }_k{\lambda }_{y,k}}& k=1,\dots ,M\\ 0& k=M+1,\dots ,K\end{array}& \left(41\right)\end{array}$

[0063] which satisfies (39). For this Q, we have

E[|u k # r y |2]=λy,k(qkk−1)2  (42)

[0064] Now for the non-zero constraints in (37) to hold with equality, we must have

q kk =1−{square root}{square root over (αk)}  (43)

[0065] and 14$\begin{array}{cc}{\mu }_k=\frac{{\sigma }_w^2}{{\lambda }_{y,k}\sqrt{{\alpha }_k}}\left(1-\sqrt{{\alpha }_k}\right)& \left(44\right)\end{array}$

[0066] Since we see from (44) that μk≧0, this proposed solution satisfies the Kuhn-Tucker necessary conditions for the constrained minimization.

[0067] We conclude that H is given by 15$\begin{array}{cc}H={\mathrm{UQU}}^{\#}Q=\mathrm{diag}\left(q_{11},\dots ,q_{\mathrm{KK}}\right)q_{\mathrm{kk}}=\{\begin{array}{cc}1-\sqrt{{\alpha }_k}& k=1,\dots ,M\\ 0& k=M+1,\dots ,K\end{array}& \left(45\right)\end{array}$

[0068] Thus the reverse spectral domain constrained estimator has a form very similar to that of our previous signal subspace estimators. The implementation of (45) is given in FIG. 3 with the gains

g (m)=1−{square root}{square root over (αk)}  (46)

SUMMARY OF THE INVENTION

[0069] According to an exemplary embodiment of the invention, a speech enhancement system receives noisy speech and produces enhanced speech. The noisy speech is characterized by spectral coefficients spanning a plurality of frequency bins and contains an original noise. The speech enhancement system includes a noise adaptation module. The noise adaptation module receives the noisy speech, and segments the noisy speech into noise-only frames and signal-containing frames. The noise adaptation module determines a noise estimate and a probability of signal absence in each frequency bin. A signal-to-noise ratio (SNR) estimator is coupled to the noise adaptation module. The signal-to-noise ratio estimator determines a first signal-to-noise ratio and a second signal-to-noise ratio based on the noise estimate. A core estimator coupled to the signal-to-noise ratio estimator receives the noisy speech. The core estimator applies to the spectral coefficients of the noisy speech one of a first set of gains for each frequency bin in the frequency domain without discarding the noise-only frames. The core estimator outputs noisy speech having a residual noise.

[0070] Each one of the first set of gains is determined based on the second signal-to-noise ratio, a level of aggression, the probability of signal absence in each frequency bin, or combinations thereof. The core estimator constrains the spectral density of the spectral coefficients of the residual noise to be below a constant proportion of the spectral density of the spectral coefficients of the original noise. A soft decision module coupled to the core estimator and to the signal-to-noise ratio estimator determines a second set of gains that is based on the first signal-to-noise ratio, the second signal-to-noise ratio and the probability of signal absence in each frequency bin. The soft decision module applies the second set of gains to the spectral coefficients of the noisy speech containing the residual noise and outputs enhanced speech.

[0071] According to an aspect of the invention, noisy speech that is characterized by spectral coefficients spanning a plurality of frequency bins and that contains an original noise is enhanced by segmenting the noisy speech into noise-only frames and signal-containing frames and determining a noise estimate and a probability of signal absence in each frequency bin. A first signal-to-noise ratio and a second signal-to-noise ratio are determined based on the noise estimate. A first set of gains is determined based on the second signal-to-noise ratio, a level of aggression, the probability of signal absence in each frequency bin, or combinations thereof. The first set of gains is applied to the spectral coefficients of the noisy speech without discarding the noise-only frames to produce noisy speech containing a residual noise, such that the spectral density of the spectral coefficients of the residual noise is maintained below a constant proportion of the spectral density of the spectral coefficients of the original noise. A second set of gains is applied to the noisy speech containing the residual noise to produce enhanced speech. The spectral amplitude of the noisy speech is modified without affecting the phase of the noisy speech. During a noise-only frame, a constant gain is applied to the noise to avoid noise structuring.

[0072] Other features and advantages of the invention will become apparent from the following detailed description, taken in conjunction with the accompanying drawings, which illustrate, by way of example, the features of the invention.

BRIEF DESCRIPTION OF THE DRAWINGS

[0073] FIG. 1 illustrates a speech enhancement setup for N noise sources for a single-channel system;

[0074] FIG. 2 is a block diagram of a modified MMSE-LSA speech enhancement system;

[0075] FIG. 3 is a block diagram of a signal subspace estimator;

[0076] FIG. 4 is a block diagram of a speech enhancement system in accordance with the principles of the invention;

[0077] FIG. 5 is a block diagram of a first embodiment of the core estimator of the speech enhancement system illustrated in FIG. 4; and

[0078] FIG. 6 is a block diagram of a second embodiment of the core estimator of the speech enhancement system illustrated in FIG. 4.

DETAILED DESCRIPTION

III. Hybrid Speech Enhancement System

[0079] Ephraim's signal subspace approach (see Section II.) and Malah's modified MMSE-LSA algorithm (see Section I.) have very different strengths and weaknesses.

[0080] Ephraim's signal subspace approach provides a simple but powerful framework for trading-off between the degree of noise suppression and signal distortion. This framework is general enough to incorporate many different criteria, including perceptual measures for general applications. This provides a good deal of flexibility when attempting to specialize an enhancement algorithm for a specific application. However, the technique offers no means for controlling noise distortion and handling non-stationary noise. Noise can be so severely distorted that the enhanced signal is less desirable than the original noisy signal, even though the noise energy has been suppressed. This forces one to operate the signal subspace algorithm in a very aggressive mode, so that the noise is practically eliminated but signal distortion may be high.

[0081] Malah's modified MMSE-LSA approach was carefully designed to reduce noise distortion and adapt to non-stationary noise. The approach is quite robust when presented with different types and levels of noise. The main difficulty is that the trade-off between the degree of noise suppression and signal distortion is awkward and is best performed by varying α in (16), which has undesirable side effects on the noise distortion. This provides very little flexibility when trying to adapt the algorithm to fit a particular application.

[0082] The present invention combines the strengths of these two approaches in order to generate a robust and flexible speech enhancement system that performs just as well. FIG. 4 schematically illustrates a speech enhancement system in accordance with the principles of the invention. The speech enhancement system shown in FIG. 4 receives noisy speech and produces enhanced speech. The speech enhancement system includes a noise adaptation processor 34 that receives the noisy speech that contains an original noise. A signal-to-noise ratio (SNR) estimator 36 is coupled to the noise adaptation processor 34 and receives the noisy speech containing the original noise. A core estimator 38 is coupled to the SNR estimator 36 and receives the noisy speech containing the original noise. The core estimator 38 applies a first set of gains in the frequency domain to the noisy speech containing the original noise without discarding noise-only frames, and outputs noisy speech containing a residual noise. A soft decision module 40 is coupled to the core estimator 38 and to the SNR estimator 36. The soft decision module 40 applies a second set of gains to the noisy speech and outputs the enhanced speech.

[0083] The noise adaptation processor 34 acts independently from the remainder of the modules. It is essential for many STSA speech enhancement algorithms to have an accurate estimate of the noise. Malah's modified MMSE-LSA approach, for example, is particularly effective in tracking non-stationary noise, especially noise with varying intensity levels. The decision directed estimation approach is buried in the SNR estimator 36, which smoothes estimates between frames when the SNR becomes poor. We have seen that the effect is to reduce noise distortion when the gain applied depends heavily on these SNR estimates. The soft decision module 40 has broad applicability, and could be considered part of the core estimator 38. Since this technique has proven most effective in handling the uncertainty of signal presence in certain frequency bands for different estimators, we consider the soft decision module 40 to be a separate entity.

III. A. Signal Subspace as a Core Estimator

[0084] Our first insight is that we can substitute anything we desire in the core estimator 38 block of FIG. 4 and take advantage of the supporting structure as long as the effective gain depends heavily on the SNR estimates provided. Our intuition is that this choice of core estimator 38 might depend on the desired application. For our present purpose, however, we will use the spectral domain constrained version of the signal subspace approach as the core estimator 38 in an effort to take advantage of its aggressive noise suppression properties and flexibility.

[0085] We modify the signal subspace approach so as to satisfy our constraints on the core estimator 38. The first modification to the signal subspace approach is using a Discrete Fourier Transform (DFT) in place of the KLT (24, FIG. 3). Since the first step of the signal subspace approach is to decompose the noisy speech into a noise-only subspace and a speech-plus-noise subspace and throw away the noise-only subspace, the approach takes advantage of the uncertainty of signal presence. When the KLT used in the signal subspace estimator is approximated with a Discrete Fourier Transform (DFT), this step is precisely a hard decision with zero gain applied to the frequency bins that contain pure noise. Such an approach leads to unpleasant noise distortion properties. The second modification to the signal subspace approach is to skip this noise-only subspace cancellation step.

[0086] Adapting the signal subspace approach to be a function of our SNR estimates is a bit more troublesome. The first difficulty is that the signal subspace approach assumes the noise is white, and to be a function of SNR's for each frequency bin implies that the noise model must be generalized. We have approximated the KLT with the DFT, and will now consider applying the signal subspace approach to a whitened version of the noisy speech. Say W is the whitening filter for the noise w. Then, after applying H to the whitened noisy speech Wy we obtain an estimate of Wx. Solving for {circumflex over (x)}, we have

{circumflex over (x)}=W−1HWy  (47)

[0087] where

H=UQU#  (48)

[0088] W=UWFU#  (49)

[0089] Since we are using a DFT approximation to the KLT, U# is the DFT matrix operator and U is the inverse DFT matrix operator. In (49), WF is the frequency domain implementation of the whitening filter. Therefore WF is a diagonal matrix, and Q is diagonal as derived in Section II.B. Substituting (48) and (49) into (47) and simplifying, we obtain 16$\begin{array}{cc}\begin{array}{c}\hat{x}={\mathrm{UW}}_F^{-1}{\mathrm{QW}}_F{U}^{\#}y\\ ={\mathrm{UQU}}^{\#}y\\ =\mathrm{Hy}\end{array}& \left(50\right)\end{array}$

[0090] We have shown that whitening the signal, applying the signal subspace technique, and then applying the inverse of the whitening filter is equivalent to applying the signal subspace technique to the colored noise directly. The constraint, however, is modified. For the whitened noisy input, we now have

E[|u k # {tilde over (r)} w |2]≦αk{tilde over (σ)}w2  (51)

[0091] where

{tilde over (r)}w=HWw  (52)

[0092] {tilde over (σ)}w2=E[|uk#Ww|2]  (53)

[0093] So {tilde over (r)}w given in (52) is the residual whitened noise, and {tilde over (σ)}w2 given in (53) is the variance of this whitened noise. Since, according to the principles of the invention, we are using the DFT approximation to the KLT, the expectations in (51) and (53) are energy spectral density coefficients of the residual whitened noise and the whitened noise respectively. Therefore, dividing the kth constraint given in (51) by the magnitude squared of the kth component of the whitening filter in the frequency domain |WFk|2, we obtain our new constraint:

S r w r w (k)≦αkSww(k)  (54)

[0094] Here Srwrw(k) and Sww(k) are the kth spectral coefficients of the residual noise and original noise, respectively.

[0095] The final step is to choose the constant constraints αk in (54). For white noise, Ephraim found that αk=exp{−νσw2/λx(k)} was a good selection for aggressive noise suppression. For the DFT approximation to the KLT, we have λx(k)=Sxx(k). To extend the technique to colored noise, we have determined to try 17$\begin{array}{cc}\begin{array}{c}{\alpha }_k=\exp \left\{-\upsilon ·S_{\mathrm{ww}}\left(k\right)/S_{\mathrm{xx}}\left(k\right)\right\}\\ =\exp \left\{-\upsilon /{\xi }_k\right\}\end{array}& \left(55\right)\end{array}$

[0096] In (55), we have ensured that the resulting gain depends heavily on the estimate of the a-priori SNR 86k. In this manner, we heavily base our core estimator on the decision-directed estimate of ξk and benefit from the resulting reduction in musical noise.

[0097] A first embodiment of our new core estimator 38 (FIG. 4) for the hybrid speech enhancement system is illustrated in FIG. 5 along with a DFT 44. The first embodiment of the core estimator 38 is coupled to the DFT 44. The DFT 44 receives the noisy signal and converts it into DFT coefficients in the frequency domain. The core estimator 38 includes a set of gains in accordance with (55), which is applied in the frequency domain to the DFT spectral coefficients of the noisy signal. One of the set of gains is applied to each DFT coefficient of the noisy speech by the core estimator 38. The DFT coefficients of the noisy signal are passed from the core estimator 38 to the soft decision module 40 (FIG. 4) for further enhancement.

III.B. Differences with the Modified MMSE-LSA

[0098] The gain that is applied to the noisy signal in the frequency domain in the hybrid speech enhancement system according to the principles of the invention is different than the gain that is applied in the frequency domain according to the modified MMSE-LSA technique developed by Malah.

[0099] In the modified MMSE-LSA approach developed by Malah, we consider clean speech x[n] that has been contaminated with uncorrelated additive noise w[n] to produce noisy speech y[n]:

y[n]=x[n]+w[n]   (56)

[0100] In the frequency domain, we have

Y k =X k +W k   (57)

[0101] where

Xk=AkeJφk  (58)

Yk=RkeJθk  (59)

[0102] We now estimate Ak by minimizing the log-spectral amplitude in a MMSE sense: 18$\begin{array}{cc}{\hat{A}}_k=\underset{B}{\arg \min }E\left[{\left(\log A_k-\log B\right)}^2\right]& \left(60\right)\end{array}$

[0103] so the enhanced signal (in the frequency domain) becomes

{circumflex over (X)} k =Â k e J θ k   (61)

[0104] It turns out that Ak can be computed by simply applying a gain in the frequency domain:

[0105] Â k =G (εk,γk)·Rk  (62)

[0106] where G(εk,γk) is a complicated function of the a-priori and a-posteriori SNR's εk and γk.

[0107] On the other hand, the gain applied in the frequency domain by the hybrid speech enhancement system in accordance with the principles of the invention is closer to that used in the signal subspace approach developed by Ephraim, but is still fundamentally different. We begin in vector notation with

y=x+w   (63)

[0108] and estimate the clean speech by filtering the noisy speech with a linear filter H:

{circumflex over (x)}=Hy  (64)

[0109] We can decompose the residual error into a term solely dependent on the clean speech, called the signal distortion rx, and a term solely dependent on the noise, called the residual noise rw: 19$\begin{array}{cc}\begin{array}{c}r=\hat{x}-x\\ =\left(H-I\right)x+\mathrm{Hw}\\ =r_x+r_w\end{array}& \left(65\right)\end{array}$

[0110] H is chosen so as to minimize the signal distortion energy while keeping the residual noise constrained in the frequency domain:

H=arg min{overscore (ε)}x2 such that Srwrw(k)≦αkSww(k)  (66)

[0111] Here {overscore (ε)}x2=tr E[rxrx#] is the signal distortion energy, Srwrw(k) is the kth spectral coefficient of the residual noise rw, Sww(k) is the kth spectral coefficient of the noise w, and the αk are constants. H turns out to (approximately) apply a gain to each frequency component of the noisy speech:

 k =G k ·R k   (67)

[0112] where

G k ={square root}{square root over (αk)}  (68)

III.C. Modular Structure

[0113] Referring to FIG. 4, the hybrid speech enhancement system includes the core estimator 38 along with the support modules that perform the noise adaptation 34, SNR estimation 36, and soft decision gain calculation 40 tasks. The core estimator 38 of the hybrid speech enhancement system performs a short-time spectral amplitude (STSA) speech enhancement process in the frequency domain by modifying the spectral amplitude of the noisy speech without touching the phase (i.e. using the noisy phase). According to the principles of the invention, the purpose of the core estimator 38 in the hybrid speech enhancement system shown in FIG. 4 is to provide a gain for each frequency bin of the spectral amplitude of the noisy speech. The core estimator 38 is constructed to take advantage of the other modules (for example, by making direct use of the estimated SNR's from the SNR estimator 36).

[0114] The noise adaptation processor 34 segments the noisy speech into noise-only and signal-containing frames, and is responsible for maintaining a current estimate of the noise spectrum as well as an estimate of the probability of signal presence in each frequency bin. These parameters are used when estimating the SNR's, and also impact the core estimator and soft decision gains directly. For example, during a noise-only frame a constant gain is applied to the noise in order to avoid noise structuring.

[0115] Given the noise estimate λw(k), two SNR's are computed. The a-posteriori SNR, γk, is directly measured, while the a-priori SNR, ξk, is estimated using the decision-directed approach.

[0116] A second embodiment of the core estimator 38 (FIG. 4) is illustrated in FIG. 6, along with a DFT 52. The core estimator 38 is coupled to the DFT 52. The DFT 52 receives the noisy speech signal containing an original amount of noise. The DFT 52 transforms the noisy signal containing the original noise into DFT coefficients in the frequency domain. After the noisy signal is transformed into the frequency domain, the core estimator applies a set of gains, Gk={square root}{square root over (αk)}, to the DFT coefficients in the frequency domain and outputs noisy speech containing a residual noise. Here the energy of the signal distortion is minimized with the residual noise constrained by the αk's. We developed a set of constraints for the αk's: 20$\begin{array}{cc}{G}_k=\sqrt{\exp \left(-\upsilon /{\eta }_k\right)},\mathrm{where}{\eta }_k=\frac{{\xi }_k}{1-q_k}& \left(69\right)\end{array}$

[0117] and ν is some constant indicating the level of aggression of the speech enhancement. In the second embodiment of the core estimator 38 depicted in FIG. 6, these gains described by (69) are applied to the DFT coefficients received from the DFT 52. After the core estimator 38 applies the gains to the DFT coefficients of the noisy speech, the noisy signal is passed to the soft decision module 40 (FIG. 4) for further enhancement.

[0118] In the hybrid speech enhancement system, the soft decision module 40 of FIG. 4 operates in the frequency domain to apply a second set of gains to further enhance the noisy signal. For each frequency bin, the soft decision module 40 computes a gain that is applied to the spectral amplitude of the noisy speech in the frequency domain. The gain for each frequency bin is based on the a-posteriori SNR, the a-priori SNR and the probability of signal absence in each frequency bin, qk.

[0119] The hybrid speech enhancement system illustrated by FIGS. 4, 5 and 6 provides the ability to place constraints on the signal distortion or residual noise energy in the frequency domain yielding a greater flexibility than the modified MMSE-LSA approach developed by Malah. Some of the constraints which can be placed include using soft decision rather than removing noise-only subspace, which results in a less artificial sounding noise. More specifically, the power spectral density of the residual noise is constrained to be below a constant proportion of the original noise power spectral density. The constraints are manipulated so as to fit into the decision-directed approach. The gain applied can depend on signal presence uncertainty, or not.

[0120] An important advantage of the hybrid speech enhancement system as compared to the signal subspace approach developed by Ephraim is the improved performance gained from making use of the modified MMSE-LSA framework. The noise adaptation processor, decision-directed SNR estimator, and soft decision module all help in reducing noise distortion and providing a better trade-off between speech distortion and noise reduction than obtainable with the signal subspace approach alone.

[0121] While several particular forms of the invention have been illustrated and described, it will also be apparent that various modifications can be made without departing from the spirit and scope of the invention.

Claims

1. A speech enhancement system, comprising: a noise adaptation module receiving noisy speech, the noisy speech being characterized by spectral coefficients spanning a plurality of frequency bins and containing an original noise, the noise adaptation module segmenting the noisy speech into noise-only frames and signal-containing frames, and the noise adaptation module determining a noise estimate and a probability of signal absence in each frequency bin; a signal-to-noise ratio estimator coupled to the noise adaptation module, the signal-to-noise ratio estimator determining a first signal-to-noise ratio and a second signal-to-noise ratio based on the noise estimate; and a core estimator coupled to the signal-to-noise ratio estimator and receiving the noisy speech, the core estimator applying to the spectral coefficients of the noisy speech a first set of gains in the frequency domain without discarding the noise-only frames to produce speech that contains a residual noise, wherein the first set of gains is determined based, at least in part, on the second signal-to-noise ratio and a level of aggression, and wherein the core estimator is operative to maintain the spectral density of the spectral coefficients of the residual noise below a proportion of the spectral density of the spectral coefficients of the original noise.

2. The system of claim 1, wherein: each one of the first set of gains is also based on the probability of signal absence in each frequency bin.

3. The system of claim 1, wherein: the system modifies the spectral amplitude of the noisy speech without affecting the phase of the noisy speech.

4. The system of claim 1, wherein: during a noise-only frame, a constant gain is applied to the noise in order to avoid noise structuring.

5. The system of claim 1, wherein: the core estimator applies to the spectral coefficients of the noisy speech one of the first set of gains for each frequency bin.

6. The system of claim 1, further comprising: a soft decision module coupled to the signal-to-noise ratio estimator and to the core estimator, the soft decision module applying a second set of gains to the spectral coefficients of the speech that contains a residual noise.

7. The system of claim 6, wherein: the soft decision module determines the second set of gains based on the first signal-to-noise ratio, the second signal-to-noise ratio and the probability of signal absence in each frequency bin.

8. A method for enhancing speech, comprising the steps of: receiving noisy speech, wherein the noisy speech is characterized by spectral coefficients spanning a plurality of frequency bins and contains an original noise; segmenting the speech into noise-only frames and signal-containing frames; determining a noise estimate and a probability of signal absence in each frequency bin; determining a first signal-to-noise ratio and a second signal-to-noise ratio based on the noise estimate; determining a first set of gains based, at least in part, on the second signal-to-noise ratio and a level of aggression; and applying the first set of gains to the spectral coefficients of the noisy speech without discarding the noise-only frames to produce speech that contains a residual amount of noise, such that the spectral density of the spectral coefficients of the residual noise is maintained below a proportion of the spectral density of the spectral coefficients of the original noise.

9. The method of claim 8, wherein: the first set of gains is also based on the probability of signal absence in each frequency bin.

10. The method of claim 8, further comprising the step of: modifying the spectral coefficients of the noisy speech without affecting the phase of the noisy speech.

11. The method of claim 8, further comprising the step of: during a noise-only frame, applying a constant gain to the noise.

12. The method of claim 8, wherein: one of the first set of gains is applied to the spectral coefficients of the noisy speech for each frequency bin.

13. The method of claim 8, further comprising the step of: applying a second set of gains to the spectral coefficients of the speech that contains a residual noise.

14. The method of claim 13, further comprising the step of: determining the second set of gains based on the first signal-to-noise ratio, the second signal-to-noise ratio and the probability of signal absence in each frequency bin.