WO2001024166A1

WO2001024166A1 - G.723.1 audio encoder

Info

Publication number: WO2001024166A1
Application number: PCT/SG1999/000096
Authority: WO
Inventors: Wenshun Tian
Original assignee: Stmicroelectronics Asia Pacific Pte Ltd
Priority date: 1999-09-30
Filing date: 1999-09-30
Publication date: 2001-04-05
Also published as: DE69926019D1; US6738733B1; EP1221162A1; EP1221162B1

Abstract

The invention relates to a speech encoder and corresponding encoding method for a dual rate coding system generally in accordance with the ITU-T G.723.1 recommendation, the dual rate coding system using Multi-Pulse Maximum Likelihood Quantization (MP-MLQ) and Algebraic Code-Excited Linear Prediction (ACELP) implementations for respective higher and lower transmission rates, wherein if the higher transmission rate is applicable, the number of gain scale factors searched in quantization steps is reduced, and if the lower transmission rate is applicable, the coding system has a threshold test for entry into the third processing loop of the ACELP implementation. The decoder may be a normal decoder. For example, a G.723.1 decoder may also be used to decode the bit streams encoded by a low complexity G.723.1 encoder in accordance with the invention.

Description

G.723.1 AUDIO ENCODER

The present invention relates to low complexity encoders, and more particularly, to low complexity encoders for implementing recommendation G.723.1 of the International Telecommunication Union (TTU-T).

Lower complexity compressors/decompressors (codecs) may be preferred for some computationally intensive applications. If the complexity of the codec is the bottleneck in a system, complexity reduction is desirable and can result in a significant reduction in millions of instructions per second (MIPS) required to be executed by the encoder.

The ITU-T recommendation G.723.1, incorporated herein by reference, relates to dual rate speech coding for multimedia communications transmitting at 5.3 and 6.3 Kbps. The recommendation prescribes certain methods of implementation for each of these transmission rates. The 6.3 Kbps codec has better quality and uses Multi-Phase Maximum Likelihood Quantization (MP-MLQ) for fixed codebook excitation. The 5.3 Kbps codec uses Algebraic Code-Excited Linear Prediction (ACELP). A functional module of the codec which executes these two encoding methods bears almost half of the computational load of the entire G.723.1 speech coder. If the methods executed by the functional module are made to have a decreased computational load, the G.723.1 speech coder will have an increased efficiency.

The present invention provides a method of reducing the computational load of a dual rate encoding system, the encoding system being configured to transmit at a first transmission rate using a Multi-Pulse Maximum Likelihood Quantization (MP-MLQ) process or at a second transmission rate using an Algebraic Code-Excited Linear Prediction (ACELP) process, wherein the MP-MLQ process normally searches subframes of excitation signals according to a nominal number of gain scale factors in the execution of quantization steps for encoding the speech signals and the ACELP process normally imposes a first correlation threshold test for entry into an embedded signal processing loop, the method including the step of: for the MP-MLQ process, reducing the number of gain scale factors employed in the quantization steps, thereby reducing the number of searches, which in turn reduces the computational load; or for the ACELP process, imposing a second correlation threshold test for entry into a previous signal processing loop in which the embedded signal processing loop is embedded, thereby reducing the number of times the previous signal processing loop and the embedded signal processing loop are entered, which in turn reduces the computational load.

The present invention further provides a dual rate speech coding system having a reduced computational load, the encoding system having Multi-Pulse Maximum Likelihood Quantization (MP-MLQ) processing means for transmitting at a first transmission rate and Algebraic Code-Excited Linear Prediction (ACELP) processing means for transmitting at a second transmission rate, wherein the MP-MLQ processing means normally searches subfram.es of excitation signals according to a nominal number of gain scale factors in quantization of the speech signals, and the ACELP processing means normally uses a first correlation threshold test for allowing entry into an embedded signal processing loop, wherein: the MP-MLQ processing means has a reduced number of gain scale factors for reducing the number of searches and thereby reducing the computational load; the ACELP processing means uses a second correlation threshold test for allowing entry into a previous signal processing loop in which the embedded signal processing loop is embedded, tliereby reducing the number of times the previous signal processing loop and the embedded signal processing loop are entered, which in turn reduces the computational load.

