This invention generally relates to wireless code division multiple access (CDMA) communication systems. In particular, the invention relates to scaling using gain factors during data detection in such systems.
In wireless CDMA communication systems, multiple communications are transmitted over a shared spectrum. The individual communications are distinguished by a code used to transmit each communication. In frequency division duplex (FDD) CDMA systems, uplink and downlink communications are separated by frequency spectrum. In time division duplex (TDD) CDMA or time division synchronous code division multiple access (TD-SCDMA) systems, uplink and downlink communications are separated by time.
One approach to receive communications in such systems is referred to as single user detection (SUD). In general, SUD is typically performed in a two part process. First, the received signal is equalized to compensate for the response of the wireless channel. Second, the equalized signal is despread using the codes of the user or users. SUD is typically utilized when all communications experience a same channel response. The preferred uses for the SUD is in the downlink and in the uplink when a single user monopolizes a frequency spectrum (in FDD) or time slot (in TDD or TD-SCDMA). Another approach to receive communications in such systems is referred to as multi-user detection (MUD). In MUD, the data from all users' data is estimated simultaneously.
In such systems, the individual communications, transmitted simultaneously, may have varying transmission power levels. To illustrate, in the uplink, a single User Equipment (UE) may transmit multiple coded composite transport channels (CCTrCHs). Each CCTrCH may have a differing quality of service (QOS), requiring a different transmit power level to reach the required QOS. In the downlink, each UE in addition to each individual UE's CCTrCH may require a differing transmission power level.
Due to the varying power levels between communications, the orthogonality between the received codes is degraded, reducing the performance of the data detection. Accordingly, it is desirable to have alternate approaches to data detection.
Data is estimated from a received vector comprising a plurality of communications. A received wireless signal is converted to a baseband signal, which is sampled to produce a received vector. Channel responses are estimated for the received communications. The noise variance is estimated. The noise variance is scaled by a scaling factor. Samples are processed using the estimated channel responses and the scaled noise variance to produce a spread data vector. The spread data vector is despread to recover the data of the received wireless signal.
After passing through the wireless air interface 30, an antenna 32 or antenna array at the base station 22 receives the communications. The received communications are converted to a baseband signal, such as by a demodulation device 34. A sampling device 36 samples the baseband signal, such as at the chip rate or a multiple of the chip rate, to produce a received vector r. A channel estimation device 38 estimates the channel responses for the received communications, as a channel response matrix H. All of the H matrixes in the following equations are presumed to be normalized. The channel estimation is typically performed using a reference signal, such as the midambles of communication bursts or a pilot code. The channel estimation device 38 also preferably estimates a variance of the noise σ2, although the estimated noise variance may be derived from other sources. The noise variance σ2 is scaled by a scaling device, such as by a factor of 1/P. Preferred derivations for the value of P are explained in detail later. Although the preferred embodiment scales the noise variance σ2, the channel response matrix H may be scaled by a scaling device 61 instead of the noise variance, producing H′, as shown in
A channel equalizer 42 uses the channel response matrix H and the scaled noise variance σ2/P to produce a spread data vector s. Preferably, the spread data vector s is scaled by a scaling device 64, such as by a factor 1/P ΛH, although this scaling device 64 may not be used. (·)H is the conjugate transpose operation. When the scaling device 64 is not used, the spread data vector s is passed directly to the despreader 46. ΛH is a diagonal matrix, which preferred derivations are explained in detail later. The despreader 46 despreads the scaled spread data or spread data using the codes C of the communications to recover the data d.
In
At a UE 20, an antenna 28 or antenna array receives the communications sent through the wireless air interface 30. The received communications are demodulated by a demodulator 54 to baseband. A sampling device 56 samples the baseband signal to produce a received vector r. A channel estimation device 58 estimates the channel response matrix H and the noise variance σ2. The noise variance is scaled by a scaling device 60, such as by a factor of 1/P. Alternately, as shown in
Discussion of preferred algorithms for scaling in data detection are explained in conjunction with a preferred wideband CDMA (WCDMA) TDD system, although the invention is applicable to other systems, such as FDD/WCDMA, TD-SCDMA, CDMA 2000 as well as others.
The received vector r can be modeled per Equation 1.
r=AΛd+n Equation 1
A is the normalized symbol response matrix. Λ is the signal amplitude gain matrix and is a diagonal matrix. n is the noise vector.
When K codes are transmitted at the same time (within the same timeslot), A=[A1, . . . , AK] is the normalized symbol response matrix of the K codes. Λ is preferably of size K·Ns. Ns is the number of data symbols in a data field. d=[d1, . . . , dK] is the data sequence carried by the K codes. Λ is per Equation 2.
Each amplitude gain sub-matrix for a kth code of the K codes, Λk, is a diagonal matrix of size of Ns with all the diagonals equal to the signal amplitudes, gk, of the kth code, per Equation 3.
To estimate data {circumflex over (d)} in joint detection, a minimum mean square error (MMSE) approach per Equation 4 and 5 may be used.
{circumflex over (d)}=Λ−1(AHA+Λ−2σ2I)−1AHr Equation 4
{circumflex over (d)}=Λ−1AH(AAH+Λ−2σ2I)−1r Equation 5
{circumflex over (d)} is the estimated data vector. Equation 4 and 5 are interchangeable by the matrix inversion lemma.
When all the communications pass through the same propagation channel H, the symbol response matrix A is per Equation 6.
A=[A1, . . . , AK]=[HC1, . . . , HCK]=H[C1, . . . , CK]=HC Equation 6
H is preferably of size Q·Ns+W−1 by Q·Ns. Ns is the number of data symbols per data block that the data detection is performed. Q is the spreading factor of the communications and W is the length of the delay spread. C is preferably of size Q·Ns by K·Ns.
The received signal for a common channel can be modeled per Equation 7.
r=HCΛd+n Equation 7
Applying a MMSE solution to determine d is per Equation 8.
{circumflex over (d)}=ΛHCHHH(HMHH+σ2I)−1r Equation 8
The matrix M is the code correlation matrix per Equation 9.
M=CGCH Equation 9
M is preferably of size Q·Ns. The matrix G is the code power matrix, per Equation 10.
The code correlation matrix M is a diagonal dominant matrix with all the diagonals having a same value, such as P. One value for P is the total power of all bursts per Equation 11.
In general, P is referred to as the code power scale factor. By ignoring the edge effects of the non-diagonal part of the matrix, an approximation for M, {circumflex over (M)}, is per Equation 12.
{circumflex over (M)}=P·I Equation 12
By substituting Equation 12 into Equation 8, Equation 13 or 14 results.
Equation 13 and 14 are interchangeable by the matrix inversion lemma.
Equation 14 can be broken down in two stages per Equations 15 and 16.
In Equation 15, the channel equalization is performed. Preferably, the scaling in Equation 15 is performed by scaling devices 40, 60. Preferably, to reduce the complexity in solving equation 15 by the channel equalizer 42, 62, an approximate Cholesky or fast Fourier transform based solution is used, although other techniques may be used.
In Equation 16, the despreading is performed. The scaling,
during or prior to despreading may or may not be performed. Such scaling tends to improve the decoding performance of turbo-encoded signals.
By careful selection of the code power scaling factor, the overall performance of the receiver can be improved. The performance of the receiver can be modeled by Equation 17.
ŝ=s+If+w Equation 17
If represents the residual interference after equalization due to the incomplete equalization of channel distortion that results in the inter-symbol interference (ISI) and multiple access interference (MAI). ŝest of Equation 18 represents the impact of If on the estimated spread data.
w represents the noise after equalization and is per Equation 19.
The interference If and noise w cannot be reduced simultaneously. Decreasing one tends to increase the other. For a large P, the residual interference is reduced, but the noise is increased. For a small P, the noise is reduced but at the cost of increased residual interference.
Two preferred embodiments for power scaling, determining P, are total power scaling and selective scaling. In total power scaling, P is determined per Equation 20.
M is the number of UEs and N is the number of CCTrCHs per UE. Km,n is the total number of codes in the nth CCTrCH of the mth UE and gm,n is the gain factor for the nth CCTrCH of the mth UE.
Total power scaling tends to optimize performance over all connections equally and not to optimize any connection over the others. To illustrate, the code power matrix G is per Equation 21.
Λm,n2 is the code power sub-matrix corresponding to the nth CCTrCH of the mth UE. The code power gm,n2 in the matrix G and in G's sub-matrix can be approximated by one single common power in a least square error approach by minimizing Equation 22.
The solution of least-square-error-power is the average power of all codes per Equation 23.
K is the total number of codes transmitted in the time period of interest for the system and is per Equation 24.
The code power scale factor is determined by Equation 25.
gavg2 is the average code power. A measure of the mismatch between code powers Δm,n is per Equation 26.
Δm,n=|gm,n2−gavg2| Equation 26
As illustrated, total power scaling is optimal over all connections by minimizing the code power mismatch over all connections.
In selective code power scaling, P is determined to optimize a particular UE connection. To optimize a jth UE's connection, Equation 27 is used.
P=αjPT Equation 27
αj is a factor for the jth UE that is based on the interference and noise level. αj should be adaptively adjusted based on the interference and noise level for optimum performance of the data detection. Two preferred equations for deriving αj are per Equations 28 and 29.
Kj is the total number of codes carried by that jth UE. Ki,j is the number of codes for the ith CCTrCH of the jth UE. gi,j is the gain factor for the ith CCTrCH for the jth UE. I is the number of CCTrCHs of the UE.
Selective code power scaling may also be used to optimize a particular CCTrCH of a UE. To optimize the ith CCTrCH of the jth UE, Equation 30 is used.
P=αi,jPT Equation 30
αi,j is a factor for the ith CCTrCH of the jth UE.
Selective code power scaling may also be used to optimize a particular code of a particular CCTrCH of a UE. To optimize the nth code of the ith CCTrCH of the jth UE, Equation 31 is used.
P=αn,i,jPT Equation 31
αn,i,j is a factor for the nth code of the ith CCTrCH of the jth UE. Two preferred equations for determining αn,i,j are Equations 32 and 33.
Two special cases of selective code power scaling are maximum code power and minimum code power scaling. In maximum code power scaling, the maximum code power is used for the scaling. Maximum code power scaling is most applicable when the over-scaling of code power degrades less than the under-scaling of code power. In minimum code power scaling, the minimum code power is used for scaling. Minimum code power scaling is most applicable when the under-scaling of code power degrades less than the over-scaling of code power.
In some cases P is not necessarily determined from the gain factors. For instance, when a common midamble is used in a downlink allocation, the estimated channel response, H′, has the total power information. Accordingly, the total power is embedded in the estimated channel response at the output of the channel estimator, i.e. H′=H·√{square root over (P)}. In this alternative, gain scaling is not required and steps 78 and 80 of
In
In an alternate approach as shown in
Alternately, the gain factor scaling can be performed in conjunction with the multiuser detection. For such gain scaling, the multiuser detection device 108, 120 performs the gain factor scaling.
This application is a continuation of U.S. patent application Ser. No. 11/962,371, filed Dec. 21, 2007, which is a continuation of U.S. patent application Ser. No. 11/408,411, filed Apr. 21, 2006, which issued as U.S. Pat. No. 7,313,172 on Dec. 25, 2007, which is a continuation of U.S. patent application Ser. No. 11/175,662, filed Jul. 6, 2005, which issued as U.S. Pat. No. 7,042,929 on May 9, 2006, which is a continuation of U.S. patent application Ser. No. 10/327,299, filed Dec. 20, 2002, which issued as U.S. Pat. No. 6,928,104 on Aug. 9, 2005, which in turn claims priority from U.S. Provisional Application No. 60/396,823, filed Jul. 18, 2002, which are incorporated by reference as if fully set forth herein.
Number | Name | Date | Kind |
---|---|---|---|
5583886 | Rohani et al. | Dec 1996 | A |
5619525 | Wiedeman et al. | Apr 1997 | A |
5737327 | Ling et al. | Apr 1998 | A |
5991636 | Won et al. | Nov 1999 | A |
6085104 | Kowalski et al. | Jul 2000 | A |
6151358 | Lee et al. | Nov 2000 | A |
6175588 | Visotsky et al. | Jan 2001 | B1 |
6356639 | Ishito et al. | Mar 2002 | B1 |
6590932 | Hui et al. | Jul 2003 | B1 |
6801565 | Bottomley et al. | Oct 2004 | B1 |
6928104 | Pan et al. | Aug 2005 | B2 |
6952561 | Kumar et al. | Oct 2005 | B1 |
20020057735 | Piirainen | May 2002 | A1 |
20030099306 | Nilsson et al. | May 2003 | A1 |
20030099309 | Fimoff et al. | May 2003 | A1 |
20060121946 | Walton et al. | Jun 2006 | A1 |
20100166100 | Medvedev et al. | Jul 2010 | A1 |
Number | Date | Country |
---|---|---|
2311702 | Oct 1997 | GB |
2362523 | Nov 2001 | GB |
H11-266232 | Sep 1999 | JP |
2001-313588 | Nov 2001 | JP |
10-0301284 | Jun 2001 | KR |
0182511 | Nov 2001 | WO |
03025513 | Mar 2003 | WO |
Number | Date | Country | |
---|---|---|---|
20090316762 A1 | Dec 2009 | US |
Number | Date | Country | |
---|---|---|---|
60396823 | Jul 2002 | US |
Number | Date | Country | |
---|---|---|---|
Parent | 11962371 | Dec 2007 | US |
Child | 12549912 | US | |
Parent | 11408411 | Apr 2006 | US |
Child | 11962371 | US | |
Parent | 11175662 | Jul 2005 | US |
Child | 11408411 | US | |
Parent | 10327299 | Dec 2002 | US |
Child | 11175662 | US |