Method and device for estimating the tonality of a sound signal

<u>Field</u><u>of</u><u>the</u><u>Invention</u>

0001The present invention relates to sound activity detection, background noise estimation and sound signal classification where sound is understood as a useful signal. The present invention also relates to corresponding sound activity detector, background noise estimator and sound signal classifier.

0002In particular but not exclusively: <ul id="ul0001" list-style="dash"><li>The sound activity detection is used to select frames to be encoded using techniques optimized for inactive frames.</li><li>The sound signal classifier is used to discriminate among different speech signal classes and music to allow for more efficient encoding of sound signals, i.e. optimized encoding of unvoiced speech signals, optimized encoding of stable voiced speech signals, and generic encoding of other sound signals.</li><li>An algorithm is provided and uses several relevant parameters and features to allow for a better choice of coding mode and more robust estimation of the background noise.</li><li>Tonality estimation is used to improve the performance of sound activity detection in the presence of music signals, and to better discriminate between unvoiced sounds and music. For example, the tonality estimation may be used in a super-wideband codec to decide the codec model to encode the signal above 7 kHz.</li></ul>

<u>Background of</u><u>the</u><u>Invention</u>

0003Demand for efficient digital narrowband and wideband speech coding techniques with a good trade-off between the subjective quality and bit rate is increasing in various application areas such as teleconferencing, multimedia, and wireless communications. Until recently, telephone bandwidth constrained into a range of 200-3400 Hz has mainly been used in speech coding applications (signal sampled at 8 kHz). However, wideband speech applications provide increased intelligibility and naturalness in communication compared to the conventional telephone bandwidth. In wideband services the input signal is sampled at 16 kHz and the encoded bandwidth is in the range 50-7000 Hz. This bandwidth has been found sufficient for delivering a good quality giving an impression of nearly face-to-face communication. Further quality improvement is achieved with so-called super-wideband, in which the signal is sampled at 32 kHz and the encoded bandwidth is in the range 50-15000 Hz. For speech signals this provides a face-to-face quality since almost all energy in human speech is below 14000 Hz. This bandwidth also gives significant quality improvement with general audio signals including music (wideband is equivalent to AM radio and super-wideband is equivalent to FM radio). Higher bandwidth has been used for general audio signals with the full-band 20-20000 Hz (CD quality sampled at 44.1 kHz or 48 kHz).

0004A sound encoder converts a sound signal (speech or audio) into a digital bit stream which is transmitted over a communication channel or stored in a storage medium. The sound signal is digitized, that is, sampled and quantized with usually 16-bits per sample. The sound encoder has the role of representing these digital samples with a smaller number of bits while maintaining a good subjective quality. The sound decoder operates on the transmitted or stored bit stream and converts it back to a sound signal.

0005<i>Code-Excited Linear Prediction</i> (CELP) coding is one of the best prior techniques for achieving a good compromise between the subjective quality and bit rate. This coding technique is a basis of several speech coding standards both in wireless and wireline applications. In CELP coding, the sampled speech signal is processed in successive blocks of <i>L</i> samples usually called <i>frames,</i> where <i>L</i> is a predetermined number corresponding typically to 10-30 ms. A linear prediction (LP) filter is computed and transmitted every frame. The L-sample frame is divided into smaller blocks called <i>subframes.</i> In each subframe, an excitation signal is usually obtained from two components, the past excitation and the innovative, fixed-codebook excitation. The component formed from the past excitation is often referred to as the adaptive codebook or pitch excitation. The parameters characterizing the excitation signal are coded and transmitted to the decoder, where the reconstructed excitation signal is used as the input of the LP filter.

0006The use of source-controlled variable bit rate (VBR) speech coding significantly improves the system capacity. In source-controlled VBR coding, the codec uses a signal classification module and an optimized coding model is used for encoding each speech frame based on the nature of the speech frame (e.g. voiced, unvoiced, transient, background noise). Further, different bit rates can be used for each class. The simplest form of source-controlled VBR coding is to use voice activity detection (VAD) and encode the inactive speech frames (background noise) at a very low bit rate. Discontinuous transmission (DTX) can further be used where no data is transmitted in the case of stable background noise. The decoder uses comfort noise generation (CNG) to generate the background noise characteristics. VAD/DTX/CNG results in significant reduction in the average bit rate, and in packet-switched applications it reduces significantly the number of routed packets. VAD algorithms work well with speech signals but may result in severe problems in case of music signals. Segments of music signals can be classified as unvoiced signals and consequently may be encoded with unvoiced-optimized model which severely affects the music quality. Moreover, some segments of stable music signals may be classified as stable background noise and this may trigger the update of background noise in the VAD algorithm which results in degradation in the performance of the algorithm. Therefore, it would be advantageous to extend the VAD algorithm to better discriminate music signals. In the present disclosure, this algorithm will be referred to as Sound Activity Detection (SAD) algorithm where sound could be speech or music or any useful signal. The present disclosure also describes a method for tonality detection used to improve the performance of the SAD algorithm in case of music signals.

0007Another aspect in speech and audio coding is the concept of embedded coding, also known as layered coding. In embedded coding, the signal is encoded in a first layer to produce a first bit stream, and then the error between the original signal and the encoded signal from the first layer is further encoded to produce a second bit stream. This can be repeated for more layers by encoding the error between the original signal and the coded signal from all preceding layers. The bit streams of all layers are concatenated for transmission. The advantage of layered coding is that parts of the bit stream (corresponding to upper layers) can be dropped in the network (e.g. in case of congestion) while still being able to decode the signal at the receiver depending on the number of received layers. Layered encoding is also useful in multicast applications where the encoder produces the bit stream of all layers and the network decides to send different bit rates to different end points depending on the available bit rate in each link.

0008Embedded or layered coding can be also useful to improve the quality of widely used existing codecs while still maintaining interoperability with these codecs. Adding more layers to the standard codec core layer can improve the quality and even increase the encoded audio signal bandwidth. Examples are the recently standardized ITU-T Recommendation G.729.1 where the core layer is interoperable with widely used G.729 narrowband standard at 8 kbit/s and upper layers produces bit rates up to 32 kbit/s (with wideband signal starting from 16 kbit/s). Current standardization work aims at adding more layers to produce a super-wideband codec (14 kHz bandwidth) and stereo extensions. Another example is ITU-T Recommendation G.718 for encoding wideband signals at 8, 12, 16, 24 and 32 kbit/s. The codec is also being extended to encode super-wideband and stereo signals at higher bit rates.

0009The requirements for embedded codecs usually ask for good quality in case of both speech and audio signals. Since speech can be encoded at relatively low bit rate using a model based approach, the first layer (or first two layers) is (or are) encoded using a speech specific technique and the error signal for the upper layers is encoded using a more generic audio encoding technique. This delivers a good speech quality at low bit rates and good audio quality as the bit rate is increased. In G.718 and G.729.1, the first two layers are based on ACELP (Algebraic Code-Excited Linear Prediction) technique which is suitable for encoding speech signals. In the upper layers, transform-based encoding suitable for audio signals is used to encode the error signal (the difference between the original signal and the output from the first two layers). The well known MDCT (Modified Discrete Cosine Transform) transform is used, where the error signal is transformed in the frequency domain. In the super-wideband layers, the signal above 7 kHz is encoded using a generic coding model or a tonal coding model. The above mentioned tonality detection can also be used to select the proper coding model to be used.

0010An example of a known method and apparatus for determining the tonality of an input audio signal is disclosed in the patent document <patcit id="pcit0001" dnum="US2004181393A1"><text>US 2004/181393 A1</text></patcit>.

<u>Summary of the Invention</u>

0011According to a first aspect of the present invention, there is provided a method for estimating a tonality of a sound signal. The method comprises: calculating a current residual spectrum of the sound signal; detecting peaks in the current residual spectrum; calculating a correlation map between the current residual spectrum and a previous residual spectrum for each detected peak; and calculating a long-term correlation map based on the calculated correlation map, the long-term correlation map being indicative of a tonality in the sound signal.

0012According to a further aspect of the present invention, there is provided a device for estimating a tonality of a sound signal. The device comprises: a calculator a current residual spectrum of the sound signal; a detector for detecting peaks in the current residual spectrum; a calculator for calculating a correlation map between the current residual spectrum and a previous residual spectrum for each detected peak; and a calculator for calculating a long-term correlation map based on the calculated correlation map, the long-term correlation map being indicative of a tonality in the sound signal.

0013The foregoing and other objects, advantages and features of the present invention will become more apparent upon reading of the following non restrictive description of an illustrative embodiment thereof, given by way of example only with reference to the accompanying drawings.

<u>Brief</u><u>Description of</u><u>the</u><u>Drawings</u>

0014In the appended drawings: <ul id="ul0002" list-style="none"><li><figref idref="f0001">Figure 1</figref> is a schematic block diagram of a portion of an example of sound communication system including sound activity detection, background noise estimation update, and sound signal classification;</li><li><figref idref="f0002">Figure 2</figref> is a non-limitative illustration of windowing in spectral analysis;</li><li><figref idref="f0003">Figure 3</figref> is a non-restrictive graphical illustration of the principle of spectral floor calculation and the residual spectrum;</li><li><figref idref="f0004">Figure 4</figref> is a non-limitative illustration of calculation of spectral correlation map in a current frame;</li><li><figref idref="f0005">Figure 5</figref> is an example of functional block diagram of a signal classification algorithm; and</li><li><figref idref="f0006">Figure 6</figref> is an example of decision tree for unvoiced speech discrimination.</li></ul>

<u>Detailed</u><u>description</u>

0015In the non-restrictive, illustrative embodiment of the present invention, sound activity detection (SAD) is performed within a sound communication system to classify short-time frames of signals as sound or background noise/silence. The sound activity detection is based on a frequency dependent signal-to-noise ratio (SNR) and uses an estimated background noise energy per critical band. A decision on the update of the background noise estimator is based on several parameters including parameters discriminating between background noise/silence and music, thereby preventing the update of the background noise estimator on music signals.

0016The SAD corresponds to a first stage of the signal classification. This first stage is used to discriminate inactive frames for optimized encoding of inactive signal. In a second stage, unvoiced speech frames are discriminated for optimized encoding of unvoiced signal. At this second stage, music detection is added in order to prevent classifying music as unvoiced signal. Finally, in a third stage, voiced signals are discriminated through further examination of the frame parameters.

