Derivation of eigenvectors for spatial processing in MIMO communication systems

[0001] The present Application for Patent claims priority to Provisional Application No. 60/432,760 entitled Derivation of Eigenvectors for Spatial Processing in MIMO Communication Systems filed Dec. 11, 2002, and assigned to the assignee hereof and hereby expressly incorporated by reference herein.

[0002] I. Field

[0003] The present invention relates generally to data communication, and more specifically to techniques for deriving eigenvectors based on steered reference and used for spatial processing in multiple-input multiple-output (MIMO) communication systems.

[0004] II. Background

[0005] A MIMO system employs multiple (N_T) transmit antennas and multiple (N_R) receive antennas for data transmission. A MIMO channel formed by the N_T transmit and N_R receive antennas may be decomposed into N_S independent or spatial channels, where N_SminN_T, N_R). Each of the N_S independent channels corresponds to a dimension. The MIMO system can provide improved performance (e.g., increased transmission capacity and/or greater reliability) if the additional dimensionalities created by the multiple transmit and receive antennas are effectively utilized.

[0006] In a wireless communication system, data to be transmitted is typically processed (e.g., coded and modulated) and then upconverted onto a radio frequency (RF) carrier signal to generate an RF modulated signal that is more suitable for transmission over a wireless channel. For a wireless MIMO system, up to N_T RF modulated signals may be generated and transmitted simultaneously from the N_T transmit antennas. The transmitted RF modulated signals may reach the N_R receive antennas via a number of propagation paths in the wireless channel. The characteristics of the propagation paths typically vary over time due to various factors such as, for example, fading, multipath, and external interference. Consequently, the RF modulated signals may experience different channel conditions (e.g., different fading and multipath effects) and may be associated with different complex gains and signal-to-noise ratios (SNRs).

[0007] To achieve high performance, it is often necessary to estimate the response of the wireless channel between the transmitter and the receiver. For a MINO system, the channel response may be characterized by a channel response matrix H, which includes N_TN_R complex gain values for N_TN_R different transmit/receive antenna pairs (i.e., one complex gain for each of the N_T transmit antennas and each of the N_R receive antennas). Channel estimation is normally performed by transmitting a pilot (i.e., a reference signal) from the transmitter to the receiver. The pilot is typically generated based on known pilot symbols and processed in a known manner (i.e., known a priori by the receiver). The receiver can then estimate the channel gains as the ratio of the received pilot symbols over the known pilot symbols.

[0008] The channel response estimate may be needed by the transmitter to perform spatial processing for data transmission. The channel response estimate may also be needed by the receiver to perform spatial processing (or matched filtering) on the received signals to recover the transmitted data. Spatial processing needs to be performed by the receiver and is typically also performed by the transmitter to utilize the N_S independent channels of the MIMO channel.

[0009] For a MIMO system, a relatively large amount of system resources may be needed to transmit the pilot from the N_T transmit antennas such that a sufficiently accurate estimate of the channel response can be obtained by the receiver in the presence of noise and interference. Moreover, extensive computation is normally needed to process the channel gains to obtain eigenvectors needed for spatial processing. In particular, the receiver is typically required to process the channel gains to derive a first set of eigenvectors used for spatial processing for data reception on one link and may further be required to derive a second set of eigenvectors used for spatial processing for data transmission on the other link. The derivation of the eigenvectors and the spatial processing for data transmission and reception are described below. The second set of eigenvectors typically needs to be sent back to the transmitter for its use. As can be seen, a large amount of resources may be needed to support spatial processing at the transmitter and receiver.

[0010] There is therefore a need in the art for techniques to more efficiently derive eigenvectors used for spatial processing in MIMO systems.

[0011] Techniques are provided herein for deriving eigenvectors based on steered reference and used for spatial processing for data reception and transmission. A steered reference is a pilot transmission on only one spatial channel or eigenmode of a MIMO channel for a given symbol period, which is achieved by performing spatial processing with a steering vector for that eigenmode, as described below. The steered reference is used by a receiver to derive estimates of both a diagonal matrix of singular values and a unitary matrix U of left eigenvectors of the channel response matrix H, without having to estimate the MIMO channel response or perform singular value decomposition of H.

[0012] The estimates of and U may be used for matched filtering of data transmission received via a first link (e.g., the uplink). For a time division duplex (TDD) system, which is characterized by downlink and uplink channel responses that are reciprocal of one another, the estimate of U may also be used for spatial processing of data transmission on a second link (e.g., the downlink).

[0013] In another aspect, a matrix with orthogonal columns is derived based on the estimates of and U. The orthogonalization of the columns of may be achieved by various techniques such as QR factorization, minimum square error computation, and polar decomposition, all of which are described below. An orthogonal matched filter matrix {tilde over (M)} may then be derived based on the matrix and the estimate of . The matrix {tilde over (M)} may be used for matched filtering for the first link, and the matrix may be used for spatial processing for the second link.

[0014] Various aspects and embodiments of the invention are described in further detail below.

[0015] The features, nature, and advantages of the present invention will become more apparent from the detailed description set forth below when taken in conjunction with the drawings in which like reference characters identify correspondingly throughout and wherein:

[0016] FIG. 1 shows a flow diagram of a process for deriving an orthogonal matched filter matrix {tilde over (M)} based on a steered reference;

[0017] FIG. 2 shows a wireless communication system;

[0018] FIG. 3 shows a frame structure for a TDD MIMO-OFDM system;

[0019] FIG. 4 shows transmission of steered reference and data on the downlink and uplink for an exemplary transmission scheme;

[0020] FIG. 5 shows a block diagram of an access point and a user terminal; and

[0021] FIG. 6 shows a block diagram of the spatial processing performed by the access point and user terminal for data transmission on the downlink and uplink.

[0022] The word exemplary is used herein to mean serving as an example, instance, or illustration. Any embodiment or design described herein as exemplary is not necessarily to be construed as preferred or advantageous over other embodiments or designs.

[0023] The techniques described herein for deriving eigenvectors may be used for various MIMO communication systems. For example, these techniques may be used for single-carrier MIMO systems as well as multi-carrier MIMO systems. For clarity, these techniques are described below for a single-carrier MIMO system.

[0024] The model for a single-carrier MIMO system may be expressed as:

[0025] where x is a transmit vector with N_T entries for the symbols sent from the N_T transmit antennas (i.e., xx₁ x₂ . . . x_NT^T);

[0026] r is a receive vector with N_R entries for the symbols received via the N_R receive antennas (i.e., rr₁ r₂ . . . r_NR^T);

[0027] H is an (N_RN_T) channel response matrix;

[0028] n is a vector of additive white Gaussian noise (AWGN); and

[0029] ^T denotes the transpose.

[0030] The noise vector n is assumed to have components with zero mean and a covariance matrix of _n²I, where I is the identity matrix and ² is the noise variance.

[0031] The channel response matrix H may be expressed as: $\begin{array}{cc}\underset{\_}{H}=\left[\begin{array}{cccc}{h}_{1,1}& h_{1,2}& & h_{1,N_T}\\ h_{2,1}& h_{2,2}& & h_{2,N_T}\\ & & & \\ h_{N_R,1}& h_{N_R,2}& & h_{N_R,N_T}\end{array}\right],& \mathrm{Eq}\left(2\right)\end{array}$

[0032] where entry h_i,j, for i1 . . . N_R and j1 . . . N_T, is the coupling (i.e., complex gain) between the j-th transmit antenna and the i-th receive antenna. For simplicity, the channel response is assumed to be flat across the entire system bandwidth, and the channel response for each transmit/receive antenna pair can be represented by a single complex value h_i,j. Also for simplicity, the following description assumes that N_RN_T, the channel response matrix H has full rank, and N_SN_TN_R.

