CA2169424C - Method and apparatus for noise reduction by filtering based on a maximum signal-to-noise ratio and an estimated noise level - Google Patents

Method and apparatus for noise reduction by filtering based on a maximum signal-to-noise ratio and an estimated noise level Download PDF

Info

Publication number
CA2169424C
CA2169424C CA002169424A CA2169424A CA2169424C CA 2169424 C CA2169424 C CA 2169424C CA 002169424 A CA002169424 A CA 002169424A CA 2169424 A CA2169424 A CA 2169424A CA 2169424 C CA2169424 C CA 2169424C
Authority
CA
Canada
Prior art keywords
noise
frame
value
level
signal
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Lifetime
Application number
CA002169424A
Other languages
French (fr)
Other versions
CA2169424A1 (en
Inventor
Joseph Chan
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sony Corp
Original Assignee
Sony Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Sony Corp filed Critical Sony Corp
Publication of CA2169424A1 publication Critical patent/CA2169424A1/en
Application granted granted Critical
Publication of CA2169424C publication Critical patent/CA2169424C/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10LSPEECH ANALYSIS OR SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING; SPEECH OR AUDIO CODING OR DECODING
    • G10L21/00Processing of the speech or voice signal to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
    • G10L21/02Speech enhancement, e.g. noise reduction or echo cancellation
    • G10L21/0208Noise filtering
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10LSPEECH ANALYSIS OR SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING; SPEECH OR AUDIO CODING OR DECODING
    • G10L21/00Processing of the speech or voice signal to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
    • G10L21/02Speech enhancement, e.g. noise reduction or echo cancellation
    • G10L21/0208Noise filtering
    • G10L21/0216Noise filtering characterised by the method used for estimating noise
    • G10L2021/02161Number of inputs available containing the signal or the noise to be suppressed
    • G10L2021/02163Only one microphone
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10LSPEECH ANALYSIS OR SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING; SPEECH OR AUDIO CODING OR DECODING
    • G10L21/00Processing of the speech or voice signal to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
    • G10L21/02Speech enhancement, e.g. noise reduction or echo cancellation
    • G10L21/0208Noise filtering
    • G10L21/0216Noise filtering characterised by the method used for estimating noise
    • G10L21/0232Processing in the frequency domain
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10LSPEECH ANALYSIS OR SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING; SPEECH OR AUDIO CODING OR DECODING
    • G10L25/00Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
    • G10L25/27Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the analysis technique

Abstract

A method for reducing the noise in an speech signal by removing the noise from an input speech signal is disclosed. The noise reducing method includes converting the input speech signal into a frequency spectrum, determining filter characteristics based upon a first value obtained on the basis of the ratio of a level of the frequency spectrum to an estimated level of the noise spectrum contained in the frequency spectrum and a second value as found from the maximum value of the ratio of the frame- based signal level of the frequency spectrum to the estimated noise level and the estimated noise level, and reducing the noise in the input speech signal by filtering responsive to the filter characteristics. A corresponding apparatus for reducing the noise is also disclosed.

Description

TITLE OF THE INVENTION

Method and Apparatus for Noise Reduction by Filtering Based on a Maximum Signal-to-Noise Ratio and an Estimated Noise Level BACKGROUND OF THE INVENTION

This invention relates to a method for removing the noise contained in a speech signal for suppressing or reducing the noise therein.

In the fi el d of a portable telephone set or speech recognition, it is felt to be necessary to suppress the noise such as background noise or environmental noise contained in the collected speech signal for emphasizing its speech components.

As a technique for emphasizing the speech or reducing the noise, a technique of employing a conditional probability function for attenuation factor adjustment is disclosed in R.J. McAulay and M.L. Maplass, "Speech Enhancement using a Soft-Decision noise Suppression Filter, in IEEE
Trans. Acoust., Speech Signal Processing, Vol.28, pp.137 to 145, April 1980.

In the above noise-suppression technique, it is a frequent occurrence that unspontaneous sound tone or distorted speech be produced due to an inappropriate suppression filter or an operation based upon an inappropriate fixed signal-to-noise ratio (sNR). It is not desirable for the user to have to adjust the SNR, as one of the parameters of a noise suppression device, in actual operation for realizing an optimum performance. In addition, it is difficult with the conventional speech signal enhancement technique to eliminate the noise sufficiently without generating_ distortion in the speech signal susceptible to significant variation in the SNR in short time.

Such speech enhancement or noise reducing technique employs a technique of discriminating a noise domain by comparing the input power or level to a pre-set threshold value. However, if the time constant of the threshold value is increased with this technique for prohibiting the threshold value from tracking the speech, a changing noise level, especially an increasing noise level, cannot be followed appropriately, thus leading occasionally to mistaken discrimination.

For overcoming this drawback, the present inventors have proposed in JP Patent Application Hei-6-99869 (1994) a noise reducing method for reducing the noise in a speech signal.

With this noise reducing method for the speech si gnal , noise suppression is achieved by adaptively controlling a maximum likelihood filter configured for calculating a speech component based upon the SNR derived from the input speech signal and the speech presence probability. This method employs a signal corresponding to the input speech spectrum less the estimated noise spectrum in calculating the speech presence probability.

With this noise reducing method for the speech si gnal , since the maximum likelihood filter is adjusted to an optimum suppression filter depending upon the SNR of the input speech signal, sufficient noise reduction for the input speech signal may be achieved.
However, since complex and voluminous processing operations are required for calculating the speech presence probability, it has been desired to simplify the processing operations.
SUMMARY OF THE INVENTION

It is therefore an object of the present invention to provide a noise reducing method for an input speech signal whereby the processing operations for noise suppression for the input speech signal may be simplified.

In one aspect, the present invention provides a method for reducing the noise in an input speech signal for noise suppression including converting the input speech signal into a frequency spectrum, determining filter characteristics based upon a first value obtained on the basis of the ratio of a level of the frequency spectrum to an estimated level of the noise spectrum contained in the frequency spectrum and a second value as found from the maximum value of the ratio of the frame-based signal level of the frequency spectrum to the estimated noise level and from the estimated noise level, and reducing the noise in the input speech signal by filtering responsive to the filter characteristics.

In another aspect, the present invention provides an apparatus for reducing the noise in an input speech signal for noise suppression including means for converting the input speech signal into a frequency spectrum, means for determining filter characteristics based upon a first value obtained on the basis af the rati_o of a level of the frequency spectrum to an estimated level of the noise spectrum contained in the frequency spectrum and a second value as found from the maximum value of the ratio of the frame-based signal level of the frequency spectrum to the estimated noise level and from the estimated noise level, and means for reducing the noise in the input speech signal by filtering responsive to the filter characteristics.

With the method and apparatus for reducing the noise in the speech signal, according to the present invention, the first value is a value calculated on the basis of the ratio of the input signal spectrum obtained by transform from the input speech signal to the estimated noise spectrum contained in the input signal spectrum, and sets an initial value of filter characteristics determining the noise reduction amount in the filtering for noise reduction. The second value is a value calculated on the basis of the maximum value of the ratio of the signal level of the input signa spectrum to the estimated noise level, that is the maximum SNR, and the estimated noise level, and is a value for variably controlling the filter characteristics. The noise may be removed in an amount corresponding to the maximum SNR from the input speech signal by the filtering conforming to the filter characteristics variably controlled by the first and second values.

Since a table having pre-set levels of the input signal spectrum and the estimated levels of the noise spectrum entered -therein may be used for finding the first value, the processing volume may be advantageously reduced.

Also, the second value is obtained responsive to the maximum SNR and the frame-based noise level, the filter characteristics may be adjusted so that the maximum noise reduction amount by the filtering will be changed substantially linearly in a dB area responsive to the maximum SN ratio.

With the above-described noise reducing method of the present invention, the first and the second value are used for controlling the filter characteristics for filtering for removing the noise from the input speech signal, whereby the noise may be removed from the input speech signal by filtering conforming to the maximum SNR in the input speech signal, in particular, the distortion in the speech signal caused by the filtering at the high SN ratio may be diminished and the volume of the processing operations for achieving the filter characteristics may also be reduced.

In addition, according to the present invention, the first value for controlling the filter characteristics may be calculated using a table having the levels of the input signal spectrum and the levels of the estimated noise spectrum entered therein for reducing the processing volume for achieving the filter characteristics.

Also, according to the present invention, the second value obtained responsive to the maximum SN ratio and to the frame-based noise level may be used for controlling the f i l te r characte ri sti cs for reducing the processing volume for a c h i e v i ng the fi 1 te r characte ri sti cs . The maximum n o i s e r e d u c t i o n amount achieved by the f i 1 te r characteri sti cs may be changed r e s p o n s i v e to the N r a t i o of the input speech signal.

BRIEF DESCRIPTION OF THE DRAWINGS

Fig.1 illustrates a first embodiment of the noise reducing method for the speech signal of the present invention, as applied to a noise reducing apparatus.

Fig.2 illustrates a specific example of the energy E[k] and the decay energy EdecaY[k] in the embodiment of Fig.l.

Fig.3 illustrates specific examples of an RMS value RMS[k], an estimated noise level value MinRMS[k] and a maximum RMS value MaxRMS[k] in the embodiment of Fig.l.

Fig.4 illustrates specific examples of the relative energy Brel[k], a maximum SNR MaxSNR[k] in dB, and a value dBth resre, [ k ], as one of th reshol d values for noise discrimination, in the embodiment shown in Fig.l.

Fig. 5 is a graph showing NR_ level [ k] as a function defined with respect to the maximum SNR
MaxSNR[k], in the embodiment shown in Fig.l.

Fig.6 shows the relation between NR[w,k] and the maximum noise reduction amount in dB, in the embodiment shown in Fig.l.

Fig.7 shows the relation between the ratio of Y[w,k]/N[w,k]
and Hn[w,k] responsive to NR[w,k] in dB, in the embodiment siaown in F i g. 1.

Fig.8 illustrates a second embodiment of the noise reducing method for the speech signal of the present invention, as applied to a noise reducing apparatus.

Figs. 9 and 10 are graphs showing the distortion of segment portions of the speech signal obtained on noise suppression by the noise reducing apparatus of Figs.1 and 8 with respect to the SN ratio of the segment portions.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

Referring to the drawings, a method and apparatus for reducing the noise in the speech signal according to the present invention will be explained in detail.

Fig.1 shows an embodiment of a noise reducing apparatus for reducing the noise in a speech signal according to the present invention.

The noise reducing apparatus includes, as main components, a fast Fourier transform unit 3 for converting the input speech signal into a frequency domain signal or frequency spectra, an Hn value calculation unit 7 for controlling filter characteristics during removing the noise portion from the input speech signal by filtering, and a spectrum correction unit 10 for reducing the noise in the input speech signal by filtering responsive to filtering characteristics produced by the Hn value calculation unit 7.

An input speech signal y[t], entering a speech signal input terminal 13 of the noise reducing apparatus, is provided to a f rami ng unit 1. A framed signal y_ f ramej k outputted by the framing unit 1, is provided to a wi_ndowing unit 2, a root mean square (RMS) calculationunit2lwithin a noise estimation unit 5, and a filtering unit 8.

An output of the windowing unit 2 is provided to the fast fourier transform unit 3, an output of which is provided to both the spectrum correction unit 10 and a band-splitting unit 4. An output of the band-splitting unit 4 is provided to the spectrum correction unit 10, a noise spectrum estimation unit 26 within the noise estimation unit 5 and to the Hn value calculation unit 7. An output of the spectrum correction unit 10 is provided to a speech signal output terminal 14 via the fast Fourier transform unit 11. and an overlap-and-add unit 12.

An output of the RMS calculation unit 21 is provided to a relative energy calculation unit 22, a maximum RMS calculation unit 23, an estimated noise level calculation unit 24 and to a noise spectrum estimation unit 26. An output of the maximum RMS
calculation unit 23 is provided to an estimated noise level calculation unit 24 and to a maximum S.NR calculation unit 25. An output of the relative energy calculation unit 22 is provided to a noise spectrum estimation unit 26. An output of the estimated noise level calculation unit 24 is provided to the filtering unit 8, maximum SNR calculation unit 25, noise spectrum estimation unit 26 and to the NR value calculation unit 6. An output of the maximum SRR calculation unit 25 is provided to the NR value calculation unit 6 and to the noise spectrum estimation unit 26, an output of which is provided to the Hn value calculation unit 7.

An output of the NR value calculation unit 6 is again provided to the NR value calculation unit 6, while being also provided to the Hn value calculation unit 7.

An output of the Hn value calculation unit 7 is provided via the filtering unit 8 and a band conversion unit 9 to the spectrum correction unit 10.

The operation of the above-described first embodiment of the noise reducing apparatus is explained.

To the speech signal input terminal 13 is supplied an input speech signal y[t] containing a speech component and a noise component. The input speech signal y[t], which is a digital signal sample at, for example, a sampling frequency FS, is provided to the framing unit 1 where it is split into plural frames each having a frame length of FL sampl es . The input speech signal y[t], thus split, is then processed on the frame basis.
The frame interval, which is an amount of displacement of the frame along the time axis, is FI samples, so that the (k+1)st rame begins after FI samples as from the k'th frame. By way of illustrative examples of the sampling frequency and the number f samples, if the sampling frequency FS is 8 kHz, the frame interval FI of 80 samples corresponds to 10 ms, while the frame length FL of 160 samples corresponds to 20 ms.

Prior to orthogonal transform calculations by the fast Fourier transform unit 2, the windowing unit 2 multiplies each framed signal y_f ramejk from the f rami ng unit 1 with a wi ndowi ng function winput= Following the inverse FFI, performed at the terminal stage of the frame-based signal processing operations, as will be explained later, an output signal is multiplied with a wi ndowi ng funct i on woutput = The wi ndowi ng funct i ons winput and woutput may be respect i vel y exempl i f i ed by the fol l owi ng equat i ons (1) and (2):

Wjnput (J) =(1 -1 cos ( 27tJ )) c, FL 0 s j s FL

...(1) Wou.'aut [J] =( 1 - 4.cos( 27[FLj ))~, 0 S j 5 FL

(2) The fast Fourier transform unit 3 then performs 256-point fast Fourier transform operations to produce frequency spectral amplitude values, which then are split by the band splitting unit 4 into, for example, 18 bands. The frequency ranges of these bands are shown as an example in Table 1:

band numbers frequency ranges 0 0 to 125 Hz 1 125 to 250 Hz 2 250 to 275 Hz 3 375 to 563 Hz 4 563 to 750 Hz 750 to 938 Hz 6 938 to 1125 Hz 7 1125 to 1313 Hz 8 1313 to 1563 Hz 9 1563 to 1813 Hz 1813 to 2063 Hz 11 2063 to 2313 Hz 12 2313 to 2563 Hz 13 2563 to 2813 Hz 14 2813 to 3063 hz 3063 to 3375 hz 16 3375 to 3688 Hz 17 3688 to 4000 Hz The amplitude values of the frequency bands, resulting from frequency spectrum splitting, become amplitudes Y[w,k] of the input signal spectrum, which are outputted to respective portions, as explained previously.

The above frequency ranges are based upon the fact that the higher the frequency, the less becomes the perceptual resolution of the human hearing mechanism. As the amplitudes of the respective bands, the maximum FFT amplitudes in the pertinent frequency ranges are employed.

In the noise estimation unit 5, the noise of the framed signal y_framejk is separated from the speech and a frame presumed to be noisy is detected, while the estimated noise level value and the maximum SN ratio are provided to the NR value I

calculatio-n unit 6. The noisy domain estimation or the noisy frame detection is performed by combination of, for example, three detection operations. An illustrative example of the noisy domain estimation is now explained.

The RMS calculation unit 21 calculates RMS values of signals every frame and outputs the calculated RMS values. The RMS value of the k'th frame, or RMS[k], is calculated by the following equation (3):

RMS[k] ,E (y_framei,k) z FLj=o (3) In the relative energy calculation unit 22, the relative energy of the k'th frame pertinent to the decay energy from the previous frame, or dBfel[k], is calculated, and the resulting value is outputted. The relative energy in dB, that is dBfel[k], is found by the following equation (4):

dBzel [k] = 101ogiol Edecay [k] }
\ E[k] 1 (4) while the energy value E[k] and the decay energy value E
decay [ k ]
are found from the following equations (5) and (6):

