Method for reducing noise in an input speech signal

This invention relates to a method for reducing the noise in speech signals and a method for detecting the noise domain. More particularly, it relates to a method for reducing the noise in the speech signals in which noise suppression is achieved by adaptively controlling a maximum likelihood filter for calculating speech components based upon the speech presence probability and the SN ratio calculated on the basis of input speech signals.

In a portable telephone or speech recognition, it is thought to be necessary to suppress environmental noise or background noise contained in the collected speech signals and to enhance the speech components.

As techniques for enhancing the speech or reducing the noise, those employing a conditional probability function for adjusting attenuation factor are shown in R.J. McAulay and M.L. Malpass, Speech Enhancement Using a Soft-Decision Noise Suppression Filter, IEEE Trans. Acoust, Speech, Signal Processing, Vol.28, pp.137-145, April 1980, and J. Yang, Frequency Domain Noise Suppression Approach in Mobile Telephone System, IEEE ICASSP, vol.II, pp. 363-366, April 1993.

With these noise suppression techniques, it may occur frequently that unnatural speech tone or distorted speech be produced due to the operation based on an inappropriate fixed signal-to-noise (S/N) ratio or to an inappropriate suppression factor. In actual application, it is not desirable for the user to adjust the S/N ratio, which is among the parameters of the noise suppression system for achieving an optimum performance. In addition, it is difficult with the conventional speech signal enhancement techniques to remove the noise sufficiently without by-producing the distortion of the speech signals susceptible to considerable fluctuations in the short-term S/N ratio.

It would be desirable to provide a method for reducing the noise in speech signals whereby the suppression factor is adjusted to a value optimized with respect to the S/N ratio of the actual input responsive to the input speech signals and sufficient noise removal may be achieved without producing distortion as secondary effect or without the necessity of preadjustment by the user.

According to one aspect of the invention, there is provided a method for reducing the noise in an input speech signal in which noise suppression is done by adaptively controlling a maximum likelihood filter adapted for calculating speech components based on the probability of speech occurrence and the S/N ratio calculated based on the input speech signal,    characterised by the characteristics of the maximum likelihood filter being smoothing filtered along the frequency axis and along the time axis.

Preferably, a median value of characteristics of the maximum likelihood filter in neighbouring left and right frequency ranges is used for smoothing filtering along the frequency axis.

Preferably the spectrum of an input signal less an estimated noise spectrum is employed in calculating the probability of speech occurring of speech occurrence.

Preferably the value of the above difference or a pre-set value, whichever is larger, is employed for calculating the probability of speech occurrence.

Preferably the value of the above difference or a pre-set value, whichever is larger, is found for the current frame and a previous frame, the value for the previous frame is multiplied with a pre-set decay coefficient, and the value for the current frame or the value for the previous frame multiplied by a pre-set decay coefficient, whichever is larger, is employed for calculating the probability of speech occurrence.

Preferably, noise domain is detected for finding the probability of speech occurrence by comparing the frame-based RMS values to a threshold value Th₁, a value th for finding the threshold value Th₁ is found responsive to the RMS value for the current frame or the value th of the previous frame multiplied a coefficient α, whichever is smaller, and the coefficient α is changed over depending on the RMS value for the current frame.

Preferably the value th for finding the threshold value Th₁ is found by employing a smaller one of the RMS value of the current frame and the value th of a previous frame multiplied by a coefficient α, whichever is smaller, or the minimum value of the RMS values over plural frames, whichever is larger.

Preferably, the noise domain detection is done by discriminating the relative energy of the current frame using a threshold value Th₂ calculated using the maximum SN ratio of the input speech signal.

According to another aspect, there is provided a method for detecting a noise domain by dividing an input speech signal on the frame basis, finding an RMS value on the frame basis and comparing the RMS values to a threshold value Th₁ for detecting the noise domain, wherein the improvement comprises    calculating a value th for finding the threshold Th₁ using the RMS value for the current frame and a value th of the previous frame multiplied by a coefficient α, whichever is smaller, and changing over the coefficient α depending on an RMS value of the current frame.

Preferably, the method comprises calculating a value th for finding the threshold Th₁ using a smaller one of the RMS value for the current frame and a value th of the previous frame multiplied by a coefficient α, or the smallest RMS value over plural frames, whichever is larger.

Preferably, the noise domain is detected based upon the results of discrimination of the relative energy of the current frame using the threshold value Th₂ calculated using the maximum SN ratio of the input speech signal and the results of comparison of the RMS value to the threshold value Th₁.

The invention will be further described by way of non-limitative example, with reference to the accompanying drawings, in which:-

Fig.1 is a block circuit diagram for illustrating a circuit arrangement for carrying out the noise reducing method for speech signals according to an embodiment of the present invention.

Fig.2 is a block circuit arrangement showing an illustrative example of a noise estimating circuit employed in the embodiment shown in Fig.1.

Fig.3 is a graph showing illustrative examples of an energy E[k] and a decay energy E_decay[k] in the embodiment shown in Fig.1.

Fig.4 is a graph showing illustrative examples of the short-term RMS value RMS[k], minimum noise RMS values MinNoise[k] and the maximum signal RMS values MaxSignal[k] in the embodiment shown in Fig.1.

Fig.5 is a graph showing illustrative examples of the relative energy in dB dB_rel[k], maximum SNR value MaxSNR[k] and dBthres_rel[k] as one of threshold values for noise discrimination.

Fig.6 is a graph for illustrating NR level[k] as a function defined with respect to the maximum SNR value MaxSNR[k] in the embodiment shown in Fig.1.

In Fig.1, a schematic arrangement of the noise reducing device for carrying out the noise reducing method for speech signals according to the preferred embodiment of the present invention is shown in a block circuit diagram.