[0033] The channel response matrix H may be diagonalized to obtain the N_T independent channels, which are also referred to as spatial channels or eigenmodes. This diagonalization may be achieved by performing either singular value decomposition of the channel response matrix H or eigenvalue decomposition of the correlation matrix of H, which is H^H H, where ^H denotes the conjugate transpose. For clarity, singular value decomposition is used for the following description.

[0034] The singular value decomposition of the channel response matrix H may be expressed as:

[0035] where U is an (N_RN_R) unitary matrix whose columns are left eigenvectors of H;

[0036] is an (N_RN_T) diagonal matrix of singular values of H, which is

[0037] V is an (N_TN_T) unitary matrix whose columns are right eigenvectors of H.

[0038] A unitary matrix M is characterized by the property M^HMI, which means that the columns of the unitary matrix are orthogonal to one another and the rows of the matrix are also orthogonal to one another. The columns of the matrix V are also referred to as steering vectors. Singular value decomposition is described in further detail by Gilbert Strang in a book entitled Linear Algebra and Its Applications, Second Edition, Academic Press, 1980.

[0039] Spatial processing may be performed by both the transmitter and the receiver to transmit data on the N_T spatial channels of the MIMO channel. The spatial processing at the transmitter may be expressed as:

[0040] where s is a data vector with up to N_T non-zero entries for data symbols to be transmitted on the N_T spatial channels. The transmit vector x is further processed and then transmitted over the MIMO channel to the receiver.

[0041] The received transmission at the receiver may be expressed as:

[0042] where all the terms are defined above.

[0043] The spatial processing at the receiver to recover the data vector s may be expressed as:

[0044] where s is the data vector;

[0045] is an estimate of the data vector s;

[0046] M is an (N_TN_R) matched filter matrix, which is M^TU^H;

[0047] G is an (N_TN_T) scaling matrix, which is

[0048] {circumflex over (n)} is the post-processed noise, which is {circumflex over (n)}G^TU^Hn.

[0049] The spatial processing by the receiver is often referred to as matched filtering. Since M^TU^H and since the columns of U are left eigenvectors of H. the columns of M^T are conjugated left eigenvectors of H scaled by the singular values in .

[0050] As shown in equation (6), the receiver needs good estimates of the matrices and U in order to perform the matched filtering to recover the data vector s. The matrices and U may be obtained by transmitting a pilot from the transmitter to the receiver. The receiver can then estimate the channel response matrix H based on the received pilot and perform the singular value decomposition of this estimate, as shown in equation (3), to obtain the matrices and U. However, as noted above, a large amount of resources may be needed to transmit this pilot and to perform the singular value decomposition.

[0051] I. Steered Reference

[0052] In an aspect, a steered reference is transmitted by the transmitter and used by the receiver to derive estimates of the matrices and U, which are needed for matched filtering. The steered reference is a pilot transmission on only one spatial channel or eigenmode for a given symbol period, which is achieved by performing spatial processing with a steering vector for that eigenmode. The receiver can then estimate the matrices and U based on the steered reference, without having to estimate the MINO channel response or perform the singular value decomposition.

[0053] A steered reference sent by the transmitter may be expressed as:

[0054] where x_sr,m is the transmit vector for the steered reference for the m-th eigenmode;

[0055] v_m is the right eigenvector of H for the m-th eigenmode; and

[0056] p is a pilot symbol transmitted for the steered reference.

[0057] The eigenvector v_m is the m-th column of the matrix V, where Vv₁ v₂ . . . v_NT.

[0058] The received steered reference at the receiver may be expressed as:

[0059] where r_sr,m is the receive vector for the steered reference for the m-th eigenmode; and

[0060] _m is the singular value for the m-th eigenmode.

[0061] As shown in equation (8), at the receiver, the received steered reference in the absence of noise is equal to u_mmp, which is the known pilot symbol p transformed by u_mm. The eigenvector u_m is the m-th column of the matrix U, and the singular value _m is the m-th diagonal element of the matrix . The receiver can thus obtain an estimate of u_mm based on the steered reference sent by the transmitter.

[0062] Various techniques may be used to process the received steered reference to obtain estimates of u_m and _m. In one embodiment, to obtain an estimate of u_mm, the receive vector r_sr,m for the steered reference sent on the m-th eigenmode is first multiplied with the complex conjugate of the pilot symbol, p*. The result may then be integrated over multiple steered reference symbols received for each eigenmode m to obtain the estimate of u_mm. A row vector {circumflex over (m)}_m may be defined to be equal to the conjugate transpose of the estimate of u_mm (i.e., {circumflex over (m)}_m{circumflex over ()}_mm^H) Each of the N_R entries of the vector {circumflex over (m)}_m is obtained based on a corresponding one of the N_R entries of the vector r_sr,m.

[0063] The row vector {circumflex over (m)}_m for the m-th eigenmode includes estimates of both u_m and _m, and may thus be referred to as a scaled vector. Since eigenvectors have unit power, the singular value _m may be estimated based on the received power of the steered reference, which can be measured for each eigenmode. In particular, the singular value estimate {circumflex over ()}_m may be set equal to the square root of the power for the vector r_sr,m, divided by the magnitude of the pilot symbol p. The vector {circumflex over (m)}_m may be scaled by 1/{circumflex over ()}_m to obtain the eigenvector _m.

[0064] In another embodiment, a minimum mean square error (MMSE) technique is used to obtain an estimate of u_m based on the receive vector r_sr,m for the steered reference. Since the pilot symbol p is known, the receiver can derive an estimate of u_m such that the mean square error between the recovered pilot symbol {circumflex over (p)} (which is obtained after performing matched filtering on the receive vector r_sr,m) and the transmitted pilot symbol p is minimized.

[0065] The steered reference is transmitted for one eigenmode at a time (i.e., one eigenmode for each symbol period of steered reference transmission). The steered reference for all N_T eigenmodes may be transmitted in various manners. In one embodiment, the steered reference is transmitted for one eigenmode for each frame, where a frame is an interval of data transmission for the system and is defined to be of a particular time duration (e.g., 2 msec). For this embodiment, the steered reference for multiple eigenmodes may be transmitted in multiple frames. In another embodiment, the steered reference is transmitted for multiple eigenmodes within one frame. This may be achieved by cycling through the N_T eigenmodes in N_T symbol periods. For both embodiments, the n-th steered reference symbol may be expressed as:

[0066] where n is an index for either symbol period or frame number and L is the number of steered reference symbols to be transmitted. Multiple steered reference symbols may be transmitted for each eigenmode m to allow the receiver to obtain more accurate estimate of u_mm.

[0067] The receiver is able to obtain the row vector {circumflex over (m)}_m for each of the N_T eigenmodes based on the received steered reference for that eigenmode. The row vectors {circumflex over (m)}_m for all N_T eigenmodes may be used to form an initial matched filter matrix {circumflex over (M)}, where {circumflex over (M)}{circumflex over (m)}₁ {circumflex over (m)}₂ . . . {circumflex over (m)}_NT^T and {circumflex over (M)}{circumflex over ()}^TH. The matrix {circumflex over (M)} may be used for matched filtering by the receiver, as shown in equation (6), to recover the transmitted data vector s.

[0068] The steered reference is sent for one eigenmode at a time and may be used by the receiver to obtain the matched filter vector {circumflex over (m)}_m for that eigenmode. Since the N_T matched filter vectors {circumflex over (m)}_m of the matrix {circumflex over (M)} are obtained individually and over different symbol periods, and due to noise and other sources of degradation in the wireless channel, the N_T vectors {circumflex over (m)}_m of the matrix {circumflex over (M)} are not likely to be orthogonal to one another. If the N_T vectors {circumflex over (m)}_m are thereafter used for matched filtering of a received data transmission, then any errors in orthogonality among these vectors would result in cross-talk between the individual symbol streams sent on the N_T eigenmodes. The cross-talk may degrade performance.

[0069] II. Eigenvector Orthogonalization

