EP0135512A4 - Improved method and means of determining coefficients for linear predictive coding. - Google Patents

Improved method and means of determining coefficients for linear predictive coding.

Info

Publication number
EP0135512A4
EP0135512A4 EP19840900575 EP84900575A EP0135512A4 EP 0135512 A4 EP0135512 A4 EP 0135512A4 EP 19840900575 EP19840900575 EP 19840900575 EP 84900575 A EP84900575 A EP 84900575A EP 0135512 A4 EP0135512 A4 EP 0135512A4
Authority
EP
European Patent Office
Prior art keywords
coefficients
signal
ram
sequencer
lattice
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
EP19840900575
Other languages
German (de)
French (fr)
Other versions
EP0135512B1 (en
EP0135512A1 (en
Inventor
Ira Alan Gerson
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.)
Motorola Solutions Inc
Original Assignee
Motorola Inc
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
Application filed by Motorola Inc filed Critical Motorola Inc
Priority to AT84900575T priority Critical patent/ATE36206T1/en
Publication of EP0135512A1 publication Critical patent/EP0135512A1/en
Publication of EP0135512A4 publication Critical patent/EP0135512A4/en
Application granted granted Critical
Publication of EP0135512B1 publication Critical patent/EP0135512B1/en
Expired 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

Definitions

  • This invention relates to linear predictive coding.
  • it is an improved method and means of determining coefficients for linear predictive coding.
  • Linear predictive coding is a method of analyzing a speech signal and characterizing that signal in terms of coefficients which can be encoded, broadcast, received, and decoded to recover a close approximation to the original signal.
  • the existence of redundancies in speech makes it possible to use encoded descriptions of the speech that can be carried in a communication channel having a bandwidth that is less than the bandwidth of the speech. This is in distinct contrast to many well-known forms of converting speech into digital signals. Most of these methods require a bandwidth that is greater than the bandwidth of the speech.
  • Linear predictive coding (LPC) of speech begins conceptually with a model of the human speech-producing system.
  • the model has a source of sound that is analogous to the vocal cords. That source is coupled acoustically to a stacked array of hollow cylindrical tubes that is analogous to the cavities of the throat and mouth in a human speaker.
  • speech is characterized by four types of quantities. The first of these is a measure of whether speech is voiced or unvoiced.
  • a voiced signal begins with an input from the vocal cords, while an unvoiced signal is produced by the action of the rest of the vocal tract on moving air alone. This produces the differences in sound between "s" which is unvoiced and its voiced equivalent "z", for example.
  • a second characteristic of the sound is the pitch, the fundamental frequency produced by the vocal cords in making a voiced sound.
  • a third characteristic is the energy.
  • the effect of the throat and mouth on either voiced or unvoiced sound is summarized by obtaining some measure of the transfer function of the vocal tract. Such a measure might be the reflection coefficients of the structure, the poles of the transfer function of the structure, the logarithmic area ratios (LAR's) of the structure, or any of several other well-known functions of such resonances.
  • various mathematical transforms of these functions may have utility for particular purposes. The functions are interrelated so that any one set can be determined from a knowledge of any other set.
  • the present invention is an improved method and means of determining quantities corresponding to reflection coefficients.
  • the reflection coefficients are coefficients with a specific all pole filter structure known as a lattice filter.
  • the reflection coefficients discussed here are the electrical analog of acoustical reflection coefficients. While the reflection coefficient is in general a complex quantity in which the imaginary portion is a measure of loss at a discontinuity, when determining the reflection coefficients that characterize a human vocal tract, the portions of reflection coefficients associated with loss are small enough compared to the real part that denotes lossless reflection that it is adequate to use a lossless model in which the reflection coefficients are real.
  • One way of obtaining the reflection coefficients that characterize a sample of a signal such as speech is to determine the characteristics of an inverse lattice filter that will reproduce the signal when excited by an impulse. Examples of several approaches of performing such an analysis are given in an article entitled “Stable and Efficient Lattice Methods for Linear Prediction” by Makhoul in “IEEE Transactions on Acoustics, Speech, and Signal Processing, Vol. ASSP-25, No. 5, October, 1977. This article points out that an analysis based on the all-pole lattice filter is stable without windowing but at a computational cost that is several times the cost of the auto-correlation and covariance methods of calculation. What Makhoul refers to as the computational costs of computer analysis are proportional to the equipment costs of realizing a circuit for comparable analysis on a semiconductor chip.
  • An improved method and means of determining reflection coefficients that characterize an electrical signal obtains characteristics of an all-zero inverse lattice filter.
  • the reflection coefficients are obtained by filtering the signal, sampling the filtered signal. obtaining the elements of a correlation matrix from the samples, initializing values of matrices of forward residual autocorrelations, backward residual autocorrelations, and cross correlation of residuals, combining matrix elements to obtain a first reflection coefficient, removing from the forward, backward and cross-correlation matrices the effect of the first reflection coefficient, calculating from the revised matrices a second reflection coefficient, and repeating the calculations to the desired order.
  • Fig. 1 is an overall block diagram of a circuit for the practice of the present invention.
  • Fig. 2 is an expanded circuit diagram of the address calculator of Fig. 1.
  • Fig. 3 is an expanded circuit diagram of the update calculator of Fig. 1.
  • Fig. 4 is an expanded circuit diagram of the reflection-coefficient calculator of Fig. 1.
  • Fig. 5 is a flowchart of a method for practicing the invention.
  • Fig. 6 is a block diagram of a section of a typical lattice filter for the practice of the present invention.
  • Fig. 1 is an overall block diagram of a circuit for the practice of the invention.
  • an electrical analog signal that is to be analyzed is applied at terminal 10. That signal will typically be an electrical analog of a voice signal, although it may be any electrical signal that exhibits redundancies analogous to those of speech. Examples of such other signals include video scans and seismic analysis records.
  • the signal at terminal 10 is applied to a filter 12, if necessary, to limit its bandwidth. If the bandwidth is already adequate, filter 12 may comprise a direct wire connection, or filter 12 may combine a bandpass filter with any of a number of systems of pre-emphasis that are commonly used in radio broadcasting.
  • ADC analog-to-digital converter
  • the digital output of ADC 14 will be separated into frames of a convenient length, of the order of tens of milliseconds. That function is here indicated as being performed by framer 16 as a means of insuring that the digital input to correlator 18 establishes correlation among samples in the same frame.
  • the function of framer 16 could also be combined into correlator 18 or ADC 14. It should be noted that the ADC 14 is not necessary if the signal is already in digital form.
  • the correlation parameters are obtained in the circuit of Fig. 1 because that circuit will be used to determine reflection coefficients.
  • the output of correlator 18 is applied, if necessary, to a normalizer 20 to normalize output values to a common level.
  • the output of normalizer 20 is taken to random-access memory (RAM) 22 where it is stored in an address that is directed by a signal from multiplexer 23.
  • Normalizer 20 also generates a signal indicating the completion of the operation of correlator 18 for one frame. That signal is taken as an input to sequencer 24. Signals from sequencer 24 are coupled out to control an address calculator 26, an update calculator 28, a reflection coefficient calculator 30, multiplexer 23, and multiplexer 32. Normalizer 20 determines the appropriate addresses in RAM 22 for storing the coefficients of a correlation matrix so that their serial readout will be in a desired order.
  • the output of normalizer 20 is applied to multiplexer 23 to apply first an initial condition to RAM 22 that is determined from normalizer 20, and then accessed by address instructions from address calculator 26.
  • the output of RAM 22 is applied through multiplexer 32 to RAM 34.
  • Multiplexer 32 selects as an input to RAM 34 either the output of RAM 22 or the output of update calculator 28.
  • Multiplexer 32 is controlled by a signal from sequencer 24.
  • the location of storage elements in RAM 34 is controlled by a signal from address calculator 26.
  • the output of RAM 34 is read into update calculator 26 and reflection coefficient calculator 30.
  • update calculator 28 also receives as an input the calculated reflection coefficients from reflection coefficient calculator 30, which coefficients are taken as the output of the circuit of Fig. 1.
  • Fig. 2 is a circuit diagram of the address calculator 26 of Fig. 1.
  • Address calculator 26 produces the addresses that direct the retrieval of matrix elements in RAMs 22 and the storage and retrieval of matrix elements in RAM 34.
  • sequencer 24 supplies addresses to read-only memory (ROM) 40.
  • ROM 40 is read out into multiplexers 42 and 44.
  • the output of multiplexer 42 is taken to RAM 46 which is loaded under the control of signals from sequencer 24.
  • the output of RAM 46 is taken directly as the input to RAM 34 of Fig. 1, and it is also applied to latches 48 and 50, as well as the add/subtract unit 42.
  • the output of latch 48 is taken to multiplexer 23 of Fig. 1, and then to RAM 22.
  • the output of latch 50 is multiplexed in multiplexer 44 with the output of ROM 40 under the control of a signal from sequencer 24.
  • the output of multiplexer 44 is applied to add-subtract unit 52 where it is supplied as an input to be either added or subtracted from the output of RAM 46, depending upon a control signal from sequencer 24.
  • the add/subtract unit 52 also has an output connected to sequencer 24 which indicates a zero result from the add/subtract operation.
  • Fig. 3 is a circuit diagram of the update calculator 28 of Fig. 1.
  • data inputs are applied at terminals 60 and 62.
  • the input at terminal 60 is the output of RAM 34 of Fig. 1.
  • Terminal 60 is connected to a shifter 64, the output of which is taken to a RAM 66 and to a multiplexer 68.
  • the output of RAM 66 is taken to multiplexer 70 which, in turn, is connected to multiplier 72.
  • the output of multiplier 72 is taken as an input to multiplexer 68, as an input to register 74, and as an input to multiplexer 76.
  • the output of register 74 is coupled through multiplexer 78 to supply a second input to multiplier 72.
  • Terminal 62 is connected to reflection-coefficient calculator 30 of Fig.
  • OR gate 90 The output of OR gate 90 is taken to register 92 which supplies an output that is taken both as a second input to OR gate 90 and as an input to priority encoder 94.
  • the output of priority encoder 94 controls register 96 that operates shifter 64.
  • the combination of shifter 64, register 96, OR gate 90, register 92 and priority encoder 94 comprises a normalizer that normalizes the output of RAM 34 of Fig. 1. If it is not desired or is considered unnecessary to normalize coefficients, then terminal 60 can supply an input directly to RAM 66 and the elements just described could be removed from the circuit of Fig. 3.
  • Fig. 4 is a circuit diagram of the reflection coefficient calculator 30 of Fig. 1. In Fig.
  • 4 terminal 100 receives a signal from RAM 34 that is applied to a subtractor 102 and a multiply-by-two circuit 104.
  • the output of subtractor 102 is connected to register 106, which produces one output that is connected as a input to subtractor 102, and another output that is taken to divider 108.
  • the output of divider 108 may be taken directly as a reflection coefficient.
  • the signal at terminal 112 is a reflection coefficient of the equivalent filter of the original input signal.
  • the reflection coefficients may be used directly as the coefficients in LPC or they may be transformed into a different function, as described earlier.
  • Fig. 1 If it is desired to convert reflection coefficients into LAR's, this can be combined readily with the quantization that is performed by quantizer 110.
  • the choice of the particular functions is one of design.
  • Filter 12 has already been described as optional. If the circuit of Fig. 1 is to be applied to speech for use in a typical radio speech bandwidth, then the input signal may need to be subjected to bandpass filtering. It may also be desirable to combine some form of pre-emphasis with filter 12. Filter 12 may also be used to prevent aliasing when the filtered signal is applied to ADC 14.
  • the sample rate produced by ADC 14 is a design parameter.
  • sampling theorem it will be necessary to sample the input at a rate that is at least twice the frequency of the highest component contained in the input.
  • sampling rates for speech are typically between 6.4 kHz and 10 kHz.
  • Each sample is then encoded into a number of bits that typically ranges between 8 and 15, with 12 as a typical number. It follows that a typical frame will have of the order of one hundred samples.
  • Correlator 18 takes these samples and determines from them the elements of a covariance matrix. This matrix is symmetric.
  • p is the prediction order, a number that is typically between 8 and 12 for speech.
  • the index k is kept less than or equal to i to avoid recalculating equal terms on both sides of the axis of the symmetric matrix.
  • Normalizer 20 of Fig. 1 performs a function that is here indicated separately but that might also be included in correlator 18. Normalizer 20 shifts the elements of the correlation matrix so that the maximum value of any element in the array is between one-half and one in magnitude. Normalizer 20 then truncates the values of the elements thus shifted to a number of bits equal to the word length of the system.
  • the covariance matrix thus has a number of elements equal to (p + 1) 2 .
  • the covariance matrix is symmetric, it can be described completely by storing the elements of the diagonal and the elements below the diagonal, a total of (p + 1)(p + 2)/2 elements. These elements are stored in RAM 22 in a location that is controlled by address calculator 26. A convenient method of loading RAM 22 is to load the diagonal elements, beginning with the element of highest order and proceeding to the diagonal element of lowest order, and then repeat in sequence down paths in the matrix that are parallel to the main diagonal. It should be noted that the F and B matrices are also symmetric and may be stored in a similar fashion.
  • an operations block 122 directs the determination of correlation coefficients. Such a determination is well known. For discrete or sampled components, it is normally accomplished by a calculation such as that of the following equation:
  • the next step in the flowchart of Fig. 5 is to initialize the matrices of F, B and C as indicated in operations block 126.
  • the quantities F, B and C are intermediate quantities used in the determination of LPC coefficients. Their initial values are determined as follows:
  • Operations box 1 28 next directs that the value of j equal 1 .
  • Operations box 1 30 then calls for the determination of the quantity k j , the jth reflection coefficient. This is determined as follows:
  • Bj (i,k) B j-1 (i+l,k+1)+k j [C j-1 (i+l,k+1)+C j-1 (k+1,i+1)]+
  • the value of j is increased by 1 in operations block 140, and control returns to operations block 130 to continue the calculations.
  • Fig. 6 is a block diagram of a lattice filter that provides a further explanation of the process by which the circuit of Fig. 1 obtains reflection coefficients.
  • a terminal 150 receives as an input the sampled signal. This signal is applied to an upper leg 152 which will calculate a forward residual, and it is applied to a lower leg 154 which will apply the signal to the first of a sequence of delay elements 156 to calculate a backward residual.
  • Fig. 6 comprises a cascade of elements, each of which applies the forward signal to a multiplier 158 and a summer 165. The backward signals are applied to a multiplier 162 and a summer 164. Both multipliers 158 and 162 receive as additional inputs the current reflection coefficient. Thus, the current forward residual is multiplied by the current reflection coefficient and added to the current backward residual in summer 164 to generate as an output the next backward residual.
  • the current backward residual is multiplied in multiplier 162 with the current reflection coefficient and added to the current forward residual in summer 160 to generate the next forward residual.
  • the process just described continues through a number of sections of the lattice filter that is determined by the designer as a number adequate to characterize the particular signal in question. This is typically a number of stages equal to 8, 10 or 12. The last such stage is shown here as receiving a current forward residual signal on terminal 166 and a current backward residual signal on terminal 168.
  • the current forward residual signal is applied to a multiplier 170 and a summer 172, while the current backward residual is delayed in delay element 174, and the delayed signal is applied to multiplier 176 and summer 178. Both multipliers 170 and 176 receive as additional inputs the current reflection coefficient.
  • the result of applying a signal to this circuit of Fig. 6 is the production in forward line 152 of a sequence of elements of a forward residual vector and to produce in line 156 the elements of a backward residual vector. Individual elements are combined to produce an autocorrelation matrix of the forward residual elements, and autocorrelation matrix of the backward residual elements, and a cross-correlation matrix between forward and backward residual elements. These matrices are used as described earlier to calculate the reflection coefficients.
  • Lattice methods of determining coefficients for linear predictive coding have been used in the past.
  • circuits and programs used to determine the lattice coefficients have used intermediate variables that varied in magnitude over a wide range. This required a wide range of quantize values to characterize the intermediate variables, and thus took more computational time.
  • the circuit arrangement and method of the present invention uses only variable and intermediate variables which are bounded in magnitude by unity. This permits operations and calculations to be performed in a fixed-point fractional implementation.
  • the input signal is windowed so that it is stationary in a statistical sense, then it can be shown that the number of computations necessary. to determine lattice coeficients is reduced still further.
  • the method and means of the present invention has been used with a frame length of approximately 15 milloseconds to determine 12 lattice coefficients in real time. By completing the calculations for the data of one frame before the end of the next succeeding frame. What is claimed is:
  • NP MIN0(NP,KK-1) .
  • AK( II ) QU(AK( II) , II,AR) C SEE IF DONE
  • FR(I1,I2) FR(I1,I2)+ 1Z*(CR(I1,I2)+CR(I2,I1))+ZS*BR(I1,I2)
  • BR(I1,I2) BR(I1+1,I2+1)+ 1Z*(CR(I1+1,I2+1)+CR(I2+1,I1+1))+ZS*FR(I1+1,I2+1)

Landscapes

  • Engineering & Computer Science (AREA)
  • Multimedia (AREA)
  • Health & Medical Sciences (AREA)
  • Audiology, Speech & Language Pathology (AREA)
  • Computational Linguistics (AREA)
  • Physics & Mathematics (AREA)
  • Signal Processing (AREA)
  • Acoustics & Sound (AREA)
  • Human Computer Interaction (AREA)
  • Compression, Expansion, Code Conversion, And Decoders (AREA)
  • Medicines Containing Material From Animals Or Micro-Organisms (AREA)
  • Investigating Or Analysing Biological Materials (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
  • Filters That Use Time-Delay Elements (AREA)
  • Radar Systems Or Details Thereof (AREA)
  • Transmission Systems Not Characterized By The Medium Used For Transmission (AREA)

Abstract

An improved method and means of determining reflection coefficients that characterize an electrical signal obtains characteristics of an all-zero inverse lattice filter. The reflection coefficients are obtained by filtering the signal, sampling the filtered signal, obtaining the elements of a correlation matrix from the samples, initializing values of matrices of forward residuals, backward residuals, and cross correlation of residuals, combining matrix elements to obtain a first reflection coefficient, removing from the forward, backward an cross-correlation matrices the effect of the first reflection coefficient, calculating from the revised matrices a second reflection coefficient, and repeating the calculations to the desired order.

Description

IMPROVED METHOD AND MEANS OF
DETERMINING COEFFICIENTS FOR
LINEAR PREDICTIVE CODING
Background of the Invention
This invention relates to linear predictive coding. In particular, it is an improved method and means of determining coefficients for linear predictive coding.
Linear predictive coding is a method of analyzing a speech signal and characterizing that signal in terms of coefficients which can be encoded, broadcast, received, and decoded to recover a close approximation to the original signal. The existence of redundancies in speech makes it possible to use encoded descriptions of the speech that can be carried in a communication channel having a bandwidth that is less than the bandwidth of the speech. This is in distinct contrast to many well-known forms of converting speech into digital signals. Most of these methods require a bandwidth that is greater than the bandwidth of the speech.
Linear predictive coding (LPC) of speech begins conceptually with a model of the human speech-producing system. The model has a source of sound that is analogous to the vocal cords. That source is coupled acoustically to a stacked array of hollow cylindrical tubes that is analogous to the cavities of the throat and mouth in a human speaker. From the model, speech is characterized by four types of quantities. The first of these is a measure of whether speech is voiced or unvoiced. A voiced signal begins with an input from the vocal cords, while an unvoiced signal is produced by the action of the rest of the vocal tract on moving air alone. This produces the differences in sound between "s" which is unvoiced and its voiced equivalent "z", for example. A second characteristic of the sound is the pitch, the fundamental frequency produced by the vocal cords in making a voiced sound. A third characteristic is the energy. Finally, the effect of the throat and mouth on either voiced or unvoiced sound is summarized by obtaining some measure of the transfer function of the vocal tract. Such a measure might be the reflection coefficients of the structure, the poles of the transfer function of the structure, the logarithmic area ratios (LAR's) of the structure, or any of several other well-known functions of such resonances. In addition, various mathematical transforms of these functions may have utility for particular purposes. The functions are interrelated so that any one set can be determined from a knowledge of any other set.
The present invention is an improved method and means of determining quantities corresponding to reflection coefficients. The reflection coefficients are coefficients with a specific all pole filter structure known as a lattice filter. The reflection coefficients discussed here are the electrical analog of acoustical reflection coefficients. While the reflection coefficient is in general a complex quantity in which the imaginary portion is a measure of loss at a discontinuity, when determining the reflection coefficients that characterize a human vocal tract, the portions of reflection coefficients associated with loss are small enough compared to the real part that denotes lossless reflection that it is adequate to use a lossless model in which the reflection coefficients are real.
One way of obtaining the reflection coefficients that characterize a sample of a signal such as speech is to determine the characteristics of an inverse lattice filter that will reproduce the signal when excited by an impulse. Examples of several approaches of performing such an analysis are given in an article entitled "Stable and Efficient Lattice Methods for Linear Prediction" by Makhoul in "IEEE Transactions on Acoustics, Speech, and Signal Processing, Vol. ASSP-25, No. 5, October, 1977. This article points out that an analysis based on the all-pole lattice filter is stable without windowing but at a computational cost that is several times the cost of the auto-correlation and covariance methods of calculation. What Makhoul refers to as the computational costs of computer analysis are proportional to the equipment costs of realizing a circuit for comparable analysis on a semiconductor chip.
Summary of the Invention
It is an object of the present invention to provide an improved method and means of determining reflection coefficients of a signal.
It is a further object of the present invention to provide an improved method and means of determining parameters of an all-pole lattice filter that characterize an electrical signal.
Other objects will become apparent in the course of a detailed description of the invention.
An improved method and means of determining reflection coefficients that characterize an electrical signal obtains characteristics of an all-zero inverse lattice filter. The reflection coefficients are obtained by filtering the signal, sampling the filtered signal. obtaining the elements of a correlation matrix from the samples, initializing values of matrices of forward residual autocorrelations, backward residual autocorrelations, and cross correlation of residuals, combining matrix elements to obtain a first reflection coefficient, removing from the forward, backward and cross-correlation matrices the effect of the first reflection coefficient, calculating from the revised matrices a second reflection coefficient, and repeating the calculations to the desired order.
Brief Description of the Drawings
Fig. 1 is an overall block diagram of a circuit for the practice of the present invention.
Fig. 2 is an expanded circuit diagram of the address calculator of Fig. 1. Fig. 3 is an expanded circuit diagram of the update calculator of Fig. 1.
Fig. 4 is an expanded circuit diagram of the reflection-coefficient calculator of Fig. 1.
Fig. 5 is a flowchart of a method for practicing the invention.
Fig. 6 is a block diagram of a section of a typical lattice filter for the practice of the present invention.
Detailed Description of the Invention
Fig. 1 is an overall block diagram of a circuit for the practice of the invention. In Fig. 1 an electrical analog signal that is to be analyzed is applied at terminal 10. That signal will typically be an electrical analog of a voice signal, although it may be any electrical signal that exhibits redundancies analogous to those of speech. Examples of such other signals include video scans and seismic analysis records. Whatever the source, the signal at terminal 10 is applied to a filter 12, if necessary, to limit its bandwidth. If the bandwidth is already adequate, filter 12 may comprise a direct wire connection, or filter 12 may combine a bandpass filter with any of a number of systems of pre-emphasis that are commonly used in radio broadcasting. The output of filter 12, treated as described, is applied to analog-to-digital converter (ADC) 14. The digital output of ADC 14 will be separated into frames of a convenient length, of the order of tens of milliseconds. That function is here indicated as being performed by framer 16 as a means of insuring that the digital input to correlator 18 establishes correlation among samples in the same frame. The function of framer 16 could also be combined into correlator 18 or ADC 14. It should be noted that the ADC 14 is not necessary if the signal is already in digital form. The correlation parameters are obtained in the circuit of Fig. 1 because that circuit will be used to determine reflection coefficients.
The output of correlator 18 is applied, if necessary, to a normalizer 20 to normalize output values to a common level. The output of normalizer 20 is taken to random-access memory (RAM) 22 where it is stored in an address that is directed by a signal from multiplexer 23. Normalizer 20 also generates a signal indicating the completion of the operation of correlator 18 for one frame. That signal is taken as an input to sequencer 24. Signals from sequencer 24 are coupled out to control an address calculator 26, an update calculator 28, a reflection coefficient calculator 30, multiplexer 23, and multiplexer 32. Normalizer 20 determines the appropriate addresses in RAM 22 for storing the coefficients of a correlation matrix so that their serial readout will be in a desired order. The output of normalizer 20 is applied to multiplexer 23 to apply first an initial condition to RAM 22 that is determined from normalizer 20, and then accessed by address instructions from address calculator 26. The output of RAM 22 is applied through multiplexer 32 to RAM 34. Multiplexer 32 selects as an input to RAM 34 either the output of RAM 22 or the output of update calculator 28. Multiplexer 32 is controlled by a signal from sequencer 24. The location of storage elements in RAM 34 is controlled by a signal from address calculator 26. The output of RAM 34 is read into update calculator 26 and reflection coefficient calculator 30. update calculator 28 also receives as an input the calculated reflection coefficients from reflection coefficient calculator 30, which coefficients are taken as the output of the circuit of Fig. 1.
Fig. 2 is a circuit diagram of the address calculator 26 of Fig. 1. Address calculator 26 produces the addresses that direct the retrieval of matrix elements in RAMs 22 and the storage and retrieval of matrix elements in RAM 34. In Fig. 2 sequencer 24 supplies addresses to read-only memory (ROM) 40. ROM 40 is read out into multiplexers 42 and 44. The output of multiplexer 42 is taken to RAM 46 which is loaded under the control of signals from sequencer 24. The output of RAM 46 is taken directly as the input to RAM 34 of Fig. 1, and it is also applied to latches 48 and 50, as well as the add/subtract unit 42. The output of latch 48 is taken to multiplexer 23 of Fig. 1, and then to RAM 22. The output of latch 50 is multiplexed in multiplexer 44 with the output of ROM 40 under the control of a signal from sequencer 24. The output of multiplexer 44 is applied to add-subtract unit 52 where it is supplied as an input to be either added or subtracted from the output of RAM 46, depending upon a control signal from sequencer 24. The add/subtract unit 52 also has an output connected to sequencer 24 which indicates a zero result from the add/subtract operation.
Fig. 3 is a circuit diagram of the update calculator 28 of Fig. 1. In Fig. 3 data inputs are applied at terminals 60 and 62. The input at terminal 60 is the output of RAM 34 of Fig. 1. Terminal 60 is connected to a shifter 64, the output of which is taken to a RAM 66 and to a multiplexer 68. The output of RAM 66 is taken to multiplexer 70 which, in turn, is connected to multiplier 72. The output of multiplier 72 is taken as an input to multiplexer 68, as an input to register 74, and as an input to multiplexer 76. The output of register 74 is coupled through multiplexer 78 to supply a second input to multiplier 72. Terminal 62 is connected to reflection-coefficient calculator 30 of Fig. 1 to receive calculated reflection coefficients which are then coupled to register 80 of Fig. 3. The outputs of registers 74 and 80 are applied to multiplexer 78 to be selected under the control of a signal from sequencer 24 of Fig. 1. The output of register 80 of Fig. 3 is also taken as an input to multiplexer 70. Multiplexer 76 produces an output that is taken to registers 82 and 84, the outputs of which supply inputs to multiplexer 86. The output of multiplier 86 is taken as an input to multiplexer 70 and as an input to summer 88. The output of summer 88 is taken as input to multiplexer 76 as a data input to RAM 34 of Fig. 1 and as an input to OR gate 90 of Fig. 3. The output of OR gate 90 is taken to register 92 which supplies an output that is taken both as a second input to OR gate 90 and as an input to priority encoder 94. The output of priority encoder 94 controls register 96 that operates shifter 64. The combination of shifter 64, register 96, OR gate 90, register 92 and priority encoder 94 comprises a normalizer that normalizes the output of RAM 34 of Fig. 1. If it is not desired or is considered unnecessary to normalize coefficients, then terminal 60 can supply an input directly to RAM 66 and the elements just described could be removed from the circuit of Fig. 3. Fig. 4 is a circuit diagram of the reflection coefficient calculator 30 of Fig. 1. In Fig. 4 terminal 100 receives a signal from RAM 34 that is applied to a subtractor 102 and a multiply-by-two circuit 104. The output of subtractor 102 is connected to register 106, which produces one output that is connected as a input to subtractor 102, and another output that is taken to divider 108. The output of divider 108 may be taken directly as a reflection coefficient. In the alternative, it may be desirable to quantize that value as in quantizer 110. In either case, the signal at terminal 112 is a reflection coefficient of the equivalent filter of the original input signal. The reflection coefficients may be used directly as the coefficients in LPC or they may be transformed into a different function, as described earlier. If it is desired to convert reflection coefficients into LAR's, this can be combined readily with the quantization that is performed by quantizer 110. The choice of the particular functions is one of design. Consider now the elements of Fig. 1 in terms of their functions. Filter 12 has already been described as optional. If the circuit of Fig. 1 is to be applied to speech for use in a typical radio speech bandwidth, then the input signal may need to be subjected to bandpass filtering. It may also be desirable to combine some form of pre-emphasis with filter 12. Filter 12 may also be used to prevent aliasing when the filtered signal is applied to ADC 14. The sample rate produced by ADC 14 is a design parameter. According to the well-known sampling theorem, it will be necessary to sample the input at a rate that is at least twice the frequency of the highest component contained in the input. Thus, sampling rates for speech are typically between 6.4 kHz and 10 kHz. Each sample is then encoded into a number of bits that typically ranges between 8 and 15, with 12 as a typical number. It follows that a typical frame will have of the order of one hundred samples. Correlator 18 takes these samples and determines from them the elements of a covariance matrix. This matrix is symmetric. Element (i,k) of the covariance matrix is obtained by summing, from n = p to the total number of samples, the product of the (n - i)th and the (n - k)th samples for zero less than or equal to k, less than or equal to i , less than or equal to p. In this statement, p is the prediction order, a number that is typically between 8 and 12 for speech. The index k is kept less than or equal to i to avoid recalculating equal terms on both sides of the axis of the symmetric matrix.
Normalizer 20 of Fig. 1 performs a function that is here indicated separately but that might also be included in correlator 18. Normalizer 20 shifts the elements of the correlation matrix so that the maximum value of any element in the array is between one-half and one in magnitude. Normalizer 20 then truncates the values of the elements thus shifted to a number of bits equal to the word length of the system. The covariance matrix thus has a number of elements equal to (p + 1)2.
Since the covariance matrix is symmetric, it can be described completely by storing the elements of the diagonal and the elements below the diagonal, a total of (p + 1)(p + 2)/2 elements. These elements are stored in RAM 22 in a location that is controlled by address calculator 26. A convenient method of loading RAM 22 is to load the diagonal elements, beginning with the element of highest order and proceeding to the diagonal element of lowest order, and then repeat in sequence down paths in the matrix that are parallel to the main diagonal. It should be noted that the F and B matrices are also symmetric and may be stored in a similar fashion.
The operation of the circuits of Figs. 1, 2, 3 and 4 is explained further in the flowchart of Fig. 5 which shows the sequence of operations that are performed by those circuits. In Fig. 5, after a start 120, an operations block 122 directs the determination of correlation coefficients. Such a determination is well known. For discrete or sampled components, it is normally accomplished by a calculation such as that of the following equation:
N-1
∅(i,k) = Σ s(n-i)s(n-k) , O≤K≤i≤p n=p
where s(n) , O≤i ≤ N-1 are samples of the audio signal during a frame and p is the order of the filter.
The next step in the flowchart of Fig. 5 is to initialize the matrices of F, B and C as indicated in operations block 126. The quantities F, B and C are intermediate quantities used in the determination of LPC coefficients. Their initial values are determined as follows:
FO(i,k) = ∅(i,k) BO(i,k) = ∅(i+1,k+1) CO(i,k) = ∅(i, k+1) for O ≤ i, k ≤ p-1
Operations box 1 28 next directs that the value of j equal 1 . Operations box 1 30 then calls for the determination of the quantity kj , the jth reflection coefficient. This is determined as follows:
A variety of techniques are available for determining the reflection coefficient and the above technique is intended as an illustration of one embodiment of the present invention.
After the value of kj is determined, it is quantized according to the instruction from operations block 132 at a quantization level that is determined by the number of bits reserved to broadcast the particular LPC coefficient in question. Decision block 134 next tests to see whether j = p. If it does, all of the LPC coefficients have been calculated, and exit is to operations block 136 to end the calculation. If j does not equal p, the calculated value of kj is used in operations block 138 to update values of F, B and C according to the following relations:
Fj(i,k)=Fj-1(i,k)+kj[Cj-1(i,k)+Cj-1(k,i)]+kj 2Bj-1(i,k)
Bj (i,k)=Bj-1(i+l,k+1)+kj[Cj-1(i+l,k+1)+Cj-1(k+1,i+1)]+
Kj 2Fj-1(i+1,k+1)
Cj(i,k)=Cj-1(i,k+1)+kj[Bj-1(i,k+1)+Fj-1(i,k+1)]+kj 2Cj-1(k+1,i)
O≤i,k≤(p-j-1)
The value of j is increased by 1 in operations block 140, and control returns to operations block 130 to continue the calculations.
The preceding description shows a circuit for determining LPC coefficients and a flowchart of the steps performed by that circuit. The flowchart of Fig. 5 is also useful in directing the steps of a program for determining LPC coefficients by the same method using a computer. This has been done, and a computer program for determining LPC coefficients using this method is included as an appendix to this application. The computer program has the particular advantage that, using the method of the present invention, calculations may be performed efficiently in integer arithmetic. This speeds calculations and makes it possible to determine LPC coefficients in real time using the present invention. Fig. 6 is a block diagram of a lattice filter that provides a further explanation of the process by which the circuit of Fig. 1 obtains reflection coefficients. In Fig. 6 a terminal 150 receives as an input the sampled signal. This signal is applied to an upper leg 152 which will calculate a forward residual, and it is applied to a lower leg 154 which will apply the signal to the first of a sequence of delay elements 156 to calculate a backward residual. Fig. 6 comprises a cascade of elements, each of which applies the forward signal to a multiplier 158 and a summer 165. The backward signals are applied to a multiplier 162 and a summer 164. Both multipliers 158 and 162 receive as additional inputs the current reflection coefficient. Thus, the current forward residual is multiplied by the current reflection coefficient and added to the current backward residual in summer 164 to generate as an output the next backward residual. The current backward residual is multiplied in multiplier 162 with the current reflection coefficient and added to the current forward residual in summer 160 to generate the next forward residual. The process just described continues through a number of sections of the lattice filter that is determined by the designer as a number adequate to characterize the particular signal in question. This is typically a number of stages equal to 8, 10 or 12. The last such stage is shown here as receiving a current forward residual signal on terminal 166 and a current backward residual signal on terminal 168. The current forward residual signal is applied to a multiplier 170 and a summer 172, while the current backward residual is delayed in delay element 174, and the delayed signal is applied to multiplier 176 and summer 178. Both multipliers 170 and 176 receive as additional inputs the current reflection coefficient. If the lattice filter has been designed with an adequate number of sections to approximate the input signal sufficiently well, then the output of the final forward residual signal at terminal 180 will be close to zero and so will the final backward residual signal at terminal 182. The result of applying a signal to this circuit of Fig. 6 is the production in forward line 152 of a sequence of elements of a forward residual vector and to produce in line 156 the elements of a backward residual vector. Individual elements are combined to produce an autocorrelation matrix of the forward residual elements, and autocorrelation matrix of the backward residual elements, and a cross-correlation matrix between forward and backward residual elements. These matrices are used as described earlier to calculate the reflection coefficients.
Lattice methods of determining coefficients for linear predictive coding have been used in the past. However, circuits and programs used to determine the lattice coefficients have used intermediate variables that varied in magnitude over a wide range. This required a wide range of quantize values to characterize the intermediate variables, and thus took more computational time. The circuit arrangement and method of the present invention uses only variable and intermediate variables which are bounded in magnitude by unity. This permits operations and calculations to be performed in a fixed-point fractional implementation. In addition if the input signal is windowed so that it is stationary in a statistical sense, then it can be shown that the number of computations necessary. to determine lattice coeficients is reduced still further. The method and means of the present invention has been used with a frame length of approximately 15 milloseconds to determine 12 lattice coefficients in real time. By completing the calculations for the data of one frame before the end of the next succeeding frame. What is claimed is:
APPENDIX FORTRAN PROGRAM FOR COEFFICIENTS
C
C FLAT
C MODIFIED COVAR LATTICE LPC ANALYSIS
C
∅ ∅ ∅1 LOGICAL*1 IDATE( 9 ) ,ITIME(8)
∅ ∅ ∅2 DIMENSION AR( 2 ∅, 3 ) , ISAMP( 512 ) ,F( 256 ) ,AK( 19 ) ,BA( 19 ) 1,C(2 ∅,2 ∅) ,S(2 ∅),BR( 19,19) ,FR( 19,19) ,CR( 19, 19 ) ,DR(19,19)
∅ ∅ ∅3 DATA APAST/ ∅./BA/19* ∅./ISAMP/512* ∅/AK/19* ∅./S/2 ∅*1 1/IFLAG/ ∅/NPT/1/IPO/1/IPI/1/IPP/1/JY/Υ'/
∅ ∅ ∅4 CALL INITT(96 ∅)
∅ ∅ ∅5 CALL TERM (3,1 ∅23)
∅ ∅ ∅6 CALL DATE(IDATE)
∅ ∅ ∅7 CALL TIME(ITIME)
∅ ∅ ∅8 TYPE 5,ITIME,IDATE
∅ ∅ ∅9 5 FORMAT (' COVAR LATTICE LPC ANALYSIS ' ,
15X,8A1,5X,9A1,/)
∅ ∅1 ∅ TYPE 1 ∅
∅ ∅11 1 ∅ FORMAT (' ENTER INPUT FILE NAME',/)
∅ ∅12 CALL ASSIGN (1,,-1)
∅ ∅13 TYPE 2 ∅
∅ ∅14 2 ∅ FORMAT (' ENTER OUTPUT FILE NAME',/)
∅ ∅15 CALL ASSIGN (2,,-1)
∅ ∅16 DEFINE FILE 1 ( ∅,256 ,U, IPI)
∅ ∅17 DEFINE FILE 2( ∅, 128,U,IPO)
∅ ∅18 TYPE 3 ∅
∅ ∅19 3 ∅ FORMAT (' ENTER PREEMPHASIS COEFFICIENT')
∅ ∅2 ∅ ACCEPT 4 ∅,APRE
∅ ∅21 4 ∅ FORMAT (F1 ∅. ∅)
∅ ∅22 TYPE 5 ∅
∅ ∅23 5 ∅ FORMAT (' ENTER NUMBER OF SAMPLES PER FRAME')
∅ ∅24 ACCEPT 4 ∅,AA
∅ ∅25 LFRAME=AA
∅ ∅26 TYPE 6 ∅
∅ ∅27 6 ∅ FORMAT (' ENTER ORDER OF PREDICTOR')
∅ ∅28 ACCEPT 4 ∅,AA
∅ ∅29 NP=AA
∅ ∅3 ∅ 69 TYPE 7 ∅
∅ ∅31 7 ∅ FORMAT (' ENTER ANALYSIS LENGTH')
∅ ∅32 ACCEPT 4 ∅,AA
∅ ∅33 LEN=AA
∅ ∅34 TYPE 8 ∅
∅ ∅35 8 ∅ FORMAT (' ENTER MAX # OF FRAMES')
∅ ∅36 ACCEPT 4 ∅,AA ∅∅37 NMAX=AA
∅ ∅38 TYPE 9 ∅
∅ ∅39 9 ∅ FORMAT ( ' ENTER GAIN' )
∅ ∅4 ∅ ACCEPT 4 ∅,GAIN
∅ ∅41 TYPE 92
∅ ∅42 92 FORMAT (' ENTER N FOR BINOMIAL SMOOTHINGN' )
∅ ∅43 ACCEPT 4 ∅,AA
∅ ∅44 NB=AA
∅ ∅45 IF (NB.EQ. ∅) GO TO 97
∅ ∅ 47 DO 93 19=1 ,NP
∅ ∅48 93 S(I9+1)=S(I9)*(NB-I9+1)/(NB+I9)
∅ ∅49 TYPE 94,S
∅ ∅5 ∅ 94 FORMAT (5F15.1 ∅)
∅ ∅51 97 TYPE 95
∅ ∅52 95 FORMAT (' DO YOU WANT QUANT IN RECURSION? Y Y/N' )
∅ ∅53 ACCEPT 63,IIY
∅ ∅54 63 FORMAT (A1)
∅ ∅55 IF (IIY.NE.JY) GO TO 99
∅ ∅57 TYPE 96
∅ ∅58 96 FORMAT (' ENTER QUANT FILE',/)
∅ ∅59 CALL ASSIGN (3,,-1)
∅ ∅6 ∅ IFLAG=1
∅ ∅61 READ (3,4 ∅)XXX
∅ ∅62 KK=XXX+1. ∅1
∅ ∅63 NP=MIN0(NP,KK-1) .
∅ ∅64 READ (3,4 ∅)((AR(K1,K2) ,K2=1,3) ,K1=1 ,KK)
∅ ∅65 99 READ (1 ,IPI)(ISAMP(IV),IV=257,512)
∅ ∅66 NPT=(LFRAME-LEN)/2+257
∅ ∅67 1 ∅ ∅ IF (IPO.GT.NMAX) GO TO 1 ∅ ∅ ∅
∅ ∅69 IF (NPT.LE.257) GO TO 115
∅ ∅71 DO 11 ∅ ILO=1,256
∅ ∅72 11 ∅ ISAMP(ILO)=ISAMP(ILO+256)
∅ ∅73 READ (1 'IPI,END=1 ∅ ∅ ∅) (ISAMP(IV) ,IV=257,512)
∅ ∅74 NPT=NPT-256
∅ ∅75 115 APAST=ISAMP(NPT-1)
∅ ∅76 12 ∅ NLEN=LEN
∅ ∅77 125 DO 13 ∅ IL2=1,NLEN
∅ ∅78 BAST=ISAMP(NPT+IL2-1)
∅ ∅79 F( IL2 )=( BAST-APRE*APAST) /32768.
∅ ∅80 13 ∅ APAST=BAST
C COMPUTE COVARIANCE MATRIX (EQ. 41)
∅ ∅81 DO 135 I1=1,NP+1
∅ ∅82 DO 135 12=1,I1
∅ ∅83 C(I1,I2)= ∅
∅ ∅84 DO 135 I3=NP+1,NLEN
∅ ∅85 135 C(I1,I2)=C(I1,I2)+F(I3-I2+1)*F(I3-I1+1)
∅ ∅86 DO 138 11=1.NP+1
∅ ∅87 DO 138 12=1,I1 ∅∅88 C(I1,I2)=C(I1,I2)*S(I1-I2+1)
∅ ∅89 138 C(I2,I1)=C(I1,I2)
∅ ∅9 ∅ R ∅=( C(1,1)+C(NP+1 ,NP+1 ) ) *GAIN/( 2.* (NLEN-NP) ) C INIT CR,BR AND FR
∅ ∅91 DO 14 ∅ 11=1 ,NP
∅ ∅92 DO 14 ∅ I2=1,NP
∅ ∅93 FR I1,I2)=C(I1,I2)
∅ ∅94 BR(I1,I2)=C(I1+1,I2+1)
∅ ∅95 14 ∅ CR(I1,I2)=C(I1,I2+1)
∅ ∅96 II=1
C LOOP STARTS HERE
∅ ∅97 15 ∅ IND=NP-II+1
C COMPUTE AK(II)
∅ ∅98 DENOM=FR(1,1)+BR(1,1)+FR(IND,IND)+BR(IND,IND)
∅ ∅99 IF (DENOM.EQ. ∅.) GO TO 165
∅1 ∅1 AK(II)=-2.*(CR(1,1)+CR(IND,IND))/DENOM
∅1 ∅2 IF (ABS(AK(II)) .GE.1.) TYPE 155
∅1 ∅4 155 FORMAT (' ERROR IN COMPUTATION')
∅1 ∅5 165 IF (IFLAG.EQ.1) AK( II )=QU(AK( II) , II,AR) C SEE IF DONE
∅1 ∅7 IF (II.EQ.NP) GO TO 2 ∅ ∅
∅1 ∅9 Z=AK(II)
∅11 ∅ ZS=Z*Z
C MODIFY FR,BR AND CR ARRAYS
∅111 DO 17 ∅ I1=1,NP-II
∅112 DO 17 ∅ I2=1,NP-II
∅113 FR(I1,I2)=FR(I1,I2)+ 1Z*(CR(I1,I2)+CR(I2,I1))+ZS*BR(I1,I2)
∅114 BR(I1,I2)=BR(I1+1,I2+1)+ 1Z*(CR(I1+1,I2+1)+CR(I2+1,I1+1))+ZS*FR(I1+1,I2+1)
∅115 17 ∅ DR(I1,I2)=CR(I1,I2+1)+
1Z*(BR(I1,I2+1)+FR(I1,I2+1))+ZS*CR(I2+1,I1)
∅116 DO 18 ∅ I1=1,NP-II
∅117 DO 18 ∅ I2=1,NP-II
∅118 18 ∅ CR(I1,I2)=DR(I1,I2)
∅119 II=II+1
∅12 ∅ GO TO 15 ∅
∅121 2 ∅ ∅ PO=FR(1 , 1 )+2*AK(NP) *CR( 1 , 1 ) +AK(NP) *AK(NP) *BR( 1 ,1)
∅122 PO=PO*GAIN/(NLEN-NP)
∅123 IF( IFLAG.EQ.1) PO=QU(PO, ∅,AR)
∅125 WRITE (2'IPO)NP,LFRAME,PO,APRE,BA,AK,R ∅
∅126 ISW=IPEEK(" 17757 ∅)
∅127 IF(ISW.AND." 1 ∅ ∅ ∅ ∅ ∅) TYPE 16 ∅, IPI-1 , IPO-1
∅129 16 ∅ FORMAT ( ' ',218)
∅13 ∅ -NPT=NPT+LFRAME
∅131 GO TO 1 ∅ ∅
∅132 10 ∅ ∅ WRITE (2'IPO)(-1)
∅133 CALL CLOSEO) ∅134 CALL CLOSE (2)
∅135 CALL CLOSE(3)
∅136 CALL EXIT
∅137 END

Claims

Claims
1. A method of processing an analog electrical signal to obtain lattice coefficient of an inverse lattice filter that characterizes the analog electrical signal, the method comprising the steps of: a) converting the analog electrical signal to a digital signal; b) selecting samples of the digital signal for a time period equal to a predetermined frame length; c) determining correlation coefficients of a selected samples; d) initializing values of F, B, and C matrices; e) determining a first lattice coefficient from initial values of F, B, and C matrices; f) adjusting values of F, B, and C matrices from the first lattice coefficient; and g) determining a second lattice coefficient from the adjusted F, B, and C matrices and repeating the steps of adjusting and determining until a predetermined number of lattice coefficients has been obtained.
2. The method of claim 1 comprising in addition the step of filtering the analog signal before it is converted to a digital signal .
3. The method of claim 1 wherein the step of determining correlation coefficients comprises the step of obtaining the sum of products of sampled signals according to the relationship
N-1 ∅(i,k) = ∑ s(n-i)s(n-k) , O≤K≤i≤p n=p
where s(n) , O<n < N-1 are samples of the audio signal during a frame and p is the order of the filter.
4. The method of claim 1 wherein the F matrice is defined according to the relation and
Fj(i,k)=Fj-1(i,k)+kj[Cj-1(i,k)+Cj-1(k,i)]+kj 2Bj-1(i,k) O≤i,k≤(p-j-l)
5. The method of claim 1 wherein the B matrice is determined according to
Bj (i,k)=Bj-1(i+1,k+1)+kj[Cj-1(i+1,k+1)+Cj-1(k+l,i+1)]+
kj 2Fj-1(i+1,k+1) O≤i,k≤(p-j-1)
6. The method of claim 1 wherein the C matrice is determined according to,
Cj (i,k)=Cj-1(i,k+1)+kj[Bj-1(i,k+1)+Fj-1(i,k+1)]+kj 2Cj-1(k+1,i) O≤i,k≤(p-j-1)
7. The method of claim 1 wherein converting the analog electrical signal is not necessary if the signal is already in digital form.
8. An apparatus for processing an analog electrical signal to obtain reflection coefficients of an inverse lattice filter that characterizes the analog electrically signal, the apparatus comprising: a) a filter receiving the analog electrical signal and producing a filtered analog electrical signal; b) analog-to-digital converter connected to the filter to produce digital signal , c) a framer to select time interval s of a predetermined amount; d) corel ator connected to the framer and to the analog-to-digital converter to produce corelation coefficients of the digital signals during each framed interval ; e) a normalizer connected to the corelator to normalize the corelation coefficients of the correlator, f) a first RAM connected to the normalizer to receive and store normalized corelation coefficients, g) a sequencer that generates a plurality of control signals, h) an address calculator connected to the sequencer and to the first RAM to control addresses at which normalized corelation coefficients are stored and the first RAM; i) a multiplexer connected to the first RAM and receiving as an input data from the first RAM j) a second RAM connected to the multiplexer and to the address calculator, the second RAM receiving input from the multiplexer and storing that input in a location controlled by the address calculator; k) a reflection coefficient calculator connected to the second RAM and to the sequencer, the reflection coefficient calculator receiving data from the second RAM and caculating reflection coefficients under control from the sequencer; and 1) an update circuit connected to the reflection coefficient calculator the second RAM and the sequencer to calculate current coefficients under the control of the sequencer and producing an output that is taken as an input to the multiplier.
EP84900575A 1983-01-03 1983-12-22 Improved method and means of determining coefficients for linear predictive coding Expired EP0135512B1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
AT84900575T ATE36206T1 (en) 1983-01-03 1983-12-22 IMPROVED METHOD AND MEANS FOR DETERMINING COEFFICIENTS FOR LINEAR PREDICTIVE CODING.

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US45492983A 1983-01-03 1983-01-03
US454929 1983-01-03