Referring to Fig.1, an input signal y[t] containing a speech component and a noise component is supplied to an input terminal 11. The input signal y[t], which is a digital signal having the sampling frequency of FS, is fed to a framing/windowing circuit 12 where it is divided into frames each having a length equal to FL samples so that the input signal is subsequently processed on the frame basis. The framing interval, which is the amount of frame movement along the time axis, is FI samples, such that the (k+1)th sample is started after FL samples as from the K'th frame. Prior to processing by a fast Fourier transform (FFT) circuit 13, the next downstream side circuit, the framing/ windowing circuit 12 preforms windowing of the frame-based signals by a windowing function W_input. Meanwhile, after inverse FFT or IFFT at the final stage of signal processing of the frame-based signals, an output signal is processed by windowing by a windowing function W_output. Examples of the windowing functions W_input and W_output are given by the following equations (1) and (2):<maths id="math0001" num="(1)"><math display="block"><mrow><msub><mrow><mtext mathvariant="italic">W</mtext></mrow><mrow><mtext mathvariant="italic">input</mtext></mrow></msub><mtext>[</mtext><mtext mathvariant="italic">j</mtext><mtext>] = (</mtext><mfrac><mrow><mtext>1</mtext></mrow><mrow><mtext>2</mtext></mrow></mfrac><mtext>-</mtext><mfrac><mrow><mtext>1</mtext></mrow><mrow><mtext>2</mtext></mrow></mfrac><mtext>·cos(</mtext><mfrac><mrow><mtext>2·π·</mtext><mtext mathvariant="italic">j</mtext></mrow><mrow><mtext mathvariant="italic">FL</mtext></mrow></mfrac><mtext>))</mtext><msup><mrow><mtext>​</mtext></mrow><mrow><mfrac><mrow><mtext>1</mtext></mrow><mrow><mtext>4</mtext></mrow></mfrac></mrow></msup><mspace linebreak="newline" /><mtext>0 ≤ j ≤ FL</mtext></mrow></math><img file="EP1065656A2_D0001.tif" /></maths><maths id="math0002" num="(2)"><math display="block"><mrow><msub><mrow><mtext mathvariant="italic">W</mtext></mrow><mrow><mtext mathvariant="italic">output</mtext></mrow></msub><mtext>[</mtext><mtext mathvariant="italic">j</mtext><mtext>] = (</mtext><mfrac><mrow><mtext>1</mtext></mrow><mrow><mtext>2</mtext></mrow></mfrac><mtext>-</mtext><mfrac><mrow><mtext>1</mtext></mrow><mrow><mtext>2</mtext></mrow></mfrac><mtext>·cos(</mtext><mfrac><mrow><mtext>2·π·</mtext><mtext mathvariant="italic">j</mtext></mrow><mrow><mtext mathvariant="italic">FL</mtext></mrow></mfrac><mtext>))</mtext><msup><mrow><mtext>​</mtext></mrow><mrow><mfrac><mrow><mtext>3</mtext></mrow><mrow><mtext>4</mtext></mrow></mfrac></mrow></msup><mspace linebreak="newline" /><mtext>0 ≤ j ≤ FL</mtext></mrow></math><img file="EP1065656A2_D0002.tif" /></maths>

If the sampling frequency FS is 8000 Hz = 8 kHz, and the framing interval FI is 80 and 160 samples, the framing interval is 10 msec and 20 msec, respectively.

The FFT circuit 13 performs FFT at 256 points to produce frequency spectral amplitude values which are divided by a frequency dividing circuit 14 into e.g., 18 bands. The following Table 1 shows examples of the frequency ranges of respective bands. <tables id="tabl0001" num="0001"><table frame="all"><title>TABLE 1</title><tgroup cols="2" colsep="1" rowsep="0"><colspec colnum="1" colname="col1" colwidth="78.75mm" /><colspec colnum="2" colname="col2" colwidth="78.75mm" /><thead valign="top"><row rowsep="1"><entry namest="col1" nameend="col1" align="center">Band Number</entry><entry namest="col2" nameend="col2" align="center">Frequency Ranges</entry></row></thead><tbody valign="top"><row><entry namest="col1" nameend="col1" align="center">0</entry><entry namest="col2" nameend="col2" align="center">0 - 125 Hz</entry></row><row><entry namest="col1" nameend="col1" align="center">1</entry><entry namest="col2" nameend="col2" align="center">125 - 250 Hz</entry></row><row><entry namest="col1" nameend="col1" align="center">2</entry><entry namest="col2" nameend="col2" align="center">250 - 375 Hz</entry></row><row><entry namest="col1" nameend="col1" align="center">3</entry><entry namest="col2" nameend="col2" align="center">375 - 563 Hz</entry></row><row><entry namest="col1" nameend="col1" align="center">4</entry><entry namest="col2" nameend="col2" align="center">563 - 750 Hz</entry></row><row><entry namest="col1" nameend="col1" align="center">5</entry><entry namest="col2" nameend="col2" align="center">750 - 938 Hz</entry></row><row><entry namest="col1" nameend="col1" align="center">6</entry><entry namest="col2" nameend="col2" align="center">938 - 1125 Hz</entry></row><row><entry namest="col1" nameend="col1" align="center">7</entry><entry namest="col2" nameend="col2" align="center">1125 - 1313 Hz</entry></row><row><entry namest="col1" nameend="col1" align="center">8</entry><entry namest="col2" nameend="col2" align="center">1313 - 1563 Hz</entry></row><row><entry namest="col1" nameend="col1" align="center">9</entry><entry namest="col2" nameend="col2" align="center">1563 - 1813 Hz</entry></row><row><entry namest="col1" nameend="col1" align="center">10</entry><entry namest="col2" nameend="col2" align="center">1813 - 2063 Hz</entry></row><row><entry namest="col1" nameend="col1" align="center">11</entry><entry namest="col2" nameend="col2" align="center">2063 - 2313 Hz</entry></row><row><entry namest="col1" nameend="col1" align="center">12</entry><entry namest="col2" nameend="col2" align="center">2313 - 2563 Hz</entry></row><row><entry namest="col1" nameend="col1" align="center">13</entry><entry namest="col2" nameend="col2" align="center">2563 - 2813 Hz</entry></row><row><entry namest="col1" nameend="col1" align="center">14</entry><entry namest="col2" nameend="col2" align="center">2813 - 3063 Hz</entry></row><row><entry namest="col1" nameend="col1" align="center">15</entry><entry namest="col2" nameend="col2" align="center">3063 - 3375 Hz</entry></row><row><entry namest="col1" nameend="col1" align="center">16</entry><entry namest="col2" nameend="col2" align="center">3375 - 3688 Hz</entry></row><row rowsep="1"><entry namest="col1" nameend="col1" align="center">17</entry><entry namest="col2" nameend="col2" align="center">3688 - 4000 Hz</entry></row></tbody></tgroup></table></tables>

These frequency bands are set on the basis of the fact that the perceptive resolution of the human auditory system is lowered towards the higher frequency side. As the amplitudes of the respective ranges, the maximum FFT amplitudes in the respective frequency ranges are employed.

A noise estimation circuit 15 distinguishes the noise in the input signal y[t] from the speech and detects a frame which is estimated to be the noise. The operation of estimating the noise domain or detecting the noise frame is performed by combining three kinds of detection operations. An illustrative example of noise domain estimation is hereinafter explained by referring to Fig.2.

In this figure, the input signal y[t] entering the input terminal 11 is fed to a root-mean-square value (RMS) calculating circuit 15A where short-term RMS values are calculated on the frame basis. An output of the RMS calculating circuit 15A is supplied to a relative energy calculating circuit 15B, a minimum RMS calculating circuit 15C, a maximum signal calculating circuit 15D and a noise spectrum estimating circuit 15E. The noise spectrum estimating circuit 15E is fed with outputs of the relative energy calculating circuit 15B, minimum RMS calculating circuit 15C and the maximum signal calculating circuit 15D,while being fed with an output of the frequency dividing circuit 14.

