EP0786872A2 - Viterbi decoder with reduced metric computation - Google Patents

Viterbi decoder with reduced metric computation Download PDF

Info

Publication number
EP0786872A2
EP0786872A2 EP96309397A EP96309397A EP0786872A2 EP 0786872 A2 EP0786872 A2 EP 0786872A2 EP 96309397 A EP96309397 A EP 96309397A EP 96309397 A EP96309397 A EP 96309397A EP 0786872 A2 EP0786872 A2 EP 0786872A2
Authority
EP
European Patent Office
Prior art keywords
components
circuit
cross correlation
autocorrelation
terms
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.)
Withdrawn
Application number
EP96309397A
Other languages
German (de)
French (fr)
Inventor
R.A. Cesari
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.)
Nokia of America Corp
Original Assignee
Lucent Technologies 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 Lucent Technologies Inc filed Critical Lucent Technologies Inc
Publication of EP0786872A2 publication Critical patent/EP0786872A2/en
Withdrawn legal-status Critical Current

Links

Images

Classifications

    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03MCODING; DECODING; CODE CONVERSION IN GENERAL
    • H03M13/00Coding, decoding or code conversion, for error detection or error correction; Coding theory basic assumptions; Coding bounds; Error probability evaluation methods; Channel models; Simulation or testing of codes
    • H03M13/63Joint error correction and other techniques
    • H03M13/6343Error control coding in combination with techniques for partial response channels, e.g. recording
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03MCODING; DECODING; CODE CONVERSION IN GENERAL
    • H03M13/00Coding, decoding or code conversion, for error detection or error correction; Coding theory basic assumptions; Coding bounds; Error probability evaluation methods; Channel models; Simulation or testing of codes
    • H03M13/37Decoding methods or techniques, not specific to the particular type of coding provided for in groups H03M13/03 - H03M13/35
    • H03M13/39Sequence estimation, i.e. using statistical methods for the reconstruction of the original codes
    • H03M13/41Sequence estimation, i.e. using statistical methods for the reconstruction of the original codes using the Viterbi algorithm or Viterbi processors
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03MCODING; DECODING; CODE CONVERSION IN GENERAL
    • H03M13/00Coding, decoding or code conversion, for error detection or error correction; Coding theory basic assumptions; Coding bounds; Error probability evaluation methods; Channel models; Simulation or testing of codes
    • H03M13/37Decoding methods or techniques, not specific to the particular type of coding provided for in groups H03M13/03 - H03M13/35
    • H03M13/39Sequence estimation, i.e. using statistical methods for the reconstruction of the original codes
    • H03M13/41Sequence estimation, i.e. using statistical methods for the reconstruction of the original codes using the Viterbi algorithm or Viterbi processors
    • H03M13/4161Sequence estimation, i.e. using statistical methods for the reconstruction of the original codes using the Viterbi algorithm or Viterbi processors implementing path management
    • H03M13/4169Sequence estimation, i.e. using statistical methods for the reconstruction of the original codes using the Viterbi algorithm or Viterbi processors implementing path management using traceback
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03MCODING; DECODING; CODE CONVERSION IN GENERAL
    • H03M13/00Coding, decoding or code conversion, for error detection or error correction; Coding theory basic assumptions; Coding bounds; Error probability evaluation methods; Channel models; Simulation or testing of codes
    • H03M13/65Purpose and implementation aspects
    • H03M13/6502Reduction of hardware complexity or efficient processing
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/03Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
    • H04L25/03006Arrangements for removing intersymbol interference
    • H04L25/03178Arrangements involving sequence estimation techniques
    • H04L25/03184Details concerning the metric
    • H04L25/03197Details concerning the metric methods of calculation involving metrics

