Method and device for coding an audio signal by “forward” and “backward” LPC analysis

[0002] The invention involves a procedure and a device for coding an audio-frequency signal, such as a speech signal, by means of “forward” and “backward” LPC analysis.

[0003] At present, the aim of coding techniques for audio-frequency signals, particularly speech signals, is to allow for the transmission of these signals in digital form, within the conditions of reduction of the transmission output, in order, particularly, to ensure a management adapted to the networks for transmitting these signals, taking into consideration the considerable growth in transactions between users.

[0004] Of the coding techniques used, that designated by LPC analysis, Linear Predictive Coding in English, consists of carrying out a linear prediction of the audio-frequency signal to be encoded, the coding being carried out temporarily by means of a linear filtering prediction applied to the successive blocks of this signal.

[0005] Of the aforementioned techniques, that known as CELP coding, Code Excited Linear Prediction, is the most widespread and provides some of the best performance. Other techniques, such as the technique designated by MP-LPC, Multi Pulse Linear Predictive Coding, or the VSELP technique, Vector Sum Excited Linear Prediction in English, are relatively similar to CELP coding.

[0006] The aforementioned coding techniques are known as “analysis by synthesis”. They have enabled in particular, for audio-frequency signals belonging to the telephonic frequency bandwidth, the transmission output of these signals to be reduced from 64 kb/s (MIC coding) to 16 kb/s with the help of the CELP coding technique and even to 8 kb/s where these encoders use the most recent developments of this coding technique, without any perceptible reduction in the quality of the voice reconstituted after transmission and decoding.

[0007] A particularly important area of application for these coding techniques is, in particular, that of mobile telephony. Within this area of application, the necessary limitation of the frequency bandwidth granted to each mobile-telephony operator and the extremely rapid increase in the number of subscribers makes necessary the corresponding reduction of the coding output, while user demands in terms of speech quality continue to grow. Other areas of application of these coding techniques concern, for example, the storage of digital data which represent these signals on memory supports, high-quality telephony for video or audio conference applications, multimedia or digital transmissions via satellite.

[0008] The linear prediction filters used in the aforementioned techniques are obtained with the help of an analysis module called “LPC analysis” operating on successive digital signal blocks. These filters are capable, according to the order of analysis, that is, according to the number of filter coefficients, of modeling more or less reliably the contours of the spectrum of frequencies of the signal to be coded. In the case of a speech signal, these contours are called formants.

[0009] However, for good quality coding, required by most current applications, the filter thus defined is not sufficient for perfectly modeling the signal. It is therefore essential to code the residue of the linear prediction. One such operating mode relating to linear prediction residue is particularly used by the coding technique, LD-CELP, Low Delay CELP in English, previously mentioned in the description. In this case, the residual signal is modeled by a waveform taken from a stochastic codepage and multiplied by a gain value. The MP-LPC coding technique, for example, models this residue with the help of variable position pulses modified by respective gain values, whereas the VSELP coding technique carries out this modeling by means of a linear combination of pulse vectors taken from appropriate lists.

[0010] An explanatory recap of the operating method of LPC analysis and especially “backward” LPC analysis and “forward” LPC analysis will be given below.

[0011] The general envelope of the frequency spectrum is modeled by means of a short-term synthesis filter, constituting the LPC filter, the coefficients of which are modeled by means of a linear prediction of the speech signal to be coded. This LPC filter, an autoregressive filter, has a transfer function of the form, equation (1): $A\left(z\right)=1-\sum _{i=1}^p{a}_i{z}^{-i}$

[0012] where p designates the name of coefficients, ai of the filter and the order of the linear prediction applied, z designating the transformed variable z of the space of the frequencies.

[0013] One method of evaluating the coefficients a_iconsists of applying a criterion of minimization of the energy of the error prediction signal of the speech signal over the analysis length of this latter.

[0014] The analysis length for a digital speech signal formed of successive samples is, in practical terms, a number N of these samples, constituting a coding frame. The energy of the error prediction signal thus confirms equation (2): $\mathrm{Ep}=\sum _{n=1}^N{\left(s\left(n\right)-\sum _{i=1}^p{a}_i·s\left(n-1\right)\right)}^2$

[0015] where s(n) designates the sample of row n in the frame of N samples.

[0016] In a block-by-block coding process, the coding frame can be advantageously divided into several subframes or adjacent LPC blocks. The analysis length N then exceeds the length of each block in order to make it possible to take into account a certain number of past or, if applicable, future samples, by means of and at the cost of delaying the appropriate coding.

[0017] The analysis is called “forward” LPC when the LPC analysis process is carried out on the block of the current frame of the speech signal to be coded, with the coding taking place at encoder level “in real time”, that is, during the block of the current frame, with the only processing delay introduced by the calculation of the filter coefficients. This analysis involves transmitting the calculated values of the filter coefficients to the decoder.

[0018] “Backward” LPC analysis, used in the LD-CELP encoder at 16 kb/s is the object of the standard UIT-T G728. This analysis technique consists of carrying out the LPC analysis not on the block of the current frame of the speech signal to be coded, but on the synthesis signal. It is understood that this LPC analysis is actually performed on the synthesis signal of the block preceding the current block, as this signal is available simultaneously at encoder and decoder level. This simultaneous operation in the encoder and decoder thus makes it possible to avoid transmitting from the encoder to the decoder the value obtained in the encoder of the LPC filter coefficients. For this reason, “backward” LPC analysis makes it possible to free up transmission output and the output thus freed can be used, for example to enrich the excitation codepages in the case of CELP coding. “Backward” LPC analysis furthermore allows an increase in the order of analysis; the number of LPC filter coefficients may be as much as 50 in the case of an LD-CELP encoder, compared to 10 coefficients for most encoders using “forward” LPC analysis.

[0019] Thus, correct operation of “backward” LPC analysis requires the following conditions:

[0020] good quality synthesis signal, very close to the speech signal to be coded, which involves a sufficiently high coding output, higher than 13 kb/s, taking into account the quality of current CELP encoders;

[0021] reduced frame and block length due to the delay of one block between the analyzed signal and the signal to be coded. The length of the frame and block should therefore be low in comparison to the mean stationary time of the speech signal to be coded;

[0022] reliability of the transmission and conservation of the integrity of the data transmitted between the encoder and the decoder, by introducing few transmission errors. As soon as the synthesis signals differ significantly from the speech signal to be coded, the encoder and decoder cease to calculate the same filter and large divergences may occur, without being able to return to a noticeable similarity of the filters calculated in the encoder or decoder.

[0023] Due to the respective advantages and disadvantages of the aforementioned “backward” and “forward” types of LPC analysis, one technique consisting of selectively associating “backward” and “forward” LPC analysis was proposed in the article titled “Dual Rate Low Delay CELP Coding (8 kbits/s/16 kbits/s) using a Mixed Backward/Forward Adaptive LPC Prediction”, published by S. PROUST, C. LAMBLIN and D. MASSALOUX, Proc. IEEE Workshop Speech Co. Telecomm., September. 1995, pp 37-38.

[0024] The conditions mentioned above, regarding the correct functioning of “backward” LPC analysis, show that this type of analysis alone presents the limitations mentioned when operating at transmission outputs appreciably below 16 kb/s. Besides the reduction in the quality of the synthesis signal, which reduces the performance of the LPC filter, it is very often necessary, in order to reduce the transmission output, to operate with a greater LPC frame length, of the order of 10 to 30 ms. It can therefore be seen that, under these conditions, the degradation occurs especially during transitions of the frequency spectrum and, more generally, in the not so stationary areas, since for generally very stationary signals, such as music signals, “backward” LPC analysis holds a considerable advantage over “forward” LPC analysis.

[0025] The association of the two aforementioned types of LPC analysis aims to reduce these disadvantages and increase the advantages inherent in each one:

[0026] “forward” LPC analysis for the coding of the transitions and the non-stationary areas;

[0027] “backward” LPC analysis, to a greater extent, for the coding of the stationary areas.

[0028] Furthermore, the introduction of LPC frames coded by “forward” LPC analysis into LPC frames coded by “backward” analysis allows the encoder and decoder to re-converge towards the same synthesis signal in the case of a transmission error and therefore offers far greater error protection than coding by “backward” LPC analysis alone.