The RMS calculating circuit 15A calculates RMS values of the frame-based signals. The RMS value RMS[k] of the k'th frame is calculated by the following equation:<maths id="math0003" num=""><img file="EP1065656A2_D0003.tif" /></maths>

The relative energy calculating circuit 15B calculates the relative energy dB_rel[k] of the k'th frame pertinent to the decay energy from a previous frame. The relative energy dB_rel[k] in dB is calculated by the following equation (4):<maths id="math0004" num="(4)"><math display="block"><mrow><msub><mrow><mtext mathvariant="italic">dB</mtext></mrow><mrow><mtext mathvariant="italic">rel</mtext></mrow></msub><mtext>[</mtext><mtext mathvariant="italic">k</mtext><msub><mrow><mtext>] = 10log</mtext></mrow><mrow><mtext>10</mtext></mrow></msub><mtext>(</mtext><mfrac><mrow><msub><mrow><mtext mathvariant="italic">E</mtext></mrow><mrow><mtext mathvariant="italic">decay</mtext></mrow></msub><mtext>[</mtext><mtext mathvariant="italic">k</mtext><mtext>]</mtext></mrow><mrow><mtext mathvariant="italic">E</mtext><mtext>[</mtext><mtext mathvariant="italic">k</mtext><mtext>]</mtext></mrow></mfrac><mtext>)</mtext></mrow></math><img file="EP1065656A2_D0004.tif" /></maths>

In the above equation (4), the energy value E[k] and the decay energy value E_decay[k] may be found respectively by the equations (5) and (6):<maths id="math0005" num=""><img file="EP1065656A2_D0005.tif" /></maths><maths id="math0006" num="(6)"><math display="block"><mrow><msub><mrow><mtext mathvariant="italic">E</mtext></mrow><mrow><mtext mathvariant="italic">decay</mtext></mrow></msub><mtext>[</mtext><mtext mathvariant="italic">k</mtext><mtext>] = max(</mtext><mtext mathvariant="italic">E</mtext><mtext>[</mtext><mtext mathvariant="italic">k</mtext><mtext>], </mtext><mtext mathvariant="italic">e</mtext><msup><mrow><mtext>​</mtext></mrow><mrow><mfrac><mrow><mtext>-</mtext><mtext mathvariant="italic">FL</mtext></mrow><mrow><mtext>0.65·</mtext><mtext mathvariant="italic">FS</mtext></mrow></mfrac></mrow></msup><mtext></mtext><msub><mrow><mtext mathvariant="italic">E</mtext></mrow><mrow><mtext mathvariant="italic">decay</mtext></mrow></msub><mtext>[</mtext><mtext mathvariant="italic">k</mtext><mtext>-1])</mtext></mrow></math><img file="EP1065656A2_D0006.tif" /></maths>

Sine the equation (5) may be represented by FL·(RMS[k])², an output RMS[k] of the RMS calculating circuit 15A may be employed. However, the value of the equation (5), obtained in the course of calculation of the equation (3) in the RMS calculating circuit 15A, may be directly transmitted to the relative energy calculating circuit 15B. In the equation (6), the decay time is set to 0.65 sec only by way of an example.

Fig.3 shows illustrative examples of the energy E[k] and the decay energy E_decay[k].

