This application claims benefit of priority under 35 U.S.C. § 119 from prior United Kingdom Patent Application No. 0620071.1, filed Oct. 10, 2006, the entire contents of which are incorporated herein by reference.
The present invention is in the field of wireless communication, and particularly, but not exclusively, the field of multiple input, multiple output (MIMO) communications systems.
Conventional communication systems can be represented mathematically as:
y=Hx+v
in which, for a MIMO communication system, y is an n-by-1 vector representing the received signal, H is an n-by-m channel matrix modelling the transmission characteristics of the communications channel, x is an m-by-1 vector representing transmit symbols, v is an n-by-1 noise vector and wherein m and n denote the number of transmit and receive antennas respectively.
It will be understood by the skilled reader that the same representation can be used for multi-user detection in CDMA systems.
Recent publications have demonstrated how the use of a technique called Lattice Reduction can improve the performance of MIMO detection methods.
For example, “Lattice-Reduction-Aided Detectors for MIMO Communication Systems”, (H. Yao and G. W. Wornell, Proc. IEEE Globecom, November 2002, pp. 424-428) describes Lattice-reduction (LR) techniques for enhancing the performance of multiple-input multiple-output (MIMO) digital communication systems.
In addition, “Low-Complexity Near-Maximum-Likelihood Detection and Precoding for MIMO Systems using Lattice Reduction”, (C. Windpassinger and R. Fischer, in Proc. IEEE Information Theory Workshop, Paris, March, 2003, pp. 346-348) studies the lattice-reduction-aided detection scheme proposed by Yao and Wornell. It extends this with the use of the well-known LLL algorithm, which enables the application to MIMO systems with arbitrary numbers of dimensions.
“Lattice-Reduction-Aided Receivers for MIMO-OFDM in Spatial Multiplexing Systems”, (I. Berenguer, J. Adeane, I. Wassell and X. Wang, in Proc. Int. Symp. on Personal Indoor and Mobile Radio Communications, September 2004, pp. 1517-1521, hereinafter referred to as “Berenguer et al.”) describes the use of Orthogonal Frequency Division Multiplexing (OFDM) to significantly reduce receiver complexity in wireless systems with Multipath propagation, and notes its proposed use in wireless broadband multi-antenna (MIMO) systems.
Finally, “MMSE-Based Lattice-Reduction for Near-ML Detection of MIMO Systems”, (D. Wubben, R. Bohnke, V. Kuhn and K. Kammeyer, in Proc. ITG Workshop on Smart Antennas, 2004, hereinafter referred to as “Wubben et al.”) adopts the lattice-reduction aided schemes described above to the MMSE criterion.
The techniques used in the publications described above use the concept that mathematically, the columns of the channel matrix, H, can be viewed as describing the basis of a lattice. An equivalent description of this lattice (a so-called ‘reduced basis’) can therefore be calculated so that the basis vectors are close to orthogonal. If the receiver then uses this reduced basis to equalise the channel, noise enhancement can be kept to a minimum and detection performance will improve (such as, as illustrated in
yr xr and Hr are defined to be the real-valued representations of y, x, and H respectively, such that:
where Re( ) and Im( ) denote the real and imaginary components of their arguments.
It will be noted that Berenguer et al. describes the equivalent method in the complex plane, though for the purpose of clarity the Real axis representation of the method is used herein.
A number of lattice reduction algorithms exist in the art. One suitable lattice reduction algorithm is the Lenstra-Lenstra-Lovasz (LLL) algorithm referred to above, which is disclosed in Wubben et al., and also in “Factoring Polynomials with Rational Coefficients”, (A. Lenstra, H. Lenstra and L. Lovasz, Math Ann., Vol. 261, pp. 515-534, 1982, hereinafter referred to as “Lenstra et al.”), and in “An Algorithmic Theory of Numbers, Graphs and Convexity”, (L. Lovasz, Philadelpia, SIAM, 1980, hereinafter referred to as “Lovasz”).
Any one of these can be used to calculate a transformation matrix, T, such that a reduced basis, {tilde over (H)}r, is given by
{tilde over (H)}r=HrT
The matrix T contains only integer entries and its determinant is +/−1.
After lattice reduction, the system is re-expressed as:
where zr=T−1xr. The received signal, yr, in this redefined system is then equalised to obtain an estimate of zr. This equalisation process then employs, for example, a linear ZF technique, which obtains:
{tilde over (z)}r=({tilde over (H)}r*{tilde over (H)}r)−1{tilde over (H)}r*yr
Since {tilde over (H)}r is close to orthogonal, {tilde over (z)}r should suffer much less noise enhancement than if the receiver directly equalised the channel Hr.
Of course, other equalisation techniques could be used. For example, MMSE techniques, or more complex successive interference cancellation based methods, such as in the published prior art identified above, could be considered for use.
A receiver in accordance with the above operates in the knowledge that the transmitted symbols contained in x are obtained from an M-QAM constellation. With this constraint, xr can be formed as xr=αs+β, where s is taken from the set of integers limited by the dimension of the lattice, and α and β are scalar values obtained from the definition of the M-QAM constellation in use. α is equal to the minimum distance between two constellation points and β is equal to the offset from the origin when s=0. In the present example, a 16-QAM constellation is used, having real and imaginary components of {+/−1, +/−3}. Therefore, as shown in
{tilde over (z)}r can then be quantised as {circumflex over (z)}r, in accordance with the method indicated in Wubben et al.:
where Q { } is a quantisation function which rounds each element of its argument to the nearest integer, and where 1 is a column vector of ones.
It will be understood from the above that, the quantisation function apart, the remaining operations are a result of M-QAM constellations being scaled and translated versions of the integer lattice. The integer quantisation therefore requires the same simple scaling and translation operations.
Finally, the estimate {tilde over (x)}r of xr is obtained by this method as
{circumflex over (x)}r=T{circumflex over (z)}r
Occasionally, if errors are present in the estimate of {circumflex over (z)}r, then it is possible that some of the symbol estimates in {circumflex over (x)}r, may not be valid symbols. In such cases, these symbols are mapped to the nearest valid symbol. For example, for the present example employing 16-QAM, the values +/−1, +/−3 may define the valid entries in {circumflex over (x)}r. Therefore if a component of {circumflex over (x)}r were, for example, equal to +5, then this would be mapped to a value of +3.
Such reduced lattice detectors (e.g. for MIMO systems) usually output hard decisions. The only mention in the literature of a technique that could be employed for obtaining soft-output is “From Lattice-Reduction-Aided Detection Towards Maximum-Likelihood Detection in MIMO Systems”, (C. Windpassinger, L. Lampe and R. Fischer, in Proc. Int. Conf. on Wireless and Optical Communications, Banff, Canada, July 2003, hereinafter referred to as “Windpassinger et al.”). The method that Windpassinger et al. proposes is complex, and the performance of this technique was not validated in the publication.
U.S. Pat. No. 6,724,843 describes a detector in which received symbols are decoded in a multiple-antenna communication system using lattice-based decoding. The symbols are generated using a modulation constellation, e.g., a diagonal modulation constellation, and the constellation is characterized as a lattice for decoding purposes. For example, if a given communication link of the multiple-antenna communication system includes M transmitter antennas and a single receiver antenna, the diagonal modulation constellation can be characterized as a lattice in M dimensions. A differential decoding operation for received differential symbols involves a determination of the closest point in the lattice corresponding to the constellation. This determination may be made in an efficient manner using a basis reduction algorithm which generates an approximately orthogonal basis for the lattice, and then utilizes component-wise rounding to determine the closest point. The lattice-based decoding has a complexity which is polynomial rather than exponential in the particular number of antennas and the system rate, but is nonetheless able to deliver error rate performance which approximates that of maximum likelihood decoding.
Therefore, with the exception of Windpassinger et al., the reduced lattice detection schemes described in all other references only output hard decisions for the estimate of the transmitted symbol vector, {circumflex over (x)}. When used in a system with an outer channel code (i.e. in any practical system) the performance of the code can be substantially improved if it is supplied with soft-information, e.g. a log-likelihood ratio (LLR) for each bit.
UK patent application 0518036.9 (as yet unpublished) describes an approach which makes use of the above and further presents a method of obtaining soft estimates using candidate vectors in a reduced basis and selecting the most likely candidate on the basis of probabilities determined for the candidate vectors.
An aspect of the invention provides a method for determining soft estimates of transmitted bit values from a received signal in a lattice-reduction-aided receiver based wireless communications system, the method comprising obtaining an estimate of the channel response, applying lattice reduction to said channel response and equalising said received signal in accordance with the reduced basis channel, and determining probabilities of transmitted bits having particular values by means of selecting a set of candidate symbols by selecting a core candidate vector in the reduced basis, determining a set of mapping functions between vectors in the reduced basis and transmitted symbol vectors from which a set of candidate transmitted symbol vectors can be generated, generating the set of candidate transmitted signal vectors on the basis of the mapping functions and, on the basis of the received signal, determining the probability of each transmitted bit value estimate having been transmitted.
Preferably, the method provides applying lattice reduction in accordance with the LLL algorithm.
The method is suitable for use in a MIMO wireless communications system. Further, it can be used with any other system wherein a received signal is the result of transmission from a plurality of antennas, which may or may not be collocated. Further, the method is applicable to CDMA systems, such as multi-user detection (MUD).
Yet another aspect of the invention provides a receiver for use in a lattice-reduction-aided receiver based wireless communications system, the receiver comprising means for obtaining an estimate of the channel response and a detector operable to process a received signal so as to determine soft estimates of transmitted bit values, the detector comprising means for applying lattice reduction to the channel response estimate and equalising the received signal in accordance with the reduced basis channel, and soft information determining means for determining probabilities of transmitted bits having particular values, the soft information determining means comprising of means for selecting a core candidate symbol vector in the reduced basis, means for determining a set of mapping functions between vectors in the reduced basis and transmitted symbol vectors from which a set of candidate transmitted symbol vectors can be generated, means for generating a set of candidate transmitted signal vectors on the basis of said mapping functions and, means for determining the probability of said transmitted bit value having been transmitted, on the basis of the received signal.
Another aim of the present invention is to provide a method for obtaining log-likelihood ratios (LLRs) for bits at the output of a lattice-reduction-aided MIMO receiver.
To this end, the probabilities determined in either the method above or by the detector above can be converted into LLRs.
The reader will appreciate that aspects of the invention can also be embodied by means of computer implemented means, for instance under the direction of a suitable software product. Such a software product can be in the form of a storage medium, for instance a computer readable optical or magnetic disk, or by way of a signal borne on a computer receivable carrier.
The present invention will now be described with reference to an implementation thereof for the equalization of a wireless communication system.
The communications system 10 comprises a transmitter device 12 and a receiver device 14. It will be appreciated that, in many circumstances, a wireless communications device will be provided with the facilities of a transmitter and a receiver in combination but, for this example, the devices have been illustrated as one way communications devices for reasons of simplicity.
The transmitter device 12 comprises a data source 16, which provides data (comprising information bits or symbols) to a channel encoder 18. The channel encoder 18 is followed by a channel interleaver 20 and, in the illustrated example, a space-time encoder 22. The space-time encoder 22 encodes an incoming symbol or symbols as a plurality of code symbols for simultaneous transmission from a transmitter antenna array 24 comprising a plurality of transmit antennas 25. In this illustrated example, three transmit antennas 25 are provided, though practical implementations may include more, or fewer antennas depending on the application.
The encoded transmitted signals propagate through a MIMO channel 28 defined between the transmit antenna array 24 and a corresponding receive antenna array 26 of the receiver device 14. The receive antenna array 26 comprises a plurality of receive antennas 27 which provide a plurality of inputs to a lattice-reduction-aided decoder 30 of the receiver device 14. In this specific embodiment, the receive antenna array 26 comprises three receive antennas 27 although, again, this number is not essential for performance of the invention.
The lattice-reduction-aided decoder 30 has the task of removing the effect of the MIMO channel 28. The output of the lattice-reduction-aided decoder 30 comprises a plurality of signal streams, one for each transmit antenna 25, each carrying so-called soft or likelihood data on the probability of a transmitted bit having a particular value. This data is provided to a channel de-interleaver 32 which reverses the effect of the channel interleaver 20, and the de-interleaved bits output by this channel de-interleaver 32 are then presented to a channel decoder 34, in this example a Viterbi decoder, which decodes the convolutional code. The output of channel decoder 34 is provided to a data sink 36, for further processing of the data in any desired manner.
The specific function of the lattice-reduction-aided decoder 30 will be described in due course.
Audio/video output hardware devices 120 are further connected to the general purpose bus 114, for the output of information to a user. Audio/video output hardware devices 120 can include a visual display unit, a speaker or any other device capable of presenting information to a user.
Communications hardware devices 122, connected to the general purpose bus 114, are connected to the antenna 26. In the illustrated embodiment in
Communications facilities 132 in accordance with the specific embodiment are also stored in the working memory 112, for establishing a communications protocol to enable data generated in the execution of one of the applications 130 to be processed and then passed to the communications hardware devices 122 for transmission and communication with another communications device. It will be understood that the software defining the applications 130 and the communications facilities 132 may be partly stored in the working memory 112 and the mass storage device 116, for convenience. A memory manager could optionally be provided to enable this to be managed effectively, to take account of the possible different speeds of access to data stored in the working memory 112 and the mass storage device 116.
On execution by the processor 110 of processor executable instructions corresponding with the communications facilities 132, the processor 110 is operable to establish communication with another device in accordance with a recognised communications protocol.
The function of the lattice-reduction aided decoder 30 will now be described in further detail in accordance with
In step S1-2, the vector {circumflex over (z)}r is taken as the core candidate vector. In step S1-4 a set of mapping functions is generated. Each mapping function is a means for the generation of another candidate transmitted symbol vector from the core transmitted symbol vector. The manner in which this is achieved, in accordance with the present embodiment, will now be described.
Expressed in the reduced basis, the additional list of candidate vectors to be considered as a transmitted symbol vector consists of all vectors that vary from the core candidate vector by modification of one or more elements of the vector {circumflex over (z)}r. Whilst any of these additional candidate vectors may differ from {circumflex over (z)}r in more than one element, the example described herein generates candidates by only allowing these to vary one element of {circumflex over (z)}r. Creating candidate vectors by allowing perturbations to multiple elements of {circumflex over (z)}r can slightly improve performance, but at the expense of increasing the length of the candidate list and hence increasing complexity.
For subsequent use in the decoder, candidate vectors need to be expressed in terms of the original lattice basis. To do this, the above mentioned candidate vectors would need to be converted from the reduced basis into the original basis. Such conversions require the use of matrix calculations, and this action would consume considerable computational resources.
By generating a set of mapping functions between the core candidate and additional candidates, it is only necessary to carry out the above conversion on the core candidate. Mapping functions that do not require full matrix calculations can be created, and so such functions will require less computation than would otherwise be possible.
In the present embodiment, {circumflex over (x)}r, is a 2-by-1 vector, so the set of 5 candidate vectors, c(1) expressed in the reduced lattice, are as follows:
In the original lattice, these candidate vectors are expressed as:
{circumflex over (x)}r(2) for example, can be expressed as:
Hence by the same method,
{circumflex over (x)}r(3)={circumflex over (x)}r(1)−αT(:,1)
{circumflex over (x)}r(4)={circumflex over (x)}r(1)+αT(:,2)
{circumflex over (x)}r(5)={circumflex over (x)}r(1)−αT(:,2)
where T(:,i) represents the ith column of matrix T.
Thus it is only necessary to perform the full matrix multiplication once, calculating for example {circumflex over (x)}r(1), in order to generate the full set of candidate transmitted symbol vectors.
Although in the present embodiment, {circumflex over (x)}r is a 2-by-1 vector, it will be appreciated that a similar approach could be taken in other embodiments of the invention wherein vectors of different dimension are used.
It will be appreciated that other methods of generating suitable candidates are possible, and would provide a suitable list for further processing in accordance with the method. One alternative method would be to perturb not just one element at a time, but to perturb two simultaneously. Therefore, for each pair of elements, the decoder would cycle through the four combinations of perturbing by +/−α. Then, the decoder would pick the next pair of elements and repeat in order to generate more candidates.
Candidate transmitted symbol vectors are then obtained in step S1-6 by applying the mapping functions calculated in step S1-4 to the core candidate transmitted symbol vector.
Just as for the hard-output detector outlined in the introduction and discussion of the prior art above, it is possible that some of the elements of the vector {circumflex over (x)}r(i) may not be valid symbols. Therefore, step S1-8 seeks to determine if this is the case, and, if so, in step S1-10, these symbols are mapped to the nearest valid symbol. (e.g. for 16-QAM, if the values +/−1, +/−3 define the valid entries as illustrated in
For each candidate symbol vector {circumflex over (x)}r(i) (corrected, if required), in step S1-12, the detector calculates its probability of being transmitted, as:
These probabilities are then used to calculate, in step S1-14, the probability of symbol x′ having been transmitted from antenna k, where x′εX and X defines the set of symbols in the chosen constellation.
Depending on the list of candidates, according to the above definition P may not be specified for all values of k and x′. In these cases P is set to a default (small) value. This default can be a fixed value or it could varied according to a method such as that described in “Adaptive Selection of Surviving Symbol Replica Candidates Based on Maximum Reliability in QRM-MLD for OFCDM MIMO Multiplexing” (K. Higuchi, H. Kawai, N. Maeda and M. Sawahashi, in Proc. IEEE Globecom, Dallas, December 2004), or by any other appropriate method.
Now that the receiver has information on the probability of different symbols having been transmitted, these are processed in a conventional manner in step S1-16 to obtain a log-likelihood ratio for each transmitted bit. In the present example, this is done as follows:
where L(bk,i) is the log-likelihood ratio of bit bk,i, k indicates the transmit antenna, i=1, . . . , M where M is the number of bits per symbol, and where X(1) and X(0) are the sets of symbols for which bk,i=1 and bk,i=0 respectively.
The graph of
It will be appreciated that the foregoing disclosure of specific embodiments of the invention can be applied to any communications product employing MIMO transmission techniques, to take advantage of the benefits of the invention. Further, the invention is applicable to any circumstance in which the detection of symbols, which may be based on multiple input, is required. This could arise in systems where a plurality of antennas are provided in separate locations. Further, CDMA MUD may be a suitable basis for use of the method of the present invention.
The invention has been described by way of a software implementation. This software implementation can be introduced as a stand alone software product, such as borne on a storage medium, e.g. an optical disk, or by means of a signal. Further, the implementation could be by means of an upgrade or plug-in to existing software, or by introduction of solid state memory means storing the relevant executable instructions.
Whereas the invention can be so provided, it could also be by way exclusively by hardware, such as on an ASIC.
The reader will appreciate that the foregoing is but one example of implementation of the present invention, and that further aspects, features, variations and advantages may arise from using the invention in different embodiments. The scope of protection is intended to be provided by the claims appended hereto, which are to be interpreted in the light of the description with reference to the drawings and not to be limited thereby.
Number | Date | Country | Kind |
---|---|---|---|
0620071.1 | Oct 2006 | GB | national |