APPARATUS, METHOD AND COMPUTER PROGRAM PRODUCT PROVIDING
AMIMO RECEIVER
TECHNICALFffiLD:
The exemplary and non-limiting embodiments of this invention relate generally to wireless communications systems, methods, devices and computer programs and, more specifically, relate to multiple input, multiple output (MlMO) wireless communications systems.
BACKGROUND:
The following abbreviations that may appear in the description are defined as follows:
AWGN additive white Gaussian noise
APP a posterior probability
MIMO multiple-input multiple-output
MISO multiple-input single-output
SISO single-input single-output BP belief propagation
SNR signal-to-noise ratio
SAP serial-to-parallel
QF quasistatic fading
IF independent fading FIR finite impulse response
E-UTRAN evolved universal terrestrial radio access network
OFDM orthogonal frequency division multiplexing
WCDMA wideband code division multiple access
BICM bit interleaved coded modulation
CM coded modulation LDPC low density parity check
ML maximum likelihood r.v. random variable
QAM quadrature amplitude modulation
QPSK quadrature phase shift keying
The use of an efficient detection method is important in order to accommodate multiple antenna transmissions and/or high order constellations. In such cases, detection by the use of an exhaustive search is prohibitive due to the large number of valid signal combinations across multiple transmit antennas (Cartesian product of individual antenna constellations). A reduced search method, that can nonetheless perform very close to optimal (maximum likelihood or ML) detection, will be required for future generations of wireless systems (cellular and non-cellular). In addition, it would be desirable to have available soft information from the detector in order to improve the performance of a decoder, or to enable iterations between detection and decoding. Further, modularization is desirable, as is any simplification over traditionaralgorithms, such as sphere decoders.
Most good signal constellations exhibit a lattice structure, which can be used to greatly simplify the search. Sphere detectors (as well as decoders) have been developed primarily in response to the need to alleviate the complexity of ML estimation- for a large number of hypotheses. Since the underlying problem is a search, any complexity reduction would need to implement a reduced search procedure. Traditionally, the result of a reduced search algorithm is a hard decision. Upon recognizing the importance of soft information at the detector's output, some researchers began to explore sphere detecting algorithms that are capable of providing a soft information output.
A hard decision sphere detector algorithm was based on an algorithm devised by Pohst
(M. Pohst, "On the computation of lattice vectors of minimal length, successive minima, and reduced basis with applications", ACM SIGSAM Bull., vol. 15,pp.37-44, 1981, and
U. Finclce and M. Pohst, "Improved methods for calculating vectors of short length in a lattice, including a complexity analysis", Math. Comput, vol. 44, pp. 463-471, Apr. 1985), and were described by Viterbo and Boutros (E. Viterbo and J. Boutros, "A universal lattice decoder for fading channels", IEEE Trans. Inform. Theory, vol. 45, No. 5, pp. 1639-1642, July 1999). An improved algorithm for closest point search in a lattice was proposed by Schnorr and Euchnerr (C. P. Schnorr and M. Euchnerr, "Lattice basis reduction: improved practical algorithms and solving subset sum problems," Math. Programming, vol. 66, pp. 181—191, Apr. 1994), which starts at the center of the valid range of lattice points, and has better efficiency. Agrell et al. (E. Agrell, T. Eriksson, A. Vardy, and K. Zeger, "Closest point search in lattices", IEEE Trans. Inform. Theory, vol. 48, No. 2, pp. 2201-2214, Aug.2002) devised another algorithm that shows a slight gain at a low signal to noise ratio (SNR).
Seeking to obtain a soft information output, Boutros et al. proposed an approach to soft output sphere detecting (J. Boutros, N. Gresset, L. Brunei, and M. Fossorier, "Soft-input soft-output lattice sphere decoder for linear channels", Proc. IEEE Conf. Globecom'03, pp. 1583-1587, 2003) without resorting to basis conversions (boundaries of search regions are difficult to determine); instead they take advantage of the finite structure of the constellation (finite modulation alphabet).
Other soft information approaches allow for complex non-lattice modulation alphabets, or replace the sphere with a different body, or implement a list (see, for example, S . Baro,
J. Hagenauer, and M. Witzke, "Iterative detection of MIMO transmission using a list-sequential (LISS) detector", IEEE International Conf. Commun., ICC'03, vol. 4, pp.
2653-2657, Anchorage, May 2003; Y. de Jong and T. Willink, "Iterative tree search detection for MIMO wireless systems", IEEE 56th Vehicular Technology Conf. VTCO2 Proceedings, vol. 2, pp. 1041-1045, Sept. 2002; and B. M. Hochwald and S. ten Brink,
"Achieving near-capacity on a multiple-antenna channel", IEEE Trans. Commun., vol.
51, pp. 389-399, March 2003.)
Regarding the soft-input soft-output lattice sphere detector from J. Boutros et al., it can be noted that it requires three reduced search passes: first a (Schnorr-Euchnerr) reduced sphere detector search is needed to obtain the initial hard decision ML point, a second
search pass then enumerates all lattice points in a sphere centered on the newly found ML point, and a third search pass then evaluates the squared distances needed for generating the soft information. This triple-pass strategy to computing a soft-information detector output increases the detector's complexity, especially since in the final pass the Euclidean distances are computed with respect to both the ML point (from the hard decision pass) and the received point.
Current trends in modern and proposed wireless communication systems, such as one known as Evolved Universal Terrestrial Radio Access Network (E-UTRAN), aim at achieving high data rates at relatively low costs, and mandate multi-carrier designs, high spectral efficiencies and MIMO technology.
SUMMARY
The foregoing and other problems are overcome, and other advantages are realized, through the use of the exemplary embodiments of this invention.
hi a first aspect thereof the exemplary embodiments of this invention provide a method that includes receiving a plurality of signals through a plurality of antennas, the plurality of signals being modulated with a space-time lattice code; removing an effect of a channel matrix from the received signals to provide an equalized received signal; and lattice detecting the equalized received signal based on a Tanner graph representation of the lattice.
In another aspect thereof the exemplary embodiments of this invention provide a computer program product that is embodied in a computer readable medium and that includes instructions, the execution of which result in performing operations that comprise: in response to receiving a plurality of signals through a plurality of antennas, the plurality of signals being modulated with a space-time lattice code, removing an effect of a channel matrix from the received signals to provide an equalized received signal; and lattice detecting the equalized received signal based on a Tanner graph representation of the lattice.
In a still further aspect thereof the exemplary embodiments of this invention provide an apparatus that includes an equalizer configured to respond to a plurality of signals received through a plurality of receive antennas to remove an effect of a channel matrix from the received signals to provide an equalized received signal, the plurality of signals being transmitted from a plurality of transmit antennas modulated with a space-time lattice code. The apparatus further includes a detector configured to operate on the equalized received signal in accordance with a Tanner graph representation of the lattice to perform lattice detection and to output soft information concerning real coordinates of complex symbols from modulation constellations used at the plurality of transmit antennas.
In another aspect thereof the exemplary embodiments of this invention provide an integrated circuit that includes an equalizer circuit configured to respond to a plurality of signals received through a plurality of receive antennas to remove an effect of a channel matrix from the received signals to provide an equalized received signal, the plurality of signals being transmitted from a plurality of transmit antennas modulated with a space- time lattice code; and a detector circuit configured to operate on the equalized received signal in accordance with a Tanner graph representation of the lattice to perform lattice detection and to output soft information concerning real coordinates of complex symbols from modulation constellations used at the plurality of transmit antennas.
In a further aspect thereof the exemplary embodiments of this invention provide an apparatus that includes means for equalizing a plurality of signals received through a plurality of receive antennas to remove an effect of a channel matrix from the received signals to provide an equalized received signal, the plurality of signals being transmitted from a plurality of transmit antennas modulated with a space-time lattice code; and means for operating on the equalized received signal in accordance with a Tanner graph representation of the lattice to perform lattice detection and to output soft information concerning real coordinates of complex symbols from modulation constellations used at the plurality of transmit antennas.
BRIEF DESCRIPTION OF THE DRAWINGS
In the attached Drawing Figures:
Figure 1 is an example of a Tanner graph.
Figure 2 illustrates a projection of a point.
Figure 3 shows a state transition diagram for a Markov process representing a sequence of lattice points.
Figure 4 is a block diagram of an iterative receiver for a super-orthogonal space-time lattice code in the presence of a coordinate interleaver, in accordance with exemplary embodiments of this invention.
Figure 5 is a graph that plots FER versus Eb/N0 for a super-orthogonal space-time lattice code, with MMSE followed by BP.
Figure 6 is a graph that plots FER versus Eb/No of iterative decoding based on IC-MMSE plus BP for a super-orthogonal space-time lattice code with coordinate interleaver.
Figure 7 shows a simplified block diagram of one non-limiting embodiment of a MIMO system that is suitable for use in practicing the exemplary embodiments of this invention.
Figure 8 is a logic flow diagram that is illustrative of a method and/or execution of a computer program product, in accordance with exemplary embodiments of this invention.
DETAILED DESCRIPTION
The exemplary embodiments of this invention relate in general directly or indirectly to transmit antenna diversity, MIMO systems, lattice constellations, lattice detection and decoding, soft information, sphere decoding, iterative receivers, belief propagation,
Tanner graphs, multipath channels, closed loop schemes, channel estimation, orthogonal
frequency division multiplexing (OFDM), space-time coding, spatial precoding, spatial redundancy, beamf orming, transmission parameter adaptation and multi-carrier systems, as non-limiting examples.
In accordance with the exemplary embodiments of this invention there is provided an efficient detection method that accommodates multiple antenna transmissions with high order signal constellations.
Further in accordance with the exemplary embodiments of this invention there is provided a reduced search method that is capable of performing in a manner that closely approaches optimal (ML) detection, and that is suitable for use in implementing future generations of wireless systems (both cellular and non-cellular).
Further hi accordance with the exemplary embodiments of this invention there is provided an ability to output soft information from the detector in order to improve decoder performance and/or to enable iterations between detection and decoding.
Further in accordance with the exemplary embodiments of this invention there is provided a modularization capability. Related to the foregoing, exemplary embodiments of this invention employ belief propagation, and this functionality may be exploited in receiver implementations that include some form of, by example, a Low Density Parity Check (LDPC) decoder, wherein the architecture may be designed so that the belief propagation module is reusable.
Still further in accordance with the exemplary embodiments of this invention there is provided a technique to avoid a step back artifact that is a result of operation of conventional sphere decoders.
Figure 7 is a block diagram of an exemplary MIMO system 10 that is suitable for practicing this invention. The MDVIO system 10 includes a transmitter 12 and at least one receiver 14. The transmitter 12 has a plurality of transmit (T) antennas (T1-TJyn.) and associated transmit amplifiers 12A5 and atransmit control function 12B. The receiver 14 has one or more receive (R) antennas (RΓRMR) an<^ associated receive amplifiers 14A,
and a receive control function 14B. In general, the number of transmit antennas may or may not equal the number of receive antennas, and both are preferably greater than one. The transmit control function 12B is assumed to include one or more sources of data, as well as an encoder and modulator, and any other circuitry needed to transmit data, such as packet data (control and/or traffic data packets), to the receiver 14. The receive control function 14B is assumed to include one or more data sinks, as well as a complementary data decoder and demodulator, and any other circuitry needed to receive data, such as packet data, from the transmitter 12.
The transmit control function 12B may include at least one data processor (DP) 12C that is operable to execute program code in order to operate as a MIMO transmitter. The receive control function 14B may include at least one data processor (DP) 14C that is operable to execute program code in order to operate as a MIMO receiver, in particular one that operates in accordance with the exemplary embodiments of this invention. Further in this regard the receiver 14 implements a novel iterative receiver, as shown in block diagram form as having, for example, an Inner Decoder (IC-MMSE) 15 A, an Outer Decoder 15B and a Soft Estimator 15 C, as is also shown in Figure 4 and described in detail below. The DPs 12C, 14C may be embodied as one or more digital signal processor (DSP) and/or other integrated circuits, or in any form that is suitable for implementing the exemplary embodiments of this invention.
In general, the exemplary embodiments of this invention may be implemented by computer software executable by at least the DP 14C, or by hardware, or by a combination of software and hardware, and also firmware.
Embodiments of the receiver 14 may be realized in, but are not limited to, cellular phones, personal digital assistants (PDAs) having wireless communication capabilities, portable computers having wireless communication capabilities, image capture devices such as digital cameras having wireless communication capabilities, gaming devices having wireless communication capabilities, music storage and playback appliances having wireless communication capabilities, Internet appliances permitting wireless Internet access and browsing, as well as portable units or terminals that incorporate combinations of such functions.
Note that at least one or more of the Inner Decoder (IC-MMSE) 15 A, the Outer Decoder 15B and the Soft Estimator 15C may be embodied in one or more integrated circuits.
The exemplary embodiments of this invention beneficially employ an efficient, low-complexity, soft-information detector for MIMO channels and lattice constellations, and are based on, as one example, Tanner graph representations of lattices. Due to the coding gain associated with a lattice, structural relations exist between certain lattice points, which can be associated via an equivalence relation for detection purposes. The detector algorithm in accordance with the exemplary embodiments of this invention is capable of generating both total and extrinsic a posteriori probability at the output of the detector. The step-back artifact (characteristic of traditional sphere decoders) is eliminated. The algorithm applies to general lattices and enables the provisioning of the iterative type of receiver 14.
As non-limiting and exemplary implementations the operation of a novel algorithm was simulated for an uncoded transmission for a super-orthogonal constellation in two scenarios, hi a quasi-static (block) fading scenario it was found to achieve ML performance even with one 'surviving' label (out of six). In an independent fading scenario, with coordinate (component) interleaving and iterations between equalization and detection, the algorithm was found to perform close to interference-free transmission. The coordinate interleaved scenario was furthermore found to outperform the former scenario, despite the absence of forward error correction coding, hi that the presence of one of six labels have been shown to be sufficient in the exemplary implementation, it can be expected that complexity may be reduced to about approximately 17% to 20% of that of the exhaustive (optimal) search, thereby enabling very efficient implementations to be provided.
A more detailed description of the foregoing exemplary aspects of this invention, including the simulations referred to above, is found below.
SUMMARY OF SPHERE DECODING RESEARCH
Sphere detectors (as Well as decoders) arose primarily from the need to alleviate the complexity of ML estimation for a large number of hypotheses. Since the underlying problem is a search, complexity reduction must come from a reduced search. Traditionally, the result of the reduced search algorithm was a hard decision; later, upon recognizing the importance of soft information at the detector's output, researchers began looking at sphere detecting algorithms capable of providing a soft information output.
The hard decision sphere detector algorithm was based on an algorithm devised by Pohst [1], [2], and described by Viterbo and Boutros [4] (there was an earlier paper by Viterbo and Biglieri, in 1993). An improved algorithm for closest point search in a lattice was proposed by Schnorr and Euchnerr [3], which starts at the center of the valid range of lattice points, and has better efficiency. Agrell et al. [31] devised another algorithm that shows a slight gain al low SNR.
Seeking to obtain a soft information output, Boutros et al. proposed a clean, elegant approach to soft output sphere detecting [5] without resorting to basis conversions (boundaries of search regions are difficult to determine); instead they take advantage of the finite structure of the constellation (finite modulation alphabet).
Other soft information approaches allow for complex non-lattice modulation alphabets, or replace the sphere with a different body, or implement a list [7], [6]. ^
Regarding the soft-input soft-output lattice sphere detector from [5], we note that it requires three reduced search passes — first a (Schnorr-Euchnerr) reduced sphere detector search needed to obtain the initial hard decision ML point, second to enumerate all lattice points in a sphere centered on the newly found ML point, and third to evaluate the squared distances needed for generating the soft information. This triple-pass strategy to computing a soft-information detector output increases the detector's complexity — especially since in the last pass Euclidean distances are computed with respect to both the ML point (from the hard decision pass) and the received point.
I. INTRODUCTION
Multiple input multiple output (MIMO) transmission has emerged as a strong scenario for future high-speed wireless communications due to the large capacity potential of MIMO channels. Space-time codes that exploit both spatial diversity and time diversity have been widely proposed as MIMO modulation in the past decade to achieve reliable transmission.
Recently, the importance of lattice MIMO constellations in constructing space-time lattice codes was recognized by El-Gamal et al. [26] from a diversity-multiplexing tradeoff perspective. Superorthogonal space-time codes — first reported in [18], then in [19], [20], [21], [22] (where they were dubbed 'superorthogonal') — are in fact lattice space time codes (see [23, Section III] and Example 2). As lattices, such constellations lend themselves to efficient detection algorithms, e.g. sphere decoding. Classic sphere decoding (see [311, and references therein) use hard decision, and a step-back provision; soft-output versions have been imagined, but rely on a list of important candidates, and retain the step-back provision. In [17], lattice partitioning is used to divide the lattice into a finite number of cosets. Each coset is then labeled by a codeword of a finite Abelian group block code. In [29], a Tanner graph (TG) representation for the label code was developed; this opens an opportunity for using belief propagation on the lattice labels.
The sequel takes a qualitatively different approach to soft output closest point search in lattices, via a form of belief propagation on a lattice. Due to the coding gain associated with a lattice, structural relations exist between certain lattice points, which can be associated via an equivalence relation for detection purposes. The algorithm can generate both total and extrinsic a posteriori probability (APP) at the detector's output. The step-back feature is eliminated. For each channel use, a filter bank for interference cancellation with minimum mean square error (IC-MMSE) is used to remove the channel effects. Then, a reduced-complexity lattice decoder based on TG lattice representation is proposed for computing total APP and extrinsic APP. The capability of calculating the extrinsic APP enables decoding schemes that iterate between detection and decoding. This novel lattice detection algorithm is applied to detecting superorthogonal space-time lattice codes [23] in quasistatic fading, and to a coordinate interleaved [34] scenario. The following notation is adhered to. Vectors are denoted by lowercase bold letters; βi denotes the i-th element of vector α. Matrices are denoted by uppercase bold letters. The i-th column vector and the y-th element of a matrix, say A, are denoted by Ci1 and a^, respectively. The superscripts T and H are used to denote transposition and complex conjugated transposition, respectively.
II. PROBLEM DEFINITION AND SYSTEM MODEL
Complex and real transmission models are described; a general formulation for lattice constellations for MTMO channels is then introduced, followed by two examples pertaining respectively to linear dispersion and superorthogonal codes.
A. Rayleigh flat fading MIMO channels
Consider MIMO wireless transmission with Nt transmit antennas and JV7. receive antennas in Rayleigh flat fading. The channel coefficients are assumed to be constant over a block of T MIMO channel uses and change independently from block to block. The transmission of each block is then given by
Y = ^IfN
tSH + N (1)
If detection and decoding can be performed iteratively, then soft information about x can be fed back from the FEC decoder and made available to the filter bank in the. form of probabilities of valid realizations of transmitted vectors x, or its elements XiI i.e. either at the vector level x, {Pr(a; = φ(C
τ)) φ(C
τ S C(A, UQ, TV) }, or at the coordinate level — e.g. in the case when coordinate interleaving [34] is used to scramble the coordinates of several vectors x prior to trasmission: In the latter case the structure present in the different multidimensional lattice points is destroyed during transmission through the channel; not only does this mean that the coordinate probabilities supplied by the decoder have to be unscrambled before being fed back to the LMMSE filter for interference cancellation (IC — see Fig. 4), but the performance can be improved (over the non-interleaved scenario) even in an uncoded system (see Section IV-B).
An iterative receiver aims at iteratively canceling the interference prior to filtering by forming as soft interference estimator in one of two ways:
1) Vector level feedback.
2) Coordinate level feedback: If Ki is the ith. coordinate alphabet, the average interference value at position i is
Let X
1Q f denote the vector obtained by setting the i-th element of asic to zero, i.e., x
loj = [
• • - , xic,j-i, 0, -cic,
t+i,
• • -]
τ, the interference cancellation is performed for the i-th branch
and the soft estimate $n of the i-th. branch after IC is
subject to a unit power constraint like (24). The estimation (30) is referred to IC-MMSE. The covariance matrix of y
i5 denoted as J2IG,», is
■wi
- Substituting JRi
G,» of (31) for R in (25), (26), yields the IC-MMSE solution mi and the corresponding MSE σf, respectively. Note that the IC-MMSE filter bank is a more general solution than a MMSE filter bank for removing channel effects in a MIMO scenarios. After IC-MMSE filtering the soft estimate of the ith branch is
with ήi ~ Λ/"(0, σf ), or written in a matrix form as
B. Belief propagation detector for lattice code based on Tanner graph representation
After IC-MMSE equalization, the soft estimate x of a lattice point is obtained. Recall that in lattice space-time schemes, the codebook of transmitted vectors a; is a lattice code C(A, uo, TV), where the generator matrix of Λ is TB. For simplicity, bet B be a generic lattice generator matrix. Lattice detection is to either decide which lattice point inside the shaping region has the minimum distance to &, or calculate the soft information (e.g., in the form of probability or log-likelihood ratio) about each candidate lattice point. The first detection criterion leads to hard decision detectors — e.g., maximum likelihood (ML). The second decoding criterion leads to soft decision detectors, which can be used in iterations between detection and decoding. In this section, a novel Tanner graph based lattice decoding algorithm is introduced. For simplicity, assume an m-dimensional lattice code, i.e., x e Rm.
The novel lattice decoding algorithm introduced below relies on Tanner graph representations of lattices [29], which are enabled by lattice partitioning; all lattice points (those inside the shaping region are of interest) are partitioned into several subgroups (cosets). Each subgroup includes several different lattice points, and is labelled by a well-defined Abelian group block codeword. Then, a reduced-complexity soft-output lattice detector can be obtained by operating on the smaller number of cosets instead of lattice points. The labels of all cosets form an Abelian block code, which can be represented by a Tanner graph similar to low-density-parity-check (LDPC) codes. Belief propagation on a lattice is performed on its non-binary label Tanner graph to yield the total and extrinsic APP of the labels and their coordinates, as described in the following subsections. The APPs of individual lattice points are obtained in a final step described in Section III-D.
A somewhat subtler point is that lattice partitioning revolves around an orthogonal sublattice Λ' of Λ, and the quotient group A/Λ'; |Λ/A'| is finite iff A and Λ' have the same dimensionality. The most straightforward way of obtaining Λ' is by G-S orthogonalization of Λ's generator matrix, whereby all orthogonal G-S directions intercept A and the intersection naturally forms a sublattice of the same dimensionality as Λ; in all other cases the orthogonal sublattice will have to be obtained by some means other than G-S orthogonalization.
where ABS(α) = ∑ |αj|. Consequently, the probability of hypotheses Hx, H
2 catrbe obtained from L(H)
For each hypothesis one can apply the lattice detection algorithm developed in Section IH for detecting χ. We treat the information-bearing vector χ as a lattice with generator matrix B, i.e.,χ = Bu. For example, the equivalent model for detecting lattice point x is χ = Bu + ή1, where % is the output of matched filtering of hypothesis Hx. Since χ is from a X>4 lattice, its generator matrix B is given in (35). The APPs can be obtained according to Section EI.
B. Iterative receiver for coordinate interleaving in fast fading
Coordinate interleaving, along with the outer iteration loop in Fig. 4, is now considered; the real and imaginary parts of all complex symbols in a frame are collectivelly scrambled before transmission [34]. Y — {y
t, y
2, ■ ■ • , y
N} denotes a frame spanning JV MIMO channel uses at the MMO channel output (before deinterleaving). Note that the structure of the superorthogonal lattice code is removed during transmission, and has to be recovered before detection. The applicable receive equation is (6) rather than (19); the iterative IC-MMSE attempts to iteratively remove the cross-antenna interference, Le. to undo the channel H on a per MIMO channel use basis. During the first iteration, the soft feedback from the detector/decoder is null. The output of IC-MMSE is always deinterleaved, thus restoring the superorthogonal structure and yielding the soft output X = -(SC
1, X
2, ■
■ ■ , £Djv} with
Since the information-bearing vector χ
θ t is a direct sum of two D
4 lattices, and the effective channel gain matrix T is unitary, the equalization approach in Section IV-A applies to eq. (62).
k = 1, 2, are associated with the following transmission models upon removing F
1, T
2 respectively:
wher
e generator matrix B is given in (35). For each hypothesis, the lattice decoding algorithm can be applied to compute the extrinsic APPs P(u; O) and P(c; O).
Inner-loop iterative decoding between SISO and BP, as shown in Fig. 4, can further improve the overall performance, especially in the presence of forward error correction coding, when decoding follows detection. Herein, only an uncoded system is considered in order to illustrate the concept. Even in an uncoded system it is possible to perform inner loop iterations between PBP(c; O) from the belief propagation module and P(u;I) from the SISO block; more benefit is derived however when a decoder is part of the inner-loop.
V. SIMULATIONS
Simulation results for a superorthogonal space-time lattice code with 4PSK constellation (Example 2), in both quasistatic and fast fading channels, are discussed. Each half of the superorthogonal constellation belongs to a D4 lattice, implicitly defining a shaping region; only six of the twelve L(A.) labels listed in Example 3 (first four, last two) are needed to cover the lattice points in the shaping region. In order to test the algorithm's efficiency, only the most likely label (or two labels) — post belief propagation — are retained; the others receive zero probabilities (re-normalization is performed after setting to zero the probabilities of discarded labels).
A. Quasistatic fading
The channel is constant over T = 2 symbol periods. In our simulations, each data packet includes 500 super-orthogonal codewords. Each point on the curves plotted in Fig. 5 and Fig. 6 is obtained by testing 2000 independent data packets.
Fig. 5 shows the FER (frame error ratio)8 vs. Eb/ NQ for super-orthogonal space-time code when the coordinate interleaver is absent. QPSK modulation is employed and the channel spectral efficiency is 2.5 bits/channel use. The performance of the ML algorithm that exhaustively searches all possible valid codewords and picks the one with the ML is plotted as reference. For the MMSE-BP algorithm, we run one iteration for the Tanner graph and collect the probability of the coordinate of label. Then we consider choosing one surviving label and two surviving labels. The simulation result shows that the MMSE-BP algorithm with one surviving label and two surviving labels have the same performance as that of the ML algorithm. The MMSE-BP with simplified initialization that reduces the overall complexity is also examined. In this case, we consider two surviving labels, the results show that it is about 0.5 dB away from the ML performance in low SNR region. As SNR increases, the MMSE-BP with simplified initialization approaches the ML performance asymptotically.
B. Fast fading
Fast fading simulations include a coordinate interleaver. In our simulations, a depth-eight traditional block interleaver is considered. QPSK is used and the channel spectral efficiency is 2.5 bits/channel use. Two inner iterations are. run between the SISO block and the BP block; one iteration is run on the lattice Tanner graph inside the BP block. We simulate different scenarios where different number of surviving labels are considered, hi addition, iterative interference cancellation scheme is considered to improve the overall performance. The soft estimator computes the soft estimates of the coordinates of lattice point based on the output from the BP ( P(u; O)). Fig. 6 shows the FER vs.Et/No for different number of surviving labels and different number of iterations between the IC-MMSE and the outer decoder.
VI. CONCLUSION
A soft output closest point search in lattices was introduced, via a form of belief propagation on a lattice. Due to the coding gain associated with a lattice, structural relations exist between certain lattice points, which can be associated via an equivalence relation for detection purposes. This leads to a soft output detection algorithm, which can generate both total and extrinsic a posteriori probability at the detector's output. The step-back feature of classic sphere decoding is eliminated.
8One frame is meant to be one super-orthogonal space-time codeword
COMPUTING THE EXTRINSIC A POSTERIORI PROBABILITIES AFTER BELIEF PROPAGATION
Herein, the expressions for extrinsic a posteriori probabilities (46), (47), at the belief propagation detector's output, are derived; the extrinsic probabilities are needed in iterative receivers. Here, the goal of detection is to provide soft information about valid channel alphabet symbols, i.e. real coordinates of the complex symbols from the modulation constellations used on various transmit antennas; this information about coordinates can be used to revert the effect of a coordinate interleaver, or can be forwarded directly to a soft decoder for some coded modulation encoder. Alternatively, it can be used for soft or hard demodulation, e.g. in the' case of bit interleaved coded modulation, or with plain uncoded transmission.
When a lattice is represented by a Tanner graph, it is possible to associate a Markov process with the model for soft detection of lattice points in a natural way. This is enabled by first viewing the sequence of lattice points passed through the channel as a Markov source. Another observation is that, 'in general, simple detection (with or without soft information) is by itself memoryless; thereby, one should expect the Markov process to be somehow degenerated, in order to reflect the memoryless nature of simple (non-iterative) detection. The objective of detection is to determine the aposteriori (total or extrinsic) probabilities of the output of the Markov source. In order to leverage off of known results — even in the case of plain, unencoded transmission (no forward error correcting redundancy added by encoding) — one can view the output c of the Markov source (a lattice point, i.e. a vector of lattice coordinates) as the result of mapping with rate one (i.e. no additional redundancy) an identical replica of the input u*= c; this is a degenerated Markov process where even the dependence of the future on the present is removed. The only remaining structure to be captured for the Markov source, in the case when the candidate points are from a lattice, must reflect the partitioning in labeled cosets, as discussed in Section III-B. To this end, note that the labels themselves can be associated with states having integer values by virtue of the following convention: the state Sk-i at time k — 1 is the index of the label that contains the most recent lattice point output by the Markov source, i.e. at time k — 1; when the Markov source outputs a new point at time k it transitions into state S). equal to the integer indexing the label that contains the new point. Alternatively, with respect to the mapping u-& • c and omitting the time index, when M = A G A occurs at the rate-one block input, the Markov process transitions into the state whose (integer) value indexes the label containing λ. This is represented in Fig. 3, where e denotes an edge between starting state ss(e) and ending state sB(e). Formally, for any edge e, at any time, if ω(e)=λeA(Zj) cA, where ie {l, . . . |L(A)|} indexes one of the |£(A)| labels, then the ending state sB(e) = i and the Markov source outputs c(e) =u(e). There is a bijective mapping I between integer states and labels S-B- I3 such that, for any' integer state sg{l, . . .|L(Λ)|}, £(s) = l3 is the label associated with s.
The Markov sequence of random points selected from-the lattice can be thus viewed as triggered by state transitions triggered by u = λ 6 A; although the realizations of u on the lattice grid are random, a state model arises as a result of partitioning the lattice in equivalence classes. That is, there exist certain structural relations between certain points, which can be associated via an equivalence relation. The state probabilities, used in a posteriori probability calculations, are seen to be associated with the probabilities of these equivalence classes (or their labels), which can be obtained separately from belief propagation on the lattice's Tanner graph, as shown next.
In general, for a Markov process generated by triggering state transitions via some input (e.g. a classical convolutional code), the new state depends on the current input and several previous inputs; in the case at hand the new state depends only on the current input. This illustrates the degenerated nature of the Markov process at hand, seen thereby to be memoryless.
The memoryless nature of the Markov process is also apparent in the fact that any state can be reached in one transition from any state, and the probability distribution of the states does not depend on time; it depends only .on the probability distribution for u, and so does the probability distribution of the output of the Markov process. The output of the Markov process does not depend on the current state, but rather on the input u; the input determines both the new output and the new state, which implies that the output any time does not depend on any previous state.
The remainder of this appendix will use the state transition diagram in Fig. 3 for the Markov process that forms the object of detection; the results in [32], [33] apply. Following [32], the extrinsic APPs P^(C3I O) and PξF(uj; O) during the ft-th transition between states have the general expressions
wher re the probabilities of the current state and the new state that are associated with edge e.
Following the well-known results and notation in [33] and using the memoryless nature of the Markd&^iόiϊdsέiirϋFi^. B,
Aφ]
where, following [33], yj denotes the observations of the relevant Markov process, as taken at the output of a discrete memoriless channel at time instants 0, 1, . . . , T. Most importantly, the factor KQ does not depend on the state s, and is thereby cancelled out during the normalization step that enforces ∑
s A
k[s] = 1. Due to the isomorphism between states and labels it follows that Pr{5
fc = s; y'
k} is the label probability Pr (t(s)) — Pr (l
s) calculated as in (45). From [33] and the properties of the degenerated Markov process,
■
which does not depend on the state s and behaves as a constant that is cancelled out during the normalization step enforcing
REFERENCES
[1] M. Pohst, "Oa the computation of lattice vectors of minimal length, successive minima, and reduced basis with applications," ACM SIGSAM Bull., vol.
15, pp. 37-44, 1981 [2] U. Fincke and M. Pohst, "Improved methods for calculating vectors of short length in a lattice, including a complexity analysis," Math. Comput., vol.
44, pp. 463-471, Apr. 1985 [3] C. P. Schnoir and M. Euchnerr, "Lattice basis reduction: improved practical algorithms and solving subset sum problems," Math. Programming, vol. 66, pp. 181-191, Apr. 1994
[4] E. Viterbo and J. Boutros, "A -universal lattice decoder for fading channels," IEEE Trans. Inform. Theory, vol. 45, No. 5, pp. 1639-1642, July 1999. . [5] J. Boutros, N. Gresset, L. Brunei, and M. Fossorier, "Soft-input soft-output lattice sphere decoder for linear channels," Proc. IEEE Conf. Globecom '03, pp. 1583-1587, 2003. [6] S. Baro, J. Hagenauer, and M. Witzke, "Iterative detection of MIMO transmission using a list-sequential (LISS) detector," IEEE International Conf.
Commun., ICC '03, vol. 4, pp. 2653-2657, Anchorage, May 2003 [7] Y. de Jong and T. Willink, "Iterative tree search detection for MIMO wireless systems," IEEE 56th Vehicular Technology Conf. VTC'02 Proceedings, vol. 2, pp. 1041-1045, Sept 2002 [8] M. O. Damen, H. El-Gamal, and G. Caire, "On Maximum-Likelihood Detection and the Search for the Closest Lattice Point," IEEE Trans. Inform.
Theory, vol. 49, No. 10, pp. 2389-24"θ2, Oct. 2003. [9] A. H. Banihashemi and I. E Blake, "On the Trellis Complexity of Root Lattices and Their Duals," IEEE Trans. Inform. Theory, vol. 45, No. 6 , pp.
2168-2172, Sept. 1999. '
[10] S. Liu and Z. Tian, "Near-Optimum Soft Decision Equalization for Frequency Selective MIMO Channels," IEEE Trans. Signal Proc, vol. 52, No. 3, pp. 721-733, Mar. 2004. Λ
[11] B. Steingrimsson, Z.-Q. Luo, and K. M. Wong, "Soft quasi-maximum-likelihood detection for multiple-antenna wireless channels," IEEE Trans. Signal
Proc, vol. 51, No. 11, pp. 2710-2719, Nov. 2003. [12] B. Dong, X. Wang, and A Doucet, "A New Class of Soft MIMO Demodulation Algorithms," IEEE Trans. Signal Proc, vol. 51, No. 11, pp. 2752-2763,
Nov. 2003. [13] D. Seethaler, H. Artέs, and F. Hlawatsch, "Efficient approximate-ML detection for MMO spatial multiplexing systems by using a 1-D nearest neighbor search," , pp. 290-293.
[14] D. Le Ruyet, T. Bertozzi, and B. Ozbek, "Breadth first algorithms for APP detectors over MIMO channels," , pp. 926-930, 2004. [15] M. E. Ramirez Montalvo, "Sub-optimal non-linear detection methods for high data rate systems using multiple antennas," Master's Thesis, Electrical and Communications Engineering Department, Helsinki Univ. Technology, Nov. 20, 2003.
[16] A. K. Lenstra, H. W. Lenstra, and L. Lovasz, "Factoring polynomials with rational coefficients," Math. Ann., vol. 261, pp. 515-534, Apr. 1982 [17] A. H. Banihashemi and I. F. Blake, "Trellis structure and decoding complexity of lattices," Ph. D. dissertation, ECE Dept, Univ. Waterloo, Waterloo,
ON, Canada, 1997 [18] D. M. Ionescu, K. K. Mukkavilli, Z. Yan, and J. Lilleberg, "Improved 8- and 16-State Space-Time codes for 4PSK with Two Transmit Antennas," IEEE
Commun. letters, vol. 5, pp. 301-303, July 2001. [19] S. Siwamogsatham and M. P. Fite, "Improved High-Rate Space-Time Codes via Concatenation of Expanded Orthogonal Block Code and M-TCM,"
Proceedings of 2002 ICC, vol. 1, pp. 636-640, May 2002. [20] S. Siwamogsatham and M. P. Fitz, "Improved High-Rate Space-Time Codes via Orthogonality and Set Partitioning," Proceedings of 2002 IEEE WCNC, vol. 1, pp. 264-270, March 2002.
[21] N. Seshadri and H. Jafarkhani, "Super-Orthogonal Space-Time Trellis Codes," Proc. ICC'02, May 2002, Vol. 3, pp. 1439-1443. [22] H. Jafarkhani, N. Seshadri, "Super-orthogonal space-time trellis codes," IEEE Trans. Inform. Theory, vol. 49, pp. 937-950, Apr. 2003. ■ [23] D. M. Ionescu and Z. Yan, "Fading-Resilient Super-Orthogonal Space-Time Signal Sets: Can Good Constellations Survive in Fading?," submitted to
IEEE Trans. Infom. Theory; available on the ArXiv server http://www.arxiv.org/abs/cs.IT/0505049 [24] O. Tirkkonβn and A. Hottinen, "Square-matrix embeddable space-time block codes for complex signal constellations," IEEE Trans. Inform. Theory, vol.
48, pp. 384-395, Feb. 2002. [25] Z. Yan and D. M. Ionescu, "Geometrical Uniformity of a Class of Space-Time Trellis Codes," IEEE Trans. Inform. Theory, vol. 50, pp. 3343-3347,
Dec. 2004. [26] H. El Gamal, G. Caire, M. O. Damen, "Lattice coding and decoding achieve optimal diversity-multiplexing tradeoff of MIMO channels," IEEE Trans.
Inform. Theory, vol. 50, no. 6, pp. 968-985, June 2004.
[27] B. Hassibi and B. M. Hochwald, "High-rate codes that are linear in space and time," JEEE Trans. Inform. Theory, vol. 48, no. 7, pp. 1804-1824, July - 2002.
[28] B. Farhang-Boroujeny, Adaptive Filters: theory and applications. Chichester, West Sussex, England: Wiley, 2000. [29] A. H. BaniBashemi and F. R. Kschischang, 'Tanner graphs for group block codes and lattices: construction and complexity," IEEE Trans. Inform. Theory, vol. 47, Feb 2001
[30] M. C. Davey and D. MacKay, 'Xow-Density Parity Check Codes over GF(g)," ZEEE Commun. Letters, vol. 2, June 1998 [31] E. Agrell, T. Eriksson, A. Vardy, and K. Zeger, "Closest point search in lattices," ZEEE Trans. Inform. Theory, vol. 48, No. 2, pp. 2201-2214, Aug. 2002.
[32] S. Benedetto, D. Divsalar, G. Montorsi, and F. Pollara, "A soft-input soft-output APP module for iterative decoding of concatenated codes," IEEE
Commun. Lett., vol. 1, pp. 22-24-, Jan. 1997. [33] L. R. Bahl, J. Cocke, F. Jelinek, and J. Raviv, "Optimal decoding of linear codes for minimizing symbol error rate," IEEE Trans. Inform. Theory, vol.
IT-20, pp. 284-287, Mar. 1974. [34] D. M. Ionescu, D. Doan, S. Gray, "On Interleaving Techniques for MIMO Channels and Limitations of Bit Interleaved Coded Modulation," submitted to ZEEB Trans. Infom. Theory; available on the ArXiv server http://www.arxiv.org/abs/cs.IT/0510072.
Referring to Figure 8, an aspect of the exemplary embodiments of this invention resides in a method, such as one that may be used in a MIMO receiver. The method includes: (Block 8A) receiving a plurality of signals through a plurality of antennas, the plurality of signals being modulated with a space-time lattice code; (Block 8B) removing an effect of a channel matrix from the received signals to provide an equalized received signal; and (Block 8C) lattice detecting the equalized received signal based on a Tanner graph representation of the lattice.
The use of the exemplary embodiments of this invention enables and provides at least the advantages of soft output detection, no step back artifact generation, modularization of receiver implementation, wherein all practical constellations may be viewed as lattices (in a sense that they may be, e.g., degenerated lattices or cubic lattices). The use of the exemplary embodiments of this invention enables and provides a practical and efficient technique and means for decoding large constellations from multiple transmit antennas.
The exemplary embodiments of this invention can be applied to and used in, as non-limiting examples, E-UTRAN systems, OFDM-based systems, WCDMA systems, multi-carrier systems, so-called 3.9G (3.9 generation) systems and so-called 4G (fourth generation) systems, as well as in multi-band and multi-mode user equipment and terminals.
In general, the various embodiments may be implemented in hardware or special purpose circuits, software, logic or any combination thereof. For example, some aspects may be implemented in hardware, while other aspects may be implemented in firmware or software which may be executed by a controller, microprocessor or other computing device, although the invention is not limited thereto. While various aspects of the invention may be illustrated and described as block diagrams, flow charts, or using some other pictorial representation, it is well understood that these blocks, apparatus, systems, techniques or methods described herein may be implemented in, as non-limiting examples, hardware, software, firmware, special purpose circuits or logic, general purpose hardware or controller or other computing devices, or some combination thereof. Embodiments of the inventions may be practiced in various components such as integrated circuit chips and modules. The design of integrated circuits is by and large a
highly automated process. Complex and powerful software tools are available for converting a logic level design into a semiconductor circuit design ready to be etched and formed on a semiconductor substrate. Commercially available programs and systems may automatically route conductors and locate components on a semiconductor chip using well established rules of design, as well as libraries of pre-stored design modules. Once the design for a semiconductor circuit has been completed, the resultant design, in a standardized electronic format (e.g., Opus, GDSII, or the like) may be transmitted to a semiconductor fabrication facility or "fab" for fabrication.
Various modifications and adaptations may become apparent to those skilled in the relevant arts in view of the foregoing description, when read in conjunction with the accompanying drawings. However, any and all modifications of the teachings of this invention will still fall within the scope of the non-limiting embodiments of this invention.
Furthermore, some of the features of the various non-limiting embodiments of this invention may be used to advantage without the corresponding use of other features. As such, the foregoing description should be considered as merely illustrative of the principles, teachings and exemplary embodiments of this invention, and not in limitation thereof.