Definitions

  • This invention relates generally to viterbi decoders such as would be used in digital receivers, and particularly to a viterbi decoder with reduced metric computations.
  • Digital systems transmit information either as zeros or ones. While accurate transmission of information is generally not at issue, accurate detection of information received depends on being able to detect transmitted ones as ones upon being received, and being able to detect transmitted zeros as zero upon being received. Accurate detection is not a simple matter due to many factors. These factors include low received signal power at a receiver, movement of the transmitter relative to the receiver which introduces changing channel characteristics. Multi-path interference, or inter-symbol interference, occurs when reflections from various objects, including atmospheric layers, cause the presence of multiple transmission paths from the transmitter to the receiver, with varying transit times for a transmitted signal, causing several images of a transmitted signal to arrive at the receiver at different instants in time.
  • a viterbi decoder is a maximum likelihood decoder that provides forward error correction. Viterbi decoders are used to decode a sequence of encoded symbols, such as a bit stream.
  • the bit stream can represent encoded information in a telecommunications system. Such encoded information can be transmitted through various media with each bit (or set of bits) representing a symbol instant.
  • a viterbi decoder works back through a sequence of possible received bit sequences at each symbol instant to determine which one bit sequence was most likely to have been transmitted. The possible transition from a bit at one symbol instant, or state, to a bit at a next or a subsequent, symbol instant or state, is limited.
  • Each possible transition from one state to a next state can he shown graphically and is defined as a branch. All possible branches constitute a trellis. A sequence of interconnected branches is defined as a path. Each state can transition only to a limited number (for example two or four) of next states upon receipt of the next bit (or set or bits) in the bit branch. Thus, some paths survive and other paths do not survive during the decoding process. By eliminating those transitions that are not permissible, computation efficiency can be achieved in determining the most likely path to survive.
  • the viterbi decoder typically defines and calculates a branch metric associated with each branch and employs this branch metric to determine which paths survive and which paths do not survive.
  • a branch metric is calculated at each symbol instant for each possible branch.
  • Each path has an associated metric, accumulated cost, that is updated at each symbol instant.
  • the accumulated cost for the next state is calculated as the extremum of the sum of the branch metric for the possible transmission and the path accumulated cost at the respective previous, or originating, state.
  • Calculation at each symbol instant of the branch metric involves an add, compare and select operation for each state. This can result in computation-intensive operation.
  • a method and apparatus for generating branch metrics that are combinations of autocorrelation and cross correlation terms are constructed by first calculating and storing components of the autocorrelation and cross correlation terms at a symbol instant. Once calculated and stored, predetermined ones of the autocorrelation components and the cross correlation components are selected. The selected autocorrelation components and cross correlation components, or their inverse as predetermined, are combined to produce a branch metric. Other predetermined combinations of the stored components of autocorrelation and cross correlation terms, or their inverse, are combined to produce other branch metrics at the same symbol instant. All branch metrics associated with the symbol instant can be calculated in this manner. The components of the autocorrelation terms and cross correlation terms associated with the next symbol instant are calculated and stored, and the process of generating the branch metrics repeated.
  • Transceiver 10 shown in the block diagram of Figure 1 is an illustrative embodiment of the invention.
  • Transceiver 10 is a communication system or portion of a communication system such as a mobile digital cellular telephone.
  • Transceiver 10 is comprised of transmitter 12 and receiver 14.
  • a microphone 121 for converting a voice signal into an electrical signal is coupled to a vocoder 122.
  • the vocoder is coupled to the channel encoder 123, which includes reordering and partitioning 1231, interleaving and formatting 1232, and symbol converter 1234.
  • the output of the channel encoder is coupled to a transmitter filter 124.
  • the filter output is coupled to a continuous balanced modulator 125, which in turn is coupled to another filter 126.
  • Filter 126 is coupled to an up converter and power amplifier 127 which is turn is coupled to an antenna 128.
  • antenna 141 is coupled to a low noise amplifier and down conversion 142, which is coupled to filtering 143.
  • Filtering 143 is coupled to automatic gain control 144 which conditions the signal for input to the Viterbi decoder 145.
  • the Viterbi decoder is coupled to a channel decoder 146 which includes deinterleaving, deformatting, deordering and departitioning 1461.
  • the channel decoder is coupled to a vocoder 147.
  • the vocoder is coupled to a speaker 148.
  • Figure 2 shows a trellis diagram useful in understanding the branch metrics used in a Viterbi decoder.
  • a Viterbi decoder reconstructs a transmitted bit stream from a received bit stream. Perfect digital signal symbols in a sequence are transmitted. However, the received signal is distorted due to noise and interference.
  • the mobile digital cellular telephone moves through its environment, the characteristics of the channel over which the signals are transmitted and received change.
  • states shown vertically
  • C-l -1 where C is the constraint length.
  • These 4 C-l states are referred to as individual states.
  • the portion of a trellis diagram shown in Figure 2 is constructed one symbol instant at a time, as is known in the art.
  • each individual state of the illustrative embodiment shown in Figure 2 can only transition to four possible individual states at the next symbol instant. Concomitantly, each individual state at a symbol instant has only four possible previous individual states from which a transition to it can originate. Other decoders could have more or less than four transitions from a previous state and to a next state.
  • a branch metric is calculated at each symbol instant for each possible transition from an individual state at one symbol instant, known as a source state, to an individual state at the next symbol instant.
  • Various methods for calculating branch metrics are known in art.
  • the branch metric of all branches transitioning to a given next individual state are calculated, then added to the accumulated cost of the respective originating individual state. This results in four potential accumulated cost sums.
  • a comparison is made of the four potential accumulated cost sums.
  • An extremum branch metric such as the least of the four sums, is selected as the next state accumulated cost.
  • a transition corresponding to the least sum is the most likely transition from the four possible originating individual states to the given next individual state.
  • the surviving branch data of the most likely transition is stored forming the trellis of Figure 2 as the surviving branch to the given next individual state as is known in the art.
  • the extremum branch metric is the accumulated cost of the next individual state. This process is repeated for each of the next individual states and repeated also for each symbol instant in the sequence.
  • Accumulated cost and the source state for the selected branch can be stored in an array.
  • a decoded symbol at each symbol instant can be obtained by a traceback operation. Starting with the individual state in the last symbol instant having the minimum associated accumulated cost, a single trace back through the trellis will provide a decoded symbol at each symbol instant. At each symbol instant of this traceback operation, the source state is translated into a decoded symbol for that symbol instant, as is known in the art.
  • a viterbi decoder is used with a partial log likelihood estimator to determine the most likely symbol sequence taken from a set of possible sequences over an interval of time. Minimizing the integral squared error between the received sequence and the expected sequence results in the most likely sequence to have been transmitted through a channel.
  • the predicted squared error or log-likelihood between the received signal and an expected sequence (or complex functions is the auto correlation of the predicted error functions, multiplied by a constant. The constant may be ignored when evaluating relative likelihoods.
  • Equation (2) represents the integral of the signal energy minus the difference between the cross correlation of the expected sequence.
  • the three terms of equation (2) represent energy, cross-correlation, and autocorrelation respectively.
  • log likelihood (m) S(i-l) m . y(n) * . h n-l (k+l) + S(i) m . y(n)* . h n (k) +S(i+l) m . y(n)* . h n+l (k-l) + [S(i-l) . S(i)] m . h n-l (k+l) . h n (k+l) + [S(i-l)* . S(i+l)] m . h n-l (k)* . h (n+l) (k) + [S(i)* . S(i+l)] m . h n-l (k)* . h (n+l) (k) + [S(i)* . S(i+l)] m . h n (k-l)* . h (n+l)
  • Odd indexed symbols are expected to be one of four odd phases on the unit circle [2 ⁇ n/8; n odd]
  • the MLSE operation remains identical for both even and odd indexed symbols, and eliminates the need for performing a scalar operation for each state in a symbol interval.
  • the - ⁇ /4 phase rotation requires calculating the sum and the difference of the real and imaginary value of the expected values and multiplying then by ⁇ 2, and re-ordering these scaled sums and differences in four possible arrangements with the appropriate sign changes.
  • Table I summarizes the arithmetic for eight phase rotations on the unit circle.
  • the r sample sum-of-products for each term is calculated only once.
  • Each term, then reduced to a single complex value assumes the real value or one of four rotations. This operation simply involves selecting either the real of imaginary value of the complex number and multiplying it by + l or - 1.
  • the partial log-likelihoods are assigned as metrics to a 16 state trellis shown in Figure 2. Four paths with unique metrics enter each state. Four paths with unique metrics also originate at each state. For example, a channel having a multipath spread of three symbols values affect each partial log-likelihood window. Viterbi decoding is applied to the trellis to determine the path with the largest accumulated partial log likelihood.
  • a path is traced back to reveal the most likely symbol sequence for the interval by decoding a symbol at each symbol instant in the traceback operation.
  • the symbol sequence is then converted to differential symbol values for final DQPSK bit detection.
  • a circuit useful in generating the 64 partial log likelihood values, which are referred to as incremental metrics of branch metrics, for each symbol instant employing equation (2) in accordance with the invention is shown in Figure 5.
  • the past and present channel estimates, as well as an approximation of a future channel estimate are provided at 500, 502 and 504.
  • the reconstructed symbol is created for a single integration window.
  • the 64 incremental metrics are developed from the integration window.
  • the three taps of the channel impulse response are provided in each of blocks 500, 502 and 504.
  • Both the in-phase I, and quadrature, Q, components of a different one of the channel impulse response taps provides an input to each rotate block.
  • the odd indexed channel estimates are rotated back by ⁇ /4, in respective rotate blocks 506, 508 and 510, to provide an expected value at 0, 2 ⁇ /4, etc. Rotation is accomplished by substitution using TABLE I. Since only the odd channel estimates are rotated, at one symbol instant ROTATE blocks 506 and 510 will rotate channel estimates, and at previous and subsequent symbol instants only ROTATE block 508 will rotate the channel estimate.
  • Each rotate block generates in-phase and quadrature outputs of the partial response.
  • Outputs 512 and 514 from rotate block 506 are provided to both complex conjugate generator 524 and complex multiply 526.
  • Complex conjugate generator 524 produces the complex conjugate of its inputs 512 and 514 as its outputs 526 and 528 and provides the outputs as two inputs to complex multiply 530.
  • the other two inputs to complex multiply 530 are outputs 516 and 518 from rotate block 508.
  • Complex multiply 530 computes the complex product of its inputs and provides in-phase and quadrature outputs 532 and 534, which are stored in registers 564 and 566.
  • Outputs 516 and 518 from rotate block 508 are provided to both complex conjugate generator 536 and to complex multiply 538.
  • the complex conjugate of inputs 516 and 518 are provided as outputs 540 and 542 from complex conjugate generator 536 which are two inputs to complex multiply 544.
  • the other two inputs to complex multiply 544 are outputs 520 and 522 from rotate block 510.
  • Complex multiply 544 computes the complex product of its inputs and provides in-phase and quadrature outputs 546 and 548 which are stored in registers 568 and 570.
  • Outputs 520 and 522 from rotate block 510 are provided to both complex conjugate generator 550 and to complex multiply 552.
  • the complex conjugate of inputs 520 and 522 are provided as outputs 554 and 556 from complex conjugate generator 550 which are two inputs to complex multiply 558.
  • the other two inputs to complex multiply 558 are outputs 512 and 514 from rotate block 506.
  • Complex multiply 558 computes the complex product of its inputs and provides in-phase and quadrature outputs 560 and 562 which are stored in registers 572 and 574.
  • Complex multiplys 526, 538 and 552 each multiply the inputs from respective rotate blocks by sampled received signal y n * , and respectively store the complex product in registers 576 and 578, 580 and 582, and 584 and 586.
  • complex multiply function represented by complex multiply blocks 526, 530, 538, 544, 552 and 558 can be performed by a single multiplier. Similarly, other functions represented in Figure 5 need not be performed by separate hardware.
  • State table control 588 shown in Table II, controls multiplexers 590-600 by select inputs 602-612. Each multiplexer output 614-624 is provided to a respective selector 626-636. Selectors 626-636 are also controlled by state table control 588 by select input 638-648. Each selector 626-636 provides an output 650-660 to adder 662 which adds each of the inputs to produce output 664.
  • the control of multiplexers 590-600 and selectors 626-636 is represented in Table II.
  • the first column is the metric counter.
  • the second, third and fourth columns are the symbol indices.
  • the six columns at the right represent the twelve controls 602-612 and 638-648.
  • the I or Q in the first of the six right-most columns represents whether the in-phase or quadrature output from registers 576 or 578 is selected as the input to multiplexer 590.
  • the sign preceding the entry in the first of the six right-most columns represents whether the selector 626 selects the output of the multiplexer taken as positive for a +, or inverted for a - .
  • the other five of the six right-most columns each control one of the other multiplexer-selector combinations.
  • control information in the six right-most columns of Table II to control multiplexers 590-600 and selectors 626-636 can be stored in an array with each branch metric requiring only twelve control bits, two for each term of equation (2).
  • the data in Table II is generated from the constellation of Figure 3, the 0, ⁇ /4, 2 ⁇ /4 and 6 ⁇ /4 rotations in Table I, and the real part of the possible symbol product combination of S 1 * S 2 .
  • Each factor of each partial response in the auto correlation is constructed by a one of the symbol values times a channel estimate.
  • components of the cross correlation terms and auto correlation terms of a group of branch metrics at a symbol instant are calculated once and stored in respective registers 564-586.
  • Selective ones of the components of the cross correlation terms and auto correlation terms are combined by multiplexers 590-600 and selectors 626-636 to form the cross correlation terms and auto correlation terms.
  • State table control 588 determines which components from registers 564-586 are combined and whether the components are positive or negative.
  • Outputs 650, 652 and 654 are the cross correlation terms in equation (2) and outputs 656, 658 and 660 are the autocorrelation terms in equation (2).
  • Adder 662 combines the cross correlation terms and auto correlation terms to produce a branch metric as output 664.
  • the metric counter is incremented and another branch metric is generated from the components stored in registers 564-586. This process continues until all of the branch metrics are generated that are based on the components stored in registers 564-586.
  • the branch metrics for all branches into states at a symbol instant are calculated based on the components stored in registers 564-586.
  • the channel impulse response in blocks 500, 502 and 504 are updated as is y * n representing the next symbol instant. Updated components are calculated and stored in registers 564-586, such as by overwriting.
  • the branch metrics for the updated symbol instant are then computed in the above manner. The process is repeated for all symbol instants represented by the received data.
  • Other aspects of the viterbi decoding, including traceback operations and determining final decoded symbols may be done in any known manner.
  • 64 branch metrics are produced as various combinations of the values stored in registers 564-586.
  • 64 branch metrics are produced at output 664, one at a time. Calculations are updated for the next symbol instant, and 64 branch metrics at the new symbol instant are calculated.
  • equation (2) There are six terms in equation (2).
  • the number of terms in equation (2) is independent of the value of r, the number of samples per symbol. Thus, the number of registers 564-586 will remain the same regardless of the number of samples per symbol.
  • the circuit of Figure 5 can be simplified.
  • the simplification is possible because the expected values of received signals take on only two real values.
  • One of the ways the circuit of Figure 5 is simplified is that multiplexers 590-600 are not required.
  • Table III is a tabulations of the number of operations required to calculate branch metrics on a per symbol basis using the Euclidean distance method.
  • L represents the channel length, that is the number of taps in the channel impulse response, and r represents the number of samples per symbol.
  • TABLE III Operation Adds/Symbol Multiples/Symbol Combine complex channel impulse responses 4 ⁇ r ⁇ 4 L Not applicable Calculate error vector 2 ⁇ r ⁇ 4 L Not applicable Square error vectors Not applicable 2(4) L Combine squared error vectors 4 L Not applicable Total (6 ⁇ r+l)4 L 2(4) L
  • Table IV is a tabulation of the number of operations required to calculate branch metrics on a per symbol basis in accordance with the present invention.
  • the computation and power savings of the present invention are beneficial up to about eight channel impulse response tape at one sample per symbol.
  • the savings increase for multiple samples per symbol, and thus the present invention becomes more advantageous even with a greater number of taps in the channel impulse response.
  • the power savings are greater for the present invention.
  • Multiply operations consume approximately five times as much power as an odd operation using the same operands.
  • the power savings are accomplished by performing relative few numbers of multiply operations and performing relative more add operations.
  • Table V illustrates the savings in operations and power for the illustrated embodiment.
  • An energy unit is that amount of energy consumed in an add operation.
  • the branch metric operation is the most redundant operation in a receiver, thus power saved in each operation, here about 50% reduction, is a significant savings.
  • the invention is particularly useful in communications systems and equipment employing integrated circuits including this reduced computation technique.
  • Such communication systems and equipment have the advantage of reduced complexity of circuitry to accomplish the signal processing. Furthermore, the reduced computation time reduces power consumption.
  • control bits can be generated in other ways, such as but not limited to a counter and logic circuitry, a control processing unit (CPU), a computer program, or stored in read-only-memory (ROM).
  • CPU control processing unit
  • ROM read-only-memory
  • the invention has been described as useful in computing branch metrics, the invention is not limited thereto.
  • the invention may also be used to calculate the difference between branch or path metrics entering each state.
  • the branch or path metric differences cab be used to calculate "soft metrics" for subsequent error correction decoding.
  • the invention is not limited thereto.
  • the invention is useful in any digital transmission and detection application employing forward error correction such a viterbi decoder.
  • other applications include reading digital data from a disk or CD ROM, and communications between earth stations and satellites.
  • Pipelining is achieved by initiating computation with a new dataset before completing computation with a previous set of data. The more latches used in pipelining, the greater the depth of pipelining. Pipelining causes an initial latency in computation time required to fill the pipeline, but maximizes usage of resources.

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Probability & Statistics with Applications (AREA)
  • Theoretical Computer Science (AREA)
  • Power Engineering (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Error Detection And Correction (AREA)
  • Complex Calculations (AREA)
  • Editing Of Facsimile Originals (AREA)

Abstract

There is disclosed a method and appartus for computing branch metrics that are combinations of autocorrelation and cross correlation terms are constructed by first calculating and storing components of the autocorrelation and cross correlation terms at a symbol instant. Once calculated and stored, predetermined ones of the autocorrelation components and the cross correlation components are selected. The selected autocorrelation components and cross correlation components, or their inverse as predetermined, are combined to produce a branch metric. Other predetermined combinations of the stored components of autocorrelation and cross correlation terms, or their inverse, are combined to produce other branch metrics at the same symbol instant. All branch metrics associated with the symbol instant can be calculated in this manner. The components of the autocorrelation terms and cross correlation terms associated with the next symbol instant are calculated and stored, and the process of generating the branch metrics repeated.

Description

    Technical Field
  • This invention relates generally to viterbi decoders such as would be used in digital receivers, and particularly to a viterbi decoder with reduced metric computations.
  • Background of the Invention
  • Although telecommunications technology originally developed using copper wires as the transmission media, today there is much interest and commercial activity in wireless technology using free space as the transmission media. Wireless technology frees the user from the constraint of needing a physical connection to the transmission media.
  • Initial commercial activity developing wireless technology used analog transmission technology. In the interest of increasing efficiency of spectrum utilization, digital transmission systems have been developed. One digital standard (Group Special Mobile, or GSM) has been adopted for use in Europe. A different standard (IS54) will likely be adopted for use in North America: either a time-division multiple access (TDMA) standard, or a competing code-division multiple access (CDMA) standard. Other parts of the world such as Japan are considering yet another digital standard, such as the Japan Digital Standard.
  • Digital systems transmit information either as zeros or ones. While accurate transmission of information is generally not at issue, accurate detection of information received depends on being able to detect transmitted ones as ones upon being received, and being able to detect transmitted zeros as zero upon being received. Accurate detection is not a simple matter due to many factors. These factors include low received signal power at a receiver, movement of the transmitter relative to the receiver which introduces changing channel characteristics. Multi-path interference, or inter-symbol interference, occurs when reflections from various objects, including atmospheric layers, cause the presence of multiple transmission paths from the transmitter to the receiver, with varying transit times for a transmitted signal, causing several images of a transmitted signal to arrive at the receiver at different instants in time.
  • Error correction techniques have been developed to increase accuracy of digital communication systems. A viterbi decoder is a maximum likelihood decoder that provides forward error correction. Viterbi decoders are used to decode a sequence of encoded symbols, such as a bit stream. The bit stream can represent encoded information in a telecommunications system. Such encoded information can be transmitted through various media with each bit (or set of bits) representing a symbol instant. In the decoding process, a viterbi decoder works back through a sequence of possible received bit sequences at each symbol instant to determine which one bit sequence was most likely to have been transmitted. The possible transition from a bit at one symbol instant, or state, to a bit at a next or a subsequent, symbol instant or state, is limited. Each possible transition from one state to a next state can he shown graphically and is defined as a branch. All possible branches constitute a trellis. A sequence of interconnected branches is defined as a path. Each state can transition only to a limited number (for example two or four) of next states upon receipt of the next bit (or set or bits) in the bit branch. Thus, some paths survive and other paths do not survive during the decoding process. By eliminating those transitions that are not permissible, computation efficiency can be achieved in determining the most likely path to survive. The viterbi decoder typically defines and calculates a branch metric associated with each branch and employs this branch metric to determine which paths survive and which paths do not survive.
  • A branch metric is calculated at each symbol instant for each possible branch. Each path has an associated metric, accumulated cost, that is updated at each symbol instant. For each possible transmission, the accumulated cost for the next state is calculated as the extremum of the sum of the branch metric for the possible transmission and the path accumulated cost at the respective previous, or originating, state. Calculation at each symbol instant of the branch metric involves an add, compare and select operation for each state. This can result in computation-intensive operation.
  • A need remains in the art for all efficient way to calculate the branch metrics that substantially reduces the number of calculations and concomitantly reduces the power and time consumed in the calculation of branch metrics.
  • Summary of the Invention
  • In accordance with the invention, a method and apparatus for generating branch metrics that are combinations of autocorrelation and cross correlation terms are constructed by first calculating and storing components of the autocorrelation and cross correlation terms at a symbol instant. Once calculated and stored, predetermined ones of the autocorrelation components and the cross correlation components are selected. The selected autocorrelation components and cross correlation components, or their inverse as predetermined, are combined to produce a branch metric. Other predetermined combinations of the stored components of autocorrelation and cross correlation terms, or their inverse, are combined to produce other branch metrics at the same symbol instant. All branch metrics associated with the symbol instant can be calculated in this manner. The components of the autocorrelation terms and cross correlation terms associated with the next symbol instant are calculated and stored, and the process of generating the branch metrics repeated.
  • Brief Description of the Drawings
  • The invention will now be described by way of example in which:
    • Figure 1 is a block diagram of a transceiver, in accordance with an illustrative embodiment of the present invention;
    • Figure 2 is a diagram showing three symbol instants and the possible branches among the individual states;
    • Figure 3 is a graphical representation of an odd symbol phase constellation for π/4- shifted differentially encoded quadrature phase shift keying;
    • Figure 4 is a graphical representation of an even symbol phase constellation for π/4- shifted differentially encoded quadrature phase shift keying; and
    • Figure 5 is a block diagram of a circuit for calculating branch metrics in a viterbi decoder, using reduced computations in accordance with the present invention; and
    • Figure 6 is a graphical comparison of the calculation of branch metrics in accordance with the present invention, as compared to the prior art.
    Detailed Description
  • Transceiver 10 shown in the block diagram of Figure 1 is an illustrative embodiment of the invention. Transceiver 10 is a communication system or portion of a communication system such as a mobile digital cellular telephone. Transceiver 10 is comprised of transmitter 12 and receiver 14.
  • As shown in Figure 1, a microphone 121 for converting a voice signal into an electrical signal is coupled to a vocoder 122. The vocoder is coupled to the channel encoder 123, which includes reordering and partitioning 1231, interleaving and formatting 1232, and symbol converter 1234. The output of the channel encoder is coupled to a transmitter filter 124. The filter output is coupled to a continuous balanced modulator 125, which in turn is coupled to another filter 126. Filter 126 is coupled to an up converter and power amplifier 127 which is turn is coupled to an antenna 128.
  • On the receiver side, antenna 141 is coupled to a low noise amplifier and down conversion 142, which is coupled to filtering 143. Filtering 143 is coupled to automatic gain control 144 which conditions the signal for input to the Viterbi decoder 145. The Viterbi decoder is coupled to a channel decoder 146 which includes deinterleaving, deformatting, deordering and departitioning 1461. The channel decoder is coupled to a vocoder 147. The vocoder is coupled to a speaker 148.
  • Figure 2 shows a trellis diagram useful in understanding the branch metrics used in a Viterbi decoder. A Viterbi decoder reconstructs a transmitted bit stream from a received bit stream. Perfect digital signal symbols in a sequence are transmitted. However, the received signal is distorted due to noise and interference. As the mobile digital cellular telephone moves through its environment, the characteristics of the channel over which the signals are transmitted and received change. At each state or symbol instant (shown horizontally), there are a number of states (shown vertically) that range from zero up to 4C-l-1 where C is the constraint length. These 4C-l states are referred to as individual states. The portion of a trellis diagram shown in Figure 2 is constructed one symbol instant at a time, as is known in the art. The trellis entries at each symbol instant comprise a cell. At each symbol instant, each individual state of the illustrative embodiment shown in Figure 2 can only transition to four possible individual states at the next symbol instant. Concomitantly, each individual state at a symbol instant has only four possible previous individual states from which a transition to it can originate. Other decoders could have more or less than four transitions from a previous state and to a next state.
  • A branch metric is calculated at each symbol instant for each possible transition from an individual state at one symbol instant, known as a source state, to an individual state at the next symbol instant. Various methods for calculating branch metrics are known in art. The branch metric of all branches transitioning to a given next individual state are calculated, then added to the accumulated cost of the respective originating individual state. This results in four potential accumulated cost sums. A comparison is made of the four potential accumulated cost sums. An extremum branch metric, such as the least of the four sums, is selected as the next state accumulated cost. A transition corresponding to the least sum is the most likely transition from the four possible originating individual states to the given next individual state. The surviving branch data of the most likely transition is stored forming the trellis of Figure 2 as the surviving branch to the given next individual state as is known in the art. The extremum branch metric is the accumulated cost of the next individual state. This process is repeated for each of the next individual states and repeated also for each symbol instant in the sequence.
  • Accumulated cost and the source state for the selected branch can be stored in an array. When the trellis is complete, a decoded symbol at each symbol instant can be obtained by a traceback operation. Starting with the individual state in the last symbol instant having the minimum associated accumulated cost, a single trace back through the trellis will provide a decoded symbol at each symbol instant. At each symbol instant of this traceback operation, the source state is translated into a decoded symbol for that symbol instant, as is known in the art.
  • A viterbi decoder is used with a partial log likelihood estimator to determine the most likely symbol sequence taken from a set of possible sequences over an interval of time. Minimizing the integral squared error between the received sequence and the expected sequence results in the most likely sequence to have been transmitted through a channel. The predicted squared error or log-likelihood between the received signal and an expected sequence (or complex functions is the auto correlation of the predicted error functions, multiplied by a constant. The constant may be ignored when evaluating relative likelihoods.
    Figure imgb0001
    where:
    • y(n) is a received sample sequence
    • S(i) is the predicted symbol at symbol instant or symbol index i
    • h is a channel estimate, estimated at the symbol rate (with index i)
    • i is the symbol index
    • l is the number of symbols in the interval
    • n is I/r
    • N is the number of samples in a sequence, N=l/r
      Figure imgb0002
    • m is the metric, or expected value index
    • r is the number of samples per symbol
    • * represents a complex conjugate.
    Recognizing and taking advantage of conjugate symmetry, and substituting the index j for the index i in the conjugate error factor, the log-likelihood can be rewritten as
    Figure imgb0003
    The expected value (the right term in the right bracket of equation (1)) can be considered a sum of partial responses, each of which is equal to S(i)m multiplied by the channel response h(n) and placed in a time location r·i. For a channel filter, L symbols long, where L is the number of symbols in the channel impulse response, L partial responses overlap. The received signal is processed r samples at a time. The total squared error, which is proportional to the log-likelihood, is the sum of the partial windowed squared error values. For a 4-ary case, each symbol can take on one of four possible values. With three symbols affecting each one-symbol window, each window can have 64 possible partial log-likelihoods. For a rapidly changing channel, it may be necessary to update the channel estimate, h(n), once per symbol.
  • Equation (2) represents the integral of the signal energy minus the difference between the cross correlation of the expected sequence. The three terms of equation (2) represent energy, cross-correlation, and autocorrelation respectively. Terms not dependent on S(i)m need not be included when calculating the relative probabilities of the symbol sequences. Such terms include the energy term and the inner products of the autocorrelation term where i=j. It is only necessary to find the sequence which maximizes the absolute value of the difference of the last two terms in the above equation.
  • The equation for evaluating the partial log-likelihood for a three-tap channel filter, 16-state case are: log likelihood (m) = S(i-l) m . y(n) * . h n-l (k+l) + S(i) m . y(n)* . h n (k) +S(i+l) m . y(n)* . h n+l (k-l) + [S(i-l) . S(i)] m . h n-l (k+l) . h n (k+l) + [S(i-l)* . S(i+l)] m . h n-l (k)* . h (n+l) (k) + [S(i)* . S(i+l)] m . h n (k-l)* . h n+l (k-l)
    Figure imgb0004
    where:
    • h(n-l) (k) is the past channel estimate,
    • h(n) (k) is the present channel estimate, and
    • h(n+l) (k) is an approximation of a future channel estimate.
    When performing the initial timing and channel estimate, a known symbol sequence is compared with the received signal to obtain the most likely channel impulse response. For synchronous detection with the differential quadrature phase shift keying (DQPSK) channel, the known preamble sequence of differential symbols should be first converted into a set of absolute phase values from which an expected value is generated for use in the channel estimation process. The first symbol in the preamble used for the channel estimate should be an even number of symbol back from the first unknown data symbol, and the value of this first phase value prior to the first known symbol used for the channel estimate should be arbitrarily set to zero phase. The DQPSK phase differential symbol values following the first symbol preamble used in the channel estimation should then be converted to absolute phases accordingly. The phase offset of the first received symbol can then be expected to be one of the even indexed phases. [2π n/8; n even]
    Figure imgb0005
    on the unit circle (when using the channel estimate derived from the above-mentioned preample sequence). The expected complex symbol value for the first unknown symbol is then one of the four complex values; 1+0j 0-1j, -1+0j, and 0+1j.
  • Odd indexed symbols, are expected to be one of four odd phases on the unit circle [2π n/8; n odd]
    Figure imgb0006
    By rotating the odd indexed expect partial responses by -π/4, the MLSE operation remains identical for both even and odd indexed symbols, and eliminates the need for performing a scalar operation for each state in a symbol interval. The - π/4 phase rotation requires calculating the sum and the difference of the real and imaginary value of the expected values and multiplying then by √2, and re-ordering these scaled sums and differences in four possible arrangements with the appropriate sign changes. Table I summarizes the arithmetic for eight phase rotations on the unit circle. TABLE I
    Θ in radians Re[e(Rc A + j Im A)] Im[e(Re A + j Im A)]
    0 Re A Im A
    π/4 √2 (Rc A-Im A) √2 (Re A+Im A)
    2π/4 -Im A Re A
    3π/4 √2 (-Re A-Im A) √2 (Re A +Im A)
    4π/4 -Re A -Im A
    5π/4 √2 (-Re A+Im A) √2 (-Re A-Im A)
    6π/4 Im A -Re A
    7π/4 √2 (Re A+Im A) √2 (-Rc A+Im A)
    The combination of appropriately chosen references greatly simplifies the calculation of the partial log-likelihoods. For example, the r sample sum-of-products for each term is calculated only once. Each term, then reduced to a single complex value, assumes the real value or one of four rotations. This operation simply involves selecting either the real of imaginary value of the complex number and multiplying it by + l or - 1. The partial log-likelihoods are assigned as metrics to a 16 state trellis shown in Figure 2. Four paths with unique metrics enter each state. Four paths with unique metrics also originate at each state. For example, a channel having a multipath spread of three symbols values affect each partial log-likelihood window. Viterbi decoding is applied to the trellis to determine the path with the largest accumulated partial log likelihood. Starting with the state in the last symbol instant of the interval having the largest accumulated partial log-likelihood, a path is traced back to reveal the most likely symbol sequence for the interval by decoding a symbol at each symbol instant in the traceback operation. The symbol sequence is then converted to differential symbol values for final DQPSK bit detection.
  • A circuit useful in generating the 64 partial log likelihood values, which are referred to as incremental metrics of branch metrics, for each symbol instant employing equation (2) in accordance with the invention is shown in Figure 5. For a channel impulse response having three terms, that is a channel filter having three taps, and one sample per symbol, in equation (2) r=l, i=N and l=N.
  • The past and present channel estimates, as well as an approximation of a future channel estimate are provided at 500, 502 and 504. At each symbol instant, the reconstructed symbol is created for a single integration window. The 64 incremental metrics are developed from the integration window. The three taps of the channel impulse response are provided in each of blocks 500, 502 and 504. Both the in-phase I, and quadrature, Q, components of a different one of the channel impulse response taps provides an input to each rotate block. In a preferred embodiment, the odd indexed channel estimates are rotated back by π/4, in respective rotate blocks 506, 508 and 510, to provide an expected value at 0, 2π/4, etc. Rotation is accomplished by substitution using TABLE I. Since only the odd channel estimates are rotated, at one symbol instant ROTATE blocks 506 and 510 will rotate channel estimates, and at previous and subsequent symbol instants only ROTATE block 508 will rotate the channel estimate.
  • Each rotate block generates in-phase and quadrature outputs of the partial response. Outputs 512 and 514 from rotate block 506 are provided to both complex conjugate generator 524 and complex multiply 526. Complex conjugate generator 524 produces the complex conjugate of its inputs 512 and 514 as its outputs 526 and 528 and provides the outputs as two inputs to complex multiply 530. The other two inputs to complex multiply 530 are outputs 516 and 518 from rotate block 508. Complex multiply 530 computes the complex product of its inputs and provides in-phase and quadrature outputs 532 and 534, which are stored in registers 564 and 566.
  • Outputs 516 and 518 from rotate block 508 are provided to both complex conjugate generator 536 and to complex multiply 538. The complex conjugate of inputs 516 and 518 are provided as outputs 540 and 542 from complex conjugate generator 536 which are two inputs to complex multiply 544. The other two inputs to complex multiply 544 are outputs 520 and 522 from rotate block 510. Complex multiply 544 computes the complex product of its inputs and provides in-phase and quadrature outputs 546 and 548 which are stored in registers 568 and 570.
  • Outputs 520 and 522 from rotate block 510 are provided to both complex conjugate generator 550 and to complex multiply 552. The complex conjugate of inputs 520 and 522 are provided as outputs 554 and 556 from complex conjugate generator 550 which are two inputs to complex multiply 558. The other two inputs to complex multiply 558 are outputs 512 and 514 from rotate block 506. Complex multiply 558 computes the complex product of its inputs and provides in-phase and quadrature outputs 560 and 562 which are stored in registers 572 and 574.
  • Complex multiplys 526, 538 and 552 each multiply the inputs from respective rotate blocks by sampled received signal yn *, and respectively store the complex product in registers 576 and 578, 580 and 582, and 584 and 586.
  • The complex multiply function represented by complex multiply blocks 526, 530, 538, 544, 552 and 558 can be performed by a single multiplier. Similarly, other functions represented in Figure 5 need not be performed by separate hardware.
  • State table control 588, shown in Table II, controls multiplexers 590-600 by select inputs 602-612. Each multiplexer output 614-624 is provided to a respective selector 626-636. Selectors 626-636 are also controlled by state table control 588 by select input 638-648. Each selector 626-636 provides an output 650-660 to adder 662 which adds each of the inputs to produce output 664.
    Figure imgb0007
    Figure imgb0008
    Figure imgb0009
  • The control of multiplexers 590-600 and selectors 626-636 is represented in Table II. The first column is the metric counter. The second, third and fourth columns are the symbol indices. The six columns at the right represent the twelve controls 602-612 and 638-648. The I or Q in the first of the six right-most columns represents whether the in-phase or quadrature output from registers 576 or 578 is selected as the input to multiplexer 590. The sign preceding the entry in the first of the six right-most columns represents whether the selector 626 selects the output of the multiplexer taken as positive for a +, or inverted for a - . The other five of the six right-most columns each control one of the other multiplexer-selector combinations.
  • The control information in the six right-most columns of Table II to control multiplexers 590-600 and selectors 626-636 can be stored in an array with each branch metric requiring only twelve control bits, two for each term of equation (2). The data in Table II is generated from the constellation of Figure 3, the 0, π/4, 2π/4 and 6π/4 rotations in Table I, and the real part of the possible symbol product combination of S1 * S2. Each factor of each partial response in the auto correlation is constructed by a one of the symbol values times a channel estimate.
  • In accordance with the invention, components of the cross correlation terms and auto correlation terms of a group of branch metrics at a symbol instant are calculated once and stored in respective registers 564-586. Selective ones of the components of the cross correlation terms and auto correlation terms are combined by multiplexers 590-600 and selectors 626-636 to form the cross correlation terms and auto correlation terms. State table control 588 determines which components from registers 564-586 are combined and whether the components are positive or negative. Outputs 650, 652 and 654 are the cross correlation terms in equation (2) and outputs 656, 658 and 660 are the autocorrelation terms in equation (2). Adder 662 combines the cross correlation terms and auto correlation terms to produce a branch metric as output 664. The metric counter is incremented and another branch metric is generated from the components stored in registers 564-586. This process continues until all of the branch metrics are generated that are based on the components stored in registers 564-586. Typically, the branch metrics for all branches into states at a symbol instant are calculated based on the components stored in registers 564-586. The channel impulse response in blocks 500, 502 and 504 are updated as is y * n
    Figure imgb0010
    representing the next symbol instant. Updated components are calculated and stored in registers 564-586, such as by overwriting. The branch metrics for the updated symbol instant are then computed in the above manner. The process is repeated for all symbol instants represented by the received data. Other aspects of the viterbi decoding, including traceback operations and determining final decoded symbols may be done in any known manner.
  • Once the calculations are updated and stored in registers 564-586, 64 branch metrics are produced as various combinations of the values stored in registers 564-586. By stepping through the 64 controls of state table control 588, 64 branch metrics are produced at output 664, one at a time. Calculations are updated for the next symbol instant, and 64 branch metrics at the new symbol instant are calculated.
  • There are six terms in equation (2). The number of terms in equation (2) is independent of the value of r, the number of samples per symbol. Thus, the number of registers 564-586 will remain the same regardless of the number of samples per symbol.
  • For GSM applications, the circuit of Figure 5 can be simplified. The simplification is possible because the expected values of received signals take on only two real values. One of the ways the circuit of Figure 5 is simplified is that multiplexers 590-600 are not required.
  • Table III is a tabulations of the number of operations required to calculate branch metrics on a per symbol basis using the Euclidean distance method. L represents the channel length, that is the number of taps in the channel impulse response, and r represents the number of samples per symbol. TABLE III
    Operation Adds/Symbol Multiples/Symbol
    Combine complex channel impulse responses 4·r·4L Not applicable
    Calculate error vector 2·r·4L Not applicable
    Square error vectors Not applicable 2(4)L
    Combine squared error vectors 4L Not applicable
    Total (6·r+l)4L 2(4)L
  • Table IV is a tabulation of the number of operations required to calculate branch metrics on a per symbol basis in accordance with the present invention.
    Figure imgb0011
  • The total number of operations to complete a metric calculation at one symbol instant for the 4-ary case using each method is tabulated in the bottom row of each table.
  • As compared to the Euclidean distance method, the computation and power savings of the present invention are beneficial up to about eight channel impulse response tape at one sample per symbol. The savings increase for multiple samples per symbol, and thus the present invention becomes more advantageous even with a greater number of taps in the channel impulse response. Even where the number of operations are about the same for the two methods compared in Tables III and IV, the power savings are greater for the present invention. Multiply operations consume approximately five times as much power as an odd operation using the same operands. The power savings are accomplished by performing relative few numbers of multiply operations and performing relative more add operations. TABLE V
    Operation Euclidean Method Reduced Metric Computation Method
    Two Samples per Symbol One Sample per Symbol Two Samples per Symbol One Sample per Symbol
    Adds 832 448 326 320
    Multiplies 128 128 48 24
    Total operations 960 576 374 344
    Energy Units Consumed 1472 1088 566 440
  • Table V illustrates the savings in operations and power for the illustrated embodiment. An energy unit is that amount of energy consumed in an add operation. The branch metric operation is the most redundant operation in a receiver, thus power saved in each operation, here about 50% reduction, is a significant savings.
  • The invention is particularly useful in communications systems and equipment employing integrated circuits including this reduced computation technique. Such communication systems and equipment have the advantage of reduced complexity of circuitry to accomplish the signal processing. Furthermore, the reduced computation time reduces power consumption.
  • While the control of multiplexers 590-600 and selectors 626-636 has been described as being controlled by state table control 588, the invention is not limited thereto. The control bits can be generated in other ways, such as but not limited to a counter and logic circuitry, a control processing unit (CPU), a computer program, or stored in read-only-memory (ROM).
  • Although the invention has been described as useful in computing branch metrics, the invention is not limited thereto. The invention may also be used to calculate the difference between branch or path metrics entering each state. The branch or path metric differences cab be used to calculate "soft metrics" for subsequent error correction decoding.
  • While the illustrative embodiment has been described as applicable to branch metric calculations for a trellis having a constraint length of C, the invention is also applicable to a "reduced constraint" trellis viterbi decoder.
  • While the illustrative embodiment of the invention has been described as useful in digital cellular telephony, the invention is not limited thereto. The invention is useful in any digital transmission and detection application employing forward error correction such a viterbi decoder. By way of example and not limitation, other applications include reading digital data from a disk or CD ROM, and communications between earth stations and satellites.
  • While the illustrative embodiment of the invention has not been described as incorporating pipelining, one skilled in the art would recognize the enhanced computation efficiency available by utilizing pipelining in the design. Pipelining is achieved by initiating computation with a new dataset before completing computation with a previous set of data. The more latches used in pipelining, the greater the depth of pipelining. Pipelining causes an initial latency in computation time required to fill the pipeline, but maximizes usage of resources.

Claims (21)

  1. A method for determining a plurality of branch metrics at a symbol instant that are comprised of autocorrelation terms and cross correlation terms, including the steps of generating components of the autocorrelation terms, storing the components of the autocorrelation terms, generating components of the cross correlation terms, storing the components of the cross correlation terms,
    characterized by:
    selecting predetermined ones of the autocorrelation components and the cross correlation components, and
    combining the selected autocorrelation components and cross correlation components in a predetermined manner to produce a branch metric.
  2. The method as recited in claim 1, further characterized the step of:
    repeating the selecting and combining steps to produce another branch metric.
  3. An integrated circuit including a branch metric calculation circuit, including a circuit (526, 538, 552) for generating components of autocorrelation terms, at least one register (576, 578, 580, 582, 584, 586) for receiving and storing components of the autocorrelation terms, a circuit (530, 544, 558) for generating components of cross correlation terms, at least one register (564, 566, 568, 570, 572, 574) for receiving and storing components of the cross correlation terms,
       characterized by:
    a selection circuit (590, 592, 594, 596, 598, 600) for selecting predetermined ones of the autocorrelation components and the cross correlation components stored in said at least one registers (576, 578, 580, 582, 584, 586);
    a selective inverter circuit (626, 628, 630, 632, 634, 636) for selectively inverting predetermined ones of the selected auto correlation components and cross correlation components;
    a combining circuit (662) for receiving the selected auto correlation components and cross correlation components and for combining the selected autocorrelation components and cross correlation components to produce a branch metric.
  4. The integrated circuit as recited in claim 3, further characterized in that the selection circuit includes a look-up table (588).
  5. The integrated circuit as recited in claim 3 or 4, further characterized in that the combining circuit is an adder (662).
  6. The integrated circuit as recited in any of claims 3-5 further characterized in that the selective inverter circuit (626, 628, 630, 632, 634, 636) provides as an output a positive or a negative version of the input depending upon a control select input (632, 634, 636, 638, 640, 642).
  7. The integrated circuit as recited in any of claims 3-6, further characterized in that the selection circuit selects all of the autocorrelation components and all of the cross correlation components.
  8. The integrated circuit as recited in any of claims 3-7, further characterized by at least one rotate circuit (506, 508, 510).
  9. An integrated circuit for having a branch metric calculation circuit including a circuit for generating components of autocorrelation terms, a circuit for generating components of cross correlation terms, and a circuit for combining autocorrelation components and cross correlation components characterized by:
    at least two registers (576, 578, 580, 582, 584, 586) for receiving and storing components of the autocorrelation terms;
    at least two registers (564, 566, 568, 570, 572, 574) for receiving and storing components of the cross correlation terms;
    a selection circuit (590, 592, 594, 596, 598, 600) for selecting predetermined ones of the autocorrelation components and the cross correlation components stored in said at least two registers;
    a selective circuit (626, 628, 630, 632, 634, 636) for selectively inverting predetermined ones of the selected auto correlation components and cross correlation components:
    a combining circuit (662) for receiving the selected auto correlation components and cross correlation components and for combining the selected autocorrelation components and cross correlation components to produce a branch metric.
  10. The integrated circuit as recited in claim 9, further characterized in that the selection circuit includes a look-up table (588).
  11. The integrated circuit as recited in claim 9 or 10, further characterized in that the combining circuit is an adder (662).
  12. The integrated circuit as recited in any of claims 9-11 further characterized in that the selective inverter circuit (626, 628, 630, 632, 634, 636) provides as an output a positive or a negative version of the input depending upon a control select input (632, 634, 636, 638, 640, 642).
  13. The integrated circuit as recited in any of claim 9-12, further characterized in that the selection circuit selects all of the autocorrelation components and all of the cross correlation components.
  14. The integrated circuit as recited in any of claims 9-13, further characterized by at least one rotate circuit (506, 508, 510).
  15. A receiver (14), the receiver (14) including a decoder having a branch metric calculation circuit, including a circuit (526, 538, 552) for generating components of autocorrelation terms, at least one register (576, 578, 580, 582, 584, 586) for receiving and storing components of the autocorrelation terms, a circuit (530, 544, 558) for generating components of cross correlation terms, at least one register (564, 566, 568, 570, 574, 586) for receiving and storing components of the cross correlation terms, a selection circuit (590, 592, 594, 596, 598, 600) for selecting predetermined ones of the autocorrelations components and the cross correlation components stored in said at least one registers (576, 578, 580, 582, 584, 586),
       characterized by:
    a selective inverter (626, 628, 630, 632, 634, 636) circuit for selectively inverting predetermined ones of the selected auto correlation components and cross correlation components, and
    a combining circuit (662) for receiving the selected auto correlation components, and cross correlation components and for combining the selected autocorrelation components and cross correlation components to produce a branch metric.
  16. A receiver as recited in claim 15, further characterized in that the selection circuit includes a look-up table (588).
  17. A receiver as recited in claim 15 or 16, further characterized in that the combining circuit is an adder (662).
  18. A receiver as recited in any of claims 15-17 further characterized in that the selective inverter circuit (626, 628, 630, 632, 634, 636) provides as an output a positive or a negative version of the input depending upon a control select input (632, 634, 636, 638, 640, 642).
  19. A receiver as recited in any of claims 15-18, further characterized in that the selection circuit selects all of the autocorrelation components and all of the cross correlation components.
  20. A receiver as recited in any of claims 15-19, further characterized by at least one circuit (506, 508, 510).
  21. A receiver (14) as recited in any of claims 15-20, further characterized by a transmitter (12).
EP96309397A 1995-12-29 1996-12-23 Viterbi decoder with reduced metric computation Withdrawn EP0786872A2 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US953795P 1995-12-29 1995-12-29
US9537 1995-12-29

Publications (1)

Publication Number Publication Date
EP0786872A2 true EP0786872A2 (en) 1997-07-30

Family

ID=21738272

Family Applications (1)

Application Number Title Priority Date Filing Date
EP96309397A Withdrawn EP0786872A2 (en) 1995-12-29 1996-12-23 Viterbi decoder with reduced metric computation

Country Status (3)

Country Link
EP (1) EP0786872A2 (en)
JP (1) JPH104362A (en)
TW (1) TW312070B (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0901234A2 (en) * 1997-09-08 1999-03-10 Lucent Technologies Inc. Viterbi compare-select operation for two-bit traceback coding
US6353913B2 (en) 1998-03-23 2002-03-05 Nokia Networks Oy Modulation detection method and apparatus
WO2002039689A2 (en) * 2000-11-13 2002-05-16 Telefonaktiebolaget Lm Ericsson (Publ) Methods and systems for accumulating metrics generated by a sequence estimation algorithm
EP1432190A1 (en) * 2002-12-18 2004-06-23 Texas Instruments Incorporated RSSE using hardware acceleration

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
None

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0901234A2 (en) * 1997-09-08 1999-03-10 Lucent Technologies Inc. Viterbi compare-select operation for two-bit traceback coding
EP0901234A3 (en) * 1997-09-08 2004-06-16 Lucent Technologies Inc. Viterbi compare-select operation for two-bit traceback coding
US6353913B2 (en) 1998-03-23 2002-03-05 Nokia Networks Oy Modulation detection method and apparatus
WO2002039689A2 (en) * 2000-11-13 2002-05-16 Telefonaktiebolaget Lm Ericsson (Publ) Methods and systems for accumulating metrics generated by a sequence estimation algorithm
WO2002039689A3 (en) * 2000-11-13 2002-08-01 Ericsson Telefon Ab L M Methods and systems for accumulating metrics generated by a sequence estimation algorithm
US6970520B1 (en) 2000-11-13 2005-11-29 Telefonaktiebolaget Lm Ericsson (Publ) Methods and systems for accumulating metrics generated by a sequence estimation algorithm
EP1334562B1 (en) * 2000-11-13 2012-01-04 TELEFONAKTIEBOLAGET LM ERICSSON (publ) Simplifications of the metric computation within a maximum likelihood sequence estimation algorithm
EP1432190A1 (en) * 2002-12-18 2004-06-23 Texas Instruments Incorporated RSSE using hardware acceleration
US7372922B2 (en) 2002-12-18 2008-05-13 Texas Instruments Incorporated RSSE optimization using hardware acceleration

Also Published As

Publication number Publication date
JPH104362A (en) 1998-01-06
TW312070B (en) 1997-08-01

Similar Documents

Publication Publication Date Title
US5844947A (en) Viterbi decoder with reduced metric computation
US5471500A (en) Soft symbol decoding
US5550870A (en) Viterbi processor
US5774504A (en) Equalization and decoding for digital communication channel
JP3822249B2 (en) Method and means for detecting communication signals with non-uniform error protection
EP1290821B1 (en) Methods and systems for decoding a received signal having a transmitter or channel induced coupling between bits
US5502735A (en) Maximum likelihood sequence detector
US5375129A (en) Maximum likelihood sequence detector
US8320441B2 (en) Efficient high performance demodulation of low BT value gaussian minimum shift keying incorporating turbo equalization
EP1424781A1 (en) Information transmission apparatus; information transmission method; information reception apparatus; and information reception method
EP0632623A2 (en) Adaptive maximum likelihood sequence estimator
EP0653715B1 (en) Integrated circuit comprising a coprocessor for Viterbi decoding
US20030115061A1 (en) MPSK equalizer
WO2017039558A1 (en) Method and apparatus for simplified generation of continuous phase modulation, cpm, waveforms
US5724390A (en) MLSE before derotation and after derotation
US7277506B1 (en) Maximum likelihood sequence estimator which computes branch metrics in real time
JP3728171B2 (en) Reliability information calculation method
FI102230B (en) Reception method and receiver
US5533065A (en) Decreasing length tracebacks
EP0786872A2 (en) Viterbi decoder with reduced metric computation
CN108055106B (en) Recognition method for code word boundary based on RM code partial decoding auxiliary differential detection
JPH05152894A (en) Adaptive equalizer
CN114598423B (en) Method, device and medium for demodulation and decoding combined GMSK and LDPC
RU184852U1 (en) A device for generating a soft solution of a signal demodulator of a Quaternary phase shift keying with smoothing and shift
JPH09266458A (en) Adaptive equalizer

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

AK Designated contracting states

Kind code of ref document: A2

Designated state(s): DE ES FR GB IT

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

Free format text: STATUS: THE APPLICATION HAS BEEN WITHDRAWN

18W Application withdrawn

Withdrawal date: 19990614