CA2618002C - Algebraic codebook with signal-selected pulse amplitudes for fast coding of speech - Google Patents

Algebraic codebook with signal-selected pulse amplitudes for fast coding of speech Download PDF

Info

Publication number
CA2618002C
CA2618002C CA002618002A CA2618002A CA2618002C CA 2618002 C CA2618002 C CA 2618002C CA 002618002 A CA002618002 A CA 002618002A CA 2618002 A CA2618002 A CA 2618002A CA 2618002 C CA2618002 C CA 2618002C
Authority
CA
Canada
Prior art keywords
amplitude
speech signal
pulse
codebook
zero
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Lifetime
Application number
CA002618002A
Other languages
French (fr)
Other versions
CA2618002A1 (en
Inventor
Jean-Pierre Adoul
Claude Laflamme
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Universite de Sherbrooke
Original Assignee
Universite de Sherbrooke
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Priority claimed from US08/508,801 external-priority patent/US5754976A/en
Application filed by Universite de Sherbrooke filed Critical Universite de Sherbrooke
Publication of CA2618002A1 publication Critical patent/CA2618002A1/en
Application granted granted Critical
Publication of CA2618002C publication Critical patent/CA2618002C/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10LSPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L19/00Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis
    • G10L19/04Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis using predictive techniques
    • G10L19/08Determination or coding of the excitation function; Determination or coding of the long-term prediction parameters
    • G10L19/10Determination or coding of the excitation function; Determination or coding of the long-term prediction parameters the excitation function being a multipulse excitation
    • G10L19/107Sparse pulse excitation, e.g. by using algebraic codebook

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computational Linguistics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Algebra (AREA)
  • General Physics & Mathematics (AREA)
  • Signal Processing (AREA)
  • Health & Medical Sciences (AREA)
  • Audiology, Speech & Language Pathology (AREA)
  • Human Computer Interaction (AREA)
  • Acoustics & Sound (AREA)
  • Multimedia (AREA)
  • Compression, Expansion, Code Conversion, And Decoders (AREA)

Abstract

The present invention relates to a method and device for conducting a search in a codebook. This codebook consists of a set of pulse amplitude/position combinations each defining a number L of positions p and comprising both zero-amplitude pulses and non-zero-amplitude pulses assigned to respective positions p = 1, 2, ...L of the combination. Also, each non-zero-amplitude pulse assumes one of q possible amplitudes. According to the method, a subset of combinations is pre-selected from the codebook, and the search is limited to this subset to reduce complexity thereof. To pre-select the subset, an amplitude/position function is pre-established in relation to the sound signal. Pre-establishin g the amplitude/position function includes pre-assigning one of the q possible amplitudes to each position p by (i) processing the sound signal to produce a backward-filtered target signal D and a pitch-removed residual signal R', (ii) calculating an amplitude estimate vector B in response to the signals D and R', and (iii) for each position p, quantizing an amplitude estimate B p of the vector B to obtain the amplitude to be selected for that particular position p.

Description

ALGEBRAIC CODEBOOK WITH SIGNAL-SELECTED
PULSE AMPLITUDES FOR FAST CODING OF SPEECH

BACKGROUND OF THE INVENTION
1. Field of the invention:

The present invention relates to an improved technique for digitally encoding a sound signal, in particular but not exclusively a speech signal, in view of transmitting and synthesizing this sound signal.
2. Brief description of the prior art:

The demand for efficient digital speech encoding techniques with a good subjective quality/bit rate trade-off is increasing for numerous applications such as voice transmission over satellites, land mobile, digital radio or packed network, voice storage, voice response and wireless telephony.

One of the best prior art techniques capable of achieving a good quality/bit rate trade-off is the so-called Code Excited Linear Prediction (CELP) technique. According to this technique, the speech signal is sampled and processed in blocks of L samples (i.e. vectors), where L is some predetermined number. The CELP technique makes use of a codebook.

A codebook, in the CELP context, is an indexed set of L-sample-long sequences which will be referred to as L-dimensional codevectors (pulse combinations defining L different positions and comprising both zero-amplitude pulses and non-zero-amplitude pulses assigned to respective positions p=1, 2, ...L of the combination). The codebook comprises an index k ranging from 1 to M, where M represents the size of the codebook sometimes expressed as a number of bits b:
M=2b A codebook can be stored in a physical memory (e.g. a look-up table), or can refer to a mechanism for relating the index to a corresponding codevector (e.g. a formula).

To synthesize speech according to the CELP technique, each block of speech samples is synthesized by filtering the appropriate codevector from the codebook through time varying filters modelling the spectral characteristics of the speech signal. At the encoder end, the synthetic output is computed for all or a subset of the candidate codevectors from the codebook (codebook search). The retained codevector is the one producing the synthetic output that is the closest to the original speech signal according to a perceptually weighted distortion measure.
A first type of codebooks is the so-called "stochastic"
codebooks. A drawback of these codebooks is that they often involve substantial physical storage. They are stochastic, i.e. random in the sense that the path from the index to the associated codevector involves look-up tables that are the result of randomly generated numbers or statistical techniques applied to large speech training sets. The size of stochastic codebooks tends to be limited by storage and/or search complexity.

A second type of codebooks are the algebraic codebooks.
By contrast with the stochastic codebooks, algebraic codebooks are not random and require no storage. An algebraic codebook is a set of indexed codevectors in which the amplitudes and positions of the pulses of the e codevector can be derived from its index k through a rule requiring no, or minimal, physical storage. Therefore, the size of an algebraic codebook is not limited by storage requirements. Algebraic codebooks can also be designed for efficient search.

OBJECTS OF THE INVENTION

An object of the present invention is therefore to provide a method and device for drastically reducing the complexity of the codebook search upon encoding a sound signal, these method and device being applicable to a large class of codebooks.
SUMMARY OF THE INVENTION

More particularly, in accordance with the present invention, there is provided a method of conducting a search in a codebook during encoding of a sound signal. The codebook contains a set of pulse amplitude/position combinations defining a number L of different positions p and comprising both zero-amplitude pulses and non-zero-amplitude pulses assigned to respective positions p = 1, 2, ...L of the pulse amplitude/position combination. Each non-zero-amplitude pulse assumes one of q possible amplitudes. This codebook search conducting method comprises pre-selecting from the codebook a subset of pulse amplitude/position combinations in relation to the sound signal, and searching only this subset of pulse amplitude/position combinations during encoding of the sound signal whereby complexity of the search is reduced as only a subset of the pulse amplitude/position combinations of the codebook is searched. Pre-selecting a subset of pulse amplitude/position combinations comprises pre-establishing, in relation to the sound signal, an amplitude/position function between the positions p= 1, 2, ... L and the q possible amplitudes. Pre-establishing an amplitude/position function comprises pre-assigning one of the q possible amplitudes as valid amplitude to each position p. Pre-assigning one of the q possible amplitudes to each position p comprises (a) processing the sound signal to produce a backward-filtered target signal D and a pitch-removed residual signal R', (b) calculating an amplitude estimate vector B in response to the backward-filtered target signal D and to the pitch-removed residual signal R', and (c) for each position p, quantizing an amplitude estimate BP of the vector B to obtain the amplitude to be selected for that position p. Finally, searching the subset of pulse amplitude/position combinations comprises limiting the search to the pulse amplitude/position combinations of the codebook having non-zero-amplitude pulses which satisfy the pre-established function.

The present invention also relates to a device for 5 conducting a search in a codebook during encoding of a sound signal. The codebook contains a set of pulse amplitude/position combinations each defining a number L of different positions p and comprising both zero-amplitude pulses and non-zero-amplitude pulses assigned to respective positions p = 1, 2, ...L of the pulse amplitude/position combination. Each non-zero-amplitude pulse assumes one of q possible amplitudes. This codebook search conducting device comprises means for pre-selecting from the codebook a subset of pulse amplitude/position combinations in relation to the sound signal, and means for searching only the subset of pulse amplitude/position combinations during encoding of the sound signal whereby complexity of the search is reduced as only a subset of the pulse amplitude/position combinations of the codebook is searched. The pre-selecting means comprises means for pre-establishing, in relation to the sound signal, an amplitude/position function between the positions p=1, 2, ...L and the q possible amplitudes, and the pre-establishing means comprises means for pre-assigning one of the q possible amplitudes as valid amplitude to each position p_ The means for pre-assigning one of the q possible amplitudes to each position p comprises (a) means for processing the sound signal to produce a backward-filtered target signal D
and a pitch-removed residual signal R', ( b) means for calculating an amplitude estimate vector B in response to the backward-filtered target signal D and to the pitch-removed residual signal R', and ( c) means for quantizing, for each of the positions p, an amplitude estimate BA of the vector B to obtain the amplitude to be selected for the position p. Finally, the searching means comprises means for limiting the search to the pulse amplitude/position combinations of the codebook having non-zero-amplitude pulses which satisfy the pre-established function.

The present invention further relates to a cellular communication system for servicing a large geographical area divided into a plurality of cells, comprising:
mobile transmitter/receiver units;
cellular base stations respectively situated in the cells;
means for controlling communication between the cellular base stations;
a bidirectional wireless communication sub-system between each mobile unit situated in one cell and the cellular base station of that cell, this bidirectional wireless communication sub-system comprising in both the mobile unit and the cellular base station (a) a transmitter including means for encoding a speech signal and means for transmitting the encoded speech signal, and (b) a receiver including means for receiving a transmitted encoded speech signal and means for decoding the received encoded speech signal.
The speech signal encoding means comprises means responsive to the speech signal for producing speech signal encoding parameters, and these speech signal encoding parameter producing means comprises a device as described hereinabove, for conducting a search in a codebook in view of producing at least one of the speech signal encoding parameters, wherein the speech signal constitutes the sound signal.
According to the invention, there are further provided:

- A cellular network element comprising (a) a transmitter including means for encoding a speech signal and means for transmitting the encoded speech signal, and (b) a receiver including means for receiving a transmitted encoded speech signal and means for decoding the received encoded speech signal. The speech signal encoding means comprises means responsive to the speech signal for producing speech signal encoding parameters, and these speech signal encoding parameter producing means comprises a device as described hereinabove, for conducting a search in a codebook in view of producing at least one of the speech signal encoding parameters, wherein the speech signal constitutes the sound signal.

- A cellular mobile transmitter/receiver unit comprising (a) a transmitter including means for encoding a speech signal and means for transmitting the encoded speech signal, and (b) a receiver including means for receiving a transmitted encoded speech signal and means for decoding the received encoded speech signal. The speech signal encoding means comprises means responsive to the speech signal for producing speech signal encoding parameters, and these speech signal encoding parameter producing means comprises the above described device for conducting a search in a codebook in view of producing at least one of the speech signal encoding parameters, wherein the speech signal constitutes the sound signal.

- In a cellular communication system for servicing a large geographical area divided into a plurality of cells, and comprising: mobile transmitter/receiver units; cellular base stations respectively situated in these cells; and means for controlling communication between the cellular base stations;
a bidirectional wireless communication sub-system between each mobile unit situated in one cell and the cellular base station of that cell, this bidirectional wireless communication sub-system comprising in both the mobile unit and the cellular base station (a) a transmitter including means for encoding a speech signal and means for transmitting the encoded speech signal, and (b) a receiver including means for receiving a transmitted encoded speech signal and means for decoding the received encoded speech signal. The speech signal encoding means comprises means responsive to the speech signal for producing speech signal encoding parameters, and these speech signal encoding parameter producing means comprises the above described device for conducting a search in a codebook in view of producing at least one of the speech signal encoding parameters, wherein the speech signal constitutes the sound signal.

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

BRIEF DESCRIPTION OF THE DRAWINGS

In the appended drawings:
Figure 1 is a schematic block diagram of a sound signal-encoding device comprising an amplitude selector and an optimizing controller in accordance with the present invention;
Figure 2 is a schematic block diagram of a decoding device associated with the encoding device of Figure 1;

Figure 3a is a sequence of basic operations for the fast codebook search in accordance with the present invention, based on signal-selected pulse amplitudes;

Figure 3b is a sequence of operations for pre-assigning one of the q amplitudes to each position p of the pulse amplitude/position combinations;

Figure 3c is a sequence of operations involved in the N-embedded loop search in which the innermost loop is skipped whenever the contribution of the first N-1 pulses to the numerator DAkT is deemed insufficient;

Figure 4 is a schematic representation of the N-nested loops used in the codebook search; and Figure 5 is a schematic block diagram illustrating the infrastructure of a typical cellular communication system.

DETAILED DESCRIPTION OF THE ILLUSTRATIVE EMBODIMENT
Figure 5 illustrates the infrastructure of a typical cellular 5 communication system 1.

Although application of the search conducting method and device according to the invention to a cellular communication system is disclosed as a non limitative example in the present specification, it should 10 be kept in mind that these method and device can be used with the same advantages in many other types of communication systems in which sound signal encoding is required.

In a cellular communication system such as 1, a telecommunications service is provided over a large geographic area by dividing that large area into a number of smaller cells. Each cell has a cellular base station 2 (Figure 5) for providing radio signalling channels, and audio and data channels.

The radio signalling channels are utilized to page mobile radio telephones (mobile transmitter/receiver units) such as 3 within the limits of the cellular base station's coverage area (cell), and to place calls to other radio telephones either inside or outside the base station's cell, or onto another network such as the Public Switched Telephone Network (PSTN) 4.

Once a radio telephone 3 has successfully placed or received a call, an audio or data channel is set up with the cellular base station 2 corresponding to the cell in which the radio telephone 3 is situated, and communication between the base station 2 and radio telephone 3 occurs over that audio or data channel. The radio telephone 3 may also receive control or timing information over the signalling channel whilst a call is in progress.

If a radio telephone 3 leaves a cell during a call and enters another cell, the radio telephone hands over the call to an available audio or data channel in the new cell. Similarly, if no call is in progress a control message is sent over the signalling channel such that the radio telephone logs onto the base station 2 associated with the new cell. In this manner mobile communication over a wide geographical area is possible.

The cellular communication system 1 further comprises a terminal 5 to control communication between the cellular base stations 2 and the Public Switched Telephone Network 4, for example during a communication between a radio telephone 3 and the PSTN 4, or between a radio telephone 3 in a first cell and a radio telephone 3 in a second cell.

Of course, a bidirectional wireless radio communication sub-system is required to establish communication between each radio telephone 3 situated in one cell and the cellular base station 2 of that cell.
Such a bidirectional wireless radio communication system typically comprises in both the radio telephone 3 and the cellular base station 2 (a) a transmitter for encoding the speech signal and for transmitting the encoded speech signal through an antenna such as 6 or 7, and (b) a receiver for receiving a transmitted encoded speech signal through the same antenna 6 or 7 and for decoding the received encoded speech signal. As well known to those of ordinary skill in the art, voice encoding is required in order to reduce the bandwidth necessary to transmit speech across the bidirectional wireless radio communication system, i.e. between a radio telephone 3 and a base station 2.
The aim of the present invention is to provide an efficient digital speech encoding technique with a good subjective quality/bit rate trade-off for example for bidirectional transmission of speech signals between a cellular base station 2 and a radio telephone 3 through an audio or data channel. Figure 1 is a schematic block diagram of a digital speech-encoding device suitable for carrying out this efficient technique.

The speech encoding device of Figure 1 is the same encoding device as illustrated in Figure 1 of U.S. parent patent No. 5 444 816 to which an amplitude selector 112 in accordance with the present invention has been added. U.S. parent patent No. 5 444 816 issued on August 22, 1995 for an invention entitled "DYNAMIC CODEBOOK FOR
EFFICIENT SPEECH CODING BASED ON ALGEBRAIC CODES".

The analog speech signal is sampled and block processed. It should be understood that the present invention is not limited to an application to speech signal. Encoding of other types of sound signal can also be contemplated.

In the illustrated example, the block of input sampled speech S (Figure 1) comprises L consecutive samples. In the CELP
literature, L is designated as the "subframe" length and is typically situated between 20 and 80. Also, the blocks of L samples are referred to as L-dimensional vectors. Various L-dimensional vectors are produced in the course of the encoding procedure. A list of these vectors, which appear in Figures 1, and 2, as well as a list of transmitted parameters is given hereinbelow:

List of the main L-dimensional vectors:
S Input speech vector;
R' Pitch-removed residual vector;
X Target vector;
D Backward-filtered target vector;
Ak Codevector of index k from the algebraic codebook; and Ck Innovation vector (filtered codevector).
List of transrnitted parameters:
k Codevector index (input of the algebraic codebook);
g Gain;
STP Short term prediction parameters (defining A(z)); and LTP Long term prediction parameters (defining a pitch gain b and a pitch delay T).

I
DECODING PRINCIPLE:

It is believed preferable to describe first the speech decoding device of Figure 2 illustrating the various steps carried out between the digital input (input of demultiplexer 205) and the output sampled speech (output of synthesis filter 204).

The demultiplexer 205 extracts four different parameters from the binary information received from a digital input channel, namely the index k, the gain g, the short term prediction parameters STP, and the long term prediction parameters LTP. The current L-dimensional vector S
of speech signal is synthesized on the basis of these four parameters as will be explained in the following description.

The speech decoding device of Figure 2 comprises a dynamic codebook 208 composed of an algebraic code generator 201 and an adaptive prefilter 202, an amplifier 206, an adder 207, a long term predictor 203, and a synthesis filter 204.

In a first step, the algebraic code generator 201 produces a codevector Ak in response to the index k.

In a second step, the codevector Ak is processed by an adaptive prefilter 202 supplied with the long term prediction parameters LTP to produce an output innovation vector Ck. The purpose of the adaptive prefilter 202 is to dynamically control the frequency content of the output innovation vector Ck so as to enhance speech quality, i.e. to reduce the audible distortion caused by frequencies annoying the human ear.

Typical transfer functions F(z) for the adaptive prefilter 202 are given below:

Fe (z) = A(z/Y,) A(Z/r2) Fb (Z) 1- bo zT

Fe(z) is a formant prefilter in which 0< Y, < yz < 1 are constants. This prefilter enhances the formant regions and works very effectively specially at coding rate below 5 kbit/s.

10 Fb(z) is a pitch prefllter where T is the time varying pitch delay and bo is either constant or equal to the quantized long-term pitch prediction parameter from the current or previous subframes. Fb(z) is.very effective to enhance pitch harmonic frequencies at all rates. Therefore, F(z) typically includes a pitch prefilter sometimes combined with a formant 15 prefilter, namely:

F(z) = Fa(z) Fc(z) In accordance with the CELP technique, the output sampled speech signal S is obtained by first scaling the innovation vector Ck from the codebook 208 by the gain g through the amplifier 206. The adder 207 then adds the scaled waveform gCk to the output E (the long term prediction component of the signal excitation of the synthesis filter 204) of a long term predictor 203 supplied with the LTP parameters, placed in a feedback loop and having a transfer function B(z) defined as follows:
B(z) = bz T

where b and T are the above defined pitch gain and delay, respectively.
The predictor 203 is a filter having a transfer function being in accordance with the last received LTP parameters b and T to model the pitch periodicity of speech. It introduces the appropriate pitch gain b and delay T of samples. The composite signal E + gCk constitutes the signal excitation of the synthesis filter 204 which has a transfer function 1/A(z) (A(z) being defined in the following description). The filter 204 provides the correct spectrum shaping in accordance with the last received STP parameters. More specifically, the filter 204 models the resonant frequencies (formants) of speech. The output block S is the synthesized sampled speech signal, which can be converted into an analog signal with proper anti-aliasing filtering in accordance with a technique well known in the art.

There are many ways to design an algebraic code generator 201. An advantageous method, disclosed in the above-mentioned U.S. patent No. 5 444 816, consists of using at least one N-interleaved single-pulse permutation code.

This concept will be illustrated by way of a simple algebraic code generator 201. In this example, L = 40 and the set of 40-dimensional codevectors contains only N = 5 non-zero-amplitude pulses that will be called Sp,, SPZ, Sp,, SP,, Sps -!n this more thorough notation, p;
stands for the location of the i"' pulse within the subframe (i.e., p; ranges from 0 to L-1). Suppose that pulse Sp, is constrained to eight possible positions p, as follows:

p, = 0,5,10,15,20,25,30,35 = 0+8m, ;m, = 0,1...7 Within these eight positions, which can be called "track" #1, SP, and seven zero-amplitude pulses can freely permute. This is a "single-pulse permutation code". Let us now interleave five such "single pulse permutation codes" by also constraining the positions of the remaining pulses in a similar fashion (i.e. track #2, track #3, track #4, and track #5).
p, = 0,5,10,15,20,25,30,35 = 0+8m, p2 = 1,6,11,16,21,26,31,36 = 1+8m2 p3 = 2,7,12,17,22,27,32,37 = 2+8m3 p4 = 3,8,13,18,23,28,33,38 = 3+8m4 p5 = 4,9,14,19,24,29,34,39 = 4+8m5 Note that the integers m, = 0, 1, ..., 7 fully define the position p, of each pulse Sp,. Thus, a simple position index kp can be derived through straightforward multiplexing of the m; s using the following relation:

kp=4096m,+512m2+64m3+8m4 +ms It should be pointed out that other codebooks can be derived using the above pulse tracks. For instance, only 4 pulses can be used, where the first three pulses occupy the positions in the first three tracks, respectively, while the fourth pulse occupies either the fourth or the fifth track with one bit to specify which track. This design gives rise to a bit position codebook.

In the prior art, the non-zero-amplitude pulses were assumed to have a fixed amplitude for all practical purposes for reasons of codevector search complexity. Indeed, if pulse gP, may assume one of q possible amplitudes, as many as q" puise-amplitude combinations will have to be considered in the search. For instance, if the five pulses of the first example are allowed to take one of q = 4 possible amplitudes, for example gP; =+ 1, -1, + 2, - 2 instead of a fixed amplitude, the algebraic codebook size jumps from 15 to 15+(5x2) bits = 25 bits; that is, a search a thousand time more complex.

The present specification discloses the surprising fact that very good performance can be achieved with q-amplitude pulses without paying a heavy price. The solution consists of limiting the search to a restrained subset of codevectors. The method of selecting the codevectors is related to the input speech signal as will be described in the following description.
A practical benefit is to enable an increase of the size of the dynamic algebraic codebook 208 by allowing individual pulses to assume different possible amplitudes without increasing the codevector search complexity.
ENCODING PRINCIPLE:
The sampled speech signal S is encoded on a block by block basis by the encoding system of Figure 1 that is broken down into 11 modules numbered from 102 to 112. The function and operation of most of these modules are unchanged with respect to the description of U.S.
parent patent No. 5 444 816. Therefore, although the following description will at least briefly explain the function and operation of each module, it will concentrate on the matter that is new with respect to the disclosure of U.S.
parent patent No. 5 444 816.

For each block of L samples of speech signal, a set of Linear Predictive Coding (LPC) parameters, called short term prediction (STP) parameters, is produced in accordance with a prior art technique through an LPC spectrum analyser 102. More specifically, the analyser 102 models the spectral characteristics of each block S of L samples.
The input block S of L-sample is whitened by a whitening filter 103 having the following transfer function based on the current values of the STP parameters:

A(z) ai z~
j__Q
where ao = 1, and z is the usual variable of the so-called z-transform. As illustrated in Figure 1, the whitening filter 103 produces a residual vector R.

A pitch extractor 104 is used to compute and quantize the LTP parameters, namely the pitch delay T and the pitch gain g. The initial state of the extractor 104 is also set to a value FS from an initial state extractor 110. A detailed procedure for computing and quantizing the LTP

parameters is described in U.S. parent patent No. 5 444 816 and is believed to be well known to those of ordinary skill in the art. Accordingly, it will not be further described in the present disclosure.

5 A filter responses characterizer 105 (Figure 1) is supplied with the STP and LTP parameters to compute a filter responses characterization FRC for use in the later steps. The FRC information consists of the following three components where n = 1, 2, ... L.

- f(n): response of F(z) Note that F(z) typically includes the pitch prefilter.

- h(n): response of /1_~l to f(n) AI zy I

\ /

where y is a perceptual factor. More generally, h(n) is the impulse response of F(z)W(z)/A(z) which is the cascade of prefilter F(z), perceptual weighting filter W(z) and synthesis fiiter 1/A(z). Note that F(z) and 1/A(z) are the same filters as used in the decoder of Figure 2.

- U(i,j): autocorrelation of h(n) according to the following expression:
u(i,j)-k' h(k-i+1)h(k- j+1) for 1 s i s L and i s j s L; h(n)=0 for n<1 The long term predictor 106 is supplied with the past excitation signal (i.e. E+ gCk of the previous subframe) for form the new E
component using proper pitch delay T and gain b.

The initial state of the perceptual filter 107 is set to the value FS supplied from the initial state extractor 110. The pitch removed residual vector R'= R-E calculated by a subtractor 121 (Figure 1) is then supplied to the perceptual filter 107 to obtain at the output of the tatter filter a target vector X. As illustrated in Figure 1, the STP parameters are applied to the filter 107 to vary its transfer function in relation to these parameters. Basically, X = R' - P where P represents the contribution of the long term prediction (LTP) including "ringing" from the past excitations.
The MSE criterion which applies to A can now be stated in the following matrix notations:

minIIAII2 = minR''-S'1I2 = minIiS'+-9AkHTjj2 = minlIX-gAkHTIl2 k where H is an L x L lower-triangular Toeplitz matrix formed from the h(n) response as follows. The term h(0) occupies the matrix diagonal and h(1), h(2), ...h(L-1) occupy the respective lower diagonals.

A backward fiitering step is performed by the filter 108 of Figure 1. Setting to zero the derivative of the above equation with respect to the gain g yields to the optimum gain as follows:

~9=0 ag g = X(Ak HT )T
ii&HTiI

With this value for g, the minimization becomes:

minllo112 = min {IIXIIZ' (x(AkHTY}

IIAkH TIlThe objective is to find the particular index k for which the minimization is achieved. Note that because IIXlf is a fixed quantity, the same index can be found by maximizing the following quantity:

max (X(Ak HT )T )2 = max ((XH) A''T )2 = max (D AkT )2 k iIAkHTII2 k akZ k ak2 where D = (XH) and qk2 = (IAkHT,I2 In the backward filter 108, a backward filtered target vector D = (XH) is computed. The term "backward filtering" for this operation comes from the interpretation of (XH) as the filtering of time-reversed X.

Only an amplitude selector 112 has been added to Figure I of the above mentioned U.S. parent patent No. 5 444 816. The function of the amplitude selector 112 is to restrain the codevectors Ak being searched by the optimizing controller 109 to the most promising codevectors Ak to thereby reduce the codevector search complexity. As described in the foregoing description, each codevector Ak is a pulse amplitude/position combination waveform defining L different positions p and comprising both zero-amplitude pulses and non-zero-amplitude pulses assigned to respective positions p = 1, 2, ...L of the combination, wherein each non-zero-amplitude pulse assumes at least one of q different possible amplitudes.

Referring now to Figure 3a, 3b and 3c, the purpose of the amplitude selector 112 is to pre-establish a function Sp between the positions p of the codevector waveform and the q possible values of the pulse amplitudes. The pre-established function Sp is derived in relation to the speech signal prior to the codebook search. More specifically, pre-establishing this function consists of pre-assigning, in relation to the speech signal, at least one of the q possible amplitudes to each position p of the waveform (step 301 of Figure 3a).

To pre-assign one of the q amplitudes to each position p of the waveform, an amplitude estimate vector B is calculated in response to the backward-filtered target vector D and to the pitch-removed residual vector R'. More specifically, the amplitude estimate vector B is calculated by summing (substep 301-1 of Figure 3b) the backward-filtered target vector D in normalized form:
(1-a)~p~

i and the pitch-removed residual vector R' in normalized form:
R' IIR'Il to thereby obtain an amplitude estimate vector B of the form:
, JIDII +R IIR'll where {3 is a fixed constant having a typical value of 1/2 (the value of R is chosen between 0 and I depending on the percentage of non-zero-amplitude pulses used in the algebraic code).

For each position p of the waveform, the amplitude Sp to be pre-assigned to that position p is obtained by quantizing a corresponding amplitude estimate Bp of vector B. More specifically, for each position p of the waveform, a peak-normalized amplitude estimate Bp of the vector B is quantized (substep 301-2 of Figure 3b) using the following expression:

Sp=Q(BP/ maxIB.{) n wherein Q() is the quantization function and max ( Bn I
n is a normalisation factor representing a peak amplitude of the non-zero-amplitude pulses.

In the important special case in which:

- q = 2, that is the pulse amplitudes can assume only two values (i.e. ~~ = t 1); and - the non-zero-amplitude pulse density N/L is lower than 10 or equal to 15%;

- the value of 0 can be equal to zero; then the amplitude estimate vector B reduces simply to the backward-filtered target vector D
and consequently Sp = sign(DP) .

The purpose of the optimizing controller 109 is to select the best codevector Ak from the algebraic codebook. The selection criterion is given in the form of a ration to be calculated for each codevector Ak and to be maximized over all codevectors (step 303):

max (D Akr )z k ak2 where D=(XH) and ak = IIAkHTIr =
Since Ak is an algebraic codevector having N non-zero-amplitude pulses of respective amplitudes Sp,, the numerator is the square of N
D Ak = iEDp,Sp, and the denominator is an energy term which can be expressed as:

ak = ~ESp,U (Pi,pi)+ ~~,Sp,Sp, U (Pi,P

where U(p,,p,) is the correlation associated with two unit-amplitude pulses, one at position p, and the other at position pj. This matrix is computed in accordance with the above equation in the filter response characterizer 105 and included in the set of parameters referred to as FRC in the block diagram of Figure 1.

A fast method for computing this denominator (step 304) involves the N-nested loops illustrated in Figure 4 in which the trim lined notation S(i) and SS(i,j) is used in the place of the respective quantities "Sp," and "Sp,Sp,"= Computation of the denominator ak2 is the most time consuming process. The computations contributing to akz which are performed in each loop of Figure 4 can be written on separate lines from the outermost loop to the innermost loop as follows:

ak = Sp,ZU(P,'P,) + SPzZU(P2P2)+ 2SP,SPzU(P,I P2) + 'SPa2 u(P9 ~ Ps) + 2[SP,SP,U(p, IPs) +SpzSP,U(P2 ' P3 )l + SPN2 U (PN ,PN ) + 2[Sp,SpNV (P1 RPN ) *SpzSPNU (P2 IpN ) + ... i'SPN.,SPNU
(PN-t IPN)J

where p; is the position of the i"' non-zero-amplitude pulse. Note that the N-nested loops of Figure 4 enables constraining the non-zero-amplitude pulses of codevectors Ak in accordance with N interleaved single-pulse permutation codes.

In the present invention search complexity is drastically reduced by restraining the subset of codevectors Akbeing searched to codevectors of which the N non-zero-amplitude pulses respect the function pre-established in step 301 of Figure 3a. The pre-established function is respected when the N non-zero-amplitude pulses of a codevector Akeach have an amplitude equal to the amplitude pre-assigned to the position p of the non-zero-amplitude pulse.
Said restraining the subset of codevectors is preformed by first combining the pre-established function SP with the entries of matrix U(i,j) (step 302 of Figure 3a) then, by using the N-nested loops of Figure 4 with all pulses S(i) assumed to be fixed, positive and of unit amplitude (step 303). Thus, even though the amplitude of non-zero pulses can take any of q possible values in the aigebraic codebook, the search complexity is reduced to the case of fixed pulse amplitudes. More precisely, the matrix U(i,j) which is supplied by the filter response characterizer 105 is combined with the pre-established function in accordance with the following relation (step 302):

U'(i.j) = Ss S; U (i,j) where S, results from the selecting method of amplitude selector 112, namely S, is the amplitude selected for an individual position i following quantization of the corresponding amplitude estimate.

With this new matrix, the computation for each loop of the fast algorithm can be written on a separate line, from the outermost to the innermost loop, as follows:

ak = U'(p,-P,) + U'(ps-p2)+ 2U'(prpz) + U'(ps.pa)+ 2U'(p,p3)+ 2U'(p2,p3) + U'(pN,PN)+ 2U'(p,,pN)+ 2U'(p2'P,)+..,+ 2U'(pN-,IPN) where pX is the position of the x'h non-zero-amplitude pulse of the waveform, and where U'(px,py) is a function dependent on the amplitude Spx pre-assigned to a position px amongst the positions p and the amplitude Spy pre-assigned to a position py amongst the positions p.
To still further reduce the search complexity, one may skip (cf Figure 3c) in particular, but not exclusively, the innermost loop whenever the following inequality is true:

W-, ~ SP. Dp~ < To where SP~ is the amplitude pre-assigned to position põ, Dp. is the p,"
component of the target vector D, and To is a threshold related to the backward-filtered target vector D.

The global signal excitation signal E + gCk is computed by an adder 120 (Figure 1) from the signal gCk from the controller 109 and the output E from the predictor 106. The initial state extractor module 110, constituted by a perceptual filter with a transfer function 1/A(zy'') varying in relation to the STP parameters, subtracts from the residual signal R the signal excitation signal E + gCk for the sole purpose of obtaining the final filter state FS for use as initial state in filter 107 and pitch extractor 104.
The set of four parameters k, g, LTP and STP are converted into the proper digital channel format by a multiplexer 111 completing the procedure for encoding a block S of samples of speech signal.

Although the present invention has been described hereinabove with reference to an illustrative embodiment thereof, this embodiment can be modified at will, within the scope of the appended claims, without departing from the spirit and nature of the subject invention.

Claims (20)

WHAT IS CLAIMED IS:
1. A method of conducting a search in a codebook during encoding of a sound signal, in which:
- the codebook contains a set of pulse amplitude/position combinations;
- each pulse amplitude/position combination defines a number L of different positions p and comprises both zero-amplitude pulses and non-zero-amplitude pulses assigned to respective positions p = 1, 2, ...L of the pulse amplitude/position combination;
- each non-zero-amplitude pulse assumes one of q possible amplitudes;
and - said codebook search conducting method comprises:
pre-selecting from said codebook a subset of pulse amplitude/position combinations in relation to the sound signal; and searching only said subset of pulse amplitude/position combinations during encoding of the sound signal whereby complexity of the search is reduced as only a subset of the pulse amplitude/position combinations of the codebook is searched;
wherein:
pre-selecting a subset of pulse amplitude/position combinations comprises pre-establishing, in relation to the sound signal, an amplitude/position function S p between the positions p = 1, 2, ...L and the q possible amplitudes;
pre-establishing an amplitude/position function S p comprises pre-assigning one of the q possible amplitudes as valid amplitude to each position p; and pre-assigning one of the q possible amplitudes to each position p comprises:

processing the sound signal to produce a backward-filtered target signal D and a pitch-removed residual signal R';
calculating an amplitude estimate vector B in response to the backward-filtered target signal D and to the pitch-removed residual signal R'; and for each of said positions p, quantizing an amplitude estimate B p of said vector B to obtain the amplitude to be selected for said position p; and searching said subset of pulse amplitude/position combinations comprises limiting the search to the pulse amplitude/position combinations of said codebook having non-zero-amplitude pulses which satisfy the pre-established function S p.
2. ~The method of claim 1, wherein the pre-established function S p is satisfied when the non-zero-amplitude pulses of a pulse amplitude/position combination each have an amplitude equal to the amplitude pre-assigned by the pre-established function S p to the position p of said non-zero-amplitude pulse.
3. ~The method of claim 1 or 2, in which calculating the amplitude estimate vector B comprises summing the backward-filtered target signal D
in normalized form:

to the pitch-removed residual signal R' in normalized form:

to thereby obtain the amplitude estimate vector B of the form:

where .beta. is a fixed constant.
4. ~The method of claim 3, wherein .beta. is a fixed constant having a value situated between 0 and 1.
5. ~The method of one of claims 1 to 4, in which for each of said positions p, the quantizing of the amplitude vector estimate comprises quantizing a peak-normalized amplitude estimate B p of said vector B using the following expression:

B p/max|B n|
n wherein the denominator max|B n|
n is a normalizing factor representing a peak amplitude of the non-zero-amplitude pulses.
6. ~The method of claim 1, further comprising restraining the positions p of the non-zero-amplitude pulses of the combinations of the codebook in accordance with a set of tracks of pulse positions.
7. ~The method of one of claims 1 to 6, wherein said pulse amplitude/position combinations each comprise a number N of non-zero-amplitude pulses, and wherein searching said subset of pulse amplitude/position combinations comprises maximizing a given ratio having a denominator .alpha. k2 computed by means of N nested loops in accordance with the following relation:

.alpha. k 2 = U'(p1,p1) + U'(p2,p2) + 2U'(p1,p2) + U'(p3,p3) + 2U'(p1,p3) + 2U'(p2,p3) + U'(p N,p N)+ 2U'(p1,p N)+ 2U'(P2,p N)+...+ 2U'(p N-1,p N) where computation for each loop is written in a separate line from an outermost loop to an innermost loop of the N nested loops, where p n is the position of the n th non-zero-amplitude pulse of the combination, and where U'(p x,p y) is a function dependent on the amplitude S p x pre-assigned to a position p x amongst the positions p and the amplitude S p y pre-assigned to a position p y amongst the positions p.
8. ~The method of claim 7, wherein maximizing said given ratio comprises skipping at least the innermost loop of the N nested loops whenever the following inequality is true where S p n is the amplitude pre-assigned to position p n, D p n is the p n th component of the target vector D, and T D is a threshold related to the backward-filtered target vector D.
9. ~A device for conducting a search in a codebook during encoding of a sound signal, in which:
- the codebook contains a set of pulse amplitude/position combinations;
- each pulse amplitude/position combination defines a number L of different positions p and comprises both zero-amplitude pulses and non-zero-amplitude pulses assigned to respective positions p = 1, 2, ...L of the pulse amplitude/position combination;

- each non-zero-amplitude pulse assumes one of q possible amplitudes;
and - said codebook search conducting device comprises:
means for pre-selecting from said codebook a subset of pulse amplitude/position combinations in relation to the sound signal; and means for searching only said subset of pulse amplitude/position combinations during encoding of the sound signal whereby complexity of the search is reduced as only a subset of the pulse amplitude/position combinations of the codebook is searched;
wherein:
the pre-selecting means comprises means for pre-establishing, in relation to the sound signal, an amplitude/position function S p between the positions p = 1, 2, ...L and the q possible amplitudes;
the pre-establishing means comprises means for pre-assigning one of the q possible amplitudes as valid amplitude to each position p;
and the means for pre-assigning one of the q possible amplitudes to each position p comprises:
means for processing the sound signal to produce a backward-filtered target signal D and a pitch-removed residual signal R';
means for calculating an amplitude estimate vector B in response to the backward-filtered target signal D and to the pitch-removed residual signal R'; and means for quantizing, for each of said positions p, an amplitude estimate B p of said vector B to obtain the amplitude to be selected for said position p; and the searching means comprises means for limiting the search to the pulse amplitude/position combinations of said codebook having non-zero-amplitude pulses which satisfy the pre-established function S p.
10. ~The device of claim 9, wherein the pre-established function S p is satisfied when the non-zero-amplitude pulses of a pulse amplitude/position combination each have an amplitude equal to the amplitude pre-assigned by the pre-established function S p to the position p of said non-zero-amplitude pulse.
11. ~The device of claim 9 or 10, in which the means for calculating the amplitude estimate vector B comprises means for summing the backward-filtered target signal D in normalized form:

to the pitch-removed residual signal R' in normalized form:

to thereby obtain the amplitude estimate vector B of the form:

where .beta. is a fixed constant.
12. ~The device of claim 11, wherein .beta. is a fixed constant having a value situated between 0 and 1.
13. ~The device of one of claims 9 to 12, in which the means for quantizing the amplitude vector estimate comprises means for quantizing, for each of said positions p, a peak-normalized amplitude estimate B p of said vector B using the following expression:

B p/max | B n |
n wherein the denominator max | B n |
n is a normalizing factor representing a peak amplitude of the non-zero-amplitude pulses.
14. ~The device of claim 9, further comprising means restraining the positions p of the non-zero-amplitude pulses of the combinations of the codebook in accordance with a set of tracks of pulse positions.
15. ~The device of one of claims 9 to 14, wherein said pulse amplitude/position combinations each comprise a number N of non-zero-amplitude pulses, and wherein the means for searching said subset of pulse amplitude/position combinations comprises means for maximizing a given ratio having a denominator a k 2 computed by means of N nested loops in accordance with the following relation:

.alpha. k 2 = U'(p1,p1) + U'(p2,p2)+ 2U'(p1,p2) + U'(p3,p3) + 2U'(p1,p3)+ 2U'(p2,p3) + U(p N,p N)+ 2U'(p1,p N)+ 2U'(p2,p N)+...+ 2U'(p N-1,p N) where computation for each loop is written in a separate line from an outermost loop to an innermost loop of the N nested loops, where p n is the position of the n th non-zero-amplitude pulse of the combination, and where U'(p x,p y) is a function dependent on the amplitude S p x pre-assigned to a position p x amongst the positions p and the amplitude S p y pre-assigned to a position p y amongst the positions p.
16. ~The device of claim 15, wherein the means for maximizing said given ratio comprises means for skipping at least the innermost loop of the N nested loops whenever the following inequality is true where S p n is the amplitude pre-assigned to position p n, D p n is the p n th component of the target vector D, and T D is a threshold related to the backward-filtered target vector D.
17. ~A cellular communication system for servicing a large geographical area divided into a plurality of cells, comprising:
mobile transmitter/receiver units;
cellular base stations respectively situated in said cells;
means for controlling communication between the cellular base stations;
a bidirectional wireless communication sub-system between each mobile unit situated in one cell and the cellular base station of said one cell, said bidirectional wireless communication sub-system comprising in both the mobile unit and the cellular base station (a) a transmitter including means for encoding a speech signal and means for transmitting the encoded speech signal, and (b) a receiver including means for receiving a transmitted encoded speech signal and means for decoding the received encoded speech signal;
-wherein said speech signal encoding means comprises means responsive to the speech signal for producing speech signal encoding parameters, and wherein said speech signal encoding parameter producing means comprises a device as recited in any of claims 9 to 16, for conducting a search in a codebook in view of producing at least one of said speech signal encoding parameters, wherein the speech signal constitutes said sound signal.
18. ~A cellular network element comprising (a) a transmitter including means for encoding a speech signal and means for transmitting the encoded speech signal, and (b) a receiver including means for receiving a transmitted encoded speech signal and means for decoding the received encoded speech signal;
-wherein said speech signal encoding means comprises means responsive to the speech signal for producing speech signal encoding parameters, and wherein said speech signal encoding parameter producing means comprises a device as recited in any of claims 9 to 16, for conducting a search in a codebook in view of producing at least one of said speech signal encoding parameters, wherein the speech signal constitutes said sound signal.
19. ~A cellular mobile transmitter/receiver unit comprising (a) a transmitter including means for encoding a speech signal and means for transmitting the encoded speech signal, and (b) a receiver including means for receiving a transmitted encoded speech signal and means for decoding the received encoded speech signal;
-wherein said speech signal encoding means comprises means responsive to the speech signal for producing speech signal encoding parameters, and wherein said speech signal encoding parameter producing means comprises a device as recited in any of claims 9 to 16, for conducting a search in a codebook in view of producing at least one of said speech signal encoding parameters, wherein the speech signal constitutes said sound signal.
20. ~In a cellular communication system for servicing a large geographical area divided into a plurality of cells, and comprising: mobile transmitter/receiver units; cellular base stations respectively situated in said cells; and means for controlling communication between the cellular base stations;
- a bidirectional wireless communication sub-system between each mobile unit situated in one cell and the cellular base station of said one cell, said bidirectional wireless communication sub-system comprising in both the mobile unit and the cellular base station (a) a transmitter including means for encoding a speech signal and means for transmitting the encoded speech signal, and (b) a receiver including means for receiving a transmitted encoded speech signal and means for decoding the received encoded speech signal;
-wherein said speech signal encoding means comprises means responsive to the speech signal for producing speech signal encoding parameters, and wherein said speech signal encoding parameter producing means comprises a device as recited in any of claims 9 to 16, for conducting a search in a codebook in view of producing at least one of said speech signal encoding parameters, wherein the speech signal constitutes said sound signal.
CA002618002A 1995-02-06 1996-02-02 Algebraic codebook with signal-selected pulse amplitudes for fast coding of speech Expired - Lifetime CA2618002C (en)

Applications Claiming Priority (5)

Application Number Priority Date Filing Date Title
US38396895A 1995-02-06 1995-02-06
US08/383,968 1995-02-06
US08/508,801 1995-07-25
US08/508,801 US5754976A (en) 1990-02-23 1995-07-28 Algebraic codebook with signal-selected pulse amplitude/position combinations for fast coding of speech
PCT/CA1996/000069 WO1996024925A1 (en) 1995-02-06 1996-02-02 Algebraic codebook with signal-selected pulse amplitudes for fast coding of speech

Publications (2)

Publication Number Publication Date
CA2618002A1 CA2618002A1 (en) 1997-08-03
CA2618002C true CA2618002C (en) 2001-08-21

Family

ID=39264252

Family Applications (1)

Application Number Title Priority Date Filing Date
CA002618002A Expired - Lifetime CA2618002C (en) 1995-02-06 1996-02-02 Algebraic codebook with signal-selected pulse amplitudes for fast coding of speech

Country Status (1)

Country Link
CA (1) CA2618002C (en)

Also Published As

Publication number Publication date
CA2618002A1 (en) 1997-08-03

Similar Documents

Publication Publication Date Title
EP1225568B1 (en) Algebraic codebook with signal-selected pulse amplitudes for fast coding of speech
KR100299408B1 (en) Depth-first algebraic codebook search for fast coding of speech
US7774200B2 (en) Method and apparatus for transmitting an encoded speech signal
US5729655A (en) Method and apparatus for speech compression using multi-mode code excited linear predictive coding
US5359696A (en) Digital speech coder having improved sub-sample resolution long-term predictor
KR19980080463A (en) Vector quantization method in code-excited linear predictive speech coder
KR20020077389A (en) Indexing pulse positions and signs in algebraic codebooks for coding of wideband signals
EP0801377B1 (en) Apparatus for coding a signal
CA2618002C (en) Algebraic codebook with signal-selected pulse amplitudes for fast coding of speech
CA2210765E (en) Algebraic codebook with signal-selected pulse amplitudes for fast coding of speech
US5943644A (en) Speech compression coding with discrete cosine transformation of stochastic elements
JP2002073097A (en) Celp type voice coding device and celp type voice decoding device as well as voice encoding method and voice decoding method
KR100389898B1 (en) Method for quantizing linear spectrum pair coefficient in coding voice
NO322594B1 (en) Algebraic codebook with signal-selected pulse amplitudes for fast speech encoding
MXPA97005997A (en) Algebraic coding book with amplitude signal deimpulse selected for a rapidacodification of
SE524905C2 (en) Pulse combination codebook searching method for sound signal encoding - preselecting sound related pulse combination subset and searching within that subset

Legal Events

Date Code Title Description
EEER Examination request
MKEX Expiry

Effective date: 20160202