[0029] In general, the above-mentioned mixed “forward”-“backward” LPC analysis consists of carrying out two LPC analyses, a “forward” LPC analysis of the speech signal or audio frequency to be coded and a “backward” LPC analysis of the synthesis signal.

[0030] Two filters are calculated for each LPC block, these filters being designated by “forward” LPC filter and “backward” LPC filter, respectively. A procedure of choosing the filter applied to the LPC block, depending on whether the signal is stationary, is therefore applied. This procedure requires two different criteria:

[0031] a first criterion based on the prediction gains of the filters;

[0032] a second criterion based on a distance parameter between the “forward” LPC filters calculated successively.

[0033] For each of these two criteria, the threshold values are established.

[0034] First Criterion:

[0035] The choice of “backward” LPC filter is made if the distance between the prediction gain of the “backward” and “forward” LPC filters is greater than a first threshold value.

[0036] Second Criterion:

[0037] For a current analysis in “backward” LPC analysis mode, prohibition of switching from “backward” LPC analysis mode to “forward” LPC analysis mode if the distance calculated on the vectors of the parameters representing two consecutive “forward” LPC filters is lower than a second threshold value, a distance which is too small characterizing a more or less stationary area, for which reason it is appropriate to avoid changing the LPC analysis mode. The calculated distance is a Euclidean distance between the spectral lines of the speech or audio-frequency signal to be coded.

[0038] A more detailed description of the aforementioned mixed LPC analysis method can be found in the article published by S. PROUST, C. LAMBLIN and D. MASSALOUX, mentioned above.

[0039] In-depth studies on the above-mentioned mixed analysis operating method have shown the following important disadvantages:

[0040] for certain signals, the prediction gain values of the “forward” and “backward” LPC filters may oscillate above and below the first threshold value. This phenomenon leads to sudden and frequent changes from “backward” LPC filter to “forward” LPC filter or vice versa. The discontinuity of filtering thus introduced constitutes a source of considerable degradation of the synthesis signal and is not, most of the time, linked to the real spectral modifications of the speech or audio-frequency signal to be coded;

[0041] the optimal value of the first threshold which should be established varies considerably according to whether the signal to be coded is stationary, more so when the coding output is low. For a coding delay corresponding to an LPC frame of 10 to 30 ms, or when the transmission output falls, there is a clear divergence between the coding mode of musical signals and speech signals; “forward” LPC analysis is mainly used.

[0042] Since music signals are quite stationary, “backward” LPC analysis is used even for long LPC frames. In the case of speech signals, however, the highly stationary areas have a very short duration and their passage in “backward” LPC analysis mode is therefore brief, thus leading to unwanted filter transitions which reduce the quality of the coding. The encoder can thus no longer correct the phenomena generated by the discontinuity introduced by the switching of the filters.

[0043] The LPC filter which gives the best subjective quality and which therefore best models the spectrum of the signal to be coded is not always that which has the best prediction gain. Certain switchings from one mode of LPC analysis to another, linked to an instantaneous decision, are therefore useless.

[0044] The object of the present invention is to resolve the aforementioned disadvantages by employing a procedure and device for coding a digital audio-frequency signal by means of specific “forward” and “backward” LPC analysis.

[0045] Another object of the present invention is also to employ a process for dynamically adapting the function of choice between “forward” LPC analysis and “backward” LPC analysis according to how stationary the signal to be coded is.

[0046] A further object of the present invention is also to employ a process for dynamically adapting the aforementioned choice function on the basis of discrimination between highly stationary signals, such as music or background noise, and other signals, such as speech, in order to allow the most appropriate code processing by “backward” LPC analysis and “forward” LPC analysis, respectively.

[0047] A further object of the present invention is, once the aforementioned most appropriate choice of coding has been made, for a signal to be coded of a given type or with given characteristics, to prevent any sudden switching to the LPC analysis mode not chosen and, therefore, to prevent the appearance of transitions from “forward” LPC filters to “backward” LPC filters and vice versa, which tend to reduce the quality of the reproduced synthesis signal.

[0048] A further object of the present invention is to employ a dynamic adaptation process of the aforementioned choice function by which the change in the LPC analysis mode corresponds reliably to a change in the stationarity of the signal to be coded, thus having a far lower chance of being linked to a simple crossover effect of the first and second threshold values.

[0049] The method and device for coding a digital audio-frequency signal, which are the object of the present invention, employ a double analysis based on the criterion of choice between “forward” and “backward” LPC analysis, respectively, to create a transmitted coded signal consisting of LPC filtering parameters accompanied by analysis decision information and a non-transmitted coding residue signal. The digital audio-frequency signal is subdivided into frames, succession of blocks of a determined number of samples, and the coding of this digital audio-frequency signal is carried out on this signal using a “forward” LPC filter for the non-stationary areas and a synthesis signal, respectively. This synthesis signal is obtained from the coding residue signal, using “backward” LPC filtering for the stationary areas.

[0050] They are notable insofar as they consist of and allow for, respectively:

[0051] determining the degree of stationarity of the digital audio-frequency signal according to a stationarity parameter whose value is between a maximum stationarity value and a minimum stationarity value;

[0052] establishing, based on the stationarity parameter, an analysis choice value, based on a decision function;

[0053] applying the analysis choice value to the LPC filtering in order to code the digital audio-frequency signal by means of “forward” LPC filtering on the non-stationary areas of the digital audio-frequency and by means of “backward” LPC filtering on the stationary areas of the synthesis signal.

[0054] This operating method makes it possible to prioritize remaining in either the “forward” or “backward” LPC filtering mode, according to the degree of stationarity of the digital audio-frequency signal and to limit the number of switchings from one mode of filtering to another and vice versa.

[0055] The method and the device which are the object of the present invention have an application not only in the area of mobile telephony, but also in the sector of creation and reproduction of phonograms, satellite transmission and high-quality telephony for multimedia video or audio conference applications.

[0056] Understanding will be facilitated by reading the description and examining the design below, where:

[0057] FIG. 1 shows, in the form of a general flow chart, an explanatory diagram of the stages which allow the performance of the coding which is the object of the present invention;

[0058] FIG. 2ashows a general flow chart of the stages of calculating the stationarity parameter for each current LPC block;

[0059] FIG. 2bshows a particularly advantageous method of carrying out the essential stages of the calculation of the stationarity parameter, according to FIG. 2a;

[0060] FIG. 2cshows a detail of the execution of FIG. 2band, more particularly, a detail of the process of tuning the value of the intermediate stationarity parameter in order to obtain the stationarity parameter;

[0061] FIGS. 2dand 2eshow, respectively, a first and second example of the application of a tuning function, allowing for the calculation of a tuning value for the intermediate stationarity function according to the comparative values of the “forward” and “backward” LPC filter gain;

[0062] FIG. 2fshows as an explanatory example a flow chart of the stages making it possible to employ the decision function and the “forward” or “backward” LPC analysis choice value;

[0063] FIG. 3 shows, in the form of functional blocks, the general diagram of an encoder which makes it possible to code an audio-frequency signal according to the object of the present invention;

[0064] FIG. 4 shows, in the form of functional blocks, the general diagram of a decoder which makes it possible to decode an audio-frequency signal which has been coded by using an encoder as shown in FIG. 3.

[0065] A more detailed description of the method for coding a digital audio-frequency signal, employing a double analysis based on the criterion of choice between “forward” and “backward” LPC analysis, respectively, of a transmitted coded signal, which is the object of the present invention, is now given in connection with FIG. 1.

[0066] In general terms, it is shown that the transmitted coded signal, written as s_c_n(t), consists in part of the LPC filtering parameters accompanied by LPC analysis decision information. Furthermore, a non-transmitted coding residue signal, res_n(t), is available for performing the coding procedure.

[0067] The digital audio-frequency signal is subdivided into LPC frames, a succession of LPC blocks, each block, for the sake of convenience of the description, being written as Bn and having a determined number of samples, N.

[0068] One aspect of the coding procedure which is the object of the present invention consists of carrying out the aforementioned coding of the digital audio-frequency signal as described above using “forward” LPC filtering for the non-stationary areas and for a synthesis signal obtained from the coding residue signal using “backward” LPC filtering for the stationary areas.