Publications (3)

Publication Number Publication Date
EP0135512A1 EP0135512A1 (en) 1985-04-03
EP0135512A4 true EP0135512A4 (en) 1985-06-10
EP0135512B1 EP0135512B1 (en) 1988-08-03

Family

ID=23806651

Family Applications (1)

Application Number Title Priority Date Filing Date
EP84900575A Expired EP0135512B1 (en) 1983-01-03 1983-12-22 Improved method and means of determining coefficients for linear predictive coding

Country Status (9)

Country Link
EP (1) EP0135512B1 (en)
JP (1) JPS60500274A (en)
AT (1) ATE36206T1 (en)
AU (1) AU566370B2 (en)
DE (1) DE3377599D1 (en)
DK (1) DK417184A (en)
FI (1) FI843448A (en)
NO (1) NO843498L (en)
WO (1) WO1984002814A1 (en)

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB1318985A (en) * 1970-02-07 1973-05-31 Nippon Telegraph & Telephone Audio response apparatus
GB2026289B (en) * 1978-04-12 1982-04-21 Secr Defence Self-adaptive linear prediction filters
JPS597120B2 (en) * 1978-11-24 1984-02-16 日本電気株式会社 speech analysis device
JPS55164700U (en) * 1979-05-14 1980-11-26
JPS5853358B2 (en) * 1980-03-31 1983-11-29 株式会社東芝 speech analysis device
US4378469A (en) * 1981-05-26 1983-03-29 Motorola Inc. Human voice analyzing apparatus

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
COMPCON 79 (19TH IEEE COMPUTER SOCIETY INTERNATIONAL CONFERENCE), Washington, D.C., 4th-7th September 1979, pages 203-206, IEEE, New York, US; A.J. GOLDBERG et al.: "Microprocessor implementation of a linear predictive coder" *
ICASSP 81 (IEEE INTERNATIONAL CONFERENCE ON ACOUSTICS, SPEECH AND SIGNAL PROCESSING), Atlanta, 30th March - 1st April 1981, pages 355-353, IEEE, New York, US; F. SUGIYAMA et al.: "An LPC vocoder for efficient implementation" *