Advantageously, embodiments of the invention simplify the ACELP and MP-MLQ methods by reducing the number of recursions which make less contribution to the metrics. This is achieved by selecting less gain levels or putting an extra threshold to decrease the chance to enter the most computational intensive loops.

Advantageously, the proposed encoder scheme is applicable for both ITU-T recommendations G.723.1 and G.723.1A. For ACELP excitation, further complexity reduction is possible by adjusting the thresholds. This complexity reduction for ACELP excitation is also applicable for G.729 and its annexes. The present invention will now be described in further detail, by way of example only, with reference to the accompanying drawing.

Figure 1 is a block diagram of the G.723.1 speech coder.

Reference is also made to the following procedures which are appended to this description.

Procedure 1 is a pseudocode representation of the standard MP-MLQ procedure of the G.723.1 speech coder;

Procedure 2 is a pseudocode representation of the MP-MLQ procedure of an embodiment of the present invention;

Procedure 3 is a pseudocode representation of the standard ACELP procedure of the G.723.1 speech coder; Procedure 4 is a pseudocode representation of the ACELP procedure of an embodiment of the present invention.

A MP-MLQ/ACELP block 10 for implementing the MP-MLQ and ACELP excitation methods is shown in Figure 1. These methods take up almost half of computational load of the whole codec. Since embodiments of the present invention only relate to these two fixed codebook excitation methods, the description relates only to these excitation techniques and not to other parts of the G.723.1 speech coder. Apart from the fixed codebook excitation part (i.e. block 10), all other modules are the same for the dual rate coders. The decoding scheme, for decoding bit streams encoded with the low complexity encoder, remains the same as for the normal ITU-T G.723.1 recommendation.

MP-MLQ Excitation (normal complexity)

The object of the quantization procedure is to find the optimized excitation e_u(n) which makes the mean square error minimum, based on an analysis by synthesis method. The excitation signal is given by

fc-0

(1) where G_u is the gain factor, δ(n) is a Dirac function, (±1)

and positions of the Dirac functions respectively, and N„ is the number of pulses, which is 5 for odd subframes and 6 for even subframes. The pulse positions in are either all odd or all even. This is indicated by a grid bit.

The scalar gain quantizer consists of 24 steps, of 3.2 dB each, Around the quantized value, G_u, additional gain values are selected within the range [G_u - 6.4dB; G_u + 3.2dB]. The optimal combination of pulse locations and gains are then transmitted to the remaining encoder modules.

To improve the quality of speech with a short pitch period, the following additional procedure is used. If the pitch lag is less than 58 samples for a particular subframe, a train of Dirac functions with a period of the pitch index is used for each location ξ_k instead of a single Dirac function in the above quantization procedure. The choice between a train of Dirac functions or a single Dirac function to represent the residual signal is made based on the mean square error computation. The configuration which yields the lowest mean square error is selected.

Based on the above brief description of MP-MLQ, the optimization procedure is represented in pseudocode as shown in Procedure 1. The symbols InsCI inside the brackets are the cycles needed for a given processor; and the number of cycles if using, for example, a D950 processor. The D950 is a normal 16-bit fixed-point digital signal processor (DSP) made by STMicroelectronics. Other 16-bit fixed-point DSPs are the ADSP-2181 by Analog Devices and the TMS320C54x series by Texas Instruments. Although the number of instructions required to execute the same function may vary among different DSPs, the invention will still achieve a significant savings in MIPS for each appropriate DSP.

The worst case for MP-MLQ is that above optimization procedure is conducted twice when the pitch is less than 58 samples. The total number of cycles per subframe is given by Total cycles per subframe =2 x{rπsCl+2x[lnsC2+InsC6+4x(InsC3+IτιsC4+InsC5)]}