[0069] A particularly notable aspect of the method which is the object of the present invention consists of, in order to establish the “forward” or “backward” LPC filter choice criteria for each current block of the succession of current blocks forming the current frame, as shown in FIG. 1, each current block, written B_n, being available in an initial stage 10, to determine in stage 11 the degree of stationarity of the digital audio-frequency signal, according to a stationarity parameter, written STAT(n). This stationarity parameter presents a digital value between a maximum stationarity value, written STAT_M, and a minimum stationarity value, written STAT_m.

[0070] By way of convention and without prejudice to the degree of generality of the coding procedure which is the object of the present invention, the stationarity parameter presents the maximum value STAT_Mfor an extremely stationary signal, whereas this stationarity parameter presents the minimum value STAT_mfor a highly non-stationary signal.

[0071] After the aforementioned stage 11, the coding method which is the object of the present invention consist of establishing, in stage 12, using the stationarity parameter STAT(n), an LPC analysis choice value. This analysis choice value corresponds, logically, to either the “forward” LPC analysis choice or the “backward” LPC analysis choice. The value of the choice of analysis is written d_n(n) and is obtained from a specific decision function, written D_n.

[0072] The aforementioned stage 12 is then followed by a test stage 13 which allows the application of the analysis choice value d_n(n), represented by C, to the LPC filtering in order to carry out the coding of the digital audio-frequency signal by means of “forward” LPC filtering for the non-stationary areas of the digital audio-frequency signal and by means of “backward” LPC filtering for the stationary areas of the synthesis signal.

[0073] The execution of the decision function D_nand the aforementioned analysis choice values d_n(n) form a particularly advantageous aspect of the coding procedure which is the object of the present invention, as they make it possible to prioritize remaining in one of the LPC filtering modes, either “forward” or “backward”, according to the degree of stationarity of the audio-frequency signal and to limit the number of switchings from one to other of the filtering modes, and vice versa.

[0074] In general terms, it is mentioned that the decision function executed in stage 12 and indicated as D_nis an adaptive function, updated for each current block B_nfrom the stationarity parameter.

[0075] Updating the adaptive function makes it possible to prioritize remaining in one of the LPC filtering modes, either “forward” or “backward”, according to the degree of stationarity of the digital audio-frequency signal and to hence limit the number of switchings from one to other of the filtering modes, and vice versa.

[0076] More specifically, the analysis choice value d_n(n) established according to the aforementioned decision function D_ncorresponds to a priority value of the LPC filtering mode, either “forward” or “backward”, and to another priority value representing in fact a value of absence of priority for returning to the “backward” or “forward” LPC filtering mode.

[0077] As a priority value for the LPC filtering mode, it is mentioned that the analysis choice value d_n(n) can, for example, correspond to a logical value, the true value of this logical value, value 1, for example, corresponding to a choice of “backward” LPC filtering, whereas the complementary value of this true value, the value zero, corresponds to a choice of “forward” LPC filtering. It can thus be seen that the test function in stage 13 can be summarized as a test value of the logical value of the aforementioned analysis choice value, to ensure in stage 14 “backward” LPC filtering for the stationary areas of the signal to be coded or “forward” LPC filtering in stage 15 for the non-stationary areas, the aforementioned stages 14 and 15 being thus followed by stages 14aand 15aback to the next block, written as B_n+1for n=n+1.

[0078] Although the analysis choice value d_n(n) is represented by a logical value, it is understood that this logical value may be associated with a value of priority and probability of the mode of filtering specifically established by the decision function D_n. It can particularly be seen that this probability value may correspond, for each current block B_n, to the true logical value for a range of probability values between zero and 1 for “backward” LPC filtering while the complementary value, the logical value zero, for example, may correspond to the complement of the aforementioned range of probability values between zero and 1 for the first aforementioned range. This probability depends on a number of successive filtering decisions within the same filtering mode.

[0079] The operating mode of the decision function D_nmakes it possible in fact to associate with the logical variable d_n(n) the filtering mode priority and is adaptive over time for each current block B_n.

[0080] In general terms, it is mentioned that the aim of adapting the decision function D_nis to progressively prioritize the “backward” LPC filtering mode or, in contrast, the “forward” LPC filtering mode, whichever works better, taking into account the overall stationarity of the signal to be coded, in order to avoid as far as possible any unnecessary switching from one mode of filtering to another.

[0081] More specifically:

[0082] the more stationary the signal to be coded, the more the decision function D_nprioritizes “backward” LPC analysis, limiting as far as possible switching to “forward” LPC filtering mode;

[0083] in contrast, the less stationary the signal to be coded, the more the decision function D_nprioritizes “forward” LPC analysis, limiting as far as possible any switching to “backward” LPC filtering mode.

[0084] A more detailed description of the execution of the specific decision function which makes it possible to adapt this decision function, according to the value of the stationarity parameter STAT(n), is given later in the description.

[0085] A method of preferential calculation of the stationarity parameter STAT(n) relating to each current LPC block B_nis now given and described in connection with FIG. 2a.

[0086] According to the aforementioned figure, stage 11, consisting of determining the degree of stationarity of each current block B_nof the digital audio-frequency signal consists, starting with an arbitrary initial value of the stationarity parameter, as shown in stage 110FIG. 2a, this arbitrary value being written STAT(O), of calculating in stage 111 for this current block B_n, an intermediate stationarity parameter value, written STAT*(n), as a function of a determined number of successive analysis choice values, these LPC analysis choice values, written d_n−1(n−1), . . . , to d_n−p(n−p), being obtained for different successive blocks prior to the current block B_nof the succession of LPC blocks and a value of the stationarity parameter of the block preceding the current block, this stationarity value being written STAT(n−1). In stage 111 shown in FIG. 2a, the function of the determined number of previous analysis choice values is given in relation to these previous values, written d_n−1(n−1) to d_n−p(n−p). The initial arbitrary value for the stationarity parameter STAT(O) can, for example, be the same as the mean value between the maximum value and the minimum value of the stationarity parameter mentioned above in the description, STAT_Mand STAT_m.

[0087] The aforementioned stage 111 is then followed by stage 112 which consists of tuning the intermediate stationarity parameter value according to the value of the prediction gains of the “forward” and “backward” LPC filters or analysis mode of the frame preceding the current frame. In stage 112 of FIG. 2a, the aforementioned function is written g(STAT*(n), Gpf, Gpb) where Gpf is the prediction gain of the “forward” LPC filter and Gpb is the prediction gain of the “backward” LPC filter for the frame preceding the current frame. In stage 112, that is, the stage which consists of tuning the intermediate stationarity parameter value, the stationarity parameter value STAT(n) of the current LPC block B_nis given the value, equation (3):

[0088] STAT(n)=g(STAT*(n), Gpf, Gpb)

[0089] corresponding to the tuned intermediate stationarity parameter value.

[0090] A more detailed description of the calculation stage 111 of the intermediate stationarity parameter STAT*(n) and of stage 112 which consists of tuning this parameter value is now given in connection with FIG. 2b.

[0091] According to the aforementioned figure, stage 111, starting with an initialization stage 1110 in which the value of the stationarity parameter STAT(n−1) and the analysis choice value d_n−1(n−1) relating to the LPC block B_n−1prior to the current block B_nis available, consists to carry out in stage 1111 a stage which consists of discriminating the “forward” or “backward” LPC analysis mode of the block B_n−1preceding the current block B_n. This discrimination stage 1111, as shown in FIG. 2b, may consist of a test stage for the analysis choice value d_n−1(n−1) in relation to the symbolic value “fwd” or the logical value zero, corresponding to the complementary value of the true logical value.

[0092] On a negative response to the aforementioned test 1111, that is, for any block B_n−1preceding the current block B_n, analyzed in “backward” LPC analysis mode, the stage which calculates the intermediate stationarity parameter value consists, in stage 1113, of determining the number of previous frames consecutively analyzed in “backward” LPC analysis mode, written N_BWD; then, in stage 1114, it consists of comparing the superiority of the number of previous frames to an initial arbitrary value, written Na, representing a number of successive frames analyzed in “backward” LPC mode.

[0093] On a positive response to the superiority comparison of test 1114, the calculation stage consists of attributing in stage 1114b, to the intermediate stationarity parameter value STAT*(n), the value of the stationarity parameter of the block preceding the current block, STAT(n−1), increased by a determined value which depends on the first arbitrary value representing a number of successive analyzed frames, that is, the number of previous frames N_BWD, analyzed consecutively in “backward” LPC analysis mode. In stage 1114b, the determined value which depends on the first arbitrary value is written f_n(N_BWD). During the aforementioned stage, it can be seen that the intermediate stationarity parameter value STAT*(n) for the current LPC bloc B_nis thus increased in relation to the value corresponding to the same stationarity parameter for the preceding block B_n−1.

