LAYERED SPACE TIME PROCESSING IN A MULTIPLE ANTENNA SYSTEM
This invention was made with United States government support awarded by the following agencies: NSF CCR-9875805, NSF ECS-9979408, and DARPA F30602-00-2-
0555. The United States has certain rights in this invention.
BACKGROUND OF THE INVEN ION
1. Field of Invention The present invention is directed to wireless communication systems. In particular, the present invention is directed to wireless communication systems utilizing multiple antenna arrays.
2. Description of Related Art Presently, multiple antenna arrays are used for transmitting data in wireless communication systems. For example, multiple antennas are used at both the transmitter and at the receiver for transmitting data. These multiple antenna arrays can increase wireless channel capacity linearly by the number of transmit antennas, providing the number of receive antennas is greater or equal to the number of transmit antennas. Unfortunately, the relatively high dimensional nature of multiple antenna array systems has high computational complexity in practical systems.
One example of a wireless communication system for generating a high data rate using a multiple antenna array is disclosed in U.S. Patent No. 6,097,771, issued August 1, 2000. This system is otherwise known as a Bell Laboratories Layered Space-Time (BLAST) system. The BLAST system utilizes multiple transmit antennas and multiple receive antennas along with signal processing for interference cancellation.
An example of such systems utilizes, for example, four transmitter antennas and four receiver antennas. This system can create four independent sub-channels within a single bandwidth by coding the sub-channels individually as opposed to coding them jointly. Thus,
the system increases capacity within a bandwidth by four-fold by exploiting the idea of diversity such that each channel corresponds with transmitting over many Raleigh fading channels. Accordingly, even if one channel is bad, it is combined with the other good channels to reduce the fluctuation of signal strength. This is because each channel then has enough diversity so that each one ends up appearing as a traditional additive white Gaussian noise channel from clear signal processing. However, the creation of the sub-channels requires signal processing to remove interference between the sub-channels.
In operation, the BLAST system performs signal processing to create four parallel sub-channels in a sequential procedure. This sequential method involves removing channel interference by sequentially subtracting the signal of one sub-channel from the other subchannels. Unfortunately, the system ends up starting with the sub-channel with the lowest diversity. This is a problem because the sub-channel with the lowest diversity order is, in essence, the worst channel. The higher diversity sub-channels are created by subtracting decisions from the worst sub-channel to create more reliable sub-channels. Therefore, this method is backward in the sense that it starts with the least reliable decision to create successive, more reliable sub-channels. Thus, problems with the worst sub-channel can create problems with the more reliable sub-channels. Also, this method is excessively complex and thus requires a significant amount of processing power. Furthermore, this method has the problem in that it creates unequal levels of diversity. For example, in a four sub-channel system, one sub-channel has a first order diversity, another sub-channel has a second order diversity, another sub-channel has a third order diversity, and the last subchannel has a fourth order diversity. This creates more unreliability because the lower order diversity sub-channels are more unreliable.
SUMMARY OF THE INVENTION
The present invention provides a more reliable wireless communication system utilizing a multiple antenna array. The system is more reliable at least because it begins with a higher order diversity sub-channel. Additionally, the present invention provides a less complex method of determining initial decisions on sub-channels. Also, the present
invention is very efficient. Furthermore, the present invention provides dramatic performance improvements over traditional systems.
These features, along with other features described below, are provided for by a system and method for performing space-time processing method for removing interference for a multiple antenna array system. A decision of a desired sub-channel of the signal vector is generated by nulling out the interference of a first set of sub-channels by multiplying the received signal vector by a unitary matrix generated from a QR decomposition of a channel matrix. An improved symbol decision is generated by successively canceling channel interference due to a second set of sub-channels. The symbol represents transmitted bits, transmitted coded bits, or other useful transmitted information. For example, the symbol represents the decoded or coded bits of a digital transmission.
According to another embodiment, the system and method performs an extended space-time processing. An improved symbol decision of a desired sub-channel of the signal vector is generated by generating a baseline decision for the sub-channel. This baseline decision can be generated by the baseline method, by the BLAST system, or by any other method useful for generating a decision for a sub-channel. Next, a contribution of a strongest sub-channel is subtracted from the signal vector to generate a modified signal vector. The modified signal vector is multiplied by a unitary matrix generated from a QR decomposition of another channel matrix. Channel interference of the remaining sub-channels of the modified signal vector is successively cancelled from a remaining sub-channel.
BRIEF DESCRIPTION OF THE DRAWINGS The preferred embodiments of the present invention will be described with reference to the following figures, wherein like numerals designate like elements, and wherein: Fig. 1 is an exemplary block diagram of a wireless system according to one embodiment;
Fig. 2 is an exemplary flowchart outlining the operation of the processor according to a baseline method;
Fig. 3 is an exemplary illustration of a space-time processing schematic for a six transmitter and six receiver system according to one baseline embodiment;
Fig. 4 is an exemplary illustration of a schematic for processing similar to BLAST for a (6, 6) system according to one embodiment; Fig. 5 is an exemplary illustration of performance analysis of the baseline method according to one embodiment;
Fig. 6 is an exemplary flowchart outlining the operation of the processor according to the extended method according to one embodiment; and
Fig. 7 is an exemplary illustration of performance analysis of the extended method according to one embodiment.
DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS Fig. 1 is an exemplary block diagram of a wireless system 100 according to a first embodiment. The wireless system 100 includes a transmitter 110 including an array of transmit antennas 120, a receiver 130 including an array of receive antennas 140, and a controller or processor 150. Preferably, the number of transmit antennas np is equal to the number of receive antennas nR. In operation, the transmitter 110 transmits a signal vector via the array of transmit antennas 120. The signal vector is received by the receiver 130 via the array of receive antennas 140. For example, the transmitter 110 transmits different bit streams on different sub-channels on different transmit antennas 120. The number of subchannels preferably corresponds to the number of receive antennas n .
The signal vector received by the receive antennas 140 can be represented as:
rt = AsHtst + nt
where rt represents the received signal vector, As represents a power normalization constant, Ht represents the channel matrix, st represents the transmitted signal vector, nt represents a noise vector, and t represents a discrete time instant. For simplification, the time index may be ignored in the following descriptions.
The received signal vector values can further be represented as:
The noise vector can be assumed to be an nr dimensional complex Gaussian vector with zero mean and covariance matrix I. Noise is assumed to be independent at different time instants. A k-dimensional complex (real) Gaussian random vector with mean m and covariance matrix R can be denoted by:
CNk(m, R)(Nk(m, R)).
Each entry hj j in the channel matrix H represents the path gain between the i-th receive antenna and the j-th transmit antenna. In a rich scattering environment, all entries are modeled as independent identically distributed (iid) 0^(0,1) random variables. A chi- squared variable with K degrees of freedom denoted by ^ is defined as:
Σ ι=l «?
where a\ represents independent identically distributed N1 (0,1/2) random variables.
For further analysis, E[«] denotes expectation, t means transpose for real matrices where the transpose is Hermitian for complex matrices, and _L stands for "orthogonal to."
Generally, transmitted power is equally distributed between transmit antennas and the same signal constellation is used at all transmit antennas. The power normalization constant As normalizes the signal constellation to have unit average power. Since the noise has unit power in the channel model, the signal to noise ratio (SNR) is represented as:
SNR = nTAs
The transmitted power is generally fixed independent of the number of transmit antennas.
Generally, the channel matrix H changes over time due to channel fading. A quasi- static approximation of the fading channel is adopted. For example, the channel remains unchanged during a coherent channel period which lasts TCOh time instants. However, the channel changes independently from one coherence period to another. Also, H is generally known to the receiver 130, but not to the transmitter 110.
The ergodic capacity of the multiple antenna system 100 can be represented as:
C = E[C(H)] = E[log det(I + (SNR/nT)HHt)]
In particular, C(H) = log(l + SNR χ2 2\ ) for a k-order diversity (l,k) system
According to one embodiment, the processor performs baseline space-time processing method or baseline method for removing interference for nj = nR = n. In this method, the processor generates a decision of a desired sub-channel of the signal vector by nulling out the interference of a first set of sub-channels by multiplying the received signal vector by a unitary matrix generated from a QR decomposition of a channel matrix and the processor further generates the improved symbol decision by successively canceling channel interference due to a second set of sub-channels. The symbol represents transmitted bits, transmitted coded bits, or other useful transmitted information. For example, the symbol represents the decoded or coded bits of a digital transmission. A more detailed description of the baseline method is described below.
For initial definitions, a tight lower bound of C(H) for a n transmitter and n receiver, (n,n), system is found according to:
C(H) > ∑log(l + (SNi?/n)^ ', 2k k=\
While this baseline space-time processing analysis is calculated for a (n,n) system, it can also apply to a p transmitter and n
R receiver, (nτ,n ), system. The channel matrix H can be written as (hi, h , . . . , h
n) where h
j with 1 < j < n denotes the j-th column of H. Then, let H
j denote the linear space spanned by hi, h
2, . . . , h
j. Accordingly, Hj - Span(h
l5 h
2, . . . , hj). Fig. 2 is an exemplary flowchart 200 outlining the more detailed operation of the processor 150 according to the baseline method. In step 205, the operation begins. In step 210, the processor 150 performs a QR decomposition of the channel matrix H. hi particular, H=QG, where Q is a unitary matrix and G is an upper triangular matrix (the variable G is used instead of the classic R because R is already used as a variable in the process). Accordingly,
Where Q is a unitary matrix with the following properties:
Hj = Span(q1,...,qy) l ≤ j ≤ n HjΛ ± qj , l ≤ j ≤ n.
The entries in the upper triangular matrix G have particular distributions summarized as:
The signal vector received by the receive antennas 140 can be rewritten as:
Symbols from different antennas interfere with each other. Accordingly, in step 215, the processor 150 performs interference nulling in one step by multiplying r by Qt. This is a significant simplification over present systems that perform extensive processing to null out interfering sub-channels. In particular, present systems perform interference nulling by decoding sn by projecting r onto H„_ , thus avoiding interference from symbols sl .., sn. The signal due to sn is then subtracted from r and sn-1 is decoded by nulling out interference from s to sn-2. This processing proceeds until the signal is decoded.
According to the present invention, the processor 150 determines the symbol y=Qtr, where r is the received signal and Qt is the Hermitian transpose of Q. In other words, y=Gs+w, where s is the transmitted signal and w is the noise. Therefore, the interference nulling pattern is created directly by the unitary transformation. In step 220, the processor 150 decodes the symbol yn of the n transmit antenna according to:
Again, where As is the power normalization constant, gn>n is the upper triangular matrix, sn is the transmitted signal, and wn is the noise. The corresponding decision is denoted sn .
In the following steps, the processor 150 successively decodes the symbols from n-1 to the first transmit antenna. Accordingly, in step 225, the processor sets j=n-l. In step 230, the processor determines if j<l. If true, then the processor exits the operation in step 235. If false, the processor proceeds to step 240 where it subtracts interference from the next symbol according to:
which can be rewritten as:
n
■ AsSj, j + ∑ gjΛsk -sk)+ w, k=j+l
h step 245, the processor 150 decodes the symbol decision Sj based on yj . In step 250, the processor decrements j by one and returns to step 230. This process effectively creates n non-interfering one-dimensional sub-channels where the j sub-channel is described as:
yj,t = As§j ,tsj,t + wj,t^ ≤ j ≤ n
where
° j,j,ι • X: 2(«~7+l)
In particular, the jth sub-channel corresponds to n-j+1 order diversity.
Fig. 3 is an exemplary illustration of a space-time processing schematic for a six transmitter and six receiver (6, 6) system according to one baseline embodiment. Tx represents the transmitter index and D represents the diversity order associated with the corresponding sub-channel. The horizontal direction indicates time instants and the vertical direction indicates spatial sub-channels created by the baseline method. As shown, each subchannel is associated with a particular transmit antenna. The space-time processing is visualized by a two dimensional stack of rectangles where STj;t denotes the rectangle corresponding to the jth sub-channel at the tth time instant. This symbolizes a space-time processing unit aimed at retrieving information from the j* sub-channel at the tth time instant by subtracting interference from other transmitters and forming decisions Sj according to the particular coding scheme used. Similarly, a row of rectangles represent the processing with respect to the corresponding spatial sub-channel at different time instants and a column of rectangles represent the processing at the corresponding time instant for all spatial subchannels. Accordingly, the system equation for STμ can be written as
which shows that interference from the transmit antennas 120 corresponding to rectangles above the unit have been nulled out and that interference from transmit antennas 120
corresponding to rectangles below the unit must be cancelled by using decision feedbacks. Therefore, the baseline method provides a framework to construct various layered space-time structures.
Fig. 4 is an exemplary illustration of a schematic for processing similar to BLAST for a (6, 6) system. According to this processing, a layer is composed of six units in the diagonal direction where the arrow indicates the processing from left to right. Accordingly, each of the six diagonal unit in the schematic comprises a layer upon which codes can be applied. Different layers are processed from left to right for interference cancellation. The baseline method improves on BLAST at least in its simplicity and efficiency of operation. For performance analysis, a D order diversity channel is described as y=hx+n where ~ χ1D , n ~ CN1 (O, σ2 j and As normalizes the signal constellation to have unit power.
The average SΝR per diversity branch is denoted p = -^ A2 . The probability of bit error, Pe, of σ the binary phase shift keying (BPSK) is then determined according to:
where
The term ^T AsgJ k (sk —sk ) + Wj is the interference term due to imperfect decision
A=7+l feedbacks. In order to quantify the effect of interference, it is approximated as a Gaussian random variable, Wj . Pe is defined for S to be pk. That is, Pr(-?A ≠ sk) = pk . Given an equal probability of 1 and -1 under BPSK modulation, it is easy to verify that
E[s* -£ = 0
and
E s - sk = 4Pk
Therefore, the mean and variance of w . are given by:
and it is independent of s
k -s
k . Then, P
e is calculated for all sub-channels under a Gaussian approximation of interference. P
e is calculated by first calculating p of the nx sub-channel by using the equation for P
e with a diversity order of D=n
R-nτ+l and
As I . Second, p is calculated successively from nχ-1 to 1. For a given k , the interference power σ1 is calculated. Then, p is calculated by utilizing the equation for Pe
with D=nR-k+l and p ■ l2 Third, the overall probability of bit error, Pe, is obtained by
1 + σ2 averaging across all the sub-channels. Accordingly,
∑Pk P (SNR) = -^ —
Fig. 5 is an exemplary illustration of performance analysis of the baseline method. This performance analysis illustrates the lower bound for the probability of bit error. This is based on a simulation of baseline processing for a 6 transmitter and 6 receiver (6,6) system with TCO = 10,000 time instants. BPSK was used at the transmitter. The estimated performance was averaged over 10,000 channel coherence periods, or a total of 108 time instants. The nominal Pe performance, assuming perfect interference cancellation was computed using the above equations and served as a lower bound on the actual Pe performance. As described in the Background, since feedback from sub-channels with low diversity are used to create sub-channels with high diversity, imperfect decision feedbacks can compromise the performance of the sub-channels with a high diversity order. In Fig. 5, the nominal and the actual Pe is plotted for several sub-channels in the (6,6) system. In comparing nominal performance with actual simulation data, it is shown that imperfect decision feedbacks can affect sub-channels with large diversity order. Also, the system performance can be limited by the worst sub-channel. Thus, the worst sub-channel in a layered space-time structure can be a bottleneck in limiting system performance.
According to another embodiment, the processor 150 performs an extended space- time processing method or extended method.
As shown above, the strongest sub-channel with the largest diversity order is formed by successive decision feedbacks from the rest of the sub-channels. Accordingly, the strongest sub-channel has better performance than those with less diversity. Unfortunately, this fact is not utilized in a BLAST-type space-time processing structure. The present inventors have recognized that a direct way to improve system performance is to subtract the strongest sub-channel from the total received signal using the strongest sub-channel's decision. This is defined as a "loopback" process because it forms a feedback flow in a reverse order. Loopback operation can effectively remove the contribution of the transmit antenna corresponding to the strongest sub-channel.
A (6,6) system can be used to further explain the loopback procedure. First, decisions are generated for the sub-channels using the baseline operation, BLAST, or the like. The resulting diversity order for each respective sub-channel is 6, 5, 4, 3, 2, 1. From the
generated decisions, a signal due to, for example, the 1st transmit antenna is effectively reconstructed. This signal is subtracted from the received signal. The received signal is now effectively a 5 transmitter and 6 receiver (5,6) system. If the baseline operation is performed, the diversity order of all of the sub-channels from 2 to 6 is improved by 1 diversity order. These improved sub-channels produce improved decisions which thus assist in the decoding of the 1st sub-channel. Thus, the first iteration of a loopback operation results in a diversity order for each respective sub-channel being 6, 6, 5, 4, 3, 2. This accordingly improves the diversity order of the second sub-channel to full diversity. The loopback operation can be continued by subtracting successively improving sub-channels. For example, both the 1st and 2nd sub-channels are next subtracted, thus effectively forming a (4,6) system.
The depth of the loopback operation is defined as Iioopback- Therefore, sub-channels 1 to Iioop ack are used successively in the loopback cancellation. Since processing structure evolves after each loopback operation, the final structure has the property that the first Iloopback + 1 sub-channels have full diversity order, while the rest have 5 down to 1 + Iιoopback- For example, the baseline operation corresponds to Iioopback = 0. Also, the loopback processing can iterate at the final stage Ifmaι times without looping back to the next sub-channel.
In the extended method, the processor 150 generates an improved symbol decision of a desired sub-channel of the signal vector by first generating a baseline decision for the subchannel. This baseline decision can be generated by the baseline method, by the BLAST system, or by any other method useful for generating a decision for a sub-channel. Next, the processor 150 subtracts a contribution of a strongest sub-channel from the signal vector to generate a modified signal vector. Then, the processor 150 multiplies the modified signal vector by a unitary matrix generated from a QR decomposition of another channel matrix. Finally, the processor 150 successively cancels channel interference of the remaining sub- channels of the modified signal vector from a remaining sub-channel.
A more detailed description of the extended method is described below: Fig. 6 is an exemplary flowchart 600 outlining the operation of the processor 150 according to the extended method. In step 605, the operation begins. In step 610, the processor 150 generates the initial symbol decisions (JSJ , i?2 ••• sn J. The processor 150
can generate these symbol decisions by utilizing the baseline method, a method such as BLAST, or any like method. For example, the processor 150 decodes each sub-channel utilizing the baseline method for (nχ5 ΏR). The decision at the kth sub-channel is denoted sk .
In step 615, the processor 150 sets i=l to begin a loop for i=l to Iioopback to use sub-channels 1 to Iioopback successively in loopback cancellation. In step 620, the processor 150 determines if i > Iioopback- If so, the processor advances to step 660. If i < Iioopback, the processor 150 advances to step 625.
In step 625, the processor subtracts the signals from the sub-channels 1 to i to generate an improved received signal r according to:
= AsUisi + ∑AshJ(sj -sJ) + n
7=1
= A.Hls, + n
In step 630, the processor 150 decodes the received signal r to generate symbol decisions
Sj = [sM , • • • , sn j , where t stands for the transpose. The processor 150 can decode the received signal r by utilizing the baseline method, a method such as BLAST, or any like method.
In step 635, the processor 150 sets k=l to begin a loop for k=i to 1 to cancel interference from the other sub-channels. In step 640, the processor determines whether to exit the loop based on k<l. If k≥l, the processor advances to step 645. In step 645, the processor cancels the interference from the other sub-channels according to:
y
and then updates symbol decision s
kby decoding r
k . In step 650, the processor 150 decrements k and returns to step 640.
If k<l in step 640, the processor 150 increments i in step 665 and returns to step 620. If i>Iioo
Pback in step 620, the processor 150 advances to step 660 where it fixes
In step 665, the processor 150 performs the functions of steps 625-650 for I
fmai-1 iterations. In step 670, the operation ends.
Fig. 7 is an exemplary illustration of performance analysis of the extended method. As discussed above, the performance can be studied analytically by approximating the residual interference after the decision feedback as Gaussian random variables. Pe in each step of the extended algorithm is obtained by calculating the variance of residual interference.
The performance algorithm is given by setting p = [p • ••, p„τ f to be a vector of Pe's of all sub channels at each step of the extended algorithm with a specific Iioopback and If,naι. Then set
p = (p1 • ■ -, p )'. Given a SNR, \A \ = . Then, p is calculated in accordance with τ nτ the baseline structure. Next, a loop is set from i = 1 to Iι00pback- By Gaussian approximation of residual interference, the noise n is distributed as CN"R (O, σ2 J with
σa =l + 4 ,a ^
7=1
Then, pi is generated as described above with respect to the baseline operation. Next, from k = i to 1, pk is updated using the above equations with a diversity order of D = nR and
Fig. 7 illustrates the performance improvement with loopback cancellation in the extended method for a (6,6) system. As before, bits are sent using BPSK at each transmit antenna. L-maι is fixed to be 1 while varying the loopback depth Iioop ack from 1 to 5. The
performance of a baseline or BLAST-type algorithm is included for comparison. As illustrated, large performance gain is achieved by using loopback cancellation. For example, the system with full loopback at 4dB already achieves Pe of the BLAST system at 9dB, thus resulting in a 5dB savings. Performance gain increases as SNR increases, which projects more power savings at higher SNR's. Also, a few loopback cancellations can be sufficient.
As illustrated in Fig. 7, large performance improvements are achieved with only 1 or 2 levels of loopback.
While this invention has been described with specific embodiments thereof, it is evident that many alternatives, modifications, and variations will be apparent to those skilled in the art. For example, various elements and various steps of different aspects and embodiments may be combined with various elements and various steps of other aspects and embodiments. Accordingly, the preferred embodiments of the invention as set forth herein are intended to be illustrative, not limiting. Various changes may be made without departing from the spirit and scope of the invention.