(2)

Therefore the total number of cycles per subframe for the procedure of Procedure 1 is 64368 if using the D950 processor.

ACELP excitation (normal complexity)

For the ACELP technique for fixed codebook excitation, a 17-bit algebraic codebook is used for the stochastic codebook excitation e_ω(n). Each fixed codevector contains four non-zero pulses which can assume the signs and positions given in the following table.

Table 1 ACELP excitation codebook

Sign Positions

±1 0 8 16 24 32 40 48 56

±1 2 10 18 26 34 42 50 58

±1 4 12 20 28 36 44 52 (60)

±1 6 14 22 30 38 46 54 (62)

In the table, all pulses are in the even positions but the positions of all pulses can be simultaneously shifted by one (to occupy odd positions) when requiring one extra bit. Note that the last position of each of the last two pulses falls outside the subframe boundary, which signifies that the pulse is not present. Each pulse position is encoded with 3 bits and each pulse sign is encoded with 1 bit. This gives a total of 16 bits for the 4 pulses. Further, an extra bit is used to encode the shift. The excitation sequence is defined as

^e«.m(ⁿ) = θoδ(n-ξo) ^{+ α}ι^δ(ⁿ-ξt) + θiδ(n-ξ₂) + a₃b(n-ξj

(3)

where ξ_k is the position of the k¹'¹ pulse and a_k is its sign (± 1). A focused search approach is used to simplify the search procedure. To limit the number of times entering the last loop, a threshold is applied and the last loop is entered only if this threshold is exceeded. The maximum number of times the loop can be entered is fixed so that a low percentage of the codebook is searched. The maximum absolute correlation C,^ and the average correlation C^_j due to the contribution of the first three pulses are found prior to the codebook search. The threshold is given by:

tΛr3 = C_ϊv3+(C_mιtt3-C„_V3)/2

(4) The fourth loop is entered only if the absolute correlation (of the three pulses) exceeds thrβ.

To further control the search, the number of times the last loop is entered (for the 4 subframes) is not allowed to exceed 600. (The average worst case per subframe is 150 times).

Based on the above brief description of ACELP, the optimization procedure is represented in pseudocode as shown in Procedure 3. In Procedure 3, InsCi is the number of instruction cycles, foDowed by an example number of cycles for the D950 implementation. The total cycles are calculated by

InsCl +InsCl2+Sx(InsC2+InsCU) + 8^lx(InsC3 +InsC10) + 8²x(InsC4+InsC9)+ time₃x(JnsC5 +InsCl + &xInsC6) + (8³-ήme )xInsCS

=27207 + time _jx238 [if using a D950 processor] =62907 [if time_z = 150]

(5)

where ιime₃ is the number of times entering the last loop. At the worst case, the maximum number of time^ is set to 150. Therefore the worst case cycles per 7.5 ms subframe are 62907 if using a D950 processor, which equates to 8.4 MIPS.

Lower Complexity Implementations

In embodiments of the invention the modules (codes) may be shared by both G.723.1 and the lower complexity implementation of the G.723.1 coder (LC-G.723.1). Preferably, the coding system is selectable between bit-exact G.723.1 and LC-G.723.1 coders, leading to an embedded system. This is shown by the procedure as follows:

If (lower complexity enabled)

If (6.3kbit/s)

6.3kbit s LC-G723.1 encoding Else 5.3kbit/s LC-G723.1 encoding Else If (6.3kbit/s)

6.3 Kbit/s encoding Else 5.3 Kbit/s encoding

For the low-complexity encoding of 6.3 Kbps and 5.3 Kbps codecs in accordance with the present invention, the operation procedures are shown in Procedure 2 and Procedure 4 respectively.

One of the characteristics of MP-MLQ is that the latter pulse contribution will be added upon the previous one and all pulses are scaled by one gain. For each new found pulse, the gain is further fine tuned within the range [-6.4dB;-3.2dB; 0; +3.2dB]. Since all pulses share one gain, the observation is that the gain level decreases as the number of found pulses increases. Due to the characteristic of MP-MLQ, the additional higher gain levels (0 and +3.2dB) are rarely selected. In this simplification, we only use two gain levels, i.e. -6.4 dB and -3.2 dB around the previous quantized gain. Therefore the number of instructions inside the gain searching loop can be decreased by about half for each subframe when the pitch lag is less than 58 samples.

The worst case number of cycles for MP-MLQ is calculated as:

Total cycles per subfraine +2x{.InsC2+lnsC6+2x(InsC3+ InsCi +InsC5)]}