[0094] On a negative response to the superiority comparison in the comparison test 1114, the value of the stationarity parameter STAT(n−1) of the block preceding the current block B_n, is attributed, in stage 1114a, to the intermediate stationarity parameter value STAT*(n).

[0095] However for every preceding bloc B_n−1analyzed in “forward” LPC analysis mode, that is, on a positive response to test 1111, the stage for calculating the intermediate stationarity parameter, 111, as shown in FIG. 2b, consists of determining in stage 1112, on the test criterion, the occurrence of a transition from “backward” LPC analysis mode to “forward” LPC analysis mode between the block before the block before the current bloc B_n−1, of row n−2, that is, the existence of an LPC analysis choice value d_n−2(n−2)=symbolic value “bwd”, whose logical value is zero, as mentioned above. The positive response to test 1112 indicates the existence of such a transition from the “backward” analysis mode by the LPC block B_n−1preceding the block preceding the current block B_n−1, whereas a negative response to the aforementioned test 1112 indicates the absence of such a transition.

[0096] On a positive response to the aforementioned occurrence test 1112, the calculation stage 111 then consists of comparing using an inferiority comparison criterion, the number of previous aforementioned N_BWD frames with a second arbitrary value N_bwhich represents a number of frames successively analyzed in “backward” LPC mode preceding the bloc B_n−1preceding the current bloc.

[0097] On a positive response to the comparison performed in test 1118, this test is followed by stage 1118a, which consists of attributing to the intermediate stationarity parameter value STAT*(n) the stationarity parameter value of the block preceding the current block, STAT(n−1), reduced by a determined value which depends on the second arbitrary value N_b; this determined value is written f₂(N_BWD). It can be seen that in the attribution stage 1118a, the intermediate stationarity parameter value is thus reduced as a result.

[0098] However on a negative response to the inferiority comparison carried out in test 1118, stage 111 consists then in allocating, in a stage 1118b, to the value of the intermediate stationary parameter STAT*(n) the value of the stationarity parameter of the block preceding the current block, i.e. STAT(n−1).

[0099] In FIG. 2bit will be noticed that the allocation stages 1118aand 118bare then followed by a stage of replacing with zero the number of successive blocks processed in the “backward” LPC analysis mode, this stage of making to zero carrying the reference 1118cand enabling the updating the whole of the calculation process of the value of the intermediate stationarity parameter.

[0100] On a negative response to the comparison test 1112, no “forward” LPC transition analysis occurring, the value of the stationarity parameter STAT(n−1) of the preceding block B_n−1is attributed to the value of the intermediate stationarity parameter STAT*(n) in a stage 1119.

[0101] At the end of stage 111, the value of the intermediate stationarity parameter STAT*(n) is set for the current block B_n.

[0102] As far as the stage 112 consisting in tuning the value of the aforementioned intermediate stationarity parameter is concerned, it is noted, by reference to FIG. 2b, that it consists to advantage, of a stage 1120, in distinguishing the prediction gains of the “backward” LPC filtering and the “forward” LC filtering, these gain values being noted Gpb and Gpf respectively. It is understood that the aforementioned discrimination consists simply in memorizing and reading the gain values calculated for the respectively aforementioned “forward” and “backward” filtering. As well as the aforementioned gain values, the stage 1120 may consist of calculating the comparative value of the prediction gains, noted DGfb, as the difference or the ratio between the aforementioned “forward” and “backward” prediction gains.

[0103] As has been shown furthermore in FIG. 2b, the stage 112 of FIG. 2aincludes behind the aforementioned stage 1120 a stage 1121 consisting of modifying the value of the intermediate stationarity parameter STAT*(n) with a refining value ΔS, this refining value according to a particularly noteworthy characteristic of the method which is the object of the present invention being a function of the comparative value of the “forward” and “backward” LPC filtering prediction gains.

[0104] As a general rule, it is indicated that the function representative of the refining value ΔS is noted:

[0105] ΔS=f_r(Gpf, Gpb)

[0106] where Gpf and Gpb designate as previously the “forward” and “backward” LPC filtering prediction gains respectively.

[0107] As a general rule, it is indicated that the function f_r(GPf, Gpb) enabling the setting up of the refining value ΔS is a function respectively increasing and decreasing with this comparative value, according to the direction in which this comparative value is considered. When the comparative value designates the value of the “backward” LPC filtering gain comparative to the “forward” LPC filtering gain, this choice may be arbitrarily retained without any damage in the general nature of the method, the object of the invention, to the aforementioned comparative value DGfb, the function fr is then increasing. It is decreasing in the opposite case.

[0108] In other terms, the modification, by increasing or decreasing, the value of the intermediate stationarity parameter of the refining value ΔS is proportional to this comparative value of the gains. As a general rule, this modification is written STAT(n)=STAT*(n)+kΔS. In practice k is taken as equal to 1. In more specific terms, it is shown that the refining value ΔS increases in algebraic value when the gap between the “forward” and “backward” LPC filtering prediction gains increases, the function f_r(GPf, Gpb) being then an increasing function, whereas this refining value ΔS decreases in algebraic value when this same aforementioned gap decreases, the aforementioned gap being defined between the prediction gain of the LPC “backward” filtering and the prediction gain of the LPC “forward” filtering. In fact, this function is increasing or decreasing according to the definition of this gap.

[0109] Consequently, at the end of stage 1121 as shown in FIG. 2b, the value of the intermediate stationarity parameter STAT*(n) can then, for k=1, be corrected by the algebraic value of the aforementioned refining value ΔS in order to calculate the value of the stationarity parameter STAT(n).

[0110] Following stage 1121, the value of the stationarity parameter STAT(n) is set in stage 1122.

[0111] A more detailed description of the stage 1121 of FIG. 2bwill be now given in connection with FIG. 2cin a preferential version in which several test criteria are applied as much as to the refining value as to the values of the LPC “forward” and “backward” prediction gain in view of optimizing the calculation process of the stationarity parameter.

[0112] As is shown in the aforementioned FIG. 2c, the stage 1121 can consist of a first stage 1121aenabling the calculation of the refining value ΔS from the previously quoted function f_r(Gpf, Gpb). Different examples of useable functions will be given later in the description.

[0113] In the first place, the refining value ΔS is subject to a superiority comparison test with the value 0, in a stage 1121b, this comparison test enabling in fact to determine the increase of this refining value ΔS.

[0114] On a positive response to the aforementioned test 1121b, the refining value ΔS being positive and corresponding to an increase in the comparative value of the “forward” and “backward” LPC filtering prediction gains, the stage of increasing the value of the intermediate stationarity parameter from the refining value ΔS is moreover subjected to a superiority condition of the gain value of “backward” LPC filtering, in comparison with a first positive value determined in a superiority comparison stage of the value of the “backward” LPC filtering gain value Gpb in comparison with this first determined positive value, called S_i.

[0115] On a negative response to the aforementioned test 1121c, the value of the intermediate stationarity parameter STAT*(n) is attributed to the value of the stationarity parameter STAT(n) in a stage 1121g.

[0116] On a positive response to the aforementioned test 1121c, the increase of the value of the intermediate stationarity parameter of the refining value ΔS is furthermore subjected to an inferiority condition of the value of the intermediate stationarity parameter STAT*(n) in comparison with a second determined positive value STAT_irepresenting of course a stationarity value. This inferiority test condition is carried out in the stage 1121e.

[0117] On a negative response to the aforementioned test 1121e, the value of the intermediate stationarity parameter STAT*(n) in the aforementioned stage 1121gis attributed to the value of the intermediate stationarity parameter STAT(n).

[0118] On a positive response to the inferiority test condition 1121ethe value of the intermediate stationarity parameter STAT*(n) increased by the positive value ΔS of the refining value in the stage 1121iis attributed to the value of the intermediate stationarity parameter STAT(n).

[0119] In contrast, on a negative response to the aforementioned test 1121b, the refining value ΔS being negative, the reduction stage of the intermediate stationarity parameter with the refining value ΔS, this value being negative, is furthermore subject to an inferiority test condition of the “backward” LPC filtering gain value Gpb in comparison with a determined third positive value called S_din a comparison stage 1121d. This third determined positive value is of course representative of an LPC filtering gain value.