[0070] In another aspect, to improve performance, an enhanced matched filter matrix {tilde over (M)} is derived based on the steered reference and has row vectors that are forced to be orthogonal to one other. The orthogonalization of the row vectors of {tilde over (M)} may be achieved by various techniques such as QR factorization, minimum square error computation, and polar decomposition. All of these orthogonalization techniques are described in detail below. Other orthogonalization techniques may also be used and are within the scope of the invention.

[0071] 1. QR Factorization

[0072] QR factorization decomposes the transposed initial matched filter matrix, {circumflex over (M)}^T, into an orthogonal matrix Q_F and an upper triangle matrix R_F. The matrix Q_F forms an orthogonal basis for the columns of the matrix {circumflex over (M)}^T (i.e., the rows of {circumflex over (M)}), and the diagonal elements of the matrix R_F give the length of the components of the columns of {circumflex over (M)}^T in the directions of the respective columns of Q_F. The matrices Q_F and R_F may be used to derive an enhanced matched filter matrix {tilde over (M)}_F.

[0073] The QR factorization may be performed by various methods, including a Gram-Schmidt procedure, a householder transformation, and so on. The Gram-Schmidt procedure is recursive and may be numerically unstable. Various variants of the Gram-Schmidt procedure have been devised and are known in the art. The classical Gram-Schmidt procedure for orthogonalizing the matrix {circumflex over (M)}^T is described below.

[0074] For QR factorization, the matrix {circumflex over (M)}^T may be expressed as:

[0075] where Q_F is an (N_RN_R) orthogonal matrix; and

[0076] R_F is an (N_RN_T) upper triangle matrix with zeros below the diagonal and possible non-zero values along and above the diagonal.

[0077] The Gram-Schmidt procedure generates the matrices Q_F and R_F column-by-column. The following notations are used for the description below:

[0078] Q_Fq₁ q₂ . . . q_NR, where q_j is thej-th column of Q_F;

[0079] q_i,j is the entry in the i-th row and j-th column of Q_F;

[0080] {tilde over (Q)}_F{tilde over (q)}₁ {tilde over (q)}₂ . . . {tilde over (q)}_NR, where {tilde over (q)}_j is the j-th column of {tilde over (Q)}_F;

[0081] r_i,j is the entry in the i-th row and j-th column of R_F;

[0082] {circumflex over (M)}^T{circumflex over (m)}₁ {circumflex over (m)}₂ . . . {circumflex over (m)}_NT, where {circumflex over (m)}_j is the j-th column of {circumflex over (M)}^T; and

[0083] {circumflex over (m)}_i,j is the entry in the i-th row and j-th column of {circumflex over (M)}^T.

[0084] The first column of Q_F and R_F may be obtained as: $\begin{array}{cc}{r}_{1,1}={\underset{\_}{\hat{m}}}_1={\left[\underset{i=1}{\overset{N_R}{}}{{\hat{m}}_{i,1}}^2\right]}^{1/2},\mathrm{and}& \mathrm{Eq}\left(11\right)\\ {\underset{\_}{q}}_1=\frac{1}{r_{1,1}}{\underset{\_}{\hat{m}}}_1.& \mathrm{Eq}\left(12\right)\end{array}$

[0085] The first column of R_F includes one non-zero value r_1,1 for the first row and zeros elsewhere, where r_1,1 is the 2-norm of {circumflex over (m)}₁. The first column of Q_F is a normalized version of the first column of {circumflex over (M)}^T, where the normalization is achieved by scaling each entry of {circumflex over (m)}₁ with the inverse of r_1,1.

[0086] Each of the remaining columns of Q_F and R_F may be obtained as follows: $\begin{array}{cc}\mathrm{FOR}j=2,3N_T\mathrm{FOR}i=1,2j-1r_{i,j}={\underset{\_}{q}}_i^H{\underset{\_}{\hat{m}}}_j& \mathrm{Eq}\left(13\right)\\ {\underset{\_}{\tilde{q}}}_j={\underset{\_}{\hat{m}}}_j-\underset{i=1}{\overset{j-1}{}}{r}_{i,j}{\underset{\_}{q}}_i& \mathrm{Eq}\left(14\right)\\ r_{j,j}={\underset{\_}{\tilde{q}}}_j& \mathrm{Eq}\left(15\right)\\ {\underset{\_}{q}}_j=\frac{1}{r_{j,j}}{\underset{\_}{\tilde{q}}}_j.& \mathrm{Eq}\left(16\right)\end{array}$

[0087] The Gram-Schmidt procedure generates one column at a time for the matrix Q_F. Each new column of Q_F is forced to be orthogonal to all prior-generated columns to the left of the new column. This is achieved by equations (14) and (16), where the j-th column of Q_F (or q_j) is generated based on {tilde over (q)}_j, which in turn is generated based on the j-th column of {circumflex over (M)}^T (or {circumflex over (m)}_j) and subtracting out any components in {circumflex over (m)}_j pointing in the direction of the other (j-1) columns to the left of {circumflex over (m)}_j. The diagonal elements of R_F are computed as the 2-norm of the columns of {tilde over (Q)}_F (where {tilde over (q)}₁{circumflex over (m)}₁), as shown in equation (15).

[0088] Improved performance may be attained by ordering the matrix {circumflex over (M)}^T based on the singular value estimates before performing the QR factorization. The initial singular value estimates {tilde over ()}_m, for m1 . . . N_T, for diagonal matrix {tilde over ()} may be computed as the 2-norm of the columns of {circumflex over (M)}^T, as described below. The initial singular value estimates may then be ordered such that {tilde over ()}₁{tilde over ()}₂ . . . {tilde over (_NT, where {tilde over ()})}₁ is the largest singular value estimate and _NT is the smallest singular value estimate. When the initial singular value estimates for the diagonal matrix {tilde over ()} are ordered, the columns of the matrix {circumflex over (M)}^T are also ordered correspondingly. The first or left-most column of {circumflex over (M)}^T would then be associated with the largest singular value estimate and the highest received SNR, and the last or right-most column of {circumflex over (M)}^T would be associated with the smallest singular value estimate and the lowest received SNR . For the QR factorization, the initial singular value estimates may be obtained as the 2-norm of the columns of {circumflex over (M)}^T and used for ordering the columns of {circumflex over (M)}^T. The final singular value estimates are obtained as the 2-norm of the columns of {tilde over (Q)}_F, as described above. The steered reference may also be transmitted in order (e.g., from the largest eigenmode to the smallest eigenmode), so that the singular value estimates are effectively ordered by the transmitter.

[0089] If the columns of {circumflex over (M)}^T are ordered based on decreasing values of their associated singular value estimates, then the columns/eigenvectors of Q_F are forced to be orthogonal to the first column/eigenvector with the best received SNR. This ordering thus has the beneficial effect of rejecting certain noise components of each of the remaining eigenvectors of Q_F. In particular, the j-th column of Q_F (or q_j) is generated based on the j-th column of {circumflex over (M)}^T (or {circumflex over (m)}_j), and noise components in {circumflex over (m)}_j that point in the direction of the j-1 eigenvectors to the left of q_j (which are associated with higher received SNRs) are subtracted from {circumflex over (m)}_j to obtain q_j. The ordering also has the beneficial effect of improving the estimates of eigenvectors associated with smaller singular values. The overall result is improved performance, especially if the orthogonalized eigenvectors of Q_F are used for spatial processing for data transmission on the other link, as described below.

[0090] The enhanced orthogonal matched filter {tilde over (M)}_F obtained based on QR factorization may then be expressed as:

[0091] where {tilde over (R)}_F includes only the diagonal elements of R_F (i.e., the elements above the diagonal are set to zeros). The diagonal elements of {tilde over (R)}_F and R_F are estimates of the singular values of H. Since M^TU^H and {tilde over (M)}_F{tilde over (R)}_F^TQ_F^T, the following substitutions may be made: {tilde over (R)}_F and Q_FU*, where * denotes the complex conjugate.