Also Published As

Publication number Publication date
JPS60500274A (en) 1985-02-28
NO843498L (en) 1984-09-03
DK417184A (en) 1984-10-04
DE3377599D1 (en) 1988-09-08
DK417184D0 (en) 1984-08-31
AU566370B2 (en) 1987-10-15
ATE36206T1 (en) 1988-08-15
EP0135512B1 (en) 1988-08-03
AU2439184A (en) 1984-08-02
FI843448A0 (en) 1984-09-03
WO1984002814A1 (en) 1984-07-19
EP0135512A1 (en) 1985-04-03
FI843448A (en) 1984-09-03

Similar Documents

Publication Publication Date Title
US4544919A (en) Method and means of determining coefficients for linear predictive coding
US4301329A (en) Speech analysis and synthesis apparatus
JP3508146B2 (en) Digital signal encoding / decoding device, digital signal encoding device, and digital signal decoding device
US6484140B2 (en) Apparatus and method for encoding a signal as well as apparatus and method for decoding signal
JP2904472B2 (en) Method, data processing system and apparatus for efficiently compressing digital audio signals
EP0910067B1 (en) Audio signal coding and decoding methods and audio signal coder and decoder
JPH0697837A (en) Digital signal decoding device
EP0424121A2 (en) Speech coding system
EP0673014A2 (en) Acoustic signal transform coding method and decoding method
EP0700032A2 (en) Methods and apparatus for reconstructing non-quantized adaptively transformed voice signals
US4922539A (en) Method of encoding speech signals involving the extraction of speech formant candidates in real time
US5926785A (en) Speech encoding method and apparatus including a codebook storing a plurality of code vectors for encoding a speech signal
US4809331A (en) Apparatus and methods for speech analysis
US5511094A (en) Signal processor for a sub-band coding system
JP3248668B2 (en) Digital filter and acoustic encoding / decoding device
Chen et al. A fixed-point 16 kb/s LD-CELP algorithm
AU566370B2 (en) Improved method and means of determining coefficients for linear predictive coding
US6006177A (en) Apparatus for transmitting synthesized speech with high quality at a low bit rate
Schafer et al. Parametric representations of speech
JP3082625B2 (en) Audio signal processing circuit
JP4645869B2 (en) DIGITAL SIGNAL PROCESSING METHOD, LEARNING METHOD, DEVICE THEREOF, AND PROGRAM STORAGE MEDIUM
Buzo et al. A two-step speech compression system with vector quantizing
JP2897940B2 (en) Speech linear prediction parameter coding method
JP3010655B2 (en) Compression encoding apparatus and method, and decoding apparatus and method
GB2199215A (en) A stochastic coder