[0120] On a negative response to the aforementioned test 1121dthe value of the intermediate stationarity parameter STAT*(n) is attributed to the value of the stationarity parameter STAT(n) in the stage 1121g.

[0121] In contrast, on a positive response to the aforementioned test 1121d, the reduction stage of the value of the intermediate stationarity parameter with the refining value ΔS is furthermore subject to a superiority condition of the value of the intermediate stationarity parameter STAT*(n) in comparison with a fourth determined positive value, called STATd in a comparison test called 1121f. Of course, the fourth determined positive value is representative of a chosen stationarity parameter value.

[0122] On a negative response to the aforementioned test 1121f, the value of the intermediate stationarity parameter STAT*(n) is attributed to the stationarity parameter STAT(n) in the stage 1121g.

[0123] On a positive response to the aforementioned test 1121f, the value of the intermediate stationarity parameter STAT*(n) increased by the algebraic value of the refining value ΔS, negative, is attributed to the stationarity parameter STAT(n), the value of the intermediate stationarity parameter being thus reduced in order to set up the value of the stationarity parameter STAT(n) in the stage 1121h.

[0124] At the end of the stages 1121g, 1121hand 1121i, the stationarity parameter STAT(n) is thus set in the stage 1122 of FIG. 2b.

[0125] As regards the function f_r(Gpf, Gpb), it is shown that it may consist of a non linear function of the comparative value of the “forward” and “backward” LPC filtering gains in which the comparative value of the “forward” and “backward” LPC filtering prediction gains may themselves consist either in the ratio of, or in the difference of the “forward” and “backward” LPC filtering prediction gains. Other types of functions, such as linear functions, may be used.

[0126] A first example of the non linear function f_r(Gpf, Gpb) is shown in FIG. 2d.

[0127] In the version example of FIG. 2d, value pairs of the “backward” LPC filtering prediction gain Gpb in the ordinate and the “forward” LPC filtering gain Gpf enable allocating the positive refining values ΔS, ΔS>0 or negative ΔS<0 for a value of the ratio ρ=Gpb/Gpf corresponding to a respectively greater or lesser slope than that of the straight line ΔS=0.

[0128] In FIG. 2e, has been shown the case where the relative value of the “forward” and “backward” filtering prediction gains no longer correspond to the ratio of the gains ρ but to the difference of the aforementioned gains. In this case, the relative value of the “forward” and “backward” LPC filtering prediction gains can also be a non linear function enabling allocating to the refining value ΔS for the values of this difference corresponding to the value pairs Gpb, Gpf corresponding to the straight lines for which the abscissa origin is respectively less or greater, in algebraic value, than the abscissa origin of the straight line ΔS=0. In the case of FIG. 2e, the straight lines delimiting the zones as a function of the sign of the refining value ΔS are parallel to each other.

[0129] According to another particular aspect of the procedure which is the object of the invention, it is recommended furthermore that it is accepted not to adapt the stationarity index of the current block B_nduring the silence frames, when for example the audio frequency signal is constituted by a speech signal comprising silences. In such a case, the stage 1111 of the stage 111 shown in FIG. 2bcan be preceded by a stage 1111aconsisting, for each successive current block, in determining the mean energy of the audio frequency digital signal and comparing in this same stage, on inferiority comparison criterion, this mean energy with a determined threshold value representative of a silence frame. In FIG. 2b, this threshold value is called ENER_SIL. On a positive response to the aforementioned test, the value of the stationarity parameter of the preceding block STAT(n−1) in the allocation stage 1111bshown in FIG. 2bis attributed to the value of the stationarity parameter of the current block STAT(n). The stages 1111aand 1111bare, in the aforementioned figure, shown as a dotted line, because it is reserved for example to the coding of a speech signal.

[0130] A more detailed description of the implementation decision function D_nenabling the decision values d_n(n) to be obtained will be now given in connection with FIG. 2f. This description is given in a preferential version in which this decision function, being able to be compared with that which is described in the previously mentioned article by the description, published by S. PROUST, C. LAMBLIN and D. MASSALOUX, is however temporally adapted, according to the object of the present invention in order to obtain the successive choice analysis d_n(n) values.

[0131] Starting with a stage 120, for the current block B_n, in the first place a distance, called d_LPC, between the LPC filter of the current block and that of the preceding block B_n−1is calculated. This distance calculation is carried out for example by using the LSP frequency parameters as previously mentioned in the description relating to the procedure described in the aforementioned article.

[0132] It is noted:

[0133] the values of the thresholds S_PRED(n) and S_TRANS, S_STAT and G₁being reached in the criterion justified on the prediction gains of the “backward” and “forward” LPC filters;

[0134] the threshold values S_LSP_L and S_LSP_H being reached in the criterion justified on the distances between LSP frequency vectors representing two “forward” LPC filters comparative to two consecutive blocks B_n−1and B_n;

[0135] the prediction gain Gpf of the “forward” LPC filter;

[0136] the prediction gain Gpb of the “backward” filter; and

[0137] the prediction gain Gpi of the “forward” filter interpolated according to the method explained in the published article, previously mentioned in the description.

[0138] The criterion for establishing the decision function, in relation to FIG. 2f, is established in the manner below:

[0139] if the consecutive LPC filters are very stationary, i.e. for d_LPC<S_LSP_L, then, no switching of the “backward” LPC filtering with the “forward” LPC filtering is carried out if it is in the “backward” LPC filtering mode, on condition that the prediction gain of the “backward” LPC filter is greater than the prediction gain of the “forward” LPC filter reduced by a S_STAT value. It is mentioned that the S_STAT value is chosen so as to favor the choice of a “backward” LPC filter in the presence of a large stationarity of the spectrum measured by means of the distance d_LPC;

[0140] if the consecutive LPC filters have a significant transition, i.e. for d_LPC>S_LSP_H and if Gpf>Gpb−S_TRANS, then the chosen filtering mode is the “forward” LPC filtering, i.e. d_n(n)=0, symbolic value “fwd”, otherwise, d_n(n) is almost equal to 1, symbolic value “bwd”. It is mentioned that the value of S_TRANS is chosen so as to strongly favor the choice of the “forward” LPC filter in the presence of a spectrum transition measured by means of the distance d_LPC;

[0141] otherwise, in all other cases, if Gpb>Gpf-S_PRED and Gpi>Gpf-S_PRED, then, the LCP filter retained is the interpolated “backward” LPC filter, on condition that the gain of this latter and that of the pure “backward” LPC filter exceeds the threshold value G_ipreviously mentioned. If the condition on the values of the aforementioned prediction gain is not fulfilled, then, the “forward” LPC filtering is chosen.

[0142] In order to increase the number of transmitted “forward” LPC filters and thus to increase the strength of the coding system to the transmission errors, the “forward” LPC filtering mode may be chosen with advantage as soon as the energy signal to be coded E_n, i.e. the energy of the corresponding block B_n, becomes less than the value of the energy of a silence frame ENER_SIL, this value of energy corresponding to the minimum audible level.

[0143] The set of the conditions enabling the establishment of the decision function D_nand the obtaining of the corresponding chosen analysis values d_n(n), is illustrated in FIG. 2fwith temporal adaptation of the decision function D_n.

[0144] The value of the stationarity parameter STAT(n) can for example be located on a scale of 0, corresponding to the non-stationary STAT_nvalue, to 100, corresponding to the very stationary STAT_(n)value.

[0145] According to the value of the stationarity parameter STAT(n), the decision function D_nis modified by adaptation of the value of the thresholds.

[0146] The more the stationarity of the signal increases, the more the “backward” LPC filtering mode is favored: the thresholds S_PRED, S_LSP and S_LSP_H are increased.

[0147] As a non-limited example, the modification functions for each current LPC block B_nof the aforementioned threshold values have been shown:

[0148] S_PRED(n)=f_sPRED(STAT(n)) with the function f_sPREDincreasing with the value of STAT(N);

[0149] S_LSP_L(n)=f_SLPCL(STAT(n)) with the function f_SLPCLincreasing;

[0150] S_LSP_L(n)=f_sLPCR(STAT(n)) with the function f_sLPCRincreasing.