[0092] 2. Mean Square Error Computation and Polar Decomposition

[0093] The initial matched filter matrix {circumflex over (M)} may also be orthogonalized based on a particular optimality criterion. One possible criterion is to minimize a measure of the squared error between the matrix {circumflex over (M)} and an optimum matched filter with the desired orthogonality properties. This may be expressed as:

[0094] where X_F is the Frobenius norm of X, and is given as: $\begin{array}{cc}{\underset{\_}{X}}_F={\left[\underset{i,j}{}{x_{\mathrm{ij}}}^2\right]}^{1/2}& \mathrm{Eq}\left(19\right)\end{array}$

[0095] The condition Q_P^HQ_PI ensures that Q_P is a unitary matrix, which would mean that the columns of Q_P are orthogonal to one another and the rows of Q_P are also orthogonal to one another. Equation (18) results in an optimum matched filter ^TQ_P that is the best fit to the measured data given by the matrix {circumflex over (M)}.

[0096] The solution to equation (18) can be obtained from the known solution to the orthogonal Procrustes problem. This problem asks the questiongiven two known matrices A and B, can a unitary matrix Q_P be found that rotates B into A. The problem may be expressed as:

[0097] The solution to the Procrustes problem can be obtained as follows. First, a matrix C_P is defined as C_PB^HA. The singular value decomposition of C_P is then given as C_PU_PPV_P^H or U_P^HC_PV_PP. The unitary matrix Q_P that solves the minimization problem shown in equation (20) is then given as:

[0098] The derivation and proof for equation (21) is described by G. H. Golub and C. F. Van Loan in Matrix Computation, Third Edition, Johns Hopkins University Press, 1996.

[0099] The solution for equation (20), which is shown in equation (21), is related to the polar decomposition of the matrix C. This polar decomposition is given as:

[0100] where Z_P is a unitary matrix, which is given as Z_PPV_P^H;

[0101] _P is a matrix of left eigenvectors of C_P that spans the column space of C_P (i.e., _P is equal to U_P or a sub-matrix of U_P depending on the dimension of C_P);

[0102] P_P is a Hermitian symmetric positive semi-definite matrix, which is given as

[0103] {tilde over ()}_P is a square matrix of singular values of C_P with dimension equal to the number of columns of C_P.

[0104] Polar decomposition can thus be performed on the matrix C_P to obtain the unitary matrix Z_P, which may be equal to either Q_P or a sub-matrix of Q_P depending on the dimension of C_P. It can be shown that the matrix Z_P is the optimal result to the minimization problem shown in equation (20).

[0105] Algorithms for direct computation of polar decomposition are described by P. Zielinski and K. Zietak in The Polar DecompositionProperties, Applications and Algorithms, Annals of the Polish Mathematical Society, 38 (1995), and by A. A. Dubrulle in An Optimum Iteration for the Matrix Polar Decomposition, Electronic Transactions on Numerical Analysis, Vol. 8, 1999, pp. 21-25.

[0106] The solution for the optimum matched filter expressed in equation (18) may be obtained based on the solution to the orthogonal Procrustes problem described above. This may be achieved by equating {circumflex over (M)} to A and ^T to B. For the computation, an estimate of the singular values, {tilde over ()}, may be obtained as the 2-norm of the columns of {circumflex over (M)}^T and used in place of . The diagonal elements of {tilde over ()} may be expressed as: $\begin{array}{cc}{\stackrel{~}{}}_{i,i}={\widehat{\underset{\_}{m}}}_i={\left[\underset{j=1}{\overset{N_R}{}}{{\hat{m}}_{i,j}}^2\right]}^{1/2},\mathrm{for}i\left\{1N_T\right\}.& \mathrm{Eq}\left(23\right)\end{array}$

[0107] It can be shown that the use of {tilde over ()} in the computation for Q_P results in nearly un-measurable degradation in performance relative to the use of the exact singular values in .

[0108] A matrix C_M may then be defined as:

[0109] The singular value decomposition of the matrix C_M is then given as:

[0110] The unitary matrix Q_M that solves the minimization problem shown in equation (18) is then given as:

[0111] An enhanced orthogonal matched filter e,otl M_M, which is the solution to the minimization problem in equation (18), may then be expressed as:

[0112] Alternatively, the polar decomposition of C_M may be performed as described above, which may be expressed as:

[0113] The unitary matrix Q_M that solves the minimization problem shown in equation (18) may then be given as:

[0114] The enhanced orthogonal matched filter {tilde over (M)}_M may then be expressed as:

[0115] It can be shown that the matrix Z_M from the polar decomposition is the optimal result for the matrix Q_M for the minimum square error computation (i.e., Q_MZ_M). Thus, the polar decomposition and minimum square error computation both yield the same orthogonal matched filter {tilde over (M)}_M.

[0116] FIG. 1 shows a flow diagram of an embodiment of a process 100 for deriving an orthogonal matched filter matrix {tilde over (M)} based on a steered reference. Initially, the receiver receives and processes the steered reference to obtain an estimate of u_mm or for each of multiple eigenmodes of H (step 112). This processing may be performed as described above. An initial matched filter matrix {circumflex over (M)} is then formed whose rows {circumflex over (m)}_m, for m1 . . . N_T, are derived based on the estimates of u_mm. The orthogonal matched filter matrix {tilde over (M)} may then be obtained from the initial matched filter matrix {circumflex over (M)} using any one of the orthogonalization techniques described above.

[0117] For the QR factorization technique, the matrix {circumflex over (M)} is factorized to obtain the matrices Q_F and R_F (step 122). The orthogonal matched filter matrix {tilde over (M)} is then obtained as shown in equation (17) (step 124) and the singular value estimates {tilde over ()} are obtained as the diagonal elements of R_F (step 126).

[0118] For the minimum square error technique, estimates of the singular values, {tilde over ()}, are obtained as the 2-norm of the columns of {circumflex over (M)}^T (step 132). The matrix C_M is then computed as shown in equation (24) (step 134). Singular value decomposition of C_M is next computed as shown in equation (25) (step 136). The orthogonal matched filter matrix {tilde over (M)} is then obtained as shown in equation (27) (step 138).

[0119] For the polar decomposition technique, estimates of the singular values, {tilde over ()}, are obtained as the 2-norm of the columns of {circumflex over (M)}^T (step 142). The matrix C_M is then computed as shown in equation (24) (step 144). Polar decomposition of C_M is next computed as shown in equation (28) (step 146). The orthogonal matched filter matrix {tilde over (M)} is then obtained as shown in equation (30) (step 148).

[0120] The orthogonal matched filter matrix {tilde over (M)} may thereafter be used to perform matched filtering of a received data transmission (step 150).

[0121] The orthogonalization of the matched filter matrix provides several benefits. First, the use of an orthogonal matched filter matrix {tilde over (M)} avoids cross-talk between the eigenmodes of H. The derivation of the initial matched filter matrix {circumflex over (M)} piecemeal based on the steered reference does not guarantee that the eigenvectors of {circumflex over (M)}^T are orthogonal. The lack of orthogonality results in performance degradation. The orthogonalization of the matched filter matrix avoids this performance degradation.

[0122] Second, QR factorization can improve the quality of the eigenvectors associated with smaller singular values. Without QR factorization, the quality of the estimates of the eigenvectors is not constant, and the estimates of the eigenvectors associated with smaller singular values are likely to be lower in quality. QR factorization can improve the quality of the eigenvectors associated with smaller singular values by rejecting certain noise components, as described above. Polar decomposition may have similar effect, but not in the direct way as QR factorization.

[0123] Third, orthogonalization may reduce the amount of resources needed to transmit the steered reference. If orthogonalization is not performed, then high quality estimates of and U would be needed to ensure low cross-talk among the eigenmodes. A longer transmission period would then be needed for the steered reference for the eigenvectors associated with smaller singular values to ensure that the desired quality is obtained. High quality estimates of and U would thus require a longer period of transmission for the steered reference (which would consume more valuable system resources) and a longer integration period for the steered reference at the receiver (which may result in longer delay for data transmission). Orthogonalization can provide the desired performance without the need for high quality estimates of and U.