The minimum RMS calculating circuit 15C finds the minimum RMS value suitable for evaluating the background noise level. The frame-based minimum short-term RMS values on the frame-basis and the minimum long-term RMS values, that is the minimum RMS values over plural frames, are found. The long-term values are used when the short-term values cannot track or follow significant changes in the noise level. The minimum short-term RMS noise value MinNoise_short is calculated by the following equation (7):<maths id="math0007" num="(7)"><math display="block"><mrow><msub><mrow><mtext mathvariant="italic">MinNoise</mtext></mrow><mrow><mtext mathvariant="italic">short</mtext></mrow></msub><mtext>[</mtext><mtext mathvariant="italic">k</mtext><mtext>] = min(</mtext><mtext mathvariant="italic">RMS</mtext><mtext>[</mtext><mtext mathvariant="italic">k</mtext><mtext>], max(α(</mtext><mtext mathvariant="italic">k</mtext><mtext>) </mtext><mtext mathvariant="italic">e</mtext><msup><mrow><mtext>​</mtext></mrow><mrow><mfrac><mrow><mtext mathvariant="italic">FL</mtext></mrow><mrow><mtext>0.8·</mtext><mtext mathvariant="italic">FS</mtext></mrow></mfrac></mrow></msup><mtext></mtext><msub><mrow><mtext mathvariant="italic">MinNoise</mtext></mrow><mrow><mtext mathvariant="italic">short</mtext></mrow></msub><mtext>[</mtext><mtext mathvariant="italic">k</mtext><mtext>-1]; </mtext><mtext mathvariant="italic">MinN</mtext></mrow></math><img file="EP1065656A2_D0007.tif" /></maths><dl id="dl0001" compact="compact"><dt>α(k) = 1</dt><dd>RMS[k] &lt; MAX_NOISE_RMS, and RMS[k] &lt; 3 MinNoise_short[k-1]</dd><dt>0</dt><dd>otherwise</dd></dl>

The minimum short-term RMS noise value MinNoise_short is set so as to be increased for the background noise, that is the surrounding noise free of speech. While the rate of rise for the high noise level is exponential, a fixed rise rate is employed for the low noise level for producing a higher rise rate.

The minimum long-term RMS noise value MinNoise_long is calculated for every 0.6 second. MinNoise_long is the minimum over the previous 1.8 second of frame RMS values which have dB_rel &gt; 19 dB. If in the previous 1.8 second, no RMS values have dB_rel &gt; 19 dB, then MinNoise_long is not used because the previous 1 second of signal may not contain any frames with only background noise. At each 0.6 second interval, if MinNoise_long &gt; MinNoise_short, then MinNoise_short at that instance is set to MinNoise_long.

The maximum signal calculating circuit 15D calculates the maximum RMS value or the maximum value of SNR (S/N ratio). The maximum RMS value is used for calculating the optimum or maximum SNR value. For the maximum RMS value, both the short-term and long-term values are calculated. The short-term maximum RMS value MaxSignal_short is found from the following equation (8):<maths id="math0008" num="(8)"><math display="block"><mrow><msub><mrow><mtext mathvariant="italic">MaxSignal</mtext></mrow><mrow><mtext mathvariant="italic">short</mtext></mrow></msub><mtext>[</mtext><mtext mathvariant="italic">k</mtext><mtext>] = </mtext><mtext mathvariant="italic">max</mtext><mtext>(</mtext><mtext mathvariant="italic">RMS</mtext><mtext>[</mtext><mtext mathvariant="italic">k</mtext><mtext>], </mtext><mtext mathvariant="italic">e</mtext><msup><mrow><mtext>​</mtext></mrow><mrow><mtext>-</mtext><mfrac><mrow><mtext mathvariant="italic">FL</mtext></mrow><mrow><mtext>3.2·</mtext><mtext mathvariant="italic">FS</mtext></mrow></mfrac></mrow></msup><mtext>·</mtext><msub><mrow><mtext mathvariant="italic">MaxSignal</mtext></mrow><mrow><mtext mathvariant="italic">short</mtext></mrow></msub><mtext>[</mtext><mtext mathvariant="italic">k</mtext><mtext>-1])</mtext></mrow></math><img file="EP1065656A2_D0008.tif" /></maths>

The maximum long-term RMS noise value MaxSignal_long is calculated at an interval of e.g., 0.4 second. This value MaxSignal_long is the maximum value of the frame RMS value during the term of 0.8 second temporally forward of the current time point. If, during each of the 0.4 second domains, MaxSignal_long is smaller than MaxSignal_short, MaxSignal_short is set to a value of (0.7·MaxSignal_short + 0.3·MaxSignal_long).

Fig.4 shows illustrative values of the short-term RMS value RMS[k], minimum noise RMS value MinNoise[k] and the maximum signal RMS value MaxSignal[k]. In Fig.4, the minimum noise RMS value MinNoise[k] denotes the short-term value of MinNoise_short which takes the long-term value MinNoise_long into account. Also, the maximum signal RMS value MaxSignal[k] denotes the short-term value of MaxSignal_short which takes the long-term value MaxSignal_long into account.

The maximum signal SNR value may be estimated by employing the short-term maximum signal RMS value MaxSignal_short and the short-term minimum noise RMS value MinNoise_short. The noise suppression characteristics and threshold value for noise domain discrimination are modified on the basis of this estimation for reducing the possibility of distorting the noise-free clean speech signal. The maximum SNR value MaxSNR is calculated by the equation:<maths id="math0009" num="(9)"><math display="block"><mrow><mtext mathvariant="italic">MaxSNR</mtext><mtext>[</mtext><mtext mathvariant="italic">k</mtext><msub><mrow><mtext>] = 20.0·log</mtext></mrow><mrow><mtext>10</mtext></mrow></msub><mtext> (</mtext><mfrac><mrow><mtext>max(1000.0, </mtext><msub><mrow><mtext mathvariant="italic">MaxSignal</mtext></mrow><mrow><mtext mathvariant="italic">short</mtext></mrow></msub><mtext>[</mtext><mtext mathvariant="italic">k</mtext><mtext>])</mtext></mrow><mrow><mtext>max(0.5, </mtext><msub><mrow><mtext mathvariant="italic">MinNoise</mtext></mrow><mrow><mtext mathvariant="italic">short</mtext></mrow></msub><mtext>[</mtext><mtext mathvariant="italic">k</mtext><mtext>])</mtext></mrow></mfrac><mtext>-1.0)</mtext></mrow></math><img file="EP1065656A2_D0009.tif" /></maths>

From the value MaxSNR, the normalized parameter NR_level in a range of from 0 to 1 indicating the relative noise level is calculated. The following NT_level function is employed.<maths id="math0010" num=""><img file="EP1065656A2_D0010.tif" /></maths>

The operation of the noise spectrum estimation circuit 15E is explained. The values calculated by the relative energy calculating circuit 15B, minimum RMS calculating circuit 15C and by the maximum signal calculating circuit 15D are used for distinguishing the speech from the background noise. If the following conditions are met, the signal in the k'th frame is classified as being the background noise.<maths id="math0011" num="(11)"><math display="block"><mrow><msub><mrow><mtext>((RMS[k] &lt; NoiseRMS</mtext></mrow><mrow><mtext>thres</mtext></mrow></msub><mtext>[k])</mtext><mspace linebreak="newline" /><msub><mrow><mtext> or (dB</mtext></mrow><mrow><mtext>rel</mtext></mrow></msub><msub><mrow><mtext>[k] &gt; dBthres</mtext></mrow><mrow><mtext>rel</mtext></mrow></msub><mtext>[k])) and (RMS[k] &lt; RMS[k-1] + 200)</mtext></mrow></math><img file="EP1065656A2_D0011.tif" /></maths> where NoiseRMS_rel[k] = min(1.05 + 0.45·NR_level[k]) MinNoise[k], MinNoise[k] + Max_Δ_NOISE_RMS) dBthres_rel[k] = .ax(MaxSNR[k] - 4.0,0.9·MaxSNR[k])

Fig.5 shows illustrative values of the relative energy dB_rel[k], maximum SNR value MaxSNR[k] and the value of dBthres_rel[k], as one of the threshold values of noise discrimination, in the above equation (11).

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

If the k'th frame is classified as being the background noise or the noise, the time averaged estimated value of the noise spectrum Y[w, k] is updated by the signal spectrum Y[w, k] of the current frame, as shown in the following equation (12):<maths id="math0012" num="(12)"><math display="block"><mrow><mtext>N[w,k] = α·max(N[w,k-1], Y[w,k])</mtext><mspace linebreak="newline" /><mtext> + (1-α)·min(N[w,k-1], Y[w,k])</mtext></mrow></math><img file="EP1065656A2_D0012.tif" /></maths><maths id="math0013" num=""><math display="block"><mrow><mtext>α = </mtext><mtext mathvariant="italic">e</mtext><msup><mrow><mtext>​</mtext></mrow><mrow><mtext>-</mtext><mfrac><mrow><mtext mathvariant="italic">FL</mtext></mrow><mrow><mtext>0.5·</mtext><mtext mathvariant="italic">FS</mtext></mrow></mfrac></mrow></msup></mrow></math><img file="EP1065656A2_D0013.tif" /></maths> where <u>w</u> denotes the band number for the frequency 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].

An output of the noise estimation circuit 15 shown in Fig.2 is transmitted to a speech estimation circuit 16, a Pr(Sp) calculating circuit 17, a Pr(Sp|Y) calculating circuit 18 and to a maximum likelihood filter 19.

In carrying out arithmetic-logical operations in the noise spectrum estimation circuit 15E of the noise estimation circuit 15, the arithmetic-logical operations may be carried out using at least one of output data of the relative energy calculating circuit 15B, minimum RMS calculating circuit 15C and the maximum signal calculating circuit 15D. Although the data produced by the estimation circuit 15E is lowered in accuracy, a smaller circuit scale of the noise estimation circuit 15 suffices. Of course, high-accuracy output data of the estimation circuit 15E may be produced by employing all of the output data of the three calculating circuits 15B, 15C and 15D. However, the arithmetic-logical operations by the estimation circuit 15E may be carried out using outputs of two of the calculating circuits 15B, 15C and 15D.

The speech estimation circuit 16 calculates the SN ratio on the band basis. The speech estimation circuit 16 is fed with the spectral amplitude data Y[w, k] from the frequency band splitting circuit 14 and the estimated noise spectral amplitude data from the noise estimation circuit 15. The estimated speech spectral data S[w, k] is derived based upon these data. A rough estimated values of the noise-free clean speech spectrum may be employed for calculating the probability Pr(Sp|Y) as later explained. This value is calculated by taking the difference of spectral values in accordance with the following equation (13).<maths id="math0014" num="(13)"><math display="block"><mrow><mtext mathvariant="italic">S'</mtext><mtext>[</mtext><mtext mathvariant="italic">w,k</mtext><mtext>] = </mtext><msqrt><mtext>max(0,</mtext><mtext mathvariant="italic">Y</mtext><mtext>[</mtext><mtext mathvariant="italic">w,k</mtext><msup><mrow><mtext>]</mtext></mrow><mrow><mtext>2</mtext></mrow></msup><mtext>-ρ·</mtext><mtext mathvariant="italic">N</mtext><mtext>[</mtext><mtext mathvariant="italic">w,k</mtext><msup><mrow><mtext>]</mtext></mrow><mrow><mtext>2</mtext></mrow></msup><mtext>)</mtext></msqrt></mrow></math><img file="EP1065656A2_D0014.tif" /></maths>

Then, using the rough estimated value S'[w, k] of the speech spectrum as calculated by the above equation (13), an estimated value S[w, k] of the speech spectrum, time-averaged on the band basis, is calculated in accordance with the following equation (14):<maths id="math0015" num="(14)"><math display="block"><mrow><mtext>S[w,k] = max(S'[w,k], S'[w,k-1]·decay_rate)</mtext><mspace linebreak="newline" /><msub><mrow><mtext mathvariant="italic">decay</mtext></mrow><mrow><mtext mathvariant="italic">rate</mtext></mrow></msub><mtext> = </mtext><mtext mathvariant="italic">e</mtext><msup><mrow><mtext>​</mtext></mrow><mrow><mfrac><mrow><mtext>-</mtext><mtext mathvariant="italic">FL</mtext></mrow><mrow><mtext>(4-0.5·</mtext><msub><mrow><mtext mathvariant="italic">nr</mtext></mrow><mrow><mtext mathvariant="italic">level</mtext></mrow></msub><mtext>)·</mtext><mtext mathvariant="italic">FS</mtext></mrow></mfrac></mrow></msup></mrow></math><img file="EP1065656A2_D0015.tif" /></maths>

In the equation (14), the decay_rate shown therein is employed.

The band-based SN ratio is calculated in accordance with the following equation (15):<maths id="math0016" num="(15)"><math display="block"><mrow><mtext mathvariant="italic">SNR</mtext><mtext>[</mtext><mtext mathvariant="italic">w,k</mtext><msub><mrow><mtext>] = 20·log</mtext></mrow><mrow><mtext>10</mtext></mrow></msub><mtext> (</mtext><mfrac><mrow><mtext>0.2·</mtext><mtext mathvariant="italic">S</mtext><mtext>[</mtext><mtext mathvariant="italic">w</mtext><mtext>-1,</mtext><mtext mathvariant="italic">k</mtext><mtext>]+0.6·</mtext><mtext mathvariant="italic">S</mtext><mtext>[</mtext><mtext mathvariant="italic">w</mtext><mtext>,</mtext><mtext mathvariant="italic">k</mtext><mtext>]+0.2</mtext><mtext mathvariant="italic">S</mtext><mtext>[</mtext><mtext mathvariant="italic">w</mtext><mtext>+1,</mtext><mtext mathvariant="italic">k</mtext><mtext>]</mtext></mrow><mrow><mtext>0.2·</mtext><mtext mathvariant="italic">N</mtext><mtext>[</mtext><mtext mathvariant="italic">w</mtext><mtext>+1,</mtext><mtext mathvariant="italic">k</mtext><mtext>]+0.6·</mtext><mtext mathvariant="italic">N</mtext><mtext>[</mtext><mtext mathvariant="italic">w</mtext><mtext>,</mtext><mtext mathvariant="italic">k</mtext><mtext>]+0.2</mtext><mtext mathvariant="italic">N</mtext><mtext>[</mtext><mtext mathvariant="italic">w</mtext><mtext>-1,</mtext><mtext mathvariant="italic">k</mtext><mtext>]</mtext></mrow></mfrac><mtext>)</mtext></mrow></math><img file="EP1065656A2_D0016.tif" /></maths> where the estimated value of the noise spectrum N[ ] and the estimated value of the speech spectrum may be found from the equations (12) and (14), respectively.

The operation of the Pr(Sp) calculating circuit 17 is explained. The probability Pr(Sp) is the probability of the speech signals occurring in an assumed input signal. This probability was hitherto fixed perpetually to 0.5. For a signal having a high SN ratio, the probability Pr(Sp) can be increased for prohibiting sound quality deterioration. Such probability Pr(Sp) may be calculated in accordance with the following equation (16):<maths id="math0017" num="(16)"><math display="block"><mrow><mtext>Pr(Sp) = 0.5 + 0.45·(1.0 - NR_level)</mtext></mrow></math><img file="EP1065656A2_D0017.tif" /></maths> using the NR_level function calculated by the maximum signal calculating circuit 15D.

The operation of the Pr(Sp|Y) calculating circuit 18 is now explained. The value Pr(Sp|Y) is the probability of the speech signal occurring in the input signal y[t], and is calculated using Pr(Sp) and SNR[w, k]. The value Pr(Sp|Y) is used for reducing the speech-free domain to a narrower value. For calculations, the method disclosed in R.J. McAulay and M.L. Malpass, Speech Enhancement Using a Soft-Decision Noise Suppression Filter, IEEE Trans. Acoust, Speech, and Signal Processing, Vo. ASSP-28, No.2, April 1980, which is now explained by referring to equations (17) to (20), was employed.<maths id="math0018" num="(17)"><math display="block"><mrow><mtext mathvariant="italic">Pr</mtext><mtext>(</mtext><mtext mathvariant="italic">H1</mtext><mtext>|</mtext><mtext mathvariant="italic">Y</mtext><mtext>) [</mtext><mtext mathvariant="italic">w,k</mtext><mtext>] = </mtext><mfrac><mrow><mtext mathvariant="italic">Pr</mtext><mtext>(</mtext><mtext mathvariant="italic">H1</mtext><mtext>)·</mtext><mtext mathvariant="italic">p</mtext><mtext>(</mtext><mtext mathvariant="italic">Y</mtext><mtext>|</mtext><mtext mathvariant="italic">H1</mtext><mtext>)</mtext></mrow><mrow><mtext mathvariant="italic">Pr</mtext><mtext>(</mtext><mtext mathvariant="italic">H1</mtext><mtext>)·</mtext><mtext mathvariant="italic">p</mtext><mtext>(</mtext><mtext mathvariant="italic">Y</mtext><mtext>|</mtext><mtext mathvariant="italic">H1</mtext><mtext>)+</mtext><mtext mathvariant="italic">Pr</mtext><mtext>(</mtext><mtext mathvariant="italic">H0</mtext><mtext>)·</mtext><mtext mathvariant="italic">p</mtext><mtext>(</mtext><mtext mathvariant="italic">Y</mtext><mtext>|</mtext><mtext mathvariant="italic">H0</mtext><mtext>)</mtext></mrow></mfrac><mspace linebreak="newline" /><mtext>(Bayes Rule)</mtext></mrow></math><img file="EP1065656A2_D0018.tif" /></maths><maths id="math0019" num="(18)"><math display="block"><mrow><mtext mathvariant="italic">p</mtext><mtext>(</mtext><mtext mathvariant="italic">Y</mtext><mtext>|</mtext><mtext mathvariant="italic">H0</mtext><mtext>) = </mtext><mfrac><mrow><mtext>2·</mtext><mtext mathvariant="italic">Y</mtext></mrow><mrow><mtext>σ</mtext></mrow></mfrac><mtext>·</mtext><mtext mathvariant="italic">e</mtext><msup><mrow><mtext>​</mtext></mrow><mrow><mtext>-</mtext><mfrac><mrow><msup><mrow><mtext mathvariant="italic">Y</mtext></mrow><mrow><mtext>2</mtext></mrow></msup></mrow><mrow><mtext>σ</mtext></mrow></mfrac></mrow></msup><mspace linebreak="newline" /><mtext>(Rayleigh pdf)</mtext></mrow></math><img file="EP1065656A2_D0019.tif" /></maths><maths id="math0020" num="(19)"><math display="block"><mrow><mtext mathvariant="italic">p</mtext><mtext>(</mtext><mtext mathvariant="italic">Y</mtext><mtext>|</mtext><mtext mathvariant="italic">H1</mtext><mtext>) = </mtext><mfrac><mrow><mtext>2·</mtext><mtext mathvariant="italic">Y</mtext></mrow><mrow><mtext>σ</mtext></mrow></mfrac><mtext>·</mtext><mtext mathvariant="italic">e</mtext><msup><mrow><mtext>​</mtext></mrow><mrow><mtext>-</mtext><mfrac><mrow><msup><mrow><mtext mathvariant="italic">Y</mtext></mrow><mrow><mtext>2</mtext></mrow></msup><mtext>+</mtext><msup><mrow><mtext mathvariant="italic">S</mtext></mrow><mrow><mtext>2</mtext></mrow></msup></mrow><mrow><mtext>σ</mtext></mrow></mfrac></mrow></msup><mtext>·</mtext><msub><mrow><mtext mathvariant="italic">I</mtext></mrow><mrow><mtext>0</mtext></mrow></msub><mtext>(</mtext><mfrac><mrow><mtext>2·</mtext><mtext mathvariant="italic">S</mtext><mtext>·</mtext><mtext mathvariant="italic">Y</mtext></mrow><mrow><mtext>σ</mtext></mrow></mfrac><mtext>)</mtext><mspace linebreak="newline" /><mtext>(Rician pdf)</mtext></mrow></math><img file="EP1065656A2_D0020.tif" /></maths><maths id="math0021" num=""><img file="EP1065656A2_D0021.tif" /></maths>

In the above equations (17) to (20), HO denotes a non-speech event, that is the event that the input signal y(t) is the noise signal n(t), while H1 denotes a speech event, that is the event that the input signal y(t) is a sum of the speech signal s(t) and the noise signal n(t) and s(t) is not equal to 0. In addition, <u>w</u>, <u>k</u>, Y, S and σ denote the band number, frame number, input signal [w, k], estimated value of the speech signal S[w, k] and a square value of the estimated noise signal N[w, k]², respectively.

Pr(H1∼Y)[w, k] is calculated from the equation (17), while p(Y|HO) and p(Y|H1) in the equation (17) may be found from the equation (19). The Bessel function I₀(|X|) is calculated from the equation (20).

The Bessel function may be approximated by the following function (21):<maths id="math0022" num="(21)"><math display="block"><mrow><msub><mrow><mtext>I</mtext></mrow><mrow><mtext>0</mtext></mrow></msub><mtext>(|x|) =</mtext><mspace linebreak="newline" /><mfrac><mrow><mtext>1</mtext></mrow><mrow><msqrt><mtext>2π|</mtext><mtext mathvariant="italic">x</mtext><mtext>|</mtext></msqrt></mrow></mfrac><mtext>·</mtext><msup><mrow><mtext mathvariant="italic">e</mtext></mrow><mrow><mtext>|</mtext><mtext mathvariant="italic">x</mtext><mtext>|</mtext></mrow></msup><mtext>+0.07, </mtext><mtext mathvariant="italic">if</mtext><mtext> |</mtext><mtext mathvariant="italic">x</mtext><mtext>|≥0.5</mtext><mspace linebreak="newline" /><mtext> 1   otherwise</mtext></mrow></math><img file="EP1065656A2_D0022.tif" /></maths>

Heretofore, a fixed value of the SN ratio, such as SNR = 5, was employed for deriving Pr(H1|Y) without employing the estimated speech signal value S[w, k]. Consequently, p(Y|H1) was simplified as shown by the following equation (22):<maths id="math0023" num="(22)"><math display="block"><mrow><mtext mathvariant="italic">p</mtext><mtext>(</mtext><mtext mathvariant="italic">Y</mtext><mtext>|</mtext><mtext mathvariant="italic">H1</mtext><mtext>) = </mtext><mfrac><mrow><mtext>2</mtext></mrow><mrow><mtext>σ</mtext></mrow></mfrac><mtext>·</mtext><mtext mathvariant="italic">e</mtext><msup><mrow><mtext>​</mtext></mrow><mrow><mtext>-</mtext><mfrac><mrow><msup><mrow><mtext mathvariant="italic">Y</mtext></mrow><mrow><mtext>2</mtext></mrow></msup></mrow><mrow><mtext>σ</mtext></mrow></mfrac><mtext>-</mtext><msup><mrow><mtext mathvariant="italic">SNR</mtext></mrow><mrow><mtext>2</mtext></mrow></msup></mrow></msup><mtext>·</mtext><msub><mrow><mtext mathvariant="italic">I</mtext></mrow><mrow><mtext>0</mtext></mrow></msub><mtext> (2·</mtext><mtext mathvariant="italic">SNR</mtext><mtext>·</mtext><mfrac><mrow><mtext mathvariant="italic">Y</mtext></mrow><mrow><msqrt><mtext>σ</mtext></msqrt></mrow></mfrac><mtext>)</mtext></mrow></math><img file="EP1065656A2_D0023.tif" /></maths>

A signal having an instantaneous SN ratio lower than the value SNR of the SN ratio employed in the calculation of p(Y|H1) is suppressed significantly. If it is assumed that the value SNR of the SN ratio is set to an excessively high value, the speech corrupted by a noise of a lower level is excessively lowered in its low-level speech portion, so that the produced speech becomes unnatural. Conversely, if the value SNR of the SN ratio is set to an excessively low value, the speech corrupted by the larger level noise is low in suppression and sounds noisy even at its low-level portion. Thus the value of p(Y|H1) conforming to a wide range of the background/speech level is obtained by using the variable value of the SN ratio SNR_new[w, k] as in the present embodiment instead of by using the fixed value of the SN ratio. The value of SNR_new[w, k] may be found from the following equation (23):<maths id="math0024" num=""><img file="EP1065656A2_D0024.tif" /></maths> in which the value of MIN_SNR is found from the equation (24):<maths id="math0025" num="(24)"><math display="block"><mrow><mtext mathvariant="italic">MIN_SNR</mtext><mtext>(</mtext><mtext mathvariant="italic">x</mtext><mtext>) =</mtext><mspace linebreak="newline" /><mtext>3,</mtext><mtext mathvariant="italic">x</mtext><mtext>&lt;10</mtext><mspace linebreak="newline" /><mtext>3-</mtext><mfrac><mrow><mtext mathvariant="italic">x</mtext><mtext>-10</mtext></mrow><mrow><mtext>35</mtext></mrow></mfrac><mtext>·1.5, 10≤</mtext><mtext mathvariant="italic">x</mtext><mtext>≤45</mtext><mspace linebreak="newline" /><mtext>1.5, </mtext><mtext mathvariant="italic">otherwise</mtext></mrow></math><img file="EP1065656A2_D0025.tif" /></maths>

The value SNR_new[w, k] is an instantaneous SNR in the k'th frame in which limitation is placed on the minimum value. The value of SNR_new[w, k] may be decreased to 1.5 for a signal having the high SN ratio on the whole. In such case, suppression is not done on segments having low instantaneous SN ratio. The value SNR_new[w, k] cannot be lowered to below 3 for a signal having a low instantaneous SN ratio as a whole. Consequently, sufficient suppression my be assured for segments having a low instantaneous S/N ratio.

The operation of the maximum likelihood filter 19 is explained. The maximum likelihood filter 19 is one of pre-filters provided for freeing the respective bands of the input signal of noise signals. In the most likelihood filter 19, the spectral amplitude dataY[w, k] from the frequency band splitting filter 14 is converted into a signal H[w, k] using the noise spectral amplitude data N[w, k] from the noise estimation circuit 15. The signal H[w, k] is calculated in accordance with he following equation (25):<maths id="math0026" num="(25)"><math display="block"><mrow><mtext mathvariant="italic">H</mtext><mtext>[</mtext><mtext mathvariant="italic">w,k</mtext><mtext>] =</mtext><mspace linebreak="newline" /><mtext mathvariant="italic">a</mtext><mtext>+(1-</mtext><mtext mathvariant="italic">a</mtext><mtext>)·</mtext><mfrac><mrow><mtext>(</mtext><msup><mrow><mtext mathvariant="italic">Y</mtext></mrow><mrow><mtext>2</mtext></mrow></msup><mtext>-</mtext><msup><mrow><mtext mathvariant="italic">N</mtext></mrow><mrow><mtext>2</mtext></mrow></msup><mtext>)</mtext><msup><mrow><mtext>​</mtext></mrow><mrow><mfrac><mrow><mtext>1</mtext></mrow><mrow><mtext>2</mtext></mrow></mfrac></mrow></msup></mrow><mrow><mtext mathvariant="italic">Y</mtext></mrow></mfrac><mtext>, </mtext><mtext mathvariant="italic">Y</mtext><mtext>&gt;0 </mtext><mtext mathvariant="italic">and Y</mtext><mtext>≥</mtext><mtext mathvariant="italic">N</mtext><mspace linebreak="newline" /><mtext mathvariant="italic">a, otherwise</mtext></mrow></math><img file="EP1065656A2_D0026.tif" /></maths> where α = 0.7 - 0.4·NR_level[k].

Although the value α in the above equation (25) is conventionally set to 1/2, the degree of noise suppression may be varied depending on the maximum SNR because an approximate value of the SNR is known.

The operation of a soft decision suppression circuit 21 is now explained. The soft decision suppression circuit 20 is one of pre-filters for enhancing the speech portion of the signal. Conversion is done by the method shown in the following equation (26) using the signal H[w, k] and the value Pr(H1|Y) from the Pr(Sp|Y) calculating circuit 18:<maths id="math0027" num="(26)"><math display="block"><mrow><mtext>H[w,k] ← Pr(H1|Y)[w,k]·H[w,k]+(1-Pr(H1|Y[w,k]·MIN_GAIN</mtext></mrow></math><img file="EP1065656A2_D0027.tif" /></maths>

In the above equation (26), MIN_GAIN is a parameter indicating the minimum gain, and may be set to, for example, 0.1, that is -15 dB.

The operation of a filter processing circuit 21 is now explained. The signal H[w, k] from the soft decision suppression circuit 20 is filtered along both the frequency axis and the time axis. The filtering along the frequency axis has the effect of shortening the effective impulse response length of the signal H[w, k]. This eliminates any circular convolution aliasing effects associated with filtering by multiplication in the frequency domain. The filtering along the time axis has the effect of limiting the rate of change of the filter in suppressing noise bursts.

The filtering along the frequency axis is now explained. Median filtering is done on the signals H[w, k] of each of 18 bands resulting from frequency band division. The method is explained by the following equations (27) and (28): <ul id="ul0001" list-style="none" compact="compact"><li>Step 1 :<maths id="math0028" num="(27)"><math display="block"><mrow><mtext>H1[w, k] = max(median(H[w-1, k], H[w, k], H[w+1, k], H[w, k]</mtext></mrow></math><img file="EP1065656A2_D0028.tif" /></maths>    where H1[w, k] = H[w, k] if (w-1) or (w+1) is absent</li><li>Step 2 :<maths id="math0029" num="(27)"><math display="block"><mrow><mtext>H2[w, k] = min(median(H[w-1, k], H[w, k], H[w+1, k], H[w, k]</mtext></mrow></math><img file="EP1065656A2_D0029.tif" /></maths>    where H2[w, k] = H1[w, k] if (w-1) or (w+1) is absent.</li></ul>

In the step 1, H1[w, k] is H[w, k] without single band nulls. In the step 2, H2[w, k] is H1[w,k] without sole band spikes. The signal resulting from filtering along the frequency axis is H2[w, k].

Next, the filtering along the time axis is explained. The filtering along time axis considers three states of the input speech signal, namely the speech, the background noise and the transient which is the rising portion of the speech. The speech signal ia smoothed along the time axis as shown by the following equation (29).<maths id="math0030" num="(29)"><math display="block"><mrow><msub><mrow><mtext>H</mtext></mrow><mrow><mtext>speech</mtext></mrow></msub><mtext>[w, k] = 0.7·H2[w, k] + 0.3·H2[w, k-1]</mtext></mrow></math><img file="EP1065656A2_D0030.tif" /></maths>

The background noise signal is smoothed along the time axis as shown by the following equation (30):<maths id="math0031" num="(30)"><math display="block"><mrow><msub><mrow><mtext>H</mtext></mrow><mrow><mtext>noise</mtext></mrow></msub><mtext>[w, k] = 0.7·Min_H + 0.3·Max_H</mtext></mrow></math><img file="EP1065656A2_D0031.tif" /></maths> where Min_H and Max_H are: <ul id="ul0002" list-style="none" compact="compact"><li>Min_H = min(H2[w, k], H2[w, k-1]</li><li>Max_H = max(H2[w, k], H2[w, k-1]</li></ul> For transient signals, no smoothing on time axis is not performed. Ultimately, calculations are carried out for producing the smoothed output signal H_{t_smooth}[w, k] by the following equation (31):<maths id="math0032" num="(31)"><math display="block"><mrow><msub><mrow><mtext>H</mtext></mrow><mrow><mtext>t_smooth</mtext></mrow></msub><mtext>[w, k]</mtext><mspace linebreak="newline" /><msub><mrow><mtext> = (1-α</mtext></mrow><mrow><mtext>tr</mtext></mrow></msub><msub><mrow><mtext>)(α</mtext></mrow><mrow><mtext>sp</mtext></mrow></msub><msub><mrow><mtext>·H</mtext></mrow><mrow><mtext>speech</mtext></mrow></msub><msub><mrow><mtext>[w, k] + (1-α</mtext></mrow><mrow><mtext>sp</mtext></mrow></msub><msub><mrow><mtext>)·H</mtext></mrow><mrow><mtext>noise</mtext></mrow></msub><msub><mrow><mtext>[w, k] + α</mtext></mrow><mrow><mtext>tr</mtext></mrow></msub><mtext>·H2[w, k]</mtext></mrow></math><img file="EP1065656A2_D0032.tif" /></maths>    α_sp and α_tr in the equation (31) are respectively found from the equations (32) and (33):<maths id="math0033" num="(32)"><math display="block"><mrow><msub><mrow><mtext>α</mtext></mrow><mrow><mtext mathvariant="italic">sp</mtext></mrow></msub><mtext> =</mtext><mspace linebreak="newline" /><mtext>1.0, </mtext><msub><mrow><mtext mathvariant="italic">SNR</mtext></mrow><mrow><mtext mathvariant="italic">inst</mtext></mrow></msub><mtext>&gt;4.0</mtext><mspace linebreak="newline" /><mtext>(</mtext><msub><mrow><mtext mathvariant="italic">SNR</mtext></mrow><mrow><mtext mathvariant="italic">inst</mtext></mrow></msub><mtext>-1)·</mtext><mfrac><mrow><mtext>1</mtext></mrow><mrow><mtext>3</mtext></mrow></mfrac><mtext>, 1.0&lt;</mtext><msub><mrow><mtext mathvariant="italic">SNR</mtext></mrow><mrow><mtext mathvariant="italic">inst</mtext></mrow></msub><mtext>&lt;4.0</mtext><mspace linebreak="newline" /><mtext>0, </mtext><mtext mathvariant="italic">otherwise</mtext></mrow></math><img file="EP1065656A2_D0033.tif" /></maths> where<maths id="math0034" num="(33)"><math display="block"><mrow><msub><mrow><mtext mathvariant="italic">SNR</mtext></mrow><mrow><mtext mathvariant="italic">inst</mtext></mrow></msub><mtext> = </mtext><mfrac><mrow><mtext mathvariant="italic">RMS</mtext><mtext>[</mtext><mtext mathvariant="italic">k</mtext><mtext>]</mtext></mrow><mrow><mtext mathvariant="italic">MinNoise</mtext><mtext>[</mtext><mtext mathvariant="italic">k</mtext><mtext>]</mtext></mrow></mfrac><mspace linebreak="newline" /><msub><mrow><mtext> α</mtext></mrow><mrow><mtext mathvariant="italic">tr</mtext></mrow></msub><mtext> =</mtext><mspace linebreak="newline" /><msub><mrow><mtext>1.0, δ</mtext></mrow><mrow><mtext mathvariant="italic">rms</mtext></mrow></msub><mtext>&gt;3.5</mtext><mspace linebreak="newline" /><msub><mrow><mtext>(δ</mtext></mrow><mrow><mtext mathvariant="italic">rms</mtext></mrow></msub><mtext>-2)·</mtext><mfrac><mrow><mtext>2</mtext></mrow><mrow><mtext>3</mtext></mrow></mfrac><msub><mrow><mtext>, 2.0&lt;δ</mtext></mrow><mrow><mtext mathvariant="italic">rms</mtext></mrow></msub><mtext>&lt;3.5</mtext><mspace linebreak="newline" /><mtext>0, </mtext><mtext mathvariant="italic">otherwise</mtext></mrow></math><img file="EP1065656A2_D0034.tif" /></maths> where<maths id="math0035" num=""><img file="EP1065656A2_D0035.tif" /></maths>

The operation in a band conversion circuit 22 is explained. The 18 band signals H_{t_smooth}[w, k] from the filtering circuit 21 is interpolated to e.g., 128 band signals H₁₂₈[w, k]. The interpolation is done in two stages, that is, the interpolation from 18 to 64 bands is done by zero-order hold and the interpolation from 64 to 128 bands is done by a low-pass filter interpolation.

The operation in a spectrum correction circuit 23 is explained. The real part and the imaginary part of the FFT coefficients of the input signal obtained at the FFT circuit 13 are multiplied with the above signal H₁₂₈[w, k] to carry out spectrum correction. The result is that the spectral amplitude is corrected, while the spectrum is not modified in phase.

An IFFT circuit 24 executes inverse FFT on the signal obtained at the spectrum correction circuit 23.

An overlap-and-add circuit 25 overlap and adds the frame boundary portions of the frame-based IFFT output signals. A noise-reduced output signal is obtained at an output terminal 26 by the procedure described above.

The output signal thus obtained is transmitted to various encoders of a portable telephone set or to a signal processing circuit of a speech recognition device. Alternatively, decoder output signals of a portable telephone set may be processed with noise reduction according to the present invention.

The present invention is not limited to the above embodiment. For example, the above-described filtering by the filtering circuit 21 may be employed in the conventional noise suppression technique employing the maximum likelihood filter. The noise domain detection method by the filter processing circuit 15 may be employed in a variety of devices other than the noise suppression device.