Legal Events

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

Free format text: ORIGINAL CODE: 0009012

17P Request for examination filed

Effective date: 19840830

AK Designated contracting states

Designated state(s): AT CH DE FR GB LI LU NL SE

17Q First examination report despatched

Effective date: 19860704

D17Q First examination report despatched (deleted)
GRAA (expected) grant

Free format text: ORIGINAL CODE: 0009210

AK Designated contracting states

Kind code of ref document: B1

Designated state(s): AT CH DE FR GB LI LU NL SE

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

Ref country code: NL

Effective date: 19880803

Ref country code: FR

Free format text: THE PATENT HAS BEEN ANNULLED BY A DECISION OF A NATIONAL AUTHORITY

Effective date: 19880803

Ref country code: AT

Effective date: 19880803

REF Corresponds to:

Ref document number: 36206

Country of ref document: AT

Date of ref document: 19880815

Kind code of ref document: T

REF Corresponds to:

Ref document number: 3377599

Country of ref document: DE

Date of ref document: 19880908

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

Ref country code: GB

Effective date: 19881222

EN Fr: translation not filed
PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: LU

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

Effective date: 19881231

Ref country code: LI

Effective date: 19881231

Ref country code: CH

Effective date: 19881231

NLV1 Nl: lapsed or annulled due to failure to fulfill the requirements of art. 29p and 29m of the patents act
PLBE No opposition filed within time limit

Free format text: ORIGINAL CODE: 0009261

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

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

26N No opposition filed
GBPC Gb: european patent ceased through non-payment of renewal fee
REG Reference to a national code

Ref country code: CH

Ref legal event code: PL

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

Ref country code: DE

Effective date: 19890901

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

Ref country code: SE

Effective date: 19891223

EUG Se: european patent has lapsed

Ref document number: 84900575.6

Effective date: 19900104