FL
E[k] _ E (Y-frameJ, k) Z
1=~

...(5) ~
Edecay [k] = max(E [k] , (exp -FI ~ *Edecay [k - 11 1 0.65*FS

(6) The equation (5) may be expressed from the equation (3) as FL*(RMS[k])2. Of course, the value of the equation (5), obtained during calculations of the equation (3) by the RMS calculation unit 21, may be directly provided to the relative energy calculation unit 21. In the equation (6), the decay time is set to 0.65 second.

Fig.2 shows illustrative examples of the energy value E[k]
and the decay energy Edecay[k] =

The maximum RMS calculation unit 23 finds and outputs a maximum RMS value necessary for estimating the maximum value of the ratio of the signal level to the noise level, that is the maximum SN ratio. This maximum RMS value MaxRMS[k] may be found by the equation (7):

MaxRMS[k] = max(4000,RMS[k] ,6*MaxRMS[k-1] +(1-6) *RMS[k) ) ...(7) where 0 is a decay constant. For 0, such a value for which the maximum RMS value is decayed by 1/e at 3.2 seconds, that is 0 = 0.993769, is employed.

The estimated noise level calculation unit 24 finds and outputs a minimum RMS value suited for evaluating the background noise level. This estimated noise level value minRMS[k] is the smallest value of five local minimum values previous to the current time point, that is five values satisfying the equation (8):

(RMS[k] < 0.6*MaxRMS[k] and RMS[k] < 4000 and RMS[k] < RMS[k+1] and RMS[k] < RMS[k-1] and RMS[k] < RMS[k-2]) or (RMS[k] < MinRMS) (8) The estimated noise level value minRMS[k] is set so as to r i se for the background noise freed of speech. The r i se rate for the high noise level is exponential, while a fixed rise rate is used for the low noise level for +-ealizing a more outstanding rise.

Fig.3 shows illustrative examples of the RMS values RMS[k], estimated noise level value minRMS[k] and the maximum RMS values MaxRMS[k].

The maximum SNR calculation unit 25 estimates and calculates the maximum SN ratio MaxSNR[k], using the maximum RMS value and the estimated noise level value, by the following equation (9);
MaxSNR[k] = 201og (MaxRMS[k] )-loMi.rnRMS [k]
... 1 (9) From the maximum SNR value MaxSNR, a normalization parameter NR level i.n a range from 0 to 1, representing the relative noise level, is calculated. For NR_level, the following function is employed:

( 2+ 2 cosTt MaxSN20[k] -30x (1-0. 002 (MaxSNR [k] -30) Z;
N1~1eve1 [k] = 30<MaxSNR[k] S50 0.0 MaxSNR[k] >50 1.0 MaxSNR [k] : o the ...(10) The operation of the noise spectrum estimation unit 26 is explained. The respective values found in the relative energy calculation unit 22, estimated noise level calculation unit 24 and the maximum SNR calculation unit 25 are used for discriminating the speech from the background noise. If the following conditions:

((RMS[k] < NoiseRMSthres[k]) or ( dBrel [ k ] > dBthres [ k l)) and (RMS[k] < RMS[k-1]+200) .(11) where NoiseRMSthres[kl = 1.05+0.45*NR_level[k]xMinRMS[k]
dBthres rel [k] = max(MaxSNR[k]-4.0, 0.9*MaxSNR[k]

are valid, the signal in the k'th frame is classified as the background noise. The amplitude of the ba.ckground noise, thus classified, is calculated and outputted as a time averaged estimated value N[w,k] of the noise spectrum.

2 i 69424 Fig.4. shows illustrative examples of the relative energy in dB, shown in Fig.11, that is dBfel[k], the maximum SNR[k] and dBthresrel , as one of the threshold values for noise discrimination.

Fig.6 shows NR_level[k], as a function of MaxSNR[k] in the equation (10).

If the k'th frame is classified as the background noise or as the noise, the time averaged estimated value of the noise spectrum N[w,k] is updated by the amplitude Y[w,k] of the input signal spectrum of the signal of the current frame by the following equation (12):

N[w,k] = a*max(N[w,k-1], Y[w,k]) + (1 - a)*min(N[w,k-1], Y[w,k]) (12) -a =exp FI
(0.5*FS
where w specifies the band number in the band splitting.

If the k'th frame is classified as the speech, the value of N[w,k-1] is directly used for N[w,k].

The NR value calculation unit 6 calculates NR[w,k], which is a value used for prohibiting the filter response from being changed abruptly, and outputs the produced value NR[w,k]. This NR [w, k] i s a val ue rangi ng f rom 0 to 1 and i s def i ned by the equation (13):

(13) adj [w, k] NR [w, k-1] -8NR<adj [ w, k] <NR [ w, k-1J +8NR
NR[w,k] = NR[w,k-1] -bl,,R NR[w,k-1] -Bj,,R>-adi [w,k]
NR[w,k-1] +8NR NR[w,k-1] +81,Rsadj [w,k]
SNR = 0.004 adj[w,k] = min(adj1[k],adj2[k])-adj3[w,k]

In the equation (13), adj[w,k] is a parameter used for taking into account the effect as explained below and is defined by the equation (14):

bNR = 0.004 and adj[w,k] = min(adjl[k],adj2[k])-adj3[w,k] ...(14) In the equation (14), adjl[k] is a value having the effect of suppressing the noise suppressing effect by the filtering at the high SNR by the filtering described below, and is defined by thr-: following equation (15):

1 MaxSNR [ k] < 2 9 adjl [k] = 1 ',1axSNR [k] -29 29 sMaxSNR [k] <43 0 MaxSNR [k] : o therwi se ...(15) In the equation (14), adj2[k] is a value having the effect of suppressing the noise suppression rate with respect to an extremely low noise level or an extremely high noise level, by the above-described filtering operation, and is defined by the following equation (16):

(16) In the above equation (14), adj3[k] is a value having the 0 MinRMS[k] <20 MinRMS[k] -20 20sMinRMS[k] <60 adj2[k] = 1 60sMinRMS[k] <1000 1_ MinRMS[k] -1000 1000sMinRMS[k] <1800 0.2 1800sMinRMS[k]
effect of suppressing the maximum noise reduction amount from 18 dB to 15 dB between 2375 Hz and 4000 Hz, and is defined by the following equation (17):

0 w<2375Hz adj3 [w, k] = 0. 059415 (w-2375) w: otherwise (17) Meanwhile, it is seen that the relation between the above values of NR[w,k] and the maximum noise reduction amount in dB
is substantially linear in the dB region, as shown in Fig.6.

The Hn value calculation unit 7 generates, from the amplitude Y[w,k] of the input signal spectrum, split into frequency bands, the time averaged estimated value of the noise spectrum N[w,k] and the value NR[w,k], a value Hn[w,k] which determines filter characteristics configured for removing the noise portion from the input speech signal. The value Hn[w,k]
is calculated based upon the following equation (18):

Hn[w,k] = 1-(2*NR[w,k]-NR2[w,k])*(1-H[w][S/N=y]) (18) The value H[w][S/N=r] in the above equation (18) is equivalent to optimum characteristics of a noise suppression filter when the SNR is fixed at a value r, and is found by the following equation (19):

H[w] [S/N=y] = 2 (1+1_2[1kJJ*p(H1:yw) [s/N=r]+G.:n*POi Yw) [s/r7=T]
' (19) Meanwhile, this value may be found previously and listed in a table in accordance with the value of Y[w,k]/N[w.k].
Meanwhile, x[w,k] in the equation (19) is equivalent to Y[w,k]/N
[w,k], while Gmin is a parameter indicating the minimum gain of H[w][S/N=r]. On the other hand, P(HijYW)[S/N=r] and p(HOjYW[S/N
=r] are parameters specifying the states of the amplitude Y[w, k] while P(H1lYW)[S/N=r] is a parameter specifying the state in which the speech component and the noise component are mixed together in Y[w,k] and P(HOIYW)[S/N=r] is a parameter specifying that only the noise component is contained in Y[w,k]. These values are calculated in accordance with the equation (20):

P(H1 i Yw) [s/tv-Y] = 1-P(HO1 Yw) [s/lv=tr]

P P(H1) * (exp (-y2) ) *10 (2*Y*x[w.k] ) +P(HO) * (exp (-x2 g{ {. . . (20) where P(hl) = P(HO) = 0.5 It is seen from the equation (20) that P(H1I YW)[S/N=r] and P(HOIYW)[S/N=r] are functions of x[w,k], while Io(2*r*x [w,k]) is a Bessel function and is found responsive to the values of r and 2ib9424 1w,k]. Both P(H1) and P(HO) are fixed at 0.5. The processing volume may be reduced to approximately one-fifth of that with the conventional method by simplifying the parameters as described above.

The relation between the Hn[w,k] value produced by the Hn value calculation unit 7, and the x[w,k] value, that is the ratio Y[w,k]/N[w,k], is such that, for a higher value of the ratio Y
[w,k]/N[w,k], that is for the speech component being higher than the noisy component, the value Hn[w,k] is increased, that is the suppression is weakened, whereas, for a lower value of the ratio Y[w,k]/N[w,k], that is for the speech component being lower than the noisy component, the value Hn[w,k] is decreased, that is the suppression is intensified. In the above equation, a solid line curve stands for the case of r = 2.7, Gmin =-18 dB and NR[w,k]
= 1. It is also seen that the curve specifying the above relation is changed within a range L depending upon the NR[w,k]
value and that respective curves for the value of NR[w,k] are changed with the same tendency as for NR[w,k] = 1.

The filtering unit 8 performs filtering for smoothing the Hn[w,k] along both the frequency axis and the time axis, so that a smoothed signal Ht_smooth[w,k] is produced as an output signal.
The filtering in a direction along the frequency axis has the effect of reducing the effective impulse response length of the signal Hn[w,k]. This prohibits the aliasing from being produced due to cyclic convolution resulting from realization of a filter by multiplication in the frequency domain. The filtering in a direction along the time axis has the effect of limiting the rate of change in filter characteristics in suppressing abrupt noise generation.

The filtering in the direction along the frequency axis is first explained. Median filtering is performed on Hn[w,k] of each band. This method is shown by the following equations (21) and (22):

step 1: H1[w,k] = max(median(Hn[w-i,k], Hn[w,k]
,Hn[w+1,k],Hn[w,k]) ...(21) step 2: H2[w,k] = min(median(H1[w-i,k],H1[w,k]

,H1[w+1,k],H1[w,k]) ...(22) If, in the equations (21) and (22), (w-1) or (w+1) is not present, H1[w,k] = Hn[w,k] and H2[w,k] = H1[w,k], respectively In the step 1, H1[w,k] is Hn[w,k] devoid of a sole or lone zero (0) band, whereas, in the 2, H2[w,k] H1[w,k] devoid of a sole, lone or protruding band. In this manner, Hn[w,k] is converted into H2[w,k].

Next, filtering in a direction along the time axis is explained. For filtering in a direction along the time axis, the fact that the input signal contains three components, namely the speech, background noise and the transient state representing the transient state of the rising portion of the speech, is taken into account . The speech signal Hspeech[w, k] is smoothed along the time axis, as shown by the equation (23):

Hspeech[w,k] = 0.7*H2[w,k]+0.3*H2[w,k-1] ...(23) The background noise is smoothed in a direction along the axis as shown in the equation (24):

Hnoise[w, k]= 0.7*Min_H+0.3*Max_H ,,.(24) In the above equation (24), Min_H and Max_H may be found by Min H= min(H2[w,k], H2[w,k-1]) and Max_H = max(H2[w,k],H2[w,k-1]), respectively.

The signals in the transient state are not smoothed in the direction along the time axis.

Using the above-described smoothed signals, a smoothed output signal Ht smooth is produced by the equat i on (25) :
Ht-smooth [ w , k] _ ( 1-atr ) (asp*Hspeech[w, k]

+ (1-asp)*Hnoise[w,k])+atr*H2[w,k]

...(25) In the above equation (25), asp and atr may be respectively found from the equation (26):

1.0 S'NRinst>4 . 0 C[sp = 3 (SNRinst-1) 1.0<SNRinsc<4.0 0 SNRirst : otherwise ...(26) where SNR t - RMS[k]
ins - Mi12RMS[k-1]
and from the equation (27):

(27) 1.0 a rms> 3. 5 aBp = 3 (8rms-2) 1.0<8rms<3.5 0 8rmB : o terwi se where RMSlocai [k]
bzma - RMSlocal [k-1]

Then, at the band conversion unit 9, the smoothing signal Ht smooth[w,k] for 18 bands from the filtering unit 8 is expanded 141KS1oca1 [k] = 1 E(y_frameJ,k) 2 FI FI
.i= Z

by interpolation to, for example, a 128-band signal H128[w,k], which is outputted. This conversion s performed by, for example, two stages, whi l e the expansion from 18 to 64 bands and that from 64 bands to 128 bands are performed by zero-order holding and by low pass filter type interpolation, respectively.

The spectrum correction unit 10 then multiplies the real and imaginary parts of FFT coefficients obtained by fast Fourier transform of the framed signal y_ frame j, k obtained by FFT unit 3 with the above signal H128[w,k] by way of performing spectrum correction, that is noise component reduction. The resulting signal is outputted. The result is that the spectral amplitudes are corrected without changes in phase.

The i.nverse FFT unit 11 then performs inverse FFT on the output signal of the spectrum correction unit 10 in order to output the resultant IFFTed signal.

The overlap-and-add unit 12 overlaps and adds the frame boundary portions of the frame-based IFFted signals. The resulting output speech signals are outputted at a speech signal output terminal 14.

Fig.8 shows another embodiment of a noise reduction apparatus for carrying out the noise reducing method for a speech signal according to the present invention. The parts or components which are used in common with the noise reduction apparatus shown in Fig.1 are represented by the same numerals and the description of the operation is omitted for simplicity.

The noise reduction appara'_.:s has a fast Fourier transform unit 3 for t ransformi ng the input speech signal into a frequency-domain signal, an Hn value calculation unit 7 for controlling filter characteristics of the filtering operation of removing the noise component from the input speech signal, and a spectrum correction unit 10 for reducing the noise in the input speech signal by the filtering operation conforming to filter characteristics obtained by the Hn value calculation unit 7.

In the noise suppression filter characteristic generating unit 35, having the Hn calculation unit 7, the band splitting portion 4 splits the amplitude of the frequency spectrum outputted from the FFT unit 3 into, for example, 18 bands, and outputs the band-based amplitude Y[w,k] to a calculation unit 31 for calculating the RMS, estimated noise level and the maximum SNR, a noise spectrum estimating unit 26 and to an initial filter response calculation unit 33.

The calculation unit 31 calculates, from y_framej k, outputted from the f rami ng unit 1 and Y[w,k] outputted by the band splitting unit 4, the frame-based RMS value RMS[k], an estimated noise level value MinRMS[k] and a maximum RMS value Max [k], and transmits these values to the noise spectrum estimating unit 26 and an adj1, adj2 and adj3 calculation unit 32.

The initial filter response calculation unit 33 provides the time-averaged noise value N[w,k] outputted from the noise spectrum estimation unit 26 and Y[w,k] outputted from the band splitting unit 4 to a filter suppression curvcs'table unit 34 -F*or finding out the value of H[w,k] corresponding to Y[w,k] and N [w, k] stored in the filter suppression curve table unit 34 to transmit the value thus found to the Hn value calculation unit 7. In the filter suppression curve table unit 34 is stored a table for H[w,k] values.

The output speech signals obtained by the noise reduction apparatus shown in Figs.1 and 8 are provided to a signal processing circuit, such as a variety of encoding circuits for a portable telephone set or to a speech recognition apparatus.
Alternatively, the noise suppression may be performed on a decoder output signal of the portable telephone set.

2} 69424 Figs.9 and 10 illustrate the distortion in the speech signals obtained on noise suppression by the noise reduction method of the present invention, shown in black color, and the distortion in the speech signals obtained on noise suppression by the conventional noise reduction method , shown in white color, respectively. In the graph of Fig.9, the SNR values of segments sampled every 20 ms are plotted against the distortion for these segments. In the graph of Fig.10, the SNR values for the segments are plotted against distortion of the entire input speech signal. In Figs.9 and 10, the ordinate stands for distortion which becomes smaller with the height from the origin, while the abscissa stands for the SN ratio of the segments which becomes higher towards right.

It is seen from these figures that, as compared to the speech signals obtained by noise suppression by the conventional noise reducing method, the speech signal obtained on noise suppression by the noise reducing method of the present invention undergoes distortion to a lesser extent especially at a high SNR
value exceeding 20.

Claims (4)

1. A method for reducing noise in an input speech signal, the method comprising the steps of:
converting the input speech signal into a frequency spectrum; determining filter characteristics by obtaining a first value representing a ratio of a signal level of the frequency spectrum to an estimated noise level of a noise spectrum contained in the frequency spectrum from a table containing a plurality of pre-set signal levels of the frequency spectrum of the input speech signal and a plurality of pre-set estimated noise levels of the noise spectrum in order to determine an initial value of the filter characteristics, and obtaining a second value representing a maximum value of a ratio of a frame-based signal level of the frequency spectrum to a frame-based estimated noise level and the frame-based estimated noise level for variably controlling the filter characteristics; and reducing noise in the input speech signal by noise filtering using the determined filter characteristics, including decreasing the noise filtering when the frame-based signal level is greater than the frame-based estimated noise level, and increasing the noise filtering when the frame-based signal level is less than the frame-based estimated noise level.
2. The method for noise reduction as claimed in claim 1, wherein the step of obtaining the second value includes obtaining a value by adjusting a maximum noise reduction amount by noise filtering based on the determined filter characteristics so that a maximum noise reduction amount changes substantially linearly in a dB domain.
3. The method for noise reduction as claimed in claim 1, further comprising the steps of:
obtaining the frame-based estimated noise level based on a root mean square value of an amplitude of the frame-based signal level and a maximum value of root mean square values; and calculating the maximum value of the ratio of the frame-based signal level to the frame-based estimated noise level based on the maximum value of the root mean square values and the frame-based estimated noise level, wherein the maximum value of the root mean square values is a maximum value among root mean square values of amplitudes of the frame-based signal level and a value obtained based on the maximum value of the mean root mean square values of a directly previous frame and a pre-set value.
4. An apparatus for reducing noise in an input speech signal and for performing noise suppression, the apparatus comprising:
means for converting the input speech signal into a frequency spectrum;
means for determining filter characteristics based upon a first value representing a ratio of a signal level of the frequency spectrum to an estimated noise level of a noise spectrum contained in the frequency spectrum obtained from a table containing a plurality of pre-set signal levels of the frequency spectrum of the input speech signal and a plurality of pre-set estimated noise levels of the noise spectrum in order to determine an initial value of the filter characteristics, and a second value representing a maximum value of a ratio of a frame-based signal level of the frequency spectrum to a frame-based estimated noise level of the noise spectrum and the frame-based estimated noise level of the noise spectrum for variably controlling the filter characteristics; and means for reducing noise in the input speech signal by noise filtering responsive to the determined filter characteristics, wherein the noise filtering is decreased when the frame-based signal level is greater than the frame-based estimated noise level, and the noise filtering is increased when the frame-based signal level is less than the frame-based estimated noise level.
CA002169424A 1995-02-17 1996-02-13 Method and apparatus for noise reduction by filtering based on a maximum signal-to-noise ratio and an estimated noise level Expired - Lifetime CA2169424C (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP02933695A JP3484801B2 (en) 1995-02-17 1995-02-17 Method and apparatus for reducing noise of audio signal
JPP07-029336 1995-02-17

Publications (2)

Publication Number Publication Date
CA2169424A1 CA2169424A1 (en) 1996-08-18
CA2169424C true CA2169424C (en) 2007-07-10

Family

ID=12273403

Family Applications (1)

Application Number Title Priority Date Filing Date
CA002169424A Expired - Lifetime CA2169424C (en) 1995-02-17 1996-02-13 Method and apparatus for noise reduction by filtering based on a maximum signal-to-noise ratio and an estimated noise level

Country Status (17)

Country Link
US (1) US6032114A (en)
EP (1) EP0727769B1 (en)
JP (1) JP3484801B2 (en)
KR (1) KR100414841B1 (en)
CN (1) CN1140869A (en)
AT (1) ATE209389T1 (en)
AU (1) AU696187B2 (en)
BR (1) BR9600761A (en)
CA (1) CA2169424C (en)
DE (1) DE69617069T2 (en)
ES (1) ES2163585T3 (en)
MY (1) MY121575A (en)
PL (1) PL184098B1 (en)
RU (1) RU2127454C1 (en)
SG (1) SG52253A1 (en)
TR (1) TR199600132A2 (en)
TW (1) TW297970B (en)

Families Citing this family (39)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3484757B2 (en) * 1994-05-13 2004-01-06 ソニー株式会社 Noise reduction method and noise section detection method for voice signal
JP3591068B2 (en) * 1995-06-30 2004-11-17 ソニー株式会社 Noise reduction method for audio signal
JP4104659B2 (en) * 1996-05-31 2008-06-18 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ Device for suppressing disturbing components of input signals
EP0992978A4 (en) * 1998-03-30 2002-01-16 Mitsubishi Electric Corp Noise reduction device and a noise reduction method
JP3454206B2 (en) 1999-11-10 2003-10-06 三菱電機株式会社 Noise suppression device and noise suppression method
US7139711B2 (en) * 2000-11-22 2006-11-21 Defense Group Inc. Noise filtering utilizing non-Gaussian signal statistics
US6985859B2 (en) * 2001-03-28 2006-01-10 Matsushita Electric Industrial Co., Ltd. Robust word-spotting system using an intelligibility criterion for reliable keyword detection under adverse and unknown noisy environments
JP3457293B2 (en) * 2001-06-06 2003-10-14 三菱電機株式会社 Noise suppression device and noise suppression method
JP3427381B2 (en) * 2001-06-20 2003-07-14 富士通株式会社 Noise cancellation method and apparatus
US20040148166A1 (en) * 2001-06-22 2004-07-29 Huimin Zheng Noise-stripping device
US6985709B2 (en) * 2001-06-22 2006-01-10 Intel Corporation Noise dependent filter
US6950036B2 (en) * 2001-09-20 2005-09-27 Honeywell International Inc. Station identification for a local area augmentation system on a visual display
AU2003209821B2 (en) * 2002-03-13 2006-11-16 Hear Ip Pty Ltd A method and system for controlling potentially harmful signals in a signal arranged to convey speech
AUPS102902A0 (en) * 2002-03-13 2002-04-11 Hearworks Pty Ltd A method and system for reducing potentially harmful noise in a signal arranged to convey speech
RU2206960C1 (en) * 2002-06-24 2003-06-20 Общество с ограниченной ответственностью "Центр речевых технологий" Method and device for data signal noise suppression
US7016651B1 (en) * 2002-12-17 2006-03-21 Marvell International Ltd. Apparatus and method for measuring signal quality of a wireless communications link
US6909759B2 (en) * 2002-12-19 2005-06-21 Texas Instruments Incorporated Wireless receiver using noise levels for postscaling an equalized signal having temporal diversity
US6920193B2 (en) * 2002-12-19 2005-07-19 Texas Instruments Incorporated Wireless receiver using noise levels for combining signals having spatial diversity
US7065166B2 (en) 2002-12-19 2006-06-20 Texas Instruments Incorporated Wireless receiver and method for determining a representation of noise level of a signal
GB2398913B (en) * 2003-02-27 2005-08-17 Motorola Inc Noise estimation in speech recognition
CN100417043C (en) * 2003-08-05 2008-09-03 华邦电子股份有限公司 Automatic gain controller and its control method
JP4519169B2 (en) * 2005-02-02 2010-08-04 富士通株式会社 Signal processing method and signal processing apparatus
JP4836720B2 (en) * 2006-09-07 2011-12-14 株式会社東芝 Noise suppressor
GB2450886B (en) * 2007-07-10 2009-12-16 Motorola Inc Voice activity detector and a method of operation
CA2715432C (en) * 2008-03-05 2016-08-16 Voiceage Corporation System and method for enhancing a decoded tonal sound signal
JP5012995B2 (en) * 2008-03-24 2012-08-29 株式会社Jvcケンウッド Audio signal processing apparatus and audio signal processing method
KR101475864B1 (en) 2008-11-13 2014-12-23 삼성전자 주식회사 Apparatus and method for eliminating noise
KR101615766B1 (en) * 2008-12-19 2016-05-12 엘지전자 주식회사 Impulsive noise detector, method of detecting impulsive noise and impulsive noise remover system
FR2944640A1 (en) * 2009-04-17 2010-10-22 France Telecom METHOD AND DEVICE FOR OBJECTIVE EVALUATION OF THE VOICE QUALITY OF A SPEECH SIGNAL TAKING INTO ACCOUNT THE CLASSIFICATION OF THE BACKGROUND NOISE CONTAINED IN THE SIGNAL.
EP2673777B1 (en) * 2011-02-10 2018-12-26 Dolby Laboratories Licensing Corporation Combined suppression of noise and out-of-location signals
US8712076B2 (en) 2012-02-08 2014-04-29 Dolby Laboratories Licensing Corporation Post-processing including median filtering of noise suppression gains
US9173025B2 (en) 2012-02-08 2015-10-27 Dolby Laboratories Licensing Corporation Combined suppression of noise, echo, and out-of-location signals
US9231740B2 (en) 2013-07-12 2016-01-05 Intel Corporation Transmitter noise in system budget
US10504538B2 (en) 2017-06-01 2019-12-10 Sorenson Ip Holdings, Llc Noise reduction by application of two thresholds in each frequency band in audio signals
CN107786709A (en) * 2017-11-09 2018-03-09 广东欧珀移动通信有限公司 Call noise-reduction method, device, terminal device and computer-readable recording medium
CN111199174A (en) * 2018-11-19 2020-05-26 北京京东尚科信息技术有限公司 Information processing method, device, system and computer readable storage medium
CN111477237B (en) * 2019-01-04 2022-01-07 北京京东尚科信息技术有限公司 Audio noise reduction method and device and electronic equipment
CN111429930B (en) * 2020-03-16 2023-02-28 云知声智能科技股份有限公司 Noise reduction model processing method and system based on adaptive sampling rate
CN113035222B (en) * 2021-02-26 2023-10-27 北京安声浩朗科技有限公司 Voice noise reduction method and device, filter determination method and voice interaction equipment

Family Cites Families (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS60140399A (en) * 1983-12-28 1985-07-25 松下電器産業株式会社 Noise remover
US4628529A (en) * 1985-07-01 1986-12-09 Motorola, Inc. Noise suppression system
US4630305A (en) * 1985-07-01 1986-12-16 Motorola, Inc. Automatic gain selector for a noise suppression system
US4630304A (en) * 1985-07-01 1986-12-16 Motorola, Inc. Automatic background noise estimator for a noise suppression system
IL84948A0 (en) * 1987-12-25 1988-06-30 D S P Group Israel Ltd Noise reduction system
US5007094A (en) * 1989-04-07 1991-04-09 Gte Products Corporation Multipulse excited pole-zero filtering approach for noise reduction
US5212764A (en) * 1989-04-19 1993-05-18 Ricoh Company, Ltd. Noise eliminating apparatus and speech recognition apparatus using the same
US5097510A (en) * 1989-11-07 1992-03-17 Gs Systems, Inc. Artificial intelligence pattern-recognition-based noise reduction system for speech processing
CA2032765C (en) * 1989-12-21 1995-12-12 Hidetaka Yoshikawa Variable rate encoding and communicating apparatus
AU633673B2 (en) * 1990-01-18 1993-02-04 Matsushita Electric Industrial Co., Ltd. Signal processing device
JP2797616B2 (en) * 1990-03-16 1998-09-17 松下電器産業株式会社 Noise suppression device
CA2040025A1 (en) * 1990-04-09 1991-10-10 Hideki Satoh Speech detection apparatus with influence of input level and noise reduced
DE69124005T2 (en) * 1990-05-28 1997-07-31 Matsushita Electric Ind Co Ltd Speech signal processing device
DE4137404C2 (en) * 1991-11-14 1997-07-10 Philips Broadcast Television S Method of reducing noise
FI92535C (en) * 1992-02-14 1994-11-25 Nokia Mobile Phones Ltd Noise reduction system for speech signals
JPH05344010A (en) * 1992-06-08 1993-12-24 Mitsubishi Electric Corp Noise reduction device for radio communication equipment
JPH06140949A (en) * 1992-10-27 1994-05-20 Mitsubishi Electric Corp Noise reduction device
US5479560A (en) * 1992-10-30 1995-12-26 Technology Research Association Of Medical And Welfare Apparatus Formant detecting device and speech processing apparatus
EP0707763B1 (en) * 1993-07-07 2001-08-29 Picturetel Corporation Reduction of background noise for speech enhancement
US5617472A (en) * 1993-12-28 1997-04-01 Nec Corporation Noise suppression of acoustic signal in telephone set
JP3484757B2 (en) * 1994-05-13 2004-01-06 ソニー株式会社 Noise reduction method and noise section detection method for voice signal
US5544250A (en) * 1994-07-18 1996-08-06 Motorola Noise suppression system and method therefor

Also Published As

Publication number Publication date
JPH08221093A (en) 1996-08-30
CA2169424A1 (en) 1996-08-18
EP0727769A3 (en) 1998-04-29
KR100414841B1 (en) 2004-03-10
KR960032294A (en) 1996-09-17
SG52253A1 (en) 1998-09-28
TW297970B (en) 1997-02-11
EP0727769B1 (en) 2001-11-21
EP0727769A2 (en) 1996-08-21
ATE209389T1 (en) 2001-12-15
JP3484801B2 (en) 2004-01-06
CN1140869A (en) 1997-01-22
PL312845A1 (en) 1996-08-19
DE69617069T2 (en) 2002-07-11
DE69617069D1 (en) 2002-01-03
RU2127454C1 (en) 1999-03-10
AU4444496A (en) 1996-08-29
BR9600761A (en) 1997-12-23
TR199600132A2 (en) 1996-10-21
ES2163585T3 (en) 2002-02-01
AU696187B2 (en) 1998-09-03
MY121575A (en) 2006-02-28
PL184098B1 (en) 2002-08-30
US6032114A (en) 2000-02-29

Similar Documents

Publication Publication Date Title
CA2169424C (en) Method and apparatus for noise reduction by filtering based on a maximum signal-to-noise ratio and an estimated noise level
US5752226A (en) Method and apparatus for reducing noise in speech signal
US5771486A (en) Method for reducing noise in speech signal and method for detecting noise domain
US5550924A (en) Reduction of background noise for speech enhancement
JP2995737B2 (en) Improved noise suppression system
US6487257B1 (en) Signal noise reduction by time-domain spectral subtraction using fixed filters
JP5203933B2 (en) System and method for reducing audio noise
US7155385B2 (en) Automatic gain control for adjusting gain during non-speech portions
EP1806739A1 (en) Noise suppressor
JP3961290B2 (en) Noise suppressor
JP2004341339A (en) Noise restriction device
US6507623B1 (en) Signal noise reduction by time-domain spectral subtraction
US20030065509A1 (en) Method for improving noise reduction in speech transmission in communication systems

Legal Events

Date Code Title Description
EEER Examination request
MKEX Expiry

Effective date: 20160215