[0124] III. MIMO-OFDM System

[0125] The techniques for deriving eigenvectors used for spatial processing are now described for an exemplary wideband MIMO communication system that employs orthogonal frequency division multiplexing (OFDM). OFDM effectively partitions the overall system bandwidth into a number of (N_F) orthogonal subbands, which are also referred to as tones, frequency bins, or frequency subchannels. With OFDM, each subband is associated with a respective subcarrier upon which data may be modulated. For a MIMO-OFDM system, each subband may be associated with multiple eigenmodes, and each eigenmode of each subband may be viewed as an independent transmission channel.

[0126] For OFDM, the data or pilot to be transmitted on each usable subband is first modulated (i.e., mapped to modulation symbols) using a particular modulation scheme. One modulation symbol may be transmitted on each usable subband in each symbol period. A signal value of zero may be sent for each unused subband. For each OFDM symbol period, the modulation symbols for the usable subbands and zero signal values for the unused subbands (i.e., the modulation symbols and zeros for all Nsubbands) are transformed to the time domain using an inverse fast Fourier transform (IFFT) to obtain a transformed symbol that comprises N_F time-domain samples. To combat inter-symbol interference (ISI) caused by frequency selective fading, a portion of each transformed symbol is often repeated (which is often referred to as adding a cyclic prefix) to form a corresponding OFDM symbol. The OFDM symbol is then processed and transmitted over the wireless channel. An OFDM symbol period, which is also referred to as a symbol period, corresponds to the duration of one OFDM symbol.

[0127] For this exemplary system, the downlink and uplink share a single frequency band using time-division duplex (TDD). For a TDD MIMO-OFDM system, the downlink and uplink channel responses may be assumed to be reciprocal of one another. That is, if H(k) represents a channel response matrix from antenna array A to antenna array B for subband k, then a reciprocal channel implies that the coupling from array B to array A is given by H^T(k).

[0128] FIG. 2 shows a wireless communication system 200 that includes a number of access points (APs) 210 that communicate with a number of user terminals (UTs) 220. (For simplicity, only one access point is shown in FIG. 2.) An access point may also be referred to as a base station or some other terminology. Each user terminal may be a fixed or mobile terminal and may also be referred to as an access terminal, a mobile station, a remote station, a user equipment (UE), a wireless device, or some other terminology. Each user terminal may communicate with one or possibly multiple access points on the downlink and/or the uplink at any given moment. The downlink (i.e., forward link) refers to transmission from the access point to the user terminal, and the uplink (i.e., reverse link) refers to transmission from the user terminal to the access point. The channel response between each access point and each user terminal may be characterized by a set of channel response matrices H(k), for kK where K represents the set of all subbands of interest (e.g., the usable subbands).

[0129] In the following description for a pair of communicating access point and user terminal, it is assumed that calibration has been performed to account for differences between the transmit and receive chains of the access point and the user terminal. The results of the calibration are diagonal matrices {circumflex over (K)}_ap(k) and {circumflex over (K)}_ut(k), for kK, to be used at the access point and the user terminal, respectively, on the transmit path. A calibrated downlink channel response, H_cdn(k), observed by the user terminal and a calibrated uplink channel response, H_cup(k), observed by the access point may then be expressed as:

[0130] where H_dn(k)R_ut(k)H(k)T_ap(k) is the effective downlink channel response, which includes the responses of the transmit chain T_ap(k) at the access point and the receive chain R_ut(k) at the user terminal;

[0131] H_up(k)R_ap(k)H^T(k)T_ut(k) is the effective uplink channel response, which includes the responses of the transmit chain T_ut(k) at the user terminal and the receive chain R_ap(k) at the access point; and

[0132] H(k) is an (N_utN_ap) channel response matrix between the N_ap antennas at the access point and the N_ut antennas at the user terminal.

[0133] If calibration is not performed, then the matrices {circumflex over (K)}_ap(k) and {circumflex over (K)}_ut(k), for kK, are each set to the identity matrix I.

[0134] FIG. 3 shows an embodiment of a frame structure 300 that may be used for a TDD MIMO-OFDM system. Data transmission occurs in units of TDD frames, with each TDD frame covering a particular time duration (e.g., 2 msec). Each TDD frame is partitioned into a downlink phase and an uplink phase. The downlink phase is further partitioned into multiple segments for multiple downlink transport channels. In the embodiment shown in FIG. 3, the downlink transport channels include a broadcast channel (BCH), a forward control channel (FCCH), and a forward channel (FCH). Similarly, the uplink phase is partitioned into multiple segments for multiple uplink transport channels. In the embodiment shown in FIG. 3, the uplink transport channels include a reverse channel (RCH) and a random access channel (RACH).

[0135] In the downlink phase, a BCH segment 310 is used to transmit one BCH protocol data unit (PDU) 312, which includes a beacon pilot 314, a MIMO pilot 316, and a BCH message 318. The beacon pilot is transmitted from all access point antennas and is used by the user terminals for timing and frequency acquisition. The MIMO pilot is transmitted from all access point antennas with different orthogonal codes and is used by the user terminals for channel estimation. The BCH message carries system parameters for the user terminals in the system. An FCCH segment 320 is used to transmit one FCCH PDU, which carries assignments for downlink and uplink resources and other signaling for the user terminals. An FCH segment 330 is used to transmit one or more FCH PDUs 332. Different types of FCH PDU may be defined. For example, an FCH PDU 332a includes only a data packet 336a, and an FCH PDU 332b includes a downlink steered reference 334b and a data packet 336b.

[0136] In the uplink phase, an RCH segment 340 is used to transmit one or more RCH PDUs 342 on the uplink. Different types of RCH PDU may also be defined. For example, an RCH PDU 342a includes an uplink steered reference 344a and a data packet 346a. An RACH segment 350 is used by the user terminals to gain access to the system and to send short messages on the uplink. An RACH PDU 352 may be sent within RACH segment 350 and includes an uplink steered reference 354 and a message 356.

[0137] For the embodiment shown in FIG. 3, the beacon and MIMO pilots are sent on the downlink in the BCH segment in each TDD frame. A steered reference may or may not be sent in any given FCH/RCH PDU. A steered reference may also be sent in an RACH PDU to allow the access point to estimate pertinent vectors during system access.

[0138] For simplicity, the following description is for a communication between one access point and one user terminal. The MIMO pilot is transmitted by the access point and used by the user terminal to obtain an estimate of the calibrated downlink channel response, _cdn(k), for kK. The calibrated uplink channel response may then be estimated as _cup(k)_cdn^T(k). Singular value decomposition may be performed to diagonalized the matrix _cup(k) for each subband, which may be expressed as:

[0139] where _ap(k) is an (N_apN_ap) unitary matrix of left eigenvectors of _cup(k);

[0140] {circumflex over ()}(k) is an (N_apN_ut) diagonal matrix of singular values of _cup(k); and

[0141] {circumflex over (V)}_ut(k) is an (N_utN_ut) unitary matrix of right eigenvectors of _cup(k).

[0142] Similarly, the singular value decomposition of the estimated calibrated downlink channel response matrix, _cdn(k), may be expressed as:

[0143] where the matrices {circumflex over (V)}_ut*(k) and _ap*(k) are unitary matrices of left and right eigenvectors, respectively, of _cdn(k).