[0151] In the adaptation of the aforementioned threshold values, it has been shown that the increasing functions mentioned are for example functions for that which concerns the functions f_SLPCLand f_SLPCH. The function f_SPREDis a refined function of the variable stationarity parameter, of the form:

[0152] S_PRED(n)=α.STAT(n)+β

[0153] where α and β are two real values between 0 and 1 and where the value of S_PRED(n) is limited in the interval [S_PRED_m, S_PRED_M], S_PRED_mand S_PRED_Mrepresenting two experimentally determined values.

[0154] In order to limit yet again the risk of switching filters, it is then possible to choose, when the stationarity parameter STAT(n) is less than a given threshold value S_FWD, to require the “forward” LPC filtering mode.

[0155] On the other hand, the S_TRANS, S_STAT and G₁threshold values retain a fixed value, these values being able for example to be equal to −1 dB, 5 dB and 0 dB respectively.

[0156] The establishment of the decision function D_nand the obtaining of the analysis choice values d_n(n) are illustrated in the following way in FIG. 2f: following the aforementioned stage 120, carrying out a test stage 121 relative to the energy of the current LPC block B_n, by an inferiority comparison with the silence energy value ENER_SIL or with the value of the stationarity parameter STAT(n), compared by an inferiority comparison with the value S_FWDquoted previously in the description. On a positive response to the aforementioned test 121, the choice analysis value d_n(n) is taken as equal to 0, i.e. a symbolic value “fwd” in the stage 122.

[0157] On a negative response to the aforementioned test 121, a new test is carried out relative to the choice analysis value d_n−1(n−1) with the logical value 1, i.e. with the symbolic value “bwd”.

[0158] On a positive response to the aforementioned test 123, a new test is carried out on the aforementioned LPC filtering distance d_LPC, in a stage 124, in comparison with the threshold value S_LSP_H(n) by superiority comparison with this threshold value.

[0159] On a positive response to the aforementioned test 124, a new test 126ais carried out, consisting of comparing the “forward LPC filtering prediction gain, Gpf, with the “backward” LPC filtering prediction gain, Gpb, reduced by the threshold value S_TRANS.

[0160] On a positive response to the aforementioned test 126a, the logical value 0, symbolic value “fwd”, is attributed to the choice analysis value d_n(n), and on a negative response to the aforementioned test 126a, the same value of choice analysis is attributed the value 1, symbolic value “bwd”. The corresponding stages are called 128 and 129.

[0161] On a negative response to the previously mentioned test 124, a new test 125 is carried out. The test 125 consists in carrying out a comparison of the distance of the LPC filtering, d_LPC, by inferiority comparison with the threshold value S_LSP_L(n).

[0162] On a positive response to the test 125, a new test 126bis carried out by superiority comparison of the “backward” LPC filtering prediction gain with the “forward” LPC filtering prediction gain reduced by the previously mentioned value S_STAT.

[0163] On a positive response to the test 126b, the logical value 1 is attributed to the value of the choice analysis d_n(n) in the stage 129, i.e. the symbolic value “bwd”.

[0164] On a negative response to the test 126b, the logical value 0 is attributed to the value of the choice analysis d_n(n), i.e. the symbolic value “fwd”, stage 128.

[0165] In contrast, on a negative response to the test 125, a new test is carried out, in a stage 127, this test consisting of verifying the comparison conditions of the “backward” LPC filtering gain Gpb with the “forward” LPC filtering prediction gain reduced by the threshold value S_PRED(n), by superiority comparison of the intermediate LPC filtering prediction gain Gpi with the “forward” LPC filtering prediction gain value reduced by the aforementioned threshold value S_PRED(n) and by superiority comparison of the “backward” filtering prediction gain Gpb with the threshold value G₁, as well as comparison of the value of the intermediate filtering prediction gain Gpi with the threshold value G₁.

[0166] It is mentioned that the negative response to the test 123 previously mentioned in the description leads also to the carrying out of the aforementioned test 127.

[0167] On a positive response to the previously mentioned test 127, the logical value 1 is attributed to the value of the choice analysis d_n(n), i.e. the symbolic value “bwd” in the stage 129, whereas with a negative response to the aforementioned test 127, the logical value 0 is on the contrary attributed to the value of the choice analysis d_n(n), i.e. the symbolic value “fwd” in the stage 128.

[0168] Thus is set, by means of using the decision function D_n, the value of the choice analysis d_n(n) obtained with the aforementioned logical values 1 or 0, these logical values being however connected to a priority or absence of priority value of returning to the “backward” or “forward” filtering mode as a function of the value of the stationarity parameter.

[0169] A more detailed description of a coding device of an audio frequency digital signal by double analysis on the criterion of respectively “forward” or “backward” LPC choice analysis in a transmitted coded signal, according to the object of the present invention, will now be given in connection with FIG. 3.

[0170] In a practical manner, it is mentioned that the digital signal to be coded is subdivided into frames constituted by successive blocks of samples, each block comprising a given number N of samples for example.

[0171] In FIG. 3, constitution mode of the audio frequency digital signal to be coded in successive blocks of samples B_nhas not been shown for this operating mode is well known in the state of the technical art and can be carried out form a simple memory buffer, for example addressed to periodically read the frame frequency and the block frequency.

[0172] As shown furthermore in the aforementioned FIG. 3, the coding device which is the object of this invention includes a “forward” LPC analysis filter, carrying the reference 1A, and a “backward” analysis filter, carrying the reference 1B, in order to enable the delivery of a transmitted coded signal consisting of LPC filtering parameters accompanied by an analysis decision indication, as well as Pr parameters, relative to the harmonic analysis and to the excitation signal CELP.

[0173] Generally, it is shown that the analysis decision indication corresponds of the value of choice analysis d_n(n) as mentioned previously in the description. In so far as the LPC filtering parameters, it is mentioned that these correspond to specific parameters, according to the mode used of the coding method which is the object of the present invention as will be described later in the description.

[0174] In FIG. 3, also has been shown, in the coding device according to the invention, the existence of an adaptive filter operating as a function of the value of the stationarity parameter, this adaptive filter carrying the reference 1E. This adaptive filter 1E receives, it is understood of course, the original digital signal called s_n(t), i.e. the current block B_n. The filter 1E uses the filtering LPC parameter in order to calculate the residual signal which in turn is coded by the module 1F. These LPC parameters, as well as the filtering decision indication constitute a part of the coded signal which is transmitted to the decoder.

[0175] Furthermore, as shown in FIG. 3, the coding device which is the object of the present invention includes a coding means, carrying the reference 1F, of a non transmitted residue coding signal, the residue coding signal, designated by res_n(t)is directly available at the output of the adaptive filter 1E, this signal being thus delivered to the input with the audio frequency digital signal at the coding module of the not transmitted residue coding signal, in order to generate a synthesis residue signal, res_syn_n(t).

[0176] A reverse filtering module, carrying the reference 1G, receives the synthesis residue signal and enables the delivery of a synthesis signal referenced s_syn_n(t).

[0177] A memorization module 1H receives the aforementioned synthesis signal s_syn_n(t)in order to deliver the aforementioned synthesis signal for the previous block to the current block B_n, the synthesis signal thus obtained being designated by s_syn_n−1(t). This synthesis signal is delivered to the “backward” LPC analysis filter carrying the reference 1B in the aforementioned FIG. 3.

[0178] The coding device, which is the object of the present invention, as shown in FIG. 3, enables carrying out a coding of the audio frequency digital signal on the aforementioned audio frequency digital signal from the “forward” LPC filter for the non-stationary zones and on the aforementioned synthesis signal s_syn_n−1(t) from the “backward” LPC filter 1B for the stationary zones, as will be described below.

[0179] As will be observed in the aforementioned FIG. 3, the device which is the object of the invention comprises in this aim, for each current LPC block B_n, a calculation module 1C of the degree of stationarity of the audio frequency digital signal according to a stationarity parameter the value of which is between a maximum stationarity value and a minimum stationarity value. Of course, the stationarity parameter is the parameter STAT(n) previously described in the description according to the coding procedure which is the object of the present invention. The maximum and minimum stationarity values are also defined previously.

[0180] As has been shown furthermore in FIG. 3, the coding device which is the object of the invention includes a module, called 1D, for establishing from the aforementioned stationarity parameter STAT(n) a decision function and an LPC choice analysis value, the decision function being called D_nas previously mentioned in the description, and the LPC choice analysis value being of course corresponding to the value of the LPC choice analysis called d_n(n) previously mentioned in the description. It will be recalled that the value of the choice analysis d_n(n) can take the values 0 or 1, logical values, which correspond to the choice analysis symbolic value “fwd” and “bwd” for the “forward” and backward” LPC analysis respectively.