0017The herein disclosed techniques can be deployed with either narrowband (NB) sound signals sampled at 8000 sample/s or wideband (WB) sound signals sampled at 16000 sample/s, or at any other sampling frequency. The encoder used in the non-restrictive, illustrative embodiment of the present invention is based on AMR-WB [<nplcit id="ncit0001" npl-type="s" url="http://www.3gpp.org"><text>AMR Wideband Speech Codec: Transcoding Functions, 3GPP Technical Specification TS 26.190 (http://www.3gpp.org</text></nplcit>)] and VMR-WB [<nplcit id="ncit0002" npl-type="s" url="http://www.3gpp2.org"><text>Source-Controlled Variable-Rate Multimode Wideband Speech Codec (VMR-WB), Service Options 62 and 63 for Spread Spectrum Systems, 3GPP2 Technical Specification C.S0052-A v1.0, April 2005 (http://www.3gpp2.org</text></nplcit>)] codecs which use an internal sampling conversion to convert the signal sampling frequency to 12800 sample/s (operating in a 6.4 kHz bandwidth). Thus the sound activity detection technique in the non-restrictive, illustrative embodiment operates on either narrowband or wideband signals after sampling conversion to 12.8 kHz.

0018<figref idref="f0001">Figure 1</figref> is a block diagram of a sound communication system 100 according to the non-restrictive illustrative embodiment of the invention, including sound activity detection.

0019The sound communication system 100 of <figref idref="f0001">Figure 1</figref> comprises a pre-processor 101. Preprocessing by module 101 can be performed as described in the following example (high-pass filtering, resampling and pre-emphasis).

0020Prior to the frequency conversion, the input sound signal is high-pass filtered. In this non-restrictive, illustrative embodiment, the cut-off frequency of the high-pass filter is 25 Hz for WB and 100 Hz for NB. The high-pass filter serves as a precaution against undesired low frequency components. For example, the following transfer function can be used: <maths id="math0001" num=""><math display="block"><msub><mi>H</mi><mrow><mi>h</mi><mo>⁢</mo><mn>1</mn></mrow></msub><mfenced><mi>z</mi></mfenced><mo>=</mo><mfrac><mrow><msub><mi mathvariant="normal">b</mi><mn mathvariant="normal">0</mn></msub><mo>+</mo><msub><mi mathvariant="normal">b</mi><mn mathvariant="normal">1</mn></msub><mo>⁢</mo><msup><mi>z</mi><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup><mo>+</mo><msub><mi mathvariant="normal">b</mi><mn mathvariant="normal">2</mn></msub><mo>⁢</mo><msup><mi>z</mi><mrow><mo>-</mo><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow><mn mathvariant="normal">1</mn><mo>+</mo><msub><mi mathvariant="normal">a</mi><mn mathvariant="normal">1</mn></msub><mo>⁢</mo><msup><mi>z</mi><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup><mo>+</mo><msub><mi mathvariant="normal">a</mi><mn mathvariant="normal">2</mn></msub><mo>⁢</mo><msup><mi>z</mi><mrow><mo>-</mo><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac></math><img file="EP2162880B1_D0001.tif" /></maths> where, for WB, <i>b</i>₀ = 0.9930820, <i>b</i>₁ = -1.98616407, <i>b</i>₂ = 0.9930820, <i>a</i>₁ = -1.9861162, <i>a</i>₂ = 0.9862119292 and, for NB, <i>b</i>₀ = 0.945976856, <i>b</i>₁ = -1.891953712, <i>b</i>₂ = 0.945976856, <i>a</i>₁ = -1.889033079, <i>a</i>₂ = 0.894874345. Obviously, the high-pass filtering can be alternatively carried out after resampling to 12.8 kHz.

0021In the case of WB, the input sound signal is decimated from 16 kHz to 12.8 kHz. The decimation is performed by an upsampler that upsamples the sound signal by 4. The resulting output is then filtered through a low-pass FIR (Finite Impulse Response) filter with a cut off frequency at 6.4 kHz. Then, the low-pass filtered signal is downsampled by 5 by an appropriate downsampler. The filtering delay is 15 samples at a 16 kHz sampling frequency.

0022In the case of NB, the sound signal is upsampled from 8 kHz to 12.8 kHz. For that purpose, an upsampler performs on the sound signal an upsampling by 8. The resulting output is then filtered through a low-pass FIR filter with a cut off frequency at 6.4 kHz. A downsampler then downsamples the low-pass filtered signal by 5. The filtering delay is 16 samples at 8 kHz sampling frequency.

0023After the sampling conversion, a pre-emphasis is applied to the sound signal prior to the encoding process. In the pre-emphasis, a first order high-pass filter is used to emphasize higher frequencies. This first order high-pass filter forms a pre-emphasizer and uses, for example, the following transfer function: <maths id="math0002" num=""><math display="block"><msub><mi>H</mi><mrow><mi>pre</mi><mo>-</mo><mi>emph</mi></mrow></msub><mfenced><mi>z</mi></mfenced><mo>=</mo><mn>1</mn><mo>-</mo><mn>0.68</mn><mo>⁢</mo><msup><mi>z</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup></math><img file="EP2162880B1_D0002.tif" /></maths>

0024Pre-emphasis is used to improve the codec performance at high frequencies and improve perceptual weighting in the error minimization process used in the encoder.

0025As described hereinabove, the input sound signal is converted to 12.8 kHz sampling frequency and preprocessed, for example as described above. However, the disclosed techniques can be equally applied to signals at other sampling frequencies such as 8 kHz or 16 kHz with different preprocessing or without preprocessing.

0026In the non-restrictive illustrative embodiment of the present invention, the encoder 109 (<figref idref="f0001">Figure 1</figref>) using sound activity detection operates on 20 ms frames containing 256 samples at the 12.8 kHz sampling frequency. Also, the encoder 109 uses a 10 ms look ahead from the future frame to perform its analysis (<figref idref="f0002">Figure 2</figref>). The sound activity detection follows the same framing structure.

0027Referring to <figref idref="f0001">Figure 1</figref>, spectral analysis is performed in spectral analyzer 102. Two analyses are performed in each frame using 20 ms windows with 50% overlap. The windowing principle is illustrated in <figref idref="f0002">Figure 2</figref>. The signal energy is computed for frequency bins and for critical bands [<nplcit id="ncit0003" npl-type="s"><text>J. D. Johnston, "Transform coding of audio signal using perceptual noise criteria," IEEE J. Select. Areas Commun., vol. 6, pp. 314-323, February 1988</text></nplcit>].

0028Sound activity detection (first stage of signal classification) is performed in the sound activity detector 103 using noise energy estimates calculated in the previous frame. The output of the sound activity detector 103 is a binary variable which is further used by the encoder 109 and which determines whether the current frame is encoded as active or inactive.

0029Noise estimator 104 updates a noise estimation downwards (first level of noise estimation and update), i.e. if in a critical band the frame energy is lower than an estimated energy of the background noise, the energy of the noise estimation is updated in that critical band.

0030Noise reduction is optionally applied by an optional noise reducer 105 to the speech signal using for example a spectral subtraction method. An example of such a noise reduction scheme is described in [<nplcit id="ncit0004" npl-type="s"><text>M. Jelínek and R. Salami, "Noise Reduction Method for Wideband Speech Coding," in Proc. Eusipco, Vienna, Austria, September 2004</text></nplcit>].

0031Linear prediction (LP) analysis and open-loop pitch analysis are performed (usually as a part of the speech coding algorithm) by a LP analyzer and pitch tracker 106. In this non-restrictive illustrative embodiment, the parameters resulting from the LP analyzer and pitch tracker 106 are used in the decision to update the noise estimates in the critical bands as performed in module 107. Alternatively, the sound activity detector 103 can also be used to take the noise update decision. According to a further alternative, the functions implemented by the LP analyzer and pitch tracker 106 can be an integral part of the sound encoding algorithm.

0032Prior to updating the noise energy estimates in module 107, music detection is performed to prevent false updating on active music signals. Music detection uses spectral parameters calculated by the spectral analyzer 102.

0033Finally, the noise energy estimates are updated in module 107 (second level of noise estimation and update). This module 107 uses all available parameters calculated previously in modules 102 to 106 to decide about the update of the energies of the noise estimation.

0034In signal classifier 108, the sound signal is further classified as unvoiced, stable voiced or generic. Several parameters are calculated to support this decision. In this signal classifier, the mode of encoding the sound signal of the current frame is chosen to best represent the class of signal being encoded.

0035Sound encoder 109 performs encoding of the sound signal based on the encoding mode selected in the sound signal classifier 108. In other applications, the sound signal classifier 108 can be an automatic speech recognition system.

<i>Spectral analysis</i>

0036The spectral analysis is performed by the spectral analyzer 102 of <figref idref="f0001">Figure 1</figref>.

0037Fourier Transform is used to perform the spectral analysis and spectrum energy estimation. The spectral analysis is done twice per frame using a 256-point Fast Fourier Transform (FFT) with a 50 percent overlap (as illustrated in <figref idref="f0002">Figure 2</figref>). The analysis windows are placed so that all look ahead is exploited. The beginning of the first window is at the beginning of the encoder current frame. The second window is placed 128 samples further. A square root Harming window (which is equivalent to a sine window) has been used to weight the input sound signal for the spectral analysis. This window is particularly well suited for overlap-add methods (thus this particular spectral analysis is used in the noise suppression based on spectral subtraction and overlap-add analysis/synthesis). The square root Harming window is given by: <maths id="math0003" num="(1)"><math display="block"><msub><mi>w</mi><mi mathvariant="italic">FFT</mi></msub><mfenced><mi>n</mi></mfenced><mo>=</mo><msqrt><mn>0.5</mn><mo>-</mo><mn>0.5</mn><mo>⁢</mo><mi>cos</mi><mfenced><mfrac><mrow><mn>2</mn><mo>⁢</mo><mi mathvariant="italic">πn</mi></mrow><msub><mi>L</mi><mi mathvariant="italic">FFT</mi></msub></mfrac></mfenced></msqrt><mo>=</mo><mi>sin</mi><mfenced><mfrac><mi mathvariant="italic">πn</mi><msub><mi>L</mi><mi mathvariant="italic">FFT</mi></msub></mfrac></mfenced><mo>,</mo><mspace width="1em" /><mi>n</mi><mo>=</mo><mn>0</mn><mo>,</mo><mo>…</mo><mo>,</mo><msub><mi>L</mi><mi mathvariant="italic">FFT</mi></msub><mo>-</mo><mn>1</mn></math><img file="EP2162880B1_D0003.tif" /></maths> where <i>L_FFT</i>=256 is the size of the FTT analysis. Here, only half the window is computed and stored since this window is symmetric (from 0 to <i>L_FFTl</i>2).

0038The windowed signals for both spectral analyses (first and second spectral analyses) are obtained using the two following relations: <maths id="math0004" num=""><math display="block"><msubsup><mi>x</mi><mi>w</mi><mfenced><mn>1</mn></mfenced></msubsup><mfenced><mi>n</mi></mfenced><mo>=</mo><msub><mi>w</mi><mi mathvariant="italic">FFT</mi></msub><mfenced><mi>n</mi></mfenced><mo>⁢</mo><mi mathvariant="italic">sʹ</mi><mfenced><mi>n</mi></mfenced><mo>,</mo><mspace width="5em" /><mi>n</mi><mo>=</mo><mn>0</mn><mo>,</mo><mo>…</mo><mo>,</mo><msub><mi>L</mi><mi mathvariant="italic">FFT</mi></msub><mo>-</mo><mn>1</mn></math><img file="EP2162880B1_D0004.tif" /></maths><maths id="math0005" num=""><math display="block"><msubsup><mi>x</mi><mi>w</mi><mfenced><mn>2</mn></mfenced></msubsup><mfenced><mi>n</mi></mfenced><mo>=</mo><msub><mi>w</mi><mi mathvariant="italic">FFT</mi></msub><mfenced><mi>n</mi></mfenced><mo>⁢</mo><mi mathvariant="italic">sʹ</mi><mo>⁢</mo><mfenced separators=""><mi>n</mi><mo>+</mo><msub><mi>L</mi><mi mathvariant="italic">FFT</mi></msub><mo>/</mo><mn>2</mn></mfenced><mo>,</mo><mspace width="2em" /><mi>n</mi><mo>=</mo><mn>0</mn><mo>,</mo><mo>…</mo><mo>,</mo><msub><mi>L</mi><mi mathvariant="italic">FFT</mi></msub><mo>-</mo><mn>1</mn></math><img file="EP2162880B1_D0005.tif" /></maths> where <i>s</i>'(0) is the first sample in the current frame. In the non-restrictive, illustrative embodiment of the present invention, the beginning of the first window is placed at the beginning of the current frame. The second window is placed 128 samples further.

0039FFT is performed on both windowed signals to obtain following two sets of spectral parameters per frame: <maths id="math0006" num=""><math display="block"><msup><mi>X</mi><mfenced><mn>1</mn></mfenced></msup><mfenced><mi>k</mi></mfenced><mo>=</mo><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover></mstyle><msubsup><mi>x</mi><mi>w</mi><mfenced><mn>1</mn></mfenced></msubsup><mfenced><mi>n</mi></mfenced><mo>⁢</mo><msup><mi>e</mi><mrow><mo>-</mo><mi>j</mi><mo>⁢</mo><mn>2</mn><mo>⁢</mo><mi>π</mi><mo>⁢</mo><mfrac><mi mathvariant="italic">kn</mi><mi>N</mi></mfrac></mrow></msup><mo>,</mo><mspace width="2em" /><mi>k</mi><mo>=</mo><mn>0</mn><mo>,</mo><mo>…</mo><mo>,</mo><msub><mi>L</mi><mi mathvariant="italic">FFT</mi></msub><mo>-</mo><mn>1</mn></math><img file="EP2162880B1_D0006.tif" /></maths><maths id="math0007" num=""><math display="block"><msup><mi>X</mi><mfenced><mn>2</mn></mfenced></msup><mfenced><mi>k</mi></mfenced><mo>=</mo><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><mrow><mi>N</mi><mo>-</mo><mn>1</mn></mrow></munderover></mstyle><msubsup><mi>x</mi><mi>w</mi><mfenced><mn>2</mn></mfenced></msubsup><mfenced><mi>n</mi></mfenced><mo>⁢</mo><msup><mi>e</mi><mrow><mo>-</mo><mi>j</mi><mo>⁢</mo><mn>2</mn><mo>⁢</mo><mi>π</mi><mo>⁢</mo><mfrac><mi mathvariant="italic">kn</mi><mi>N</mi></mfrac></mrow></msup><mo>,</mo><mspace width="2em" /><mi>k</mi><mo>=</mo><mn>0</mn><mo>,</mo><mo>…</mo><mo>,</mo><msub><mi>L</mi><mi mathvariant="italic">FFT</mi></msub><mo>-</mo><mn>1</mn></math><img file="EP2162880B1_D0007.tif" /></maths> where <i>N</i>= <i>L_FFT.</i>

0040The FFT provides the real and imaginary parts of the spectrum denoted by <i>X_R</i>(<i>k</i>), <i>k</i>=0 to 128, and <i>X_I</i>(<i>k</i>), <i>k=</i>1 to 127. <i>X_R</i>(0) corresponds to the spectrum at 0 Hz (DC) and <i>X_R</i>(128) corresponds to the spectrum at 6400 Hz. The spectrum at these points is only real valued.

0041After FFT analysis, the resulting spectrum is divided into critical bands using the intervals having the following upper limits [<nplcit id="ncit0005" npl-type="s"><text>M. Jelinek and R. Salami, "Noise Reduction Method for Wideband Speech Coding," in Proc. Eusipco, Vienna, Austria, September 2004</text></nplcit>] (20 bands in the frequency range 0-6400 Hz): Critical bands = {100.0, 200.0, 300.0, 400.0, 510.0, 630.0, 770.0, 920.0, 1080.0, 1270.0, 1480.0, 1720.0, 2000.0, 2320.0, 2700.0, 3150.0, 3700.0, 4400.0, 5300.0, 6350.0} Hz.

0042The 256-point FFT results in a frequency resolution of 50 Hz (6400/128). Thus after ignoring the DC component of the spectrum, the number of frequency bins per critical band is <i>M_CB</i> = {2, 2, 2, 2, 2, 2, 3, 3, 3, 4, 4, 5, 6, 6, 8, 9, 11, 14, 18, 21}, respectively.

0043The average energy in a critical band is computed using the following relation: <maths id="math0008" num="(2)"><math display="block"><msub><mi>E</mi><mi mathvariant="italic">CB</mi></msub><mfenced><mi>i</mi></mfenced><mo>=</mo><mfrac><mn>1</mn><mrow><msup><mfenced separators=""><msub><mi>L</mi><mi mathvariant="italic">FFT</mi></msub><mo>/</mo><mn>2</mn></mfenced><mn>2</mn></msup><mo>⁢</mo><msub><mi>M</mi><mi mathvariant="italic">CB</mi></msub><mfenced><mi>i</mi></mfenced></mrow></mfrac><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>k</mi><mo>=</mo><mn>0</mn></mrow><mrow><msub><mi>M</mi><mi mathvariant="italic">CB</mi></msub><mfenced><mi>i</mi></mfenced><mo>-</mo><mn>1</mn></mrow></munderover></mstyle><mfenced separators=""><msubsup><mi>X</mi><mi>R</mi><mn>2</mn></msubsup><mo>⁢</mo><mfenced separators=""><mi>k</mi><mo>+</mo><msub><mi>j</mi><mi>i</mi></msub></mfenced><mo>+</mo><msubsup><mi>X</mi><mi>I</mi><mn>2</mn></msubsup><mo>⁢</mo><mfenced separators=""><mi>k</mi><mo>+</mo><msub><mi>j</mi><mi>i</mi></msub></mfenced></mfenced><mo>,</mo><mspace width="2em" /><mi>i</mi><mo>=</mo><mn>0</mn><mo>,</mo><mo>…</mo><mo>,</mo><mn>19</mn></math><img file="EP2162880B1_D0008.tif" /></maths> where <i>K_R</i>(<i>k</i>) and <i>X_l</i>(<i>k</i>) are, respectively, the real and imaginary parts of the <i>k</i>^th frequency bin and <i>j_i</i> is the index of the first bin in the <i>i</i>^th critical band given by <i>j_i</i> ={1, 3, 5, 7, 9, 11, 13, 16, 19, 22, 26, 30, 35, 41, 47, 55, 64, 75, 89, 107}.

0044The spectral analyzer 102 also computes the normalized energy per frequency bin, <i>E_BIN</i>(<i>k</i>)<i>,</i> in the range 0-6400 Hz, using the following relation: <maths id="math0009" num="(3)"><math display="block"><msub><mi>E</mi><mi mathvariant="italic">BIN</mi></msub><mfenced><mi>k</mi></mfenced><mo>=</mo><mfrac><mn>4</mn><msubsup><mi>L</mi><mi mathvariant="italic">FFT</mi><mn>2</mn></msubsup></mfrac><mo>⁢</mo><mfenced separators=""><msubsup><mi>X</mi><mi>R</mi><mn>2</mn></msubsup><mfenced><mi>k</mi></mfenced><mo>+</mo><msubsup><mi>X</mi><mi>I</mi><mn>2</mn></msubsup><mfenced><mi>k</mi></mfenced></mfenced><mo>,</mo><mspace width="2em" /><mi>k</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>…</mo><mo>,</mo><mn>127</mn></math><img file="EP2162880B1_D0009.tif" /></maths> Furthermore, the energy spectra per frequency bin in both analyses are combined together to obtain the average log-energy spectrum (in decibels), i.e. <maths id="math0010" num="(4)"><math display="block"><msub><mi>E</mi><mi mathvariant="italic">dB</mi></msub><mfenced><mi>k</mi></mfenced><mo>=</mo><mn>10</mn><mo>⁢</mo><mi>log</mi><mo>⁢</mo><mfenced open="[" close="]" separators=""><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>⁢</mo><mfenced separators=""><msubsup><mi>E</mi><mi mathvariant="italic">Bin</mi><mfenced><mn>1</mn></mfenced></msubsup><mfenced><mi>k</mi></mfenced><mo>+</mo><msubsup><mi>E</mi><mi mathvariant="italic">Bin</mi><mfenced><mn>2</mn></mfenced></msubsup><mfenced><mi>k</mi></mfenced></mfenced></mfenced><mo>,</mo><mspace width="2em" /><mi>k</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>…</mo><mo>,</mo><mn>127</mn><mo>,</mo></math><img file="EP2162880B1_D0010.tif" /></maths> where the superscripts (1) and (2) are used to denote the first and the second spectral analysis, respectively.

0045Finally, the spectral analyzer 102 computes the average total energy for both the first and second spectral analyses in a 20 ms frame by adding the average critical band energies <i>E_CB.</i> That is, the spectrum energy for a certain spectral analysis is computed using the following relation: <maths id="math0011" num="(5)"><math display="block"><msub><mi>E</mi><mi mathvariant="italic">frame</mi></msub><mo>=</mo><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mn>19</mn></munderover></mstyle><msub><mi>E</mi><mi mathvariant="italic">CB</mi></msub><mfenced><mi>i</mi></mfenced></math><img file="EP2162880B1_D0011.tif" /></maths> and the total frame energy is computed as the average of spectrum energies of both the first and second spectral analyses in a frame. That is <maths id="math0012" num="(6)"><math display="block"><msub><mi>E</mi><mi>t</mi></msub><mo>=</mo><mn>10</mn><mo>⁢</mo><mi>log</mi><mrow><mo>(</mo><mn>0.5</mn><mo>⁢</mo><mfenced separators=""><msub><mi>E</mi><mi mathvariant="italic">frame</mi></msub><mfenced><mn>0</mn></mfenced><mo>+</mo><msub><mi>E</mi><mi mathvariant="italic">frame</mi></msub><mfenced><mn>1</mn></mfenced></mfenced></mrow><mo>,</mo><mspace width="1em" /><mi>dB</mi><mn>.</mn></math><img file="EP2162880B1_D0012.tif" /></maths>

0046The output parameters of the spectral analyzer 102, that is the average energy per critical band, the energy per frequency bin and the total energy, are used in the sound activity detector 103 and in the rate selection. The average log-energy spectrum is used in the music detection.

0047In narrowband input signals sampled at 8000 sample/s, after sampling conversion to 12800 sample/s, there is no content at both ends of the spectrum, thus the first lower frequency critical band as well as the last three high frequency bands are not considered in the computation of relevant parameters (only bands from i=1 to 16 are considered). However, equations (3) and (4) are not affected.

<i>Sound activity detection (SAD)</i>

0048The sound activity detection is performed by the SNR-based sound activity detector 103 of <figref idref="f0001">Figure 1</figref>.

0049The spectral analysis described above is performed twice per frame by the analyzer 102. Let <maths id="math0013" num=""><math display="inline"><msubsup><mi>E</mi><mi mathvariant="italic">CB</mi><mfenced><mn>1</mn></mfenced></msubsup><mfenced><mi>i</mi></mfenced></math><img file="EP2162880B1_D0013.tif" /></maths>and <maths id="math0014" num=""><math display="inline"><msubsup><mi>E</mi><mi mathvariant="italic">CB</mi><mfenced><mn>2</mn></mfenced></msubsup><mfenced><mi>i</mi></mfenced></math><img file="EP2162880B1_D0014.tif" /></maths>as computed in Equation (2) denote the energy per critical band information in the first and second spectral analyses, respectively. The average energy per critical band for the whole frame and part of the previous frame is computed using the following relation: <maths id="math0015" num="(7)"><math display="block"><msub><mi>E</mi><mi mathvariant="italic">av</mi></msub><mfenced><mi>i</mi></mfenced><mo>=</mo><mn>0.2</mn><mo>⁢</mo><msubsup><mi>E</mi><mi mathvariant="italic">CB</mi><mfenced><mn>0</mn></mfenced></msubsup><mfenced><mi>i</mi></mfenced><mo>+</mo><mn>0.4</mn><mo>⁢</mo><msubsup><mi>E</mi><mi mathvariant="italic">CB</mi><mfenced><mn>1</mn></mfenced></msubsup><mfenced><mi>i</mi></mfenced><mo>+</mo><mn>0.4</mn><mo>⁢</mo><msubsup><mi>E</mi><mi mathvariant="italic">CB</mi><mfenced><mn>2</mn></mfenced></msubsup><mfenced><mi>i</mi></mfenced></math><img file="EP2162880B1_D0015.tif" /></maths> where <maths id="math0016" num=""><math display="inline"><msubsup><mi>E</mi><mi mathvariant="italic">CB</mi><mfenced><mn>0</mn></mfenced></msubsup><mfenced><mi>i</mi></mfenced></math><img file="EP2162880B1_D0016.tif" /></maths> denotes the energy per critical band information from the second spectral analysis of the previous frame. The signal-to-noise ratio (SNR) per critical band is then computed using the following relation: <maths id="math0017" num="(8)"><math display="block"><msub><mi mathvariant="italic">SNR</mi><mi mathvariant="italic">CB</mi></msub><mfenced><mi>i</mi></mfenced><mo>=</mo><msub><mi>E</mi><mi mathvariant="italic">av</mi></msub><mfenced><mi>i</mi></mfenced><mo>/</mo><msub><mi>N</mi><mi mathvariant="italic">CB</mi></msub><mfenced><mi>i</mi></mfenced><mspace width="2em" /><mi>bounded by</mi><mspace width="1em" /><msub><mi mathvariant="italic">SNR</mi><mi mathvariant="italic">CB</mi></msub><mo>≥</mo><mn>1.</mn></math><img file="EP2162880B1_D0017.tif" /></maths> where <i>N_CB(i)</i> is the estimated noise energy per critical band as will be explained below. The average SNR per frame is then computed as <maths id="math0018" num="(9)"><math display="block"><msub><mi mathvariant="italic">SNR</mi><mi mathvariant="italic">av</mi></msub><mo>=</mo><mn>10</mn><mo>⁢</mo><mi>log</mi><mfenced separators=""><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><msub><mi>b</mi><mi>min</mi></msub></mrow><msub><mi>b</mi><mi>max</mi></msub></munderover></mstyle><msub><mi mathvariant="italic">SNR</mi><mi mathvariant="italic">CB</mi></msub><mfenced><mi>i</mi></mfenced></mfenced><mo>,</mo></math><img file="EP2162880B1_D0018.tif" /></maths> where <i>b_min</i>=0 and <i>b_max</i>=19 in the case of wideband signals, and <i>b_min</i>=1 and <i>b_max</i>=16 in case of narrowband signals.

0050The sound activity is detected by comparing the average SNR per frame to a certain threshold which is a function of the long-term SNR. The long-term SNR is given by the following relation: <maths id="math0019" num="(10)"><math display="block"><msub><mi mathvariant="italic">SNR</mi><mi mathvariant="italic">LT</mi></msub><mo>=</mo><msub><mover><mi>E</mi><mo>‾</mo></mover><mi>f</mi></msub><mo>-</mo><msub><mover><mi>N</mi><mo>‾</mo></mover><mi>f</mi></msub></math><img file="EP2162880B1_D0019.tif" /></maths> where <i><o ostyle="single">E</o>_f</i> and <i><o ostyle="single">N</o>_f</i> are computed using equations (13) and (14), respectively, which will be described later. The initial value of <i><o ostyle="single">E</o>_f</i> is 45 dB.

0051The threshold is a piece-wise linear function of the long-term SNR. Two functions are used, one optimized for clean speech and one optimized for noisy speech.

0052For wideband signals, If <i>SNR_LT</i> &lt; 35 (noisy speech) then the threshold is equal to: <maths id="math0020" num=""><math display="block"><msub><mi mathvariant="italic">th</mi><mi mathvariant="italic">SAD</mi></msub><mo>=</mo><mn>0.41287</mn><mspace width="1em" /><msub><mi mathvariant="italic">SNR</mi><mi mathvariant="italic">LT</mi></msub><mo>+</mo><mn>13.259625</mn></math><img file="EP2162880B1_D0020.tif" /></maths> else (clean speech): <maths id="math0021" num=""><math display="block"><msub><mi mathvariant="italic">th</mi><mi mathvariant="italic">SAD</mi></msub><mo>=</mo><mn>1.0333</mn><mspace width="1em" /><msub><mi mathvariant="italic">SNR</mi><mi mathvariant="italic">LT</mi></msub><mo>-</mo><mn>18</mn></math><img file="EP2162880B1_D0021.tif" /></maths>

0053For narrowband signals, If <i>SNR_LT &lt;</i> 20 (noisy speech) then the threshold is equal to: <maths id="math0022" num=""><math display="block"><msub><mi mathvariant="italic">th</mi><mi mathvariant="italic">SAD</mi></msub><mo>=</mo><mn>0.1071</mn><mspace width="1em" /><msub><mi mathvariant="italic">SNR</mi><mi mathvariant="italic">LT</mi></msub><mo>+</mo><mn>16.5</mn></math><img file="EP2162880B1_D0022.tif" /></maths> else (clean speech): <maths id="math0023" num=""><math display="block"><msub><mi mathvariant="italic">th</mi><mi mathvariant="italic">SAD</mi></msub><mo>=</mo><mn>0.4773</mn><mspace width="1em" /><msub><mi mathvariant="italic">SNR</mi><mi mathvariant="italic">LT</mi></msub><mo>-</mo><mn>6.1364</mn></math><img file="EP2162880B1_D0023.tif" /></maths>

0054Furthermore, a hysteresis in the SAD decision is added to prevent frequent switching at the end of an active sound period. The hysteresis strategy is different for wideband and narrowband signals and comes into effect only if the signal is noisy.

0055For wideband signals, the hysteresis strategy is applied in the case the frame is in a "hangover period" the length of which varies according to the long-term SNR as follows: <maths id="math0024" num=""><math display="block"><msub><mi>l</mi><mi mathvariant="italic">hang</mi></msub><mo>=</mo><mn>0</mn><mspace width="4em" /><mi>if</mi><mspace width="1em" /><msub><mi mathvariant="italic">SNR</mi><mi mathvariant="italic">LT</mi></msub><mo>≥</mo><mn>35</mn></math><img file="EP2162880B1_D0024.tif" /></maths><maths id="math0025" num=""><math display="block"><msub><mi>l</mi><mi mathvariant="italic">hang</mi></msub><mo>=</mo><mn>1</mn><mspace width="4em" /><mi>if</mi><mspace width="1em" /><mn>15</mn><mo>≤</mo><msub><mi mathvariant="italic">SNR</mi><mi mathvariant="italic">LT</mi></msub><mo>&lt;</mo><mn>35.</mn></math><img file="EP2162880B1_D0025.tif" /></maths><maths id="math0026" num=""><math display="block"><msub><mi>l</mi><mi mathvariant="italic">hang</mi></msub><mo>=</mo><mn>2</mn><mspace width="4em" /><mi>if</mi><mspace width="1em" /><msub><mi mathvariant="italic">SNR</mi><mi mathvariant="italic">LT</mi></msub><mo>&lt;</mo><mn>15</mn></math><img file="EP2162880B1_D0026.tif" /></maths>

0056The hangover period starts in the first inactive sound frame after three (3) consecutive active sound frames. Its function consists of forcing every inactive frame during the hangover period as an active frame. The SAD decision will be explained later.

0057For narrowband signals, the hysteresis strategy consists of decreasing the SAD decision threshold as follows: <maths id="math0027" num=""><math display="block"><msub><mi mathvariant="italic">th</mi><mi mathvariant="italic">SAD</mi></msub><mo>=</mo><msub><mi mathvariant="italic">th</mi><mi mathvariant="italic">SAD</mi></msub><mo>-</mo><mn>5.2</mn><mspace width="3em" /><mi>if</mi><mspace width="1em" /><msub><mi mathvariant="italic">SNR</mi><mi mathvariant="italic">LT</mi></msub><mo>&lt;</mo><mn>19</mn></math><img file="EP2162880B1_D0027.tif" /></maths><maths id="math0028" num=""><math display="block"><msub><mi mathvariant="italic">th</mi><mi mathvariant="italic">SAD</mi></msub><mo>=</mo><msub><mi mathvariant="italic">th</mi><mi mathvariant="italic">SAD</mi></msub><mo>-</mo><mn>2</mn><mspace width="3em" /><mi>if</mi><mspace width="1em" /><mn>19</mn><mo>≤</mo><msub><mi mathvariant="italic">SNR</mi><mi mathvariant="italic">LT</mi></msub><mo>&lt;</mo><mn>35</mn></math><img file="EP2162880B1_D0028.tif" /></maths><maths id="math0029" num=""><math display="block"><msub><mi mathvariant="italic">th</mi><mi mathvariant="italic">SAD</mi></msub><mo>=</mo><msub><mi mathvariant="italic">th</mi><mi mathvariant="italic">SAD</mi></msub><mspace width="3em" /><mi>if</mi><mspace width="1em" /><mn>35</mn><mo>≤</mo><msub><mi mathvariant="italic">SNR</mi><mi mathvariant="italic">LT</mi></msub></math><img file="EP2162880B1_D0029.tif" /></maths> Thus, for noisy signals with low SNR, the threshold becomes lower to give preference to active signal decision. There is no hangover for narrowband signals.

0058Finally, the sound activity detector 103 has two outputs - a SAD flag and a local SAD flag. Both flags are set to one if active signal is detected and set to zero otherwise. Moreover, the SAD flag is set to one in hangover period. The SAD decision is done by comparing the average SNR per frame with the SAD decision threshold (via a comparator for example), that is: <pre listing-type="program-listing"> if SNRav > thSAD SADlocal = 1 SAD = 1 else SADlocal = 0 if in hangover period SAD = 1 else SAD = 0 end end.</pre>

<i>First level of noise estimation and update</i>

0059A noise estimator 104 as illustrated in <figref idref="f0001">Figure 1</figref> calculates the total noise energy, relative frame energy, update of long-term average noise energy and long-term average frame energy, average energy per critical band, and a noise correction factor. Further, the noise estimator 104 performs noise energy initialization and update downwards.

0060The total noise energy per frame is calculated using the following relation: <maths id="math0030" num="(11)"><math display="block"><msub><mi mathvariant="italic">N</mi><mi mathvariant="italic">tot</mi></msub><mo>=</mo><mn>10</mn><mo>⁢</mo><mi>log</mi><mfenced separators=""><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mn>19</mn></munderover></mstyle><msub><mi mathvariant="italic">N</mi><mi mathvariant="italic">CB</mi></msub><mfenced><mi>i</mi></mfenced></mfenced></math><img file="EP2162880B1_D0030.tif" /></maths> where <i>N_CB(i)</i> is the estimated noise energy per critical band.

0061The relative energy of the frame is given by the difference between the frame energy in dB and the long-term average energy. The relative frame energy is calculated using the following relation: <maths id="math0031" num="(12)"><math display="block"><msub><mi>E</mi><mi mathvariant="italic">rel</mi></msub><mo>=</mo><msub><mi>E</mi><mi>t</mi></msub><mo>-</mo><msub><mover><mi>E</mi><mo>‾</mo></mover><mi>f</mi></msub></math><img file="EP2162880B1_D0031.tif" /></maths> where <i>E_t</i> is given in Equation (6).

0062The long-term average noise energy or the long-term average frame energy is updated in every frame. In case of active signal frames (SAD flag = 1), the long-term average frame energy is updated using the relation: <maths id="math0032" num="(13)"><math display="block"><msub><mover><mi>E</mi><mo>‾</mo></mover><mi>f</mi></msub><mo>=</mo><mn>0.99</mn><mo>⁢</mo><msub><mover><mi>E</mi><mo>‾</mo></mover><mi>f</mi></msub><mo>+</mo><mn>0.01</mn><mo>⁢</mo><msub><mi>E</mi><mi>t</mi></msub></math><img file="EP2162880B1_D0032.tif" /></maths> with initial value <i><o ostyle="single">E</o>_f</i> = 45<i>dB.</i>

0063In case of inactive speech frames (SAD flag = 0), the long-term average noise energy is updated as follows: <maths id="math0033" num="(14)"><math display="block"><msub><mover><mi>N</mi><mo>‾</mo></mover><mi>f</mi></msub><mo>=</mo><mn>0.99</mn><mo>⁢</mo><msub><mover><mi>N</mi><mo>‾</mo></mover><mi>f</mi></msub><mo>+</mo><mn>0.01</mn><mo>⁢</mo><msub><mi>N</mi><mi mathvariant="italic">tot</mi></msub></math><img file="EP2162880B1_D0033.tif" /></maths>

0064The initial value of <i><o ostyle="single">N</o>_f</i> is set equal to <i>N_tot</i> for the first 4 frames. Also, in the first four (4) frames, the value of <i><o ostyle="single">E</o>_f</i> is bounded by <i><o ostyle="single">E</o>_f</i> ≥ <i><o ostyle="single">N</o>_tot</i> +10.

0065The frame energy per critical band for the whole frame is computed by averaging the energies from both the first and second spectral analyses in the frame using the following relation: <maths id="math0034" num="(15)"><math display="block"><msub><mover><mi>E</mi><mo>‾</mo></mover><mi mathvariant="italic">CB</mi></msub><mfenced><mi>i</mi></mfenced><mo>=</mo><mn>0.5</mn><mo>⁢</mo><msubsup><mi>E</mi><mi mathvariant="italic">CB</mi><mfenced><mn>1</mn></mfenced></msubsup><mfenced><mi>i</mi></mfenced><mo>+</mo><mn>0.5</mn><mo>⁢</mo><msubsup><mi>E</mi><mi mathvariant="italic">CB</mi><mfenced><mn>2</mn></mfenced></msubsup><mfenced><mi>i</mi></mfenced></math><img file="EP2162880B1_D0034.tif" /></maths>

0066The noise energy per critical band <i>N_CB</i>(<i>i</i>) is initialized to 0.03.

0067At this stage, only noise energy update downward is performed for the critical bands whereby the energy is less than the background noise energy. First, the temporary updated noise energy is computed using the following relation: <maths id="math0035" num="(18)"><math display="block"><msub><mi>N</mi><mi mathvariant="italic">tmp</mi></msub><mfenced><mi>i</mi></mfenced><mo>=</mo><mn>0.9</mn><mo>⁢</mo><msub><mi>N</mi><mi>CB</mi></msub><mfenced><mi>i</mi></mfenced><mo>+</mo><mn>0.1</mn><mo>⁢</mo><mfenced separators=""><mn>0.25</mn><mo>⁢</mo><msubsup><mi>E</mi><mi mathvariant="italic">CB</mi><mfenced><mn>0</mn></mfenced></msubsup><mfenced><mi>i</mi></mfenced><mo>+</mo><mn>0.75</mn><mo>⁢</mo><msub><mover><mi>E</mi><mo>‾</mo></mover><mi mathvariant="italic">CB</mi></msub><mfenced><mi>i</mi></mfenced></mfenced></math><img file="EP2162880B1_D0035.tif" /></maths> where <maths id="math0036" num=""><math display="inline"><msubsup><mi>E</mi><mi mathvariant="italic">CB</mi><mfenced><mn>0</mn></mfenced></msubsup><mfenced><mi>i</mi></mfenced></math><img file="EP2162880B1_D0036.tif" /></maths>denotes the energy per critical band corresponding to the second spectral analysis from the previous frame.

0068Then for <i>i</i>=0 to 19, if <i>N_tmp</i>(<i>i</i>) <i>&lt; N_CB</i>(<i>i</i>) then <i>N_CB</i>(<i>i</i>) <i>= N_tmp</i>(<i>i</i>).

0069A second level of noise estimation and update is performed later by setting <i>N_CB</i>(<i>i</i>) = <i>N_tmp</i>(<i>i</i>) if the frame is declared as an inactive frame.

<i>Second level of noise estimation and update</i>

0070The parametric sound activity detection and noise estimation update module 107 updates the noise energy estimates per critical band to be used in the sound activity detector 103 in the next frame. The update is performed during inactive signal periods. However, the SAD decision performed above, which is based on the SNR per critical band, is not used for determining whether the noise energy estimates are updated. Another decision is performed based on other parameters rather independent of the SNR per critical band. The parameters used for the update of the noise energy estimates are: pitch stability, signal non-stationarity, voicing, and ratio between the 2^nd order and 16^th order LP residual error energies and have generally low sensitivity to the noise level variations. The decision for the update of the noise energy estimates is optimized for speech signals. To improve the detection of active music signals, the following other parameters are used: spectral diversity, complementary non-stationarity, noise character and tonal stability. Music detection will be explained in detail in the following description.

0071The reason for not using the SAD decision for the update of the noise energy estimates is to make the noise estimation robust to rapidly changing noise levels. If the SAD decision was used for the update of the noise energy estimates, a sudden increase in noise level would cause an increase of SNR even for inactive signal frames, preventing the noise energy estimates to update, which in turn would maintain the SNR high in the following frames, and so on. Consequently, the update would be blocked and some other logic would be needed to resume the noise adaptation.

0072In the non-restrictive illustrative embodiment of the present invention, an open-loop pitch analysis is performed in a LP analyzer and pitch tracker module 106 in <figref idref="f0001">Figure 1</figref>) to compute three open-loop pitch estimates per frame: <i>d</i>₀, <i>d</i>₁ and <i>d</i>₂ corresponding to the first half-frame, second half-frame, and the lookahead, respectively. This procedure is well known to those of ordinary skill in the art and will not be further described in the present disclosure (e.g. VMR-WB [<nplcit id="ncit0006" npl-type="s" url="http://www.3gpp2.org"><text>Source-Controlled Variable-Rate Multimode Wideband Speech Codec (VMR-WB), Service Options 62 and 63 for Spread Spectrum Systems, 3GPP2 Technical Specification C.S0052-A v1.0, April 2005 (http://www.3gpp2.org</text></nplcit>)]). The LP analyzer and pitch tracker module 106 calculates a pitch stability counter using the following relation: <maths id="math0037" num="(19)"><math display="block"><mi mathvariant="italic">pc</mi><mo>=</mo><mfenced open="|" close="|" separators=""><msub><mi>d</mi><mn>0</mn></msub><mo>-</mo><msub><mi>d</mi><mrow><mo>-</mo><mn>1</mn></mrow></msub></mfenced><mo>+</mo><mfenced open="|" close="|" separators=""><msub><mi>d</mi><mn>1</mn></msub><mo>-</mo><msub><mi>d</mi><mn>0</mn></msub></mfenced><mo>+</mo><mfenced open="|" close="|" separators=""><msub><mi>d</mi><mn>2</mn></msub><mo>-</mo><msub><mi>d</mi><mn>1</mn></msub></mfenced></math><img file="EP2162880B1_D0037.tif" /></maths> where <i>d</i>_-1 is the lag of the second half-frame of the previous frame. For pitch lags larger than 122, the LP analyzer and pitch tracker module 106 sets <i>d</i>₂<i>= d</i>₁<i>.</i> Thus, for such lags the value of <i>pc</i> in equation (19) is multiplied by 3/2 to compensate for the missing third term in the equation. The pitch stability is true if the value of <i>pc</i> is less than 14. Further, for frames with low voicing, <i>pc</i> is set to 14 to indicate pitch instability. More specifically: <maths id="math0038" num="(20)"><math display="block"><mi>If</mi><mspace width="1em" /><mfenced separators=""><msub><mi>C</mi><mi mathvariant="italic">norm</mi></msub><mfenced><msub><mi>d</mi><mn>0</mn></msub></mfenced><mo>+</mo><msub><mi>C</mi><mi mathvariant="italic">norm</mi></msub><mfenced><msub><mi>d</mi><mn>1</mn></msub></mfenced><mo>+</mo><msub><mi>C</mi><mi mathvariant="italic">norm</mi></msub><mfenced><msub><mi>d</mi><mn>2</mn></msub></mfenced></mfenced><mo>/</mo><mn>3</mn><mo>+</mo><msub><mi>r</mi><mi>e</mi></msub><mo>&lt;</mo><msub><mi mathvariant="italic">th</mi><mi mathvariant="italic">Cpc</mi></msub><mspace width="2em" /><mi>then</mi><mspace width="2em" /><mi mathvariant="italic">pc</mi><mo>=</mo><mn>14</mn><mo>,</mo></math><img file="EP2162880B1_D0038.tif" /></maths> where <i>C_norm</i>(<i>d</i>) is the normalized raw correlation and <i>r_e</i> is an optional correction added to the normalized correlation in order to compensate for the decrease of normalized correlation in the presence of background noise. The voicing threshold <i>th_Cpc</i> = 0.52 for WB, and <i>th_Cpc =</i> 0.65 for NB. The correction factor can be calculated using the following relation: <maths id="math0039" num=""><math display="block"><msub><mi>r</mi><mi>e</mi></msub><mo>=</mo><mn>0.00024492</mn><mo>⁢</mo><msup><mi>e</mi><mrow><mn>0.1596</mn><mo>⁢</mo><mfenced separators=""><msub><mi>N</mi><mi mathvariant="italic">tot</mi></msub><mo>-</mo><mn>14</mn></mfenced></mrow></msup><mo>-</mo><mn>0.022</mn></math><img file="EP2162880B1_D0039.tif" /></maths> where <i>N</i>_tot is the total noise energy per frame computed according to Equation (11).

0073The normalized raw correlation can be computed based on the decimated weighted sound signal <i>s_wd</i>(<i>n</i>) using the following equation: <maths id="math0040" num=""><math display="block"><msub><mi>C</mi><mi mathvariant="italic">norm</mi></msub><mfenced><mi>d</mi></mfenced><mo>=</mo><mfrac><mrow><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><msub><mi>L</mi><mi>sec</mi></msub></munderover></mstyle><msub><mi>s</mi><mi mathvariant="italic">wd</mi></msub><mfenced><msub><mi>t</mi><mi mathvariant="italic">start</mi></msub></mfenced><mo>⁢</mo><msub><mi>s</mi><mi mathvariant="italic">wd</mi></msub><mo>⁢</mo><mfenced separators=""><msub><mi>t</mi><mi mathvariant="italic">start</mi></msub><mo>-</mo><mi>d</mi></mfenced></mrow><msqrt><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><msub><mi>L</mi><mi>sec</mi></msub></munderover></mstyle><msubsup><mi>s</mi><mi mathvariant="italic">wd</mi><mn>2</mn></msubsup><mfenced><msub><mi>t</mi><mi mathvariant="italic">start</mi></msub></mfenced><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><msub><mi>L</mi><mi>sec</mi></msub></munderover></mstyle><msubsup><mi>s</mi><mi mathvariant="italic">wd</mi><mn>2</mn></msubsup><mo>⁢</mo><mfenced separators=""><msub><mi>t</mi><mi mathvariant="italic">start</mi></msub><mo>-</mo><mi>d</mi></mfenced></msqrt></mfrac><mo>,</mo></math><img file="EP2162880B1_D0040.tif" /></maths> where the summation limit depends on the delay itself. The weighted signal <i>s_wd</i>(<i>n</i>) is the one used in open-loop pitch analysis and given by filtering the pre-processed input sound signal from pre-processor 101 through a weighting filter of the form <i>A</i>(<i>z</i>/<i>γ</i>)/(1<i>-µz</i>^-1)<i>.</i> The weighted signal <i>s_wd</i>(<i>n</i>) is decimated by 2 and the summation limits are given according to: <ul id="ul0003" list-style="none" compact="compact"><li><i>L</i>_sec = 40 for <i>d</i> =10,...,16</li><li><i>L</i>_sec = 40 for <i>d</i> =17, ..., 31</li><li><i>L</i>_sec = 62 for <i>d</i> = 32, ..., 61</li><li><i>L</i>_sec =115 for <i>d</i> = 62, ...,115</li></ul>

0074These lengths assure that the correlated vector length comprises at least one pitch period which helps to obtain a robust open-loop pitch detection. The instants <i>t_start</i> are related to the current frame beginning and are given by: <ul id="ul0004" list-style="none" compact="compact"><li><i>t_start</i> = 0 for first half-frame</li><li><i>t_start</i> = 128 for second half-frame</li><li><i>t_start</i> = 256 for look-ahead</li></ul> at 12.8 kHz sampling rate.

0075The parametric sound activity detection and noise estimation update module 107 performs a signal non-stationarity estimation based on the product of the ratios between the energy per critical band and the average long term energy per critical band.

0076The average long term energy per critical band is updated using the following relation: <maths id="math0041" num="(21)"><math display="block"><msub><mi>E</mi><mrow><mi mathvariant="italic">CB</mi><mo>,</mo><mi mathvariant="italic">LT</mi></mrow></msub><mfenced><mi>i</mi></mfenced><mo>=</mo><msub><mi>α</mi><mi>e</mi></msub><mo>⁢</mo><msub><mi>E</mi><mrow><mi mathvariant="italic">CB</mi><mo>,</mo><mi mathvariant="italic">LT</mi></mrow></msub><mfenced><mi>i</mi></mfenced><mo>+</mo><mfenced separators=""><mn>1</mn><mo>-</mo><msub><mi>α</mi><mi>e</mi></msub></mfenced><mo>⁢</mo><msub><mover><mi>E</mi><mo>‾</mo></mover><mi mathvariant="italic">CB</mi></msub><mfenced><mi>i</mi></mfenced><mo>,</mo><mspace width="2em" /><mi>for i</mi><mo>=</mo><msub><mi>b</mi><mi mathvariant="italic">min</mi></msub><mspace width="1em" /><mi>to</mi><mspace width="1em" /><msub><mi>b</mi><mi mathvariant="italic">max</mi></msub><mo>,</mo></math><img file="EP2162880B1_D0041.tif" /></maths> where <i>b_min</i>=0 and <i>b_max</i>=19 in the case of wideband signals, and <i>b_min</i>=1 and <i>b_max</i>=16 in case of narrowband signals, and <i><o ostyle="single">E</o>_CB(i)</i> is the frame energy per critical band defined in Equation (15). The update factor <i>α_e</i> is a linear function of the total frame energy, defined in Equation (6), and it is given as follows: <ul id="ul0005" list-style="none"><li>For wideband signals: <i>α_e</i> = 0.0245<i>E_t</i> - 0.235 bounded by 0.5 ≤ <i>α_e</i> ≤ 0.99.</li><li>For narrowband signals: <i>α_e</i> = 0.00091<i>E_t</i> + 0.3185 bounded by 0.5 ≤ <i>α_e</i> ≤ 0.999.</li></ul>

0077<i>E_t</i> is given by Equation (6).

0078The frame non-stationarity is given by the product of the ratios between the frame energy and average long term energy per critical band. More specifically: <maths id="math0042" num="(22)"><math display="block"><mi mathvariant="italic">nonstat</mi><mo>=</mo><mstyle displaystyle="true"><munderover><mo>∏</mo><mrow><mi>i</mi><mo>=</mo><msub><mi>b</mi><mi>min</mi></msub></mrow><msub><mi>b</mi><mi>max</mi></msub></munderover></mstyle><mfrac><mrow><mi>max</mi><mfenced separators=""><msub><mover><mi>E</mi><mo>‾</mo></mover><mi mathvariant="italic">CB</mi></msub><mfenced><mi>i</mi></mfenced><mo>,</mo><msub><mi>E</mi><mrow><mi mathvariant="italic">CB</mi><mo>,</mo><mi mathvariant="italic">LT</mi></mrow></msub><mfenced><mi>i</mi></mfenced></mfenced></mrow><mrow><mi>min</mi><mfenced separators=""><msub><mover><mi>E</mi><mo>‾</mo></mover><mi mathvariant="italic">CB</mi></msub><mfenced><mi>i</mi></mfenced><mo>,</mo><msub><mi>E</mi><mrow><mi mathvariant="italic">CB</mi><mo>,</mo><mi mathvariant="italic">LT</mi></mrow></msub><mfenced><mi>i</mi></mfenced></mfenced></mrow></mfrac></math><img file="EP2162880B1_D0042.tif" /></maths>

0079The parametric sound activity detection and noise estimation update module 107 further produces a voicing factor for noise update using the following relation: <maths id="math0043" num="(23)"><math display="block"><mi mathvariant="italic">voicing</mi><mo>=</mo><mfenced separators=""><msub><mi>C</mi><mi mathvariant="italic">norm</mi></msub><mfenced><msub><mi>d</mi><mn>0</mn></msub></mfenced><mo>+</mo><msub><mi>C</mi><mi mathvariant="italic">norm</mi></msub><mfenced><msub><mi>d</mi><mn>1</mn></msub></mfenced></mfenced><mo>/</mo><mn>2</mn><mo>+</mo><msub><mi>r</mi><mi>e</mi></msub></math><img file="EP2162880B1_D0043.tif" /></maths>

0080Finally, the parametric sound activity detection and noise estimation update module 107 calculates a ratio between the LP residual energy after the 2^nd order and 16^th order LP analysis using the relation: <maths id="math0044" num="(24)"><math display="block"><mi mathvariant="italic">resid_ratio</mi><mo>=</mo><mi>E</mi><mfenced><mn>2</mn></mfenced><mo>/</mo><mi>E</mi><mfenced><mn>16</mn></mfenced></math><img file="EP2162880B1_D0044.tif" /></maths> where E(2) and <i>E</i>(16) are the LP residual energies after 2^nd order and 16^th order LP analysis as computed in the LP analyzer and pitch tracker module 106 using a Levinson-Durbin recursion which is a procedure well known to those of ordinary skill in the art. This ratio reflects the fact that to represent a signal spectral envelope, a higher order of LP is generally needed for speech signal than for noise. In other words, the difference between E(2) and <i>E</i>(16) is supposed to be lower for noise than for active speech.

0081The update decision made by the parametric sound activity detection and noise estimation update module 107 is determined based on a variable <i>noise</i>_<i>update</i> which is initially set to 6 and is decreased by 1 if an inactive frame is detected and incremented by 2 if an active frame is detected. Also, the variable <i>noise</i>_<i>update</i> is bounded between 0 and 6. The noise energy estimates are updated only when <i>noise_update</i>=0<i>.</i>

0082The value of the variable <i>noise_update</i> is updated in each frame as follows: <pre listing-type="program-listing"> If (nonstat > thstat) OR (pc < 14) OR (voicing > thCnorm) OR (resid_ratio > thresid) noise_update = noise_update + 2 Else noise_update = noise_update - 1</pre> where for wideband signals, <i>th_stat</i> = <i>th_Cnorm =</i> 0.85 and <i>th_resid</i> = 1.6, and for narrowband signals, <i>th_stat</i> = 500000, <i>th_Cnorm</i> = 0.7 and <i>th_resid</i>= 10.4.

0083In other words, frames are declared inactive for noise update when <maths id="math0045" num=""><math display="block"><mfenced separators=""><mi mathvariant="italic">nonstat</mi><mo>≤</mo><msub><mi mathvariant="italic">th</mi><mi mathvariant="italic">stat</mi></msub></mfenced><mspace width="1em" /><mi>AND</mi><mspace width="1em" /><mfenced separators=""><mi mathvariant="italic">pc</mi><mo>≥</mo><mn>14</mn></mfenced><mspace width="1em" /><mi>AND</mi><mspace width="1em" /><mfenced separators=""><mi mathvariant="italic">voicing</mi><mo>≤</mo><msub><mi mathvariant="italic">th</mi><mi mathvariant="italic">Cnorm</mi></msub></mfenced><mspace width="1em" /><mi>AND</mi><mspace width="1em" /><mfenced separators=""><mi mathvariant="italic">resid_ratio</mi><mo>≤</mo><msub><mi mathvariant="italic">th</mi><mi mathvariant="italic">resid</mi></msub></mfenced></math><img file="EP2162880B1_D0045.tif" /></maths> and a hangover of 6 frames is used before noise update takes place.

0084Thus, if <i>noise_update</i>=0 then for i=0 to 19 <i>N_CB</i>(<i>i</i>) = <i>N_tmp</i>(<i>i</i>) where <i>N_tmp (i)</i> is the temporary updated noise energy already computed in Equation (18).

Improvement of noise detection for music signals

0085The noise estimation described above has its limitations for certain music signals, such as piano concerts or instrumental rock and pop, because it was developed and optimized mainly for speech detection. To improve the detection of music signals in general, the parametric sound activity detection and noise estimation update module 107 uses other parameters or techniques in conjunction with the existing ones. These other parameters or techniques comprise, as described hereinabove, spectral diversity, complementary non-stationarity, noise character and tonal stability, which are calculated by a spectral diversity calculator, a complementary non-stationarity calculator, a noise character calculator and a tonality estimator, respectively. They will be described in detail herein below.

Spectral diversity

0086Spectral diversity gives information about significant changes of the signal in frequency domain. The changes are tracked in critical bands by comparing energies in the first spectral analysis of the current frame and the second spectral analysis two frames ago. The energy in a critical band i of the first spectral analysis in the current frame is denoted as <maths id="math0046" num=""><math display="inline"><msubsup><mi>E</mi><mi mathvariant="italic">CB</mi><mfenced><mn>1</mn></mfenced></msubsup><mfenced><mi>i</mi></mfenced><mn>.</mn></math><img file="EP2162880B1_D0046.tif" /></maths>Let the energy in the same critical band calculated in the second spectral analysis two frames ago be denoted as <maths id="math0047" num=""><math display="inline"><msubsup><mi>E</mi><mi mathvariant="italic">CB</mi><mfenced separators=""><mo>-</mo><mn>2</mn></mfenced></msubsup><mfenced><mi>i</mi></mfenced><mn>.</mn></math><img file="EP2162880B1_D0047.tif" /></maths> Both of these energies are initialized to 0.0001. Then, for all critical bands higher than 9, the maximum and the minimum of the two energies are calculated as follows: <maths id="math0048" num=""><math display="block"><mtable columnalign="left"><mtr><mtd><msub><mi>E</mi><mi>max</mi></msub><mfenced><mi>i</mi></mfenced><mo>=</mo><mi>max</mi><mfenced open="{" close="}" separators=""><msubsup><mi>E</mi><mi mathvariant="italic">CB</mi><mfenced><mn>1</mn></mfenced></msubsup><mfenced><mi>i</mi></mfenced><mo>,</mo><msubsup><mi>E</mi><mi mathvariant="italic">CB</mi><mfenced separators=""><mo>-</mo><mn>2</mn></mfenced></msubsup><mfenced><mi>i</mi></mfenced></mfenced></mtd></mtr><mtr><mtd><msub><mi>E</mi><mi>min</mi></msub><mfenced><mi>i</mi></mfenced><mo>=</mo><mi>min</mi><mfenced open="{" close="}" separators=""><msubsup><mi>E</mi><mi mathvariant="italic">CB</mi><mfenced><mn>1</mn></mfenced></msubsup><mfenced><mi>i</mi></mfenced><mo>,</mo><msubsup><mi>E</mi><mi mathvariant="italic">CB</mi><mfenced separators=""><mo>-</mo><mn>2</mn></mfenced></msubsup><mfenced><mi>i</mi></mfenced></mfenced></mtd></mtr></mtable><mo>,</mo><mspace width="1em" /><mi>for i</mi><mo>=</mo><mn>10</mn><mo>,</mo><mo>…</mo><mo>,</mo><msub><mi>b</mi><mi>max</mi></msub><mn>.</mn></math><img file="EP2162880B1_D0048.tif" /></maths> Subsequently, a ratio between the maximum and the minimum energy in a specific critical band is calculated as <maths id="math0049" num=""><math display="block"><msub><mi>E</mi><mi mathvariant="italic">rat</mi></msub><mfenced><mi>i</mi></mfenced><mo>=</mo><mfrac><mrow><msub><mi>E</mi><mi>max</mi></msub><mfenced><mi>i</mi></mfenced></mrow><mrow><msub><mi>E</mi><mi>min</mi></msub><mfenced><mi>i</mi></mfenced></mrow></mfrac><mo>,</mo><mspace width="1em" /><mi>for i</mi><mo>=</mo><mn>10</mn><mo>,</mo><mo>…</mo><mo>,</mo><msub><mi>b</mi><mi>max</mi></msub><mn>.</mn></math><img file="EP2162880B1_D0049.tif" /></maths>

0087Finally, the parametric sound activity detection and noise estimation update module 107 calculates a spectral diversity parameter as a normalized weighted sum of the ratios with the weight itself being the maximum energy <i>E</i>_max(<i>i</i>). This spectral diversity parameter is given by the following relation: <maths id="math0050" num="(25)"><math display="block"><mi mathvariant="italic">spec_div</mi><mo>=</mo><mfrac><mrow><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mn>10</mn></mrow><msub><mi>b</mi><mi>max</mi></msub></munderover></mstyle><msub><mi>E</mi><mi>max</mi></msub><mfenced><mi>i</mi></mfenced><mo>⁢</mo><msub><mi>E</mi><mi mathvariant="italic">rat</mi></msub><mfenced><mi>i</mi></mfenced></mrow><mrow><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mn>10</mn></mrow><msub><mi>b</mi><mi>max</mi></msub></munderover></mstyle><msub><mi>E</mi><mi>max</mi></msub><mfenced><mi>i</mi></mfenced></mrow></mfrac></math><img file="EP2162880B1_D0050.tif" /></maths>

0088The <i>spec_div</i> parameter is used in the final decision about music activity and noise energy update. The <i>spec_div</i> parameter is also used as an auxiliary parameter for the calculation of a complementary non-stationarity parameter which is described bellow.

Complementary non-stationarity

0089The inclusion of a complementary non-stationarity parameter is motivated by the fact that the non-stationarity parameter, defined in Equation (22), fails when a sharp energy attack in a music signal is followed by a slow energy decrease. In this case the average long term energy per critical band, <i>E_CB,LT</i>(<i>i</i>), defined in Equation (21), slowly increases during the attack whereas the frame energy per critical band, defined in Equation (15), slowly decreases. In a certain frame after the attack these two energy values meet and the <i>nonstat</i> parameter results in a small value indicating an absence of active signal. This leads to a false noise update and subsequently a false SAD decision.

0090To overcome this problem an alternative average long term energy per critical band is calculated using the following relation: <maths id="math0051" num="(26)"><math display="block"><msub><mrow><mi>E</mi><mo>⁢</mo><mn mathvariant="italic">2</mn></mrow><mrow><mi mathvariant="italic">CB</mi><mo>,</mo><mi mathvariant="italic">LT</mi></mrow></msub><mfenced><mi>i</mi></mfenced><mo>=</mo><msub><mi>β</mi><mi>e</mi></msub><mo>⁢</mo><msub><mrow><mi>E</mi><mo>⁢</mo><mn mathvariant="italic">2</mn></mrow><mrow><mi mathvariant="italic">CB</mi><mo>,</mo><mi mathvariant="italic">LT</mi></mrow></msub><mfenced><mi>i</mi></mfenced><mo>+</mo><mfenced separators=""><mn>1</mn><mo>-</mo><msub><mi>β</mi><mi>e</mi></msub></mfenced><mo>⁢</mo><msub><mover><mi>E</mi><mo>‾</mo></mover><mi mathvariant="italic">CB</mi></msub><mfenced><mi>i</mi></mfenced><mo>,</mo><mspace width="2em" /><mi>for i</mi><mo>=</mo><msub><mi>b</mi><mi mathvariant="italic">min</mi></msub><mspace width="1em" /><mi>to</mi><mspace width="1em" /><msub><mi>b</mi><mi mathvariant="italic">max</mi></msub><mn>.</mn></math><img file="EP2162880B1_D0051.tif" /></maths>

0091The variable <i>E2_CB,LT</i>(<i>i</i>) is initialized to 0.03 for all i. Equation (26) closely resembles equation (21) with the only difference being the update factor <i>β_e</i> which is given as follows: <pre listing-type="program-listing"> if (spec_div > thspec_div) βe = 0 else βe = αe end,</pre> where <i>th_{spec_div}</i> = 5. Thus, when an energy attack is detected (<i>spec_div ></i> 5) the alternative average long term energy is immediately set to the average frame energy, i.e. <i>E2_CB,LT (i) = <o ostyle="single">E</o>_CB</i>(<i>i</i>)<i>.</i> Otherwise this alternative average long term energy is updated in the same way as the conventional non-stationarity, i.e. using the exponential filter with the update factor <i>α_e.</i> The complementary non-stationarity parameter is calculated in the same way as <i>nonstat,</i> but using <i>E2_CB,LT</i>(<i>i</i>), i.e. <maths id="math0052" num="(27)"><math display="block"><mi mathvariant="italic">nonstat</mi><mo>⁢</mo><mn mathvariant="italic">2</mn><mo>=</mo><mstyle displaystyle="true"><munderover><mo>∏</mo><mrow><mi>i</mi><mo>=</mo><msub><mi>b</mi><mi>min</mi></msub></mrow><msub><mi>b</mi><mi>max</mi></msub></munderover></mstyle><mfrac><mrow><mi>max</mi><mfenced separators=""><msub><mover><mi>E</mi><mo>‾</mo></mover><mi mathvariant="italic">CB</mi></msub><mfenced><mi>i</mi></mfenced><mo>,</mo><msub><mrow><mi>E</mi><mo>⁢</mo><mn mathvariant="italic">2</mn></mrow><mrow><mi mathvariant="italic">CB</mi><mo>,</mo><mi mathvariant="italic">LT</mi></mrow></msub><mfenced><mi>i</mi></mfenced></mfenced></mrow><mrow><mi>max</mi><mfenced separators=""><msub><mover><mi>E</mi><mo>‾</mo></mover><mi mathvariant="italic">CB</mi></msub><mfenced><mi>i</mi></mfenced><mo>,</mo><msub><mrow><mi>E</mi><mo>⁢</mo><mn mathvariant="italic">2</mn></mrow><mrow><mi mathvariant="italic">CB</mi><mo>,</mo><mi mathvariant="italic">LT</mi></mrow></msub><mfenced><mi>i</mi></mfenced></mfenced></mrow></mfrac></math><img file="EP2162880B1_D0052.tif" /></maths>

0092The complementary non-stationarity parameter, <i>nonstat2,</i> may fail a few frames right after an energy attack, but should not fail during the passages characterized by a slowly-decreasing energy. Since the <i>nonstat</i> parameter works well on energy attacks and few frames after, a logical disjunction of <i>nonstat</i> and <i>nonstat2</i> therefore solves the problem of inactive signal detection on certain musical signals. However, the disjunction is applied only in passages which are "likely to be active". The likelihood is calculated as follows: <pre listing-type="program-listing"> if ((nonstat > thstat) OR (tonal_stability = 1)) act_pred_LT = ka act_pred_LT + (1- ka).1 else act_pred_LT = ka act_pred_LT + (1- ka).0 end.</pre> The coefficient <i>k_a</i> is set to 0.99. The parameter <i>act_pred_LT</i> which is in the range <0:1> may be interpreted as a predictor of activity. When it is close to 1, the signal is likely to be active, and when it is close to 0, it is likely to be inactive. The <i>act_pred_LT</i> parameter is initialized to one. In the condition above, <i>tonal_stability</i> is a binary parameter which is used to detect stable tonal signal. This <i>tonal_stability</i> parameter will be described in the following description.

0093The <i>nonstat2</i> parameter is taken into consideration (in disjunction with <i>nonstat)</i> in the update of noise energy only if <i>act_pred_LT</i> is higher than certain threshold, which has been set to 0.8. The logic of noise energy update is explained in detail at the end of the present section.

Noise character

0094Noise character is another parameter which is used in the detection of certain noise-like music signals such as cymbals or low-frequency drums. This parameter is calculated using the following relation: <maths id="math0053" num="(28)"><math display="block"><mi mathvariant="italic">noise_char</mi><mo>=</mo><mfrac><mrow><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mn>10</mn></mrow><msub><mi>b</mi><mi>max</mi></msub></munderover></mstyle><msub><mi>E</mi><mi mathvariant="italic">CB</mi></msub><mfenced><mi>i</mi></mfenced></mrow><mrow><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><msub><mi>b</mi><mi>min</mi></msub></mrow><mn>9</mn></munderover></mstyle><msub><mi>E</mi><mi mathvariant="italic">CB</mi></msub><mfenced><mi>i</mi></mfenced></mrow></mfrac></math><img file="EP2162880B1_D0053.tif" /></maths>

0095The <i>noise_char</i> parameter is calculated only for the frames whose spectral content has at least a minimal energy, which is fulfilled when both the numerator and the denominator of Equation (28) are larger than 100. The <i>noise</i>_<i>char</i> parameter is upper limited by 10 and its long-term value is updated using the following relation: <maths id="math0054" num="(29)"><math display="block"><mi mathvariant="italic">noise_char_LT</mi><mo>=</mo><msub><mi>α</mi><mi>n</mi></msub><mo>⁢</mo><mi mathvariant="italic">noise_char_LT</mi><mo>+</mo><mfenced separators=""><mn>1</mn><mo>-</mo><msub><mi>α</mi><mi>n</mi></msub></mfenced><mo>⁢</mo><mi mathvariant="italic">noise_char</mi></math><img file="EP2162880B1_D0054.tif" /></maths>

0096The initial value of <i>noise_char_LT</i> is 0 and <i>α_n</i> is set equal to 0.9. This <i>noise_char_LT</i> parameter is used in the decision about noise energy update which is explained at the end of the present section.

Tonal stability

0097Tonal stability is the last parameter used to prevent false update of the noise energy estimates. Tonal stability is also used to prevent declaring some music segments as unvoiced frames. Tonal stability is further used in an embedded super-wideband codec to decide which coding model will be used for encoding the sound signal above 7 kHz. Detection of tonal stability exploits the tonal nature of music signals. In a typical music signal there are tones which are stable over several consecutive frames. To exploit this feature, it is necessary to track the positions and shapes of strong spectral peaks since these may correspond to the tones. The tonal stability detection is based on a correlation analysis between the spectral peaks in the current frame and those of the past frame. The input is the average log-energy spectrum defined in Equation (4). The number of spectral bins is denoted as <i>N_SPEC</i> (bin 0 is the DC component and <i>N_SPEC</i>=<i>L_FFT</i>/2). In the following disclosure, the term "spectrum" will refer to the average log-energy spectrum, as defined by Equation (4).

0098Detection of tonal stability proceeds in three stages. Furthermore, detection of tonal stability uses a calculator of a current residual spectrum, a detector of peaks in the current residual spectrum and a calculator of a correlation map and a long-term correlation map, which will be described hereinabelow.

0099In the first stage, the indexes of local minima of the spectrum are searched (by a spectrum minima locator for example), in a loop described by the following formula and stored in a buffer <i>i_min</i> that can be expressed as follows: <maths id="math0055" num="(30)"><math display="block"><msub><mi>i</mi><mi>min</mi></msub><mo>=</mo><mrow><mo>(</mo><mo>∀</mo><mi>i</mi><mo>:</mo><mfenced separators=""><msub><mi>E</mi><mi mathvariant="italic">dB</mi></msub><mo>⁢</mo><mfenced separators=""><mi>i</mi><mo>-</mo><mn>1</mn></mfenced><mo>></mo><msub><mi>E</mi><mi mathvariant="italic">dB</mi></msub><mfenced><mi>i</mi></mfenced></mfenced><mo>∧</mo><mfenced separators=""><msub><mi>E</mi><mi mathvariant="italic">dB</mi></msub><mfenced><mi>i</mi></mfenced><mo><</mo><msub><mi>E</mi><mi mathvariant="italic">dB</mi></msub><mo>⁢</mo><mfenced separators=""><mi>i</mi><mo>+</mo><mn>1</mn></mfenced></mfenced></mrow><mspace width="2em" /><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>…</mo><mo>,</mo><msub><mi>N</mi><mi mathvariant="italic">SPEC</mi></msub><mo>-</mo><mn>2</mn></math><img file="EP2162880B1_D0055.tif" /></maths> where the symbol ∧ means logical AND.

0100In Equation (30), <i>E_dB</i>(<i>i</i>) denotes the average log-energy spectrum calculated through Equation (4). The first index in <i>i</i>_min is 0, if <i>E_dB</i>(0) &lt; <i>E_dB</i>(1). Consequently, the last index in <i>i</i>_min is <i>N_SPEC</i>-1<i>,</i> if <i>E_dB</i>(<i>N_SPEC</i> -1) <i>&lt; E_dB</i>(<i>N_SPEC</i> - 2). Let us denote the number of minima found as <i>N_min.</i>

0101The second stage consists of calculating a spectral floor (through a spectral floor estimator for example) and subtracting it from the spectrum (via a suitable subtractor for example). The spectral floor is a piece-wise linear function which runs through the detected local minima. Every linear piece between two consecutive minima <i>i</i>_min (x) and <i>i</i>_min (<i>x</i> + 1) can be described as: <maths id="math0056" num=""><math display="block"><mi mathvariant="italic">fl</mi><mfenced><mi>j</mi></mfenced><mo>=</mo><mi>k</mi><mn>.</mn><mfenced separators=""><mi>j</mi><mo>-</mo><msub><mi>i</mi><mi>min</mi></msub><mfenced><mi>x</mi></mfenced></mfenced><mo>+</mo><mi>q</mi><mspace width="8em" /><mi>j</mi><mo>=</mo><msub><mi>i</mi><mi>min</mi></msub><mfenced><mi>x</mi></mfenced><mo>,</mo><mo>…</mo><mo>,</mo><msub><mi>i</mi><mi>min</mi></msub><mo>⁢</mo><mfenced separators=""><mi>x</mi><mo>+</mo><mn>1</mn></mfenced><mo>,</mo></math><img file="EP2162880B1_D0056.tif" /></maths> where <i>k</i> is the slope of the line and <i>q = E_dB</i>(<i>i</i>_min(<i>x</i>)). The slope <i>k</i> can be calculated using the following relation: <maths id="math0057" num=""><math display="block"><mi>k</mi><mo>=</mo><mfrac><mrow><msub><mi>E</mi><mi mathvariant="italic">dB</mi></msub><mo>⁢</mo><mfenced separators=""><msub><mi>i</mi><mi>min</mi></msub><mo>⁢</mo><mfenced separators=""><mi>x</mi><mo>+</mo><mn>1</mn></mfenced></mfenced><mo>-</mo><msub><mi>E</mi><mi mathvariant="italic">dB</mi></msub><mfenced separators=""><msub><mi>i</mi><mi>min</mi></msub><mfenced><mi>x</mi></mfenced></mfenced></mrow><mrow><msub><mi>i</mi><mi>min</mi></msub><mo>⁢</mo><mfenced separators=""><mi>x</mi><mo>+</mo><mn>1</mn></mfenced><mo>-</mo><msub><mi>i</mi><mi>min</mi></msub><mfenced><mi>x</mi></mfenced></mrow></mfrac><mn>.</mn></math><img file="EP2162880B1_D0057.tif" /></maths>

0102Thus, the spectral floor is a logical connection of all pieces: <maths id="math0058" num="(31)"><math display="block"><mtable columnalign="left"><mtr><mtd><mi mathvariant="italic">sp_floor</mi><mfenced><mi>j</mi></mfenced><mo>=</mo><msub><mi>E</mi><mi mathvariant="italic">dB</mi></msub><mfenced><mi>j</mi></mfenced><mspace width="3em" /></mtd><mtd><mi>j</mi><mo>=</mo><mn>0</mn><mo>,</mo><mo>…</mo><mo>,</mo><msub><mi>i</mi><mi>min</mi></msub><mfenced><mn>0</mn></mfenced><mo>-</mo><mn>1</mn></mtd></mtr><mtr><mtd><mi mathvariant="italic">sp_floor</mi><mfenced><mi>j</mi></mfenced><mo>=</mo><mi mathvariant="italic">fl</mi><mfenced><mi>j</mi></mfenced></mtd><mtd><mi>j</mi><mo>=</mo><msub><mi>i</mi><mi>min</mi></msub><mfenced><mn>0</mn></mfenced><mo>,</mo><mo>…</mo><mo>,</mo><msub><mi>i</mi><mi>min</mi></msub><mo>⁢</mo><mfenced separators=""><msub><mi>N</mi><mi>min</mi></msub><mo>-</mo><mn>1</mn></mfenced><mo>-</mo><mn>1.</mn></mtd></mtr><mtr><mtd><mi mathvariant="italic">sp_floor</mi><mfenced><mi>j</mi></mfenced><mo>=</mo><msub><mi>E</mi><mi mathvariant="italic">dB</mi></msub><mfenced><mi>j</mi></mfenced></mtd><mtd><mi>j</mi><mo>=</mo><msub><mi>i</mi><mi>min</mi></msub><mo>⁢</mo><mfenced separators=""><msub><mi>N</mi><mi>min</mi></msub><mo>-</mo><mn>1</mn></mfenced><mo>,</mo><mo>…</mo><mo>,</mo><msub><mi>N</mi><mi mathvariant="italic">SPEC</mi></msub><mo>-</mo><mn>1</mn></mtd></mtr></mtable></math><img file="EP2162880B1_D0058.tif" /></maths>

0103The leading bins up to <i>i</i>_min (0) and the terminating bins from <i>i</i>_min (<i>N</i>_min - 1) of the spectral floor are set to the spectrum itself. Finally, the spectral floor is subtracted from the spectrum using the following relation: <maths id="math0059" num="(32)"><math display="block"><msub><mi>E</mi><mrow><mi mathvariant="italic">dB</mi><mo>,</mo><mi mathvariant="italic">res</mi></mrow></msub><mfenced><mi>j</mi></mfenced><mo>=</mo><msub><mi>E</mi><mi mathvariant="italic">dB</mi></msub><mfenced><mi>j</mi></mfenced><mo>-</mo><mi mathvariant="italic">sp_floor</mi><mfenced><mi>j</mi></mfenced><mspace width="4em" /><mi>j</mi><mo>=</mo><mn>0</mn><mo>,</mo><mo>…</mo><mo>,</mo><msub><mi>N</mi><mi mathvariant="italic">SPEC</mi></msub><mo>-</mo><mn>1</mn></math><img file="EP2162880B1_D0059.tif" /></maths> and the result is called the residual spectrum. The calculation of the spectral floor is illustrated in <figref idref="f0003">Figure 3</figref>.

0104In the third stage, a correlation map and a long-term correlation map are calculated from the residual spectrum of the current and the previous frame. This is again a piece-wise operation. Thus, the correlation map is calculated on a peak-by-peak basis since the minima delimit the peaks. In the following disclosure, the term "peak" will be used to denote a piece between two minima in the residual spectrum <i>E_db,res.</i>

0105Let us denote the residual spectrum of the previous frame as <maths id="math0060" num=""><math display="inline"><msubsup><mi>E</mi><mrow><mi mathvariant="italic">dB</mi><mo>,</mo><mi mathvariant="italic">res</mi></mrow><mfenced separators=""><mo>-</mo><mn>1</mn></mfenced></msubsup><mfenced><mi>j</mi></mfenced><mn>.</mn></math><img file="EP2162880B1_D0060.tif" /></maths>For every peak in the current residual spectrum a normalized correlation is calculated with the shape in the previous residual spectrum corresponding to the position of this peak. If the signal was stable, the peaks should not move significantly from frame to frame and their positions and shapes should be approximately the same. Thus, the correlation operation takes into account all indexes (bins) of a specific peak, which is delimited by two consecutive minima. More specifically, the normalized correlation is calculated using the following relation: <maths id="math0061" num="(33)"><math display="block"><mtable columnalign="left"><mtr><mtd><mi mathvariant="italic">cor_map</mi><mo>⁢</mo><mfenced separators=""><msub><mi>i</mi><mi>min</mi></msub><mfenced><mi>x</mi></mfenced><mo>:</mo><msub><mi>i</mi><mi>min</mi></msub><mo>⁢</mo><mfenced separators=""><mi>x</mi><mo>+</mo><mn>1</mn></mfenced></mfenced><mo>=</mo><mfrac><msup><mfenced separators=""><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>j</mi><mo>=</mo><msub><mi>i</mi><mi>min</mi></msub><mfenced><mi>x</mi></mfenced></mrow><mrow><msub><mi>i</mi><mi>min</mi></msub><mo>⁢</mo><mfenced separators=""><mi>x</mi><mo>+</mo><mn>1</mn></mfenced><mo>-</mo><mn>1</mn></mrow></munderover></mstyle><msub><mi>E</mi><mrow><mi mathvariant="italic">dB</mi><mo>,</mo><mi mathvariant="italic">res</mi></mrow></msub><mfenced><mi>j</mi></mfenced><mo>⁢</mo><msubsup><mi>E</mi><mrow><mi mathvariant="italic">dB</mi><mo>,</mo><mi mathvariant="italic">res</mi></mrow><mfenced separators=""><mo>-</mo><mn>1</mn></mfenced></msubsup><mfenced><mi>j</mi></mfenced></mfenced><mn>2</mn></msup><mrow><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>j</mi><mo>=</mo><msub><mi>i</mi><mi>min</mi></msub><mfenced><mi>x</mi></mfenced></mrow><mrow><msub><mi>i</mi><mi>min</mi></msub><mo>⁢</mo><mfenced separators=""><mi>x</mi><mo>+</mo><mn>1</mn></mfenced><mo>-</mo><mn>1</mn></mrow></munderover></mstyle><msup><mfenced separators=""><msub><mi>E</mi><mrow><mi mathvariant="italic">dB</mi><mo>,</mo><mi mathvariant="italic">res</mi></mrow></msub><mfenced><mi>j</mi></mfenced></mfenced><mn>2</mn></msup><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>j</mi><mo>=</mo><msub><mi>i</mi><mi>min</mi></msub><mfenced><mi>x</mi></mfenced></mrow><mrow><msub><mi>i</mi><mi>min</mi></msub><mo>⁢</mo><mfenced separators=""><mi>x</mi><mo>+</mo><mn>1</mn></mfenced></mrow></munderover></mstyle><msup><mfenced separators=""><msubsup><mi>E</mi><mrow><mi mathvariant="italic">dB</mi><mo>,</mo><mi mathvariant="italic">res</mi></mrow><mfenced separators=""><mo>-</mo><mn>1</mn></mfenced></msubsup><mfenced><mi>j</mi></mfenced></mfenced><mn>2</mn></msup></mrow></mfrac></mtd></mtr><mtr><mtd><mi>x</mi><mo>=</mo><mn>0</mn><mo>,</mo><mo>…</mo><mo>,</mo><msub><mi>N</mi><mi>min</mi></msub><mo>-</mo><mn>2</mn></mtd></mtr></mtable><mo>,</mo></math><img file="EP2162880B1_D0061.tif" /></maths>

0106The leading bins of <i>cor_map</i> up to <i>i</i>_min(0) and the terminating bins <i>cor_map</i> from <i>i</i>_min(<i>N</i>_min -1) are set to zero. The correlation map is shown in <figref idref="f0004">Figure 4</figref>.

0107The correlation map of the current frame is used to update its long term value which is described by: <maths id="math0062" num="(34)"><math display="block"><mtable columnalign="left"><mtr><mtd><mi mathvariant="italic">cor_map_LT</mi><mfenced><mi>k</mi></mfenced><mo>=</mo><msub><mi>α</mi><mi mathvariant="italic">map</mi></msub><mo>⁢</mo><mi mathvariant="italic">cor_map_LT</mi><mfenced><mi>k</mi></mfenced><mo>+</mo><mfenced separators=""><mn>1</mn><mo>-</mo><msub><mi>α</mi><mi mathvariant="italic">map</mi></msub></mfenced><mo>⁢</mo><mi mathvariant="italic">cor_map</mi><mfenced><mi>k</mi></mfenced><mo>,</mo></mtd></mtr><mtr><mtd><mi>k</mi><mo>=</mo><mn>0</mn><mo>,</mo><mo>…</mo><mo>,</mo><msub><mi>N</mi><mi mathvariant="italic">SPEC</mi></msub><mo>-</mo><mn>1</mn><mo>,</mo></mtd></mtr></mtable></math><img file="EP2162880B1_D0062.tif" /></maths> where <i>α_map</i> = 0.9. The <i>cor_map_LT</i> is initialized to zero for all <i>k.</i>

0108Finally, all values of the <i>cor_map_LT</i> are summed together (through an adder for example) as follows: <maths id="math0063" num="(35)"><math display="block"><mi mathvariant="italic">cor_map_sum</mi><mo>=</mo><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>j</mi><mo>=</mo><mn>0</mn></mrow><mrow><msub><mi>N</mi><mi mathvariant="italic">SPEC</mi></msub><mo>-</mo><mn>1</mn></mrow></munderover></mstyle><mi mathvariant="italic">cor_map_LT</mi><mfenced><mi>j</mi></mfenced><mn>.</mn></math><img file="EP2162880B1_D0063.tif" /></maths>

0109If any value of the <i>cor_map_LT</i>(<i>j</i>)<i>,j=0,...,N_SPEC</i>-1, exceeds a threshold of 0.95, a flag <i>cor_strong</i> (which can be viewed as a detector) is set to one, otherwise it is set to zero.

0110The decision about tonal stability is calculated by subjecting <i>cor_map_sum</i> to an adaptive threshold, <i>thr_tonal.</i> This threshold is initialized to 56 and is updated in every frame as follows: <pre listing-type="program-listing"> if (cor_map_sum > 56) thr_tonal = thr_tonal - 0.2 else thr_tonal = thr_tonal + 0.2 end.</pre>

0111The adaptive threshold <i>thr_tonal</i> is upper limited by 60 and lower limited by 49. Thus, the adaptive threshold <i>thr_tonal</i> decreases when the correlation is relatively good indicating an active signal segment and increases otherwise. When the threshold is lower, more frames are likely to be classified as active, especially at the end of active periods. Therefore, the adaptive threshold may be viewed as a hangover.

0112The <i>tonal_stability</i> parameter is set to one whenever <i>cor_map_sum</i> is higher than <i>thr_tonal</i> or when <i>cor_strong</i> flag is set to one. More specifically: <pre listing-type="program-listing"> if ((cor_map_sum > thr_tonal) OR (cor_strong = 1)) tonal_stability = 1 else tonal_stability = 0 end.</pre>

Use of the music detection parameters in noise energy update

0113All music detection parameters are incorporated in the final decision made in the parametric sound activity detection and noise estimation update (Up) module 107 about update of the noise energy estimates. The noise energy estimates are updated as long as the value of <i>noise</i>_<i>update</i> is zero. Initially, it is set to 6 and updated in each frame as follows: <pre listing-type="program-listing"> if (nonstat > thstat) OR (pc < 14) OR (voicing > thCnorm) OR (resid_ratio > thresid) OR (tonal_stability = 1) OR (noise_char_LT > 0.3) OR ((act_pred_LT > 0.8) AND (nonstat2 > thstat)) noise_update = noise_update + 2 else noise_update = noise_update - 1 end.</pre>

0114If the combined condition has a positive result, the signal is active and the <i>noise</i>_<i>update</i> parameter is increased. Otherwise, the signal is inactive and the parameter is decreased. When it reaches 0, the noise energy is updated with the current signal energy.

0115In addition to the noise energy update, the <i>tonal_stability</i> parameter is also used in the classification algorithm of unvoiced sound signal. Specifically, the parameter is used to improve the robustness of unvoiced signal classification on music as will be described in the following section.

Sound <b>signal</b> classification (Sound <b>signal classifier 108)</b>

0116The general philosophy under the sound signal classifier 108 (<figref idref="f0001">Figure 1</figref>) is depicted in <figref idref="f0005">Figure 5</figref>. The approach can be described as follows. The sound signal classification is done in three steps in logic modules 501, 502, and 503, each of them discriminating a specific signal class. First, a signal activity detector (SAD) 501 discriminates between active and inactive signal frames. This signal activity detector 501 is the same as that referred to as signal activity detector 103 in <figref idref="f0001">Figure 1</figref>. The signal activity detector has already been described in the foregoing description.

0117If the signal activity detector 501 detects an inactive frame (background noise signal), then the classification chain ends and, if Discontinuous Transmission (DTX) is supported, an encoding module 541 that can be incorporated in the encoder 109 (<figref idref="f0001">Figure 1</figref>) encodes the frame with comfort noise generation (CNG). If DTX is not supported, the frame continues into the active signal classification, and is most often classified as unvoiced speech frame.

0118If an active signal frame is detected by the sound activity detector 501, the frame is subjected to a second classifier 502 dedicated to discriminate unvoiced speech frames. If the classifier 502 classifies the frame as unvoiced speech signal, the classification chain ends, an encoding module 542 that can be incorporated in the encoder 109 (<figref idref="f0001">Figure 1</figref>) encodes the frame with an encoding method optimized for unvoiced speech signals.

0119Otherwise, the signal frame is processed through to a "stable voiced" classifier 503. If the frame is classified as a stable voiced frame by the classifier 503, then an encoding module 543 that can be incorporated in the encoder 109 (<figref idref="f0001">Figure 1</figref>) encodes the frame using a coding method optimized for stable voiced or quasi periodic signals.

0120Otherwise, the frame is likely to contain a non-stationary signal segment such as a voiced speech onset or rapidly evolving voiced speech or music signal. These frames typically require a general purpose encoding module 544 that can be incorporated in the encoder 109 (<figref idref="f0001">Figure 1</figref>) to encode the frame at high bit rate for sustaining good subjective quality.

0121In the following, the classification of unvoiced and voiced signal frames will be disclosed. The SAD detector 501 (or 103 in <figref idref="f0001">Figure 1</figref>) used to discriminate inactive frames has been already described in the foregoing description.

0122The unvoiced parts of the speech signal are characterized by missing the periodic component and can be further divided into unstable frames, where the energy and the spectrum changes rapidly, and stable frames where these characteristics remain relatively stable. The non-restrictive illustrative embodiment of the present invention proposes a method for the classification of unvoiced frames using the following parameters: <ul id="ul0006" list-style="bullet" compact="compact"><li>voicing measure, computed as an averaged normalized correlation (<i><o ostyle="single">r</o>_x</i>);</li><li>average spectral tilt measure (<i><o ostyle="single">e</o></i>);</li><li>maximum short-time energy increase from low level (<i>dE</i>0) designed to efficiently detect speech plosives in a signal;</li><li>tonal stability to discriminate music from unvoiced signal (described in the foregoing description); and</li><li>relative frame energy (<i>E</i>_rel) to detect very low-energy signals.</li></ul>

<i>Voicing measure</i>

0123The normalized correlation, used to determine the voicing measure, is computed as part of the open-loop pitch analysis made in the LP analyzer and pitch tracker module 106 of <figref idref="f0001">Figure 1</figref>. Frames of 20 ms, for example, can be used. The LP analyzer and pitch tracker module 106 usually outputs an open-loop pitch estimate every 10 ms (twice per frame). Here, the LP analyzer and pitch tracker module 106 is also used to produce and output the normalized correlation measures. These normalized correlations are computed on a weighted signal and a past weighted signal at the open-loop pitch delay. The weighted speech signal <i>s_w</i>(<i>n</i>) is computed using a perceptual weighting filter. For example, a perceptual weighting filter with fixed denominator, suited for wideband signals, can be used. An example of a transfer function for the perceptual weighting filter is given by the following relation: <maths id="math0064" num=""><math display="block"><mi>W</mi><mfenced><mi>z</mi></mfenced><mo>=</mo><mfrac><mrow><mi>A</mi><mfenced separators=""><mi mathvariant="italic">z</mi><mo>/</mo><msub><mi mathvariant="italic">γ</mi><mn mathvariant="italic">1</mn></msub></mfenced></mrow><mrow><mn>1</mn><mo>-</mo><msub><mi mathvariant="italic">γ</mi><mn mathvariant="italic">2</mn></msub><mo>⁢</mo><msup><mi mathvariant="italic">z</mi><mrow><mo>-</mo><mn mathvariant="italic">1</mn></mrow></msup></mrow></mfrac><mo>,</mo></math><img file="EP2162880B1_D0064.tif" /></maths> where 0 &lt; γ₂ &lt; γ₁ ≤ 1 where <i>A(z)</i> is the transfer function of a linear prediction (LP) filter computed in the LP analyzer and pitch tracker module 106, which is given by the following relation: <maths id="math0065" num=""><math display="block"><mi mathvariant="italic">A</mi><mfenced><mi mathvariant="italic">z</mi></mfenced><mo>=</mo><mn>1</mn><mo>+</mo><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>p</mi></munderover></mstyle><msub><mi>a</mi><mi>i</mi></msub><mo>⁢</mo><msup><mi>z</mi><mrow><mo>-</mo><mi>i</mi></mrow></msup><mn>.</mn></math><img file="EP2162880B1_D0065.tif" /></maths>

0124The details of the LP analysis and open-loop pitch analysis will not be further described in the present specification since they are believed to be well known to those of ordinary skill in the art.

0125The voicing measure is given by the average correlation <i><o ostyle="single">C</o>_norm</i> which is defined as: <maths id="math0066" num="(36)"><math display="block"><msub><mover><mi>C</mi><mo>‾</mo></mover><mi mathvariant="italic">norm</mi></msub><mo>=</mo><mfrac><mn>1</mn><mn>3</mn></mfrac><mo>⁢</mo><mfenced separators=""><msub><mi>C</mi><mi mathvariant="italic">norm</mi></msub><mfenced><msub><mi>d</mi><mn>0</mn></msub></mfenced><mo>+</mo><msub><mi>C</mi><mi mathvariant="italic">norm</mi></msub><mfenced><msub><mi>d</mi><mn>1</mn></msub></mfenced><mo>+</mo><msub><mi>C</mi><mi mathvariant="italic">norm</mi></msub><mfenced><msub><mi>d</mi><mn>2</mn></msub></mfenced></mfenced><mo>+</mo><msub><mi>r</mi><mi>e</mi></msub></math><img file="EP2162880B1_D0066.tif" /></maths> where <i>C_norm</i>(<i>d</i>₀), <i>C_norm</i>(<i>d</i>₁) and <i>C_norm</i>(<i>d</i>₂) are respectively the normalized correlation of the first half of the current frame, the normalized correlation of the second half of the current frame, and the normalized correlation of the lookahead (the beginning of the next frame). The arguments to the correlations are the above mentioned open-loop pitch lags calculated in the LP analyzer and pitch tracker module 106 of <figref idref="f0001">Figure 1</figref>. A lookahead of 10 ms can be used, for example. A correction factor <i>r_e</i> is added to the average correlation in order to compensate for the background noise (in the presence of background noise the correlation value decreases). The correction factor is calculated using the following relation: <maths id="math0067" num="(37)"><math display="block"><msub><mi>r</mi><mi>e</mi></msub><mo>=</mo><mn>0.00024492</mn><mo>⁢</mo><msup><mi>e</mi><mrow><mn>0.1596</mn><mo>⁢</mo><mfenced separators=""><msub><mi>N</mi><mi mathvariant="italic">tot</mi></msub><mo>-</mo><mn>14</mn></mfenced></mrow></msup><mo>-</mo><mn>0.022</mn></math><img file="EP2162880B1_D0067.tif" /></maths> where <i>N</i>_tot is the total noise energy per frame computed according to Equation (11).

<i>Spectral tilt</i>

0126The spectral tilt parameter contains information about frequency distribution of energy. The spectral tilt can be estimated in the frequency domain as a ratio between the energy concentrated in low frequencies and the energy concentrated in high frequencies. However, it can be also estimated using other methods such as a ratio between the two first autocorrelation coefficients of the signal.

0127The spectral analyzer 102 in <figref idref="f0001">Figure 1</figref> is used to perform two spectral analyses per frame as described in the foregoing description. The energy in high frequencies and in low frequencies is computed following the perceptual critical bands [<nplcit id="ncit0007" npl-type="s"><text>M. Jelinek and R. Salami, "Noise Reduction Method for Wideband Speech Coding," in Proc. Eusipco, Vienna, Austria, September 2004</text></nplcit>], repeated here for convenience <maths id="math0068" num=""><math display="block"><mtable columnalign="left"><mtr><mtd><mi>Critical bands</mi><mo>=</mo><mrow><mo>{</mo><mn>100.0</mn><mo>,</mo><mn>200.0</mn><mo>,</mo><mn>300.0</mn><mo>,</mo><mn>400.0</mn><mo>,</mo><mn>510.0</mn><mo>,</mo><mn>630.0</mn><mo>,</mo><mn>770.0</mn><mo>,</mo><mn>920.0</mn><mo>,</mo></mrow></mtd></mtr><mtr><mtd><mn>1080.0</mn><mo>,</mo><mn>1270.0</mn><mo>,</mo><mn>1480.0</mn><mo>,</mo><mn>1720.0</mn><mo>,</mo><mn>2000.0</mn><mo>,</mo><mn>2320.0</mn><mo>,</mo><mn>2700.0</mn><mo>,</mo><mn>3150.0</mn><mo>,</mo><mn>3700.0</mn><mo>,</mo></mtd></mtr><mtr><mtd><mn>4400.0</mn><mo>,</mo><mn>5300.0</mn><mo>,</mo><mn>6350.0</mn><mrow><mo>}</mo></mrow><mspace width="1em" /><mi>Hz</mi><mn>.</mn></mtd></mtr></mtable></math><img file="EP2162880B1_D0068.tif" /></maths>

0128The energy in high frequencies is computed as the average of the energies of the last two critical bands using the following relations: <maths id="math0069" num="(39)"><math display="block"><msub><mover><mi>E</mi><mo>‾</mo></mover><mi>h</mi></msub><mo>=</mo><mn>0.5</mn><mo>⁢</mo><mfenced open="[" close="]" separators=""><msub><mi>E</mi><mi mathvariant="italic">CB</mi></msub><mo>⁢</mo><mfenced separators=""><msub><mi>b</mi><mi>max</mi></msub><mo>-</mo><mn>1</mn></mfenced><mo>+</mo><msub><mi>E</mi><mi mathvariant="italic">CB</mi></msub><mfenced><msub><mi>b</mi><mi>max</mi></msub></mfenced></mfenced></math><img file="EP2162880B1_D0069.tif" /></maths> where the critical band energies <i>E_CB(i)</i> are calculated according to Equation (2). The computation is performed twice for both spectral analyses.

0129The energy in low frequencies is computed as the average of the energies in the first 10 critical bands (for NB signals, the very first band is not included), using the following relation: <maths id="math0070" num="(40)"><math display="block"><msub><mover><mi>E</mi><mo>‾</mo></mover><mi>l</mi></msub><mo>=</mo><mfrac><mn>1</mn><mrow><mn>10</mn><mo>-</mo><msub><mi>b</mi><mi>min</mi></msub></mrow></mfrac><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><msub><mi>b</mi><mi>min</mi></msub></mrow><mn>9</mn></munderover></mstyle><msub><mi>E</mi><mi mathvariant="italic">CB</mi></msub><mfenced><mi>i</mi></mfenced><mn>.</mn></math><img file="EP2162880B1_D0070.tif" /></maths>

0130The middle critical bands have been excluded from the computation to improve the discrimination between frames with high energy concentration in low frequencies (generally voiced) and with high energy concentration in high frequencies (generally unvoiced). In between, the energy content is not characteristic for any of the classes and increases the decision confusion.

0131However, the energy in low frequencies is computed differently for harmonic unvoiced signals with high energy content in low frequencies. This is due to the fact that for voiced female speech segments, the harmonic structure of the spectrum can be exploited to increase the voiced-unvoiced discrimination. The affected signals are either those whose pitch period is shorter than 128 or those which are not considered as a priori unvoiced. A priori unvoiced sound signals must fulfill the following condition: <maths id="math0071" num="(41)"><math display="block"><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>⁢</mo><mfenced separators=""><msub><mi>C</mi><mi mathvariant="italic">norm</mi></msub><mfenced><msub><mi>d</mi><mn>0</mn></msub></mfenced><mo>+</mo><msub><mi>C</mi><mi mathvariant="italic">norm</mi></msub><mfenced><msub><mi>d</mi><mn>1</mn></msub></mfenced></mfenced><mo>+</mo><msub><mi>r</mi><mi>e</mi></msub><mo>&lt;</mo><mn>0.6.</mn></math><img file="EP2162880B1_D0071.tif" /></maths>

0132Thus, for the signals discriminated by the above condition, the energy in low frequencies is computed bin-wise and only frequency bins sufficiently close to the harmonics are taken into account into the summation. More specifically, the following relation is used: <maths id="math0072" num="(42)"><math display="block"><msub><mover><mi>E</mi><mo>‾</mo></mover><mi>l</mi></msub><mo>=</mo><mfrac><mn>1</mn><mi mathvariant="italic">cnt</mi></mfrac><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><msub><mi>K</mi><mi>min</mi></msub></mrow><mn>25</mn></munderover></mstyle><msub><mi>E</mi><mi mathvariant="italic">BIN</mi></msub><mfenced><mi>i</mi></mfenced><mo>⁢</mo><msub><mi>w</mi><mi>h</mi></msub><mfenced><mi>i</mi></mfenced><mn>.</mn></math><img file="EP2162880B1_D0072.tif" /></maths> where <i>K_min</i> is the first bin (<i>K_min</i>=1 for WB and <i>K_min</i>=3 for NB) and <i>E_BIN</i>(<i>k</i>) are the bin energies, as defined in Equation (3), in the first 25 frequency bins (the DC component is omitted). These 25 bins correspond to the first 10 critical bands. In the summation above, only terms close to the pitch harmonics are considered; <i>w_h</i>(<i>i</i>) is set to 1 if the distance between the nearest harmonics is not larger than a certain frequency threshold (for example 50 Hz) and is set to 0 otherwise; therefore only bins closer than 50 Hz to the nearest harmonics are taken into account. The counter <i>cnt</i> is equal to the number of non-zero terms in the summation. Hence, if the structure is harmonic in low frequencies, only high energy terms will be included in the sum. On the other hand, if the structure is not harmonic, the selection of the terms will be random and the sum will be smaller. Thus even unvoiced sound signals with high energy content in low frequencies can be detected.

0133The spectral tilt is given by the following relation: <maths id="math0073" num="(43)"><math display="block"><msub><mi>e</mi><mi>t</mi></msub><mo>=</mo><mfrac><mrow><msub><mover><mi>E</mi><mo>‾</mo></mover><mi>l</mi></msub><mo>-</mo><msub><mover><mi>N</mi><mo>‾</mo></mover><mi>l</mi></msub></mrow><mrow><msub><mover><mi>E</mi><mo>‾</mo></mover><mi>h</mi></msub><mo>-</mo><msub><mover><mi>N</mi><mo>‾</mo></mover><mi>h</mi></msub></mrow></mfrac></math><img file="EP2162880B1_D0073.tif" /></maths> where <i><o ostyle="single">N</o>_h</i> and <i><o ostyle="single">N</o>_l</i> are the averaged noise energies in the last two (2) critical bands and the first 10 critical bands (or the first 9 critical bands for NB), respectively, computed in the same way as <i><o ostyle="single">E</o>_h</i> and <i><o ostyle="single">E</o>_l</i> in Equations (39) and (40). The estimated noise energies have been included in the tilt computation to account for the presence of background noise. For NB signals, the missing bands are compensated by multiplying <i>e_t</i> by 6. The spectral tilt computation is performed twice per frame to obtain <i>e_t</i>(0) and <i>e_t</i>(1) corresponding to both the first and second spectral analyses per frame. The average spectral tilt used in unvoiced frame classification is given by <maths id="math0074" num="(44)"><math display="block"><msub><mover><mi>e</mi><mo>‾</mo></mover><mi>t</mi></msub><mo>=</mo><mfrac><mn>1</mn><mn>3</mn></mfrac><mo>⁢</mo><mfenced separators=""><msub><mi>e</mi><mi mathvariant="italic">old</mi></msub><mo>+</mo><msub><mi>e</mi><mi>t</mi></msub><mfenced><mn>0</mn></mfenced><mo>+</mo><msub><mi>e</mi><mi>t</mi></msub><mfenced><mn>1</mn></mfenced></mfenced><mo>,</mo></math><img file="EP2162880B1_D0074.tif" /></maths> where <i>e</i>_old is the tilt in the second half of the previous frame.

<i>Maximum short-time energy increase at low level</i>

0134The maximum short-time energy increase at low level <i>dE</i>0 is evaluated on the sound signal <i>s</i>(<i>n</i>), where <i>n</i>=0 corresponds to the beginning of the current frame. For example, 20 ms speech frames are used and every frame is divided into 4 subframes for speech encoding purposes. The signal energy is evaluated twice per subframe, i.e. 8 times per frame, based on short-time segments of a length of 32 samples (at a 12.8 kHz sampling rate). Further, the short-term energies of the last 32 samples from the previous frame are also computed. The short-time energies are computed using the following relation: <maths id="math0075" num="(45)"><math display="block"><msubsup><mi>E</mi><mi mathvariant="italic">st</mi><mfenced><mn>1</mn></mfenced></msubsup><mfenced><mi>j</mi></mfenced><mo>=</mo><munderover><mi>max</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mn>31</mn></munderover><mfenced separators=""><msup><mi>s</mi><mn>2</mn></msup><mo>⁢</mo><mfenced separators=""><mi>i</mi><mo>+</mo><mn>32</mn><mo>⁢</mo><mi>j</mi></mfenced></mfenced><mo>,</mo><mspace width="1em" /><mi>j</mi><mo>=</mo><mo>-</mo><mn>1</mn><mo>,</mo><mo>…</mo><mo>,</mo><mn>7</mn><mo>,</mo></math><img file="EP2162880B1_D0075.tif" /></maths> where <i>j</i>=-1 and <i>j</i>=0,...,7 correspond to the end of the previous frame and the current frame, respectively. Another set of 9 maximum energies is computed by shifting the signal indices in Equation (45) by 16 samples. That is <maths id="math0076" num="(46)"><math display="block"><msubsup><mi>E</mi><mi mathvariant="italic">st</mi><mfenced><mn>2</mn></mfenced></msubsup><mfenced><mi>j</mi></mfenced><mo>=</mo><munderover><mi>max</mi><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mn>31</mn></munderover><mfenced separators=""><msup><mi>s</mi><mn>2</mn></msup><mo>⁢</mo><mfenced separators=""><mi>i</mi><mo>+</mo><mn>32</mn><mo>⁢</mo><mi>j</mi><mo>+</mo><mn>16</mn></mfenced></mfenced><mo>,</mo><mspace width="1em" /><mi>j</mi><mo>=</mo><mo>-</mo><mn>1</mn><mo>,</mo><mo>…</mo><mo>,</mo><mn>7.</mn></math><img file="EP2162880B1_D0076.tif" /></maths>

0135For those energies that are sufficiently low, i.e. which fulfill the condition 10log(<i>E_st</i>(<i>j</i>)) &lt; 37, the following ratio is calculated: <maths id="math0077" num="(47)"><math display="block"><msup><mi mathvariant="italic">rat</mi><mfenced><mn>1</mn></mfenced></msup><mfenced><mi>j</mi></mfenced><mo>=</mo><mfrac><mrow><msubsup><mi>E</mi><mi mathvariant="italic">st</mi><mfenced><mn>1</mn></mfenced></msubsup><mo>⁢</mo><mfenced separators=""><mi>j</mi><mo>+</mo><mn>1</mn></mfenced></mrow><mrow><msubsup><mi>E</mi><mi mathvariant="italic">st</mi><mfenced><mn>1</mn></mfenced></msubsup><mfenced><mi>j</mi></mfenced><mo>+</mo><mn>100</mn></mrow></mfrac><mo>,</mo><mspace width="1em" /><mi>for</mi><mspace width="1em" /><mi>j</mi><mo>=</mo><mo>-</mo><mn>1</mn><mo>,</mo><mo>…</mo><mo>,</mo><mn>6</mn><mo>,</mo></math><img file="EP2162880B1_D0077.tif" /></maths> for the first set of indices and the same calculation is repeated for <maths id="math0078" num=""><math display="inline"><msubsup><mi>E</mi><mi mathvariant="italic">st</mi><mfenced><mn>2</mn></mfenced></msubsup><mfenced><mi>j</mi></mfenced></math><img file="EP2162880B1_D0078.tif" /></maths>to obtain two sets of ratios <i>rat</i>⁽¹⁾(<i>j</i>) and <i>rat</i>⁽²⁾(<i>j</i>). The only maximum in these two sets is searched as follows: <maths id="math0079" num="(48)"><math display="block"><mi mathvariant="italic">dE</mi><mo>⁢</mo><mn>0</mn><mo>=</mo><mi>max</mi><mfenced separators=""><msup><mi mathvariant="italic">rat</mi><mfenced><mn>1</mn></mfenced></msup><mfenced><mi>j</mi></mfenced><mo>,</mo><msup><mi mathvariant="italic">rat</mi><mfenced><mn>2</mn></mfenced></msup><mfenced><mi>j</mi></mfenced></mfenced></math><img file="EP2162880B1_D0079.tif" /></maths> which is the maximum short-time energy increase at low level.

<i>Measure on background noise spectrum flatness</i>

0136In this example, inactive frames are usually coded with a coding mode designed for unvoiced speech in the absence of DTX operation. However, in the case of a quasi-periodic background noise, like some car noises, more faithful noise rendering is achieved if generic coding is instead used for WB.

0137To detect this type of background noise, a measure of background noise spectrum flatness is computed and averaged over time. First, average noise energy is computed for first and last four critical bands as follows: <maths id="math0080" num=""><math display="block"><msub><mover><mi>N</mi><mo>‾</mo></mover><mrow><mi>l</mi><mo>⁢</mo><mn>4</mn></mrow></msub><mo>=</mo><mfrac><mn>1</mn><mn>4</mn></mfrac><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mn>0</mn></mrow><mn>3</mn></munderover></mstyle><msub><mi>N</mi><mi mathvariant="italic">CB</mi></msub><mfenced><mi>i</mi></mfenced></math><img file="EP2162880B1_D0080.tif" /></maths><maths id="math0081" num=""><math display="block"><msub><mover><mi>N</mi><mo>‾</mo></mover><mrow><mi>h</mi><mo>⁢</mo><mn>4</mn></mrow></msub><mo>=</mo><mfrac><mn>1</mn><mn>4</mn></mfrac><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mn>15</mn></mrow><mn>19</mn></munderover></mstyle><msub><mi>N</mi><mi mathvariant="italic">CB</mi></msub><mfenced><mi>i</mi></mfenced></math><img file="EP2162880B1_D0081.tif" /></maths>

0138The flatness measure is then computed using the following relation: <maths id="math0082" num=""><math display="block"><msub><mi>f</mi><mi mathvariant="italic">noise_flat</mi></msub><mo>=</mo><mfenced separators=""><msub><mover><mi>N</mi><mo>‾</mo></mover><mrow><mi>l</mi><mo>⁢</mo><mn>4</mn></mrow></msub><mo>-</mo><msub><mover><mi>N</mi><mo>‾</mo></mover><mrow><mi>h</mi><mo>⁢</mo><mn>4</mn></mrow></msub></mfenced><mo>/</mo><msub><mover><mi>N</mi><mo>‾</mo></mover><mrow><mi>l</mi><mo>⁢</mo><mn>4</mn></mrow></msub><mo>+</mo><mn>0.5</mn><mo>⁢</mo><mfenced open="[" close="]" separators=""><msub><mi>N</mi><mi mathvariant="italic">CB</mi></msub><mfenced><mn>1</mn></mfenced><mo>+</mo><msub><mi>N</mi><mi mathvariant="italic">CB</mi></msub><mfenced><mn>2</mn></mfenced></mfenced><mo>/</mo><msub><mi>N</mi><mi mathvariant="italic">CB</mi></msub><mfenced><mn>0</mn></mfenced></math><img file="EP2162880B1_D0082.tif" /></maths> and averaged over time using the following relation: <maths id="math0083" num=""><math display="block"><msubsup><mover><mi>f</mi><mo>‾</mo></mover><mi mathvariant="italic">noise_flat</mi><mfenced open="[" close="]"><mn>0</mn></mfenced></msubsup><mo>=</mo><mn>0.99</mn><mo>⁢</mo><msubsup><mover><mi>f</mi><mo>‾</mo></mover><mi mathvariant="italic">noise_flat</mi><mfenced open="[" close="]" separators=""><mo>-</mo><mn>1</mn></mfenced></msubsup><mo>+</mo><mn>0.01</mn><mo>⁢</mo><msub><mi>f</mi><mi mathvariant="italic">noise_flat</mi></msub></math><img file="EP2162880B1_D0083.tif" /></maths> where <maths id="math0084" num=""><math display="inline"><msubsup><mover><mi>f</mi><mo>‾</mo></mover><mi mathvariant="italic">noise_flat</mi><mfenced open="[" close="]" separators=""><mo>-</mo><mn>1</mn></mfenced></msubsup></math><img file="EP2162880B1_D0084.tif" /></maths>is the averaged flatness measure of the past frame and <maths id="math0085" num=""><math display="inline"><msubsup><mover><mi>f</mi><mo>‾</mo></mover><mi mathvariant="italic">noise_flat</mi><mfenced open="[" close="]"><mn>0</mn></mfenced></msubsup></math><img file="EP2162880B1_D0085.tif" /></maths> is the updated value of the averaged flatness measure of the current frame.

<i>Unvoiced signal classification</i>

0139The classification of unvoiced signal frames is based on the parameters described above, namely: the voicing measure <i><o ostyle="single">C</o>_norm</i>, the average spectral tilt <i><o ostyle="single">e</o>_t</i>, the maximum short-time energy increase at low level <i>dE</i>0 and the measure of background noise spectrum flatness, <maths id="math0086" num=""><math display="inline"><msubsup><mover><mi>f</mi><mo>‾</mo></mover><mrow><mi>noise</mi><mi>_</mi><mi>flat</mi></mrow><mn>0</mn></msubsup><mn>.</mn></math><img file="EP2162880B1_D0086.tif" /></maths> The classification is further supported by the tonal stability parameter and the relative frame energy calculated during the noise energy update phase (module 107 in <figref idref="f0001">Figure 1</figref>). The relative frame energy is calculated using the following relation: <maths id="math0087" num="(50)"><math display="block"><msub><mi>E</mi><mi mathvariant="italic">rel</mi></msub><mo>=</mo><msub><mi>E</mi><mi>t</mi></msub><mo>-</mo><msub><mover><mi>E</mi><mo>‾</mo></mover><mi>f</mi></msub></math><img file="EP2162880B1_D0087.tif" /></maths> where <i>E_t</i> is the total frame energy (in dB) calculated in Equation (6) and <i><o ostyle="single">E</o>_f</i> is the long-term average frame energy, updated in each active frame using the following relation: <maths id="math0088" num=""><math display="block"><msub><mover><mi>E</mi><mo>‾</mo></mover><mi>f</mi></msub><mo>=</mo><mn>0.99</mn><mo>⁢</mo><msub><mover><mi>E</mi><mo>‾</mo></mover><mi>f</mi></msub><mo>-</mo><mn>0.01</mn><mo>⁢</mo><msub><mi>E</mi><mi>t</mi></msub><mn>.</mn></math><img file="EP2162880B1_D0088.tif" /></maths>

0140The updating takes place only when SAD flag is set (variable <i>SAD</i> equal to 1).

0141The rules for unvoiced classification of WB signals are summarized below: <ul id="ul0007" list-style="none" compact="compact"><li>[((<i><o ostyle="single">C</o>_norm</i> < 0.695) AND (<i><o ostyle="single">e</o>_t</i> < 4.0)) OR (<i>E</i>_rel < -14)] AND</li><li>[last frame INACTIVE OR UNVOICED OR ((<i>e</i>_old < 2.4) AND</li><li>(<i>C_norm</i>(<i>d</i>₀)+<i>r_e</i> < 0.66))] AND</li><li>[<i>dE</i>0 < 250] AND</li><li>[<i>e_t</i>(1) < 2.7] AND</li><li>[(local SAD flag = 1) OR ( <maths id="math0089" num=""><math display="inline"><msubsup><mover><mi>f</mi><mo>‾</mo></mover><mi mathvariant="italic">noise_flat</mi><mfenced open="[" close="]"><mn>0</mn></mfenced></msubsup></math><img file="EP2162880B1_D0089.tif" /></maths> < 1.45) OR (<i><o ostyle="single">N</o>_f</i> < 20)] AND</li><li>NOT [(<i>tonal_stability</i> AND (((<i><o ostyle="single">C</o>_norm</i> > 0.52) AND (<i><o ostyle="single">e</o>_t</i> > 0.5)) OR (<i><o ostyle="single">e</o>_t</i> > 0.85)) AND (<i>E</i>_rel > -14) AND SAD flag set to 1]</li></ul>

0142The first line of the condition is related to low-energy signals and signals with low correlation concentrating their energy in high frequencies. The second line covers voiced offsets, the third line covers explosive segments of a signal and the fourth line is for the voiced onsets. The fifth line ensures flat spectrum in case of noisy inactive frames. The last line discriminates music signals that would be otherwise declared as unvoiced.

0143For NB signals the unvoiced classification condition takes the following form: <ul id="ul0008" list-style="none" compact="compact"><li>[local SAD flag set to 0 OR (<i>E</i>_rel < -25) OR ((<i><o ostyle="single">C</o>_norm</i> < 0.61) AND (<i><o ostyle="single">e</o>_t</i> < 7.0) AND (last frame INACTIVE OR UNVOICED OR ((<i>e</i>_old < 7.0) AND <i>(C_norm</i>(<i>d</i>₀)+<i>r_e</i> < 0.52))))] AND</li><li>[<i>dE</i>0 < 250] AND</li><li>[<i><o ostyle="single">e</o>_t</i> < 390] AND</li><li>NOT [(<i>tonal_stability</i> AND (((<i><o ostyle="single">C</o>_norm</i> > 0.52) AND (<i><o ostyle="single">e</o>_t</i> > 0.5)) OR (<i><o ostyle="single">e</o>_t</i> > 0.75)) AND (<i>E</i>_rel > -10) AND SAD flag set to 1]</li></ul>

0144The decision trees for the WB case and NB case are shown in <figref idref="f0006">Figure 6</figref>. If the combined conditions are fulfilled the classification ends by selecting unvoiced coding mode.

<i>Voiced signal classification</i>

0145If a frame is not classified as inactive frame or as unvoiced frame then it is tested if it is a stable voiced frame. The decision rule is based on the normalized correlation in each subframe (with 1/4 subsample resolution), the average spectral tilt and open-loop pitch estimates in all subframes (with 1/4 subsample resolution).

0146The open-loop pitch estimation procedure is made by the LP analyzer and pitch tracker module 106 of <figref idref="f0001">Figure 1</figref>. In Equation (19), three open-loop pitch estimates are used: <i>d</i>₀<i>, d</i>₁ and <i>d</i>₂, corresponding to the first half-frame, the second half-frame and the look ahead. In order to obtain precise pitch information in all four subframes, 1/4 sample resolution fractional pitch refinement is calculated. This refinement is calculated on the weighted sound signal <i>s_wd</i>(<i>n</i>)<i>.</i> In this exemplary embodiment, the weighted signal <i>s_wd</i>(<i>n</i>) is not decimated for open-loop pitch estimation refinement. At the beginning of each subframe a short correlation analysis (64 samples at 12.8 kHz sampling frequency) with resolution of 1 sample is done in the interval (-7,+7) using the following delays: do for the first and second subframes and <i>d</i>₁ for the third and fourth subframes. The correlations are then interpolated around their maxima at the fractional positions <i>d</i>_max - 3/4, <i>d</i>_max - 1/2, <i>d</i>_max - 1/4, <i>d</i>_max, <i>d</i>_max + 1/4, <i>d</i>_max + 1/2, <i>d</i>_max + 3/4. The value yielding the maximum correlation is chosen as the refined pitch lag.

0147Let the refined open-loop pitch lags in all four subframes be denoted as <i>T</i>(0), <i>T</i>(1), <i>T</i>(2) and <i>T</i>(3) and their corresponding normalized correlations as <i>C</i>(0), <i>C</i>(1), <i>C(2)</i> and C(3). Then, the voiced signal classification condition is given by: <ul id="ul0009" list-style="none" compact="compact"><li>[<i>C</i>(0) > 0.605] AND</li><li>[<i>C</i>(1) > 0.605] AND</li><li>[C(2) > 0.605] AND</li><li>[C(3) > 0.605] AND</li><li>[<i><o ostyle="single">e</o>_t</i> > 4] AND</li><li>[|<i>T</i>(1) - <i>T</i>(0)| < 3] AND</li><li>[|<i>T</i>(2) - <i>T</i>(1)| < 3] AND</li><li>[|<i>T</i>(3) - <i>T</i>(2)| < 3]</li></ul>

0148The condition says that the normalized correlation is sufficiently high in all subframes, the pitch estimates do not diverge throughout the frame and the energy is concentrated in low frequencies. If this condition is fulfilled the classification ends by selecting voiced signal coding mode, otherwise the signal is encoded by a generic signal coding mode. The condition applies to both WB and NB signals.

<i>Estimation of tonality in the super wideband content</i>

0149In the encoding of super wideband signals, a specific coding mode is used for sound signals with tonal structure. The frequency range which is of interest is mostly 7000 - 14000 Hz but can also be different. The objective is to detect frames having strong tonal content in the range of interest so that the tonal-specific coding mode may be used efficiently. This is done using the tonal stability analysis described earlier in the present disclosure. However, there are some aberrations which are described in this section.

0150First, the spectral floor which is subtracted from the log-energy spectrum is calculated in the following way. The log-energy spectrum is filtered using a moving-average (MA) filter, or FIR filter, the length of which is <i>L_MA</i>=15 samples. The filtered spectrum is given by: <maths id="math0090" num=""><math display="block"><mi mathvariant="italic">sp_floor</mi><mfenced><mi>j</mi></mfenced><mo>=</mo><mfrac><mn>1</mn><mrow><mn>2</mn><mo>⁢</mo><msub><mi>L</mi><mi mathvariant="italic">MA</mi></msub><mo>+</mo><mn>1</mn></mrow></mfrac><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>k</mi><mo>=</mo><mo>-</mo><msub><mi>L</mi><mi mathvariant="italic">MA</mi></msub></mrow><msub><mi>L</mi><mi mathvariant="italic">MA</mi></msub></munderover></mstyle><msub><mi>E</mi><mi mathvariant="italic">dB</mi></msub><mo>⁢</mo><mfenced separators=""><mi>j</mi><mo>+</mo><mi>k</mi></mfenced><mo>,</mo><mspace width="1em" /><mi>for</mi><mspace width="1em" /><mi>j</mi><mo>=</mo><msub><mi>L</mi><mi mathvariant="italic">MA</mi></msub><mo>,</mo><mo>…</mo><mo>,</mo><msub><mi>N</mi><mi mathvariant="italic">SPEC</mi></msub><mo>-</mo><msub><mi>L</mi><mi mathvariant="italic">MA</mi></msub><mo>-</mo><mn>1.</mn></math><img file="EP2162880B1_D0090.tif" /></maths>

0151To save computational complexity, the filtering operation is done only for <i>j</i>=<i>L_MA</i> and for the other lags, it is calculated as: <maths id="math0091" num=""><math display="block"><mtable columnalign="left"><mtr><mtd><mi mathvariant="italic">sp_floor</mi><mfenced><mi>j</mi></mfenced><mo>=</mo><mi mathvariant="italic">sp_floor</mi><mo>⁢</mo><mfenced separators=""><mi>j</mi><mo>-</mo><mn>1</mn></mfenced><mo>+</mo><mfrac><mn>1</mn><mrow><mn>2</mn><mo>⁢</mo><msub><mi>L</mi><mi mathvariant="italic">MA</mi></msub><mo>+</mo><mn>1</mn></mrow></mfrac><mo>⁢</mo><mfenced open="[" close="]" separators=""><msub><mi>E</mi><mi mathvariant="italic">dB</mi></msub><mo>⁢</mo><mfenced separators=""><mi>j</mi><mo>+</mo><msub><mi>L</mi><mi mathvariant="italic">MA</mi></msub></mfenced><mo>-</mo><msub><mi>E</mi><mi mathvariant="italic">dB</mi></msub><mo>⁢</mo><mfenced separators=""><mi>j</mi><mo>-</mo><msub><mi>L</mi><mi mathvariant="italic">MA</mi></msub><mo>-</mo><mn>1</mn></mfenced></mfenced><mo>,</mo></mtd></mtr><mtr><mtd><mi>for</mi><mspace width="1em" /><mi>j</mi><mo>=</mo><msub><mi>L</mi><mi mathvariant="italic">MA</mi></msub><mo>+</mo><mn>1</mn><mo>,</mo><mo>…</mo><mo>,</mo><msub><mi>N</mi><mi mathvariant="italic">SPEC</mi></msub><mo>-</mo><msub><mi>L</mi><mi mathvariant="italic">MA</mi></msub><mo>-</mo><mn>1.</mn></mtd></mtr></mtable></math><img file="EP2162880B1_D0091.tif" /></maths>

0152For the lags 0,.., <i>L_MA</i>-1 and <i>N_SPEC-L_MA</i>,...,<i>N_SPEC</i>-1<i>,</i> the spectral floor is calculated by means of extrapolation. More specifically, the following relation is used: <maths id="math0092" num=""><math display="block"><mtable><mtr><mtd><mi mathvariant="italic">sp_floor</mi><mfenced><mi>j</mi></mfenced><mo>=</mo><mn>0.9</mn><mo>⁢</mo><mi mathvariant="italic">sp_floor</mi><mo>⁢</mo><mfenced separators=""><mi>j</mi><mo>+</mo><mn>1</mn></mfenced><mo>+</mo><mn>0.1</mn><mo>⁢</mo><msub><mi>E</mi><mi mathvariant="italic">dB</mi></msub><mfenced><mi>j</mi></mfenced><mo>,</mo><mspace width="1em" /><mi>for</mi><mspace width="1em" /><mi>j</mi><mo>=</mo><msub><mi>L</mi><mi mathvariant="italic">MA</mi></msub><mo>-</mo><mn>1</mn><mo>,</mo><mo>…</mo><mo>,</mo><mn>0</mn><mo>,</mo></mtd></mtr><mtr><mtd><mi mathvariant="italic">sp_floor</mi><mfenced><mi>j</mi></mfenced><mo>=</mo><mn>0.9</mn><mo>⁢</mo><mi mathvariant="italic">sp_floor</mi><mo>⁢</mo><mfenced separators=""><mi>j</mi><mo>-</mo><mn>1</mn></mfenced><mo>+</mo><mn>0.1</mn><mo>⁢</mo><msub><mi>E</mi><mi mathvariant="italic">dB</mi></msub><mfenced><mi>j</mi></mfenced><mo>,</mo><mspace width="1em" /><mi>for</mi><mspace width="1em" /><mi>j</mi><mo>=</mo><msub><mi>N</mi><mi mathvariant="italic">SPEC</mi></msub><mo>-</mo><msub><mi>L</mi><mi mathvariant="italic">MA</mi></msub><mo>,</mo><mo>…</mo><mo>,</mo><msub><mi>N</mi><mi mathvariant="italic">SPEC</mi></msub><mo>-</mo><mn>1.</mn></mtd></mtr></mtable></math><img file="EP2162880B1_D0092.tif" /></maths>

0153In the first equation above the updating proceeds from <i>L_MA</i>-1 downwards to 0.

0154The spectral floor is then subtracted from the log-energy spectrum in the same way as described earlier in the present disclosure.

0155The residual spectrum, denoted as <i>E_res,dB</i>(<i>j</i>)<i>,</i> is then smoothed over 3 samples as follows using a short-time moving-average filter: <maths id="math0093" num=""><math display="block"><msubsup><mi>E</mi><mrow><mi mathvariant="italic">res</mi><mo>,</mo><mi mathvariant="italic">dB</mi></mrow><mi>ʹ</mi></msubsup><mfenced><mi>j</mi></mfenced><mo>=</mo><mn>0.33</mn><mo>⁢</mo><mfenced open="[" close="]" separators=""><msub><mi>E</mi><mrow><mi mathvariant="italic">res</mi><mo>,</mo><mi mathvariant="italic">dB</mi></mrow></msub><mo>⁢</mo><mfenced separators=""><mi>j</mi><mo>-</mo><mn>1</mn></mfenced><mo>+</mo><msub><mi>E</mi><mrow><mi mathvariant="italic">res</mi><mo>,</mo><mi mathvariant="italic">dB</mi></mrow></msub><mfenced><mi>j</mi></mfenced><mo>+</mo><msub><mi>E</mi><mrow><mi mathvariant="italic">res</mi><mo>,</mo><mi mathvariant="italic">dB</mi></mrow></msub><mo>⁢</mo><mfenced separators=""><mi>j</mi><mo>+</mo><mn>1</mn></mfenced></mfenced><mo>,</mo><mspace width="1em" /><mi>for</mi><mspace width="1em" /><mi>j</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>…</mo><mo>,</mo><msub><mi>N</mi><mi mathvariant="italic">SPEC</mi></msub><mo>-</mo><mn>1.</mn></math><img file="EP2162880B1_D0093.tif" /></maths>

0156The search of spectral minima and their indexes, the calculation of correlation map and the long term correlation map are the same as in the method described earlier in the present disclosure, using the smoothed spectrum <maths id="math0094" num=""><math display="inline"><msubsup><mi>E</mi><mrow><mi mathvariant="italic">res</mi><mo>,</mo><mi mathvariant="italic">dB</mi></mrow><mi>ʹ</mi></msubsup><mfenced><mi>j</mi></mfenced><mn>.</mn></math><img file="EP2162880B1_D0094.tif" /></maths>

0157The decision about signal tonality in the super-wideband content is also the same as described earlier in the present disclosure, i.e. based on an adaptive threshold. However, in this case a different fixed threshold and step are used. The threshold <i>thr_tonal</i> is initialized to 130 and is updated in every frame as follows: <pre listing-type="program-listing"> if (cor_map_sum > 130) thr_tonal = thr_tonal - 1.0 else thr_tonal = thr_tonal + 1.0 end.</pre>

0158The adaptive threshold <i>thr_tonal</i> is upper limited by 140 and lower limited by 120. The fixed threshold has been set with respect to the frequency range 7000 - 14000 Hz. For a different range, it will have to be adjusted. As a general rule of thumb, the following relationship may be applied <i>thr_tonal=N_SPEC</i>/2.

0159The last difference to the method described earlier in the present disclosure is that the detection of strong tones is not used in the super wideband content. This is motivated by the fact that strong tones are perceptually not suitable for the purpose of encoding the tonal signal in the super wideband content.

0160The present invention has been described in the foregoing disclosure by way of a non-restrictive, illustrative embodiment thereof. The scope of the present invention is defined by the appended claims.