DELAYED DECISION FEEDBACK LOG-MAP EQUALIZER

Abstract: A delayed decision feedback Log-MAP equalizer generates optimal soft information at the bit level for MPSK and MDPSK modulations over multipath fading channels. The algorithms have a forward-backward recursion structure similar to conventional Log-MAP algorithms. To reduce complexity, the forward recursion incorporates the techniques of delayed decision feedback sequence estimation to construct a reduced state space. The backward recursion then takes this state space as given, in that it does not hypothesize feedbacks, and computes required soft information.

Description: DELAYED DECISION FEEDBACK LOG-MAP EQUALIZER BACKGROUND OF THE INVENTION The present invention relates generally to cellular and wireless communication and, more particularly, to an equalizer hypothesizing modulation symbols received over multipath fading channels. Signals in a wireless communication system are subject to a number of phenomena that degrade signal quality. Each signal is reflected from many different man-made and natural objects. The receiver thus receives a number of signals delayed by one or more signal periods, called "multipath", as each reflected signal is received. If the period of delay is more than the time required to transmit one symbol, producing intersymbol interference (ISI), then a receiver decoder may not be able to decode the symbol. This can cause poor sound or data quality to the user. Many different algorithms are used at the receiver to attempt to compensate for such effects. These techniques are discussed, for example, by G. L. Stϋber in Principles of Mobile Communication, Chapter 6, Kluwer Academic Publishers, 1996. One such algorithm is the maximum likelihood sequence estimation (MLSE) or "Viterbi" algorithm. Digital cellular and personal communication systems based on IS- 136 and GSM (global system for mobile communications) require equalization to handle ISI arising from time dispersion. With the advent of the high date rates and the high level modulation in EDGE (enhanced data rates for global evolution) systems, the extent of ISI has increased considerably. To cover the same delay spread as today's GSM, EDGE receivers are specified to resolve a 5-tap channel. Furthermore, these systems are convolutional coded. To enable soft-decision decoding, it is, therefore, highly desirable for the equalizers to produce soft information at the bit level. For ΛJ-ary modulation, a traditional MLSE equalizer requiresΛJ 1 at each trellis stage in order to properly model an intersymbol interference (ISI) channel driven by L + 1 M- ary phase-shift keying (MPSK) symbols Because of the high computational and storage requirements this implies, the delayed decision feedback sequence estimation (DDFSE) equalizer uses a kind of hybrid of MLSE and DFE (decision feedback estimation) in which only theK+ 1 most recent symbols are hypothesized as part of the state model, while the remaining L - K symbols required to compute each path metrics are determines from the surviving paths of each state Therefore, a trellis of only states is needed in this case and the L - K "DFE" symbols constitute side information of the states DDFSE equalizers are discussed, for example, by A Duel-Hallen and C Heegard in "Delayed Decision Feedback Sequence Estimation," IEEE Transactions on Communications, vol 37, no 5, pp 428-436, May 1989 One convenient way of describing the DDFSE operation is to say that the energy in the most recent K + 1 "MLSE" symbols is used to equalize the data, while the energy in the remaining "DFE" symbols is simply canceled The performance of DDFSE equalizers can be improved by pre-filtering the received signals with minimum-phase filters to maximize the energy of the leading ISI taps Such filtering techniques are discussed, for example, in copending Application Serial No 09/378,314 entitled "Method and Apparatus for Computing Prefilter Coefficients" of K C Zangi et al and filed Aug 20, 1999 In a convolutional coded system, a Log-MAP (maximum aposterion) equalizer is preferred because of its ability to generate high-quality soft information for the Viterbi decoder Traditional MAP algorithms, as described in L R Bahl, J Cocke, F Jelinek and J Raviv (BCJR) in "Optimal Decoding of Linear Codes for Minimizing Symbol Error Rate," IEEE Transactions on Information Theory, vol 20, pp 284-287, Mar 1974, are not very suited to practical implementation because these algorithms require many multiplication and logarithm operations The Log-MAP algorithms circumvent this complexity problem by operating the signals in the logarithmic domain, as discussed, for example, by P Robertson, E Villebrun and P Hoehner, in "A Comparison of Optimal and Sub-Optimal MAP decoding algorithms operating in the log domain," Proceedings of IEEE International Communications Conference '95, pp 1009-1013, June 1995 However, both approaches are designed to generate optimal soft information for trellis with static state spaces, i e , the definitions of the states do not change over time Moreover, these definitions are known a priori and are usually given exogenously For example, the trellis of a convolutional code falls into this category the states are defined by the contents of the encoding shift-registers and the complete trellis can be constructed even before the codewords are transmitted For multipath equalization, the traditional MAP algorithms thus require the construction of the complete state space as in the case for MLSE equalizers For many practical applications, this proves to be excessively complex For example, a MAP equalizer for 8PSK with 5 channel taps (EDGE specifications) requires 4096 states Since the Log-MAP equalizers suffer from the same complexity problems as the MLSE equalizers, it would be of great practical value if the principle of DDFSE could be combined with the Log-MAP algorithms to reduce complexity An approach to this problem is discussed in the following Several problems arise when this MLSE/DFE hybrid method is applied to a forward-backward Log-MAP equalizer More specifically, the forward DDFSE recursion would use an anti-causal prefilter to create a minimum-phase channel response ADDFSE run backwards through the trellis, on the other hand, would need a causal prefilter to crease a maximum-phase channel, because in this case the MLSE symbols appear at the other end of the true state, as illustrated in Fig 1 This results in two major issues 1 The forward-backward Log-MAP algorithm breaks down The two passes are now using two sets of data which are correlated with each other (because of the two prefilters with overlapping support) Because of this correlation, the Markovian property of the data sequence, which means that the future data symbols are independent of the past data symbols given the current state, no longer holds This Markovian property plays an essential role in the derivation of the forward-backward algorithm, without which the Log-MAP algorithm is invalid 2 The state space from the two recursions do not match A complete state space is used in a Log-MAP equalizer, which guarantees the forward and backward recursions work on the same state space On the other hand, for a naive delayed decision feedback Log-MAP equalizer, only the "MLSE" symbols can be matched up The "DFE" symbols setup by the forward and backward recursions generally do not agree with each other The overlapping supports of the two prefilters makes the situation even more confusing The definitions of the states in a DDFSE trellis are constructed endogenously by the equalization recursion itself and, therefore, become unknown a priori The traditional MAP algorithms are thus not directly applicable in these situations In order to extract soft information, ad hoc post-processing methods such as the soft-output Viterbi algorithms (SOVA) are usually used instead The SOVA algorithm is discussed, for example, by J Jagenauer and P. Hoehner in "A Viterbi Algorithm with Soft-Decision Outputs and its Applications," Proceedings of IEEE Global Telecommunications Conference '89, pp. 1680- 1686, Nov. 1989. These algorithms, however, are sub-optimal and their shortcomings become more apparent in iterative processing settings. Accordingly, there is a need for an equalizer and method that generates optimal soft information at the bit level for MPSK and MDPSK modulations over multipath fading channels. SUMMARY OF THE INVENTION The present invention meets this need by providing a delayed decision feedback Log-MAP equalizer that generates optimal soft information at the bit level for MPSK and MDPSK modulations over multipath fading channels. Broadly, there is disclosed in accordance with one aspect of the invention a method of hypothesizing modulation symbols received over multipath fading channels. An estimate of a received sequence of symbols is calculated and stored using stored channel coefficients and a previous state path history of modulation symbols. For each possible current state symbol, a branch metric is calculated and stored from a prior state to a current state using the calculated estimate of the received sequence of symbols. A forward recursion is performed by calculating a forward state metric using the branch metric from the prior state to the current state. A backward recursion is performed by calculating a backward state metric using the branch metric from the prior state to the current state. It is a feature of the invention that the estimate of a received sequence of symbols calculated and stored using stored channel coefficients and a previous state path history of modulation symbols is implemented using a delayed decision feedback sequence estimation equalizer. It is another feature of the invention that the estimate of a received sequence of symbols calculated and stored using stored channel coefficients and a previous state path history of modulation symbols is implemented using a maximum likelihood sequence estimation equalizer for most recent symbols in the sequence of symbols and using a decision feedback equalizer for remaining symbols in the sequence of symbols. It is a further feature of the invention that a Log-MAP equalizer is used to calculate and store the branch metric, the forward state metric and the backward state metric. It is yet another feature of the invention that a posteriori log-likelihood ratio of each bit in the hypothesized modulation symbols is calculated using the calculated forward state metric and the calculated backward state metric. It is still another feature of the invention that the forward recursion is executed prior to the backward recursion. It is still an additional feature of the invention that the backward recursion is executed prior to the forward recursion. It is yet still another feature of the invention that the previous state path history of modulation symbols is updated using the calculated branch metric from a prior state to a current state using the calculated estimate of the received sequence of symbols. It is still a further feature of the invention that the modulation symbols received over multipath fading channels comprise M-ary modulated symbols. It is yet a further feature of the invention that the modulation symbols received over multipath fading channels comprise M-ary differential PSK modulated symbols. There is disclosed in accordance with another aspect of the invention, in a receiver receiving modulation symbols transmitted over multipath fading channels, an equalizer including a programmed controller. The programmed controller calculates and stores an estimate of a received sequence of symbols using stored channel coefficients and a previous state path history of modulation symbols. A branch metric is calculated and stored for each possible current state symbol from a prior state to a current state using the calculated estimate of the received sequence of symbols. A forward recursion is performed by calculating a forward state metric using the branch metric from the prior state to the current state. A backward recursion is performed by calculating a backward state metric using the branch metric from the prior state to the current state. There is disclosed in accordance with a further aspect of the invention a method of generating soft information at bit level of modulation symbols received over multipath fading channels. An estimate of a received sequence of symbols is calculated and stored using stored channel coefficients and a previous state path history of modulation symbols, each symbol comprising an M-ary modulated symbol representing select bits c. A branch metric for each possible current state symbol is calculated and stored from a prior state to a current state using the calculated estimate of the received sequence of symbols. A forward recursion is performed by calculating a forward state metric using the branch metric from the prior state to the current state. A backward recursion is performed by calculating a backward state metric using the branch metric from the prior state to the current state. A posteriori log-likelihood ratio of each select bit c in the estimated modulation symbols is calculated using the calculated forward state metric and the calculated backward state metric. Particularly, a DDFSE method is incorporated with a forward recursion to construct a reduced state space. Starting with a complete state space, the X states are partitioned into hX subsets Each of the subsets consists of the AX L"K> states with the same K leading hypothesized symbols while the remaining L - K hypotheses are completely different Since the minimum-phase prefilter shifts most of the energy to the leading taps, the last L - K hypotheses have little effect on the recursion metrics To reduce computation, L - K decision feedbacks are then used to select one representative state for each subset The same reduced state space is used by the backward recursion Using this structure, the decision feedback side information for the backward recursion is not associated with the surviving paths in the backward trellis, but rather is determined by the forward recursion to represent backward prediction symbols along each path That is, as shown in Fig 2, the DFE symbols appear on the left of the MLSE symbols for both forward and backward recursions, meaning that a minimum-phase channel is appropriate for both directions Therefore, both directions can use the same prefilter and the two concerns described in the previous section no longer apply The disclosed invention generalizes the traditional BCJR paradigm to this reduced-state trellis, which is referred to as a DDF-Log-MAP (delayed decision feedback Log- MAP) algorithm The disclosed method has a forward-backward recursion structure similar to the classic Log-MAP algorithms To reduce complexity, however, the forward recursion incorporates DDFSE to construct a reduced state space The backward recursion then takes this state space as given, i e , it does not hypothesize feedbacks, and computes the required soft information The description includes exemplary applications of the invention to MPSK and M-ary differential phase-shift keying (MDPSK) modulations, respectively, with a reduced state space equivalent to that of a DDFSE with 2 MLSE taps Generalizations to other modulations and other choices of reduced state spaces are straightforward Further features and advantages of the invention will be readily apparent from the specification and the drawings. BRIEF DESCRIPTION OF THE DRAWINGS Fig. 1 illustrates state configurations for forward and backward recursions using a naive approach to combining DDFSE and Log-MAP algorithms; Fig. 2 illustrates state configurations for forward and backward recursions using an approach to combining DDFSE and Log-MAP algorithms in accordance with the invention; Fig. 3 is a block diagram of a mobile terminal implementing the method according to the invention; Fig. 4 is a block diagram of a system model for a wireless communication system implementing the method according to the invention; Fig. 5 illustrates trellis notation used in explaining the method according to the invention; and Fig. 6 is a flow diagram illustration of a program run by the processor of Fig. 3 to implement the method according to the invention. DETAILED DESCRIPTION OF THE INVENTION In accordance with the invention, a DDFSE method is incorporated with the forward recursion to construct a reduced state space. Starting with a complete state space, the A/' states are partitioned into Xt subsets. Each of the subset consists of the M* 1' ^ states with the same K leading hypothesized symbols while the remaining L -K hypotheses would be completely different. Since the minimum-phase prefilter shifts most of the energy to the leading taps, the last L - K hypothesis have little effect on the recursion metrics To reduce computation, L -K decision feedbacks are then used to select one representative state for each subset The same reduced state space is used by the backward recursion The decision feedback side information for the backward recursion is not associated with the surviving paths in the backward trellis, but rather is determined by the forward recursion to represent backward prediction symbols along each path That is, as shown in Fig 2, the DFE symbols appear on the left of the MLSE symbols for both forward and backward recursions, meaning that a minimum-phase channel is appropriate for both directions Therefore, both directions can use the same prefilter Fig 3 is a block diagram of a typical mobile terminal 10 The mobile terminal 10 includes an antenna 12, a receiver 16, a transmitter, 18, a speaker 20, a processor 22, a memory 24, a user interface 26 and a microphone 32 The antenna is configured to send and receive radio signals between the mobile terminal 10 and a wireless network (not shown) The antenna 12 is connected to a duplex filter 14 which enables the receiver 16 and the transmitter 18 to receive and broadcast (respectively) on the same antenna 12 The receiver 16, in conjunction with the processor 22, demodulates, demultiplexes and decodes the radio signals into one or more channels Such channels include a control channel and a traffic channel for speech or data The speech or data are delivered to the speaker 20 (or other output device, such as a modem or fax connector) The receiver 16 delivers messages from the control channel to the processor 22. The processor 22 controls and coordinates the functioning of the mobile terminal 10 and is responsive to messages on the control channel using programs and data stored in the memory 24, so that the mobile terminal 10 can operate within the wireless network (not shown). The processor 22 also controls the operation of the mobile terminal 10 and is responsive to input from the user interface 26. The user interface 26 includes a keypad 28 as a user-input device and a display 30 to give the user information. Other devices are frequently included in the user interface 26, such as lights and special purpose buttons. The processor 22 controls the operations of the transmitter 18 and the receiver 16 over control lines 34 and 36, respectively, responsive to control messages and user input. The microphone 32 (or other data input device) receives speech signal input and converts the input into analog electrical signals. The analog electrical signals are delivered to the transmitter 18. The transmitter 18 converts the analog electrical signals into digital data, encodes the data with error detection and correction information and multiplexes this data with control messages from the processor 22. The transmitter 18 modulates this combined data stream and broadcasts the resultant radio signals to the wireless network through the duplex filter 14 and the antenna 12. While the present invention is described relative to a cellular communication system, the invention is applicable to other wireless communication systems having communication channels subject to intersymbol interference (ISI). Digital cellular and personal communication systems based on IS- 136 and GSM (global system for mobile communications) require equalization to handle ISI arising from time dispersion. With the advent of the high date rates and the high level modulation in EDGE (enhanced data rates for global evolution) systems, the extent of ISI has increased considerably. To cover the same delay spread as today's GSM, EDGE receivers are specified to resolve a 5-tap channel. Furthermore, these systems are convolutional coded, as illustrated in Fig 4 To enable soft-decision decoding, it is, therefore, highly desirable for the equalizers to produce soft information at the bit level As illustrated in the system model of Fig 4, data is input to a convolutional encoder 40 The encoded data is output through a channel interleaver 42 to a modulator 44 The modulator 44 modulates the bits in a conventional manner to produce a signal x„ to be transmitted via a multipath fading channel represented at 46 A received sequence r„ is input to an equalizer 48 in accordance with the invention to produce the soft information The output of the equalizer 48 is delivered through a channel de-interleaver 50 to a convolutional decoder 52 For n = 1, 2, , N, the received sequence r n is related to the transmitted sequence x„ by where /.,, / = 0, 1 , , L, are complex-valued multipath fading coefficients and w π is complex- valued, zero-mean white Gaussian noise As discussed above, to improve performance, minimum-phase pre-filtering can be used to maximize the energy in the leading coefficients, i e., the energy in h, for small /. For MPSK modulation, where M = 2 q , the transmitted signal x„ is uniquely 2π — specified by the input <r/-tuple, (c nl , c„ 2 , ,c n ) e {0, 1 } q , via x n = m 5 = e κl , where S = c χ * 2 q~l + C 2 * 2 q~2 +...+C . A reverse <r/-ary to binary mapping is defined as -S^s) = c b , b = l, 2, , q For MDPSK modulation, where M = 2 q , the transmitted signal x n is given by where ./„ is an MPSK symbol uniquely specified by the input g-tuple, (c nl , c n2 , . . ., c nq ) e 2π J — {0,l } q , via x n = m s - e Ki , where S = C * 2 q~ + c 2 * 2 9~ + ...+ c . A reverse q-ary to binary mapping is defined as B b (s) = c b , b = 1, 2, . . ., q. The problem is to find, for all n and b, the α posteriori log-likelihood ratio of the bit c N given the complete received sequence γ = (r x , r 2 , . . ., r N ) and the α priori log-likelihood ratios Using the traditional log-MAP algorithm, referring to the notations described in Fig. 4, the forward state metric (5), backward state metric (6), and the branch metric (7, 8) are defined as ^ β ^ gPv {s n = j, r "} (5) /UJ logPr (r k = J} (6) = log Pr {S„ = j I S„_, = /} + log Pr {r n \ S n _ = i,S n = j] , (8) where γ represents the sequence (/ • „„ /'„, + ,, . . ., /'„). The forward and backward recursions are: 1,2,... , JV (9) !,... , ! (10) where the soft combining operation COM is defined as Efficient implementation of this operation is given using COM(x,y) = maxO, ) + log(l + 6T ,X_ I ) . (12) That is, the soft combiner first finds the maximum of the two inputs, as in the case for the Viterbi algorithm, and then adds a correction term Many efficient implementations of this operation are discussed, for example, by P Robertson, E Villebrun and P Hoehner, in "A Comparison of Optimal and Sub-Optimal MAP Decoding Algorithms Operating in the Log Domain," Proceedings of IEEE International Communications Conference '95, pp 1009- 1013, June 1995 One of the approach ignores the correction term completely to obtain a very computationally efficient algorithm, which is widely referred to as the Log-Max algorithm The invention disclosed herein can be used in conjunction with any of these techniques If the trellis begins at known state B and ends at known state E, the two recursions are initialized by • (I ) where ∞ represents a very large number permitted by the processor 22 or a special symbol recognizable by the processor 22 The a posteriori state log-likelihood and the a posteriori edge log-likelihood from state / to state / at time n can be then computed as i gP γ {s„ = j\ r "} = a„ω + fJ j (14) iogPr{ , = /A = j \ r = «„_,( )+ r„(/, ) + AO) • 05) The following describes the use of DDF-Log-MAP for MPSK in accordance with the invention Under the structure of DDFSE with 2 MLSE taps, there are M states at each stage of the DDFSE trellis Associated with each state / = 0, 1, , M - 1, at the (n - l)-th stage is a path history of /.-tuple MPSK symbols corresponding to hypothesized transmitted signals from stage n - Lion - 1 P t (n- 1)= ,(Λ- L),x t (n- L+ 1)...,JC.(«- 1) m. (16) For branch metrics, since all the paths that lead to the same end statey are driven by the same hypothesized symbol m p the first term of equation (8) is determined by logPrfø, = logPrjx,, = <7 = ∑B b (j) V(c nb )+K" b=\ where V(c„ έ ) is the a priori log-likelihood ratio for the bit c nb and B b (j) is the q-aτ to binary function defined above For this equalizer, the second term of equation (8) can be written as (18) where Re means taking the real part and rJt) is defined as which is the estimated received signal after decision-feedback cancellation. The two terms, / ,, h j X^ - /) and 2Re j ? * m * r n (/ ' ) \ , can be efficiently computed by, for example, techniques discussed in copending Application Serial No. entitled "Methods, Receivers and Equalizers Having Increased Computation Efficiency" to K. C. Zangi et al. and filed Feb. 8, 2000 (Attorney Docket No. 8194-354), which is incorporated by reference herein. For the forward recursion, since COM (a + x t ) = a + COM(x ι ) , the«-th stage of the forward recursion can then be re-arranged as a H U B>UW{c«). (20) For the decision feedback update, the new path history for statey is defined as P j (n) = xfn - L + \),... ,x s (u - 2),m s ,m j (21) where s is the index of the prior state from which the surviving branch is originated and: (22) That is, the surviving branch into state/ is defined as the branch ending at j with the maximum accumulated metric V.( - >;( 2 + 2 Re(/7»;( ]) * (23) Since finding the maximum is part of the COM operation, the last expression does not incur any new complexity. For the backward recursion, the («-l)-th stage can be re-arranged as: £_,(/) - (24) That is, the backward recursion does not hypothesize feedbacks and instead re-uses the modified received sequence P n (i) computed by the forward recursion. For the soft output, the a posteriori log-likelihood ratio of the bit c„ b is then given by Λ (c„ = COM (a n (s) + β n (s)) - CqM (a n s > ) + β n (s < )) . (25) That is, for each b - 1, 2,...,q, the a posteriori state log-likelihoods (a n (s) + β n (S)), , for 5 = 0, 1, ..., -1, are first separated into two sets such that one set has the b-th position in the binary representation of 5 being 1 and the other set has that position being 0. These two sets of log-likelihoods are then processed by the soft combiner separtely. Finally, the a posteriori log-likelihood for bit c.. > is obtained by taking the difference of the two combined results. The extrinsic information for the bit c nb can be obtained as W{c nb ) = K (c nb ) - V(c nb ) . (26) The above described method is formulated to execute the forward recursion first and then backward recursion in the above. It is straightforward to generalize it to execute the backward recursion first and followed by the forward recursion. Using proper implementations for the COM operations and techniques to compute d hιx,(n - /) and 2 Qγι 0 * m * r„ (ι) , the disclosd method is comparable to the traditional or SOVA-type DDFSE in total computational complexity However, in addition to those variables generally used in MLSE or DDFSE equalizers, the disclosed method stores the forward state metrics a n (j) if the forward recursion is executed before the backward recursion If the backward recursion is executed before the forward recursion, then the backward state metrics β „(ι) are stored For systems with small blocks, e g , D-AMPS, this storage requirement is not too excessive For systems with large blocks, the disclosed method can be applied in a sliding-window style as discussed by, for example, S Benedetto, G Montorsi, D Divsalar and F Pollara in "Soft-Output Decoding Algorithms in Iterative Decoding of Turbo Codes," The Telecommunications and Data Acquisition Progress Report, Jet Propulsion Laboratory, California Institute of Technology, vol 42-124, pp 63-87, Feb 1996, to reduce the required strorage size The following describes the use of DDF-Log-MAP for MDPSK in accordance with the invention Under the structure of DDFSE with 2 MLSE taps, there are M states at each stage of the DDFSE trellis Associated with each state / = 0, 1, , M - 1, at the (n - l)-th stage is a path history of I-tuple MPSK symbols corresponding to hypothesized transmitted signals from stage n - L to n - 1 P, (n - \) = [x, (n - L),x, (n - L + \) ...,x, (n - 2),x, (n - 1) = m . (27) For branch metrics, because of the differential encoding, the transition from state i to state / is driven by //„ = m δ , where δ e {0, 1, , M - 1 } is obtained by <? = (/ - J) Λ/ Δ (/ - J) mod . (28) where mod is the modulo operation. Therefore, the first term of equation (8) can be computed as where V(c„ A ) is the a priori log-likelihood ratio for the bit c nb and B b is the q- r to binary function defined above. Note that — j) jV(c nb ) can only assume one of the following M values: (30) M which can be pre-computed using only (q-2)— + \ real additions. Since all paths that lead to the same end state j hypothesize the same transmitted symbol m p for this equalizer, the second term of equation (8) can be computed as: logPr{r„ K_, = AT- |r„(/)| 2 + 2Re{/ι 0 * w;r w (#)} , where Re means taking the real part and V \l J is defined as L } (')= r n - d h l x l {n-l) (32) 1=1 which is the estimated received signal after decision-feedback cancellation. The two terms, / , _ h l x ι (n — I) and 2 Re (/ ' ) , can be efficiently computed by, for example, techniques discussed in the aforementioned reference incorporated herein. For the forward recursion, the /?-th stage can then be formulated as ) For the decision feedback update, the new path history for state; is defined as P.(n)= L+ \),...,x s (π-2),m s , m , (34) where s is the index of the prior state from which the surviving branch is originated and: That is, the surviving branch into state; is defined as the branch ending aty with the maximum accumulated metric (36) (α„-ι(0- K | 2 )+ 2 Re{/ 0 r n (/)}+ ∑ B b ((i - j) A )v(c„ b ) b- 1 Since finding the maximum is part of the COM operation, the last expression does not incur any new complexity. For the backward recursion, since COM(a + x t ) = a + COM(x ι ) ,the(»- l)-th stage of the backward recursion can be rearranged as: A-, ( = C {'f (37) That is, the backward recursion does not hypothesize feedbacks and instead re-uses the modified received sequence P n (i) computed by the forward recursion For the soft output, the a posteriori log-likelihood ratio of the bit c nb is then given by A (c nb ) = (*„_,(/) + Λ,0,./)+ AC/)) (38) - COM (*,,_,(/) + 7„ ,J)+ / J)) >> B i, ('-Λ, =o where y n (ι,j) as defined in equation (8) is simply the sum of equations (29) and (31) That is, theM" a posteriori edge log-likelihoods can be computed as ( nΛ (ι)+ γ„(ι,j) + β n (j)), for each 7 = 0, 1, , M- 1 anάj = 0, 1, , M- 1 For each b = 1, 2, , q, these log-likelihoods are first separated into two sets such that one set has the b-th position in the binary representation of ( — j )λi being 1 and the other set has that position being 0 These two sets of log- likelihoods are then processed by the soft combiner separately Finally, the α posteriori log- likelihood for bit Cnb is obtained by taking the difference of the two combined results It should be noted that the computational complexity of equation (38) can actually be reduced by first computing the α posteriori log-likelihoods of the driving symbol u„ and then deriving the bit log-likelihoods from these symbol log-likelihoods ForΛJ= 8, this two-step approach can reduce the complexity of equation (38) by 60% This method is demonstrated in the following First, theMdriving symbol log- likelihoods are obtained as r„(*) = COM (*,,-,(') + „(/,(/ + s) u ) + β,χd + δ) A/ ) (3 ) for δ = 0, 1 , , M - 1 The a posteriori log-likelihood for bit c nb can then be computed using Λ ( = CO (Γ„(^))- CO (Γ„(^))* (40) That is, for each b = 1, 2, , , these symbol log-likelihoods are first separated into two sets such that one set has the b-th position in the binary representation of δ being 1 and the other set has that position being 0 These two sets of log-likelihoods are then processed by the soft combiner separately Finally, the a posteriori log-likelihood for bit c„ b is obtained by taking the difference of the two combined results The extrinsic value for bit c nb can be obtained as Similar to the DDF-Log-MAP method for MPSK equalization discussed above, the disclosed method for MDPSK equalization can also be generalized to execute the backward recursion first and followed by the forward recursion Similar comments on complexity from that discussion apply to the MDPSK equalization method as well Referring to Fig 6, a flow diagram illustrates a program implemented in the processor 22 of Fig 3 for implementing the method of delayed decision feedback Log-MAP equalization in accordance with the invention As will be apparent, the equalization method is done in coordination with hardware in the receiver 16 of Fig 3 In describing the method, reference will be made to equations for both MPSK and, in parenthesis, MDPSK The program begins at a start node 60 when a signal is received The program has a first loop 62 performing the forward recursion for n=l, 2, N and, for each n, j=0, 1, ...M-l A second loop 64 performs the backward recursion for n=N, N-l, 1 and, for each n, i=0, 1, . M-l The first loop 62 begins at a block 66 that calculates and stores the received signal estimate using equation 1 (32) based on the path history from equation 16 (27) Then, at a block 68, the branch metric is calculated using equations 17 and 18 (29 and 31 ) Finally, at a block 70, the forward recursion is performed The forward recursion implements equation 20 (33) and then for each state j updates the path history of equation 21 (34) using equation 22 (35) A decision block 72 loops back to the block 66 until the loop 62 completes a forward trellis for all n and j and the estimated received signal values, forward metrics and path histories are stored in memory 24, see Fig 3, for each state and time Thereafter, the second loop 64 is performed to run a backward trellis The second loop 64 begins at a block 74 that performs the backward recursion The backward recursion implements equation 24 (37) for each state i Then, at a block 76, the soft output is calculated by determining the a posteriori log likelihood ratio using equation 25 (38 or 39 and 40) Finally, extrinsic information is obtained for each bit at a block 78 using equation 26 (41) A decision block 80 loops back to the block 74 until the loop 64 completes the backward trellis for all n and i The routine is then done As is apparent, to implement the backward recursion first the functions of blocks 70 and 74 would be reversed Thus, in accordance with the invention, an efficient method generates optimal soft information at the bit level for MPSK and MDPSK modulations over multipath fading channels The algorithms have a forward-backward recursion structure similar to conventional Log-MAP algorithms To reduce complexity, the forward recursion incorporates the techniques of DDFSE to construct a reduced state space The backward recursion then takes this state space as given, in that it does not hypothesize feedbacks, and computes required soft information

Claims: CLAIMS 1 A method of hypothesizing modulation symbols received over multipath fading channels, comprising calculating and storing an estimate of a received sequence of symbols using stored channel coefficients and a previous state path history of modulation symbols, for each possible current state symbol, calculating and storing a branch metric from a prior state to a current state using the calculated estimate of the received sequence of symbols, performing a forward recursion by calculating a forward state metric using the branch metric from the prior state to the current state, and performing a backward recursion by calculating a backward state metric using the branch metric from the prior state to the current state 2 The method of claim 1 wherein calculating and storing the estimate of a received sequence of symbols using stored channel coefficients and a previous state path history of modulation symbols is implemented using a delayed decision feedback sequence estimation equalizer 3 The method of claim 1 wherein calculating and storing the estimate of a received sequence of symbols using stored channel coefficients and a previous state path history of modulation symbols is implemented using a maximum likelihood sequence estimation equalizer for most recent symbols in the sequence of symbols and using a decision feedback equalizer for remaining symbols in the sequence of symbols 4 The method of claim 1 wherein calculating and storing the branch metric from a prior state to a current state using the calculated estimate of the received sequence of symbols, calculating the forward state metric using the branch metric from the prior state to the current state, and calculating the backward state metric using the branch metric from the prior state to the current state are implemented using a Log-MAP equalizer 5 The method of claim 1 further comprising calculating a posteriori log- likelihood ratio of each bit in the estimated modulation symbols using the calculated forward state metric and the calculated backward state metric 6 The method of claim 1 wherein performing the forward recursion by calculating the forward state metric using the branch metric from the prior state to the current state is executed prior to performing the backward recursion by calculating the backward state metric using the branch metric from the prior state to the current state 7 The method of claim 1 wherein performing the backward recursion by calculating the backward state metric using the branch metric from the prior state to the current state is executed prior to performing the forward recursion by calculating the forward state metric using the branch metric from the prior state to the current state 8 The method of claim 1 further comprising updating the previous state path history of modulation symbols using the calculated branch metric from a prior state to a current state using the calculated estimate of the received sequence of symbols 9 The method of claim 1 wherein the modulation symbols received over multipath fading channels comprise M-ary modulated symbols 10 The method of claim 1 wherein the modulation symbols received over multipath fading channels comprise M-ary differential PSK modulated symbols 11 In a receiver receiving modulation symbols transmitted over multipath fading channels, an equalizer comprising a controller programmed to calculate and store an estimate of a received sequence of symbols using stored channel coefficients and a previous state path history of modulation symbols, for each possible current state symbol calculate and store a branch metric from a prior state to a current state using the calculated estimate of the received sequence of symbols; perform a forward recursion by calculating a forward state metric using the branch metric from the prior state to the current state, and perform a backward recursion by calculating a backward state metric using the branch metric from the prior state to the current state 12 The equalizer of claim 1 1 wherein the controller calculates the estimate of a received sequence of symbols using stored channel coefficients and a previous state path history of modulation symbols using a delayed decision feedback sequence estimation 13 The equalizer of claim 1 1 wherein the controller calculates the estimate of a received sequence of symbols using stored channel coefficients and a previous state path history of modulation symbols using maximum likelihood sequence estimation for most recent symbols in the sequence of symbols and using a decision feedback equalization for remaining symbols in the sequence of symbols 14 The equalizer of claim 1 1 wherein the controller calculates the branch metric from a prior state to a current state using the calculated estimate of the received sequence of symbols, calculates the forward state metric using the branch metric from the prior state to the current state, and calculates the backward state metric using the branch metric from the prior state to the current state using Log-MAP equalization 15 The equalizer of claim 1 1 wherein the controller calculates a posteriori log-likelihood ratio of each bit in the estimated modulation symbols using the calculated forward state metric and the calculated backward state metric 16 The equalizer of claim 1 1 wherein the controller performs the forward recursion by calculating the forward state metric using the branch metric from the prior state to the current state prior to performing the backward recursion by calculating the backward state metric using the branch metric from the prior state to the current state 17 The equalizer of claim 1 1 wherein the controller performs the backward recursion by calculating the backward state metric using the branch metric from the prior state to the current state prior to performing the forward recursion by calculating the forward state metric using the branch metric from the prior state to the current state 18 The equalizer of claim 1 1 wherein the controller updates the previous state path history of modulation symbols using the calculated branch metric from a prior state to a current state using the calculated estimate of the received sequence of symbols 19 The equalizer of claim 1 1 wherein the modulation symbols received over multipath fading channels comprise M-ary modulated symbols. 20 The equalizer of claim 1 1 wherein the modulation symbols received over multipath fading channels comprise M-ary differential PSK modulated symbols 21 A method of generating soft information at bit level of modulation symbols received over multipath fading channels, comprising calculating and storing an estimate of a received sequence of symbols using stored channel coefficients and a previous state path history of modulation symbols, each symbol comprising an M-ary modulated symbol representing select bits c, for each possible current state symbol, calculating and storing a branch metric from a prior state to a current state using the calculated estimate of the received sequence of symbols; performing a forward recursion by calculating a forward state metric using the branch metric from the prior state to the current state, performing a backward recursion by calculating a backward state metric using the branch metric from the prior state to the current state, and calculating a posteriori log-likelihood ratio of each select bit c in the estimated received modulation symbols using the calculated forward state metric and the calculated backward state metric 22 The method of claim 21 wherein calculating and storing the estimate of a received sequence of symbols using stored channel coefficients and a previous state path history of modulation symbols is implemented using a delayed decision feedback sequence estimation equalizer 23 The method of claim 21 wherein calculating and storing the estimate of a received sequence of symbols using stored channel coefficients and a previous state path history of modulation symbols is implemented using a maximum likelihood sequence estimation equalizer for most recent symbols in the sequence of symbols and using a decision feedback equalizer for remaining symbols in the sequence of symbols 24 The method of claim 21 wherein calculating and storing the branch metric from a prior state to a current state using the calculated estimate of the received sequence of symbols, calculating the forward state metric using the branch metric from the prior state to the current state, and calculating the backward state metric using the branch metric from the prior state to the current state are implemented using a Log-MAP equalizer 25 The method of claim 21 wherein performing the forward recursion by calculating the forward state metric using the branch metric from the prior state to the current state is executed prior to performing the backward recursion by calculating the backward state metric using the branch metric from the prior state to the current state 26 The method of claim 21 wherein performing the backward recursion by calculating the backward state metric using the branch metric from the prior state to the current state is executed prior to performing the forward recursion by calculating the forward state metric using the branch metric from the prior state to the current state 27 The method of claim 21 further comprising updating the previous state path history of modulation symbols using the calculated branch metric from a prior state to a current state using the calculated estimate of the received sequence of symbols 28 The method of claim 21 wherein the modulation symbols received over multipath fading channels comprise M-ary differential PSK modulated symbols 29 A mobile terminal receiving modulation symbols transmitted over multipath fading channels, comprising a receiver for receiving modulation symbols transmitted over multipath fading channels and including an equalizer calculating and storing an estimate of a received sequence of symbols using stored channel coefficients and a previous state path history of modulation symbols, calculating and storing a branch metric from a prior state to a current state using the calculated estimate of the received sequence of symbols for each possible current state symbol, performing a forward recursion by calculating a forward state metric using the branch metric from the prior state to the current state, and performing a backward recursion by calculating a backward state metric using the branch metric from the prior state to the current state 30 The mobile terminal of claim 29 wherein the equalizer calculates the estimate of a received sequence of symbols using stored channel coefficients and a previous state path history of modulation symbols using a delayed decision feedback sequence estimation 31 The mobile terminal of claim 29 wherein the equalizer calculates the estimate of a received sequence of symbols using stored channel coefficients and a previous state path history of modulation symbols using maximum likelihood sequence estimation for most recent symbols in the sequence of symbols and using a decision feedback equalization for remaining symbols in the sequence of symbols 32 The mobile terminal of claim 29 wherein the equalizer calculates the branch metric from a prior state to a current state using the calculated estimate of the received sequence of symbols, calculates the forward state metric using the branch metric from the prior state to the current state, and calculates the backward state metric using the branch metric from the prior state to the current state using Log-MAP equalization 33 The mobile terminal of claim 29 wherein the controller calculates a posteriori log-likelihood ratio of each bit in the hypothesized modulation symbols using the calculated forward state metric and the calculated backward state metric 34 A method of generating soft information at bit level of modulation symbols received over multipath fading channels, comprising performing a forward recursion by calculating and storing an estimate of received symbols using decision feedback cancellation, each symbol comprising an M-ary modulated symbol representing select bits c, calculating and storing branch metrics using the calculated estimate of the received symbols, and calculating forward state metrics using the branch metrics and the estimate of received symbols and, performing a backward recursion by calculating backward state metrics using the branch metrics and the estimate of received symbols stored from the forward recursion, and calculating a posteriori log-likelihood ratio of each select bit c in the estimate of received modulation symbols using the calculated forward state metric and the calculated backward state metric

Inventors: CHENG JUNG-FU, CHENG, JUNG-FU
Applicants: ERICSSON INC, ERICSSON INC.
Application year: 2001
Main classification: H04L25
IPC codes: H04B3, H04L25, H04L27