[0181] It is understood in particular that concerning the establishment of the decision function D_n, which corresponds for example to a software implementation, such as previously described in connection with FIG. 2f. Furthermore, the coding device according to the invention such as shown in FIG. 3 includes an LPC filtering analysis discrimination module, called 1D₂, this module receiving the value of the choice analysis d_n(n) and enabling delivering, for the current LPC block B_n, the value of the LPC “backward” and “forward” filtering parameters respectively as a function of the aforementioned value of choice analysis. It is clearly understood that the “backward” LPC filtering analysis as well as the “forward” LPC filtering analysis parameters are of course available in digital form at the filters carrying the reference 1B and 1A respectively in FIG. 3. These parameters are designated respectively Af_n(z) for the “forward” LPC filtering analysis parameters with regard to the “forward” LPC analysis filter, carrying the reference 1A, and by Ab_n(z) for the “backward” LPC analysis parameters with regard to the “backward” LPC analysis filter carrying the reference 1B. These parameters are delivered to the module 1D₁and the module 1D₂respectively.

[0182] As regards the creation of the equipment of the discrimination module 1D₁, it is shown that it may for example, in a non-limitative version, consist of two distinct memory zones enabling the memorization of the filtering parameters Af_n(z) and Ab_n(z) respectively, the value of choice analysis d_n(n) as a function of its current logical value, 0 or 1, enabling the addressing for reading the values of the memorized filtering parameters by the module 1D₂for example and the transmission of these filtering parameters by this latter.

[0183] Finally, as shown in FIG. 3, it has been shown that the coding device according to the object of the present invention, for the operation of the adaptive filter according to the stationarity value carrying the reference 1E , can be carried out by a filtering element the transfer function of which, called A(z), is established from the filtering parameter values delivered by the discrimination module 1D₂previously mentioned.

[0184] It is understood as well that the adaptive filtering module 1E can be achieved by a filter with adjustable coefficients, with the value of the coefficients of this latter delivered by the discrimination module 1D₂previously mentioned. The filtering carried out by the module 1E is thus of the adaptive type operating as a function of the degree of stationarity of the audio frequency digital signal to be coded. The module 1E thus delivers, from the original audio frequency digital signal sn_(t), the LPC filtering residue signal designated by res_n(t) to the coding module of the residue 1F, which enables then the delivery of the LPC synthesis residue signal designated by res_syn_n(t).

[0185] Finally, the module 1 G is a filtering module the transfer function of which is the reverse of the transfer function of the module 1E obtained form the memorized parameters of this latter. It receives the LPC synthesis residue signal res_syn_n(t) delivered by the coding module of the coding residue delivered by the module 1F. It is thus understood that the coding of the audio frequency digital signal s_n(t) is carried out in the module 1E through the LPC “forward” and “backward” analysis respectively which is carried out by the LPC “forward” and “backward” analysis filters 1A, 1B, the coded signal s_c_n(t) consisting in the transmission of the “forward” LPC filtering parameters when the value of the choice analysis d_n(n) has the symbolic value “fwd” as well as the indication of the choice analysis, i.e. of the value of the preceding quoted value of the choice analysis. This mode of operation enables carrying out the coding of the audio frequency digital signal and favoring holding it in one of the respectively “forward” and “backward” LPC filtering modes, as a function of the degree of stationarity of the digital signal, and limiting furthermore the number of switchings from one to the other of the considered filtering modes.

[0186] A decoding device of a coded audio frequency digital signal by double analysis on the criterion of respectively “forward” and “backward” LPC analysis, to a coded signal transmitted according to the coding method which is the object of the present invention, and by means of using a coding device such as shown in FIG. 3 for example, will now be described in connection with FIG. 4.

[0187] In a general manner, it is shown that the transmitted coded signal s_c_n(t) consists for each LPC analysis block of the value of the aforementioned choice analysis and, in the case where the value of choice analysis corresponds for the considered LPC analysis block to a “forward” LPC analysis, of the “forward” LPC filtering parameters as well as the coding parameters of the LPC filtering residue, Pr_nparameters, i.e. of the signal res_n(t) in a synthesis residue signal res_syn_n(t) by the residue coding module 1F.

[0188] As shown in FIG. 4, it is shown that the decoding device comprises at least a synthesis module, referenced 2A, of the filtering residue signal receiving the coding parameters of the LPC residue delivered by the module 1F. The module 2A decodes the coding parameters supplied by the module 1F and delivers consequently a synthesis residue signal, which is referenced in FIG. 4 res_syn_n(t).

[0189] The decoding device as shown in FIG. 4 comprises also a module, carrying the reference 2B, of reverse adaptive filtering as a function of the degree of stationarity, receiving the previously quoted synthesis residue signal, delivered by the module 2A, and enabling the generation of a synthesis signal s_syn_n(t) representative of the audio frequency digital signal, this signal constituting in fact the decoded signal.

[0190] It is of course understood that the reverse filtering module 2B uses the filtering parameters received by the decoder due to the fact of the transmission, are the “forward” LPC analysis parameters when these are transmitted and that the analysis decision corresponds to a “forward” LPC analysis or, in contrast, the “backward” filtering analysis parameters as will be described below.

[0191] With this aim, the decoding device which is the object of the present invention comprises of course a “backward” LPC filtering module, carrying the reference 2D, receiving the synthesis signal, i.e. the signal referenced s_syn_n(t) for the LPC block preceding the current LPC block, this synthesis signal being thus referenced s s_syn_n−1(t) in FIG. 4. With this aim, it is understood that the synthesis signal relative to the current block B_nand referenced s_syn_n(t) may then be delivered to the “backward” filtering module 2D by means of a memorization module, carrying the reference 2E, enabling in fact, by an adapted addressing for reading, to shift the reading of the synthesis signal to that corresponding to the block preceding, the current block B_n.

[0192] Finally, and to ensure the aforementioned operating mode, the decoding device which is the object of the present invention, as shown in FIG. 4, further includes a discriminator module carrying the reference 2C, enabling the carrying out of a “forward” and “backward” LPC discrimination analysis respectively. The module 2C receives, on the one hand, to control the discrimination, the value of choice analysis received, i.e. the value d_n(n), and, on the other hand, the “forward” LPC filtering parameters, i.e. the parameters Af_n(z) transmitted, as well as the “backward” LPC filtering parameters Ab_n(z) obtained by means of the module 2D. The module 2C thus enables delivering, as a function of the choice analysis value, i.e. of the value d_n(n), either the “forward” filtering parameters Af_n(z), or the “backward” filtering parameters Ab_n(z) to the reverse adaptive filtering module 2B as a function of the degree of stationarity.

[0193] As regards the material embodiment of the modules 2C and 2B, it is mentioned that these may simply consist of modules approximately identical to the modules 1D₂and 1E or, more particularly, 1G of FIG. 3.

[0194] As regards the effective embodiment of a coding device according to the object of the present invention, enabling using the procedure such as described previously in the description, two specific versions have been carried out.

[0195] Telephonic Band CELP Type Encoder, According to High Output Extension of the UIT-T Standard at 8 kb/s:

[0196] The actual encoder consisted of a telephonic band encoder from 300 to 3400 Hz, with an output of 12 kb/s of CELP type. The frames were constituted over a duration of 10 ms for an excitation supplied by algebraic codepages according to the technique called ACELP previously mentioned in the description.

[0197] The “forward” LPC analysis was an analysis of order 10 and the “backward” LPC analysis an analysis of order 30 every 80 samples.

[0198] A separation for the coding of the residue into two sub-blocks of 40 samples has been carried out. Each block B_nincluded 80 samples.

[0199] Adaptation of the Stationarity Parameter STAT(n)

[0200] The aforementioned stationarity parameter varies between two extreme values 0 and 100, the aforementioned values STAT_mand STAT_M.

