EP3020043A1 - Optimized scale factor for frequency band extension in an audiofrequency signal decoder - Google Patents

Optimized scale factor for frequency band extension in an audiofrequency signal decoder

Info

Publication number
EP3020043A1
EP3020043A1 EP14749907.3A EP14749907A EP3020043A1 EP 3020043 A1 EP3020043 A1 EP 3020043A1 EP 14749907 A EP14749907 A EP 14749907A EP 3020043 A1 EP3020043 A1 EP 3020043A1
Authority
EP
European Patent Office
Prior art keywords
frequency
filter
band
frequency band
signal
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
EP14749907.3A
Other languages
German (de)
French (fr)
Other versions
EP3020043B1 (en
Inventor
Magdalena KANIEWSKA
Stéphane RAGOT
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Koninklijke Philips NV
Original Assignee
Orange SA
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Family has litigation
First worldwide family litigation filed litigation Critical https://patents.darts-ip.com/?family=49753286&utm_source=***_patent&utm_medium=platform_link&utm_campaign=public_patent_search&patent=EP3020043(A1) "Global patent litigation dataset” by Darts-ip is licensed under a Creative Commons Attribution 4.0 International License.
Application filed by Orange SA filed Critical Orange SA
Publication of EP3020043A1 publication Critical patent/EP3020043A1/en
Application granted granted Critical
Publication of EP3020043B1 publication Critical patent/EP3020043B1/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L19/00—Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis
    • G10L19/04—Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis using predictive techniques
    • G10L19/08—Determination or coding of the excitation function; Determination or coding of the long-term prediction parameters
    • G10L19/087—Determination or coding of the excitation function; Determination or coding of the long-term prediction parameters using mixed excitation models, e.g. MELP, MBE, split band LPC or HVXC
    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L25/00—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
    • G10L25/48—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 specially adapted for particular use
    • G10L25/72—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 specially adapted for particular use for transmitting results of analysis
    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L21/00—Speech or voice signal processing techniques to produce another audible or non-audible signal, e.g. visual or tactile, in order to modify its quality or its intelligibility
    • G10L21/02—Speech enhancement, e.g. noise reduction or echo cancellation
    • G10L21/038—Speech enhancement, e.g. noise reduction or echo cancellation using band spreading techniques
    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L19/00—Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis
    • G10L19/005—Correction of errors induced by the transmission channel, if related to the coding algorithm
    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L19/00—Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis
    • G10L19/008—Multichannel audio signal coding or decoding using interchannel correlation to reduce redundancy, e.g. joint-stereo, intensity-coding or matrixing
    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L19/00—Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis
    • G10L19/02—Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis using spectral analysis, e.g. transform vocoders or subband vocoders
    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L19/00—Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis
    • G10L19/04—Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis using predictive techniques
    • G10L19/16—Vocoder architecture
    • G10L19/18—Vocoders using multiple modes
    • G10L19/24—Variable rate codecs, e.g. for generating different qualities using a scalable representation such as hierarchical encoding or layered encoding

Definitions

  • the present invention relates to the field of coding / decoding and audio-frequency signal processing (such as speech, music or other signals) for their transmission or storage.
  • the invention relates to a method and an apparatus for determining an optimized scale factor for adjusting the level of an excitation signal or, in a similar manner, a filter during a band extension. frequency in a decoder or a processor performing audio-frequency signal enhancement.
  • 3GPP AMR-WB Adaptive Multi-Rate Wideband
  • codec and decoder which operates at an input / output frequency of 16 kHz and in which the signal is divided into two sub-bands, the low band (0-6.4 kHz) which is sampled at 12.8 kHz and coded by CELP model and the high band (6.4-7 kHz) which is parametrically reconstructed by "band extension" ( or BWE for "Bandwidth Extension” with or without additional information depending on the mode of the current frame.
  • BWE Bandwidth Extension
  • the limitation of the coded band of the AMR-WB codec at 7 kHz is essentially related to the fact that the transmission frequency response of the broadband terminals has been approximated at the time of standardization (ETSI / 3GPP then ITU-T T) according to the frequency mask defined in the ITU-T P.341 standard and more precisely by using a so-called "P341" filter defined in the ITU-T G.191 standard. which cuts frequencies above 7 kHz (this filter respects the mask defined in P.341).
  • a signal sampled at 16 kHz may have a defined audio band of 0 to 8000 Hz; the AMR-WB codec thus introduces a limitation of the high band in comparison with the theoretical bandwidth of 8 kHz.
  • the 3GPP AMR-WB speech codec was standardized in 2001 mainly for circuit-mode (CS) telephony applications on GSM (2G) and UMTS (3G). This same codec was also standardized in 2003 in ITU-T as Recommendation G.722.2 "Wideband coding speech at around 16kbit / s using Adaptive Multi-Rate Wideband (AMR-WB)".
  • AMR-WB coding and decoding algorithm The details of the AMR-WB coding and decoding algorithm are not repeated here, a detailed description of this codec is found in the 3GPP specifications (TS 26.190, 26.191, 26.192, 26.193, 26.194, 26.204) and ITU-TG .722.2 (and the corresponding Appendices and Appendix) and in the article by B. Bessette et al. entitled "The adaptive multirate broadband speech coded (AMR-WB)", IEEE Transactions on Speech and Audio Processing, vol. 10, no. 8, 2002, pp. 620-636 and associated 3GPP and ITU-T standard source codes.
  • AMR-WB adaptive multirate broadband speech coded
  • the principle of band extension in the AMR-WB codec is rather rudimentary. Indeed, the high band (6.4-7 kHz) is generated by formatting a white noise through a temporal envelope (applied in the form of gains per subframe) and frequency (by the application of a linear prediction synthesis filter or LPC for "Linear Predictive Coding").
  • This band extension technique is illustrated in Figure 1.
  • a white noise, u HB1 (n), n 0, ⁇ ⁇ ⁇ , 79, is generated at 16 kHz per 5 ms subframe per linear congruent generator (block 100).
  • This noise u HB1 (n) is shaped in time by applying gains per subframe; this operation is broken down into two processing steps (blocks 102, 106 or 109):
  • a first factor is calculated (block 101) to set the white noise u HB1 (n) (block 101).
  • the normalization of the energies is done by comparing blocks of different size (64 for u (n) and 80 for u HB1 (n)), without compensation of the differences of sampling frequencies (12.8 or 16 kHz) .
  • u HB ⁇ n g HB u HB2 ⁇ n
  • w sp is a weighting function that depends on Voice Activity Detection (VAD).
  • VAD Voice Activity Detection
  • the factor g HB in the decoding AMR-WB is bounded to take values in the interval [0.1, 1.0]. Indeed, for signals whose energy increases when the frequency increases (e tilt close to -1, g sp close to 2), the gain g HB is usually underestimated.
  • correction information is transmitted by the encoder AMR-WB and decoded (blocks 107, 108) in order to refine the estimated gain per subframe (4 bits every 5ms, ie 0.8 kbit / s) .
  • Artificial excitation u HB (n) is then filtered (block 111) by a transfer function (LPC) block LPC synthesis filter (I-1) HB (z) and operating at the sampling frequency of 16 kHz. The realization of this filter depends on the rate of the current frame:
  • ⁇ ⁇ ⁇ ⁇ ( ⁇ ) ⁇ ⁇ ⁇ - '( ⁇ ⁇ ⁇ )
  • the filter 1 / A HB (z) is of order 16 and simply corresponds to:
  • a low-pass filter also FIR type (block 113) is added to the treatment to further attenuate frequencies above 7 kHz.
  • the synthesis at high frequencies (HF) is finally added (block 130) to the low frequency synthesis (BF) obtained with the blocks 120 to 122 and resampled at 16 kHz (block 123).
  • the estimate of gains per subframe is not optimal. In part, it is based on an equalization of the "absolute" energy per sub-frame (block 101) between signals at different frequencies: the artificial excitation at 16 kHz (white noise) and a signal at 12.8 kHz ( ACELP excitation decoded).
  • the 7 kHz low-pass filter (block 113) introduces an offset of nearly 1 ms between the low and high bands, which can potentially degrade the quality of some signals by slightly desynchronizing the two bands at 23.85 kbit / s - this desynchronization can also be problematic when switching from 23.85 kbit / s to other modes.
  • FIGS. 2a global diagram
  • 2b gain prediction by response level correction
  • the input signal (mono) sampled at the frequency Fs (in Hz) is divided into two disjointed frequency bands, in which two LPC filters are calculated and coded separately:
  • a HF (z) Another LPC filter, denoted A HF (z), in the spectrally folded high band (Fs / 4 Fs / 2) - its quantized version is denoted by HF (z)
  • the band extension is done in the AMR-WB + code as detailed in sections 5.4 (HF coding) and 6.2 (HF decoding) of the 3GPP specification TS 26.290.
  • the extension consists in using the decoded excitation at low frequencies (LF excit.) And in shaping this excitation by a temporal gain per subframe (block 205) and a synthesis LPC filtering (block 207); in addition, enhancement processes (post-processing of the excitation (block 206) and smoothing of the energy of the reconstructed RF signal (block 208) are implemented as illustrated in FIG. 2a.
  • the technique of coding bandwidth gains in AMR-WB +, and more specifically the level compensation of LPC filters at their junction point, is a suitable method in the context of band-based LPC bandwidth expansion. low and high, and it can be noticed that such a level compensation between LPC filters is not present in the band extension of the AMR-WB codec.
  • the direct equalization of the level between the two LPC filters at the separation frequency is not an optimal method and can cause an overestimation of high-band energy and audible artifacts in some case; It is recalled that an LPC filter represents a spectral envelope, so the principle of equalizing the level between two LPC filters for a given frequency amounts to adjusting the relative level of two LPC envelopes.
  • the gain compensation in AMR-WB + is primarily a prediction of the known gain to the encoder and the decoder and which serves to reduce the bit rate necessary for the transmission of gain information scaling the excitation signal. high band.
  • the present invention improves the situation.
  • the invention provides a method for determining an optimized scale factor to be applied to an excitation signal or a filter during a frequency band extension process of an audio frequency signal, the band extension method comprising a step of decoding or extracting, in a first frequency band, an excitation signal and parameters of the first frequency band comprising coefficients of a linear prediction filter a step of generating an extended excitation signal on at least a second frequency band and a filtering step by a linear prediction filter for the second frequency band.
  • the determination method is such that it comprises the following steps:
  • an additional filter of a lower order than the linear prediction filter of the first frequency band, the coefficients of the additional filter being obtained from the parameters decoded or extracted from the first frequency band;
  • an additional filter of a lower order than the filter of the first frequency band to be equalized makes it possible to avoid overestimations of energy in the high frequencies which could result from local fluctuations of the envelope and which can disrupt the equalization of the prediction filters.
  • the band extension method comprises a step of applying the optimized scaling factor to the extended excitation signal.
  • the application of the optimized scaling factor is combined with the filtering step in the second frequency band.
  • the coefficients of the additional filter are obtained by truncation of the transfer function of the linear prediction filter of the first frequency band to obtain a lower order.
  • the coefficients of the additional filter are modified according to a criterion of stability of the additional filter.
  • the computation of the optimized scaling factor comprises the following steps:
  • the optimized scaling factor is calculated to avoid annoying artifacts that might arise in the event that the higher order filter frequency response of the first band near the common frequency would reveal a peak or a valley of the signal.
  • the method further comprises the following steps, implemented for a predetermined decoding rate:
  • the invention also relates to a device for determining an optimized scale factor to be applied to an excitation signal or to a filter in a frequency band extension device of an audiofrequency signal, the extension device band comprising a decoding or extraction module, in a first frequency band, of an excitation signal and parameters of the first frequency band comprising coefficients of a linear prediction filter, a generation module an excitation signal extended over at least a second frequency band and a filtering module by a linear prediction filter for the second frequency band.
  • the determination device is such that it comprises:
  • a module for determining a linear prediction filter called an additional filter, of a lower order than the linear prediction filter of the first frequency band, the coefficients of the additional filter being obtained from the parameters decoded or extracted from the first band; frequency; and
  • a module for calculating the scaling factor optimized according to at least the coefficients of the additional filter a module for calculating the scaling factor optimized according to at least the coefficients of the additional filter.
  • the invention relates to a decoder comprising a device as described.
  • It relates to a computer program comprising code instructions for implementing the steps of the method of determining an optimized scale factor as described, when these instructions are executed by a processor.
  • the invention relates to a storage medium, readable by a processor, integrated or not to the device for determining an optimized scale factor, possibly removable, storing a computer program implementing a method of determining a optimized scaling factor as previously described.
  • FIG. 1 illustrates a part of an AMR-WB decoder implementing frequency band extension steps of the state of the art and as previously described;
  • FIGS. 2a and 2b show the coding of the high band in the AMR code.
  • FIG. 3 illustrates an interoperable decoder with the AMR-WB coding and integrating a band extension device used according to one embodiment of the invention
  • FIG. 4 illustrates a device for determining a subframe-optimized scale factor as a function of the flow rate, according to one embodiment of the invention
  • FIGS. 5a and 5b illustrate the frequency responses of the filters used for calculating the optimized scaling factor according to one embodiment of the invention
  • FIG. 6 illustrates in flowchart form the main steps of a method for determining an optimized scale factor according to an embodiment of the invention
  • FIG. 7 illustrates a frequency domain embodiment of an optimized scale factor determination device during a band extension
  • FIG. 8 illustrates a hardware embodiment of an optimized scale factor determination device during a band extension according to the invention.
  • FIG. 3 illustrates an exemplary decoder, compatible with the AMR-WB / G.722.2 standard, in which there is a band extension comprising a determination of an optimized scale factor according to one embodiment of the method of the invention. implemented by the tape extender illustrated by block 309.
  • AMR-WB decoding which operates with an output sampling frequency of 16 kHz
  • CELP decoding (BF for low frequencies) always operates at the internal frequency of 12.8 kHz, as in AMR-WB, and the band extension (HF for high frequencies) used for the invention operates at the same time.
  • frequency of 16 kHz the synthesis BF and HF are combined (block 312) at the frequency fs after adequate resampling (block 306 and internal processing block 311).
  • the combination of the low and high bands can be done at 16 kHz, after resampling the low band of 12.8 to 16 kHz, before resampling the combined signal at the frequency fs.
  • the decoding according to FIG. 3 depends on the mode (or bit rate) AMR-WB associated with the current frame received.
  • This excitation u ⁇ n) is used in the adaptive dictionary of the following subframe; it is then post-processed and one distinguishes as in G.718 the excitation u ⁇ n) (also noted exc) of its modified post-processed version u (n) (also noted exc2) which serves as input to the filter of synthesis, 1 / ⁇ (z), in block 303.
  • the post-treatments applied to the excitation can be modified (for example, the phase dispersion can be improved) or these post-treatments can be extended (for example, inter-harmonic noise reduction can be implemented), without affecting the nature of the band extension. It can be noted that the use of blocks 306, 308, 314 is optional.
  • the decoding of the low band described above assumes a current frame called "active" with a rate between 6.6 and 23.85 kbit / s.
  • active a current frame
  • some frames can be coded as “inactive” and in this case you can either transmit a silence descriptor (on 35 bits) or not transmit anything.
  • SID frame describes several parameters: ISF parameters averaged over 8 frames, average energy over 8 frames, "dithering" flag for the non-stationary noise reconstruction.
  • the decoder makes it possible to extend the decoded low band (50-6400 Hz while taking into account 50 Hz high-pass filtering at the decoder, 0-6400 Hz in the general case) to an extended band whose width varies, ranging from approximately 50-6900 Hz to 50-7700 Hz depending on the mode implemented in the current frame.
  • the extension of the excitation is carried out in the frequency domain in a band of 5000 to 8000 Hz, to allow bandpass filtering of 6000 to 6900 or 7700 Hz width.
  • the HF gain correction information (0.8 kbit / s) transmitted at 23.85 kbit / s is decoded here. Its use is detailed below, with reference to FIG. 4.
  • the high band synthesis part is realized in block 309 representing the band extension device used for the invention and which is detailed in FIG. production.
  • a delay (block 310) is introduced to synchronize the outputs of the blocks 306 and 307 and the synthesized high band at 16 kHz is resampled from 16 kHz to the frequency fs (output of block 311).
  • the value of the delay T depends on how to synthesize the high band signal, the frequency fs as well as the post-processing of the low frequencies. Thus, in general, the value of T in the block 310 will have to be adjusted according to the specific implementation.
  • the low and high bands are then combined (added) in block 312 and the resulting synthesis is post-processed by high-order 50 Hz (type IIR) high-pass filtering whose coefficients depend on the frequency fs (block 313) and output post-processing with optional noise gate application similar to G.718 (block 314).
  • high-order 50 Hz type IIR
  • the block 400 from a decoded excitation signal in a first frequency band u (n), performs a band extension to obtain an extended excitation signal u HB (n) on at least a second frequency band.
  • the optimized scale factor estimation according to the invention is independent of how to obtain the signal u HB (n).
  • a condition regarding its energy is important, however. Indeed, it is necessary that the energy of the high band of 6000 to 8000 Hz is at a level similar to the energy of the band 4000 to 6000 Hz of the decoded excitation signal at the output of the block 302. since the low-band signal is de-emphasized (block 305), it is also necessary to apply the de-emphasis to the high-band excitation signal, either by using an own de-emphasis filter or by multiplying by a constant factor which corresponds to a mean attenuation of mentioned filter. This condition does not apply to the 23.85 kbit / s rate that uses the additional information transmitted by the encoder. In this case, the energy of the high band excitation signal must be consistent with the signal energy corresponding to the encoder, as explained later.
  • the frequency band extension may for example be implemented in the same way as for the AMR-WB type decoder described with reference to FIG. 1 in the blocks 100 to 102, from a white noise.
  • this band extension can be performed from a combination of a white noise and a decoded excitation signal as illustrated and subsequently described for blocks 700 to 707 of FIG. .
  • the tape expansion module may also be independent of the decoder and may extend a band of an existing audio signal stored or transmitted to the extension module, with an analysis of the audio signal to extract an excitation and an LPC filter.
  • the excitation signal at the input of the extension module is no longer a decoded signal but a signal extracted after analysis, as are the coefficients of the linear prediction filter of the first frequency band used in the method of determining the optimized scale factor in an implementation of the invention.
  • this calculation is preferably carried out by subframe and it consists in equalizing the levels of the frequency responses of the LPC filters 1 / ⁇ (z) and 1 / ⁇ ( ⁇ / ⁇ ) used at low frequencies and high frequencies, as described later with reference to Figure 7, with additional precautions to avoid the cases of overestimations that can result in too much energy of the synthesized high band and thus generate audible artifacts.
  • the determination of the optimized scale factor is also performed by the determination (in 401a) of a linear prediction filter called additional filter, of a lower order than the linear prediction filter of the first frequency band 1 / M z ) , the coefficients of the additional filter being obtained from the parameters decoded or extracted from the first frequency band.
  • the optimized scaling factor is then calculated (at 401b) based on at least one of these coefficients to be applied to the extended excitation signal u HB (n).
  • FIGS. 5a and 5b The principle of the determination of the optimized scaling factor implemented in block 401 is illustrated in FIGS. 5a and 5b with concrete examples obtained from signals sampled at 16 kHz; the amplitude values of frequency responses, noted later R, P, Q, of 3 filters are calculated at the common frequency of 6000 Hz (vertical dashed line) in the current subframe, whose index m n This is not recalled here in the notation of LPC filters interpolated by sub-frame to lighten the text.
  • the value of 6000 Hz is chosen so that it is close to the Nyquist frequency of the low band, ie 6400 Hz. It is preferable not to take this Nyquist frequency to determine the optimized scale factor.
  • the energy of the signal decoded at low frequencies is typically already attenuated at 6400 Hz.
  • the band extension described here is performed on a second so-called high band frequency band which ranges from 6000 to 8000 Hz. that in variants of the invention, another frequency that 6000 Hz can be chosen, without loss of generality to determine the optimized scale factor.
  • the two LPC filters are defined for separate bands (as in AMR-WB +). In this case R, P and Q will be calculated at the separation frequency.
  • FIGS 5a and 5b illustrate how the quantities R, P, Q are defined.
  • the first step is to calculate the R and P frequency responses respectively of the linear prediction filter of the first frequency band (low band) and the second frequency band (high band) at the frequency of 6000 Hz. 'on board :
  • M 16 is the order of the decoded LPC filter 1 / ⁇ (z), and ⁇ corresponds to the normalized 6000 Hz frequency for the sampling frequency of 12.8 kHz, ie:
  • the quantities P and R are calculated according to the following pseudocode:
  • the additional prediction filter is obtained for example by appropriately truncating the polynomial ⁇ (z) to order 2.
  • the stability of the filter 1 + can be verified in different ways, here a conversion is used in the domain of the PARCOR coefficients (or reflection coefficients) by calculating:
  • the second reflection coefficient, k 2 characterizes the resonance level of the signal model at order 2; since the use of a filter of order 2 aims at eliminating the influence of such resonances around the frequency of 6000 Hz, the value of k 2 is more strongly limited, this limit is fixed at 0.6.
  • the quantity Q calculated from the first 3 decoded LPC coefficients, takes better account of the influence of the spectral slope (or tilt) in the spectrum and avoids the influence of nearby "parasitic" peaks or valleys 6000 Hz which can bias or raise the value of the quantity R, calculated from all the LPC coefficients.
  • the optimized scaling factor is derived from the precomputed quantities R, P, Q conditionally as follows:
  • a smoothing is applied to the value of R.
  • an exponential smoothing is performed with a fixed factor in time (0.5) in the form:
  • R is the value of R in the preceding sub-frame and the factor 0.5 is empirically optimized - of course, the factor 0.5 can be changed to another value and other smoothing methods are also possible. Note that smoothing reduces temporal variations and therefore avoids artifacts.
  • g HB2 (m) mm (R, P, Q) / P
  • the smoothing of R may be replaced by a smoothing of g HB2 (m) as calculated above.
  • g HB (-i) is the scale factor or gain calculated for the last subframe of the previous frame.
  • the above tilt-only condition may be extended to take into account not only the tilt parameter but also other parameters to refine the decision.
  • the calculation of g HB2 (m) can be adjusted according to these additional parameters.
  • ZCR zero crossing rate
  • the zcr parameter usually gives the results similar to the tilt.
  • a good classification criterion is the ratio between zcr s calculated for the synthesized signal s ⁇ ri) and zcr u calculated for the excitation signal u (n) at 12800 Hz. This ratio is between 0 and 1, where 0 means that the signal has a decreasing spectrum, 1 that the spectrum is increasing (which corresponds to (1 - tilt) I ' 2.
  • a ratio zcr / zcr u > 0.5 corresponds to the case tilt ⁇ 0
  • a ratio zcr / zcr u ⁇ 0.5 corresponds to tilt> 0.
  • a function of a parameter tilt h may be used where tilt hp is the calculated tilt for the synthesized signal s (n) filtered by a high-pass filter with a cut-off frequency, for example at 4800 Hz; in this case, the 1 / ⁇ (z / ⁇ ) response of 6 to 8 kHz (applied at 16 kHz) corresponds to the weighted response of 1 / A (z) of 4.8 to 6.4 kHz. Since 1 / A (z / ⁇ ) has a more flattened response, you have to compensate for this change in tilt.
  • the scale factor function according to tilt hp is then given in one embodiment by:
  • the correction gain is calculated by comparing the energy of the sampled original signal at 16 kHz and filtered by a 6-7 kHz bandpass filter, s HB (n) with white noise energy at 16 kHz filtered by a synthesis filter 1 / A (z I ⁇ ) and a bandpass filter 6-7 kHz (before filtering the noise energy is set at a level similar to that of the excitation at 12.8 kHz), s HB2 (n).
  • the gain is the root of the energy ratio of the original signal on the energy of the noise divided by two.
  • the bandpass filter can be changed for a filter with a wider band (for example from 6 to 7.6 kHz).
  • block 404 scales the excitation signal according to the following equation:
  • g HB3 (m) is a subframe gain calculated in block 403 in the form:
  • index HF gain (m) is demultiplexed from the bitstream (block 405) and decoded by block 406 as follows:
  • HP _gain (.) is the HF gain quantization dictionary defined in the AMR-WB encoding and recalled below:
  • Block 407 scales the excitation signal according to the following equation:
  • the numerator represents here the band-high signal energy that would be obtained in the 23.05 mode.
  • the energy level between the decoded excitation signal and the extended excitation signal u HB (n) must be maintained, but this constraint is not necessary in the case of the 23.85 kbit / s rate, since u HB (n) is in this case scaled by the gain g HB3 (m).
  • certain multiplication operations applied to the signal in block 400 are applied in block 402 by multiplying by g (m).
  • g (m) depends on the synthesis algorithm of u HB (n) and must be adjusted so that the energy level between the decoded low band excitation signal and the g (m) signal u HB (n) is preserved.
  • g (m) 0.6 g HB1 (m)
  • g HB1 (m) is a gain which ensures, for the signal u HB , the same ratio between energy per subframe and energy per frame as for the signal u (n) and
  • 0.6 corresponds to the average value of frequency response amplitude of the deemphasis filter from 5000 to 6400 Hz.
  • this tilt is calculated as in the AMR-WB coding according to blocks 103 and 104, however other methods tilt estimation are possible without changing the principle of the invention.
  • Equalization is preferential at a frequency different from the Nyquist frequency (6400 Hz) associated with the low band.
  • the LPC modeling implicitly represents the attenuation of the signal typically caused by the resampling operations and therefore the frequency response of an LPC filter can be subjected to the frequency of Nyquist a decrease which is not found at the frequency chosen ses.
  • ⁇ Equalization is based on a lower order filter (here of order 2) in addition to 2 filters to equalize. This additional filter makes it possible to avoid the effects of local spectral fluctuations (peak or valley) that may be present at the common frequency for calculating the frequency response of the prediction filters.
  • the advantage of the invention is that the quality of the decoded signal at 23.85 kbit / s according to the invention is improved compared to a signal decoded at 23.05 kbit / s, which is not the case in an AMR-WB decoder.
  • this aspect of the invention makes it possible to use the additional information (0.8 kbit / s) received at 23.85 kbit / s, but in a controlled manner (block 408), to improve the quality of the excitation signal extended to flow rate of 23.85.
  • the optimized scale factor determining device as illustrated by blocks 401 to 408 of Fig. 4 implements an optimized scale factor determination method now described with reference to Fig. 6.
  • an extended excitation signal u HB (n) esX. obtained during a frequency band extension method E601 which comprises a step of decoding or extracting, in a first so-called low band frequency band, an excitation signal and parameters of the first band of such as, for example, the coefficients of the linear prediction filter of the first frequency band.
  • a step E602 determines a linear prediction filter called additional filter, of a lower order than that of the first frequency band. To determine this filter, the parameters of the first decoded or extracted frequency band are used.
  • this step is performed by truncation of the transfer function of the linear prediction filter of the low band to obtain a lower filter order, for example 2. These coefficients can then be modified according to a criterion of stability as explained above with reference to FIG. 4.
  • a step E603 is implemented to calculate the optimized scaling factor to be applied to the extended excitation signal.
  • This optimized scale factor is for example calculated from the frequency response of the additional filter at a common frequency between the low band (first frequency band) and the high band (second frequency band). A minimum value that can be chosen between the frequency response of this filter and those of the low band and high band filters.
  • This step of calculating the optimized scale factor is for example described above with reference to FIG. 4 and FIGS. 5a and 5b.
  • Step E604 performed by block 402 or 409 (depending on the decoding rate) for the band extension, applies the optimized scaling factor thus calculated to the extended excitation signal so as to obtain an extended extension signal.
  • the optimized scaling factor device 708 is integrated into a tape expansion device now described with reference to FIG. 7.
  • This optimized scale factor determining device illustrated by FIG. block 708 implements the method of determining the optimized scale factor described above with reference to FIG. 6.
  • the band extension block 400 of FIG. 4 comprises the blocks 700 to 707 of FIG. 7 now described.
  • a decoded or analytically estimated low band excitation signal is received (u (n)).
  • the band extension here uses the decoded excitation at 12.8 kHz (exc2 or u (n)) at the output of the block 302 of FIG. It should be noted that in this embodiment, the generation of the over-sampled and extended excitation occurs in a frequency band ranging from 5 to 8 kHz, thus including a second frequency band (6.4-8 kHz) greater than first frequency band (0-6.4 kHz).
  • the generation of an extended excitation signal is carried out at least second frequency band but also on a part of the first frequency band.
  • this signal is transformed to obtain an excitation signal spectrum U (k) by the time-frequency transformation module 500.
  • the DCT-IV transformation is implemented by FFT according to the "Evolved ZXT (EDCT)" algorithm described in the article by D. M. Zhang, HT. Li, A Low Complexity Transform - Evolved DCT, IEEE 14th International Conference on Computational Science and Engineering (CSE), Aug. 2011, pp. 144-149, and implemented in ITU-T G.718 Annex B and G.729.1 Annex E.
  • EDCT Evolved ZXT
  • the DCT-IV transformation may be replaced by other short-term time-frequency transformations of the same length and in the field of excitation, such as an FFT (for Fast Fourier Transform "in English) or DCT-II (Discrete Cosine Transform - Type II).
  • FFT Fast Fourier Transform
  • DCT-II Discrete Cosine Transform - Type II
  • MDCT for "Modified Discrete Cosine Tranform”
  • This approach preserves the original spectrum in this band and avoids introducing distortions in the 5000-6000 Hz band during the addition of HF synthesis with BF synthesis - particularly the signal phase (implicitly represented in the DCT-IV domain) in this band is preserved.
  • the 6000-8000 Hz band of U HB1 (k) is here defined by copying the 4000-6000 Hz band of U (k) since the value of start_band is preferably fixed at 160.
  • the value of start_band can be made adaptive around the value of 160.
  • the details of the adaptation of the value start_band are not described here because they go beyond the scope of the invention without changing the scope.
  • the high band may be noisy, harmonic or have a mixture of noise and harmonics.
  • the level of harmonicity in the 6000-8000 Hz band is generally correlated with that of the lower frequency bands.
  • the noise (in the 6000-8000 Hz band) is generated pseudo-randomly with a 16-bit linear congruent generator:
  • the combination block 703 can be realized in different ways.
  • the energy of the noise is calculated in three bands: 2000-4000 Hz, 4000-6000 Hz and 6000-8000 Hz, with
  • V 4 ⁇ u (k)
  • N (fc 1 , fc 2 ) is the set of indices & for which the index coefficient k is classified as being associated with noise. This set can be obtained for example by detecting the local peaks in U + and considering that these lines are not associated with noise, ie (by applying the negation of the previous condition):
  • N (a, b) ⁇ a ⁇ k ⁇ b ⁇ U + 1)
  • the ratio of the noise energy in the 4-6 kHz and 6-8 kHz bands is set so that between the 2-4 kHz and 4-6 kHz bands:
  • the calculation of a may be replaced by other methods.
  • the linear regression could for example be estimated so supervised by estimating the factor a by giving himself the original high band in a learning base. It will be noted that the method of calculating a does not limit the nature of the invention.
  • the factors ⁇ and a may be adapted to take account of the fact that noise injected into a given band of the signal is generally perceived as stronger than a harmonic signal at the same energy in the same band.
  • the factors ⁇ and a may be adapted to take account of the fact that noise injected into a given band of the signal is generally perceived as stronger than a harmonic signal at the same energy in the same band.
  • the block 703 realizes the equivalent of the block 101 of FIG. 1 in order to normalize the white noise as a function of an excitation which is on the other hand here in the frequency domain, already extended at the rate of 16 kHz; in addition, the mix is limited to the band 6000-8000 Hz.
  • the block 704 optionally carries out a dual operation of application of bandpass filter frequency response and deemphasis filtering (or deemphasis) in the frequency domain.
  • the deemphasis filtering may be performed in the time domain, after block 705 or even before block 700; however, in this case, bandpass filtering performed in block 704 may leave some low frequency components of very low levels that are amplified by de-emphasis, which may slightly discern the decoded low band. For this reason, it is preferred here to perform the deemphasis in the frequency domain.
  • G deemph (k) is the frequency response of the filter l / (l - 0.68z _1 ) over a restricted discrete frequency band.
  • G deem h (k) as:
  • the definition of O k can be adjusted (for example for even frequencies).
  • the HF synthesis is not de-emphasized.
  • the high frequency signal is on the contrary de-emphasized so as to bring it back to a domain coherent with the low frequency signal (0-6.4 kHz) which leaves block 305 of FIG. 3. This is important for the estimation and the subsequent adjustment of the energy of the HF synthesis.
  • the de-emphasis can be performed in an equivalent way in the time domain after inverse DCT.
  • band-pass filtering is applied with two separate parts: one fixed high-pass, the other adaptive low-pass (flow-rate function).
  • This filtering is performed in the frequency domain.
  • the partial low-pass filter response in the frequency domain is calculated as follows:
  • N lp 60 to 6.6 kbit / s, 40 to 8.85 kbit / s, 20 at rates> 8.85 bit / s.
  • bandpass filtering can be adapted by defining a single filtering step combining the high-pass and low-pass filtering.
  • the bandpass filtering may be performed equivalently in the time domain (as in block 112 of FIG. 1) with different filter coefficients according to the bit rate, after an inverse DCT step.
  • it is advantageous to carry out this step directly in the frequency domain because the filtering is carried out in the field of LPC excitation and therefore the problems of circular convolution and edge effects are very limited in this field. .
  • the inverse transform block 705 performs an inverse DCT on 320 samples to find the high-frequency excitation sampled at 16 kHz. Its implementation is identical to block 700, because the DCT-IV is orthonormed, except that the length of the transform is 320 instead of 256, and we obtain:
  • This excitation sampled at 16 kHz is then optionally scaled by gains defined by subframe of 80 samples (block 707).
  • Block 707 scales the combined signal according to the following equation:
  • the realization of the block 706 differs from that of the block 101 of Figure 1, because the energy at the current frame is taken into account in addition to that of the sub-frame. This makes it possible to have the ratio of the energy of each sub-frame with respect to the energy of the frame. Energy ratios (or relative energies) are compared rather than the absolute energies between low band and high band.
  • this scaling step makes it possible to keep in the high band the energy ratio between the subframe and the frame in the same way as in the low band.
  • the block 708 then performs a scaling factor calculation by subframe of the signal (steps E602 to E 603 of FIG. 6), as previously described with reference to FIG. 6 and detailed in FIG. 5.
  • this filtering can be done in the same way as described for the block 111 of FIG. 1 of the AMR-WB decoder, however the order of the filter goes to 20 at the rate of 6.6, which does not change. not significantly the quality of the synthesized signal.
  • the step of filtering by a linear prediction filter 710 for the second frequency band is combined with the application of the optimized scaling factor, which reduces the processing complexity.
  • the filtering steps 1 / ⁇ ( ⁇ / ⁇ ) and the application of the optimized scaling factor g HB2 are combined with a single filter step g HB2 1 (z / y) to reduce the processing complexity.
  • the coding of the low band (0-6.4 kHz) may be replaced by a CELP coder other than that used in AMR-WB, for example the CELP coder in G.718 to 8. kbit / s.
  • a CELP coder other than that used in AMR-WB, for example the CELP coder in G.718 to 8. kbit / s.
  • other encoders in wide band or operating at frequencies higher than 16 kHz in which the coding of the low band operates at an internal frequency at 12.8 kHz could be used.
  • the invention can be obviously adapted to other sampling frequencies than 12.8 kHz, when a low frequency encoder operates at a sampling frequency lower than that of the original or reconstructed signal.
  • the low band decoding does not use a linear prediction, it does not have an excitation signal to be extended, in this case it will be possible to carry out an LPC analysis of the reconstructed signal in the current frame and calculate an LPC excitation. so as to be able to apply the invention.
  • the excitation (u (n)) is resampled, for example by linear interpolation or cubic "spline", from 12.8 to 16 kHz before transformation (for example DCT-IV) of length 320.
  • This variant has the defect of being more complex, because the transform (DCT-IV) of the excitation is then calculated over a greater length and the resampling is not carried out in the field of the transform.
  • all the calculations necessary for estimating the gains can be carried out in a field logarithmic.
  • the low band excitation u (n) and the LPC 1 / A (z) filter will be estimated per frame, by LPC analysis of a low band signal whose band must be extended. The low band excitation signal is then extracted by analyzing the audio signal.
  • the low band audio signal is resampled before the excitation extraction step, so that the excitation extracted from the audio signal (by linear prediction) is already resolved. sampled.
  • the band extension illustrated in Figure 7, in this case applies to a low band which is not decoded but analyzed.
  • FIG. 8 represents an exemplary embodiment of a device for determining an optimized scale factor 800 according to the invention. This may be an integral part of an audio-frequency signal decoder or equipment receiving decoded or non-decoded audio signals.
  • This type of device comprises a PROC processor cooperating with a memory block BM having a memory storage and / or work MEM.
  • Such a device comprises an input module E able to receive a decoded or extracted excitation audio signal in a first so-called low band frequency band (u (n) or U (k)) and the parameters of a filter of linear prediction synthesis (A (z)). It comprises an output module S adapted to transmit the synthesized and optimized high frequency signal (U H B '(I ") for example to a filtering module such as block 710 of FIG. 7 or to a resampling module as the module 311 of FIG.
  • the memory block may advantageously comprise a computer program comprising code instructions for implementing the steps of the method for determining an optimized scale factor to be applied to an excitation signal or to a filter within the meaning of FIG. invention, when these instructions are executed by the processor PROC, and in particular the steps of determination (E602) of a linear prediction filter called additional filter, of order less than the linear prediction filter of the first frequency band, the coefficients additional filter being obtained from the parameters decoded or extracted from the first frequency band, calculation (E603) of an optimized scale factor according to at least the coefficients of the additional filter.
  • FIG. 6 shows the steps of an algorithm of such a computer program.
  • the computer program can also be stored on a memory medium readable by a reader of the device or downloadable in the memory space thereof.
  • the memory MEM generally records all the data necessary for the implementation of the method.
  • the device thus described may also include the functions of applying the optimized scaling factor to the extended excitation signal, frequency band extension, low band decoding and other processing functions. described for example in FIGS. 3 and 4 in addition to the optimized scale factor determination functions according to the invention.

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Signal Processing (AREA)
  • Audiology, Speech & Language Pathology (AREA)
  • Human Computer Interaction (AREA)
  • Computational Linguistics (AREA)
  • Acoustics & Sound (AREA)
  • Multimedia (AREA)
  • Quality & Reliability (AREA)
  • Mathematical Physics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Compression, Expansion, Code Conversion, And Decoders (AREA)
  • Organic Low-Molecular-Weight Compounds And Preparation Thereof (AREA)

Abstract

The invention relates to a method for determining an optimized scale factor to be applied to an excitation signal or to a filter during a process for frequency band extension of an audiofrequency signal, the band extension process (E601) comprising a step of decoding or extracting, in a first frequency band, an excitation signal and parameters of the first frequency band including the coefficients of a linear prediction filter, a step of generating an excitation signal extending over at least one second frequency band, and a step of filtering by means of a linear prediction filter for the second frequency band. The determination method comprises the steps of determining (E602) a linear prediction filter referred to as an additional filter, of a lower order than that of the linear prediction filter of the first frequency band, the coefficients of the additional filter being obtained from the parameters decoded or extracted from the first frequency band and calculating (E603) the optimized scale factor as a function of at least the coefficients of the additional filter. The invention also relates to a device for determining an optimized scale factor using the method as described and to a decoder including such a device.

Description

Facteur d'échelle optimisé pour l'extension de bande de fréquence dans un décodeur de signaux audiofréquences  Scaling factor optimized for frequency band extension in an audio-frequency decoder
La présente invention se rapporte au domaine du codage/décodage et du traitement de signaux audiofréquences (comme des signaux de parole, de musique ou autres) pour leur transmission ou leur stockage. The present invention relates to the field of coding / decoding and audio-frequency signal processing (such as speech, music or other signals) for their transmission or storage.
Plus particulièrement, l'invention concerne un procédé et un dispositif de détermination d'un facteur d'échelle optimisé servant à ajuster le niveau d'un signal d'excitation ou de façon équivalent d'un filtre lors d'une extension de bande de fréquence dans un décodeur ou un processeur réalisant une amélioration de signal audiofréquence.  More particularly, the invention relates to a method and an apparatus for determining an optimized scale factor for adjusting the level of an excitation signal or, in a similar manner, a filter during a band extension. frequency in a decoder or a processor performing audio-frequency signal enhancement.
De nombreuses techniques existent pour compresser (avec perte) un signal audiofréquence comme la parole ou la musique.  Many techniques exist to compress (with loss) an audiofrequency signal such as speech or music.
Les méthodes classiques de codage pour les applications conversationnelles sont en général classifiées en codage de forme d'onde (MIC pour "Modulation par Impulsion et codage", MICDA pour "Modulation par Impulsion et Codage Différentiel Adaptatif, codage par transformée...), codage paramétrique (LPC pour "Linear Prédictive Coding" en anglais, codage sinusoïdal...) et codage hybride paramétrique avec une quantification des paramètres par "analyse par synthèse" dont le codage CELP (pour "Code Excited Linear Prédiction" en anglais) est l'exemple le plus connu.  Conventional methods of coding for conversational applications are generally classified in waveform coding (MIC for "pulse modulation and coding", ADPCM for "Pulse Modulation and Adaptive Differential Coding, transform coding, etc.), parametric coding (LPC for "Linear Predictive Coding" in English, sinus coding ...) and parametric hybrid coding with a quantification of the parameters by "analysis by synthesis" whose coding CELP (for "Code Excited Linear Prediction" in English) is the best known example.
Pour les applications non conversationnelles, l'état de l'art en codage de signal audio For non-conversational applications, the state of the art in audio signal coding
(mono) est constitué par le codage perceptuel par transformée ou en sous-bandes, avec un codage paramétrique des hautes fréquences par réplication de bande. (mono) is constituted by the perceptual coding by transform or in sub-bands, with a parametric coding of the high frequencies by band replication.
Une revue des méthodes classiques de codage de parole et audio se trouve dans les ouvrages W.B. Kleijn and K.K. Paliwal (Eds.), Speech Coding and Synthesis, Elsevier, 1995 ; M. Bosi, R.E. Goldberg, Introduction to Digital Audio Coding and Standards, Springer 2002 ; J. Benesty, M. M. Sondhi, Y. Huang (Eds.), Handbook of Speech Processing, Springer 2008. A review of conventional speech and audio coding methods can be found in W. B. Kleijn and K. K. Paliwal (Eds.), Speech Coding and Synthesis, Elsevier, 1995; M. Bosi, R. E. Goldberg, Introduction to Digital Audio Coding and Standards, Springer 2002; J. Benesty, M. Sondhi M., Y. Huang (Eds.), Handbook of Speech Processing, Springer 2008.
On s'intéresse ici plus particulièrement au codée (codeur et décodeur) normalisé 3GPP AMR-WB (pour "Adaptive Multi-Rate Wideband" en anglais) qui fonctionne à une fréquence d'entrée/sortie de 16 kHz et dans lequel le signal est divisé en deux sous-bandes, la bande basse (0-6.4 kHz) qui est échantillonnée à 12.8 kHz et codée par modèle CELP et la bande haute (6.4-7 kHz) qui est reconstruite de façon paramétrique par « extension de bande » (ou BWE pour "Bandwidth Extension" en anglais) avec ou sans information supplémentaire selon le mode de la trame courante. On peut noter ici que la limitation de la bande codée du codée AMR-WB à 7kHz est essentiellement liée au fait que la réponse en fréquence en émission des terminaux en bande élargie a été approximée au moment de la normalisation (ETSI/3GPP puis UIT-T) selon le masque fréquentiel défini dans la norme UIT-T P.341 et plus précisément en utilisant un filtre dit « P341 » défini dans la norme UIT-T G.191 qui coupe les fréquences au-dessus de 7 kHz (ce filtre respecte le masque défini dans P.341). Cependant, en théorie, il est bien connu qu'un signal échantillonné à 16 kHz peut avoir une bande audio définie de 0 à 8000 Hz ; le codée AMR-WB introduit donc une limitation de la bande haute en comparaison à la largeur de bande théorique de 8 kHz. Of particular interest here is the 3GPP AMR-WB ("Adaptive Multi-Rate Wideband") codec and decoder, which operates at an input / output frequency of 16 kHz and in which the signal is divided into two sub-bands, the low band (0-6.4 kHz) which is sampled at 12.8 kHz and coded by CELP model and the high band (6.4-7 kHz) which is parametrically reconstructed by "band extension" ( or BWE for "Bandwidth Extension" with or without additional information depending on the mode of the current frame. It can be noted here that the limitation of the coded band of the AMR-WB codec at 7 kHz is essentially related to the fact that the transmission frequency response of the broadband terminals has been approximated at the time of standardization (ETSI / 3GPP then ITU-T T) according to the frequency mask defined in the ITU-T P.341 standard and more precisely by using a so-called "P341" filter defined in the ITU-T G.191 standard. which cuts frequencies above 7 kHz (this filter respects the mask defined in P.341). However, in theory, it is well known that a signal sampled at 16 kHz may have a defined audio band of 0 to 8000 Hz; the AMR-WB codec thus introduces a limitation of the high band in comparison with the theoretical bandwidth of 8 kHz.
Le codée de parole 3GPP AMR-WB a été normalisé en 2001 principalement pour les applications de téléphonie en mode circuit (CS) sur GSM (2G) et UMTS (3G). Ce même codée a été aussi normalisé en 2003 à l'UIT-T en tant que recommandation G.722.2 "Wideband coding speech at around 16kbit/s using Adaptive Multi-Rate Wideband (AMR-WB)".  The 3GPP AMR-WB speech codec was standardized in 2001 mainly for circuit-mode (CS) telephony applications on GSM (2G) and UMTS (3G). This same codec was also standardized in 2003 in ITU-T as Recommendation G.722.2 "Wideband coding speech at around 16kbit / s using Adaptive Multi-Rate Wideband (AMR-WB)".
Il comprend neuf débits, appelés modes, de 6.6 à 23.85 kbit/s, et comprend des mécanismes de transmission continue (DTX pour "Discontinuous Transmission") avec détection d'activité vocale (VAD pour "Voice Activity Détection") et génération de bruit de confort (CNG pour "Confort Noise Génération") à partir de trames de description de silence (SID pour "Silence Insertion Descriptor"), ainsi que des mécanismes de correction de trames perdues (FEC pour "Frame Erasure Concealment", parfois appelé PLC pour "Packet Loss Concealment").  It includes nine speeds, called modes, from 6.6 to 23.85 kbit / s, and includes continuous transmission mechanisms (DTX for "Discontinuous Transmission") with voice activity detection (VAD for "Voice Activity Detection") and noise generation of comfort (CNG for "Comfort Noise Generation") from frames of silence description (SID for "Silence Insertion Descriptor"), as well as mechanisms of correction of lost frames (FEC for "Frame Erasure Concealment", sometimes called PLC for "Packet Loss Concealment").
On ne reprend pas ici les détails de l'algorithme de codage et de décodage AMR-WB, une description détaillée de ce codée se trouve dans les spécifications 3GPP (TS 26.190, 26.191, 26.192, 26.193, 26.194, 26.204) et UIT-T-G.722.2 (et les Annexes et Appendice correspondantes) ainsi que dans l'article de B. Bessette et al. intitulé « The adaptive multirate wideband speech codée (AMR-WB) », IEEE Transactions on Speech and Audio Processing, vol. 10, no. 8, 2002, pp. 620-636 et les codes sources des standards 3GPP et UIT-T associés.  The details of the AMR-WB coding and decoding algorithm are not repeated here, a detailed description of this codec is found in the 3GPP specifications (TS 26.190, 26.191, 26.192, 26.193, 26.194, 26.204) and ITU-TG .722.2 (and the corresponding Appendices and Appendix) and in the article by B. Bessette et al. entitled "The adaptive multirate broadband speech coded (AMR-WB)", IEEE Transactions on Speech and Audio Processing, vol. 10, no. 8, 2002, pp. 620-636 and associated 3GPP and ITU-T standard source codes.
Le principe de l'extension de bande dans le codée AMR-WB est assez rudimentaire. En effet, la bande haute (6.4-7 kHz) est générée en mettant en forme un bruit blanc par le biais d'une enveloppe temporelle (appliquée sous la forme de gains par sous-trame) et fréquentielle (par l'application d'un filtre de synthèse de prédiction linéaire ou LPC pour "Linear Prédictive Coding"). Cette technique d'extension de bande est illustrée à la figure 1.  The principle of band extension in the AMR-WB codec is rather rudimentary. Indeed, the high band (6.4-7 kHz) is generated by formatting a white noise through a temporal envelope (applied in the form of gains per subframe) and frequency (by the application of a linear prediction synthesis filter or LPC for "Linear Predictive Coding"). This band extension technique is illustrated in Figure 1.
Un bruit blanc, uHB1 (n) , n = 0, · · · , 79 , est généré à 16 kHz par sous-trame de 5 ms par générateur congruentiel linéaire (bloc 100). Ce bruit uHB1 (n) est mis en forme dans le temps par application de gains par sous-trame ; cette opération est décomposée en deux étapes de traitement (blocs 102, 106 ou 109) : A white noise, u HB1 (n), n = 0, · · ·, 79, is generated at 16 kHz per 5 ms subframe per linear congruent generator (block 100). This noise u HB1 (n) is shaped in time by applying gains per subframe; this operation is broken down into two processing steps (blocks 102, 106 or 109):
• Un premier facteur est calculé (bloc 101) pour mettre le bruit blanc uHB1 (n) (bloc• A first factor is calculated (block 101) to set the white noise u HB1 (n) (block
102) à un niveau semblable à celui de l'excitation, u(n) , « = (), · · - , 63 , décodée à 12.8 kHz dans la bande basse : 102) at a level similar to that of the excitation, u (n), "= (), · · -, 63, decoded at 12.8 kHz in the low band:
On peut noter ici que la normalisation des énergies se fait en comparant des blocs de taille différente (64 pour u(n) et 80 pour uHB1 (n) ), sans compensation des différences de fréquences d'échantillonnage (12.8 ou 16 kHz). It can be noted here that the normalization of the energies is done by comparing blocks of different size (64 for u (n) and 80 for u HB1 (n)), without compensation of the differences of sampling frequencies (12.8 or 16 kHz) .
• L'excitation dans la bande haute est ensuite obtenue (bloc 106 ou 109) sous la forme :  • The excitation in the high band is then obtained (block 106 or 109) in the form:
uHB {n) = gHBuHB2 {n) u HB {n) = g HB u HB2 {n)
où le gain gHB est obtenu différemment selon le débit. Si le débit de la trame actuelle est <23.85 kbit/s, le gain gHB est estimé « en aveugle » (c'est-à-dire sans information supplémentaire); dans ce cas, le bloc 103 filtre le signal décodé en bande basse par un filtre passe-haut ayant une fréquence de coupure à 400 Hz pour obtenir un signal shp (n) , « = (), · · · , 63 - ce filtre passe-haut élimine l'influence des très basses fréquences qui peuvent biaiser l'estimation faite dans le bloc 104 - puis on calcule le « tilt » (indicateur de pente spectrale) noté etilt du signal sh (n) par autocorrélation normalisée (bloc 104): where the gain g HB is obtained differently depending on the flow rate. If the rate of the current frame is <23.85 kbit / s, the gain g HB is estimated to be "blind" (ie without additional information); in this case, the block 103 filters the low-band decoded signal by a high-pass filter having a cut-off frequency at 400 Hz to obtain a signal s hp (n), "= (), · · ·, 63 - this high-pass filter eliminates the influence of the very low frequencies that can bias the estimate made in block 104 - then the "tilt" (spectral slope indicator) noted e tilt of the signal s h (n) by standard autocorrelation is calculated (block 104):
63 et enfin on calcule gHB sous la forme : 63 and finally we calculate g HB in the form:
S m = wspSsp + ^ - wSP )gBG S m = w spSsp + ^ - w SP ) g BG
où gsp = l - etilt est le gain appliqué dans les trames actives de parole (SP pour speech), gBG = l .25gsp est le gain appliqué dans les trames inactives de parole associées à un bruit de fond (BG pour Background) et wsp est une fonction de pondération qui dépend de la détection d'activité vocale (VAD). On comprend que l'estimation du tilt ( etilt ) permet d'adapter le niveau de la bande haute en fonction de la nature spectrale du signal ; cette estimation est particulièrement importante quand la pente spectrale du signal décodé CELP est telle que l'énergie moyenne décroît quand la fréquence augmente (cas d'un signal voisé où etilt est proche de 1, donc gsp = l - etilt est ainsi réduit). A noter aussi que le facteur gHB dans le décodage AMR-WB est borné pour prendre des valeurs dans l'intervalle [0.1, 1.0]. En effet, pour les signaux dont énergie croît quand la fréquence augmente ( etilt proche de -1, gsp proche de 2), le gain gHB est d'habitude sous-estimé. where g sp = l - e tilt is the gain applied in active speech frames (SP for speech), g BG = l .25g sp is the gain applied in inactive speech frames associated with background noise (BG for Background) and w sp is a weighting function that depends on Voice Activity Detection (VAD). It is understood that the estimate of the tilt (e tilt ) makes it possible to adapt the level of the high band according to the spectral nature of the signal; this estimate is particularly important when the spectral slope of the decoded signal CELP is such that the average energy decreases when the frequency increases (case of a voiced signal where e tilt is close to 1, so g sp = l - e tilt is thus reduced). Note also that the factor g HB in the decoding AMR-WB is bounded to take values in the interval [0.1, 1.0]. Indeed, for signals whose energy increases when the frequency increases (e tilt close to -1, g sp close to 2), the gain g HB is usually underestimated.
A 23.85 kbit/s, une information de correction est transmise par le codeur AMR-WB et décodée (blocs 107, 108) afin d'affiner le gain estimé par sous-trame (4 bits toutes les 5ms, soit 0.8 kbit/s). L'excitation artificielle uHB (n) est ensuite filtrée (bloc 111) par un filtre de synthèse LPC (bloc 111) de fonction de transfert \ l AHB (z) et fonctionnant à la fréquence d'échantillonnage de 16 kHz. La réalisation de ce filtre dépend du débit de la trame courante:At 23.85 kbit / s, correction information is transmitted by the encoder AMR-WB and decoded (blocks 107, 108) in order to refine the estimated gain per subframe (4 bits every 5ms, ie 0.8 kbit / s) . Artificial excitation u HB (n) is then filtered (block 111) by a transfer function (LPC) block LPC synthesis filter (I-1) HB (z) and operating at the sampling frequency of 16 kHz. The realization of this filter depends on the rate of the current frame:
• A 6.6 kbit/s, le filtre \ l AHB (z) est obtenu en pondérant par un facteur ^=0.9 un filtre LPC d'ordre 20, 1 / Aext (z) qui « extrapole » le filtre LPC d'ordre 16, 1 / Â(z) , décodé dans la bande basse (à 12.8 kHz) - les détails de l'extrapolation dans le domaine des paramètres ISF (pour "Imittance Spectral Frequency" en anglais) sont décrits dans la norme G.722.2 à la section 6.3.2.1; dans ce cas, • At 6.6 kbit / s, the filter \ l A HB (z) is obtained by weighting by a factor ^ = 0.9 an LPC filter of order 20, 1 / A ext (z) which "extrapolates" the LPC filter from order 16, 1 / Â (z), decoded in the low band (at 12.8 kHz) - details of the extrapolation in the domain of the ISF parameters (for "Imittance Spectral Frequency") are described in the standard G. 722.2 in section 6.3.2.1; in that case,
\Ι ΑΗΒ (ζ) = \ Ι Α-' (ζ Ι γ) \ Ι Α ΗΒ ( ζ ) = \ Ι Α- '( ζ Ι γ)
· Aux débits > 6.6 kbit/s, le filtre 1 / AHB (z) est d'ordre 16 et correspond simplement à :· At rates> 6.6 kbit / s, the filter 1 / A HB (z) is of order 16 and simply corresponds to:
où ^=0.6. A noter que dans ce cas le filtre 1 / A(z l Y) est utilisé à 16 kHz, ce qui résulte en un étalement (par homothétie) de la réponse en fréquence de ce filtre de [0, 6.4 kHz] à [0, 8 kHz].  where ^ = 0.6. Note that in this case the filter 1 / A (zl Y) is used at 16 kHz, which results in a spread (by homothety) of the frequency response of this filter from [0, 6.4 kHz] to [0, 8 kHz].
Le résultat, sHB (n) , est enfin traité par un filtre passe-bande (bloc 112) de type FIR ("FiniteThe result, s HB (n), is finally processed by a band-pass filter (block 112) of type FIR ("Finite
Impulse Response"), pour ne garder que la bande 6 - 7 kHz ; à 23.85 kbit/s, un filtre passe- bas également de type FIR (bloc 113) se rajoute au traitement pour atténuer encore plus les fréquences supérieures à 7 kHz. La synthèse en hautes fréquences (HF) est finalement additionnée (bloc 130) à la synthèse en basses fréquences (BF) obtenue avec les blocs 120 à 122 et ré-échantillonnée à 16 kHz (bloc 123). Ainsi même si la bande haute s'étend en théorie de 6.4 à 7 kHz dans le codée AMR-WB, la synthèse HF est plutôt comprise dans la bande 6-7 kHz avant addition avec la synthèse BF. Impulse Response "), to keep only the band 6 - 7 kHz, at 23.85 kbit / s, a low-pass filter also FIR type (block 113) is added to the treatment to further attenuate frequencies above 7 kHz. The synthesis at high frequencies (HF) is finally added (block 130) to the low frequency synthesis (BF) obtained with the blocks 120 to 122 and resampled at 16 kHz (block 123). theoretically extends from 6.4 to 7 kHz in the AMR-WB coding, the HF synthesis is rather in the band 6-7 kHz before addition with the synthesis BF.
On peut identifier plusieurs inconvénients à la technique d'extension de bande du codée AMR-WB, en particulier :  Several disadvantages can be identified with the AMR-WB coding band extension technique, in particular:
• L'estimation de gains par sous-trame (bloc 101, 103 à 105) n'est pas optimale. Pour partie, elle se base sur une égalisation de l'énergie « absolue » par sous-trame (bloc 101) entre des signaux à des fréquences différentes : l'excitation artificielle à 16 kHz (bruit blanc) et un signal à 12.8 kHz (excitation ACELP décodée). On peut noter en particulier que cette approche induit implicitement une atténuation de l'excitation bande haute (par un ratio 12.8/16=0.8) ; en fait, on notera également qu'aucune désaccentuation (ou déemphase) n'est effectuée sur la bande haute dans le codée AMR-WB, ce qui induit implicitement une amplification relative proche de 0.6 (qui correspond à la valeur de la réponse en fréquence de l / (l - 0.68z_1 ) à 6400 Hz). En fait, les facteurs de 1/0.8 et de 0.6 se compensent approximativement. • The estimate of gains per subframe (block 101, 103 to 105) is not optimal. In part, it is based on an equalization of the "absolute" energy per sub-frame (block 101) between signals at different frequencies: the artificial excitation at 16 kHz (white noise) and a signal at 12.8 kHz ( ACELP excitation decoded). It may be noted in particular that this approach implicitly induces an attenuation of the excitation band high (by a ratio 12.8 / 16 = 0.8); in fact, it will also be noted that no deemphasis is performed on the high band in the AMR-WB coding, which implicitly induces a relative amplification close to 0.6 (which corresponds to the value of the frequency response s / (l - 0.68z _1) to 6400 Hz). In fact, the factors of 1 / 0.8 and 0.6 compensate each other approximately.
• Sur la parole, les tests de caractérisation du codée 3GPP AMR-WB documentés dans le rapport 3GPP TR 26.976 ont montré que le mode à 23.85 kbit/s a une qualité moins bonne qu'à 23.05 kbit/s, sa qualité est en fait similaire à celle du mode à 15.85 kbit/s. Ceci montre en particulier que le niveau du signal HF artificiel doit être contrôlé de façon très prudente, car la qualité est dégradée à 23.85 kbit/s alors que les 4 bits par trame sont sensés permettre de mieux approcher l'énergie des hautes fréquences originales.  • On speech, the 3GPP AMR-WB coding characterization tests documented in the 3GPP TR 26.976 report showed that the 23.85 kbit / sa mode is not as good as 23.05 kbit / s, its quality is in fact similar to that of the 15.85 kbit / s mode. This shows in particular that the level of artificial RF signal must be controlled very carefully, because the quality is degraded to 23.85 kbit / s while the 4 bits per frame are supposed to better approach the energy of the original high frequencies.
• Le filtre passe-bas à 7 kHz (bloc 113) introduit un décalage de près de 1 ms entre les bandes basses et hautes, ce qui peut potentiellement dégrader la qualité de certains signaux en désynchronisant légèrement les deux bandes à 23.85 kbit/s - cette désynchronisation peut également poser problème lors d'une commutation de débit de 23.85 kbit/s à d'autres modes.  • The 7 kHz low-pass filter (block 113) introduces an offset of nearly 1 ms between the low and high bands, which can potentially degrade the quality of some signals by slightly desynchronizing the two bands at 23.85 kbit / s - this desynchronization can also be problematic when switching from 23.85 kbit / s to other modes.
Un exemple d'extension de bande par approche temporelle est décrit dans la norme 3GPP TS 26.290 décrivant le codée AMR-WB+ (normalisé en 2005). Cet exemple est illustré dans les schémas-blocs des figures 2a (schéma global) et 2b (prédiction de gain par correction de niveau de réponses) qui correspondent respectivement aux figures 16 et 10 de la spécification 3GPP TS 26.290. An example of a time approach band extension is described in the 3GPP TS 26.290 standard describing the AMR-WB + coding (standardized in 2005). This example is illustrated in the block diagrams of FIGS. 2a (global diagram) and 2b (gain prediction by response level correction) which respectively correspond to FIGS. 16 and 10 of the specification 3GPP TS 26.290.
Dans le codée AMR-WB+, le signal d'entrée (mono) échantillonné à la fréquence Fs (en Hz) est divisés en deux bandes de fréquences disjointes, dans lesquelles deux filtres LPC sont calculés et codés séparément:  In the AMR-WB + codec, the input signal (mono) sampled at the frequency Fs (in Hz) is divided into two disjointed frequency bands, in which two LPC filters are calculated and coded separately:
• un filtre LPC, noté A(z) , dans la bande basse (0-Fs/4) - sa version quantifiée est notée Â(z)  • an LPC filter, denoted A (z), in the low band (0-Fs / 4) - its quantized version is denoted ((z)
• un autre filtre LPC, noté AHF (z) , dans la bande haute repliée spectralement (Fs/4- Fs/2) - sa version quantifiée est notée ÂHF (z) Another LPC filter, denoted A HF (z), in the spectrally folded high band (Fs / 4 Fs / 2) - its quantized version is denoted by HF (z)
L'extension de bande se fait dans le codée AMR-WB+ comme détaillé dans les sections 5.4 (codage HF) et 6.2 (décodage HF) de la spécification 3GPP TS 26.290. On en résume ici le principe : l'extension consiste à utiliser l'excitation décodée en basses fréquences (LF excit.) et à mettre en forme cette excitation par un gain temporel par sous-trame (bloc 205) et un filtrage LPC de synthèse (bloc 207) ; de plus, des traitements d'améliorations (post-traitement de l'excitation (bloc 206) et lissage de l'énergie du signal HF reconstruit (bloc 208) sont mis en œuvre comme illustré à la figure 2a. Il est important de remarquer que cette extension dans AMR-WB+ nécessite la transmission d'informations supplémentaires : les coefficients du filtre ÂHF (z) en 204 et un gain de mise en forme temporelle par sous-trame (bloc 201). Une particularité de l'algorithme d'extension de bande dans AMR-WB+ est que le gain par sous-trame est quantifié par une approche prédictive ; autrement dit, on ne code pas les gains directement, mais plutôt des corrections de gain qui sont relatives à une estimation du gain notée Cette estimation, correspondant en fait à un facteur d'égalisation du niveau entre les filtres Â(z) et AHF (z) à la fréquence de séparation entre bande basse et bande haute (Fs/4). Le calcul du facteur (bloc 203) est détaillé à la figure 10 de la spécification 3GPP TS 26.290 reprise ici à la figure 2b. On ne détaillera pas plus ici cette figure. On retiendra pour résumer que les blocs The band extension is done in the AMR-WB + code as detailed in sections 5.4 (HF coding) and 6.2 (HF decoding) of the 3GPP specification TS 26.290. Here we summarize the principle: the extension consists in using the decoded excitation at low frequencies (LF excit.) And in shaping this excitation by a temporal gain per subframe (block 205) and a synthesis LPC filtering (block 207); in addition, enhancement processes (post-processing of the excitation (block 206) and smoothing of the energy of the reconstructed RF signal (block 208) are implemented as illustrated in FIG. 2a. It is important to note that this extension in AMR-WB + requires the transmission of additional information: the coefficients of the filter  HF (z) 204 and a temporal shaping gain per subframe (block 201). A peculiarity of the band extension algorithm in AMR-WB + is that subframe gain is quantified by a predictive approach; in other words, the gains are not coded directly, but rather the gain corrections that relate to an estimate of the gain noted. This estimate, in fact, corresponds to a level equalization factor between the λ (z) and A HF filters. (z) at the separation frequency between low band and high band (Fs / 4). The calculation of the factor (block 203) is detailed in FIG. 10 of the specification 3GPP TS 26.290, reproduced here in FIG. 2b. We will not detail here this figure. To summarize, the blocks
Mz) Mz)
210 a 213 servent a calculer l'énergie de la réponse impulsionnelle de -— - , 210 to 213 are used to calculate the energy of the impulse response of - -,
(\ - Q.9z-l)AHF (z) en se rappelant que le filtre AHF (z) modélise une bande haute repliée spectrale (à cause des propriétés spectrales du banc de filtre séparant les bandes basse et haute). Puisque les filtres sont interpolés par sous-trames, le gain n'est calculé qu'une fois par trame, et il est interpolé par sous-trames. (\ - Q.9z- l ) A HF (z), remembering that the HF filter (z) models a spectral folded high band (because of the spectral properties of the filter bank separating the low and high bands). Since the filters are interpolated by subframes, the gain is calculated only once per frame, and it is interpolated by subframes.
La technique de codage des gains d'extension de bande dans AMR-WB+, et plus précisément la compensation de niveaux des filtres LPC en leur point de jonction, est une méthode adaptée dans le contexte d'une extension de bande par modèles LPC en bande basse et haute, et on peut remarquer qu'une telle compensation de niveau entre filtres LPC n'est pas présente dans l'extension de bande du codée AMR-WB. Cependant, on peut vérifier dans la pratique que l'égalisation directe du niveau entre les deux filtres LPC à la fréquence de séparation n'est pas une méthode optimale et peut provoquer une surestimation d'énergie en bande-haute et des artefacts audibles dans certains cas ; on rappelle qu'un filtre LPC représente une enveloppe spectrale, ainsi le principe de l'égalisation du niveau entre deux filtres LPC pour une fréquence donnée revient à ajuster le niveau relatif de deux enveloppes LPC. Or un telle égalisation réalisée en une fréquence précise n'assure pas une complète continuité et cohérence globale de l'énergie (en fréquence) au voisinage du point d'égalisation lorsque l'enveloppe fréquentielle du signal fluctue de façon significative dans ce voisinage. Une façon mathématique de poser le problème consiste à remarquer que la continuité entre deux courbes peut être assurée en les forçant à se rejoindre en un même point, mais rien ne garantit que les propriétés locales (dérivées successives) coïncident de façon à assurer une cohérence plus globale. Le risque en assurant une continuité ponctuelle entre des enveloppes LPC bandes basse et haute est de fixer l'enveloppe de LPC en bande haute à un niveau relatif trop fort ou trop faible, le cas d'un niveau trop fort étant plus dommageable car il résulte en des artefacts plus gênants. Par ailleurs, la compensation de gain dans AMR-WB+ est avant tout une prédiction du gain connue au codeur et au décodeur et qui sert à réduire le débit nécessaire à la transmission d'information de gain mettant à l'échelle le signal d'excitation bande haute. Or, dans le contexte d'une amélioration du codage/décodage AMR-WB de façon interopérable, il n'est pas possible de modifier le codage existant des gains par sous-trames (0.8 kbit/s) de l'extension de bande dans le mode 23.85 kbit/s d'AMR-WB. De plus, pour les débits strictement inférieurs à 23.85 kbit/s, la compensation de niveaux de filtres LPC en bandes basse et haute peut être appliquée dans l'extension de bande d'un décodage compatible avec AMR-WB, cependant l'expérience montre que cette seule technique dérivée du codage AMR-WB+, appliquée sans optimisation, peut engendrer des problèmes de surestimation d'énergie de la bande haute (>6 kHz). The technique of coding bandwidth gains in AMR-WB +, and more specifically the level compensation of LPC filters at their junction point, is a suitable method in the context of band-based LPC bandwidth expansion. low and high, and it can be noticed that such a level compensation between LPC filters is not present in the band extension of the AMR-WB codec. However, it can be verified in practice that the direct equalization of the level between the two LPC filters at the separation frequency is not an optimal method and can cause an overestimation of high-band energy and audible artifacts in some case; It is recalled that an LPC filter represents a spectral envelope, so the principle of equalizing the level between two LPC filters for a given frequency amounts to adjusting the relative level of two LPC envelopes. However, such equalization performed at a precise frequency does not ensure complete continuity and overall coherence of the energy (in frequency) in the vicinity of the equalization point when the frequency envelope of the signal fluctuates significantly in this vicinity. A mathematical way of posing the problem consists in noting that the continuity between two curves can be assured by forcing them to meet at the same point, but there is no guarantee that the local properties (successive derivatives) coincide in order to ensure a more coherent overall. The risk in ensuring a one-off continuity between low and high band LPC envelopes is to set the LPC envelope in high band at a relative level that is too strong or too weak, the case of a too high level being more damaging because it results in more troublesome artifacts. Furthermore, the gain compensation in AMR-WB + is primarily a prediction of the known gain to the encoder and the decoder and which serves to reduce the bit rate necessary for the transmission of gain information scaling the excitation signal. high band. However, in the context of improved interoperable AMR-WB coding / decoding, it is not possible to modify the existing coding of subframe gains (0.8 kbit / s) of the band the 23.85 kbit / s mode of AMR-WB. In addition, for bit rates strictly below 23.85 kbit / s, high and low band LPC filter level compensation can be applied in the bandwidth of an AMR-WB compatible decoding, however the experiment shows that this technique derived from the AMR-WB + coding, applied without optimization, can cause problems of overestimation of energy of the high band (> 6 kHz).
Il existe donc un besoin pour améliorer la compensation de gains entre des filtres de prédiction linéaire de bande de fréquences différentes pour l'extension de bande de fréquence dans un codée de type AMR-WB ou une version interopérable de ce codée sans pour autant surestimer l'énergie dans une bande de fréquence et sans nécessiter d'informations supplémentaires du codeur. There is therefore a need to improve the gain compensation between linear frequency band prediction filters of different frequencies for the frequency band extension in an AMR-WB type code or an interoperable version of this codec without overestimating the energy in a frequency band and without requiring additional information from the coder.
La présente invention vient améliorer la situation.  The present invention improves the situation.
A cet effet, l'invention vise un procédé de détermination d'un facteur d'échelle optimisé à appliquer à un signal d'excitation ou à un filtre lors d'un procédé d'extension de bande de fréquence d'un signal audiofréquence, le procédé d'extension de bande comportant une étape de décodage ou d'extraction, dans une première bande de fréquence, d'un signal d'excitation et de paramètres de la première bande de fréquence comprenant des coefficients d'un filtre de prédiction linéaire, une étape de génération d'un signal d'excitation étendu sur au moins une deuxième bande de fréquence et une étape de filtrage par un filtre de prédiction linéaire pour la deuxième bande de fréquence. Le procédé de détermination est tel qu'il comporte les étapes suivantes:  For this purpose, the invention provides a method for determining an optimized scale factor to be applied to an excitation signal or a filter during a frequency band extension process of an audio frequency signal, the band extension method comprising a step of decoding or extracting, in a first frequency band, an excitation signal and parameters of the first frequency band comprising coefficients of a linear prediction filter a step of generating an extended excitation signal on at least a second frequency band and a filtering step by a linear prediction filter for the second frequency band. The determination method is such that it comprises the following steps:
- détermination d'un filtre de prédiction linéaire dit filtre additionnel, d'ordre inférieur au filtre de prédiction linéaire de la première bande de fréquence, les coefficients du filtre additionnel étant obtenus à partir des paramètres décodés ou extraits de la première bande de fréquence; et  determination of a linear prediction filter called an additional filter, of a lower order than the linear prediction filter of the first frequency band, the coefficients of the additional filter being obtained from the parameters decoded or extracted from the first frequency band; and
- calcul du facteur d'échelle optimisé en fonction au moins des coefficients du filtre additionnel.  calculating the optimized scale factor according to at least the coefficients of the additional filter.
Ainsi, l'utilisation d'un filtre additionnel d'ordre inférieur au filtre de la première bande de fréquence à égaliser, permet d'éviter les surestimations d'énergie dans les hautes fréquences qui pourraient résulter de fluctuations locales de l'enveloppe et qui peuvent perturber l'égalisation des filtres de prédiction.  Thus, the use of an additional filter of a lower order than the filter of the first frequency band to be equalized makes it possible to avoid overestimations of energy in the high frequencies which could result from local fluctuations of the envelope and which can disrupt the equalization of the prediction filters.
L'égalisation de gains entre les filtres de prédiction linéaire de la première et deuxième bande de fréquence est ainsi améliorée. Dans une application avantageuse du facteur d'échelle optimisé ainsi obtenu, le procédé d'extension de bande comprend une étape d'application du facteur d'échelle optimisé au signal d'excitation étendu. The equalization of gains between the linear prediction filters of the first and second frequency band is thus improved. In an advantageous application of the optimized scale factor thus obtained, the band extension method comprises a step of applying the optimized scaling factor to the extended excitation signal.
Dans un mode de réalisation adapté, l'application du facteur d'échelle optimisé est combinée à l'étape de filtrage dans la deuxième bande de fréquence.  In a suitable embodiment, the application of the optimized scaling factor is combined with the filtering step in the second frequency band.
Ainsi les étapes de filtrage et d'application du facteur d'échelle optimisé sont combinées à une seule étape de filtrage pour réduire la complexité de traitement.  Thus the steps of filtering and applying the optimized scaling factor are combined with a single filtering step to reduce the processing complexity.
Dans un mode de réalisation particulier, les coefficients du filtre additionnel sont obtenus par troncature de la fonction de transfert du filtre de prédiction linéaire de la première bande de fréquence pour obtenir un ordre inférieur.  In a particular embodiment, the coefficients of the additional filter are obtained by truncation of the transfer function of the linear prediction filter of the first frequency band to obtain a lower order.
Ce filtre additionnel d'ordre inférieur est donc obtenu de façon simple.  This additional low-order filter is thus obtained in a simple manner.
De plus, de façon à obtenir un filtre stable, les coefficients du filtre additionnel sont modifiés en fonction d'un critère de stabilité du filtre additionnel.  In addition, in order to obtain a stable filter, the coefficients of the additional filter are modified according to a criterion of stability of the additional filter.
Dans un mode de réalisation particulier, le calcul du facteur d'échelle optimisé comporte les étapes suivantes:  In a particular embodiment, the computation of the optimized scaling factor comprises the following steps:
- calcul des réponses en fréquence des filtres de prédiction linéaire des première et deuxième bandes de fréquence pour une fréquence commune;  calculating the frequency responses of the linear prediction filters of the first and second frequency bands for a common frequency;
- calcul de la réponse en fréquence du filtre additionnel pour cette fréquence commune; calculation of the frequency response of the additional filter for this common frequency;
- calcul du facteur d'échelle optimisé en fonction des réponses en fréquence ainsi calculées. calculation of the optimized scale factor according to the frequency responses thus calculated.
Ainsi, le facteur d'échelle optimisé est calculé de façon à éviter les artefacts gênants qui pourraient survenir dans le cas où la réponse en fréquence de filtre d'ordre supérieur de la première bande à proximité de la fréquence commune révélerait un pic ou une vallée du signal.  Thus, the optimized scaling factor is calculated to avoid annoying artifacts that might arise in the event that the higher order filter frequency response of the first band near the common frequency would reveal a peak or a valley of the signal.
Dans un mode de réalisation particulier, le procédé comprend en outre les étapes suivantes, mises en œuvre pour un débit de décodage prédéterminé:  In a particular embodiment, the method further comprises the following steps, implemented for a predetermined decoding rate:
- première mise à l'échelle du signal d'excitation étendu par un gain calculé par sous- trame fonction d'un rapport d'énergie entre le signal d'excitation décodé et le signal d'excitation étendu;  first scaling of the extended excitation signal by a calculated gain per subframe as a function of an energy ratio between the decoded excitation signal and the extended excitation signal;
-deuxième mise à l'échelle du signal d'excitation issu de la première mise à l'échelle par un gain de correction décodé;  second scaling of the excitation signal from the first scaling by a decoded correction gain;
- ajustement de l'énergie de l'excitation pour la sous-trame courante par un facteur d'ajustement calculé en fonction de l'énergie du signal obtenu après la deuxième mise à l'échelle et en fonction du signal obtenu après application du facteur d'échelle optimisé.  - adjustment of the excitation energy for the current subframe by an adjustment factor calculated as a function of the energy of the signal obtained after the second scaling and as a function of the signal obtained after application of the factor optimized scale.
Ainsi, des informations supplémentaires peuvent être utilisées pour améliorer la qualité du signal étendu pour un mode de fonctionnement prédéterminé. L'invention vise également un dispositif de détermination d'un facteur d'échelle optimisé à appliquer à un signal d'excitation ou à un filtre dans un dispositif d'extension de bande de fréquence d'un signal audiofréquence, le dispositif d'extension de bande comprenant un module de décodage ou d'extraction, dans une première bande de fréquence, d'un signal d'excitation et de paramètres de la première bande de fréquence comprenant des coefficients d'un filtre de prédiction linéaire, un module de génération d'un signal d'excitation étendu sur au moins une deuxième bande de fréquence et un module de filtrage par un filtre de prédiction linéaire pour la deuxième bande de fréquence. Le dispositif de détermination est tel qu'il comprend: Thus, additional information can be used to improve the quality of the extended signal for a predetermined mode of operation. The invention also relates to a device for determining an optimized scale factor to be applied to an excitation signal or to a filter in a frequency band extension device of an audiofrequency signal, the extension device band comprising a decoding or extraction module, in a first frequency band, of an excitation signal and parameters of the first frequency band comprising coefficients of a linear prediction filter, a generation module an excitation signal extended over at least a second frequency band and a filtering module by a linear prediction filter for the second frequency band. The determination device is such that it comprises:
- un module de détermination d'un filtre de prédiction linéaire dit filtre additionnel, d'ordre inférieur au filtre de prédiction linéaire de la première bande de fréquence, les coefficients du filtre additionnel étant obtenus à partir des paramètres décodés ou extraits de la première bande de fréquence; et  a module for determining a linear prediction filter called an additional filter, of a lower order than the linear prediction filter of the first frequency band, the coefficients of the additional filter being obtained from the parameters decoded or extracted from the first band; frequency; and
- un module de calcul du facteur d'échelle optimisé en fonction au moins des coefficients du filtre additionnel.  a module for calculating the scaling factor optimized according to at least the coefficients of the additional filter.
L'invention vise un décodeur comportant un dispositif tel que décrit.  The invention relates to a decoder comprising a device as described.
Elle vise un programme informatique comportant des instructions de code pour la mise en œuvre des étapes du procédé de détermination d'un facteur d'échelle optimisé tel que décrit, lorsque ces instructions sont exécutées par un processeur.  It relates to a computer program comprising code instructions for implementing the steps of the method of determining an optimized scale factor as described, when these instructions are executed by a processor.
Enfin l'invention se rapporte à un support de stockage, lisible par un processeur, intégré ou non au dispositif de détermination d'un facteur d'échelle optimisé, éventuellement amovible, mémorisant un programme informatique mettant en œuvre un procédé de détermination d'un facteur d'échelle optimisé tel que décrit précédemment.  Finally, the invention relates to a storage medium, readable by a processor, integrated or not to the device for determining an optimized scale factor, possibly removable, storing a computer program implementing a method of determining a optimized scaling factor as previously described.
D'autres caractéristiques et avantages de l'invention apparaîtront plus clairement à la lecture de la description suivante, donnée uniquement à titre d'exemple non limitatif, et faite en référence aux dessins annexés, sur lesquels :  Other features and advantages of the invention will appear more clearly on reading the following description, given solely by way of nonlimiting example, and with reference to the appended drawings, in which:
la figure 1 illustre une partie d'un décodeur de type AMR-WB mettant en œuvre des étapes d'extension de bande de fréquence de l'état de l'art et tel que décrit précédemment;  FIG. 1 illustrates a part of an AMR-WB decoder implementing frequency band extension steps of the state of the art and as previously described;
- les figures 2a et 2b présentent le codage de la bande haute dans le codée AMR- FIGS. 2a and 2b show the coding of the high band in the AMR code.
WB+ selon l'état de l'art et tel que décrit précédemment; WB + according to the state of the art and as described above;
la figure 3 illustre un décodeur interopérable avec le codage AMR-WB et intégrant un dispositif d'extension de bande utilisé selon un mode de réalisation de l'invention ;  FIG. 3 illustrates an interoperable decoder with the AMR-WB coding and integrating a band extension device used according to one embodiment of the invention;
- la figure 4 illustre un dispositif de détermination d'un facteur d'échelle optimisé par sous-trame en fonction du débit, selon un mode de réalisation de l'invention; et les figures 5a et 5b illustrent les réponses en fréquences des filtres utilisées pour le calcul du facteur d'échelle optimisé selon un mode de réalisation de l'invention; la figure 6 illustre sous forme d'organigramme, les étapes principales d'un procédé de détermination d'un facteur d'échelle optimisé selon un mode de réalisation de l'invention; FIG. 4 illustrates a device for determining a subframe-optimized scale factor as a function of the flow rate, according to one embodiment of the invention; and FIGS. 5a and 5b illustrate the frequency responses of the filters used for calculating the optimized scaling factor according to one embodiment of the invention; FIG. 6 illustrates in flowchart form the main steps of a method for determining an optimized scale factor according to an embodiment of the invention;
la figure 7 illustre un mode de réalisation dans le domaine fréquentiel d'un dispositif de détermination de facteur d'échelle optimisé lors d'une extension de bande;  FIG. 7 illustrates a frequency domain embodiment of an optimized scale factor determination device during a band extension;
la figure 8 illustre une réalisation matérielle d'un dispositif de détermination de facteur d'échelle optimisé lors d'une extension de bande selon l'invention.  FIG. 8 illustrates a hardware embodiment of an optimized scale factor determination device during a band extension according to the invention.
La figure 3 illustre un exemple de décodeur, compatible avec la norme AMR- WB/G.722.2 dans lequel on retrouve une extension de bande comprenant une détermination d'un facteur d'échelle optimisé selon un mode de réalisation du procédé de l'invention, mis en œuvre par le dispositif d'extension de bande illustré par le bloc 309. FIG. 3 illustrates an exemplary decoder, compatible with the AMR-WB / G.722.2 standard, in which there is a band extension comprising a determination of an optimized scale factor according to one embodiment of the method of the invention. implemented by the tape extender illustrated by block 309.
Contrairement au décodage AMR-WB qui fonctionne avec une fréquence d'échantillonnage de sortie de 16 kHz, on considère ici un décodeur qui peut fonctionner avec un signal de sortie (synthèse) à la fréquence fs = 8, 16, 32 ou 48 kHz. A noter qu'on suppose ici que le codage a été effectué selon l'algorithme AMR-WB avec une fréquence interne de 12.8 kHz pour le codage CELP en bande basse et à 23.85 kbit/s un codage de gain par sous- trame à la fréquence de 16 kHz ; même si l'invention est décrite ici au niveau du décodage, on suppose ici que le codage peut aussi fonctionner avec un signal d'entrée à la fréquence fs = 8, 16, 32 ou 48 kHz et des opérations de ré-échantillonnage adéquate, dépassant le cadre de l'invention, sont mises en œuvre au codage en fonction de la valeur de fs. On peut noter que quand fs=8 kHz, dans le cas d'un décodage compatible avec AMR-WB, il n'est pas nécessaire d'étendre la bande basse 0-6.4 kHz, car la bande audio reconstruite à la fréquence fs est limitée à 0-4000 Hz.  Unlike AMR-WB decoding which operates with an output sampling frequency of 16 kHz, we consider here a decoder that can operate with an output signal (synthesis) at the frequency fs = 8, 16, 32 or 48 kHz. Note that it is assumed here that the coding was performed according to the AMR-WB algorithm with an internal frequency of 12.8 kHz for CELP coding in the low band and at 23.85 kbit / s a sub-frame gain coding at the same time. frequency of 16 kHz; even if the invention is described here at the decoding level, it is assumed here that the coding can also operate with an input signal at the frequency fs = 8, 16, 32 or 48 kHz and adequate resampling operations, beyond the scope of the invention, are implemented at the coding according to the value of fs. It can be noted that when fs = 8 kHz, in the case of a decoding compatible with AMR-WB, it is not necessary to extend the low band 0-6.4 kHz, because the audio band reconstructed at the frequency fs is limited to 0-4000 Hz.
A la figure 3, le décodage CELP (BF pour basses fréquences) fonctionne toujours à la fréquence interne de 12.8 kHz, comme dans AMR-WB, et l'extension de bande (HF pour hautes fréquences) utilisée pour l'invention fonctionne à la fréquence de 16 kHz, les synthèses BF et HF sont combinées (bloc 312) à la fréquence fs après ré-échantillonnage adéquat (bloc 306 et traitement interne au bloc 311). Dans des variantes de réalisation, la combinaison des bandes basse et haute pourra se faire à 16 kHz, après avoir rééchantillonnée la bande basse de 12.8 à 16 kHz, avant de ré-échantillonner le signal combiné à la fréquence fs.  In Figure 3, CELP decoding (BF for low frequencies) always operates at the internal frequency of 12.8 kHz, as in AMR-WB, and the band extension (HF for high frequencies) used for the invention operates at the same time. frequency of 16 kHz, the synthesis BF and HF are combined (block 312) at the frequency fs after adequate resampling (block 306 and internal processing block 311). In alternative embodiments, the combination of the low and high bands can be done at 16 kHz, after resampling the low band of 12.8 to 16 kHz, before resampling the combined signal at the frequency fs.
Le décodage selon la figure 3 dépend du mode (ou débit) AMR-WB associé à la trame courante reçue. A titre indicatif et sans que cela impacte le bloc 309, le décodage de la partie CELP en bande basse comporte les étapes suivantes: • Démultiplexage des paramètres codés (bloc 300) en cas de trame correctement reçue ( =0 où est le « bad frame indicator » valant 0 pour une trame reçue et 1 pour une trame perdue) The decoding according to FIG. 3 depends on the mode (or bit rate) AMR-WB associated with the current frame received. As an indication and without this having an effect on block 309, the decoding of the low band CELP part comprises the following steps: • Demultiplexing the coded parameters (block 300) in the case of a correctly received frame (= 0 where is the "bad frame indicator" worth 0 for a received frame and 1 for a lost frame)
• Décodage des paramètres ISF avec interpolation et conversion en coefficients LPC (bloc 301) comme décrit dans la clause 6.1 de la norme G.722.2.  • Decoding of ISF parameters with interpolation and conversion to LPC coefficients (block 301) as described in clause 6.1 of G.722.2.
• Décodage de l'excitation CELP (bloc 302), avec une partie adaptative et fixe pour reconstruire l'excitation (exc ou u \n) ) dans chaque sous-trame de longueur 64 à 12.8 kHz:  • Decoding of the CELP excitation (block 302), with an adaptive and fixed part to reconstruct the excitation (exc or u \ n)) in each subframe of length 64 to 12.8 kHz:
u \n) = gpv(n) + gcc(n) , « = (), · · · , 63 u \ n) = g p v (n) + g c c (n) "= (), · · ·, 63
en suivant les notations de la clause 7.1.2.1 de la recommandation ITU-T G.718 d'un décodeur interopérable avec le codeur/décodeur AMR-WB, concernant le décodage CELP, où v(n) et c(n) sont respectivement les mots de code des dictionnaires adaptatif et fixe, et g et gc sont les gains décodés associés. Cette excitation u \n) est utilisée dans le dictionnaire adaptatif de la sous-trame suivante ; elle est ensuite post-traitée et on distingue comme dans G.718 l'excitation u \n) (aussi notée exc) de sa version post-traitée modifiée u(n) (aussi notée exc2) qui sert d'entrée au filtre de synthèse, 1 / Â(z) , dans le bloc 303. by following the notation of clause 7.1.2.1 of ITU-T Recommendation G.718 of an interoperable decoder with the AMR-WB encoder / decoder, concerning CELP decoding, where v (n) and c (n) are respectively the code words of the adaptive and fixed dictionaries, and g and g c are the associated decoded gains. This excitation u \ n) is used in the adaptive dictionary of the following subframe; it is then post-processed and one distinguishes as in G.718 the excitation u \ n) (also noted exc) of its modified post-processed version u (n) (also noted exc2) which serves as input to the filter of synthesis, 1 / Â (z), in block 303.
• Filtrage de synthèse par 1 / Â(z) (bloc 303) où le filtre LPC décodé Â(z) est d'ordre 16  • Synthetic filtering by 1 /  (z) (block 303) where the decoded LPC filter  (z) is of order 16
• Post-traitement bande étroite (bloc 304) selon la clause 7.3 de G.718 si fs=8 kHz. • Narrowband post-processing (block 304) according to clause 7.3 of G.718 if fs = 8 kHz.
• Désaccentuation (bloc 305) par le filtre 1 / (l - 0.68z_1 ) • Deactivation (block 305) by the filter 1 / (l - 0.68z _1 )
• Post-traitement des basses fréquences (dit « bass posfilter ») (bloc 306) atténuant le bruit inter-harmonique en basses fréquences tel que décrit à la clause 7.14.1.1 de G.718. Ce traitement introduit un retard qui est pris en compte dans le décodage de la bande haute (>6.4 kHz).  • Post-processing of low frequencies (called "bass posfilter") (block 306) attenuating the inter-harmonic noise at low frequencies as described in clause 7.14.1.1 of G.718. This processing introduces a delay which is taken into account in the decoding of the high band (> 6.4 kHz).
• Ré-échantillonnage de la fréquence interne de 12.8 kHz à la fréquence de sortie fs (bloc 307). Plusieurs réalisations sont possibles. Sans perte de généralité, on considère ici à titre d'exemple que si fs=8 ou 16 kHz, le ré-échantillonnage décrit dans la clause 7.6 de G.718 est repris ici, et si fs=32 ou 48 kHz, des filtres à réponse impulsionnelle finie (FIR) supplémentaires sont utilisés.  • Resampling of the internal frequency from 12.8 kHz to the output frequency fs (block 307). Several achievements are possible. Without loss of generality, we consider here as an example that if fs = 8 or 16 kHz, the resampling described in clause 7.6 of G.718 is repeated here, and if fs = 32 or 48 kHz, filters Finite Impulse Response (FIR) are used.
• Calcul des paramètres du "noise gâte" (bloc 308) qui est réalisé de façon préférentielle comme décrit dans la clause 7.14.3 de G.718 pour « améliorer » la qualité des silences par réduction du niveau.  • Calculation of the parameters of the "noise gate" (block 308) which is carried out preferentially as described in clause 7.14.3 of G.718 to "improve" the quality of the silences by reducing the level.
Dans des variantes qui peuvent être mises en œuvre pour l'invention, les post-traitements appliqués à l'excitation peuvent être modifiés (par exemple, la dispersion de phase peut être améliorée) ou ces post-traitements peuvent être étendus (par exemple, une réduction du bruit inter-harmonique peut être mise en œuvre), sans affecter la nature de l'extension de bande. On peut noter que l'utilisation des blocs 306, 308, 314 est optionnelle. In variants that can be implemented for the invention, the post-treatments applied to the excitation can be modified (for example, the phase dispersion can be improved) or these post-treatments can be extended (for example, inter-harmonic noise reduction can be implemented), without affecting the nature of the band extension. It can be noted that the use of blocks 306, 308, 314 is optional.
On notera également que le décodage de la bande basse décrit ci-dessus suppose une trame courante dite « active » avec un débit entre 6.6 et 23.85 kbit/s. En fait, quand le mode DTX (transmission continue en français) est activé, certaines trames peuvent être codées comme « inactives » et dans ce cas on peut soit transmettre un descripteur de silence (sur 35 bits) soit ne rien transmettre. En particulier, on rappelle que la trame SID décrit plusieurs paramètres : paramètres ISF moyennés sur 8 trames, énergie moyenne sur 8 trames, flag de "dithering" pour la reconstruction de bruit non stationnaire. Dans tous les cas, au décodeur, on retrouve le même modèle de décodage que pour une trame active, avec une reconstruction de l'excitation et d'un filtre LPC pour la trame courante, ce qui permet d'appliquer l'extension de bande même sur des trames inactives. Le même constat s'applique pour le décodage de « trames perdues » (ou FEC, PLC) dans lequel le modèle LPC est appliqué.  Note also that the decoding of the low band described above assumes a current frame called "active" with a rate between 6.6 and 23.85 kbit / s. In fact, when the DTX (Continuous Transmission in French) mode is activated, some frames can be coded as "inactive" and in this case you can either transmit a silence descriptor (on 35 bits) or not transmit anything. In particular, it is recalled that the SID frame describes several parameters: ISF parameters averaged over 8 frames, average energy over 8 frames, "dithering" flag for the non-stationary noise reconstruction. In all cases, at the decoder, we find the same decoding model as for an active frame, with a reconstruction of the excitation and an LPC filter for the current frame, which makes it possible to apply the band extension. even on inactive frames. The same applies for the decoding of "lost frames" (or FEC, PLC) in which the LPC model is applied.
Dans le mode de réalisation décrit ici et en référence à la figure 7, , le décodeur permet d'étendre la bande basse décodée (50-6400 Hz en tenant en compte du filtrage passe-haut à 50 Hz au décodeur, 0-6400 Hz dans le cas général) à une bande étendue dont la largeur varie, allant approximativement de 50-6900 Hz à 50-7700 Hz en fonction du mode mis en œuvre dans la trame courante. On peut ainsi parler d'une première bande de fréquence de 0 à 6400Hz et d'une deuxième bande de fréquence de 6400 à 8000Hz. En réalité, dans le mode de réalisation privilégié, l'extension de l'excitation est réalisée dans le domaine fréquentiel dans une bande de 5000 à 8000 Hz, pour permettre un filtrage passe- bande de largeur 6000 à 6900 ou 7700 Hz.  In the embodiment described here and with reference to FIG. 7, the decoder makes it possible to extend the decoded low band (50-6400 Hz while taking into account 50 Hz high-pass filtering at the decoder, 0-6400 Hz in the general case) to an extended band whose width varies, ranging from approximately 50-6900 Hz to 50-7700 Hz depending on the mode implemented in the current frame. We can talk about a first frequency band from 0 to 6400Hz and a second frequency band from 6400 to 8000Hz. In reality, in the preferred embodiment, the extension of the excitation is carried out in the frequency domain in a band of 5000 to 8000 Hz, to allow bandpass filtering of 6000 to 6900 or 7700 Hz width.
A 23.85 kbit/s, l'information de correction de gain HF (0.8 kbit/s) transmise à 23.85 kbit/s est ici décodée. Son utilisation est détaillée plus loin, en référence à la figure 4. La partie synthèse bande haute est réalisée dans le bloc 309 représentant le dispositif d'extension de bande utilisé pour l'invention et qui est détaillé à la figure 7 dans un mode de réalisation.  At 23.85 kbit / s, the HF gain correction information (0.8 kbit / s) transmitted at 23.85 kbit / s is decoded here. Its use is detailed below, with reference to FIG. 4. The high band synthesis part is realized in block 309 representing the band extension device used for the invention and which is detailed in FIG. production.
Afin d'aligner les bandes basses et hautes décodées, un retard (bloc 310) est introduit pour synchroniser les sorties des blocs 306 et 307 et la bande haute synthétisée à 16 kHz est ré-échantillonnée de 16 kHz à la fréquence fs (sortie de bloc 311). La valeur du retard T dépend de la façon de synthétiser le signal bande haute, de la fréquence fs de même que du post-traitement des basses fréquences. Ainsi, de façon générale la valeur de T dans le bloc 310 devra être ajustée en fonction de la mise en œuvre spécifique. Les bandes basse et haute sont ensuite combinées (ajoutées) dans le bloc 312 et la synthèse obtenue est post-traitée par filtrage passe-haut à 50 Hz (de type IIR) d'ordre 2 dont les coefficients dépendent de la fréquence fs (bloc 313) et post-traitement de sortie avec application optionnelle du "noise gate"de façon similaire à G.718 (bloc 314). In order to align the decoded low and high bands, a delay (block 310) is introduced to synchronize the outputs of the blocks 306 and 307 and the synthesized high band at 16 kHz is resampled from 16 kHz to the frequency fs (output of block 311). The value of the delay T depends on how to synthesize the high band signal, the frequency fs as well as the post-processing of the low frequencies. Thus, in general, the value of T in the block 310 will have to be adjusted according to the specific implementation. The low and high bands are then combined (added) in block 312 and the resulting synthesis is post-processed by high-order 50 Hz (type IIR) high-pass filtering whose coefficients depend on the frequency fs (block 313) and output post-processing with optional noise gate application similar to G.718 (block 314).
En référence à la figure 3, on décrit à présent un mode de réalisation d'un dispositif de détermination d'un facteur d'échelle optimisé à appliquer à un signal d'excitation lors d'un processus d'extension de bande de fréquence. Ce dispositif est inclus dans le bloc d'extension de bande 309 décrit précédemment.  Referring to FIG. 3, an embodiment of an optimized scale factor determining device to be applied to an excitation signal during a frequency band extension process is now described. This device is included in the band extension block 309 previously described.
Ainsi, le bloc 400, à partir d'un signal d'excitation décodé dans une première bande de fréquence u(n) , effectue une extension de bande pour obtenir un signal d'excitation étendu uHB (n) sur au moins une deuxième bande de fréquence. Thus, the block 400, from a decoded excitation signal in a first frequency band u (n), performs a band extension to obtain an extended excitation signal u HB (n) on at least a second frequency band.
On notera ici, que l'estimation de facteur d'échelle optimisé selon l'invention est indépendante de la façon d'obtenir le signal uHB (n) . Une condition concernant son énergie est cependant importante. En effet, il faut que l'énergie de la bande haute de 6000 à 8000 Hz soit à un niveau similaire à l'énergie de la bande de 4000 à 6000 Hz du signal d'excitation décodé à la sortie du bloc 302. De plus, puisque le signal bande basse est désaccentué (bloc 305), il faut aussi appliquer la désaccentuation au signal d'excitation bande haute, soit en utilisant un filtre de désaccentuation propre, soit en multipliant par un facteur constant qui correspond à une atténuation moyenne du filtre mentionné. Cette condition ne s'applique pas au cas du débit 23.85 kbit/s qui utilise les informations supplémentaires transmises par le codeur. Dans ce cas l'énergie du signal d'excitation bande haute doit être cohérente avec l'énergie du signal correspondant au codeur, comme expliqué plus tard. It will be noted here that the optimized scale factor estimation according to the invention is independent of how to obtain the signal u HB (n). A condition regarding its energy is important, however. Indeed, it is necessary that the energy of the high band of 6000 to 8000 Hz is at a level similar to the energy of the band 4000 to 6000 Hz of the decoded excitation signal at the output of the block 302. since the low-band signal is de-emphasized (block 305), it is also necessary to apply the de-emphasis to the high-band excitation signal, either by using an own de-emphasis filter or by multiplying by a constant factor which corresponds to a mean attenuation of mentioned filter. This condition does not apply to the 23.85 kbit / s rate that uses the additional information transmitted by the encoder. In this case, the energy of the high band excitation signal must be consistent with the signal energy corresponding to the encoder, as explained later.
L'extension de bande de fréquence peut par exemple être mise en œuvre de la même façon que pour le décodeur de type AMR-WB décrit en référence à la figure 1 dans les blocs 100 à 102, à partir d'un bruit blanc.  The frequency band extension may for example be implemented in the same way as for the AMR-WB type decoder described with reference to FIG. 1 in the blocks 100 to 102, from a white noise.
Dans un autre mode de réalisation, cette extension de bande peut s'effectuer à partir d'une combinaison d'un bruit blanc et d'un signal d'excitation décodé comme illustré et décrit ultérieurement pour les blocs 700 à 707 de la figure 7.  In another embodiment, this band extension can be performed from a combination of a white noise and a decoded excitation signal as illustrated and subsequently described for blocks 700 to 707 of FIG. .
D'autres méthodes d'extension de bande de fréquence avec conservation du niveau d'énergie entre le signal d'excitation décodé et le signal d'excitation étendu comme décrit ci- dessous, peuvent bien sur être envisagées pour le bloc 400.  Other frequency band extension methods with conservation of the energy level between the decoded excitation signal and the extended excitation signal as described below can of course be envisaged for block 400.
De plus, le module d'extension de bande peut également être indépendant du décodeur et peut effectuer une extension de bande d'un signal audio existant stocké ou transmis au module d'extension, avec une analyse du signal audio pour en extraire une excitation et un filtre LPC. Dans ce cas, le signal d'excitation en entrée du module d'extension, n'est plus un signal décodé mais un signal extrait après analyse, de même que les coefficients du filtre de prédiction linéaire de la première bande de fréquence utilisés dans le procédé de détermination du facteur d'échelle optimisé dans une mise en œuvre de l'invention. In addition, the tape expansion module may also be independent of the decoder and may extend a band of an existing audio signal stored or transmitted to the extension module, with an analysis of the audio signal to extract an excitation and an LPC filter. In this case, the excitation signal at the input of the extension module is no longer a decoded signal but a signal extracted after analysis, as are the coefficients of the linear prediction filter of the first frequency band used in the method of determining the optimized scale factor in an implementation of the invention.
Dans l'exemple illustré à la figure 4, on considère tout d'abord le cas des débits <23.85kbit/s, pour lequel la détermination du facteur d'échelle optimisé se limite au bloc 401. On calcule dans ce cas un facteur d'échelle optimisé, noté gHB2 (m) . Dans un mode de réalisation, ce calcul est effectué de façon préférentielle par sous-trame et il consiste à égaliser les niveaux des réponses en fréquences des filtres LPC 1 / Â(z) et 1 / Λ(ζ/ γ) utilisés en basses et hautes fréquences, comme décrit ultérieurement en référence à la figure 7, avec des précautions supplémentaires pour éviter les cas de surestimations qui peuvent résulter en une énergie trop grande de la bande haute synthétisée et générer donc des artefacts audibles. In the example illustrated in FIG. 4, we first consider the case of rates <23.85 kbit / s, for which the determination of the optimized scaling factor is limited to block 401. In this case, a factor of optimized scale, denoted g HB2 (m). In one embodiment, this calculation is preferably carried out by subframe and it consists in equalizing the levels of the frequency responses of the LPC filters 1 / Â (z) and 1 / Λ (ζ / γ) used at low frequencies and high frequencies, as described later with reference to Figure 7, with additional precautions to avoid the cases of overestimations that can result in too much energy of the synthesized high band and thus generate audible artifacts.
Dans un mode de réalisation alternatif, on pourra garder le filtre de synthèse HF extrapolé 1 / Aex'(z/ γ) tel que mis en œuvre dans le décodeur AMR-WB ou un décodeur interopérable avec le codeur/décodeur AMR-WB, par exemple selon la recommandation ITU-T G.718, au lieu du filtre 1 / A(z/ γ) . La compensation selon l'invention s'effectue alors à partir des filtresIn an alternative embodiment, it will be possible to keep the extrapolated HF synthesis filter 1 / A ex '(z / γ) as implemented in the AMR-WB decoder or a decoder interoperable with the AMR-WB encoder / decoder, for example according to ITU-T Recommendation G.718, instead of filter 1 / A (z / γ). The compensation according to the invention is then carried out from the filters
1 / Â(z) et I I Aext(zl γ) . 1 / (z) and IIA ext (zl γ).
La détermination du facteur d'échelle optimisé s'effectue aussi par la détermination (en 401a) d'un filtre de prédiction linéaire dit filtre additionnel, d'ordre inférieur au filtre de prédiction linéaire de la première bande de fréquence 1 / Mz) , les coefficients du filtre additionnel étant obtenus à partir des paramètres décodés ou extraits de la première bande de fréquence. Le facteur d'échelle optimisé est ensuite calculé (en 401b) en fonction au moins de ces coefficients pour être appliqué au signal d'excitation étendu uHB (n) . The determination of the optimized scale factor is also performed by the determination (in 401a) of a linear prediction filter called additional filter, of a lower order than the linear prediction filter of the first frequency band 1 / M z ) , the coefficients of the additional filter being obtained from the parameters decoded or extracted from the first frequency band. The optimized scaling factor is then calculated (at 401b) based on at least one of these coefficients to be applied to the extended excitation signal u HB (n).
Le principe de la détermination du facteur d'échelle optimisé, mis en œuvre dans le bloc 401 est illustré aux figures 5a et 5b avec des exemples concrets obtenus à partir de signaux échantillonnés à 16 kHz ; les valeurs d'amplitude de réponses en fréquence, notées plus loin R, P, Q, de 3 filtres sont calculées à la fréquence commune de 6000 Hz (ligne verticale en pointillés) dans la sous-trame courante, dont l'indice m n'est pas rappelé ici dans la notations des filtres LPC interpolés par sous-trame pour alléger le texte. La valeur de 6000 Hz est choisie de sorte qu'elle soit proche de la fréquence de Nyquist de la bande basse, soit 6400 Hz. Il est préférable de ne pas prendre cette fréquence de Nyquist pour déterminer le facteur d'échelle optimisé. En effet l'énergie du signal décodé en basses fréquences est typiquement déjà atténuée à 6400 Hz. De plus, l'extension de bande décrite ici est effectuée sur une deuxième bande de fréquence dite bande haute qui va de 6000 à 8000 Hz. A noter que dans des variantes de l'invention, une autre fréquence que 6000 Hz pourra être choisie, sans perte de généralité pour déterminer le facteur d'échelle optimisé. On pourra aussi considérer le cas où les deux filtres LPC sont définis pour les bandes séparés (comme dans AMR-WB+). Dans ce cas R, P et Q seront calculées à la fréquence de séparation. The principle of the determination of the optimized scaling factor implemented in block 401 is illustrated in FIGS. 5a and 5b with concrete examples obtained from signals sampled at 16 kHz; the amplitude values of frequency responses, noted later R, P, Q, of 3 filters are calculated at the common frequency of 6000 Hz (vertical dashed line) in the current subframe, whose index m n This is not recalled here in the notation of LPC filters interpolated by sub-frame to lighten the text. The value of 6000 Hz is chosen so that it is close to the Nyquist frequency of the low band, ie 6400 Hz. It is preferable not to take this Nyquist frequency to determine the optimized scale factor. Indeed, the energy of the signal decoded at low frequencies is typically already attenuated at 6400 Hz. In addition, the band extension described here is performed on a second so-called high band frequency band which ranges from 6000 to 8000 Hz. that in variants of the invention, another frequency that 6000 Hz can be chosen, without loss of generality to determine the optimized scale factor. We can also consider the case where the two LPC filters are defined for separate bands (as in AMR-WB +). In this case R, P and Q will be calculated at the separation frequency.
Les figures 5a et 5b illustrent comment sont définies les quantités R, P, Q. Figures 5a and 5b illustrate how the quantities R, P, Q are defined.
La première étape consiste à calculer les réponses en fréquences R et P respectivement du filtre de prédiction linéaire de la première bande de fréquence (bande basse) et de la deuxième bande de fréquence (bande haute) à la fréquence de 6000 Hz. On calcule d'abord : The first step is to calculate the R and P frequency responses respectively of the linear prediction filter of the first frequency band (low band) and the second frequency band (high band) at the frequency of 6000 Hz. 'on board :
où M = 16 est l'ordre du filtre LPC décodé 1 / Â(z) , et ^correspond à la fréquence de 6000 Hz normalisée pour la fréquence d'échantillonnage de 12.8 kHz, soit : where M = 16 is the order of the decoded LPC filter 1 / Â (z), and ^ corresponds to the normalized 6000 Hz frequency for the sampling frequency of 12.8 kHz, ie:
12800 12800
Ensuite, de façon similaire, on calcule : Then, similarly, we calculate:
ou or
6000  6000
θ' = 2π θ '= 2π
16000  16000
Dans un mode de réalisation privilégié, les quantités P et R sont calculées selon le pseudocode suivant:  In a preferred embodiment, the quantities P and R are calculated according to the following pseudocode:
px = py = 0  px = py = 0
rx = ry = 0  rx = ry = 0
fori=0 to 16  fori = 0 to 16
px = px + Ap[i]*exp_tab_p[i]  px = px + Ap [i] * exp_tab_p [i]
py = py + Ap[i]*exp_tab_p[33-i]  py = py + Ap [i] * exp_tab_p [33-i]
rx = rx + Aq[i]*exp_tab_q[i]  rx = rx + Aq [i] * exp_tab_q [i]
ry = ry + Aq[i]*exp_tab_q[33-i]  ry = ry + Aq [i] * exp_tab_q [33-i]
end for  end for
P = l/sqrt(px*px+py*py)  P = 1 / sqrt (px * px + py * py)
R = l/sqrt(rx*rx+ry*ry)  R = 1 / sqrt (rx * rx + ry * ry)
où Aq[i]= â; correspond aux coefficients de Â(z) (d'ordre 16), Ap[i]= ylât correspond aux coefficient de Λ(ζ/ γ) , sqrt() correspond à l'opération de racine carrée et les tableaux exp_tab_p et exp_tab_q de taille 34 contiennent les parties réelles et imaginaires des exponentielles complexes associée à la fréquence de 6000 Hz, avec where Aq [i] = â ; corresponds to the coefficients of  (z) (of order 16), Rv [i] = y L to T corresponds to the coefficient of Λ (ζ / γ), sqrt () is the square root operation and the tables exp_tab_p and exp_tab_q of size 34 contain the real and imaginary parts of complex exponentials associated with the frequency of 6000 Hz, with
Le filtre de prédiction additionnel est obtenu par exemple en tronquant de façon adéquate le polynôme Â(z) à l'ordre 2. The additional prediction filter is obtained for example by appropriately truncating the polynomial  (z) to order 2.
En fait la troncature directe à l'ordre conduit au filtre 1 + <31 + â2 , ce qui peut poser problème car rien ne garantit en général que ce filtre d'ordre 2 est stable. Dans un mode de réalisation privilégiée, on détecte donc la stabilité du filtre 1 + <31 + â2 et on utilise un filtre 1 + <¾ '+ <¾ ' , dont les coefficients sont tirés de 1 + en fonction de la détection d'instabilité. Plus précisément, on initialise : In fact direct truncation to the order leads to the filter 1 + <3 1 + â 2 , which can be problematic because there is no guarantee in general that this filter of order 2 is stable. In a preferred embodiment, the stability of the filter 1 + <3 1 + to 2 is thus detected and a filter 1 + <¾ '+ <¾' is used whose coefficients are derived from 1 + as a function of the detection. instability. More precisely, we initialize:
at = <¾ , i=l, 2 a t = <¾, i = 1, 2
La stabilité du filtre 1 + peut être vérifiée de différente façon, on utilise ici une conversion dans le domaine des coefficients PARCOR (ou coefficients de réflexion) en calculant :  The stability of the filter 1 + can be verified in different ways, here a conversion is used in the domain of the PARCOR coefficients (or reflection coefficients) by calculating:
kx = <¾ '/ (1 + â2 ') k x = <¾ '/ (1 + to 2 ')
La stabilité est vérifiée si < 1 , i=l,2. On modifie donc de façon conditionnelle la valeur de gavant d'assurer la stabilité du filtre, avec les étapes suivantes : Stability is checked if <1, i = 1, 2. The value of gavant is then conditionally modified to ensure the stability of the filter, with the following steps:
Γ min(0.6,fc ) L > 0  Γ min (0.6, fc) L> 0
k2 ^ \ 2 2 k 2 ^ \ 2 2
max(-0.6,fc2) k2 < 0max (-0.6, fc 2 ) k 2 <0
1 ^ [max(-0.99, ^2) ^ < 0 où min(.,.) et max(.,.) donnent respectivement le minimum et le maximum de 2 opérandes. On note que les valeurs de seuils, 0.99 pour ^ et 0.6 pour k2 , pourront être ajustées dans des variantes de l'invention. On rappelle que le premier coefficient de réflexion, , caractérise la pente spectrale (ou tilt) du signal modélisé à l'ordre 1 ; dans l'invention on sature la valeur de kx à une valeur proche de la limite de stabilité, afin de préserver cette pente et conserver un tilt similaire à celui de 1/ A(z) . On rappelle aussi que le second coefficient de réflexion, k2 , caractérise le niveau de résonance du modèle de signal à l'ordre 2 ; puisque l'utilisation d'un filtre d'ordre 2 vise à éliminer l'influence de telles résonances autour de la fréquence de 6000 Hz, on limite plus fortement la valeur de k2 , cette limite est fixée à 0.6. 1 ^ [max (-0.99, ^ 2 ) ^ <0 where min (.,.) And max (.,.) Respectively give the minimum and maximum of 2 operands. Note that the threshold values, 0.99 for ^ and 0.6 for k 2 , may be adjusted in variants of the invention. It is recalled that the first reflection coefficient,, characterizes the spectral slope (or tilt) of the signal modeled at the order 1; in the invention saturates the value of k x to a value close to the stability limit, in order to preserve this slope and maintain a tilt similar to that of 1 / A (z). It is also recalled that the second reflection coefficient, k 2 , characterizes the resonance level of the signal model at order 2; since the use of a filter of order 2 aims at eliminating the influence of such resonances around the frequency of 6000 Hz, the value of k 2 is more strongly limited, this limit is fixed at 0.6.
Les coefficients de 1 + âj '+ â2 ' sont alors obtenus par: 1 coefficients j + â 'â + 2' are then obtained by:
<¾ ' = (l + fcj)^ On calcule donc finalement la réponse en fréquence du filtre additionnel: „ „ 6000 <¾ '= (l + fc j ) ^ So we finally calculate the frequency response of the additional filter: "" 6000
avec Θ = TÎ . Cette quantité est calculée de façon préférentielle selon le pseudo-code with Θ = T. This quantity is preferably calculated according to the pseudo-code
12800  12800
suivant : next :
qx = qy = 0  qx = qy = 0
fori=0 to 2  fori = 0 to 2
qx = qx + As[i]*exp_tab_q[i];  qx = qx + As [i] * exp_tab_q [i];
qy = qy + As[i]*exp_tab_q[33-i];  qy = qy + As [i] * exp_tab_q [33-i];
end for  end for
Q = l/sqrt(qx*qx+qy*qy) où As[i]= â; ' -Q = l / sqrt (qx * qx + qy * qy) where As [i] = â ; '-
Sans perte de généralité, on pourra calculer les coefficients du filtre d'ordre 2 autrement, par exemple en appliquant au filtre LPC A(z) d'ordre 16 la procédure de réduction de l'ordreWithout loss of generality, it is possible to calculate the coefficients of the second order filter otherwise, for example by applying to the LPC filter A (z) of order 16 the order reduction procedure.
LPC dite « STEP DOWN » décrite dans J.D. Markel and A.H. Gray, Linear Prédiction of Speech, Springer Verlag, 1976 ou en effectuant deux itérations d'algorithme de Levinson- Durbin (ou STEP-UP) à partir des autocorrélations calculées sur le signal synthétisé (décodé) à 12.8 kHz et fenêtré. LPC called "STEP DOWN" described in JD Markel and AH Gray, Linear Prediction of Speech, Springer Verlag, 1976 or by performing two Levinson-Durbin (or STEP-UP) algorithm iterations from the autocorrelations calculated on the synthesized signal (decoded) at 12.8 kHz and windowed.
Pour certains signaux, la quantité Q , calculée à partir des 3 premiers coefficients LPC décodés, prend mieux en compte l'influence de la pente spectrale (ou tilt) dans le spectre et évite l'influence de pics ou de vallées « parasites » proches de 6000 Hz qui peuvent biaiser ou élever la valeur de la quantité R , calculée à partir de tous les coefficients LPC. Dans un mode de réalisation privilégié, le facteur d'échelle optimisé est déduit des quantités pré-calculées R, P, Qde façon conditionnelle, comme suit : For some signals, the quantity Q, calculated from the first 3 decoded LPC coefficients, takes better account of the influence of the spectral slope (or tilt) in the spectrum and avoids the influence of nearby "parasitic" peaks or valleys 6000 Hz which can bias or raise the value of the quantity R, calculated from all the LPC coefficients. In a preferred embodiment, the optimized scaling factor is derived from the precomputed quantities R, P, Q conditionally as follows:
Si le tilt (calculé comme dans AMR-WB dans le bloc 104, par autocorrélation normalisée sous la forme r(l)/r(0) où r(i) est l'autocorrélation) est négatif (tilt <0 comme représenté à la figure 5b), le calcul du facteur d'échelle se fait de la façon suivante: If the tilt (calculated as in AMR-WB in block 104, by normalized autocorrelation in the form r (l) / r (0) where r (i) is the autocorrelation) is negative (tilt <0 as shown in FIG. figure 5b), the calculation of the scale factor is as follows:
Pour éviter des artefacts dus à des variations trop brusques d'énergie de la bande haute, on applique un lissage à la valeur de R . Dans un mode de réalisation privilégié, un lissage exponentiel est effectué avec un facteur fixe dans le temps (0.5) sous la forme :  To avoid artefacts due to abrupt changes in energy of the high band, a smoothing is applied to the value of R. In a preferred embodiment, an exponential smoothing is performed with a fixed factor in time (0.5) in the form:
R = 0.5R + 0.5Rprev R = 0.5R + 0.5R prev
Rprev = R R prev = R
où R correspond à la valeur de R dans la sous-trame précédente et le facteur 0.5 est optimisé de façon empirique - bien entendu, le facteur 0.5 pourra être changé pour une autre valeur et d'autres méthodes de lissage sont également possibles. A noter que le lissage permet de réduire les variantes temporelles et évite donc des artéfacts. where R is the value of R in the preceding sub-frame and the factor 0.5 is empirically optimized - of course, the factor 0.5 can be changed to another value and other smoothing methods are also possible. Note that smoothing reduces temporal variations and therefore avoids artifacts.
Le facteur d'échelle optimisé est alors donné par : gHB2 (m) = max(min(R, Q), P) I P The optimized scaling factor is then given by: g HB2 (m) = max (min (R, Q), P) IP
Dans un mode de réalisation alternative, on pourra remplacer le lissage de R par un lissage de gHB2(m) tel que: In an alternative embodiment, it is possible to replace the smoothing of R by a smoothing of g HB2 (m) such that:
8HBI (m) <- °-58HB2 (m) + °-58HB2 (M - !) 8HBI ( m ) <-? - 5 8HB2 ( m ) + ° - 5 8HB2 ( M -!)
Si le tilt (calculé comme dans AMR-WB dans le bloc 104) est positif (tilt>0 comme à la figure 5a), le calcul du facteur d'échelle se fait de la façon suivante: If the tilt (calculated as in AMR-WB in block 104) is positive (tilt> 0 as in Figure 5a), the scale factor is calculated as follows:
La quantité R est lissée de façon adaptative dans le temps, avec un lissage plus fort quand R est faible - comme dans le cas précédent, ce lissage permet de réduire les variantes temporelles et évite donc des artéfacts:  The quantity R is smoothed adaptively over time, with a stronger smoothing when R is weak - as in the previous case, this smoothing makes it possible to reduce the temporal variants and thus avoids artifacts:
R = (l - a)R + aRprev avec a = l - R2 Rprev = R R = (l - a) R + aR prev with a = 1 - R 2 R prev = R
Ensuite, le facteur d'échelle optimisé est donné par :  Then the optimized scaling factor is given by:
gHB2(m) = mm{R, P, Q) / P Dans un mode de réalisation alternative, on pourra remplacer le lissage de R par un lissage de gHB2 (m) tel que calculé ci-dessus. gHB m) + SHB m = 0,..., 3 , (X = l- gH 2 B (m) g HB2 (m) = mm (R, P, Q) / P In an alternative embodiment, the smoothing of R may be replaced by a smoothing of g HB2 (m) as calculated above. gHB m ) + S H B m = 0, ..., 3, ( X = l- g H 2 B (m)
ou gHB (-ï) est le facteur d'échelle ou gain calculé pour la dernière sous-trame de la trame précédente. or g HB (-i) is the scale factor or gain calculated for the last subframe of the previous frame.
On prend ici le minimum de R, P, Q afin d'éviter de surestimer le facteur d'échelle.  Here we take the minimum of R, P, Q in order to avoid overestimating the scale factor.
Dans une variante, la condition ci-dessus dépendant uniquement du tilt pourra être étendue pour tenir compte non seulement du paramètre de tilt mais également d'autres paramètres afin d'affiner la décision. De plus, le calcul de gHB2 (m) pourra être ajusté en fonction de ces dits paramètres supplémentaires. Alternatively, the above tilt-only condition may be extended to take into account not only the tilt parameter but also other parameters to refine the decision. In addition, the calculation of g HB2 (m) can be adjusted according to these additional parameters.
Un exemple de paramètre supplémentaire est le nombre de passage par zéro (ZCR, zéro crossing rate) qui peut être défini comme :  An example of an additional parameter is the zero crossing rate (ZCR) which can be defined as:
N-l  N-l
= -∑| sgn [s(n)] - sgn [s(n - l)] | OÙ = -Σ | sgn [s (n)] - sgn [s (n - l)] | OR
Le paramètre zcr donne généralement les résultats similaires au tilt. Un bon critère de classification est le ratio entre zcrs calculé pour le signal synthétisé s{ri) et zcru calculé pour le signal d'excitation u(n) à 12800 Hz. Ce ratio est entre 0 et 1, où 0 signifie que le signal a un spectre décroissant, 1 que le spectre est croissant (ce qui correspond à (1 - tilt) I ' 2 . Dans ce cas, un ratio zcr / zcru >0.5 correspond au cas tilt <0, un ratio zcr / zcru <0.5 correspond à tilt >0. The zcr parameter usually gives the results similar to the tilt. A good classification criterion is the ratio between zcr s calculated for the synthesized signal s {ri) and zcr u calculated for the excitation signal u (n) at 12800 Hz. This ratio is between 0 and 1, where 0 means that the signal has a decreasing spectrum, 1 that the spectrum is increasing (which corresponds to (1 - tilt) I ' 2. In this case, a ratio zcr / zcr u > 0.5 corresponds to the case tilt <0, a ratio zcr / zcr u <0.5 corresponds to tilt> 0.
Dans une variante, on pourra utiliser une fonction d'un paramètre tilth où tilthp est le tilt calculé pour le signal synthétisé s(n) filtré par un filtre passe haut avec une fréquence de coupure par exemple à 4800 Hz ; dans ce cas, la réponse 1 / Â(z/ γ) de 6 à 8 kHz (appliquée à 16 kHz) correspond à la réponse pondérée de 1 / A(z) de 4.8 à 6.4 kHz. Comme 1 / A(z/ γ) a une réponse plus aplatie, il faut compenser ce changement de tilt. La fonction de facteur d'échelle selon tilthp est alors donnée dans un mode de réalisation par : In a variant, a function of a parameter tilt h may be used where tilt hp is the calculated tilt for the synthesized signal s (n) filtered by a high-pass filter with a cut-off frequency, for example at 4800 Hz; in this case, the 1 / Â (z / γ) response of 6 to 8 kHz (applied at 16 kHz) corresponds to the weighted response of 1 / A (z) of 4.8 to 6.4 kHz. Since 1 / A (z / γ) has a more flattened response, you have to compensate for this change in tilt. The scale factor function according to tilt hp is then given in one embodiment by:
(l - tilt hpf + 0.6 . On multiplie donc Q et R par mm {l,(l - tilt hp ) + 0.ôj quand tilt >0 ou par max (l,(l -tilthp ) + 0.ôj quand tilt <0. On considère maintenant le cas du débit de 23.85kbit/s, pour lequel une correction de gain est réalisée par les blocs 403 à 408. Cette correction de gain pourrait d'ailleurs faire l'objet d'une invention séparée. Dans ce mode particulier selon l'invention, l'information de correction de gain, notée gHBcorr (m) , transmise par le codage (compatible) AMR-WB avec un débit de 0.8 kbit/s est utilisée pour améliorer la qualité à 23.85 kbit/s. (l - tilt hp f + 0.6, so multiply Q and R by mm {l, (l - tilt hp ) + 0.ôj when tilt> 0 or by max (l, (l -tilt hp ) + 0.ôj when tilt <0. We now consider the case of the 23.85kbit / s rate, for which a gain correction is performed by the blocks 403 to 408. This gain correction could also be the subject of a separate invention. In this particular mode according to the invention, the gain correction information, denoted g HBcorr (m), transmitted by the encoding (compatible) AMR-WB with a bit rate of 0.8 kbit / s is used to improve the quality at 23.85. kbit / s.
On suppose ici que le codage (compatible) AMR-WB a effectué une quantification de gain de correction sur 4 bits comme décrit dans la clause UIT-T G.722.2/5.11 ou de façon équivalente dans la clause 3GPP TS 26.190/5.11.  It is assumed here that the (compatible) AMR-WB encoding performed 4-bit correction gain quantization as described in clause ITU-T G.722.2 / 5.11 or equivalently in clause 3GPP TS 26.190 / 5.11.
Dans le codeur AMR-WB, le gain de correction est calculé en comparant l'énergie du signal original échantillonné à 16 kHz et filtré par un filtre passe-bande 6-7 kHz, sHB (n) avec l'énergie du bruit blanc à 16 kHz filtré par un filtre de synthèse 1 / A(z I γ) et un filtre passe- bande 6-7 kHz (avant le filtrage l'énergie du bruit est mis à un niveau semblable à celui de l'excitation à 12.8 kHz), sHB2 (n) . Le gain est la racine du ratio d'énergie du signal original sur l'énergie du bruit divisé par deux. Dans un mode possible de réalisation on pourra changer le filtre passe-bande pour un filtre avec une bande plus large (par exemple de 6 à 7.6 kHz). In the AMR-WB encoder, the correction gain is calculated by comparing the energy of the sampled original signal at 16 kHz and filtered by a 6-7 kHz bandpass filter, s HB (n) with white noise energy at 16 kHz filtered by a synthesis filter 1 / A (z I γ) and a bandpass filter 6-7 kHz (before filtering the noise energy is set at a level similar to that of the excitation at 12.8 kHz), s HB2 (n). The gain is the root of the energy ratio of the original signal on the energy of the noise divided by two. In a possible embodiment, the bandpass filter can be changed for a filter with a wider band (for example from 6 to 7.6 kHz).
, {m) = , m = 0, ..., 3 , {m) =, m = 0, ..., 3
Pour pouvoir appliquer l'information de gain reçue à 23.85 kbit/s (dans le bloc 407), il est important de ramener l'excitation à un niveau similaire à celui attendu au codage (compatible) AMR-WB. Ainsi, le bloc 404 effectue la mise à l'échelle du signal d'excitation selon l'équation suivante:  In order to be able to apply the received gain information at 23.85 kbit / s (in block 407), it is important to reduce the excitation to a level similar to that expected at AMR-WB (compatible) coding. Thus, block 404 scales the excitation signal according to the following equation:
UHBI (n) = SHB3 (m)uHB (n) , n = 80m, · · · , 80(m + 1) - 1 U HBI ( n ) = SHB3 ( m ) u HB (n), n = 80m, · · ·, 80 (m + 1) -1
où gHB3 (m) est un gain par sous-trame calculé dans le bloc 403 sous la forme : where g HB3 (m) is a subframe gain calculated in block 403 in the form:
où le facteur 5 au dénominateur sert à compenser la différence de largeur de bande entre le signal u(n) et le signal uHB (n) , sachant qu'au codage AMR-WB l'excitation HF est un bruit blanc sur la bande 0-8000 Hz. L'indice de 4 bits par sous-trame, noté indexHF gain (m) , envoyé à 23.85 kbit/s est démultiplexé du train binaire (bloc 405) et décodé par le bloc 406 de la façon suivante :where the factor 5 at the denominator serves to compensate for the difference in bandwidth between the signal u (n) and the signal u HB (n), knowing that at the AMR-WB coding the HF excitation is a white noise on the band 0-8000 Hz. The index of 4 bits per subframe, denoted index HF gain (m), sent at 23.85 kbit / s is demultiplexed from the bitstream (block 405) and decoded by block 406 as follows:
<-> HBcorr (m) = 2.HP_gain(indexHF gain (m)) <-> HBcorr (m) = 2.HP_gain ( HF gain index (m))
où HP _gain(.) est le dictionnaire de quantification de gain HF défini dans le codage AMR- WB et rappelé ci-dessous :  where HP _gain (.) is the HF gain quantization dictionary defined in the AMR-WB encoding and recalled below:
Le bloc 407 effectue la mise à l'échelle du signal d'excitation selon l'équation suivante:Block 407 scales the excitation signal according to the following equation:
2 (n) = S HBcorr (m)uHBl (ji) , n = 80m, · · · , 80(m + 1) - 1 2 ( n ) = S HBcorr ( m ) u HB1 (ji), n = 80m, · · ·, 80 (m + 1) -1
Enfin, on ajuste l'énergie de l'excitation au niveau de la sous-trame courante avec les conditions suivantes (bloc 408). On calcule :  Finally, the excitation energy is adjusted at the level of the current subframe with the following conditions (block 408). We calculate:
∑{ g (m)gHB2 (m)uHB (n)) Σ {g (m) g HB2 (m) u HB (n))
fac(m)  fac (m)
∑½B2 («) Σ½ B2 («)
Le numérateur représente ici l'énergie de signal bande-haute qui serait obtenue dans le mode 23.05. Comme expliqué avant, pour les débits <23.85 kbit/s il faut conserver le niveau d'énergie entre le signal d'excitation décodé et le signal d'excitation étendu uHB (n) , mais cette contrainte n'est pas nécessaire dans le cas du débit de 23.85 kbit/s, puisque uHB (n) est dans ce cas mis à l'échelle par le gain gHB3 (m) . Pour éviter les doubles multiplications certaines opérations de multiplications appliqués au signal dans le bloc 400 sont appliquées dans le bloc 402 en multipliant par g (m) . La valeur de g (m) dépend de l'algorithme de synthèse de uHB (n) et doit être ajusté de telle sorte que le niveau d'énergie entre le signal d'excitation décodé en bande basse et le signal g (m)uHB (n) soit conservé. Dans un mode de réalisation particulier, qui sera décrit en détail plus tard en référence à la figure 7, g (m) = 0.6gHB1 (m) , où gHB1 (m) est un gain qui assure, pour le signal uHB , le même ratio entre énergie par sous-trame et énergie par trame que pour le signal u(n) etThe numerator represents here the band-high signal energy that would be obtained in the 23.05 mode. As explained above, for rates <23.85 kbit / s the energy level between the decoded excitation signal and the extended excitation signal u HB (n) must be maintained, but this constraint is not necessary in the case of the 23.85 kbit / s rate, since u HB (n) is in this case scaled by the gain g HB3 (m). To avoid double multiplications, certain multiplication operations applied to the signal in block 400 are applied in block 402 by multiplying by g (m). The value of g (m) depends on the synthesis algorithm of u HB (n) and must be adjusted so that the energy level between the decoded low band excitation signal and the g (m) signal u HB (n) is preserved. In a particular embodiment, which will be described in detail later with reference to FIG. 7, g (m) = 0.6 g HB1 (m), where g HB1 (m) is a gain which ensures, for the signal u HB , the same ratio between energy per subframe and energy per frame as for the signal u (n) and
0.6 correspond à la valeur moyenne d'amplitude de réponse en fréquence du filtre de désaccentuation de 5000 à 6400 Hz. 0.6 corresponds to the average value of frequency response amplitude of the deemphasis filter from 5000 to 6400 Hz.
On suppose que dans le bloc 408 on dispose d'une information sur le tilt du signal bande basse - dans un mode de réalisation privilégié ce tilt est calculé comme dans le codée AMR- WB selon les blocs 103 et 104, cependant d'autres méthodes d'estimation du tilt sont possibles sans changer le principe de l'invention.  It is assumed that in block 408 there is information on the tilt of the low band signal - in a preferred embodiment this tilt is calculated as in the AMR-WB coding according to blocks 103 and 104, however other methods tilt estimation are possible without changing the principle of the invention.
Si fac(m) >1 ou tilt<0, on prend : If fac (m)> 1 or tilt <0, we take:
UHB = UHB2 (n ' n = 80m, · · · , 80(m + 1) - 1 U HB = U HB2 ( n ' n = 80m, · · ·, 80 (m + 1) -1
Sinon : If not :
uHB \n) = max (Vl - tilt , fac(m)^ MHB2(n) , n = 80m, - · · , 80(m + 1) - 1 u HB ( n) = max (Vl - tilt, fac (m) ^ M HB2 (n), n = 80m, - · ·, 80 (m + 1) - 1
On notera que le calcul de facteur d'échelle optimisé présenté ici, notamment dans les blocs 401 et 402, se distingue de l'égalisation précitée de niveaux de filtres effectuée dans le codée AMR-WB+ par plusieurs aspects :  It will be noted that the optimized scale factor calculation presented here, in particular in the blocks 401 and 402, differs from the aforementioned equalization of filter levels performed in the AMR-WB + codec by several aspects:
• Le facteur d'échelle optimisé est calculé directement à partir des fonctions de  • The optimized scaling factor is calculated directly from the functions of
transfert des filtres LPC sans impliquer de filtrage temporel. Ceci simplifie le procédé. transfer LPC filters without involving temporal filtering. This simplifies the process.
• L'égalisation est faite de préférentielle à une fréquence différente de la fréquence de Nyquist (6400 Hz) associée à la bande basse. En effet, la modélisation LPC représente implicitement l'atténuation du signal typiquement causée par les opérations de ré-échantillonnage et donc la réponse en fréquence d'un filtre LPC peut être subir à la fréquence de Nyquist une diminution qui se retrouve pas à la fréquence commune choisie. • Equalization is preferential at a frequency different from the Nyquist frequency (6400 Hz) associated with the low band. Indeed, the LPC modeling implicitly represents the attenuation of the signal typically caused by the resampling operations and therefore the frequency response of an LPC filter can be subjected to the frequency of Nyquist a decrease which is not found at the frequency chosen commune.
· L'égalisation repose ici sur un filtre d'ordre moins élevé (ici d'ordre 2) en plus des 2 filtres à égaliser. Ce filtre additionnel permet d'éviter les effets de fluctuations spectrales locales (pic ou vallée) qui peuvent être présentes à la fréquence commune pour le calcul de la réponse en fréquence des filtres de prédiction.  · Equalization here is based on a lower order filter (here of order 2) in addition to 2 filters to equalize. This additional filter makes it possible to avoid the effects of local spectral fluctuations (peak or valley) that may be present at the common frequency for calculating the frequency response of the prediction filters.
Pour les blocs 403 à 408, l'avantage de l'invention est que la qualité du signal décodé à 23.85 kbit/s selon l'invention est améliorée par rapport à un signal décodé à 23.05 kbit/s, ce qui n'est pas le cas dans un décodeur AMR-WB. En fait, cet aspect de l'invention permet d'utiliser l'information supplémentaire (0.8 kbit/s) reçue à 23.85 kbit/s, mais de façon contrôlée (bloc 408), pour améliorer la qualité du signal d'excitation étendu au débit de 23.85. For the blocks 403 to 408, the advantage of the invention is that the quality of the decoded signal at 23.85 kbit / s according to the invention is improved compared to a signal decoded at 23.05 kbit / s, which is not the case in an AMR-WB decoder. In fact, this aspect of the invention makes it possible to use the additional information (0.8 kbit / s) received at 23.85 kbit / s, but in a controlled manner (block 408), to improve the quality of the excitation signal extended to flow rate of 23.85.
Le dispositif de détermination du facteur d'échelle optimisé tel qu'illustré par les blocs 401 à 408 de la figure 4, met en œuvre un procédé de détermination du facteur d'échelle optimisé décrit maintenant en référence à la figure 6. The optimized scale factor determining device as illustrated by blocks 401 to 408 of Fig. 4 implements an optimized scale factor determination method now described with reference to Fig. 6.
Les étapes principales sont mises en œuvre par le bloc 401. Ainsi, un signal d'excitation étendu uHB(n) esX. obtenu lors d'un procédé d'extension de bande de fréquence E601 qui comporte une étape de décodage ou d'extraction, dans une première bande de fréquence dite bande basse, d'un signal d'excitation et de paramètres de la première bande de fréquence comme par exemple les coefficients du filtre de prédiction linéaire de la première bande de fréquence. The main steps are implemented by block 401. Thus, an extended excitation signal u HB (n) esX. obtained during a frequency band extension method E601 which comprises a step of decoding or extracting, in a first so-called low band frequency band, an excitation signal and parameters of the first band of such as, for example, the coefficients of the linear prediction filter of the first frequency band.
Une étape E602 détermine un filtre de prédiction linéaire dit filtre additionnel, d'ordre inférieur à celui de la première bande de fréquence. Pour déterminer ce filtre, les paramètres de la première bande de fréquence décodés ou extraits sont utilisés.  A step E602 determines a linear prediction filter called additional filter, of a lower order than that of the first frequency band. To determine this filter, the parameters of the first decoded or extracted frequency band are used.
Dans un mode de réalisation cette étape est effectuée par troncature de la fonction de transfert du filtre de prédiction linéaire de la bande basse pour obtenir un ordre de filtre inférieur, par exemple 2. Ces coefficients peuvent ensuite être modifiés en fonction d'un critère de stabilité comme expliqué précédemment en référence à la figure 4.  In one embodiment, this step is performed by truncation of the transfer function of the linear prediction filter of the low band to obtain a lower filter order, for example 2. These coefficients can then be modified according to a criterion of stability as explained above with reference to FIG. 4.
A partir des coefficients du filtre additionnel ainsi déterminé, une étape E603 est mise en œuvre pour calculer le facteur d'échelle optimisé à appliquer au signal d'excitation étendu. Ce facteur d'échelle optimisé est par exemple calculé à partir de la réponse en fréquence du filtre additionnel à une fréquence commune entre la bande basse (première bande de fréquence) et la bande haute (deuxième bande de fréquence). Une valeur minimale pouvant être choisie entre la réponse en fréquence de ce filtre et celles des filtres bande basse et bande haute.  From the coefficients of the additional filter thus determined, a step E603 is implemented to calculate the optimized scaling factor to be applied to the extended excitation signal. This optimized scale factor is for example calculated from the frequency response of the additional filter at a common frequency between the low band (first frequency band) and the high band (second frequency band). A minimum value that can be chosen between the frequency response of this filter and those of the low band and high band filters.
Ceci évite donc les surestimations d'énergie qui pouvaient exister dans les méthodes de l'état de l'art. This avoids the overestimation of energy that could exist in state-of-the-art methods.
Cette étape de calcul du facteur d'échelle optimisé est par exemple décrite précédemment en référence à la figure 4 et aux figures 5a et 5b.  This step of calculating the optimized scale factor is for example described above with reference to FIG. 4 and FIGS. 5a and 5b.
L'étape E604 réalisée par le bloc 402 ou 409 (selon le débit de décodage) pour l'extension de bande, applique le facteur d'échelle optimisé ainsi calculée au signal d'excitation étendu de façon à obtenir un signal d'extension étendu optimisé uHB'(n). Step E604 performed by block 402 or 409 (depending on the decoding rate) for the band extension, applies the optimized scaling factor thus calculated to the extended excitation signal so as to obtain an extended extension signal. optimized u HB '(n).
Dans un mode de réalisation particulier, le dispositif de détermination du facteur d'échelle optimisé 708 est intégré dans un dispositif d'extension de bande décrit maintenant en référence à la figure 7. Ce dispositif de détermination du facteur d'échelle optimisé illustré par le bloc 708 met en œuvre le procédé de détermination du facteur d'échelle optimisé décrit précédemment en référence à la figure 6. In a particular embodiment, the optimized scaling factor device 708 is integrated into a tape expansion device now described with reference to FIG. 7. This optimized scale factor determining device illustrated by FIG. block 708 implements the method of determining the optimized scale factor described above with reference to FIG. 6.
Dans ce mode de réalisation, le bloc 400 d'extension de bande de la figure 4 comprend les blocs 700 à 707 de la figure 7 décrite maintenant.  In this embodiment, the band extension block 400 of FIG. 4 comprises the blocks 700 to 707 of FIG. 7 now described.
Ainsi, à l'entrée du dispositif d'extension de bande, un signal d'excitation bande basse décodé ou estimé par analyse est reçu ( u(n) ). L'extension de bande utilise ici l'excitation décodée à 12.8 kHz (exc2 ou u(n) ) en sortie du bloc 302 de la figure 3. On notera que dans ce mode de réalisation, la génération de l'excitation sur- échantillon né et étendu s'effectue dans une bande de fréquence allant de 5 à 8 kHz incluant donc une deuxième bande de fréquence (6.4-8kHz) supérieure à la première bande de fréquence (0-6.4 kHz). Thus, at the input of the band extender, a decoded or analytically estimated low band excitation signal is received (u (n)). The band extension here uses the decoded excitation at 12.8 kHz (exc2 or u (n)) at the output of the block 302 of FIG. It should be noted that in this embodiment, the generation of the over-sampled and extended excitation occurs in a frequency band ranging from 5 to 8 kHz, thus including a second frequency band (6.4-8 kHz) greater than first frequency band (0-6.4 kHz).
Ainsi, la génération d'un signal d'excitation étendu s'effectue au moins deuxième bande de fréquence mais aussi sur une partie de la première bande de fréquence.  Thus, the generation of an extended excitation signal is carried out at least second frequency band but also on a part of the first frequency band.
Bien évidemment, les valeurs définissant ces bandes de fréquences peuvent être différentes selon le décodeur ou le dispositif de traitement dans lequel l'invention s'applique.  Of course, the values defining these frequency bands may be different depending on the decoder or the processing device in which the invention applies.
Pour cet exemple de réalisation, ce signal est transformé pour obtenir un spectre de signal d'excitation U(k) par le module de transformation temps-fréquence 500.  For this exemplary embodiment, this signal is transformed to obtain an excitation signal spectrum U (k) by the time-frequency transformation module 500.
Dans un mode de réalisation particulier, la transformée utilise une DCT-IV (pour "Discrète Cosine Transform" - Type IV en anglais) (bloc 700) sur la trame courante de 20 ms (256 échantillons), sans fenêtrage, ce qui revient à transformer directement u(n) avec n = 0 · · · , 255 selon la formule suivante : In a particular embodiment, the transform uses a DCT-IV (for "Discrete Cosine Transform" - Type IV in English) (block 700) on the current frame of 20 ms (256 samples), without windowing, which amounts to directly transform u (n) with n = 0 · · ·, 255 according to the following formula:
où N = 256 et & = (),· · · , 255 . where N = 256 and & = (), · · ·, 255.
On note ici que la transformation sans fenêtrage (ou de façon équivalente avec une fenêtre rectangulaire implicite de la longueur de la trame) est possible car le traitement est effectué dans le domaine de l'excitation, et non le domaine du signal, si bien qu'aucun artefact (effets de bloc) n'est audible, ce qui constitue un avantage important de ce mode de réalisation de l'invention.  We note here that the transformation without windowing (or equivalently with a rectangular window implicit in the length of the frame) is possible because the processing is performed in the field of excitation, and not the domain of the signal, so that no artifact (block effects) is audible, which is an important advantage of this embodiment of the invention.
Dans ce mode de réalisation, la transformation DCT-IV est mise en œuvre par FFT suivant l'algorithme dit « Evolved ZXT(EDCT) » décrit dans l'article de D.M. Zhang, HT. Li, A Low Complexity Transform - Evolved DCT, IEEE 14th International Conférence on Computational Science and Engineering (CSE), Aug. 2011, pp. 144-149, et mis en œuvre dans les normes UIT-T G.718 Annex B et G.729.1 Annex E. In this embodiment, the DCT-IV transformation is implemented by FFT according to the "Evolved ZXT (EDCT)" algorithm described in the article by D. M. Zhang, HT. Li, A Low Complexity Transform - Evolved DCT, IEEE 14th International Conference on Computational Science and Engineering (CSE), Aug. 2011, pp. 144-149, and implemented in ITU-T G.718 Annex B and G.729.1 Annex E.
Dans des variantes de l'invention et sans perte de généralité, la transformation DCT- IV pourra être remplacée par d'autres transformations temps-fréquences court-terme de même longueur et dans le domaine de l'excitation, comme une FFT (pour "Fast Fourier Transform" en anglais) ou une DCT-II {Discrète Cosine Transform - Type II). De façon alternative, on pourra remplacer la DCT-IV sur la trame par une transformation avec recouvrement-addition et fenêtrage de longueur supérieure à la longueur de la trame courante, par exemple en utilisant une MDCT (pour "Modified Discrète Cosine Tranform" en anglais). Dans ce cas le retard dans le bloc 310 de la figure 3, devra être ajusté (réduit) de façon adéquate en fonction du retard additionnel dû à l'analyse/synthèse par cette transformée. In variants of the invention and without loss of generality, the DCT-IV transformation may be replaced by other short-term time-frequency transformations of the same length and in the field of excitation, such as an FFT (for Fast Fourier Transform "in English) or DCT-II (Discrete Cosine Transform - Type II). Alternatively, it will be possible to replace the DCT-IV on the frame by a recovery-addition and windowing transformation of length greater than the length of the current frame, for example by using an MDCT (for "Modified Discrete Cosine Tranform"). ). In this case, the delay in block 310 of FIG. 3 will have to be adjusted (reduced) by adequately depending on the additional delay due to the analysis / synthesis by this transform.
Le spectre DCT, U(k) , de 256 échantillons couvrant la bande 0-6400 Hz (à 12.8 kHz), est ensuite étendu (bloc 701) en un spectre de 320 échantillons couvrant la bande 0- 8000 Hz (à 16 kHz) sous la forme suivante :  The DCT spectrum, U (k), of 256 samples covering the band 0-6400 Hz (at 12.8 kHz), is then extended (block 701) into a spectrum of 320 samples covering the band 0-8000 Hz (at 16 kHz) in the following form:
UHBl (k) =U HBl (k) =
où on prend de façon préférentielle start_band = 160. where preferentially start_band = 160 is used.
Le bloc 701 fonctionne comme module de génération d'un signal d'excitation suréchantillonné et étendu et réalise un ré-échantillonnage de 12.8 à 16 kHz dans le domaine fréquentiel, en rajoutant ¼ d'échantillons ( k = 240, · · · , 319 ) au spectre, le ratio entre 16 et 12.8 étant de 5/4.  Block 701 functions as a module for generating an oversampled and extended excitation signal and performs a resampling of 12.8 to 16 kHz in the frequency domain, adding ¼ samples (k = 240, · · ·, 319 ) to the spectrum, the ratio between 16 and 12.8 being 5/4.
De plus, le bloc 701 réalise un filtrage passe-haut implicite dans la bande 0-5000 Hz puisque les 200 premiers échantillons de UHB1 (k) sont mis à zéro ; comme expliqué plus tard, ce filtrage passe-haut est également complété par une partie d'atténuation progressive des valeurs spectrales d'indices k = 200, · · · , 255 dans la bande 5000-6400 Hz, cette atténuation progressive est mise en œuvre dans le bloc 704 mais pourrait être réalisée séparément en dehors du bloc 704. De façon équivalente et dans des variantes de l'invention, la mise en œuvre du filtrage passe-haut séparée en blocs de coefficients d'indice k = 0, · · · ,199 mis à zéro, de coefficients k = 200, · · - , 255 atténués, dans le domaine transformé, pourra donc être effectué en une seule étape. In addition, block 701 performs high pass filtering implicit in the 0-5000 Hz band since the first 200 samples of U HB1 (k) are set to zero; as explained later, this high-pass filtering is also completed by a progressive attenuation part of the spectral values of indices k = 200, · · ·, 255 in the band 5000-6400 Hz, this progressive attenuation is implemented in block 704 but could be performed separately outside block 704. Equivalently and in variants of the invention, the implementation of high-pass filtering separated into blocks of coefficients of index k = 0, · · ·, 199 set to zero, with coefficients k = 200, · · -, 255 attenuated, in the transformed domain, can therefore be performed in a single step.
Dans cet exemple de réalisation et selon la définition de UHB1 (k) , on remarque que la bande 5000-6000 Hz de U HBl (k) (qui correspond aux indices k = 200, · · - , 239 ) est copiée à partir de la bande 5000-6000 Hz de U(k) . Cette approche permet de conserver le spectre original dans cette bande et elle évite d'introduire des distorsions dans la bande 5000-6000 Hz lors de l'addition de la synthèse HF avec la synthèse BF - en particulier la phase du signal (implicitement représentée dans le domaine DCT-IV) dans cette bande est préservée. In this exemplary embodiment and according to the definition of U HB1 (k), it is noted that the 5000-6000 Hz band of U HB1 (k) (which corresponds to the indices k = 200, · · -, 239) is copied from of the 5000-6000 Hz U (k) band. This approach preserves the original spectrum in this band and avoids introducing distortions in the 5000-6000 Hz band during the addition of HF synthesis with BF synthesis - particularly the signal phase (implicitly represented in the DCT-IV domain) in this band is preserved.
La bande 6000-8000 Hz de UHB1 (k) est ici définie en copiant la bande 4000-6000 Hz de U(k) puisque la valeur de start_band est fixée préférentiel lement à 160. The 6000-8000 Hz band of U HB1 (k) is here defined by copying the 4000-6000 Hz band of U (k) since the value of start_band is preferably fixed at 160.
Dans une variante du mode de réalisation, la valeur de start_band pourra être rendue adaptative autour de la valeur de 160. Les détails de l'adaptation de la valeur start_band ne sont pas décrits ici car ils dépassent le cadre de l'invention sans en changer la portée. Pour certains signaux en bande élargie (échantillonnés à 16 kHz), la bande haute (>6 kHz) peut être bruitée, harmonique ou comporter un mélange de bruit et d'harmoniques. De plus, le niveau d'harmonicité dans la bande 6000-8000 Hz est généralement corrélé à celui des bandes de fréquences inférieures. Ainsi le bloc 702 de génération de bruit, réalise une génération de bruit dans le domaine fréquentiel, ½BN( ) pour k = 240,· · · , 319 (80 échantillons) correspondant à une deuxième bande de fréquence dite haute fréquence afin de combiner ensuite ce bruit avec le spectre UHB1(k) dans le bloc 703. In a variant of the embodiment, the value of start_band can be made adaptive around the value of 160. The details of the adaptation of the value start_band are not described here because they go beyond the scope of the invention without changing the scope. For some broadband signals (sampled at 16 kHz), the high band (> 6 kHz) may be noisy, harmonic or have a mixture of noise and harmonics. In addition, the level of harmonicity in the 6000-8000 Hz band is generally correlated with that of the lower frequency bands. Thus, the noise generation block 702 generates frequency domain noise generation, ½ BN () for k = 240, · · ·, 319 (80 samples) corresponding to a second so-called high frequency frequency band in order to combine then this noise with the spectrum U HB1 (k) in block 703.
Dans un mode de réalisation particulier, le bruit (dans la bande 6000-8000 Hz) est généré de façon pseudo-aléatoire avec un générateur congruentiel linéaire sur 16 bits : In a particular embodiment, the noise (in the 6000-8000 Hz band) is generated pseudo-randomly with a 16-bit linear congruent generator:
avec la convention que UHBN (239) dans la trame courante correspond à la valeur UHBN (319) de la trame précédente. Dans des variantes de l'invention, on pourra remplacer cette génération de bruit par d'autres méthodes. with the convention that U HBN (239) in the current frame corresponds to the value U HBN (319) of the previous frame. In variants of the invention, this noise generation can be replaced by other methods.
Le bloc 703 de combinaison peut être réalisé de différentes façons. De façon privilégiée, on considère un mixage additif adaptatif de la forme : The combination block 703 can be realized in different ways. In a privileged way, we consider an adaptive additive mix of the form:
UHB2(k) = βϋΗΒΙ (k) + GHBNUHBN (k) , k = 240, · · · , 319 U HB2 (k) = βϋ ΗΒΙ (k) + G HBN U HBN (k), k = 240, · · ·, 319
où GHBN est un facteur de normalisation servant à égaliser le niveau d'énergie entre les deux signaux, avec ε =0.01, et le coefficient (compris entre 0 et 1) est ajusté en fonction de paramètres estimés à partir de la bande basse décodée et le coefficient β (compris entre 0 et 1) dépend de a . where G HBN is a normalization factor for equalizing the energy level between the two signals, with ε = 0.01, and the coefficient (between 0 and 1) is adjusted according to parameters estimated from the decoded low band and the coefficient β (between 0 and 1) depends on a.
Dans un mode de réalisation privilégié on calcule l'énergie du bruit dans trois bandes : 2000-4000 Hz, 4000-6000 Hz et 6000-8000 Hz, avec  In a preferred embodiment, the energy of the noise is calculated in three bands: 2000-4000 Hz, 4000-6000 Hz and 6000-8000 Hz, with
V4 = ∑ u (k) V 4 = Σ u (k)
i eN(80,159) 6 = ∑ u'2(k) i eN (80,159) 6 = Σ u ' 2 (k)
ieN(160,239)  Ien (160.239)
EN4-6 = ∑ U '2(k) E N4 - 6 = Σ U ' 2 (k)
ieN(240,319) ou Ien (240.319) or
U \k) = U(k) £ = 160,..., 239  U \ k) = U (k) £ = 160, ..., 239
239  239
fc=160  fc = 160
319 UmAk) £ = 240, ..., 319 319 U m Ak) £ = 240, ..., 319
∑ uHB1 2 (k) Σ u HB1 2 (k)
k=240 et N(fc1 , fc2) est l'ensemble des indices & pour lesquels le coefficient d'indice k est classifié comme étant associé à du bruit. Cet ensemble peut être par exemple obtenu en détectant les pics locaux dans U + et en considérant que ces raies ne sont pas associés à du bruit, soit (en appliquant la négation de la condition précédente): k = 240 and N (fc 1 , fc 2 ) is the set of indices & for which the index coefficient k is classified as being associated with noise. This set can be obtained for example by detecting the local peaks in U + and considering that these lines are not associated with noise, ie (by applying the negation of the previous condition):
N(a,b) = {a≤ k≤ b \\U + 1)|} N (a, b) = {a≤k≤ b \\ U + 1) |}
On peut noter que d'autres méthodes de calcul de l'énergie du bruit sont possibles, par exemple en prenant la valeur médiane du spectre sur la bande considérée ou en appliquant un lissage à chaque raie fréquentielle avant de calculer l'énergie par bande.  It may be noted that other methods of calculating the noise energy are possible, for example by taking the median value of the spectrum on the band in question or by applying a smoothing to each frequency line before calculating the energy per band.
On fixe a de telle sorte que le ratio entre l'énergie du bruit dans les bandes 4-6 kHz et 6-8 kHz soit le même qu'entre les bandes 2-4 kHz et 4-6 kHz : The ratio of the noise energy in the 4-6 kHz and 6-8 kHz bands is set so that between the 2-4 kHz and 4-6 kHz bands:
ou  or
EN4_6 — ΐΆάχ(ΕΝ4_6 , ΕΝ2_4), P ~ p = max(p, EN6_& ) E N4 _ 6 - ΐΆάχ (Ε Ν4 _ 6 , Ε Ν2 _ 4 ), P ~ p = max (p, E N6 _ & )
Dans des variantes de l'invention, le calcul de a pourra être remplacé par d'autres méthodes. Par exemple, dans une variante, on pourra extraire (calculer) différents paramètres (ou « features » en anglais) caractérisant le signal en bande basse, dont un paramètre « tilt » similaire à celui calculé dans le codée AMR-WB, et on estimera le facteur a en fonction d'une régression linéaire à partir de ces différents paramètres en limitant sa valeur entre 0 et 1. La régression linéaire pourra par exemple être estimée de façon supervisée en estimant le facteur a en se donnant la bande haute originale dans une base d'apprentissage. On notera que le mode de calcul de a ne limite pas la nature de l'invention. In variants of the invention, the calculation of a may be replaced by other methods. For example, in one variant, it will be possible to extract (calculate) different parameters (or "features" in English) characterizing the low band signal, including a "tilt" parameter similar to that calculated in the AMR-WB codec, and it will be estimated the factor a according to a linear regression from these different parameters by limiting its value between 0 and 1. The linear regression could for example be estimated so supervised by estimating the factor a by giving himself the original high band in a learning base. It will be noted that the method of calculating a does not limit the nature of the invention.
Dans un mode de réalisation privilégié, on prendIn a preferred embodiment, we take
afin de préserver l'énergie du signal étendu après mixage. to preserve the energy of the extended signal after mixing.
Dans une variante les facteurs β et a pourront être adaptés pour tenir compte du fait qu'un bruit injecté dans une bande donnée du signal est perçu en général comme plus fort qu'un signal harmonique à la même énergie dans la même bande. Ainsi on pourra modifier les facteurs β et comme suit:  In a variant, the factors β and a may be adapted to take account of the fact that noise injected into a given band of the signal is generally perceived as stronger than a harmonic signal at the same energy in the same band. Thus we can modify the factors β and as follows:
β - β./(α) β - β ./ (α)
a <- a.f( ) a <- a.f ()
où f(a) est une fonction décroissante de , par exemple f{ ) = b - a, [â , b = l . l , a = 1.2 , f(a) limité de 0.3 à 1. Il faut remarquer qu'après multiplication par f(a) , a2 + β2 < 1 si bien que l'énergie du signal UHB2(k) = βυ HB1(k) + GHBNU HBN(k) est plus basse que l'énergie de UHB1(k) (la différence d'énergie dépend de , plus on rajoute de bruit, plus l'énergie est atténuée). where f (a) is a decreasing function of, for example f {) = b - a, [a, b = 1. l, a = 1.2, f (a) limited from 0.3 to 1. It should be noted that after multiplication by f (a), a 2 + β 2 <1 so that the energy of the signal U HB2 (k) = βυ HB1 (k) + G HBN U HBN (k) is lower than the energy of U HB1 (k) (the energy difference depends on, the more noise is added, the more energy is attenuated).
Dans d'autres variantes de l'invention on pourra prendre :  In other variants of the invention, it will be possible to take:
β = \ - α β = \ - α
ce qui permet de préserver le niveau d'amplitude (quand les signaux combinés sont de même signe) ; cependant cette variante a le désavantage de résulter en une énergie globale (au niveau de UHB2(k) ) qui n'est pas monotone en fonction de a . which preserves the amplitude level (when the combined signals are of the same sign); however, this variant has the disadvantage of resulting in a global energy (at U HB2 (k)) which is not monotonic as a function of a.
On remarque donc ici que le bloc 703 réalise l'équivalent du bloc 101 de la figure 1 pour normaliser le bruit blanc en fonction d'une excitation qui est par contre ici dans le domaine fréquentiel, déjà étendue à la cadence de 16 kHz ; de plus, le mixage est limité à la bande 6000-8000 Hz.  It should therefore be noted here that the block 703 realizes the equivalent of the block 101 of FIG. 1 in order to normalize the white noise as a function of an excitation which is on the other hand here in the frequency domain, already extended at the rate of 16 kHz; in addition, the mix is limited to the band 6000-8000 Hz.
Dans une variante simple, on peut considérer une réalisation du bloc 703, où les spectres, UHB1(k) ou GHBNU HBN (k) , sont sélectionnés (commutés) de façon adaptative, ce qui revient à n'autoriser que les valeurs 0 ou 1 pour a ; cette approche revient à classifier le type d'excitation à générer dans la bande 6000-8000 Hz In a simple variant, it is possible to consider an embodiment of block 703, where the spectra, U HB1 (k) or G HBN U HBN (k), are selected (switched) adaptively, which amounts to allowing only the values 0 or 1 for a; this approach amounts to classifying the type of excitation to be generated in the 6000-8000 Hz band
Le bloc 704 réalise de façon optionnelle, une double opération d'application de réponse en fréquence de filtre passe-bande et de filtrage de désaccentuation (ou déemphase) dans le domaine fréquentiel. Dans une variante de l'invention, le filtrage de désaccentuation pourra être réalisé dans le domaine temporel, après le bloc 705 voire avant le bloc 700 ; cependant, dans ce cas, le filtrage passe-bande réalisé dans le bloc 704 peut laisser certaines composantes basses fréquences de très faibles niveaux qui se voient amplifier par désaccentuation, ce qui peut modifier de façon légèrement perceptible la bande basse décodée. Pour cette raison, on préfère ici réaliser la désaccentuation dans le domaine fréquentiel. Dans le mode de réalisation privilégié, les coefficients d'indice & = 0,· · · ,199 sont mis à zéro, ainsi la désaccentuation est limitée aux coefficients supérieurs. The block 704 optionally carries out a dual operation of application of bandpass filter frequency response and deemphasis filtering (or deemphasis) in the frequency domain. In a variant of the invention, the deemphasis filtering may be performed in the time domain, after block 705 or even before block 700; however, in this case, bandpass filtering performed in block 704 may leave some low frequency components of very low levels that are amplified by de-emphasis, which may slightly discern the decoded low band. For this reason, it is preferred here to perform the deemphasis in the frequency domain. In the preferred embodiment, the index coefficients & = 0, · · ·, 199 are set to zero, so the deemphasis is limited to the higher coefficients.
L'excitation est d'abord désaccentuée selon l'équation suivante : The excitation is first de-emphasized according to the following equation:
0 £ = 0,· · · ,199  0 £ = 0, · · ·, 199
U '(k) Gdeemph (k)UHB2(k) * = 200, - , 255 U '(k) G deemph (k) U HB2 (k) * = 200, -, 255
Gdeem≠(255)UHB2(k) = 256,· · · , 319 où Gdeemph (k) est la réponse en fréquence du filtre l / (l - 0.68z_1 ) sur une bande de fréquence discrète restreinte. En prenant en compte les fréquences discrètes (impaires) de la DCT-IV, on définit ici Gdeem h (k) comme: G deem 255 (255) U HB2 (k) = 256, · · ·, 319 where G deemph (k) is the frequency response of the filter l / (l - 0.68z _1 ) over a restricted discrete frequency band. Taking into account the discrete (odd) frequencies of the DCT-IV, we define here G deem h (k) as:
G W , £ = 0, · · · , 255 G W, £ = 0, · · ·, 255
e ¾ - 0.68 e ¾ - 0.68
OÙ OR
256 -80 + £ + - θ„ = 2 256 -80 + £ + - θ "= 2
256  256
Dans le cas où une autre transformation que la DCT-IV est utilisée, la définition de Ok pourra être ajustée (par exemple pour des fréquences paires). In the case where a transformation other than the DCT-IV is used, the definition of O k can be adjusted (for example for even frequencies).
On note que la désaccentuation est appliquée en deux phases pour k = 200, · · - ,255 correspondant à la bande de fréquence 5000-6400 Hz, où la réponse l / (l - 0.68z_1 ) est appliquée comme à 12.8 kHz, et pour k = 256,· · - ,319 correspondant à la bande de fréquence 6400-8000 Hz, où la réponse est étendue de 16 kHz ici à une valeur constante dans la bande 6.4-8 kHz. It is noted that the de-emphasis is applied in two phases for k = 200, · · -, 255 corresponding to the frequency band 5000-6400 Hz, where the response l / (l - 0.68z _1 ) is applied as at 12.8 kHz, and for k = 256, · · -, 319 corresponding to the 6400-8000 Hz frequency band, where the response is extended from 16 kHz here to a constant value in the 6.4-8 kHz band.
On peut noter que dans le codée AMR-WB la synthèse HF n'est pas désaccentuée. Dans le mode de réalisation présenté ici, le signal hautes fréquences est au contraire désaccentué de manière à le ramener dans un domaine cohérent avec le signal basses fréquences (0-6.4 kHz) qui sort du bloc 305 de la figure 3. Ceci est important pour l'estimation et l'ajustement ultérieur de l'énergie de la synthèse HF.  It may be noted that in the AMR-WB coding the HF synthesis is not de-emphasized. In the embodiment presented here, the high frequency signal is on the contrary de-emphasized so as to bring it back to a domain coherent with the low frequency signal (0-6.4 kHz) which leaves block 305 of FIG. 3. This is important for the estimation and the subsequent adjustment of the energy of the HF synthesis.
Dans une variante du mode de réalisation, afin de réduire la complexité, on pourra fixer Gdeem h (k) à une valeur constante indépendante de en prenant par exemple Gdeemph(k) = 0.6 ce qui correspond approximativement à la valeur moyenne de Gdeemph(k) pour k = 200,·· -,319 dans les conditions du mode de réalisation décrit ci-dessus. In a variant of the embodiment, in order to reduce the complexity, it will be possible to fix G deem h (k) at a constant value independent of, for example, by taking G deemph (k) = 0.6 which corresponds approximately to the average value of G deemph (k) for k = 200, ·· -, 319 under the conditions of the embodiment described above.
Dans une autre variante du mode de réalisation du dispositif d'extension, la désaccentuation pourra être réalisée de façon équivalente dans le domaine temporel après DCT inverse.  In another variant of the embodiment of the extension device, the de-emphasis can be performed in an equivalent way in the time domain after inverse DCT.
En plus de la désaccentuation, un filtrage passe-bande est appliqué avec deux parties séparées : l'une passe-haut fixe, l'autre passe-bas adaptative (fonction du débit).  In addition to de-emphasis, band-pass filtering is applied with two separate parts: one fixed high-pass, the other adaptive low-pass (flow-rate function).
Ce filtrage est effectué dans le domaine fréquentiel.  This filtering is performed in the frequency domain.
Dans le mode de réalisation privilégiée, on calcule la réponse partielle de filtre passe- bas dans le domaine fréquentiel comme suit:  In the preferred embodiment, the partial low-pass filter response in the frequency domain is calculated as follows:
G k) = 1-0.999—^-G k) = 1-0.999 - ^ -
Nlp-l N lp -l
où Nlp =60 à 6.6 kbit/s, 40 à 8.85 kbit/s, 20 aux débits >8.85 bit/s. where N lp = 60 to 6.6 kbit / s, 40 to 8.85 kbit / s, 20 at rates> 8.85 bit / s.
Ensuite on applique un filtre passe-bande sous la forme :  Then we apply a bandpass filter in the form:
0 £ = 0,···,199  0 £ = 0, ···, 199
Ghp(k-200)UHB2\k) £ = 200,· --,255 G hp (k-200) HB HB2 \ k) = = 200, · -, 255
UHB3(k) = U HB3 (k) =
UHB2\k) fc = 256,-,319-Nlp U HB2 \ k) fc = 256, -, 319-N lp
Glp (k - 320 -Nlv)UHB2\k) k = 320 - Nlp , · · · , 319 G lp (k-320 -N lv ) U HB2 \ k) k = 320 - N lp , · · ·, 319
La définition de Gh (k) , k = 0,··· ,55, est donnée par exemple au tableau 1 ci-dessous. The definition of G h (k), k = 0, ···, 55, is given for example in Table 1 below.
Tableau 2 On notera que dans des variantes de l'invention les valeurs de Ghp (k) pourront être modifiées tout en gardant une atténuation progressive. De même le filtrage passe-bas à largeur de bande variable, G, (k) , pourra être ajusté avec des valeurs ou un support fréquentiel différents, sans changer le principe de cette étape de filtrage. Table 2 It will be noted that in variants of the invention the values of G hp (k) may be modified while keeping a gradual attenuation. Similarly, the low-pass filtering with variable bandwidth, G, (k), may be adjusted with different values or frequency support, without changing the principle of this filtering step.
On notera aussi que le filtrage passe-bande pourra être adapté en définissant une seule étape de filtrage combinant les filtrages passe-haut et passe-bas.  Note also that the bandpass filtering can be adapted by defining a single filtering step combining the high-pass and low-pass filtering.
Dans un autre mode de réalisation, le filtrage passe-bande pourra être réalisé de façon équivalente dans le domaine temporel (comme dans le bloc 112 de la figure 1) avec différents coefficients de filtre selon le débit, après une étape de DCT inverse. Cependant, on notera qu'il est avantageux de réaliser cette étape directement dans le domaine fréquentiel car le filtrage est effectué dans le domaine de l'excitation LPC et donc les problèmes de convolution circulaire et d'effets de bord sont très limités dans ce domaine.  In another embodiment, the bandpass filtering may be performed equivalently in the time domain (as in block 112 of FIG. 1) with different filter coefficients according to the bit rate, after an inverse DCT step. However, it will be noted that it is advantageous to carry out this step directly in the frequency domain because the filtering is carried out in the field of LPC excitation and therefore the problems of circular convolution and edge effects are very limited in this field. .
On notera aussi que dans le cas du débit de 23.85 kbit/s on n'effectue pas la désaccentuation de l'excitation UHB2(k) pour rester en accord avec la façon dont le gain de correction est calculé dans le codeur AMR-WB et pour éviter les doubles multiplications. Dans ce cas block 704 réalise seulement le filtrage passe-bas. It will also be noted that in the case of the 23.85 kbit / s rate, the deemphasis of the excitation U HB2 (k) is not carried out in order to remain in agreement with the way in which the correction gain is calculated in the AMR-WB encoder. and to avoid double multiplications. In this case block 704 only performs the low-pass filtering.
Le bloc 705 de transformée inverse réalise une DCT inverse sur 320 échantillons pour trouver l'excitation haute-fréquence échantillonnée à 16 kHz. Sa mise en œuvre est identique au bloc 700, car la DCT-IV est orthonormée, sauf que la longueur de la transformée est de 320 au lieu de 256, et on obtient:The inverse transform block 705 performs an inverse DCT on 320 samples to find the high-frequency excitation sampled at 16 kHz. Its implementation is identical to block 700, because the DCT-IV is orthonormed, except that the length of the transform is 320 instead of 256, and we obtain:
où N16k = 320 et £ = (), · · · , 319 where N 16k = 320 and £ = (), · · ·, 319
Cette excitation échantillonnée à 16 kHz est ensuite de façon optionnelle mise à l'échelle par des gains définis par sous-trame de 80 échantillons (bloc 707). This excitation sampled at 16 kHz is then optionally scaled by gains defined by subframe of 80 samples (block 707).
Dans un mode de réalisation privilégié, on calcule d'abord (bloc 706) un gain gHBi(m) par sous-trame par des ratios d'énergie des sous-trames tel que dans chaque sous-trame d'indice =0, 1, 2 ou 3 de la trame courante: In a preferred embodiment, a gain g H Bi (m) per subframe is first calculated (block 706) by energy ratios of the subframes such that in each subscript of index = 0 , 1, 2 or 3 of the current frame:
où e1 (m) = u(n + 64m)2 + £ or e 1 (m) = u (n + 64m) 2 + £
79  79
e2 (m) = uHB0 (n + 80m)2 + ε e 2 (m) = u HB0 (n + 80m) 2 + ε
avec ε = 0.01. On peut écrire le gain par sous-trame gHB1 (m) sous la forme with ε = 0.01. We can write the gain by subframe g HB1 (m) in the form
ce qui montre qu'on assure dans le signal uHB le même ratio entre énergie par sous-trame et énergie par trame que dans le signal u(n) . which shows that the signal u HB has the same ratio between energy per subframe and energy per frame as in the signal u (n).
Le bloc 707 effectue la mise à l'échelle du signal combiné selon l'équation suivante:  Block 707 scales the combined signal according to the following equation:
UHB (n) = SHBI (m uHBo = 80m, · · · , 80(m + 1) - 1 UHB ( n ) = SHBI ( mu HBo = 80m, · · ·, 80 (m + 1) -1
On notera que la réalisation du bloc 706 diffère de celle du bloc 101 de la figure 1, car l'énergie au niveau de la trame courante est prise en compte en plus de celle de la sous- trame. Cela permet d'avoir le ratio de l'énergie de chaque sous-trame par rapport à l'énergie de la trame. On compare donc des ratios d'énergie (ou énergies relatives) plutôt que les énergies absolues entre bande basse et bande haute.  Note that the realization of the block 706 differs from that of the block 101 of Figure 1, because the energy at the current frame is taken into account in addition to that of the sub-frame. This makes it possible to have the ratio of the energy of each sub-frame with respect to the energy of the frame. Energy ratios (or relative energies) are compared rather than the absolute energies between low band and high band.
Ainsi, cette étape de mise à l'échelle permet de conserver dans la bande haute le ratio d'énergie entre la sous-trame et la trame de la même façon que dans la bande basse.  Thus, this scaling step makes it possible to keep in the high band the energy ratio between the subframe and the frame in the same way as in the low band.
On notera ici que dans le cas du débit 23.85 kbit/s les gains gHB1 (m) sont calculés mais appliqués dans l'étape suivante, comme expliqué en référence à la figure 4, pour éviter les doubles multiplications. Dans ce cas uHB (n) = uHB0 (n) . It will be noted here that in the case of the 23.85 kbit / s rate the gains g HB1 (m) are calculated but applied in the next step, as explained with reference to FIG. 4, to avoid double multiplications. In this case u HB (n) = u HB0 (n).
Selon l'invention, le bloc 708 effectue ensuite un calcul de facteur d'échelle par sous- trame du signal (étapes E602 à E 603 de la figure 6), comme décrit précédemment en référence à la figure 6 et détaillé en figure 4 et 5. Enfin, l'excitation corrigée uHB n) est filtrée par le module de filtrage 710 qui peut être réalisé ici en prenant comme fonction de transfert 1 / A(z/ γ) , où γ =0.9 à 6.6 kbit/s et γ =0.6 aux autres débits, ce qui limite l'ordre du filtre à l'ordre 16. According to the invention, the block 708 then performs a scaling factor calculation by subframe of the signal (steps E602 to E 603 of FIG. 6), as previously described with reference to FIG. 6 and detailed in FIG. 5. Finally, the corrected excitation u HB n) is filtered by the filtering module 710 which can be realized here by taking as transfer function 1 / A (z / γ), where γ = 0.9 to 6.6 kbit / s and γ = 0.6 at other rates, which limits the order of the filter to order 16.
Dans une variante, ce filtrage pourra être réalisé de la même façon que ce qui est décrit pour le bloc 111 de la figure 1 du décodeur AMR-WB, cependant l'ordre du filtre passe à 20 au débit de 6.6, ce qui ne change pas de façon significative la qualité du signal synthétisé. Dans une autre variante, on pourra effectuer le filtrage de synthèse LPC dans le domaine fréquentiel, après avoir calculé la réponse en fréquence du filtre mis en œuvre dans le bloc 710.  In a variant, this filtering can be done in the same way as described for the block 111 of FIG. 1 of the AMR-WB decoder, however the order of the filter goes to 20 at the rate of 6.6, which does not change. not significantly the quality of the synthesized signal. In another variant, it will be possible to carry out LPC synthesis filtering in the frequency domain, after having calculated the frequency response of the filter implemented in block 710.
Dans une variante de réalisation, l'étape de filtrage par un filtre 710 de prédiction linéaire pour la deuxième bande de fréquence est combinée à l'application du facteur d'échelle optimisé, ce qui permet de réduire la complexité de traitement. Ainsi les étapes de filtrage 1 / Λ(ζ/ γ) et d'application du facteur d'échelle optimisé gHB2 sont combinées à une seule étape de filtrage gHB2 1 (z/ y) pour réduire la complexité de traitement. In an alternative embodiment, the step of filtering by a linear prediction filter 710 for the second frequency band is combined with the application of the optimized scaling factor, which reduces the processing complexity. Thus, the filtering steps 1 / Λ (ζ / γ) and the application of the optimized scaling factor g HB2 are combined with a single filter step g HB2 1 (z / y) to reduce the processing complexity.
Dans des variantes de réalisation de l'invention, le codage de la bande basse (0-6.4 kHz) pourra être remplacé par un codeur CELP autre que celui utilisé dans AMR-WB, comme par exemple le codeur CELP dans G.718 à 8 kbit/s. Sans perte de généralité d'autres codeurs en bande élargie ou fonctionnant à des fréquences supérieurs à 16 kHz, dans lesquels le codage de la bande basse fonctionne à une fréquence interne à 12.8 kHz pourraient être utilisés. Par ailleurs, l'invention peut être adaptée de façon évidente à d'autres fréquences d'échantillonnage que 12.8 kHz, lorsqu'un codeur de basses fréquences fonctionne à une fréquence d'échantillonnage inférieure à celle du signal original ou reconstruit. Lorsque le décodage en bande basse n'utilise pas de prédiction linéaire, on ne dispose pas d'un signal d'excitation à étendre, dans ce cas on pourra réaliser une analyse LPC du signal reconstruit dans la trame courante et on calculera une excitation LPC de manière à pouvoir appliquer l'invention.  In alternative embodiments of the invention, the coding of the low band (0-6.4 kHz) may be replaced by a CELP coder other than that used in AMR-WB, for example the CELP coder in G.718 to 8. kbit / s. Without loss of generality other encoders in wide band or operating at frequencies higher than 16 kHz, in which the coding of the low band operates at an internal frequency at 12.8 kHz could be used. Moreover, the invention can be obviously adapted to other sampling frequencies than 12.8 kHz, when a low frequency encoder operates at a sampling frequency lower than that of the original or reconstructed signal. When the low band decoding does not use a linear prediction, it does not have an excitation signal to be extended, in this case it will be possible to carry out an LPC analysis of the reconstructed signal in the current frame and calculate an LPC excitation. so as to be able to apply the invention.
Enfin, dans une autre variante de l'invention, l'excitation ( u(n) ) est rééchantillonnée, par exemple par interpolation linéaire ou "spline" cubique, de 12.8 à 16 kHz avant transformation (par exemple DCT-IV) de longueur 320. Cette variante a le défaut d'être plus complexe, car la transformée (DCT-IV) de l'excitation est alors calculée sur une plus grande longueur et le ré-échantillonnage n'est pas effectué dans le domaine de la transformée.  Finally, in another variant of the invention, the excitation (u (n)) is resampled, for example by linear interpolation or cubic "spline", from 12.8 to 16 kHz before transformation (for example DCT-IV) of length 320. This variant has the defect of being more complex, because the transform (DCT-IV) of the excitation is then calculated over a greater length and the resampling is not carried out in the field of the transform.
De plus, dans des variantes de l'invention, tous les calculs nécessaires à l'estimation des gains ( GHBN , gHB1 (m) , gHB2 (m) , gHBN , ...) pourront être effectués dans un domaine logarithmique. Dans des variantes de l'extension de bande, l'excitation en bande basse u(n) et le filtre LPC 1 / A(z) seront estimés par trame, par analyse LPC d'un signal en bande basse dont la bande doit être étendue. Le signal d'excitation bande basse est alors extrait par analyse du signal audio. Moreover, in variants of the invention, all the calculations necessary for estimating the gains (G HBN , g HB1 (m), g HB2 (m), g HBN , ...) can be carried out in a field logarithmic. In variants of the band extension, the low band excitation u (n) and the LPC 1 / A (z) filter will be estimated per frame, by LPC analysis of a low band signal whose band must be extended. The low band excitation signal is then extracted by analyzing the audio signal.
Dans un mode possible de réalisation de cette variante, le signal audio bande basse est ré-échantillonné avant l'étape d'extraction de l'excitation, si bien que l'excitation extraite du signal audio (par prédiction linéaire) est déjà ré-échantillonnée.  In a possible embodiment of this variant, the low band audio signal is resampled before the excitation extraction step, so that the excitation extracted from the audio signal (by linear prediction) is already resolved. sampled.
L'extension de bande illustrée à la figure 7, s'applique dans ce cas à une bande basse qui n'est pas décodée mais analysée.  The band extension illustrated in Figure 7, in this case applies to a low band which is not decoded but analyzed.
La figure 8 représente un exemple de réalisation matérielle d'un dispositif de détermination d'un facteur d'échelle optimisé 800 selon l'invention. Celui-ci peut faire partie intégrante d'un décodeur de signal audiofréquence ou d'un équipement recevant des signaux audiofréquences décodés ou non.  FIG. 8 represents an exemplary embodiment of a device for determining an optimized scale factor 800 according to the invention. This may be an integral part of an audio-frequency signal decoder or equipment receiving decoded or non-decoded audio signals.
Ce type de dispositif comporte un processeur PROC coopérant avec un bloc mémoire BM comportant une mémoire de stockage et/ou de travail MEM.  This type of device comprises a PROC processor cooperating with a memory block BM having a memory storage and / or work MEM.
Un tel dispositif comporte un module d'entrée E apte à recevoir un signal audio d'excitation décodé ou extrait dans une première bande de fréquence dite bande basse ( u(n) ou U(k) ) et les paramètres d'un filtre de synthèse de prédiction linéaire ( A(z) ). Il comporte un module de sortie S apte à transmettre le signal haute fréquence synthétisé et optimisé (UHB'(I" ) par exemple à un module de filtrage comme le bloc 710 de la figure 7 ou à un module de ré-échantillonnage comme le module 311 de la figure 3. Such a device comprises an input module E able to receive a decoded or extracted excitation audio signal in a first so-called low band frequency band (u (n) or U (k)) and the parameters of a filter of linear prediction synthesis (A (z)). It comprises an output module S adapted to transmit the synthesized and optimized high frequency signal (U H B '(I ") for example to a filtering module such as block 710 of FIG. 7 or to a resampling module as the module 311 of FIG.
Le bloc mémoire peut avantageusement comporter un programme informatique comportant des instructions de code pour la mise en œuvre des étapes du procédé de détermination d'un facteur d'échelle optimisé à appliquer à un signal d'excitation ou à un filtre au sens de l'invention, lorsque ces instructions sont exécutées par le processeur PROC, et notamment les étapes de détermination (E602) d'un filtre de prédiction linéaire dit filtre additionnel, d'ordre inférieur au filtre de prédiction linéaire de la première bande de fréquence, les coefficients du filtre additionnel étant obtenus à partir des paramètres décodés ou extraits de la première bande de fréquence, de calcul (E603) d'un facteur d'échelle optimisé en fonction au moins des coefficients du filtre additionnel.  The memory block may advantageously comprise a computer program comprising code instructions for implementing the steps of the method for determining an optimized scale factor to be applied to an excitation signal or to a filter within the meaning of FIG. invention, when these instructions are executed by the processor PROC, and in particular the steps of determination (E602) of a linear prediction filter called additional filter, of order less than the linear prediction filter of the first frequency band, the coefficients additional filter being obtained from the parameters decoded or extracted from the first frequency band, calculation (E603) of an optimized scale factor according to at least the coefficients of the additional filter.
Typiquement, la description de la figure 6 reprend les étapes d'un algorithme d'un tel programme informatique. Le programme informatique peut également être stocké sur un support mémoire lisible par un lecteur du dispositif ou téléchargeable dans l'espace mémoire de celui-ci.  Typically, the description of FIG. 6 shows the steps of an algorithm of such a computer program. The computer program can also be stored on a memory medium readable by a reader of the device or downloadable in the memory space thereof.
La mémoire MEM enregistre de manière générale, toutes les données nécessaires à la mise en œuvre du procédé. Dans un mode possible de réalisation, le dispositif ainsi décrit peut également comporter les fonctions d'application du facteur d'échelle optimisé au signal d'excitation étendu, d'extension de bande de fréquence, de décodage bande basse et autre fonctions.de traitement décrites par exemple en figure 3 et 4 en plus des fonctions de détermination de facteur d'échelle optimisé selon l'invention. The memory MEM generally records all the data necessary for the implementation of the method. In one possible embodiment, the device thus described may also include the functions of applying the optimized scaling factor to the extended excitation signal, frequency band extension, low band decoding and other processing functions. described for example in FIGS. 3 and 4 in addition to the optimized scale factor determination functions according to the invention.

Claims

REVENDICATIONS
Procédé de détermination d'un facteur d'échelle optimisé à appliquer à un signal d'excitation ou à un filtre lors d'un procédé d'extension de bande de fréquence d'un signal audiofréquence, le procédé d'extension de bande (E601) comportant une étape de décodage ou d'extraction, dans une première bande de fréquence, d'un signal d'excitation et de paramètres de la première bande de fréquence comprenant des coefficients d'un filtre de prédiction linéaire, une étape de génération d'un signal d'excitation étendu sur au moins une deuxième bande de fréquence et une étape de filtrage par un filtre de prédiction linéaire pour la deuxième bande de fréquence, le procédé de détermination étant caractérisé en ce qu'il comporte les étapes suivantes: A method of determining an optimized scale factor to be applied to an excitation signal or a filter in a frequency band extension method of an audio frequency signal, the band extension method (E601 ) having a step of decoding or extracting, in a first frequency band, an excitation signal and parameters of the first frequency band comprising coefficients of a linear prediction filter, a step of generating a an excitation signal extended over at least a second frequency band and a filtering step by a linear prediction filter for the second frequency band, the determination method being characterized in that it comprises the following steps:
- détermination (E602) d'un filtre de prédiction linéaire dit filtre additionnel, d'ordre inférieur au filtre de prédiction linéaire de la première bande de fréquence, les coefficients du filtre additionnel étant obtenus à partir des paramètres décodés ou extraits de la première bande de fréquence; et  determination (E602) of a linear prediction filter called additional filter, of a lower order than the linear prediction filter of the first frequency band, the coefficients of the additional filter being obtained from the parameters decoded or extracted from the first band; frequency; and
- calcul (E603) du facteur d'échelle optimisé en fonction au moins des coefficients du filtre additionnel.  - calculating (E603) the scaling factor optimized according to at least the coefficients of the additional filter.
Procédé selon la revendication 1, caractérisé en ce que le procédé d'extension de bande comprend une étape d'application (E604) du facteur d'échelle optimisé au signal d'excitation étendu. Method according to claim 1, characterized in that the band extension method comprises a step of applying (E604) the optimized scaling factor to the extended excitation signal.
Procédé selon la revendication 2, caractérisé en ce que l'application du facteur d'échelle optimisé est combinée à l'étape de filtrage dans la deuxième bande de fréquence. Method according to claim 2, characterized in that the application of the optimized scaling factor is combined with the filtering step in the second frequency band.
Procédé selon la revendication 1, caractérisé en ce que les coefficients du filtre additionnel sont obtenus par troncature de la fonction de transfert du filtre de prédiction linéaire de la première bande de fréquence pour obtenir un ordre inférieur. Method according to claim 1, characterized in that the coefficients of the additional filter are obtained by truncation of the transfer function of the linear prediction filter of the first frequency band to obtain a lower order.
Procédé selon la revendication 4, caractérisé en ce que les coefficients du filtre additionnel sont modifiés en fonction d'un critère de stabilité du filtre additionnel. Method according to claim 4, characterized in that the coefficients of the additional filter are modified according to a criterion of stability of the additional filter.
Procédé selon la revendication 1, caractérisé en ce que le calcul du facteur d'échelle optimisé comporte les étapes suivantes: Method according to Claim 1, characterized in that the calculation of the optimized scaling factor comprises the following steps:
- calcul des réponses en fréquence des filtres de prédiction linéaire des première et deuxième bandes de fréquence pour une fréquence commune;  calculating the frequency responses of the linear prediction filters of the first and second frequency bands for a common frequency;
- calcul de la réponse en fréquence du filtre additionnel pour cette fréquence commune; calculation of the frequency response of the additional filter for this common frequency;
- calcul du facteur d'échelle optimisé en fonction des réponses en fréquence ainsi calculées. calculation of the optimized scale factor according to the frequency responses thus calculated.
7. Procédé selon la revendication 1, caractérisé en ce qu'il comprend en outre les étapes suivantes, mises en œuvre pour un débit de décodage prédéterminé: 7. Method according to claim 1, characterized in that it further comprises the following steps, implemented for a predetermined decoding rate:
- première mise à l'échelle du signal d'excitation étendu par un gain calculé par sous- trame fonction d'un rapport d'énergie entre le signal d'excitation décodé et le signal d'excitation étendu;  first scaling of the extended excitation signal by a calculated gain per subframe as a function of an energy ratio between the decoded excitation signal and the extended excitation signal;
-deuxième mise à l'échelle du signal d'excitation issu de la première mise à l'échelle par un gain de correction décodé;  second scaling of the excitation signal from the first scaling by a decoded correction gain;
- ajustement de l'énergie de l'excitation pour la sous-trame courante par un facteur d'ajustement calculé en fonction de l'énergie du signal obtenu après la deuxième mise à l'échelle et en fonction du signal obtenu après application du facteur d'échelle optimisé.  - adjustment of the excitation energy for the current subframe by an adjustment factor calculated as a function of the energy of the signal obtained after the second scaling and as a function of the signal obtained after application of the factor optimized scale.
8. Dispositif de détermination d'un facteur d'échelle optimisé à appliquer à un signal d'excitation ou à un filtre dans un dispositif d'extension de bande de fréquence d'un signal audiofréquence, le dispositif d'extension de bande (400) comprenant un module de décodage ou d'extraction, dans une première bande de fréquence, d'un signal d'excitation et de paramètres de la première bande de fréquence comprenant des coefficients d'un filtre de prédiction linéaire, un module de génération d'un signal d'excitation étendu sur au moins une deuxième bande de fréquence et un module de filtrage par un filtre de prédiction linéaire pour la deuxième bande de fréquence, le dispositif de détermination étant caractérisé en ce qu'il comprend: 8. A device for determining an optimized scale factor to be applied to an excitation signal or a filter in a frequency band extension device of an audio frequency signal, the band extension device (400) ) comprising a decoding or extraction module, in a first frequency band, of an excitation signal and parameters of the first frequency band comprising coefficients of a linear prediction filter, a generation module an excitation signal extended over at least a second frequency band and a filtering module by a linear prediction filter for the second frequency band, the determination device being characterized in that it comprises:
- un module de détermination (401a) d'un filtre de prédiction linéaire dit filtre additionnel, d'ordre inférieur au filtre de prédiction linéaire de la première bande de fréquence, les coefficients du filtre additionnel étant obtenus à partir des paramètres décodés ou extraits de la première bande de fréquence; et  a determination module (401a) of a linear prediction filter called an additional filter, of a lower order than the linear prediction filter of the first frequency band, the coefficients of the additional filter being obtained from the parameters decoded or extracted from the first frequency band; and
- un module de calcul (401b) du facteur d'échelle optimisé en fonction au moins des coefficients du filtre additionnel.  a computation module (401b) of the scaling factor optimized according to at least the coefficients of the additional filter.
9. Décodeur de signal audiofréquence caractérisé en ce qu'il comporte un dispositif de détermination d'un facteur d'échelle optimisé conforme à la revendication 8. 9. Audiofrequency signal decoder characterized in that it comprises an optimized scale factor determination device according to claim 8.
10. Programme informatique comportant des instructions de code pour la mise en œuvre des étapes du procédé détermination d'un facteur d'échelle optimisé selon l'une des revendications 1 à 7, lorsque ces instructions sont exécutées par un processeur. Computer program comprising code instructions for carrying out the steps of the method for determining an optimized scale factor according to one of claims 1 to 7, when these instructions are executed by a processor.
11. Support de stockage lisible par un dispositif de détermination d'un facteur d'échelle optimisé sur lequel est enregistré un programme informatique comprenant des instructions de code pour l'exécution des étapes du procédé de détermination d'un facteur d'échelle optimisé selon l'une des revendications 1 à 7. 11. A storage medium readable by an optimized scale factor determining device on which is recorded a computer program comprising code instructions for performing the steps of the method for determining an optimized scale factor according to one of claims 1 to 7.
EP14749907.3A 2013-07-12 2014-07-04 Optimized scale factor for frequency band extension in an audiofrequency signal decoder Active EP3020043B1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
FR1356909A FR3008533A1 (en) 2013-07-12 2013-07-12 OPTIMIZED SCALE FACTOR FOR FREQUENCY BAND EXTENSION IN AUDIO FREQUENCY SIGNAL DECODER
PCT/FR2014/051720 WO2015004373A1 (en) 2013-07-12 2014-07-04 Optimized scale factor for frequency band extension in an audiofrequency signal decoder

Publications (2)

Publication Number Publication Date
EP3020043A1 true EP3020043A1 (en) 2016-05-18
EP3020043B1 EP3020043B1 (en) 2017-02-08

Family

ID=49753286

Family Applications (1)

Application Number Title Priority Date Filing Date
EP14749907.3A Active EP3020043B1 (en) 2013-07-12 2014-07-04 Optimized scale factor for frequency band extension in an audiofrequency signal decoder

Country Status (11)

Country Link
US (8) US10446163B2 (en)
EP (1) EP3020043B1 (en)
JP (4) JP6487429B2 (en)
KR (4) KR102423081B1 (en)
CN (4) CN107492385B (en)
BR (4) BR122017018556B1 (en)
CA (4) CA2917795C (en)
FR (1) FR3008533A1 (en)
MX (1) MX354394B (en)
RU (4) RU2668058C2 (en)
WO (1) WO2015004373A1 (en)

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2631906A1 (en) * 2012-02-27 2013-08-28 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. Phase coherence control for harmonic signals in perceptual audio codecs
CN105976830B (en) * 2013-01-11 2019-09-20 华为技术有限公司 Audio-frequency signal coding and coding/decoding method, audio-frequency signal coding and decoding apparatus
FR3008533A1 (en) * 2013-07-12 2015-01-16 Orange OPTIMIZED SCALE FACTOR FOR FREQUENCY BAND EXTENSION IN AUDIO FREQUENCY SIGNAL DECODER
TWI557726B (en) * 2013-08-29 2016-11-11 杜比國際公司 System and method for determining a master scale factor band table for a highband signal of an audio signal
US20160323425A1 (en) * 2015-04-29 2016-11-03 Qualcomm Incorporated Enhanced voice services (evs) in 3gpp2 network
US9830921B2 (en) * 2015-08-17 2017-11-28 Qualcomm Incorporated High-band target signal control
US10825467B2 (en) * 2017-04-21 2020-11-03 Qualcomm Incorporated Non-harmonic speech detection and bandwidth extension in a multi-source environment
US20190051286A1 (en) * 2017-08-14 2019-02-14 Microsoft Technology Licensing, Llc Normalization of high band signals in network telephony communications
US10681486B2 (en) * 2017-10-18 2020-06-09 Htc Corporation Method, electronic device and recording medium for obtaining Hi-Res audio transfer information
TWI834582B (en) * 2018-01-26 2024-03-01 瑞典商都比國際公司 Method, audio processing unit and non-transitory computer readable medium for performing high frequency reconstruction of an audio signal
CN110660409A (en) * 2018-06-29 2020-01-07 华为技术有限公司 Method and device for spreading spectrum
JP2022527111A (en) * 2019-04-03 2022-05-30 ドルビー ラボラトリーズ ライセンシング コーポレイション Scalable audio scene media server
CN115136236A (en) * 2020-02-25 2022-09-30 索尼集团公司 Signal processing device, signal processing method, and program
RU2747368C1 (en) * 2020-07-13 2021-05-04 федеральное государственное казенное военное образовательное учреждение высшего образования "Военная академия связи имени Маршала Советского Союза С.М. Буденного" Министерства обороны Российской Федерации Method for monitoring and managing information security of mobile communication network
CN114333856A (en) * 2021-12-24 2022-04-12 南京西觉硕信息科技有限公司 Method, device and system for solving later-half-frame speech signal when linear prediction coefficient is given

Family Cites Families (75)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE69232202T2 (en) * 1991-06-11 2002-07-25 Qualcomm Inc VOCODER WITH VARIABLE BITRATE
US5455888A (en) * 1992-12-04 1995-10-03 Northern Telecom Limited Speech bandwidth extension method and apparatus
SE502244C2 (en) * 1993-06-11 1995-09-25 Ericsson Telefon Ab L M Method and apparatus for decoding audio signals in a system for mobile radio communication
JP3189614B2 (en) * 1995-03-13 2001-07-16 松下電器産業株式会社 Voice band expansion device
US6002352A (en) * 1997-06-24 1999-12-14 International Business Machines Corporation Method of sampling, downconverting, and digitizing a bandpass signal using a digital predictive coder
US7072832B1 (en) * 1998-08-24 2006-07-04 Mindspeed Technologies, Inc. System for speech encoding having an adaptive encoding arrangement
US6453287B1 (en) * 1999-02-04 2002-09-17 Georgia-Tech Research Corporation Apparatus and quality enhancement algorithm for mixed excitation linear predictive (MELP) and other speech coders
JP4792613B2 (en) * 1999-09-29 2011-10-12 ソニー株式会社 Information processing apparatus and method, and recording medium
FI119576B (en) * 2000-03-07 2008-12-31 Nokia Corp Speech processing device and procedure for speech processing, as well as a digital radio telephone
US6889182B2 (en) * 2001-01-12 2005-05-03 Telefonaktiebolaget L M Ericsson (Publ) Speech bandwidth extension
AUPR433901A0 (en) * 2001-04-10 2001-05-17 Lake Technology Limited High frequency signal construction method
US6732071B2 (en) * 2001-09-27 2004-05-04 Intel Corporation Method, apparatus, and system for efficient rate control in audio encoding
US7353168B2 (en) * 2001-10-03 2008-04-01 Broadcom Corporation Method and apparatus to eliminate discontinuities in adaptively filtered signals
US6895375B2 (en) * 2001-10-04 2005-05-17 At&T Corp. System for bandwidth extension of Narrow-band speech
WO2003038812A1 (en) * 2001-11-02 2003-05-08 Matsushita Electric Industrial Co., Ltd. Audio encoding and decoding device
AU2003281128A1 (en) * 2002-07-16 2004-02-02 Koninklijke Philips Electronics N.V. Audio coding
JP4676140B2 (en) * 2002-09-04 2011-04-27 マイクロソフト コーポレーション Audio quantization and inverse quantization
US7299190B2 (en) * 2002-09-04 2007-11-20 Microsoft Corporation Quantization and inverse quantization for audio
DE602004030594D1 (en) * 2003-10-07 2011-01-27 Panasonic Corp METHOD OF DECIDING THE TIME LIMIT FOR THE CODING OF THE SPECTRO-CASE AND FREQUENCY RESOLUTION
CN100507485C (en) * 2003-10-23 2009-07-01 松下电器产业株式会社 Spectrum coding apparatus, spectrum decoding apparatus, acoustic signal transmission apparatus, acoustic signal reception apparatus and methods thereof
CA2457988A1 (en) * 2004-02-18 2005-08-18 Voiceage Corporation Methods and devices for audio compression based on acelp/tcx coding and multi-rate lattice vector quantization
EP1914722B1 (en) * 2004-03-01 2009-04-29 Dolby Laboratories Licensing Corporation Multichannel audio decoding
FI119533B (en) * 2004-04-15 2008-12-15 Nokia Corp Coding of audio signals
US20070147518A1 (en) * 2005-02-18 2007-06-28 Bruno Bessette Methods and devices for low-frequency emphasis during audio compression based on ACELP/TCX
WO2006103488A1 (en) * 2005-03-30 2006-10-05 Nokia Corporation Source coding and/or decoding
AU2006232364B2 (en) * 2005-04-01 2010-11-25 Qualcomm Incorporated Systems, methods, and apparatus for wideband speech coding
TWI317933B (en) * 2005-04-22 2009-12-01 Qualcomm Inc Methods, data storage medium,apparatus of signal processing,and cellular telephone including the same
FR2888699A1 (en) * 2005-07-13 2007-01-19 France Telecom HIERACHIC ENCODING / DECODING DEVICE
US7974713B2 (en) * 2005-10-12 2011-07-05 Fraunhofer-Gesellschaft Zur Foerderung Der Angewandten Forschung E.V. Temporal and spatial shaping of multi-channel audio signals
US8332216B2 (en) * 2006-01-12 2012-12-11 Stmicroelectronics Asia Pacific Pte., Ltd. System and method for low power stereo perceptual audio coding using adaptive masking threshold
US7831434B2 (en) * 2006-01-20 2010-11-09 Microsoft Corporation Complex-transform channel coding with extended-band frequency coding
US8260620B2 (en) * 2006-02-14 2012-09-04 France Telecom Device for perceptual weighting in audio encoding/decoding
US20080004883A1 (en) * 2006-06-30 2008-01-03 Nokia Corporation Scalable audio coding
US8032371B2 (en) * 2006-07-28 2011-10-04 Apple Inc. Determining scale factor values in encoding audio data with AAC
US8260609B2 (en) * 2006-07-31 2012-09-04 Qualcomm Incorporated Systems, methods, and apparatus for wideband encoding and decoding of inactive frames
US9454974B2 (en) * 2006-07-31 2016-09-27 Qualcomm Incorporated Systems, methods, and apparatus for gain factor limiting
CN101140759B (en) * 2006-09-08 2010-05-12 华为技术有限公司 Band-width spreading method and system for voice or audio signal
KR101565919B1 (en) 2006-11-17 2015-11-05 삼성전자주식회사 Method and apparatus for encoding and decoding high frequency signal
KR100905585B1 (en) * 2007-03-02 2009-07-02 삼성전자주식회사 Method and apparatus for controling bandwidth extension of vocal signal
US8392198B1 (en) * 2007-04-03 2013-03-05 Arizona Board Of Regents For And On Behalf Of Arizona State University Split-band speech compression based on loudness estimation
RU2439721C2 (en) * 2007-06-11 2012-01-10 Фраунхофер-Гезелльшафт цур Фёрдерунг дер ангевандтен Audiocoder for coding of audio signal comprising pulse-like and stationary components, methods of coding, decoder, method of decoding and coded audio signal
KR101373004B1 (en) * 2007-10-30 2014-03-26 삼성전자주식회사 Apparatus and method for encoding and decoding high frequency signal
US8515767B2 (en) * 2007-11-04 2013-08-20 Qualcomm Incorporated Technique for encoding/decoding of codebook indices for quantized MDCT spectrum in scalable speech and audio codecs
US20090201983A1 (en) * 2008-02-07 2009-08-13 Motorola, Inc. Method and apparatus for estimating high-band energy in a bandwidth extension system
CN101281748B (en) * 2008-05-14 2011-06-15 武汉大学 Method for filling opening son (sub) tape using encoding index as well as method for generating encoding index
CA2729752C (en) * 2008-07-10 2018-06-05 Voiceage Corporation Multi-reference lpc filter quantization and inverse quantization device and method
US8577673B2 (en) * 2008-09-15 2013-11-05 Huawei Technologies Co., Ltd. CELP post-processing for music signals
US8463599B2 (en) * 2009-02-04 2013-06-11 Motorola Mobility Llc Bandwidth extension method and apparatus for a modified discrete cosine transform audio coder
JP4932917B2 (en) * 2009-04-03 2012-05-16 株式会社エヌ・ティ・ティ・ドコモ Speech decoding apparatus, speech decoding method, and speech decoding program
US8571231B2 (en) * 2009-10-01 2013-10-29 Qualcomm Incorporated Suppressing noise in an audio signal
KR101411759B1 (en) * 2009-10-20 2014-06-25 프라운호퍼 게젤샤프트 쭈르 푀르데룽 데어 안겐반텐 포르슝 에. 베. Audio signal encoder, audio signal decoder, method for encoding or decoding an audio signal using an aliasing-cancellation
CA2683983A1 (en) 2009-10-21 2011-04-21 Carbon Solutions Inc. Stabilization and remote recovery of acid gas fractions from sour wellsite gas
US8484020B2 (en) * 2009-10-23 2013-07-09 Qualcomm Incorporated Determining an upperband signal from a narrowband signal
CN102044250B (en) * 2009-10-23 2012-06-27 华为技术有限公司 Band spreading method and apparatus
US8380524B2 (en) * 2009-11-26 2013-02-19 Research In Motion Limited Rate-distortion optimization for advanced audio coding
US8455888B2 (en) * 2010-05-20 2013-06-04 Industrial Technology Research Institute Light emitting diode module, and light emitting diode lamp
US9294060B2 (en) * 2010-05-25 2016-03-22 Nokia Technologies Oy Bandwidth extender
US8600737B2 (en) * 2010-06-01 2013-12-03 Qualcomm Incorporated Systems, methods, apparatus, and computer program products for wideband speech coding
US8862465B2 (en) * 2010-09-17 2014-10-14 Qualcomm Incorporated Determining pitch cycle energy and scaling an excitation signal
US8924200B2 (en) * 2010-10-15 2014-12-30 Motorola Mobility Llc Audio signal bandwidth extension in CELP-based speech coder
CN103035248B (en) * 2011-10-08 2015-01-21 华为技术有限公司 Encoding method and device for audio signals
EP3089164A1 (en) * 2011-11-02 2016-11-02 Telefonaktiebolaget LM Ericsson (publ) Generation of a high band extension of a bandwidth extended audio signal
EP2774148B1 (en) * 2011-11-03 2014-12-24 Telefonaktiebolaget LM Ericsson (PUBL) Bandwidth extension of audio signals
US8909539B2 (en) * 2011-12-07 2014-12-09 Gwangju Institute Of Science And Technology Method and device for extending bandwidth of speech signal
CN102930872A (en) * 2012-11-05 2013-02-13 深圳广晟信源技术有限公司 Method and device for postprocessing pitch enhancement in broadband speech decoding
ES2924427T3 (en) * 2013-01-29 2022-10-06 Fraunhofer Ges Forschung Decoder for generating a frequency-enhanced audio signal, decoding method, encoder for generating an encoded signal, and encoding method using compact selection side information
FR3008533A1 (en) * 2013-07-12 2015-01-16 Orange OPTIMIZED SCALE FACTOR FOR FREQUENCY BAND EXTENSION IN AUDIO FREQUENCY SIGNAL DECODER
CN104517610B (en) * 2013-09-26 2018-03-06 华为技术有限公司 The method and device of bandspreading
CN104517611B (en) * 2013-09-26 2016-05-25 华为技术有限公司 A kind of high-frequency excitation signal Forecasting Methodology and device
US10163447B2 (en) * 2013-12-16 2018-12-25 Qualcomm Incorporated High-band signal modeling
US9542955B2 (en) * 2014-03-31 2017-01-10 Qualcomm Incorporated High-band signal coding using multiple sub-bands
US9697843B2 (en) * 2014-04-30 2017-07-04 Qualcomm Incorporated High band excitation signal generation
JP2017145792A (en) 2016-02-19 2017-08-24 株式会社ケーヒン Sensor fixing structure at intake manifold
RU2636700C1 (en) * 2016-03-18 2017-11-27 Акционерное общество "Лаборатория Касперского" Method for eliminating vulnerabilities of devices having access to internet
TWI596952B (en) * 2016-03-21 2017-08-21 固昌通訊股份有限公司 In-ear earphone

Also Published As

Publication number Publication date
CN107527628A (en) 2017-12-29
CN107492385A (en) 2017-12-19
US10943594B2 (en) 2021-03-09
JP6515157B2 (en) 2019-05-15
JP2016528539A (en) 2016-09-15
RU2751104C2 (en) 2021-07-08
CN107527629B (en) 2022-01-04
RU2017144515A3 (en) 2021-04-19
KR102343019B1 (en) 2021-12-27
US20160203826A1 (en) 2016-07-14
KR102319881B1 (en) 2021-11-02
JP6515158B2 (en) 2019-05-15
KR102315639B1 (en) 2021-10-21
US10446163B2 (en) 2019-10-15
JP2017215601A (en) 2017-12-07
US20180018983A1 (en) 2018-01-18
CA2917795C (en) 2021-11-30
CA3108924A1 (en) 2015-01-15
RU2756435C2 (en) 2021-09-30
US20190371350A1 (en) 2019-12-05
CA3108921A1 (en) 2015-01-15
RU2017144518A (en) 2019-02-15
BR122017018557B1 (en) 2021-08-03
US10783895B2 (en) 2020-09-22
CN107527629A (en) 2017-12-29
RU2668058C2 (en) 2018-09-25
US10438600B2 (en) 2019-10-08
US20190385625A1 (en) 2019-12-19
MX2016000255A (en) 2016-04-28
WO2015004373A1 (en) 2015-01-15
US10943593B2 (en) 2021-03-09
BR122017018556B1 (en) 2022-03-29
US20180018982A1 (en) 2018-01-18
RU2017144515A (en) 2019-02-15
CA3108921C (en) 2024-01-30
CA2917795A1 (en) 2015-01-15
US20180082699A1 (en) 2018-03-22
KR20170103996A (en) 2017-09-13
CA3109028A1 (en) 2015-01-15
RU2756434C2 (en) 2021-09-30
RU2017144519A (en) 2019-02-15
US20190378527A1 (en) 2019-12-12
RU2016104466A3 (en) 2018-05-28
US10438599B2 (en) 2019-10-08
RU2017144518A3 (en) 2021-05-07
KR20170103995A (en) 2017-09-13
CN107492385B (en) 2022-02-11
RU2016104466A (en) 2017-08-18
EP3020043B1 (en) 2017-02-08
CN107527628B (en) 2021-03-30
US10672412B2 (en) 2020-06-02
BR112016000337B1 (en) 2021-02-23
JP2017215619A (en) 2017-12-07
US10354664B2 (en) 2019-07-16
KR20170103042A (en) 2017-09-12
CA3109028C (en) 2024-01-30
BR122017018553B1 (en) 2022-04-19
CN105378837A (en) 2016-03-02
CN105378837B (en) 2019-09-13
RU2017144519A3 (en) 2021-04-19
FR3008533A1 (en) 2015-01-16
JP2017215618A (en) 2017-12-07
KR20160030555A (en) 2016-03-18
MX354394B (en) 2018-02-23
KR102423081B1 (en) 2022-07-21
JP6487429B2 (en) 2019-03-20
JP6515147B2 (en) 2019-05-15
US20190385626A1 (en) 2019-12-19

Similar Documents

Publication Publication Date Title
EP3020043B1 (en) Optimized scale factor for frequency band extension in an audiofrequency signal decoder
EP3014611B1 (en) Improved frequency band extension in an audio signal decoder
EP3330966B1 (en) Improved frequency band extension in an audio frequency signal decoder
WO2014009657A1 (en) Processing method for improving the quality of audio frequency signals

Legal Events

Date Code Title Description
PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

Free format text: ORIGINAL CODE: 0009012

17P Request for examination filed

Effective date: 20160210

AK Designated contracting states

Kind code of ref document: A1

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

AX Request for extension of the european patent

Extension state: BA ME

DAX Request for extension of the european patent (deleted)
RIN1 Information on inventor provided before grant (corrected)

Inventor name: RAGOT, STEPHANE

Inventor name: KANIEWSKA, MAGDALENA

GRAP Despatch of communication of intention to grant a patent

Free format text: ORIGINAL CODE: EPIDOSNIGR1

INTG Intention to grant announced

Effective date: 20161123

GRAS Grant fee paid

Free format text: ORIGINAL CODE: EPIDOSNIGR3

GRAA (expected) grant

Free format text: ORIGINAL CODE: 0009210

AK Designated contracting states

Kind code of ref document: B1

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

REG Reference to a national code

Ref country code: GB

Ref legal event code: FG4D

Free format text: NOT ENGLISH

REG Reference to a national code

Ref country code: CH

Ref legal event code: EP

Ref country code: AT

Ref legal event code: REF

Ref document number: 867179

Country of ref document: AT

Kind code of ref document: T

Effective date: 20170215

REG Reference to a national code

Ref country code: IE

Ref legal event code: FG4D

Free format text: LANGUAGE OF EP DOCUMENT: FRENCH

REG Reference to a national code

Ref country code: DE

Ref legal event code: R096

Ref document number: 602014006746

Country of ref document: DE

REG Reference to a national code

Ref country code: LT

Ref legal event code: MG4D

REG Reference to a national code

Ref country code: NL

Ref legal event code: MP

Effective date: 20170208

REG Reference to a national code

Ref country code: FR

Ref legal event code: PLFP

Year of fee payment: 4

REG Reference to a national code

Ref country code: AT

Ref legal event code: MK05

Ref document number: 867179

Country of ref document: AT

Kind code of ref document: T

Effective date: 20170208

REG Reference to a national code

Ref country code: DE

Ref legal event code: R081

Ref document number: 602014006746

Country of ref document: DE

Owner name: KONINKLIJKE PHILIPS N.V., NL

Free format text: FORMER OWNER: ORANGE, PARIS, FR

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: GR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170509

Ref country code: FI

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170208

Ref country code: LT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170208

Ref country code: HR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170208

Ref country code: NO

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170508

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: PT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170608

Ref country code: LV

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170208

Ref country code: SE

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170208

Ref country code: NL

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20170208

Ref country code: RS

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170208

Ref country code: ES

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170208

Ref country code: BG

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170508

Ref country code: AT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170208

REG Reference to a national code

Ref country code: GB

Ref legal event code: 732E

Free format text: REGISTERED BETWEEN 20170824 AND 20170830

RAP2 Party data changed (patent owner data changed or rights of a patent transferred)

Owner name: KONINKLIJKE PHILIPS N.V.

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: IT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170208

Ref country code: SK

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170208

Ref country code: EE

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170208

Ref country code: RO

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170208

Ref country code: CZ

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170208

REG Reference to a national code

Ref country code: DE

Ref legal event code: R097

Ref document number: 602014006746

Country of ref document: DE

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: PL

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170208

Ref country code: SM

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170208

Ref country code: DK

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170208

PLBE No opposition filed within time limit

Free format text: ORIGINAL CODE: 0009261

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: NO OPPOSITION FILED WITHIN TIME LIMIT

26N No opposition filed

Effective date: 20171109

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: SI

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170208

REG Reference to a national code

Ref country code: CH

Ref legal event code: PL

REG Reference to a national code

Ref country code: IE

Ref legal event code: MM4A

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: LI

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20170731

Ref country code: CH

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20170731

Ref country code: IE

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20170704

REG Reference to a national code

Ref country code: BE

Ref legal event code: MM

Effective date: 20170731

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: LU

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20170704

REG Reference to a national code

Ref country code: FR

Ref legal event code: PLFP

Year of fee payment: 5

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: BE

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20170731

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: MT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170208

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: MC

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170208

Ref country code: HU

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT; INVALID AB INITIO

Effective date: 20140704

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: CY

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170208

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: MK

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170208

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: TR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170208

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: AL

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170208

Ref country code: IS

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170608

REG Reference to a national code

Ref country code: DE

Ref legal event code: R082

Ref document number: 602014006746

Country of ref document: DE

Representative=s name: EISENFUEHR SPEISER PATENTANWAELTE RECHTSANWAEL, DE

REG Reference to a national code

Ref country code: DE

Ref legal event code: R039

Ref document number: 602014006746

Country of ref document: DE

Ref country code: DE

Ref legal event code: R008

Ref document number: 602014006746

Country of ref document: DE

REG Reference to a national code

Ref country code: GB

Ref legal event code: S72Z

Free format text: COUNTERCLAIM LODGED; COUNTERCLAIM FOR REVOCATION LODGED AT THE PATENTS COURT ON 30 SEPTEMBER 2022 (HP-2020-000010)

P01 Opt-out of the competence of the unified patent court (upc) registered

Effective date: 20230527

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: GB

Payment date: 20230725

Year of fee payment: 10

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: FR

Payment date: 20230725

Year of fee payment: 10

Ref country code: DE

Payment date: 20230726

Year of fee payment: 10

REG Reference to a national code

Ref country code: DE

Ref legal event code: R040

Ref document number: 602014006746

Country of ref document: DE