[0144] As shown in equations (32) and (33), the matrices of left and right eigenvectors for one link are the complex conjugate of the matrices of right and left eigenvectors, respectively, for the other link. For simplicity, reference to the matrices _ap(k) and {circumflex over (V)}_ut(k) in the following description may also refer to their various other forms (e.g., {circumflex over (V)}_ut(k) may refer to {circumflex over (V)}_ut(k), {circumflex over (V)}_ut*(k), {circumflex over (V)}_ut^T(k), and {circumflex over (V)}_ut^H(k)). The matrices _ap(k) and {circumflex over (V)}_ut(k) may be used by the access point and user terminal, respectively, for spatial processing and are denoted as such by their subscripts. The matrix {circumflex over ()}(k) includes singular value estimates that represent the gains for the independent channels (or eigenmodes) of the channel response matrix H(k) for the k-th subband.

[0145] The singular value decomposition may be performed independently for the channel response matrix _cup(k) for each of the usable subbands to determine the N_S eigenmodes for the subband. The singular value estimates for each diagonal matrix {circumflex over ()}(k) may be ordered such that {circumflex over ()}₁(k){circumflex over ()}₂(k) . . . {circumflex over ()}_NS(k), where {circumflex over ()}₁(k) is the largest singular value estimate and {circumflex over ()}_NS(k) is the smallest singular value estimate for subband k. When the singular value estimates for each diagonal matrix {circumflex over ()}(k) are ordered, the eigenvectors (or columns) of the associated matrices (k) and {circumflex over (V)}(k) are also ordered correspondingly. After the ordering, {circumflex over ()}₁(k) represents the singular value estimate for the best eigenmode for subband k, which is also often referred to as the principal eigenmode.

[0146] A wideband eigenmode may be defined as the set of same-order eigenmodes of all subbands after the ordering. Thus, the m-th wideband eigenmode includes the m-th eigenmodes of all subbands. Each wideband eigenmode is associated with a respective set of eigenvectors for all of the subbands. The principal wideband eigenmode is the one associated with the largest singular value estimate in the matrix {circumflex over ()}(k) for each of the subbands.

[0147] The user terminal can transmit a steered reference on the uplink. The uplink steered reference for the m-th wideband eigenmode may be expressed as:

[0148] where {circumflex over (v)}_ut,m(k) is the m-th column of the matrix {circumflex over (V)}_ut(k) for the k-th subband, with

[0149] p(k) is the pilot symbol for the k-th subband.

[0150] The received uplink steered reference at the access point may be expressed as: $\begin{array}{cc}\begin{array}{ccc}{\underset{\_}{r}}_{\mathrm{up},\mathrm{sr},m}& ={\underset{\_}{H}}_{\mathrm{up}}\left(k\right){\underset{\_}{x}}_{\mathrm{up},\mathrm{sr},m}\left(k\right)+{\underset{\_}{n}}_{\mathrm{up}}\left(k\right)& ,\mathrm{for}kK,\\ & ={\underset{\_}{H}}_{\mathrm{up}}\left(k\right){\underset{\_}{\hat{K}}}_{\mathrm{ut}}\left(k\right){\underset{\_}{\hat{v}}}_{\mathrm{ut},m}\left(k\right)p\left(k\right)+{\underset{\_}{n}}_{\mathrm{up}}\left(k\right)& \\ & {\underset{\_}{\hat{u}}}_{\mathrm{ap},m}\left(k\right){\hat{}}_m\left(k\right)p\left(k\right)+{\underset{\_}{n}}_{\mathrm{up}}\left(k\right)& \end{array}& \mathrm{Eq}\left(35\right)\end{array}$

[0151] where _ap,m(k) is the m-th column of the matrix _ap(k) for the k-th subband, with

[0152] {circumflex over ()}_m(k) is the singular value estimate for the k-th subband of the m-th wideband eigenmode.

[0153] The access point can obtain an initial matched filter matrix {circumflex over (M)}_ap(k), for kK, based on the uplink steered reference, as described above. The access point may thereafter obtain an enhanced orthogonal matched filter matrix {tilde over (M)}_ap(k), for kK, based on {circumflex over (M)}_ap(k) and using any one of the orthogonalization techniques described above.

[0154] Using QR factorization, the matrix {tilde over (M)}_ap(k) may be obtained as:

[0155] where Q_ap(k) is a unitary matrix that is the ortho-normal basis for {tilde over (M)}_ap(k);

[0156] {tilde over (R)}_ap(k) is a diagonal matrix derived based on {circumflex over (M)}_ap(k); and

[0157] {tilde over ()}_ap(k){tilde over (R)}_ap(k) and _ap*(k)Q_ap(k).

[0158] Using mean square error computation, the matrix {tilde over (M)}_ap(k) may be obtained as:

[0159] {tilde over ()}_ap(k) is the diagonal matrix whose elements are the 2-norm of the columns of {circumflex over (M)}_ap^T(k); and

[0160] Using polar decomposition, the matrix {tilde over (M)}_ap(k) may be obtained as:

[0161] The matrix {tilde over (M)}_ap(k) may be used by the access point for matched filtering of uplink data transmission from the user terminal, as described below.

[0162] The spatial processing performed by the user terminal to transmit data on multiple eigenmodes on the uplink may be expressed as:

[0163] where s_up(k) is the data vector and x_up(k) is the transmit vector for the k-th subband for the uplink. Uplink data transmission can occur on any number of wideband eigenmodes from 1 to N_S.

[0164] The received uplink data transmission at the access point may be expressed as: $\begin{array}{cc}\begin{array}{c}{\underset{\_}{r}}_{\mathrm{up}}\left(k\right)={\underset{\_}{H}}_{\mathrm{up}}\left(k\right){\underset{\_}{x}}_{\mathrm{up}}\left(k\right)+{\underset{\_}{n}}_{\mathrm{up}}\left(k\right)\\ ={\underset{\_}{H}}_{\mathrm{up}}\left(k\right){\underset{\_}{\hat{K}}}_{\mathrm{ut}}\left(k\right){\underset{\_}{\hat{V}}}_{\mathrm{ut}}\left(k\right){\underset{\_}{s}}_{\mathrm{up}}\left(k\right)+{\underset{\_}{n}}_{\mathrm{up}}\left(k\right),\\ \mathrm{for}kK,\\ ={\underset{\_}{\hat{U}}}_{\mathrm{ap}}\left(k\right)\underset{\_}{\hat{}}\left(k\right){\underset{\_}{s}}_{\mathrm{up}}\left(k\right)+{\underset{\_}{n}}_{\mathrm{up}}\left(k\right)\end{array}& \mathrm{Eq}\left(42\right)\end{array}$

[0165] where r_up(k) is the receive vector for the uplink data transmission for the k-th subband.

[0166] The matched filtering by the access point may be expressed as: $\begin{array}{cc}\begin{array}{c}{\underset{\_}{\hat{s}}}_{\mathrm{up}}\left(k\right)={\underset{\_}{G}}_{\mathrm{ap}}\left(k\right){\underset{\_}{\tilde{M}}}_{\mathrm{ap}}\left(k\right){\underset{\_}{r}}_{\mathrm{up}}\left(k\right)\\ ={\underset{\_}{G}}_{\mathrm{ap}}\left(k\right){\widetilde{\underset{\_}{}}}_{\mathrm{ap}}^T\left(k\right){\underset{\_}{\tilde{U}}}_{\mathrm{ap}}^H\left({\widehat{\underset{\_}{U}}}_{\mathrm{ap}}\left(k\right)\underset{\_}{\hat{}}\left(k\right){\underset{\_}{s}}_{\mathrm{up}}\left(k\right)+{\underset{\_}{n}}_{\mathrm{up}}\left(k\right)\right),\\ \mathrm{for}kK,\\ {\underset{\_}{s}}_{\mathrm{up}}\left(k\right)+{\underset{\_}{\tilde{n}}}_{\mathrm{up}}\left(k\right)\end{array}& \mathrm{Eq}\left(43\right)\end{array}$

[0167] For the TDD MIMO system, the access point may also use the matrices _ap(k), for kK, for spatial processing for data transmission on the downlink to the user terminal. The spatial processing performed by the access point to transmit data on multiple eigenmodes on the downlink may be expressed as:

[0168] where s_dn(k) is the data vector and x_dn(k) is the transmit vector for the k-th subband for the downlink. Downlink data transmission can similarly occur on any number of wideband eigenmodes from 1 to N_S.

[0169] The received downlink data transmission at the user terminal may be expressed as: $\begin{array}{cc}\begin{array}{c}{\underset{\_}{r}}_{\mathrm{dn}}\left(k\right)={\underset{\_}{H}}_{\mathrm{dn}}\left(k\right){\underset{\_}{x}}_{\mathrm{dn}}\left(k\right)+{\underset{\_}{n}}_{\mathrm{dn}}\left(k\right)\\ ={\underset{\_}{H}}_{\mathrm{dn}}\left(k\right){\underset{\_}{\hat{K}}}_{\mathrm{ap}}\left(k\right){\underset{\_}{\tilde{U}}}_{\mathrm{ap}}^{*}\left(k\right){\underset{\_}{s}}_{\mathrm{dn}}\left(k\right)+{\underset{\_}{n}}_{\mathrm{dn}}\left(k\right),\\ \mathrm{for}kK,\\ ={\underset{\_}{\hat{V}}}_{\mathrm{ut}}^{*}\left(k\right)\underset{\_}{\hat{}}\left(k\right){\underset{\_}{s}}_{\mathrm{dn}}\left(k\right)+{\underset{\_}{n}}_{\mathrm{dn}}\left(k\right)\end{array}& \mathrm{Eq}\left(45\right)\end{array}$

[0170] where r_dn(k) is the receive vector for the downlink data transmission for the k-th subband.

[0171] The matched filtering by the user terminal may be expressed as: $\begin{array}{cc}\begin{array}{c}{\underset{\_}{\hat{s}}}_{\mathrm{dn}}\left(k\right)={\underset{\_}{G}}_{\mathrm{ut}}\left(k\right){\underset{\_}{\hat{M}}}_{\mathrm{ut}}\left(k\right){\underset{\_}{r}}_{\mathrm{dn}}\left(k\right)\\ ={\underset{\_}{G}}_{\mathrm{ut}}\left(k\right){\underset{\_}{\hat{}}}^T\left(k\right){\underset{\_}{\hat{V}}}_{\mathrm{ut}}^T\left(k\right)\left({\underset{\_}{\hat{V}}}_{\mathrm{ut}}^{*}\left(k\right)\underset{\_}{\hat{}}\left(k\right){\underset{\_}{s}}_{\mathrm{dn}}\left(k\right)+{\underset{\_}{n}}_{\mathrm{dn}}\left(k\right)\right),\\ \mathrm{for}kK,\\ {\underset{\_}{s}}_{\mathrm{dn}}\left(k\right)+{\underset{\_}{\tilde{n}}}_{\mathrm{dn}}\left(k\right)\end{array}& \mathrm{Eq}\left(46\right)\end{array}$

[0172] where {circumflex over (M)}_ut(k){circumflex over ()}^T(k){circumflex over (V)}_ut^T(k) is the matched filter for the user terminal;

[0173] The diagonal matrix {circumflex over ()}(k) is derived from the singular value decomposition shown in equation (32).

[0174] Table 1 summarizes the spatial processing at the access point and user terminal for both data transmission and reception on multiple wideband eigenmodes.

	TABLE 1
	Downlink	Uplink
Access	Transmit:	Receive:
Point	xdn (k)	up (k)
	{circumflex over (K)}ap (k)*ap (k)sdn (k)	Gap (k){tilde over ()}T (k)Hap (k)rup (k)
User	Receive:	Transmit:
Terminal	dn (k)	xup (k)
	Gut (k){circumflex over ()}T (k){circumflex over (V)}utT (k)rdn (k)	{circumflex over (K)}ut (k){circumflex over (V)}ut (k)sup (k)

[0175] In Table 1, s(k) is the data vector, x(k) is the transmit vector, r(k) is the receive vector, and (k) is an estimate of the data vector s(k), where all vectors are for subband k. The subscripts dn and up for these vectors denote downlink and uplink transmissions, respectively.

[0176] It can be shown that the use of the matrices _ap(k), for kK, (with orthogonalized columns) for spatial processing for downlink data transmission can provide substantial improvement over the use of matrices _ap(k), for kK, (with un-orthogonalized columns) obtained from the initial matched filter matrices {circumflex over (M)}_ap(k), for kK.

[0177] FIG. 4 shows transmission of steered reference and data on the downlink and uplink for an exemplary transmission scheme. The MIMO pilot is transmitted on the downlink by the access point in each TDD frame (block 412). The user terminal receives and processes the downlink MIMO pilot to obtain an estimate the downlink channel response _cdn(k), for kK. The user terminal then estimates the uplink channel response as _cup(k)_cdn^T(k) and performs singular value decomposition of _cup(k) to obtain the matrices {circumflex over ()}(k) and {circumflex over (V)}_ut(k), for kK, as shown in equation (32) (block 414).

[0178] The user terminal then transmits the uplink steered reference on the RACH or the RCH using the matrices {circumflex over (V)}_ut(k), for kK, as shown in equation (34) and FIG. 3, during system access (step 422). The columns of {circumflex over (V)}_ut(k) are also referred to as steering vectors when used for data transmission. The access point receives and processes the uplink steered reference on the RACH or the RCH to obtain the matrices {tilde over ()}(k) and _ap(k), for kK, as described above (step 424). The columns of _ap(k) are eigenvectors that may be used for both data reception as well as data transmission. The user terminal may thereafter transmit the uplink steered reference and data on the RCH using the matrices {circumflex over (V)}_ut(k), for kK, as shown in equation (41) and FIG. 3 (step 432). The access point receives and processes the uplink steered reference on the RCH to update the matrices {tilde over ()}(k) and _ap(k), for kK (step 434). The access point also performs matched filtering for the received uplink data transmission using the matrices {tilde over ()}(k) and _ap(k) (also step 434).

[0179] The access point may thereafter transmit an optional downlink steered reference and data on the FCH using the matrices _ap(k), for kK, as shown in equation (44) and FIG. 3 (step 442). If a downlink steered reference is transmitted, then the user terminal can process the downlink steered reference to update the matrices {circumflex over ()}(k) and {circumflex over (V)}_ut(k), for kK (step 444) and may also perform orthogonalization to ensure that the columns of {circumflex over (V)}_ut(k) are orthogonal. The user terminal also performs matched filtering for the received downlink data transmission using the matrices {circumflex over ()}(k) and {circumflex over (V)}_ut(k) (also step 444).

[0180] The pilot and data transmission scheme shown in FIG. 4 provides several advantages. First, the MIMO pilot transmitted by the access point may be used by multiple user terminals in the system to estimate the response of their respective MIMO channels. Second, the computation for the singular value decomposition of _cup(k), for kK, is distributed among the user terminals (i.e., each user terminal performs singular value decomposition of its own set of estimated channel response matrices for the usable subbands). Third, the access point can obtain the matrices {tilde over ()}(k) and _ap(k), for kK, which are used for uplink and downlink spatial processing, based on the steered reference without having to estimate the MIMO channel response.

[0181] Various other transmission schemes may also be implemented for MIMO and MIMO-OFDM systems, and this is within the scope of the invention. For example, the MIMO pilot may be transmitted by the user terminal and the steered reference may be transmitted by the access point.

[0182] FIG. 5 shows a block diagram of an embodiment of an access point 210x and a user terminal 220x in MIMO-OFDM system 200. For clarity, in this embodiment, access point 210x is equipped with four antennas that can be used for data transmission and reception, and user terminal 220x is also equipped with four antennas for data transmission/reception. In general, the access point and user terminal may each be equipped with any number of transmit antennas and any number of receive antennas.