[0201] The adaptation functions previously described in the description, and in particular the functions f_a(N_BWD) and f_b(N_BWD) were such that: $f_a\left(N_{—}\mathrm{BWD}\right)=\{\begin{array}{c}1.56\mathrm{if}{N}_{—}\mathrm{BWD}>20\\ 7.81\mathrm{if}{N}_{—}\mathrm{BWD}=20\\ 0\mathrm{otherwise}\end{array}$ $f_b\left(N_{—}\mathrm{BWD}\right)=\{\begin{array}{c}0.78·\left(20-N_{—}\mathrm{BWD}\right)\mathrm{if}{N}_{—}\mathrm{BWD}\le 20\\ 0\mathrm{otherwise}\end{array}$

[0202] In these relations, x=DGfb.

[0203] As regards the function f_renabling the refining value ΔS previously mentioned in the description to be established, this is a step function of the variable x, with x=Gpb−Gpf and ΔS=f_r(x) and having the value: $f_r\left(x\right)=\{\begin{array}{c}10\mathrm{if}x\ge 4\\ 7.5\mathrm{if}x\in [3;4[\\ 5\mathrm{if}x\in [2;3[\\ 2.5\mathrm{if}x\in [1;2[\\ 1.25\mathrm{if}x\in [0;1[\\ 0.625\mathrm{if}x\in [-1;0[\\ 0\mathrm{if}x\in [-2;-1[\\ -2.5\mathrm{if}x\in [-3;-2[\\ -5\mathrm{if}x\in [-4;-3[\\ -10\mathrm{if}x\in [-5;-4[\\ -20\mathrm{if}x<-5\end{array}$

[0204] The tuning of STAT(n) is furthermore subject to the following conditions previously mentioned in relation with FIG. 2c:

[0205] If ΔS>0:

[0206] If STAT*(n)<STAT_iSTAT(n)=STAT*(n)+ΔS

[0207] Otherwise STAT(n)=STAT*(n)

[0208] Otherwise:

[0209] STAT(n)=STAT*(n)

[0210] with STAT_i=40.6

[0211] The other test conditions referenced 1121d, 1121cand 1121fin FIG. 2chave not been used in the version.

[0212] Adaptation of the Decision Thresholds

[0213] As regards the decision thresholds:

[0214] S_PRED is adapted in the following manner:

[0215] S_PRED(n)=0.03.STAT(n)+1.0

[0216] S_PRED ε[S_PRED_m, S_PRED_M], S_PRED_m=1.03 and S_PRED_M=4;

[0217] The threshold S_LSP_L is adapted using the following step function: $S_{—}{\mathrm{LPS}}_{—}L\left(n\right)=f_{S_{-}{\mathrm{LPC}}_{-}L}\left(\mathrm{STAT}\left(n\right)\right)=\{\begin{array}{c}0.015\mathrm{if}\mathrm{STAT}\left(n\right)=100\\ 0\mathrm{otherwise}\end{array}$

[0218] The threshold value S_STAT used in case of stationarity of the LPC filters measured using the threshold S_LSP_L has been fixed at 4.0 dB.

[0219] The threshold S_LSP_H has not been used in this version.

[0220] The value of the threshold G₁has been fixed at 0 db.

[0221] As regards the energy value characterizing a silence frame ENER_SIL, this value has been fixed at 40 dB measured over the 80 samples s(i) of the current block B_n: ${\mathrm{ENER}}_{—}\mathrm{SIL}=10·\mathrm{Log}\left(\sum _{i=0}^{i(80}{\left(s\left(i\right)\right)}^2\right)$

[0222] As regards the value of the previously mentioned S_FWDthreshold intended to limit still further the risk of switching by imposing the “forward” LPC filter mode when the STAT(n) value is lower than this threshold, this S_FWDvalue has been set at 40.6.

[0223] A second Version of a CELP Broadened Band Encoder with Two Sub-bands 16/24/32 kb/s was Carried Out in the Following Conditions:

[0224] a broadened band encoder of 0 to 7000 Hz in two sub-bands. A main band was encoded with the CELP technique, frame with 120 samples, excitation created by algebraic codepages, and transmission of certain energy and spectrum characteristics of a host band of between 6000 Hz and 7000 Hz.

[0225] “forward” LPC analysis with 14 coefficients and “backward” LPC analysis with 50 coefficients every 120 samples. In “forward” LPC analysis mode, separation into two 60 sample LPC sub-blocks, the filter used for the first sub-block being interpolated from the current filter and the previous filter.

[0226] Calculation of the Stationarity Parameter STAT(n)

[0227] In this version, the aforementioned stationarity parameter varies between the two extreme values 0 and 120, the aforementioned STAT_mand STAT_Mvalues.

[0228] As regards the adaptation of the stationarity parameter value STAT(n), the values of the f_a(N_BWD) and f_b(N_BWD) functions are such that: $f_a\left(N_{—}\mathrm{BWD}\right)=\{\begin{array}{c}4\mathrm{if}{N}_{—}\mathrm{BWD}>10\\ 20\mathrm{if}{N}_{—}\mathrm{BWD}=10\\ 0\mathrm{otherwise}\end{array}{f}_b\left(N_{—}\mathrm{BWD}\right)=\{\begin{array}{c}10-N_{—}\mathrm{BWD}\mathrm{if}{N}_{—}\mathrm{BWD}\le 10\\ 0\mathrm{otherwise}.\end{array}$

[0229] As regards the fr function allowing the refining value ΔS previously mentioned in the description to be set, this is a step function of the variable x, with x=Gpb/Gpf and ΔS=f_r(x) and having a value of: $f_r\left(x\right)=\{\begin{array}{c}9\mathrm{if}x\ge 1.2\\ 6\mathrm{if}x\in [1.1;1.2[\\ 3\mathrm{if}x\in [1.05;1.1[\\ 1.5\mathrm{if}x\in [1.0;1.05[\\ 0.75\mathrm{if}x\in [0.95;1.0[\\ 0\mathrm{if}x\in [0.9;0.95[\\ -1.5\mathrm{if}x\in [0.85;0.9[\\ -3\mathrm{if}x\in [0.8;0.85[\\ -6\mathrm{if}x\in [0.75;0.8[\\ -12\mathrm{if}x<0.75\end{array}$

[0230] The tuning of STAT(n) is moreover subject to the following previously mentioned conditions in relation with FIG. 2c: If ΔS>0:

[0231] The other test conditions referenced 1121hand 1121din FIG. 2chave not been used in this version.

[0232] Adaptation of Decision Thresholds

[0233] As regards the decision thresholds:

[0234] S_PRED is adapted in the following way:

[0235] S_PRED(n)=0.03 STAT(n)−0.5 limited in the interval

[0236] [S_PRED_m, S_PREDM ]

[0237] with S_PRED_m=0.5 and S_PRED_M=2.5.

[0238] The S_LSP_L threshold is adapted with the help of the following step function: $S_{—}{\mathrm{LSP}}_{—}L\left(n\right)=f_{S_{-}{\mathrm{LSP}}_{-}L}\left(\mathrm{STAT}\left(n\right)\right)=\{\begin{array}{c}0.02\mathrm{if}\mathrm{STAT}\left(n\right)>100\\ 0.01\mathrm{otherwise}\end{array}$

[0239] The S_LSP_H threshold is adapted with the help of the following step function: $S_{—}{\mathrm{LSP}}_{—}H\left(n\right)=f_{S_{-}{\mathrm{LSP}}_{-}L}\left(\mathrm{STAT}\left(n\right)\right)=\{\begin{array}{c}0.08\mathrm{if}\mathrm{STAT}\left(n\right)>100\\ 0.01\mathrm{otherwise}\end{array}$

[0240] The value of the S_TRANS threshold used in the case of transition of the LPC filters measured with the help of the S_LSP_H threshold has been set at 0 dB.

[0241] The value of the S_STAT threshold used in the case of stationarity of the LPC filters measured with the help of the S_LSP_L threshold has been set at 2.5 dB.

[0242] The value of the G threshold has been set at 0 dB.

[0243] As regards the energy value characterizing a frame of silence ENER_SIL, this value has been set at 50 dB measured over the 120 samples s(i) of the current block B_n: ${\mathrm{ENER}}_{—}\mathrm{SIL}=\mathrm{Log}\left(\sum _{i=0}^{i<120}{\left(s\left(i\right)\right)}^2\right)$

[0244] As regards the value of the previously mentioned S_FWDthreshold intended to limit still further the risk of switching by imposing the “forward” LPC filter mode when the STAT(n) value is lower than this threshold, this S_FWDvalue has been set at 60.