This invention relates to a Multiple Input, Multiple Output communications system, and to a receiver and a method of symbol detection for use in such a system.
Multiple Input, Multiple Output (MIMO) communications systems take advantage of spatial multiplexing to increase wireless bandwidth and range. Specifically, MIMO transmitters send information out using two or more antennas, and the information is received via multiple antennas as well. MIMO systems use the additional pathways to transmit more information, and then recombine the signal on the receiving end. MIMO systems provide a significant capacity gain over conventional single antenna systems, along with more reliable communication. MIMO-based transceivers can for example be employed in WLAN 802.11n, WiMax and cellular communications systems.
However, with this increase in capacity comes an increase in calculation complexity. It is a nontrivial task to separate out in the receiver the symbols transmitted from each of the transmit antennas, with the optimal joint decoding techniques requiring large numbers of calculations to be performed in order to obtain a solution. Moreover, the noise added to each signal during its transit between the transmitter and the receiver, largely in the receiver front-end amplifier, further complicates matters, meaning suboptimal, but more rapid, symbol detection techniques can have a large symbol error rate.
Digital communications systems often use signal space diagrams to represent signals, for example the signal being transmitted from a transmitter. For example, in systems using quadrature amplitude modulation (QAM), the in-phase and quadrature components of a signal represent the data being transmitted. Each point represents a symbol, a unique signal state of a modulation scheme which conveys one or more user bits to the receiver. A signal space diagram showing all the possible transmitted symbols is known as a constellation. In MIMO systems, each transmit antenna transmits a symbol, and the set of transmitted symbols at any time forms a symbol vector.
The task of the MIMO receiver is to determine the transmitted symbol vector in each time period, using the detected symbols.
The document US 2004/0066866 discloses a method of decoding space-time coded signals transmitted from a number of transmit antennas. First, a separate detection technique is used to determine initial decoding solutions corresponding to the symbols transmitted from each of a number of transmit antennas at a given time. For each initial solution, a limited area about the initial solution is defined. Each of the limited areas will correspond to regions including constellation points proximate to the initial solution. The initial solutions are used to define a limited, multi-dimensional space. Finally, a joint decoding technique is implemented within the limited space to find a final solution.
However, the multi-dimensional space defined by these initial solutions is still relatively large, meaning that a relatively large number of calculations need to be performed in the joint decoding step.
It is an object of the present invention to provide a method of symbol detection that requires a smaller number of calculations in a decoding step.
The present invention provides a method of identifying transmitted symbol vectors as part of a communications system, wherein a plurality of transmitting antennas each transmit a respective symbol during a time period, each symbol being selected from a plurality of possible transmitted symbols, and the plurality of possible transmitted symbols being represented by a constellation plane, wherein the plurality of symbols transmitted by said plurality of transmitting antennas form a transmitted symbol vector, and wherein the method comprises: receiving a signal over a channel originating from the plurality of transmitting antennas; applying a first algorithm to the received signal to obtain an initial solution of the transmitted symbol vector comprising a plurality of initial values for the transmitted symbols; hard-demapping each of said initial values to one of said possible transmitted symbols in the constellation plane, in order to form a respective estimated transmitted symbol, the set of estimated transmitted symbols comprising an estimated transmitted symbol vector; defining a selected area in the constellation plane about each estimated transmitted symbol; generating a list of candidate symbol vectors, each candidate symbol vector comprising symbols that are within the respective selected areas surrounding each estimated transmitted symbol, and each candidate symbol vector differing from the estimated transmitted symbol vector only in a subset of the symbols; and applying a decoding technique to the list of candidate symbol vectors.
This has the advantage that the decoding step requires a smaller number of calculations, as the size of the search space is reduced. Moreover, although the decoding step requires a smaller number of calculations, the decoding performance of the receiver is not seriously adversely affected.
According to a second aspect of the present invention, there is provided a receiver operating in accordance with the method of the first aspect of the invention.
According to a third aspect of the present invention, there is provided a communications system, in which the receiver operates in accordance with the method of the first aspect of the invention.
Specifically, by limiting the search space to include only vectors that differ from the initial estimate of the transmitted symbol vector in a restricted number of symbols, the number of calculations can be significantly reduced, without having excessively damaging effects on the symbol error rate.
These and other aspects of the invention will be apparent from and elucidated with reference to the embodiments described hereinafter.
The transmitter system 12 includes first transmitter circuitry 16, connected through first RF circuitry 18 to a first transmit antenna 20. The transmitter system 12 also includes second transmitter circuitry 22, connected through second RF circuitry 24 to a second transmit antenna 26. The transmitter system 12 is generally conventional, and will not be described further herein. Although
It will further be appreciated that, although
Signals from the transmitter system 12 are transmitted from the antennas 20, 26 over an air interface to the receiver system 14.
The receiver system 14 includes two receive antennas 28, 30. Again, although only two receive antennas are shown in
The first receive antenna 28 is connected to first RF receiver circuitry 32, and the output of the first RF receiver circuitry 32 is connected to a first sampling block 34, for forming digital samples of the signal received at the first receive antenna 28. The digital samples are passed to a first FFT block 36, for conversion to the frequency domain. As mentioned above, in this illustrated embodiment, the invention is applied to an OFDM system, and so this conversion to the frequency domain is required. However, the invention is equally applicable to other communications systems.
The output of the first FFT block 36 is passed to a symbol/bit detection block 38. The output demapped symbols are passed to a deinterleaver and decoder block 40, for forming a decoded output signal.
The second receive antenna 30 is connected to second RF receiver circuitry 42, and the output of the second RF receiver circuitry 42 is connected to a second sampling block 44, for forming digital samples of the signal received at the second receive antenna 30. The digital samples are passed to a second FFT block 46, and the output of the second FFT block 46 is passed to the symbol/bit detection block 38. The output demapped symbols are passed to the deinterleaver and decoder block 40.
The symbol/bit detection block 38 includes a zero-forcing detector block 100. The zero-forcing detector block 100 receives Nr inputs, r1, . . . , rN
The operation of each block in the symbol/bit detection block 38 will be described in more detail below, and further with reference to
A model for wireless communication over an Nt×Nr MIMO channel (with Nt transmit antennas, Nr receive antennas) can be given by the following equation:
r=H*s+n
where r is the received vector, H is the channel matrix, with Nr rows and Nt columns, s is the transmitted vector and n is the noise.
The simplest symbol detection technique uses the zero-forcing (ZF) algorithm. This method applies the inverse of the channel matrix H to the received vector r to obtain an output sZF:
s
ZF
=H
−1
*r=s+H
−1
*n.
The problem with zero forcing is that the inverse of the channel matrix will on average amplify the noise at the receiver, thus making symbol detection more erroneous. Zero forcing, even though simple to implement, gives noisy estimates of the received symbol vector and is a sub-optimal algorithm.
The optimal symbol detection technique is maximum-likelihood detection (MLD), but used conventionally it is computationally expensive. The technique works by computing the Euclidean distance between the received vector r and the product of the channel matrix and a possible transmitted vector si:
d
i
2
=∥r−H*s
i∥2.
However, the conventional method is to determine the solution after performing this calculation for all MN
A process flow is described for carrying out the method of the present invention.
First, a separate detection technique, such as zero forcing or minimum mean square error, is used to determine initial solutions sinit1, . . . , sinitN
Each initial solution is then hard-demapped on to a nearby one of the nearest possible transmitted symbols, generating an estimated transmitted symbol vector, sest=(sest1, sest2, . . . sestN
A selected area in the constellation plane is then defined around each of the symbol estimates, encompassing nearby possible transmitted symbols in addition to the symbol estimate itself.
A list of possible transmitted symbol vectors is then generated, including the estimated transmitted symbol vector itself, and all symbol vectors that differ from that vector only in a certain number of symbols and such that the differing symbols are within the selected area defined above.
Finally, a joint detection technique, such as maximum-likelihood detection (MLD), is applied to the vectors in the list to find a final solution. Accordingly, the initial solution is used to define a limited search space for the subsequent joint detection technique, reducing the implementation complexity considerably.
A preferred embodiment of the present method is described in
Accordingly, the MIMO channel matrix H is first determined (step S2). It is to be noted that, although here we assume H to be accurate, this is a nontrivial task and normally H will be an estimate. However, the person skilled in the art will be aware of many channel estimation techniques for determining H, and so this step is not described further. As mentioned above, the channel matrix H is an Nt×Nr matrix, where Nt is the number of transmit antennas and Nr is the number of receive antennas.
Next, the inverse of the MIMO channel matrix is computed (step S4). Preferably, for MIMO systems having different numbers of transmit and receive antennas, the pseudo-inverse, pseudoinv, will be used, as defined below:
pseudoinv(H)=(H*×H)−1×H*
The channel matrix inverse is then applied to the received vector r according to the ZF approach:
s
init
=H
−1
*r=s+H
−1
*n
to obtain an initial ZF solution symbol vector sinit, where s is the actual transmitted symbol vector and n is the noise (step S6).
Each symbol of the initial solution is then hard-demapped to the possible transmitted symbol nearest to it in the constellation plane (step S8), defining an estimated transmitted symbol (see
In step S6 of the process shown in
The value of the initial solution sinit1 for the symbol transmitted from the first transmit antenna is represented in
The action of the hard-demapper, in step S8 of the process shown in
Returning to
The estimated transmitted symbols sest1 and sest2 are represented as dots in circles 1521 and 1522. The four symbols included in the areas defined by step S10 of
In one embodiment of the invention, if sest is at a corner of the constellation plane then only the two neighbors are considered to be in the selected area, while if sest is on the edge of the constellation plane then only the three neighbors are considered to be in the selected area. However, other possibilities exist. For example, other symbols could be included in the selected area such that the selected area contains the same number of symbols irrespective of the position of sest within the constellation plane.
Returning to
Therefore, for the illustrated example having two transmit antennas with estimated transmitted symbols sest1 and sest2 represented by sj and sk, the list of candidate transmitted symbol vectors is as follows: (sj, sk), (sj+e, sk), (sj−e, sk), (sj+e.i, sk), (sj−e.i, sk), (sj, sk+e), (sj, sk−e), (sj, sk+e.i), (sj, sk-e.i), where i=√−1, such that adding or subtracting e involves moving a distance e in the horizontal direction in the complex number plane, while adding or subtracting e.i involves moving a distance e in the vertical direction in the complex number plane.
In step S14 of the process shown in
d
i
2
=∥r−H*s
i∥2.
Thus, the list of candidate transmitted symbol vectors generated in the nearby symbol list generator (block 104 in
Similarly, the symbols sj and sk+e, which together form a second of the candidate symbol vectors, are input to a second calculation block 1602, and so on, with the symbols sj−e.i and sk, which together form the nth of the candidate symbol vectors, are input to the nth calculation block 160n. As shown in
In step S16 of the process shown in
As described above, steps S14 and S16 find the symbol that has the maximum likelihood of having been the transmitted symbol. Thus, the technique described above can be regarded as maximum likelihood symbol detection.
As an alternative, it is possible to perform a maximum likelihood bit detection technique, as described below. Specifically, the symbol/bit detector 38 passes to the deinterleaver and decoder 40 a set of bit-metrics, based on the list of candidate vectors obtained in step S12 of the process of
In one embodiment of the invention, the bit-metrics can be obtained from log likelihood ratios (LLRs) for each of the bits of the symbol.
Thus, for one bit position in the symbol vector, we form a measure of the likelihood that the transmitted symbol contained a “1” in this bit position, relative to the likelihood that the transmitted symbol contained a “0” in this bit position. The LLR, L, is the log of the quotient of the probability that the transmitted bit bk was a “0” (given the received symbol vector, r), and the probability that it was a “1” (also given the received symbol vector, r). That is:
In this embodiment of the invention, the computation is simplified by considering only those possibly transmitted symbol vectors that are within the set of candidate vectors identified in step S12 of the process of
The posterior probability of each symbol vector is proportional to exp(−∥r−Hs∥2). So the total probability that the transmitted bit was a “0”, is equal to the sum over all vectors in X0 of this a posteriori probability, while the total probability that the transmitted bit was a “1” is equal to the sum over all vectors in X1 of this a posteriori probability. That is:
Again, reducing the size of the search space, by generating a list of candidate vectors as described above, greatly reduces the complexity of these calculations, surprisingly without a correspondingly adverse effect on the results.
There is thus described a method, and an apparatus, for decoding received signals in a MIMO receiver. Although it can be seen that the method and apparatus of the invention are highly efficient, in terms of the number of calculations required in the joint decoding step, and in terms of the implementation of the system in a hardware device, they achieve symbol error rates that are surprisingly close to the symbol error rates achieved by the optimal MLD symbol detector.
Number | Date | Country | Kind |
---|---|---|---|
05112141.6 | Dec 2005 | EP | regional |
PCT/IB2006/054668 | Dec 2006 | IB | international |
Filing Document | Filing Date | Country | Kind | 371c Date |
---|---|---|---|---|
PCT/IB2006/054668 | 12/7/2006 | WO | 00 | 3/11/2009 |