[0183] On the downlink, at access point 210x, a transmit (TX) data processor 514 receives traffic data from a data source 512 and signaling and other data from a controller 530. TX data processor 514 formats, codes, interleaves, and modulates the data to provide modulation symbols, which are also referred to as data symbols. A TX spatial processor 520 then receives and multiplexes the data symbols with pilot symbols, performs the required spatial processing with the matrices _ap*(k), for kK, and provides four streams of transmit symbols for the four transmit antennas. Each modulator (MOD) 522 receives and processes a respective transmit symbol stream to provide a corresponding downlink modulated signal. The four downlink modulated signals from modulators 522a through 522d are then transmitted from antennas 524a through 524d, respectively.

[0184] At user terminal 220x, four antennas 552a through 552d receive the transmitted downlink modulated signals, and each antenna provides a received signal to a respective demodulator (DEMOD) 554. Each demodulator 554 performs processing complementary to that performed by modulator 522 and provides received symbols. A receive (RX) spatial processor 560 then performs matched filtering on the received symbols from all demodulators 554a through 554d to provide recovered data symbols, which are estimates of the data symbols transmitted by the access point. An RX data processor 570 further processes (e.g., symbol demaps, deinterleaves, and decodes) the recovered data symbols to provide decoded data, which may be provided to a data sink 572 for storage and/or a controller 580 for further processing.

[0185] RX spatial processor 560 also processes the received pilot symbols to obtain an estimate of the downlink channel response, _cdn(k), for kK. Controller 580 may then decompose each matrix _cdn(k) to obtain {circumflex over ()}(k) and {circumflex over (V)}_ut(k). Controller 580 may further derive (1) the downlink matched filter matrices {circumflex over (M)}_ut(k), for kK, based on {circumflex over ()}(k) and {circumflex over (V)}_ut(k), and (2) the scaling matrices G_ut(k), for kK, based on {circumflex over ()}(k). Controller 580 may then provide {circumflex over (M)}_ut(k) to RX data processor 560 for downlink matched filtering and {circumflex over (V)}_ut(k) to a TX spatial processor 590.

[0186] The processing for the uplink may be the same or different from the processing for the downlink. Data and signaling are processed (e.g., coded, interleaved, and modulated) by a TX data processor 588, multiplexed with pilot symbols, and further spatially processed by TX spatial processor 590 with the matrices {circumflex over (V)}_ut(k), for kK. The transmit symbols from TX spatial processor 590 are further processed by modulators 554a through 554d to generate four uplink modulated signals, which are then transmitted via antennas 552a through 552d.

[0187] At access point 510, the uplink modulated signals are received by antennas 524a through 524d and demodulated by demodulators 522a through 522d to provide received symbols for the uplink steered reference and data transmission. An RX spatial processor 540 then processes the received uplink steered reference to obtain estimates of u_mm, for kK and m1 . . . N_S, which are provided to controller 530. Controller then obtains {circumflex over (M)}_ap(k) and {tilde over ()}(k) based on the estimates of u_mm, performs orthogonalization of {circumflex over (M)}_ap(k) to obtain {tilde over (M)}_ap(k) and _ap(k), and derives G_ap(k) based on {tilde over ()}(k). Controller 580 then provides {tilde over (M)}_ap(k) and G_ap(k) to RX spatial processor 540 for uplink matched filtering and _ap*(k) to TX spatial processor 520 for downlink spatial processing.

[0188] RX spatial processor 540 performs matched filtering of the received uplink data transmission with {tilde over (M)}_ap(k) and G_ap(k) to provide recovered data symbols, which are further processed by an RX data processor 542 to provide decoded data. The decoded data may be provided to a data sink 544 for storage and/or controller 530 for further processing.

[0189] Controller 530 performs the processing to obtain the matched filter matrices {tilde over (M)}_ap(k) and the scaling matrices G_ap(k), for kK, for uplink data transmission and the matrices _ap*(k), for kK, for downlink data transmission. Controller 580 performs the processing to obtain the matched filter matrices {circumflex over (M)}_ut(k) and the scaling matrices G_ut(k), for kK, for downlink data transmission and the matrices {circumflex over (V)}_ut(k), for kK, for uplink data transmission. Controllers 530 and 580 further control the operation of various processing units at the access point and user terminal, respectively. Memory units 532 and 582 store data and program codes used by controllers 530 and 580, respectively.

[0190] FIG. 6 shows a block diagram of the spatial processing performed by access point 210x and user terminal 220x to transmit data on multiple eigenmodes on the downlink and uplink.

[0191] On the downlink, within TX spatial processor 520 at access point 210x, the data vector s_dn(k) for each subband k is first multiplied with the matrix _ap*(k) by a unit 610 and further multiplied with the correction matrix {circumflex over (K)}_ap(k) by a unit 612 to obtain the transmit vector x_dn(k) for subband k. The columns of the matrix _ap*(k) are orthogonalized as described above. The transmit vectors x_dn(k), for kK, are then processed by a transmit chain 614 within modulator 522 and transmitted over the MIMO channel to user terminal 220x. Unit 610 performs the spatial processing for downlink data transmission.

[0192] At user terminal 220x, the downlink modulated signals are processed by a receive chain 654 within demodulator 554 to obtain the receive vectors r_dn(k), for kK. Within RX spatial processor 560, the receive vector r_dn(k) for each subband k is first multiplied with the matched filter matrix {circumflex over (M)}_ut(k) by a unit 656 and further multiplied with the scaling matrix G_ut(k) by a unit 658 to obtain the vector _dn(k), which is an estimate of the data vector s_dn(k) transmitted for subband k. Units 656 and 658 perform the downlink matched filtering.

[0193] On the uplink, within TX spatial processor 590 at user terminal 220x, the data vector s_up(k) for each subband k is first multiplied with the matrix {circumflex over (V)}_ut(k) by a unit 660 and then further multiplied with the correction matrix {circumflex over (K)}_ut(k) by a unit 662 to obtain the transmit vector x_up(k) for subband k. The transmit vectors x_up(k), for kK, are then processed by a transmit chain 664 within modulator 554 and transmitted over the MIMO channel to access point 210x. Unit 660 performs the spatial processing for uplink data transmission.

[0194] At access point 210x, the uplink modulated signals are processed by a receive chain 624 within demodulator 522 to obtain the receive vectors r_up(k), for kK. Within RX spatial processor 540, the receive vector r_up(k) for each subband k is first multiplied with the matched filter matrix {tilde over (M)}_ap(k) by a unit 626 and further multiplied by the scaling matrix G_ap(k) by a unit 628 to obtain the vector _up(k), which is an estimate of the data vector s_up(k) transmitted for subband k. Units 626 and 628 perform the uplink matched filtering.

[0195] The techniques described herein to derive eigenvectors for spatial processing may be implemented by various means. For example, these techniques may be implemented in hardware, software, or a combination thereof. For a hardware implementation, the elements used for these techniques may be implemented within one or more application specific integrated circuits (ASICs), digital signal processors (DSPs), digital signal processing devices (DSPDs), programmable logic devices (PLDs), field programmable gate arrays (FPGAs), processors, controllers, micro-controllers, microprocessors, other electronic units designed to perform the functions described herein, or a combination thereof.

[0196] For a software implementation, the techniques may be implemented with modules (e.g., procedures, functions, and so on) that perform the functions described herein. The software codes may be stored in a memory unit (e.g., memory units 532 and 582 in FIG. 5) and executed by a processor (e.g., controllers 530 and 580). The memory unit may be implemented within the processor or external to the processor, in which case it can be communicatively coupled to the processor via various means as is known in the art.

[0197] Headings are included herein for reference and to aid in locating certain sections. These headings are not intended to limit the scope of the concepts described therein under, and these concepts may have applicability in other sections throughout the entire specification.

[0198] The previous description of the disclosed embodiments is provided to enable any person skilled in the art to make or use the present invention. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the generic principles defined herein may be applied to other embodiments without departing from the spirit or scope of the invention. Thus, the present invention is not intended to be limited to the embodiments shown herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.