(6) For the D950 example, the total number of cycles per subframe is 39424.

For an adaptive codebook search, the worst case is when the pitch lag ≥ 58, which is just the opposite of fixed codebook excitation. If the number of gain levels decreases from 4 to 2 for fixed codebook excitation, the computational load is reduced from Equation (2) to Equation (6). To balance the computational load for all cases, the codes are also simplified for when the pitch lag ≥ 58. The number of searched gain levels is reduced from 4 to 3, i.e. -6.4, -3.2 and 0 dB. (please refer to Procedure 2).

The number of cycles per subframe for MP-MLQ with a pitch lag ≥ 58 is calculated as

Total cycles per subframe = Iπ$Cl +2x[InsC2+InsC6+3x(lnsC3+InsC4+InsC5)]

(7)

The total number of cycles per subframe would then be 19826 for the D950 processor example.

Comparing Equations (2) and (6), 24944 cycles per subframe can be saved at worst case (of MP-MLQ) if using the D950 processor. This equates to a saving of 3.3 MIPS. For the normal case, in which MP-MLQ is conducted once, the saved cycles are 12358 per subframe, which equates to 1.65 MIPS. This unbalanced complexity reduction in the fixed codebook search (MP-MLQ) corresponds to the unbalanced computational load adaptive codebook search, in which, for example, about 30,000 and 46,000 cycles are needed respectively for the worst case and normal case of MP-MLQ.

A purpose of embodiments of the invention is to reduce the complexity for the worst case scenario (i.e., under the most intensive computational load). If the complexity is reduced in the worst case, the overall MIPS requirement is reduced accordingly. At the higher bit rate, the most complex modules are the fixed codebook excitation module (MP-MLQ) and adaptive ^' excitation module. The complexity of these two modules changes depending on the pitch lag, while other modules are relatively stable in terms of computational load. Shown in Table 2 below is a comparison of the MIPS requirements for the worst case (pitch lag < 58 samples) and the normal case (pitch lag ≥58) for a D950 DSP.

Table 2 Complexity comparison 6.3kbits/s for one subframe (7.5ms)

From Procedure 3 and Equation (5), it is apparent that any instructions inside the i₂ and i₃ loops will be executed hundreds of times. It may be advantageous to further limit the numbers entering these two loops. Instead of using one threshold, two thresholds are used.

Both the maximum absolute correlation and the average correlation due to the contribution of the first two and three pulses, C_max2 and C_m2, and Q_imx3 and _vJ , are found prior to the codebook search. The thresholds are calculated by:

thr2 =(C_flv2 + (C_miu,₂-C_ev2) / 4)

(8) thr 3= (C_flVi + (C_mttrJ - C_avJ) / 3)

(9)

Now we have two thresholds. To further control the search, the average number of times the third and last loops are entered is not allowed to exceed 32 and 75 (for example), respectively for each subframe. The proposed low-complexity ACELP optimization procedure is modified as in Procedure 4. The total number of cycles per subframe is given by:

Inscl + InsC\2 + 8 x (InsC2 4- InsCU) + 64 x (7^C3 + ΛwClO) + tlme₂ x (InsQU + InsCU) +(64 - tim ^ x InsC 15 -I- time₂ x 8 x (InsCA + InsC9) + time₃ x 8 x (InsCl + 8 x //w C6) + (tim ₂ x 8 -

x 7ΛSC8 = 8373 + time₂ x 336 + time₃ x 238

(10)

where ιime₂ and t^e_, are the number of times the processor enters into the 3"¹ and 4^th loops respectively. For the worst case, the time₂ and tim are set to 32 and 75 respectively. Therefore the worst case number of cycles will become 36976. Comparing with Equation (5), 25932 cycles or 3.45 MIPS can be saved (if using the D950 processor).

It should be noted that further complexity reduction is simple to effect for this ACELP excitation by choosing smaller time₂ and tirrψ parameters and corresponding higher thresholds. The proposed parameters for this LC-G.723.1 are based on the objective that LC- G.723.1 should have similar performances to G.723.1. If further reduction of complexity is needed, the performance will be smoothly degraded. For example, by increasing the threshold levels and corresponding allowed loop entry times time₂ and time₃ to 20 and 60 respectively, a further 1.01 MIPS can be saved.

PROCEDURE 1

do twice or once depending on the pitch lag pre-scarch (InsCl, 6116) for k = 0,1 (grid) fmd the first pulse and quantize gain (InsCl, 537) for i = 1, 2, 3, 4 (gain) prepare searching (InsC3, 173) find remain pulses (InsCA, 2261) evaluate the error (7ΛSC5, 687) post process (InsC6, 13)

PROCEDURE 2

if pitch lag < 58 do twice

pre-search (7π.sCl, 6116) for k a 0,1 (grid) find the first pulse and quantize gain (InsCl, 2537) for i = 1, 2 (gain) prepare searching (InsC3, 173) find remain pulses (7ΛSC4, 2261) evaluate the error (InsC5, 687) post process (InsC6, 13)

else do once

pre-search (InsC 1 , 116) for k = 0,1 (grid) find the first pulse and quantize gain (InsCl, 2537) for i = 1, 2, 3 (gain) prepare searching (7πi'C3, 173) fmd remain pulses (7/wC4, 2261) evaluate the error (InsC5, 687) post process (InsCβ, 13) PROCEDURE 3

pre-search (InsCl ,4195) far iO = 0, ... ,7 (first pulse) partial correlation & energy calculation (InsCl, 12) for il = 0,..., 7 (second) partial correlation Sc energy calculation (InsC3, 26) for 12 = 0 7 (third) partial correlation & energy calculation (7tz.yC4, 30) if (correlation > threshold) prepare for last loop (InsCS, 14) for i3 = 0, ..., 7 Qast)

{correlation & energy calculation and evaluation (InsCβ, 25)} update (InsCl, 26) else

Update (InsC8, 2) post process (7rz5C9, 5) post process (InsCIO, 22) do post (7?uCll, 17) post process (7rcsC12, 746) PROCEDURE 4

pre-search (InsCl, 4915) fa r iO = 0.....7 (first loop) partial correlation & energy calculation (InsCl, 12) for i! = 0,..., 7 (second) partial correlation & energy calculation (InsC3, 16) if (correlation > threshold!) preparation (7τwC 13, 14) for i2 = 0,..., 7 (third) partial correlation & energy calculation (7>wC4, 30) if (correlation > threshold) preparation (InsC5, 14) for i3 = 0,..., 7 (last) {correlation & energy calculation, and evaluation

(InsC6,25)} post process else update (InsC8, 2) do post process (InsC9, 5) post process (7Λ_?C14, 26) else update (InsCl5, 2) post process (7ΛSC10, 22) post process (InsCU, 17) post process (7?LSC12, 746)

Claims

CLAIMS:

1. A method of reducing the computational load of a dual rate encoding system, the encoding system being configured to transmit at a first transmission rate using a Multi-Pulse Maximum Likelihood Quantization (MP-MLQ) process or at a second transmission rate using an Algebraic Code-Excited Linear Prediction (ACELP) process, wherein the MP-MLQ process normally searches subframes of excitation signals according to a nominal number of gain scale factors in the execution of quantization steps for encoding the speech signals and the ACELP process normally imposes a first correlation threshold test for entry into an embedded signal processing loop, the method including the step of: for the MP-MLQ process, reducing the number of gain scale factors employed in the quantization steps, thereby reducing the number of searches, which in turn reduces the computational load; or for the ACELP process, imposing a second correlation threshold test for entry into a previous signal processing loop in which the embedded signal processing loop is embedded, thereby reducing the number of times the previous signal processing loop and the embedded signal processing loop are entered, which in turn reduces the computational load.

2. The method of claim 1, wherein the second threshold test is apphcable for entry into the third of four signal processing loops.

3. The method of claim 2, wherein if the second transmission rate is applicable, the method further includes the step of substituting for the first threshold a higher threshold for entry into the fourth signal processing loop.

4. The method of claim 3, wherein if the second transmission rate is applicable, further including the step of limiting the number of times the third or fourth signal processing loop may be entered.

5. The method of claim 4, wherein the third or fourth signal processing loops may be entered up to 32 or 75 times respectively for each of the speech subframes.

6. The method of claim 5, wherein the dual rate coding system is generally in accordance with the ITU-T G.723.1 recommendation.

7. The method of claim 1, wherein if a pitch lag of the subframe is less than a predetermined parameter, the number of gain scale factors searched is reduced from four to two.

8. The method of claim 7, wherein if the pitch lag of the subframe is equal to or greater than the predetermined parameter, the number of gain scale factors searched is reduced from four to three.

9. The method of claim 8, wherein the predetermined parameter is 58.

10. The method of claim 1 or 9, wherein the quantization steps and a pre-search are executed once if the pitch lag is greater than or equal to 58 and twice if the pitch lag is less than 58.

11. A dual rate speech coding system having a reduced computational load, the encoding system having Multi-Pulse Maximum Likelihood Quantization (MP-MLQ) processing means for transmitting at a first transmission rate and Algebraic Codc-Excited Linear Prediction (ACELP) processing means for transmitting at a second transmission rate, wherein the MP- MLQ processing means normally searches subframes of excitation signals according to a nominal number of gain scale factors in quantization of the speech signals, and the ACELP processing means normally uses a first correlation threshold test for allowing entry into an embedded signal processing loop, wherein: the MP-MLQ processing means has a reduced number of gain scale factors for reducing the number of searches and thereby reducing the computational load; the ACELP processing means uses a second correlation threshold test for allowing entry into a previous signal processing loop in which the embedded signal processing loop is embedded, thereby reducing the number of times the previous signal processing loop and the embedded signal processing loop are entered, which in turn reduces the computational load.

12. The coding system of claim 11, wherein the second correlation threshold test is used to allow entry into the third of four signal processing loops.

13. The coding system of claim 12, wherein if the second transmission rate is applicable, 5 the coding system has a higher threshold, used in place of the first threshold, for entry into the fourth signal processing loop.

14. The coding system of claim 13, wherein if the second transmission rate is applicable, the number of times the third or fourth signal processing loop may be entered is limited.

10

15. The coding system of claim 14, wherein the limit for entry into the third or fourth signal processing loops is 32 or 75 times respectively for each of a plurality of speech subframes.

15 16. The coding system of claim 11 or 15 , wherein the dual rate coding system is generally in accordance with the ITU-T G.723.1 recommendation.

17. The coding system of claim 11, wherein if a pitch lag of the subframe is less than a predetermined parameter, the number of gain scale factors searched is reduced from four to 0 two.

18 . The coding system of claim 17, wherein if the pitch lag of the subframe is equal to or greater than the predeteπnined parameter, the number of gain scale factors searched is reduced from four to three. 5

19. The coding system of claim 18, wherein the predetermined parameter is 58.

20. The coding system of claim 11 or 19, wherein the quantization steps and a pre-search are executed once if the pitch lag is greater than or equal to 58 and twice if the pitch lag is 0 less than 58.