1. Field
This disclosure relates generally to wireless communication systems and, more specifically, to techniques for frequency-domain joint detection in wireless communication systems.
2. Related Art
To increase capacity and performance of wireless communication systems that employ time division-synchronous code division multiple access (TD-SCDMA) techniques, a base station (BS) receiver may employ joint (multiple user) detection. Joint detection is similar to solving a least squares (LS) problem, which may represent a significant computational effort due to the amount of data that may be involved. In general, joint detection combines knowledge about all subscriber stations (SSs) that are active in one burst in a relatively large system of equations. This knowledge has included channel impulse responses (that have been estimated from training sequences), spreading codes, and received antenna samples. Typically, designers have attempted to develop algorithms that lower computational complexity associated with joint detection without significantly degrading joint detection performance. Traditionally, joint detection has been performed using time-domain approaches (in contrast to frequency-domain approaches), due to the lower complexity traditionally associated with time-domain approaches.
At least one known joint detection frequency-domain approach has attempted to reduce computational complexity associated with joint detection by employing an algorithm that is based on extending a system matrix of a least squares (LS) problem to a block circulant matrix that is then block diagonalized. As is known, a circulant matrix is a square matrix where each column has the same elements as the column to the left of it rotated down by one position. Due to the fact that Fourier vectors are eigenvectors of circulant matrix blocks, a system of equations whose defining matrix is circulant can be solved efficiently in the frequency-domain. The transformation to and from the frequency-domain can be done efficiently with fast Fourier transforms (FFTs), which can be extended to block-circulant systems. Similar to how circulant matrix blocks can be diagonalized by FFTs, block circulant matrix blocks can be block diagonalized by block FFTs.
According to this approach, a time-domain block circulant channel matrix has been converted to frequency-domain block diagonal channel matrix by using block FFTs. In a TD-CDMA system, K CDMA codes may be simultaneously active on the same frequency band and in the same time slot. For example, a TD-SCDMA system that is designed for eight simultaneous users may include eight (K=8) simultaneously active different CDMA codes. The different spreading codes facilitate signal separation at a receiver of a base station (BS). According to the required data rate, a given user might use several different CDMA codes and/or time slots. The transmission of one block of N data symbols can be modeled by a system of linear equations with a vector containing the received samples from all antennas, an unstructured matrix that contains knowledge about estimated channel impulse responses and spreading codes, and another vector that represents temporal and spatial noise. While joint detection in the frequency-domain may be performed with reduced computational complexity using FFTs, the computational complexity of the frequency-domain approach may still be relatively high and unacceptable for many applications.
The present invention is illustrated by way of example and is not intended to be limited by the accompanying figures, in which like references indicate similar elements. Elements in the figures are illustrated for simplicity and clarity and have not necessarily been drawn to scale.
In the following detailed description of exemplary embodiments of the invention, specific exemplary embodiments in which the invention may be practiced are described in sufficient detail to enable those of ordinary skill in the art to practice the invention, and it is to be understood that other embodiments may be utilized and that logical, architectural, programmatic, mechanical, electrical and other changes may be made without departing from the spirit or scope of the present invention. The following detailed description is, therefore, not to be taken in a limiting sense, and the scope of the present invention is defined only by the appended claims and their equivalents. In particular, although various embodiments are described below in conjunction with a base station, the techniques may be employed in a subscriber station, such as a cellular handset. It will be appreciated that the present invention is not limited to a cellular handset and may be embodied in various other mobile devices, e.g., personal digital assistants (PDAs), digital cameras, portable storage devices, audio players, computer systems, and portable gaming devices.
According to one embodiment of the present disclosure, a method of operating a wireless communication device includes receiving, at a first wireless communication device, respective signals transmitted from multiple second wireless communication devices. Respective channel matrix blocks are then generated. Finally, respective channel equalizer matrix blocks for each of the respective channel matrix blocks are generated. At least one of the respective channel equalizer matrix blocks is generated without performing a matrix inversion operation on an expression that includes an associated one of the respective channel matrix blocks.
According to another embodiment of the present disclosure, a wireless communication device includes a receiver and a processor. The receiver is configured to receive respective code division multiple access signals. The processor is in communication with the receiver and is configured to generate respective channel matrix blocks for each block of the received respective signals. The processor is also configured to generate respective channel equalizer matrix blocks for each of the respective channel matrix blocks. At least one of the respective channel equalizer matrix blocks is generated without performing a matrix inversion operation on an associated one of the respective channel matrix blocks.
According to a different aspect of the present disclosure, a method of operating a wireless communication device includes receiving, at a first wireless communication device, respective code division multiple access signals transmitted from multiple second wireless communication devices. Respective channel matrix blocks are then generated. Finally, respective channel equalizer matrix blocks for each of the respective channel matrix blocks are generated. At least one of the respective channel equalizer matrix blocks is generated without performing a matrix inversion operation on an expression that includes an associated one of the respective channel matrix blocks.
According to one or more aspects of the present disclosure, a zero-force (ZF) equalizer solution for a received time-domain signal in a time division-synchronous code division multiple access (TD-SCDMA) wireless communication system is given by:
d=(CHC)−1CHx′=F(K)H(ΛHΛ)−1ΛHF(P)x′
where C is the time-domain block channel matrix; CH is the Hermitian of the time-domain block channel matrix C; x′=L−1x; C=L−1A; Rn=LLH is the space-time covariance matrix; d is the decoded received time-domain signal; Λ=diagonal[F(p)C(:,1:K)] is frequency-domain block channel matrix; F(p)=FI(p), where F is a fast Fourier transform (FFT) matrix; and I(p) is the identity matrix. A minimum mean square error (MMSE) solution for a TD-SCDMA system is given by:
where Rs=σs2I is the signal covariance matrix. While the following analysis is based on the ZF equalizer solution, it should be appreciated that the joint detection techniques disclosed herein may be readily extended to the MMSE solution.
In general, a frequency-domain solution includes the following four steps: generate the time-domain channel matrix blocks C1 and C2; calculate the vector g (g=F(P)x′=(FI(P))x′), which is the FFT of the received time-domain signal x; calculate the frequency-domain diagonalized channel matrix blocks Λd (Λd=C1+WDdC2, where
and D is the total number of blocks (e.g., twenty-four blocks when a 24-point DFT is employed); pass the frequency-domain signal g through an equalizer to provide an equalized vector h (h=(ΛHΛ)−1ΛHg) with Λ=diag[Λ1, Λ2, . . . , Λ24]; and convert the frequency-domain equalized vector h back to time-domain signal d (d=F(K)Hh). In the exact conventional solution, twenty-four matrix inversion operations have been required to convert an expression (i.e., Ωd=(ΛdHΛd)−1) including each of the respective channel matrix blocks to channel equalizer matrix blocks for a 24-point DFT.
According to various aspects of the present disclosure, techniques are disclosed that reduce the computational complexity of generating channel equalizer matrix blocks (Ωd=(ΛdHΛd)−1) in the frequency-domain for joint detection in a wireless communication system, e.g., a time division-synchronous code division multiple access (TD-SCDMA) system. The disclosed techniques reduce joint detection computation complexity by replacing direct matrix inversion of, in this case twenty-four matrix blocks (see diagram 100 of
According to one embodiment of the present disclosure, joint detection computation complexity is reduced by only determining some channel equalizer matrix blocks by performing matrix inversion operations of expressions including associated channel matrix blocks (e.g., only inverting every second, third, or fourth expression in matrix interpolation algorithm). Interpolation is then utilized to determine values for elements in remaining ones of the channel equalizer matrix blocks that are not associated with matrix inversion operations. With reference to
The filter block 206 processes the data received from the block 204, according to the channel matrix blocks (Λd). The output of the filter block 206 is provided to an equalizer block 208. The matrix blocks (Λd) (provided by the block 214) are also provided to a process 218 that may be configured to convert the matrix blocks (Λd) to channel equalizer matrix blocks (Ωd) according to various embodiments of the present disclosure. In various embodiments, the process 218 employs fewer matrix inversions than conventional approaches for determining the channel equalizer matrix blocks (Ωd). According to the conventional approach, the channel equalizer matrix blocks (Ωd) are found by inverting twenty-four (i.e., d=1, . . . , 24) matrix blocks (Ωd=(ΛdHΛd)−1) that include the channel matrix blocks (Λd) multiplied by a Hermitian (ΛdH) of the matrix blocks. The channel equalizer matrix blocks (Ωd) are provided to the equalizer block 208, which equalizes the channels based on the channel equalizer matrix blocks (Ωd).
Outputs of the equalizer block 208 are provided to a matrix transpose and inverse DFT (IDFT) block 210 (included within a DFT engine), which converts the frequency-domain signal to a decoded time-domain signal d using an IDFT, e.g., an inverse fast Fourier transform (IFFT), and provides decoded time-domain signals (dd) for each subscriber station (SS) that is in communication with the BS.
According to another embodiment of the present disclosure, a matrix power expansion is performed to determine one of the channel equalizer matrix blocks, instead of directly calculating the inverse of twenty-four matrix blocks. A matrix power expansion for one of the channel equalizer matrix blocks may be performed as follows:
A first-order approximation for the channel equalizer matrix is given by:
Ωd≈A−1−xdA−1B(1)A−1−ydA−1B(2)A−1
A=C1HC1+C2HC2; B(1)=C1HC2+C2HC1; and B(2)=j(C1HC2−C2HC1). According to this approach, only one matrix (i.e., the A matrix) requires inversion to solve for all of the channel equalizer matrix blocks.
According to another aspect of the present disclosure, coefficient interpolation is employed on a first-order matrix power expansion of a channel equalizer matrix block. In this case, a first-order expansion of the channel equalizer matrix block may be given by:
Ωd(xd,yd)=(A+xdB(1)+ydB(2))−1≈Q00+xdQ10+ydQ01+ . . .
Q00, Q10, and Q01 are then estimated from three specific values of Ωd (xd,yd), namely:
where ρ is less than one and corresponds to a radius of a circle on which the points lie. In this case, Q00, Q10, and Q01 are determined from:
The above-described approach requires three matrix inversions, instead of twenty-four matrix inversions. The values of the twenty-four channel equalizer matrix blocks are then calculated based on Q00, Q10, and Q01, as follows:
Ωd(xd,yd)=Q00+xdQ10+ydQ01
To increase accuracy of the channel equalizer matrix blocks (at the expense of increased complexity), a second-order expansion may be employed. For second-order expansion of the channel equalizer matrix blocks, estimates of Q00, Q10, Q01, Q20, Q11, and Q02 can be obtained from six specific values of Ωd (xd,yd). A second-order expansion is given by:
In this case the six Q matrix blocks (i.e., Q00, Q10, Q01, Q20, Q11, and Q02) may be estimated by choosing five complex numbers z1, . . . z5 that sum to zero (it may be convenient to choose zn=ρ·exp(2πj(n−1)/5). With zn=xn+j·yn:
Θn≡Ω(xn,yn)=[Q00Q10Q01Q20Q11Q02]·rn
where rn≡[1xnynxn2xn·ynyn2]T. Since Ω(x+j·y) is complex analytic, then by the harmonic property:
Ω(0,0)≈1/5·Σn=1 . . . 5Ω(xn,yn)
which yields:
Q00=1/5·Σ·n=1 . . . 5Θn.
To find Q10, Q01, Q20, Q11, and Q02 an R matrix is formed as follows:
R≡[r1r2r3r4r5].
Then for each matrix row ‘a’ and each matrix column ‘b’ we have:
{[Θ1[a,b],Θ2[a,b],Θ3[a,b],Θ4[a,b],Θ5[a,b]]−Q00[a,b]}=[Q10[a,b],Q01[a,b],Q20[a,b],Q11[a,b],Q02[a,b]]·R
or
{[Θ1[a,b],Θ2[a,b],Θ3[a,b],Θ4[a,b],Θ5[a,b]]−Q00[a,b]}·R−1=[Q10[a,b],Q01[a,b],Q20[a,b],Q11[a,b],Q02[a,b]]
The system of equations provides a solution for all matrix entries for all five matrix blocks. It should be noted that ‘R’ is a fixed matrix that only needs to be inverted once.
With reference to
Moving to
Moving to
Based on employing a 24-point FFT engine, analysis indicates a complex multiply and accumulate (CMAC) reduction for an intra-cell based receiver of 17.6 percent, 19.8 percent, 19.8 percent, and 22 percent for matrix interpolation on every three matrix blocks, matrix interpolation on every four matrix blocks, first-order matrix power expansion, and first-order power expansion with coefficient interpolation, respectively. Analysis also indicates a CMAC reduction for an inter-cell based receiver of 24.8 percent, 28 percent, 34.3 percent, and 38.4 percent for matrix interpolation on every three matrix blocks, matrix interpolation on every four matrix blocks, first-order matrix power expansion, and first-order power expansion with coefficient interpolation, respectively.
With reference to
As is noted above, according to various embodiments of the present disclosure, the BTS 914 is configured to perform joint detection for wireless communication devices, such as the wireless devices 902. The wireless devices 902 communicate with a base station controller (BSC) 912 of a base station subsystem (BSS) 910, via the BTS 914, to receive or transmit voice, data, or both voice and data. The BSC 912 may, for example, be configured to schedule communications for the wireless devices 902. Alternatively, the BTS 914 may schedule communications for the wireless devices 902 in which the BTS 914 is in communication. In either case, a scheduler typically employs one or more processors (that execute a software system) to schedule communications and to perform joint detection. The BTS 914 may, for example, employ one or digital signal processors (DSPs) each of which include one or more engines for performing discrete Fourier transforms (DFTs) and inverse DFTs (IDFTs) to facilitate joint detection in the frequency-domain.
The BSC 912 is also in communication with a packet control unit (PCU) 916, which is in communication with a serving general packet radio service (GPRS) support node (SGSN) 922. The SGSN 922 is in communication with a gateway GPRS support node (GGSN) 924, both of which are included within a GPRS core network 920. The GGSN 924 provides access to computer(s) 926 coupled to Internet/intranet 928. In this manner, the wireless devices 902 may receive data from and/or transmit data to computers coupled to the Internet/intranet 928. For example, when the devices 902 include a camera, images may be transferred to a computer 926 coupled to the Internet/intranet 928 or to another one of the devices 902. The BSC 912 is also in communication with a mobile switching center/visitor location register (MSC/VLR) 934, which is in communication with a home location register (HLR), an authentication center (AUC), and an equipment identity register (EIR) 932. In a typical implementation, the MSC/VLR 934 and the HLR, AUC, and EIR 932 are located within a network and switching subsystem (NSS) 930, which may also perform scheduling for the system 900. The SGSN 922 may communicate directly with the HLR, AUC and EIR 932. As is also shown, the MSC/VLR 934 is in communication with a public switched telephone network (PSTN) 942, which facilitates communication between wireless devices 902 and land telephones 940. It should be appreciated that other types of wireless systems, having different configurations, may implement various aspects of the joint detection techniques disclosed herein.
Accordingly, a number of techniques have been disclosed herein that generally reduce the complexity of performing joint detection in a TD-SCDMA system. It is contemplated that the joint detection techniques described herein may be advantageously employed in various wireless communication systems that employ a code division multiple access (CDMA) scheme on at least some channels.
As used herein, a software system can include one or more objects, agents, threads, subroutines, separate software applications, two or more lines of code or other suitable software structures operating in one or more separate software applications, on one or more different processors, or other suitable software architectures.
As will be appreciated, the processes in preferred embodiments of the present invention may be implemented using any combination of computer programming software, firmware or hardware. For example, joint detection software that implements the processes 600, 700, and 800 of
Although the invention is described herein with reference to specific embodiments, various modifications and changes can be made without departing from the scope of the present invention as set forth in the claims below. For example, the preceding techniques disclosed herein are generally broadly applicable to transmitters and receivers, irrespective of location, in wireless communication systems. Accordingly, the specification and figures are to be regarded in an illustrative rather than a restrictive sense, and all such modifications are intended to be included with the scope of the present invention. Any benefits, advantages, or solution to problems that are described herein with regard to specific embodiments are not intended to be construed as a critical, required, or essential feature or element of any or all the claims.
Unless stated otherwise, terms such as “first” and “second” are used to arbitrarily distinguish between the elements such terms describe. Thus, these terms are not necessarily intended to indicate temporal or other prioritization of such elements.
Number | Name | Date | Kind |
---|---|---|---|
20030021359 | He | Jan 2003 | A1 |
20060126753 | Jwa et al. | Jun 2006 | A1 |
20070160029 | Luo et al. | Jul 2007 | A1 |
20080151985 | Chin et al. | Jun 2008 | A1 |
Number | Date | Country | |
---|---|---|---|
20090109947 A1 | Apr 2009 | US |