US 6,549,585 B2Grant

Combined interference cancellation and maximum likelihood decoding of space-time block codes

Issue Date:2003-04-15

•24 Claims

•1 Drawing Sheets

Abstract

Block-encoded transmissions of a multi-antenna terminal unit are effectively detected in the presence of co-channel interfering transmissions from other multi-antenna terminal units, when the base station has a plurality of antennas, and interference cancellation is combined with maximum likelihood decoding. The signals received in one base station antenna are employed in processing the signals received in a second base station antenna so as to cancel the signals of one terminal unit, while decoding the signals transmitted by the other terminal unit. Zero-forcing and MMSE approaches are presented. In another embodiment of this invention the basic decoding approach is used to obtain an initial estimate for the symbols from each terminal. Assuming that the signals of the first terminal unit has been decoded correctly, the receiver employs this initial decoded signal of the first terminal unit to cancel their contribution to the signals received at the base station antennas while decoding the signals of the second terminal unit. This process is then repeated assuming that the signals of the second terminal unit has been decoded correctly, the receiver employs this initial decoded signal of the second terminal unit to cancel their contribution to the signals received at the base station antennas while decoding the signals of the first terminal unit. The above disclosed techniques are viable for any number K of terminal units concurrently transmitting over a given channel, where each terminal unit is using a space-time block code with N transmit antennas, and a base station has at least K receive antennas.

Metadata

Assignee

AT&T Corp

Inventors

Ayman F. Naguib
Nambirajan Seshadri

Application Information

Application Number:US 09/767,019

Filing Date:2001-01-22

Priority Date:1997-10-06

Art Unit:7

Classifications

IPC:

H04B 702H04B 169

Field of Search:

375267375347375148375262375341375346375144375348375349

Patent Drawings (1 sheets)

Description

[0002] This application is a continuation of application Ser. No. 09/149,163, filed Sep. 4, 1998, now U.S. Pat. No. 6,178,196, which is claims the benefit of U.S. Provisional Application No. 60/061,145, filed Oct. 6, 1997. This application is also related to U.S. application Ser. No. 09/074224, filed May 7, 1998, now U.S. Pat. No. 6,185,258, titled “Transmitter Diversity Technique for Wireless Communications”.

Background of the Invention

[0003] This invention relates to wireless communication and, more particularly, to techniques for effective wireless communication in the presence of fading, co-channel interference, and other degradations.

[0004] Rapid growth in mobile computing and other wireless data services is inspiring many proposals for high speed data services in the range of 64-144 kbps for micro cellular wide area and high mobility applications and up to 2 Mbps for indoor applications. Research challenges include the development of efficient coding and modulation, and signal processing techniques to improve the quality and spectral efficiency of wireless communications and better techniques for sharing the limited spectrum among different high capacity users.

[0005] The physical limitation of the wireless channel presents a fundamental technical challenge for reliable communications. The channel is susceptible to time-varying noise, interference, and multipaths. Power and size limitations of the communications and computing device in a mobile handset constitute another major design consideration. Most personal communications and wireless services portables are meant to be carried in a briefcase and/or pocket and must, therefore, be small and lightweight. This translates to a low power requirement since small batteries must be used. However, many of the signal processing techniques which may be used for reliable communications and efficient spectral utilization demand significant processing power, precluding the use of low power devices. Continuing advances in VLSI and integrated circuit technology for low power applications will provide a partial solution to this problem. Still, placing most of the signal processing burden on fixed locations (base stations) with relatively larger power resources than the mobile units will, likely, continue to be the trend in wireless systems design.

[0006] Perhaps the single most important parameter in providing reliable communications over wireless channels is diversity. Diversity techniques which may be used include time, frequency, and space diversity

[0007] Time diversity: channel coding in combination with limited interleaving is used to provide time diversity. However, while channel coding is extremely effective in fast fading environments (high mobility), it offers very little protection under slow fading (low mobility) unless significant interleaving delays can be tolerated.

[0008] Frequency diversity: the fact that signals transmitted over different frequencies induce different multipath structures and independent fading is. However, when the multipath delay spread is small compared to the symbol period, frequency diversity is not helpful.

[0009] Space diversity: the receiver/transmitter uses multiple antennas that are separated for reception/transmission and/or differently polarized antennas to create independent fading channels. Currently, multiple antennas at base-stations are used for receive diversity at the base. However, it is difficult to have more than two antennas at the mobile unit due to the size and cost of multiple chains of RF conversion.

[0010] Previous work on transmit diversity can be classified into three broad categories: schemes using feedback, schemes with feedforward or training information but no feedback, and blind schemes. The third category (blind schemes) relies primarily on multiple transmit antennas combined with channel coding to provide diversity. An example of this approach is disclosed in the aforementioned copending application Ser. No. 09/074224, 1998, titled “Transmitter Diversity Technique for Wireless Communications,” filed May 7.

Summary of the Invention

[0011] Improved performance is attained in an illustrative arrangement where K synchronized terminal units that transmit on N antennas to a base station having M≧K antennas, by combining interference cancellation (IC) and maximum likelihood (ML) decoding. More specifically, space-time block coding is employed in transmitters that employ N transmit antennas each, and the signals are received in a receiver that employs M receiving antennas. In accordance with the processing disclosed herein, by exploiting the structure of the space-time block code, K-1 interfering transmitting units are cancelled at the receiver, regardless of the number of transmitting antennas, N, when decoding the signals transmitted by a given mobile unit. In another embodiment of the principles of this invention, signals of a first terminal unit are decoded first, and the resulting decoded signals are employed to cancel their contribution to the signals received at the base station antennas while decoding the signals of the remaining K-1 terminal units. The process is repeated among the remaining K-1 terminal units. That is, among the remaining K-1, signals of a first terminal unit is decoded first and the resulting decoded signals are employed to cancel their contribution to the signals received at the base station antennas while decoding the signals of the remaining K-2 terminal units, and so on. This procedure is repeated K times, each time starting with decoding signals of a particular terminal unit. This successive procedure will yield additional performance improvement.

[0012] Both zero-forcing (ZF) and minimum mean-squared error (MMSE) interference cancellation (IC) and maximum likelihood (ML) techniques are disclosed.

Brief Description of the Drawing

[0013] FIG. 1 depicts an arrangement that, illustratively, includes a receiving base station (20) and two transmitting terminal units (10 and 30).

Detailed Description

[0014] FIG. 1 illustrates two transmitting units and one receiving unit that comport with the principles disclosed herein. However, this is merely illustrative, and the disclosed method is useful for more than two terminals (K>2).

[0015] Single Transmitting Unit

[0016] Transmitting unit 10 may correspond to the transmitting circuitry in a terminal unit, while receiving unit 20 may correspond to the receiving circuitry in a base station. Terminal unit 30 is shown identical to terminal unit 10. It should be understood, of course, that each terminal unit has a receiving circuit, and the base station has a transmitting circuit. The terminal units are shown to have two antennas each. Receiving unit 20 is also shown to have two receiving antennas. Here, too, it should be kept in mind that, generally, any number, M≧2, of receiving antennas can be had. Particular advantage is realized when M≧K Since the mathematical treatment below is couched in general matrix notations, the expressions are valid for any number K and/or M.

[0017] Considering terminal unit 10, the information source provides input symbols to element 13 which develops a block code. The symbols are divided into groups of two symbols each, and at a given symbol period, the two symbols in each group {c₁, c₂} are transmitted simultaneously from the two antennas. The signal transmitted from antenna 11 is c₁and the signal transmitted from antenna 12 is c₂. In the next symbol period, the signal—c₂* is transmitted from antenna 11 and the signal c₁* is transmitted from antenna 12. The symbols are modulated prior to transmission with constellation mappers 14 and 15, followed by pulse shapers 16 and 17, respectively, in a conventional manner.

[0018] In receiver 20, signals are received by antennas 21 and 22 and are applied to detector 25.

[0019] In the mathematical development of the algorithms disclosed herein, it is assumed that the channel from each of the two transmit antennas remains fixed over two consecutive symbol periods. That is,

[0020] h_i(nT)=h_i((n+1)T), i=1, 2. (1)

[0021] To ascertain the channel characteristics, the transmitter carries out a calibration session, during which pilot signals or tones are transmitted. The signals received during the calibration session are applied to channel estimator circuits 23 and 24, which are well known circuits, and the channel characteristics are thus obtained.

[0022] When only transmitter 10 is considered, the received signals at antenna 21 can be expressed as

[0023] r₁₁=h₁₁c₁+h₁₂c₂+₁ (2)

[0024] r₁₂=−h₁₁c₂*+h₁₂c₁*+η₂ (3)

[0025] where r₁₁and r₁₂are the received signals over two consecutive symbol periods, h₁₁, denotes the fading channel between transmit antenna 11 and receive antenna 21, h₁₂denotes channel between transmit antenna 12 and receive antenna 21, and η₁and η₂are noise terms, which are assumed to be complex Gaussian random variables with zero mean and power spectral density N₀/2 per dimension. Defining the vectors r=[r₁₁r₁₂*]^T, c=[c₁c₂]^T, and η=[η₁η₂*]^T, equations (2) and (3) can be rewritten in a matrix form as

[0026] R=H·c+η (4)

[0027] where the channel matrix H is defined as

\begin{matrix} H = [\begin{matrix} h_{11} & h_{12} \\ h_{12}^{*} & - h_{11}^{*} \end{matrix}] . & (5) \end{matrix}

[0028] The vector η is a complex Gaussian random vector with zero mean and covariance N₀·I. Defining C as the set of all possible symbol-pairs c={c₁,c₂}, and assuming that all symbol pairs are equi-probable, it can be easily shown that the optimum maximum likelihood (ML) decoder selects from C the symbol-pair ĉ that minimizes the expression ∥r−H·ĉ∥². This can be written as

\begin{matrix} \hat{c} = \arg \min_{\hat{c} \in C} {‖ r - H \cdot \hat{c} ‖}^{2} . & (6) \end{matrix}

[0029] It was shown by S. Alamouti in “Space Block Coding: A simple Transmitter Diversity Scheme for wireless Communications,” submitted to IEEE JSAC, September 1997 that the diversity order of the above space-time block code is equivalent to that of a two branch maximal ratio receive combining (MRRC). Because of the orthogonality of the matrix H, Alamouti also showed that this decoding rule decomposed into two separate decoding rules for c₁and c₂. The uncertainty, Δ_c, of the decoded symbols ĉ is defined as

[0030] Δ_c=∥r−H₁·ĉ∥². (7)

[0031] It should be noted that the above analysis addresses only receive antenna 21. When receiver 20 uses both antennas, i.e., antennas 21 and 22, two received signal vectors r₁and r₂can be defined for antenna 21 and 22, respectively, as

[0032] r₁=H₁·c+η₁ (8)

[0033] r₂=H₂·c+η₂ (9)

[0034] where H₁and H₂are the channel matrices to receive antennas 21 and 22, respectively, and η₁and η₂are the corresponding noise vectors. That is,

\begin{matrix} H_{1} = [\begin{matrix} h_{11} & h_{12} \\ h_{12}^{*} & - h_{11}^{*} \end{matrix}], and H_{2} = [\begin{matrix} h_{21} & h_{22} \\ h_{22}^{*} & - h_{21}^{*} \end{matrix}], & (9a) \end{matrix}

[0035] where h₂₁denotes the channel between transmit antenna 12 and receive antenna 22, and h₂₂denotes the channel between transmit antenna 11 and receive antenna 22. In this case, the ML decoding rule is

\begin{matrix} \hat{c} = \arg \min_{\hat{c} = C} ({‖ r_{1} - H_{1} \cdot \hat{c} ‖}^{2} + {‖ r_{2} - H_{2} \cdot \hat{c} ‖}^{2}), & (10) \end{matrix}

[0036] and the uncertainty of the decoded symbols is defined as

[0037] Δ_c=∥r₁−H₁·ĉ∥+∥r₂−H₂·ĉ∥². (11)

[0038] As before, both the matrices H₁and H₂are orthogonal matrices and hence the above decoding rule also decomposes to two separate decoding rules for c₁and c₂. Note that the rate of transmission (of information symbols) in the space-time block coding scheme is 1

[0039] Interference Cancellation and ML Decoding: BASIC CASE

[0040] FIG. 1, however, shows two terminal units, and the issue that needs to be addressed is the detection performance at the base station receiver when the two terminal units transmit over the same time and frequency channel.

[0041] In the notation below, g₁₁denotes the fading channel between transmit antenna 31 and receive antenna 21, g₁₂denotes the channel between antenna 31 and antenna 21, g₂₁denotes the channel between antenna 32 and antenna 22, and g₂₂denotes the channel between antenna 32 and antenna 22. Also, {c,c₂} and {s₁,s₂} denote the two symbols transmitted from terminal units 10 and 30, respectively.

[0042] At receiver 20, the received signals over two consecutive symbol periods at receive antenna 21, r₁₁and r₁₂, are

[0043] r₁₁=h₁₁c₁+h₁₂c₂+g₁₁s₁+g₁₂s₂+η₁₁ (12)

[0044] r₁₂=−h₁₁c₂*+h₁₂c₁*−g₁₁s₂*+g₁₂s₁*+η₁₂ (13)

[0045] Defining r₁=[r₁₁r₁₂*]^T, c=[c₁c₂]^T, s=[s₁s₂]^T, and n₁=[η₁₁η₁₂*]^Tequations (12) and (13) can be rewritten in matrix form as

[0046] r₁=H₁·c+G₁·s+n₁ (14)

[0047] where the channel matrices H₁and G₁between the transmitter units 10 and 30 and receive antenna 21 are given by

\begin{matrix} H_{1} = [\begin{matrix} h_{11} & h_{12} \\ h_{12}^{*} & - h_{11}^{*} \end{matrix}], and G_{1} = [\begin{matrix} g_{11} & g_{12} \\ g_{12}^{*} & - g_{11}^{*} \end{matrix}], & (15) \end{matrix}

[0048] The vector n₁is a complex Gaussian random vector with zero mean and covariance N₀·I. Similarly, the received signals over two consecutive symbol periods at receive antenna 22, r₂₁and r₂₂are

[0049] r₂₁=h₂₁c₁+h₂₂c₂+g₂₁s₁+g₂₂s₂+η₂₁ (16)

[0050] r₂₂=−h₂₁c₂*+h₂₂c₁*−g₂₁s₂*+g₂₂s₁*+η₂₂ (17)

[0051] In a similar fashion, defining r₂=[r₂₁r₂₂*]^Tand n₂=[η₂₁η₂₂*]^Tequations (16) and (17) can be rewritten as

[0052] r₂=H₂·c+G₂·s+η₂ (18)

[0053] where the channel matrices H₂and G₂between transmitter units 10 and 30 and antenna 22 are given by

\begin{matrix} H_{2} = [\begin{matrix} h_{21} & h_{22} \\ h_{22}^{*} & - h_{21}^{*} \end{matrix}], and G_{2} = [\begin{matrix} g_{21} & g_{22} \\ g_{22}^{*} & - g_{21}^{*} \end{matrix}] & (19) \end{matrix}

[0054] Equations (14) and (18) can be combined to yield the matrix form

\begin{matrix} r = [\begin{matrix} r_{1} \\ r_{2} \end{matrix}] = [\begin{matrix} H_{1} & G_{1} \\ H_{2} & G_{2} \end{matrix}] [\begin{matrix} c \\ s \end{matrix}] + [\begin{matrix} n_{1} \\ n_{2} \end{matrix}] . & (20) \end{matrix}

[0055] Zero-Forcing IC and ML Decoding Scheme: In this case, a matrix W is chosen such that

\begin{matrix} W \cdot [\begin{matrix} r_{1} \\ r_{2} \end{matrix}] = [\begin{matrix} {\tilde{r}}_{1} \\ {\tilde{r}}_{2} \end{matrix}] = [\begin{matrix} \tilde{H} & 0 \\ 0 & \tilde{G} \end{matrix}] \cdot [\begin{matrix} c \\ s \end{matrix}] + [\begin{matrix} {\tilde{n}}_{1} \\ {\tilde{n}}_{2} \end{matrix}] & (21) \end{matrix}

[0056] We can find an appropriate matrix W by realizing that

\begin{matrix} {[\begin{matrix} H_{1} & G_{1} \\ H_{2} & G_{2} \end{matrix}]}^{- 1} = [\begin{matrix} A_{1}^{- 1} & 0 \\ 0 & A_{2}^{- 1} \end{matrix}] [\begin{matrix} I & - G_{1} G_{2}^{- 1} \\ - H_{2} H_{1}^{- 1} & I \end{matrix}] & (22) \end{matrix}

[0057] where

[0058] A₁=H₁−G₁G₂⁻¹H₂and A₂=G₂−H₂H₁⁻¹G₁ (23)

[0059] Hence, if we select W as

\begin{matrix} W = [\begin{matrix} I & - G_{1} G_{2}^{- 1} \\ - H_{2} H_{1}^{- 1} & I \end{matrix}] & (24) \end{matrix}

[0060] we will have

\begin{matrix} W \cdot [\begin{matrix} r_{1} \\ r_{2} \end{matrix}] = [\begin{matrix} {\tilde{r}}_{1} \\ {\tilde{r}}_{2} \end{matrix}] = [\begin{matrix} \tilde{H} & 0 \\ 0 & \tilde{G} \end{matrix}] \cdot [\begin{matrix} c \\ s \end{matrix}] + [\begin{matrix} {\tilde{n}}_{1} \\ {\tilde{n}}_{2} \end{matrix}] & (25) \end{matrix}

[0061] where

\begin{matrix} \tilde{H} = H_{1} - G_{1} G_{2}^{- 1} H_{2} \tilde{G} = G_{2} - H_{2} H_{1}^{- 1} G_{1} . {\tilde{n}}_{1} = n_{1} - G_{1} G_{2}^{- 1} n_{2} {\tilde{n}}_{2} = n_{2} - H_{2} H_{1}^{- 1} n_{1} & (26) \end{matrix}

[0062] From equation (25) it can be easily observed that the modified received signal vector {tilde over (r)}₁contains signals only from transmitter 10 (i.e. signals from transmitter 30 have been canceled or removed) and, correspondingly, the modified received signal vector {tilde over (r)}₂contains signals only from transmitter 30 ( i.e. signals from transmitter 10 have been canceled or removed). A number of other attributes can be shown from the above to be true.

[0063] 1) the modified noise vector ñ₁is a zero mean complex Gaussian random vector with covariance

\begin{matrix} R_{n1} = N_{o} \cdot (1 + \frac{D_{g1}}{D_{g2}}) \cdot I & (27) \end{matrix}

[0064] where D_g1=|g₁₁|²+|g₁₂|²and D_g2=|g₂₁|²+|g₂₂|². Hence, the modified noise vector ñ₁is also white.

[0065] 2) the modified noise vector ñ₂is also a zero mean Gaussian random vector with covriance

\begin{matrix} R_{n2} = N_{o} \cdot (1 + \frac{D_{h2}}{D_{h1}}) \cdot I & (28) \end{matrix}

[0066] where D_h1=|h₁₁|²+|h₁₂|²and D_h2=|h₂₁|²+|h₂₂|², and hence it is also white.

[0067] 3) the matrices {tilde over (H)} and {tilde over (G)} the form

\begin{matrix} \tilde{H} = [\begin{matrix} {\tilde{h}}_{1} & {\tilde{h}}_{2} \\ {\tilde{h}}_{2}^{*} & - {\tilde{h}}_{1}^{*} \end{matrix}], and \tilde{G} = [\begin{matrix} {\tilde{g}}_{1} & {\tilde{g}}_{2} \\ {\tilde{g}}_{2}^{*} & - {\tilde{g}}_{1}^{*} \end{matrix}] . & (29) \end{matrix}

[0068] 4) Conditioned on G₁and G₂, the random variables {tilde over (h)}₁and {tilde over (h)}₂are both zero mean complex Gaussian random variables with variance σ_g²=1+D_g1/D_g2.

[0069] 5) Conditioned on H₁and H₂, the random variables {tilde over (g)}₁and {tilde over (g)}₂are both zero mean complex Gaussian random variables with variance σ_g²=1+D_h2/D_h1.

[0070] 6) the modified channel matrices {tilde over (H)} and {tilde over (G)} have a structure similar to that in equation (5), i.e. the modified channel matrices {tilde over (H)} and {tilde over (G)} are orthogonal matrices.

[0071] Considering the modified received signal vector {tilde over (r)}₁, which contains signals only from transmitter 10, i.e.,

[0072] {tilde over (r)}₁={tilde over (H)}·c+ñ₁, (30)

[0073] it is noted that this expression resembles the expression in equation (4). Hence, in accordance with the principles disclosed therein, the optimum ML decoder for the symbols transmitted by terminal unit 10 evaluates an equation that is similar to the expression of equation (6), and is given by

\begin{matrix} \hat{c} = \underset{\hat{c} \in C}{\arg \min} {‖ {\tilde{r}}_{1} - \tilde{H} \cdot \hat{c} ‖}^{2} . & (31) \end{matrix}

[0074] The corresponding decoder uncertainty is given by

[0075] Δ_c=∥{tilde over (r)}₁−{tilde over (H)}·{tilde over (c)}∥₂. (32)

[0076] Moreover, since the channel Matrix {tilde over (H)} is orthogonal, the ML decoder will also decompose into two separate rules for c₁and c₂.

[0077] In a similar fashion, considering the modified received signal vector {tilde over (r)}₂, which contains signals only from transmitter 10, i.e.,

[0078] {tilde over (r)}₂={tilde over (H)}·c+ñ₂, (33)

[0079] it is noted that this expression resembles the expression in equation (4). Hence, the optimum ML decoder for the symbols transmitted by terminal unit 30 evaluates an equation that is similar to the expression of equation (6), and is given by

\begin{matrix} \hat{s} = \underset{\hat{s} \in S}{\arg \min} {‖ {\tilde{r}}_{2} - \tilde{G} \cdot \hat{s} ‖}^{2} . & (34) \end{matrix}

[0080] The corresponding decoder uncertainty is given by

[0081] Δ_s=∥{tilde over (r)}₂−{tilde over (G)}·{tilde over (s)}∥₂. (35)

[0082] Moreover, since the channel Matrix {tilde over (G)} is also orthogonal, the ML decoder will also decompose into two separate rules for s₁and s₂.

[0083] The above-disclosed technique can be easily implemented within a detector 25 that comprises a stored program general purpose processor. Specifically, a subroutine (ĉ,Δ)=ZF.DECODE(r₁, r₂, H₁, H₂, G₁, G₂) can be installed which returns the values ĉ,Δ in response to submitted inputs r₁, r₂, H₁, H₂, G₁, and G₂, as shown below:

(\hat{c}, Δ_{c}) = ZF . DECODE (r_{1}, r_{2}, H_{1}, H_{2}, G_{1}, G_{2})

{\tilde{r} = r_{1} - G_{1} G_{2}^{- 1} r_{2} \tilde{H} = H_{1} - G_{1} G_{2}^{- 1} H_{2} \hat{c} = \underset{\hat{c} \in C}{\arg \min} {‖ \tilde{r} - \tilde{H} \cdot c ‖}^{2} Δ_{c} = {‖ \tilde{r} - \tilde{H} \cdot \hat{c} ‖}^{2}}

[0084] With such a subroutine, both ŝ and ĉ can be estimated, as follows:

[0085] (ĉ,Δ)=ZF.DECODE(r₁, r₂, H₁, H₂, G₁, G₂) (36)

[0086] (ŝ,Δ)=ZF.DECODE(r₂, r₁, G₂, G₁, H₂, H₁) (37)

[0087] It may be noted that knowledge of the noise power N₀is not required. Simulation results reveal that the performance of the FIG. 1 system which employs the principles disclosed herein and the ZF.DECODE subroutine is equivalent to that when only one terminal unit exists and the base station uses a single receive antenna. That is, the reception is equivalent to the performance of two branch MRRC diversity receiver. Advantageously, however, the disclosed technique is able to support two co-channel terminal units.

[0088] The discussion above centered on describing the technique for canceling out the signal of transmitter 10 when detecting the signal of transmitter 30, and for canceling out the signal of transmitter 30 when detecting the signal of transmitter 10. Effectively, detector 25 of receiver 20 can either comprise two processors, with one making the subroutine call of equation (31) and the other making the subroutine call of equation (32). Alternatively, the signals can be stored within detector 25 and the subroutine calls of equations 31 and 32 can be made seriatim.

[0089] Minimum Mean-Squared Error IC and ML Decoding Scheme: The above-disclosed approach for canceling the contribution of an interfering terminal unit is known as the zero-forcing (ZF) as a minimum mean-squared error technique (MMSE).

[0090] Recalling equation (20), the vector r can also be written as

\begin{matrix} r = H \cdot \tilde{c} + n where \tilde{c} = [\begin{matrix} {\tilde{c}}^{T} & {{\tilde{s}}^{T}]}^{T}, \end{matrix} r = [\begin{matrix} {\tilde{r}}_{1}^{T} & {\tilde{r}}_{2}^{T}] = {[\begin{matrix} r_{11} & r_{21}^{*} & r_{12} & r_{22}^{*} \end{matrix}]}^{T}, and \end{matrix} & (38) \\ H = [\begin{matrix} H_{1} & G_{1} \\ H_{2} & G_{2} \end{matrix}] & (39) \end{matrix}

[0091] To simplify notations, the vector r is also defined as r [r₁r₂r₃r₄]^T.

[0092] When seeking to detect and decode signals {c₁, c₂} by minimizing a mean-squared error criterion, the goal is find a linear combination of the received signals such that the mean-squared error in detecting the signals {c₁, c₂} is minimized. In general terms, this can be expressed by an error cost finction that is to be minimized, such as the function

\begin{matrix} J (α, β) = {‖ \sum_{i = 1}^{4} α_{i} r_{i} - (β_{1} c_{1} + β_{2} c_{2}) ‖}^{2} = {‖ α \cdot r - β \cdot c ‖}^{2} & (40) \end{matrix}

[0093] One may note that a minimum is certainly reached when both α and β are equal to 0, but that, of course, is not desired. Therefore, either β₁or β₂is set to 1. When β₂is set to 1, we get the following minimization criterion from equation (40)

\begin{matrix} J_{1} (α_{1}, β_{1}) = {‖ \sum_{i = 1}^{5} α_{1 i} r_{1 i} - c_{2} ‖}^{2} = {‖ {\tilde{α}}_{1} {\tilde{r}}_{r} - c_{2} ‖}^{2} & (41) \end{matrix}

[0094] where {tilde over (α)}₁=[α₁₁, α₁₂, α₁₃, α₁₄, −β₁]=[α₁−β₁] and {tilde over (r)}₁=[r^Tc₁]^T. From this it can be seen that

\begin{matrix} {\tilde{r}}_{1} = [\begin{matrix} H & 0^{T} \\ 0 & 1 \end{matrix}] [\begin{matrix} \tilde{c} \\ c_{1} \end{matrix}] + [\begin{matrix} n \\ 0 \end{matrix}] = R \cdot d + η & (42) \end{matrix}

[0095] where 0=[0 0 0 ].

[0096] What is needed is to select {tilde over (α)}₁so that the expected value of the expression in equation (41) is minimized. That is, select {tilde over (α)}×₁to minimize

[0097] E{J₁({tilde over (α)}₁)}=E{({tilde over (α)}₁{tilde over (r)}₁−c₂) ({tilde over (α)}₁{tilde over (r)}₁−c₂)*}. (43)

[0098] Taking the partial derivative with respect to {tilde over (α)}₁and setting it to zero, what results is

\begin{matrix} [\begin{matrix} M & h_{1} \\ h_{1}^{*} & 1 \end{matrix}] [\begin{matrix} α_{1}^{*} \\ - β_{1} \end{matrix}] = [\begin{matrix} h_{2} \\ 0 \end{matrix}] & (44) \end{matrix}

[0099] where M=HH* +{fraction (I/Γ)}I, Γ is the signal to noise ratio, I is the 4 by 4 identity matrix, h₁is the first column of H, and h₂is the second column of H. It follows that

[0100] α₁*=(M−h₁h₁*)⁻¹h₂and β₁*=h₁*(M−h₁h₁)⁻¹h₂. (45)

[0101] It can be shown that

\begin{matrix} {(M - h_{1} h_{1}^{*})}^{- 1} = M^{- 1} - \frac{M^{- 1} h_{1} h_{1}^{*} M^{- 1}}{1 - h_{1}^{*} M^{- 1} h_{1}}, & (46) \end{matrix}

[0102] which yields

\begin{matrix} β_{1}^{*} = \frac{h_{1}^{*} M^{- 1} h_{2}}{1 - h_{1}^{*} M^{- 1} h_{1}} & (47) \end{matrix}

[0103] From the structure of the matrix H we can easily verify that h₁and h₂are orthogonal. Using this fact and the structure of the matrix M, it can be shown that

[0104] β₁= (48)

[0105] α₁*=M⁻¹h₂ (49)

[0106] Hence, the MMSE IC solution given in equations (45) and (46) will minimize the mean-squared error in c₂without any regard to c₁. Considering the alternative cost function when β₁is set to 1, a similar analysis leads to the conclusion that

[0107] β₂= (47)

[0108] α₂*=M⁻¹h₁ (48)

[0109] In this case, the MMSE IC solution given in equations (45) and (46) will minimize the mean-squared error in c₁without any regard to c₂. Therefore, from equations (45)-(48), we can easily see that the MMSE interference canceller for signals from terminal unit 10 will consist of two different sets of weights α₁and α₂for c₂and c₁, respectively. The weights for decoding signals from terminal 30 can be obtained in a similar fashion, as expected. Thus, the decoding of signals from terminal units 10 and 30 can be performed with a single subroutine MMSE.DECODE in decoder 25 as follows:

(c, Δ_{c}) = MMSE . DECODE (r_{1}, r_{2}, H_{1}, H_{2}, G_{1}, G_{2}, Γ)

{\tilde{r} = {[\begin{matrix} r_{1}^{T} & r_{1}^{T} \end{matrix}]}^{T} \tilde{H} = [\begin{matrix} H_{1} & G_{1} \\ H_{2} & G_{2} \end{matrix}] M = {HH}^{*} + \frac{1}{Γ} I h_{1} = {[\begin{matrix} h_{11}^{T} & h_{21}^{T} \end{matrix}]}^{T} = first column of H h_{2} = {[\begin{matrix} h_{12}^{T} & h_{22}^{T} \end{matrix}]}^{T} = second column of H α_{1}^{*} = M^{- 1} h_{1}, α_{2}^{*} = M^{- 1} h_{2} c_{1} = \underset{{\hat{c}}_{1} \in C}{\arg \min} {‖ α_{1}^{*} \tilde{r} - {\hat{c}}_{1} ‖}^{2}, c_{2} = \underset{{\hat{c}}_{1} \in C}{\arg \min} {‖ α_{1}^{*} \tilde{r} - {\hat{c}}_{2} ‖}^{2} Δ_{c} = {‖ α_{1}^{*} \tilde{r} - {\hat{c}}_{1} ‖}^{2} + {‖ α_{1}^{*} \tilde{r} - {\hat{c}}_{2} ‖}^{2}}

[0110] With such a subroutine, both s and c can be estimated, as follows:

[0111] (ĉ,Δ)=MMSE.DECODE(r₁,r₂,H₁,H₂,G₁,G₂,Γ) (49)

[0112] (ŝ,Δ)=MMSE.DECODE(r₂,r₁,G₁,G₂,H₁,H₂,Γ) (50)

[0113] Similar to the zero-forcing case, simulation results reveal that the performance of the disclosed technique MMSE.DECODE is equivalent to that when only one terminal unit exists and the base station uses a single receive antenna which is equivalent to the performance of two branch MRRC diversity. However, this technique is also able to support two co-channel terminal units. In addition, when the SIR (signal-to-interference ratio, which is a ratio between the desired terminal power to the interfering terminal power) increases, the MMSE approach will have a better performance as compared to the ZF case (the ZF performance is almost the same as the performance of the MMSE approach at 0 dB SIR).

[0114] Two-Step Interference Cancellation: Actually, additional improvement can be realized by employing a two-step interference cancellation approach using either the zero-forcing or the MMSE interference cancellation techniques disclosed above. Below, we will describe this approach based on the MMSE technique. However, as one might expect there a similar approach based on the zero-forcing technique. In this two-step approach, the receiver decodes signals from both terminals using the subroutine MMSE.DECODE disclosed above. Assuming that symbols from the terminal unit 10, ĉ_o, have been decoded correctly, the receiver can, then, perfectly cancel the contribution of the terminal unit 10 in the received signal vectors r₁and r₂. The receiver then uses x₁and x₂, the received signal vectors after canceling signals from terminal unit 10, to re-decode symbols from terminal unit 30 ŝ₀using the optimum ML decoding rule in equation (10). Assuming that the symbols from terminal unit 10 has been decoded correctly, we can easily see that, the performance for terminal unit 30 will be equivalent to that with 2 transmit and 2 receive antennas (which is equivalent to 4 branch MRC diversity). If we let Δ₀=Δ_c+Δ_sdenotes the overall uncertainty for ĉ_oand ŝ_o. The receiver then repeats the above step assuming that symbols from terminal unit 30 ŝ₁have been decoded correctly using the MMSE.DECODE subroutine. As before, the receiver cancels the contribution of terminal unit 30 in the received signal vectors r₁and uses y₁and y₂, the received signal vectors after cancelling signals from terminal unit 30, to re-decode symbols from terminal unit 10 ĉ₁using the optimum ML decoding rule in equation (10). As before, assuming that symbols from terminal unit 30, the performance for terminal unit 10 will be equivalent to that with 2 transmit and 2 receive antennas. Similarly, let Δ₁=Δ_c+Δ_sdenotes the overall uncertainty for ĉ₁and ŝ₁. The receiver then compares the overall uncertainty and chooses the pair (ĉ₀, ŝ₀) if Δ₀<Δ₁and (ĉ₁, ŝ₁) otherwise. The two-step interference cancellation approach based on the MMSE technqiue disclosed above is presented below in pseudo-code subroutine II.MMSE.DECODE. As we mentioned earlier, the techniques can be also used with the zero forcing approach. Below, we also present the pseudo code subroutine II.ZF.DECODE for the two-step interference cancellation based on the zero-forcing approach.

(\hat{c}, \hat{s}) = II . MMSE . DECODE (r_{1}, r_{2}, H_{1}, H_{2}, G_{1}, G_{2}, Γ)

({\hat{c}}_{0}, Δ_{c, o}) = MMSE . DECODE (r_{1}, r_{2}, H_{1}, H_{2}, G_{1}, G_{2}, Γ)

x_{1} = r_{1} - H_{1} \cdot {\hat{c}}_{o}, x_{2} = r_{2} - H_{2} \cdot {\hat{c}}_{o}

{\hat{s}}_{o} = \underset{s \in S}{\arg \min} ({‖ x_{1} - G_{1} \cdot s ‖}^{2} + {‖ x_{2} - G_{2} \cdot s ‖}^{2})

Δ_{s, o} = {‖ x_{1} - G_{1} \cdot s ‖}^{2} + {‖ x_{2} - G_{2} \cdot s ‖}^{2}

({\hat{s}}_{1}, Δ_{s, 1}) = MMSE . DECODE (r_{1}, r_{2}, G_{1}, G_{2}, H_{1}, H_{2}, Γ)

y_{1} = r_{1} - G_{1} \cdot {\hat{s}}_{1}, {\hat{y}}_{2} = r_{2} - G_{2} \cdot {\hat{s}}_{1}

{\hat{c}}_{1} = \underset{c \in C}{\arg \min} ({‖ y_{1} - H_{1} \cdot c ‖}^{2} + {‖ y_{2} - H_{2} \cdot c ‖}^{2})

Δ_{c, 1} = {‖ y_{1} - H_{1} \cdot c ‖}^{2} + {‖ y_{2} - H_{2} \cdot c ‖}^{2}

If (Δ_{c, o} + Δ_{s, o}) < (Δ_{c, 1} + Δ_{s, 1}) \Rightarrow (\hat{c}, \hat{s}) = ({\hat{c}}_{o}, {\hat{s}}_{o})

Else (\hat{c}, \hat{s}) = ({\hat{c}}_{1}, {\hat{s}}_{1})

(\hat{c}, \hat{s}) = II . ZF . DECODE (r_{1}, r_{2}, H_{1}, H_{2}, G_{1}, G_{2})

({\hat{c}}_{0}, Δ_{c, o}) = ZF . DECODE (r_{1}, r_{2}, H_{1}, H_{2}, G_{1}, G_{2},)

x_{1} = r_{1} - H_{1} \cdot {\hat{c}}_{o}, x_{2} = r_{2} - H_{2} \cdot {\hat{c}}_{o}

{\hat{s}}_{o} = \underset{s \in S}{\arg \min} ({‖ x_{1} - G_{1} \cdot s ‖}^{2} + {‖ x_{2} - G_{2} \cdot s ‖}^{2})

Δ_{s, o} = {‖ x_{1} - G_{1} \cdot s ‖}^{2} + {‖ x_{2} - G_{2} \cdot s ‖}^{2}

({\hat{s}}_{1}, Δ_{s, 1}) = ZF . DECODE (r_{1}, r_{2}, G_{2}, G_{1}, H_{2}, H_{1},)

y_{1} = r_{1} - G_{1} \cdot {\hat{s}}_{1}, {\hat{y}}_{2} = r_{2} - G_{2} \cdot {\hat{s}}_{1}

{\hat{c}}_{1} = \underset{c \in C}{\arg \min} ({‖ y_{1} - H_{1} \cdot c ‖}^{2} + {‖ y_{2} - H_{2} \cdot c ‖}^{2})

Δ_{c, 1} = {‖ y_{1} - H_{1} \cdot c ‖}^{2} + {‖ y_{2} - H_{2} \cdot c ‖}^{2}

If (Δ_{c, o} + Δ_{s, o}) < (Δ_{c, 1} + Δ_{s, 1}) \Rightarrow (\hat{c}, \hat{s}) = ({\hat{c}}_{o}, {\hat{s}}_{o})

Else (\hat{c}, \hat{s}) = ({\hat{c}}_{1}, {\hat{s}}_{1})

[0115] Interference Cancellation and ML Decoding: GENERAL CASE

[0116] In the above basic case, we focused on the basic case where we assumed two co-channel terminals (K=2) each uses two transmit antennas (N=2). Both terminals communicate with a base station that is equipped with two transmit antennas (M=2). In this section we will consider the more general case of K≧2 co-channel terminals each is equipped with N≧2 transmitting antennas, both terminals communicate with a base that has receive M≧K antennas. We will develop similar interference cancellation and ML decoding scheme for this case.

[0117] In a paper submitted to IEEE Transactions on Information Theory, Vahid Tarokh et al. extended the above space-time block coding scheme to the case when more than two antennas are used for transmission (N≧2). There, a technique for constructing space-time block codes (with similar properties to the simple scheme described above) was developed. It was also shown that for real constellations space-time block codes with transmission rate 1 can be constructed. However, for a general complex constellation the rate of transmission for these codes will be less than 1.

[0118] In general, let us assume that the input information symbols to the transmitter are grouped into groups of Q symbols c₁, c₂, . . . , c_Q. A space-time block code in this case will map the symbols c₁, c₂, . . . , c_Qinto an N×L array C whose entries are made ±c₁, ±c₂, . . . , ±c_Qand ±c₁, ±c₂*, . . . , ±c_Q*. At time t, where 1≦t≦L, the t-th column of C is transmitted from the N antennas. In this case, the transmission rate of such code will be Q/L. In general, for a rate Q/L space-ime block code (as constructed by V. Tarokh et al.) designed for N transmit antenna, let r₁, r₂, . . . , r_Lbe the received signals at time t=1, 2, . . . , L. As before, we define the received signal vector as

[0119] r=[r₁r₂. . . r_L/2r*_L/2+1r*_L/2+2. . . r*_L]^T (51)

[0120] where the L×1 vector r can be written as

[0121] r=H·c+η (52)

[0122] and H is the L×Q channel matrix whose entries are from ±h₁, ±h₂, . . . , ±h_N, ±h*₁,±h*₂, . . . ,±h*_N, and it is an orthogonal matrix, c=[c₁c₂. . . c_Q]^T, and η is an L×1 zero-mean complex Gaussian random vector with covariance N_o·I which models the noise. The ML decoder in this case is similar to that in equation (6), that is

\begin{matrix} \hat{c} = \underset{\hat{c} \in C}{\arg \min} {‖ r - H \cdot \hat{c} ‖}^{2} & (53) \end{matrix}

[0123] and the uncertainty, Δ_c, of the decoded symbols ĉ is given by

[0124] Δ_c=∥r−H·{tilde over (c)}∥₂. (54)

[0125] As before, since the channel matrix H is orthogonal, the decoding rule in (53) decomposes into Q separate decoding rules for c₁,c₂, . . .,c_Q. For example, assuming that the terminal unit uses 4 transmit antenna, a rate {fraction (4/8)} (i.e. it is a rate ½) space-time block code is given by

[\begin{matrix} c_{1} \\ c_{2} \\ c_{3} \\ c_{4} \end{matrix}] \to [\begin{matrix} c_{1} & - c_{2} & - c_{3} & - c_{4} & c_{1}^{*} & - c_{2}^{*} & - c_{3}^{*} & - c_{4}^{*} \\ c_{2} & c_{1} & c_{4} & - c_{3} & c_{2}^{*} & c_{1}^{*} & c_{4}^{*} & - c_{3}^{*} \\ c_{3} & - c_{4} & c_{1} & c_{2} & c_{3}^{*} & - c_{4}^{*} & c_{1}^{*} & c_{2}^{*} \\ c_{4} & c_{3} & - c_{2} & c_{1} & c_{4}^{*} & c_{3}^{*} & - c_{2}^{*} & c_{1}^{*} \end{matrix}]

[0126] In this case, at time t=1 c₁, c₂, c₃, c₄are transmitted from antenna 1 through 4, respectively. At time t=2, −c₂, c₁, −c₄, c₃are transmitted from antenna 1 through 4, respectively, and so on. For this example, let r₁, r₂, . . . ,r₈be the received signals at time t=1, 2, . . . , 8. Define the received signal vector r=[r₁r₂r₃r₄r*₅r*₆r*₇r*₈]^T. In this case, we can write the received signal vector r can be written as

[0127] r=H·c+η (55)

[0128] where η is the 8×1 AWGN noise vector and H is the 8×4 channel matrix given by:

\begin{matrix} H = [\begin{matrix} h_{1} & h_{2} & h_{3} & h_{4} \\ h_{2} & - h_{1} & h_{4} & - h_{3} \\ h_{3} & - h_{4} & - h_{1} & h_{2} \\ h_{4} & h_{3} & - h_{2} & - h_{1} \\ h_{1}^{*} & h_{2}^{*} & h_{3}^{*} & h_{4}^{*} \\ h_{2}^{*} & - h_{1}^{*} & h_{4}^{*} & - h_{3}^{*} \\ h_{3}^{*} & - h_{4}^{*} & - h_{1}^{*} & h_{2}^{*} \\ h_{4}^{*} & h_{3}^{*} & - h_{2}^{*} & - h_{1}^{*} \end{matrix}] & (56) \end{matrix}

[0129] We can immediately notice that the matrix H is orthogonal, that is H*H=D_h·I where

D_{h} = \sum_{i = 1}^{4} {| h_{i} |}^{2}

[0130] and I is a 4×4 identity matrix.

[0131] Let us now assume a multi-user environment with K co-channel synchronized terminals. Each terminal uses a rate Q/L space-time block code with N transmit antenna (as constructed by V. Tarokh et al). The base station uses M≧Kantennas for reception. We can write the received signal vector at the m-th receive antenna as

\begin{matrix} r_{m} = \sum_{k = 1}^{K} H_{km} \cdot c_{k} + η_{m}, m = 1, 2, \dots, M & (57) \end{matrix}

[0132] where H_kmis the L×Q k-th user channel matrix to antenna m, c_k=[c_k1c_k2. . . c_kQ]^Tis the Q×1 information symbols vector for k-th user, and η_mis the L×1 noise vector. The entries of the k-th user channel matrix H_kmare from ±h_k,m,1, ±h_k,m,2, . . . , ±h_k,m,Nand ±h*_k,m,1, ±h*_k,m,N, where h_k,m,nis the complex channel gain between transmit antenna n of the k-th user and receive antenna m. As we stated earlier, the matrix H_kmis orthogonal.

[0133] Zero-Forcing IC and ML Decoding: Without loss of generality, let us assume that we are interested in suppressing signals from co-channel terminals 2, 3, . . . , K while decoding signals from the first terminal. This can be done in a successive manner as follows.

[0134] First, let us define r_m⁽⁰⁾=r_m. Let us assume that we start by canceling out the contributions of the K-th terminal. We can use the M-th antenna received signal vector r_Mto cancel out the contribution of the K-th terminal in the remaining M−1 received signal vectors by forming the modified received signal vectors r_m⁽¹⁾, m=1, . . . , M−1 as follows:

[0135] r_m⁽¹⁾=r_m⁽⁰⁾−H_KmH⁺_KMr_M_M⁽⁰⁾m=1, 2, . . . M−1 (58)

[0136] where H_km⁺ is the generalized inverse of the channel matrix H_kmand is given by

[0137] H_km⁺=(H*^kmH_km)⁻¹H*_km (59)

[0138] We can easily verify that H⁺_kmH_km=I, where I is the Q×Q identity matrix. We can easily verify that the modified received signal vectors r_m⁽¹⁾, m=1, . . . , M−1, do not contain any signal contribution due to the K-th user. Moreover, we can easily verify that r_m⁽¹⁾can be written as

\begin{matrix} r_{m}^{(1)} = \sum_{k = 1}^{K - 1} H_{km}^{(1)} \cdot c_{k} + η_{m}^{(1)}, m = 1, 2, \dots, M - 1 & (60) \end{matrix}

[0139] where H⁽¹⁾_kmand η⁽¹⁾_mare given by

[0140] H⁽¹⁾_km=H⁽⁰⁾_km−H⁽⁰⁾_Km(H⁽⁰⁾_KM)⁺H⁽⁰⁾_kM, m=1, 2, . . . , M−1 (61)

[0141] η⁽¹⁾_m=η⁽⁰⁾_m−H⁽⁰⁾_km(H⁽⁰⁾_KM)⁺η⁽⁰⁾_M, m=1, 2, . . . , M− (62)

[0142] Moreover, it can be shown that for codes constructed by V. Tarokh et al, the modified channel matrix H⁽¹⁾_kmwill have exactly the same structure as that of H_km. That is, the entries of the k-th user modified channel matrix H⁽¹⁾_kmare from ±h⁽¹⁾_km,1, ±h⁽¹⁾_k,m,2, . . . , ±h⁽¹⁾_k,m,Nand ±h⁽¹⁾_k,m,1, ±h^(1)*_k,m,2, . . . , ±h^(1)*_k,m,N, where h⁽¹⁾_k,m,nis the modified complex channel gain between transmit antenna n of the k-th user and receive antenna m, m=1, . . . , M−1. Hence, the modified channel matrix H⁽¹⁾_kmwill be orthogonal as well.

[0143] It can then be observed that the expression for the M−1 modified received signal vector r⁽¹⁾_min equation (60) is the same as that in equation (57) except that we now have one less interfering terminal. In a similar fashion, we can cancel out the contributions of terminal K−1 and obtain M−2 modified received signal vectors r⁽²⁾_m, m=1, . . . , M−2 that do not contain any contributions from terminals K-th and K−1. In general, after stagej, where j=1,2, . . . ,K−1 contributions from terminals K,K−1, . . . ,K−j+1 are canceled out and we are left with Mj modified received signal vectors r^(j)_m, m=1, . . . , M−j,j=1,2, . . . , K−1 that contain signals due to terminals 1, 2, . . . , K−j only. In this case, we will have

\begin{matrix} r_{m}^{(j)} = \sum_{k = 1}^{K - 1} H_{k, m}^{(j)} \cdot c_{k} + η_{m}^{(j)}, m = 1, 2, \dots, M - j & (63) \end{matrix}

[0144] where H^J_kmand η^(J)_mare given by

[0145] r^(j)_m=r^(J−1)_m−H^(J-1)_K−j,m(H^(J−1)_K−j,M−J)⁺r^(j−1)_M−jm=1, 2, . . . , M−j (64)

[0146] H^J)_k,m=H^(J−1)_k,m−H^(J−1)_K−j,m(H^J−1)_K−j,M−j)⁺H^(j−1)_k,M−1, 1≦m≦M−j,1≦k≦K−j (65)

[0147] η^(J)_m=η^(J−1)_m−H^(J−1)_K−j,m(H^(J−1)_K−j,M−J)⁺η^J−1)_M−j, m=1, 2, . . . , M−j (66)

[0148] This process is repeated until we are left with M−K+1 modified received signal vectors r^(K−1)_m, m=1, . . . , M−K+1 that contain only contributions due to the first terminal. In this case we will have

[0149] r^(K−1)_m=H^(K−1)_1,m·c₁+η^(K−1)_m, m=1, 2, . . . , M−K+1 (67)

[0150] which contains signals due to the first terminal only. Similarly, the modified channel matrix H^(K−1)_1,m, m=1, 2, . . . , M−K+1, will have a similar structure and is also orthogonal. Hence, it is straight forward to see that the ML decoder for signals from the first terminal is given by

\begin{matrix} {\hat{c}}_{1} = \underset{{\hat{c}}_{1} \in C}{\arg \min} \sum_{m = 1}^{M - K + 1} {‖ r_{m}^{(K - 1)} - H_{1, m}^{(K - 1)} \cdot {\hat{c}}_{1} ‖}^{2} & (68) \end{matrix}

[0151] and the corresponding uncertainty will be given by

\begin{matrix} Δ_{1} = \sum_{m = 1}^{M - K + 1} {‖ r_{m}^{(K - 1)} - H_{1, m}^{(K - 1)} \cdot {\hat{c}}_{1} ‖}^{2} & (69) \end{matrix}

[0152] Similarly, since the modified channel matrices H^(K−1)_1,m, 1≦m≦M−K+1 are orthogonal, as before, the decoding rule in (68) will decompose into Q separate rules for decoding c₁₁, c₁₂, . . . , c_1Q. We may observe that the basic case for zero-forcing IC and ML decoding that we discussed in detail earlier is a special case of the above approach.

({\hat{c}}_{1}, Δ_{1}) = G_ZFDECODE ({r_{m}}_{1 \leq m \leq M}, {H_{km}}_{1 \leq k \leq K, 1 \leq m \leq M})

{r_{m}^{(0)} = r_{m}, 1 \leq m \leq M H_{k, m}^{(0)} = H_{k, m}, 1 \leq m \leq M, 1 \leq k \leq K for j = 1 \to K - 1 M_{J} = M - j, K_{j} = K - j for i = 1 \to M_{j} r_{i}^{(J)} = r_{i}^{(j - 1)} - {H_{K_{j} i}^{(j - 1)} (H_{K_{j}, M_{j}}^{(j - 1)})}^{+} r_{M_{J}}^{(j - 1)} H_{k, i}^{(j)} = H_{k, i}^{(j - 1)} - {H_{K_{j}, m}^{(j - 1)} (H_{K_{j}, M_{j}}^{(j - 1)})}^{+} H_{k, M_{j}}^{(j - 1)}, 1 \leq k \leq K_{j} end end {\hat{c}}_{1} = \underset{{\hat{c}}_{1} \in C}{\arg \min} \sum_{m = 1}^{M - K + 1} {‖ r_{m}^{(K - 1)} - H_{1, m}^{(K - 1)} \cdot {\hat{c}}_{1} ‖}^{2} Δ_{1} = \sum_{m = 1}^{M - K + 1} {‖ r_{m}^{(K - 1)} - H_{1, m}^{(K - 1)} \cdot {\hat{c}}_{1} ‖}^{2}}

[0153] The above-disclosed technique can be easily implemented within a detector 25 that comprises a stored program general purpose processor. Specifically, a subroutine (c,ΔA)=G_ZFDECODE({r_m}_1≦m≦M, {H_km}_{1≦k≦K,1≦m≦M}) can be installed which returns the values c,Δ in response to submitted inputs {r_m}_1≦m≦Mand {H_km} _{1≦k≦K,1≦m≦M}, as shown above.

[0154] Minimum Mean-Squared Error IC and ML Decoding Scheme: The MMSE IC and ML decoding in the general case can be developed in a similar fashion as follows. We recall the received signal vector at the m-th receive antenna in equation (57)

\begin{matrix} r_{m} = \sum_{k = 1}^{K} H_{km} \cdot c_{k} + η_{m}, m = 1, 2, \dots, M & (70) \end{matrix}

[0155] This can be written in a matrix form as in equation (38)

[0156] r=H·{tilde over (c)}+n (71)

[0157] where r=[r^T₁r^T₂. . . r^T_M]^Tis a ML×1 vector, {tilde over (c)}=[c₁^Tc₂^T. . . c_K^T]^Tis QK×1 a vector, n=[η₁^Tη₂^T. . . η_M^T]^Tis a ML×1 vector, and

\begin{matrix} H = [\begin{matrix} H_{11} & H_{21} & \dots & H_{K1} \\ H_{12} & H_{22} & \dots & H_{K2} \\ ⋮ & ⋮ & ⋰ & ⋮ \\ H_{1 M} & H_{2 M} & \dots & H_{KM} \end{matrix}] & (72) \end{matrix}

[0158] is the ML×QK channel matrix. As before, we redefme the vector r as r=[r₁r₂. . . r_ML]^T. As before, we assume that we are interested in decoding the symbols transmitted from terminal 1 c₁₁, c₁₂, . . . , c_1Q. As before, when seeking to detect and decode signals c₁₁, c₁₂, . . . , c_1Qby minimizing a mean-squared error criterion, the goal is find a linear combination of the received signals such that the mean-squared error in detecting the signals c₁₁, c₁₂, . . . , c_1Qis minimized. In general terms, this can be expressed by an error cost function that is to be minimized, such as the function

\begin{matrix} J (α, β) = {‖ \sum_{i = 1}^{LM} α_{l} r_{l} - \sum_{j = 1}^{Q} β_{j} c_{lj} ‖}^{2} = {‖ α \cdot r - β \cdot c_{1} ‖}^{2} & (73) \end{matrix}

[0159] Similarly, as before we can see that one of the β_j, 1≦j≦Q must be set to 1 or else we get an all zero solution for α amd β. Consider the case where we set β_j=1. Hence, in this case, the criteria to be minimized is

\begin{matrix} \begin{matrix} J_{j} (α_{j}, β_{j}) = {‖ \sum_{i = 1}^{LM + Q - 1} α_{ji} r_{ji} - c_{1 j} ‖}^{2} \\ = {‖ {\tilde{α}}_{j} {\tilde{r}}_{j} - c_{1 j} ‖}^{2}, 1 \leq j \leq Q \end{matrix} & (74) \end{matrix}

[0160] where

[0161] {tilde over (α)}_j=[α_J1, α_J12. . . , α_JLM,−β₁, . . . , −β_J−1,−β_J+1, . . . ,−β_Q]=[α_j−β_j] (75)

[0162] {tilde over (r)}_J=[r^T_Jc₁₁. . . c_1j−1c_1j+1. . . c_1Q]^T (76)

[0163] If we follow the same steps as in the basic case, we arrive at the conclusion that

[0164] α_j*=M⁻1h_j,

[0165] 1≦j≦Q (78)

[0166] α_j*=M⁻¹h_j, 1≦j≦Q (78)

[0167] where h is the j-th column of the channel matrix H, and M=HH*+{fraction (1/Γ)}I, is an ML×ML matrix, Γ is the signal to noise ratio, and I is the ML×ML identity matrix.

[0168] In this case, as before, the error in decoding the j-th symbol c_1jwill be minimized without any regard to the other symbols. Hence, the MMSE-IC and ML decoder will consist of Q different combiners, one for each symbol. It should be clear now that the MMSI-IC solution for the general case is a straight forward extension to the basic case shown earlier. The MMSEBC solution for the general case can be implemented using the subroutine G₁₃MMSE.DECODE shown below.

(\hat{c}, Δ) = G_ZDECODE ({r_{m}}_{1 \leq m \leq M}, {H_{km}}_{1 \leq k \leq K, 1 \leq m \leq M})

{r = {[r_{1}^{T} r_{2}^{T} ⃛ r_{M}^{T}]}^{T} H = [\begin{matrix} H_{11} & H_{21} & \dots & H_{K1} \\ H_{12} & H_{22} & \dots & H_{K2} \\ ⋮ & ⋮ & ⋰ & ⋮ \\ H_{1 M} & H_{2 M} & \dots & H_{KM} \end{matrix}] M = {HH}^{*} + \frac{1}{Γ} I for j = 1 \to Q h_{j} = j - th column of H α_{j}^{*} = M^{- 1} h_{j} c_{j} = \underset{{\hat{c}}_{j} \in C}{\arg \min} {‖ α_{j}^{*} r - {\hat{c}}_{j} ‖}^{2} Δ_{j} = {‖ α_{j}^{*} \tilde{r} - c_{j} ‖}^{2} end \hat{c} = {[c_{1} c_{2} ⃛ c_{Q}]}^{T} Δ = \sum_{j = 1}^{Q} Δ_{j}}

Claims

We claim:

1. A method for decoding signals received at M receiving antennas, where M≧2, from K synchronized terminal units, each of which employs N transmitting antennas, comprising the steps of:

processing said M signals, received over L time intervals, where L is an integer, with signals having components related to channel coefficients between transmitting points of said terminal units and an input interface, to detect signals transmitted from each of said K terminal units by canceling interference from K−1 of said terminal units, where M and K are integers such that M≧2 and K≦M, and

identifying probable signals transmitted by said terminal units through maximum likelihood detection.

2. The method of claim 1 where transmissions of said terminal units are synchronized.

3. The method of claim 1 where each of the K terminal units employs at least two transmitting antennas, and said input interface includes M receiving antennas.

4. The method of claim 1 where each of the K terminal units employs N transmitting antennas, where N is greater than 1, and said input interface includes M receiving antennas.

5. The method of claim 4, where said interference canceling involves algebraic operations carried out on said M received signals with said signals having components related to channel coefficients between transmitting points of said terminal units, to form a plurality of signals, each of which is substantially unrelated to signals sent by all but one of the terminal units.

6. The method of claim 4 where said K=M=2 and 1≦N≦2.

7. The method of claim 6 where said processing comprises multiplying said vector of received signals by conditioning matrix

[\begin{matrix} I & - G_{1} G_{2}^{- 1} \\ - H_{2} H_{1}^{- 1} & I \end{matrix}], where

H_{1} = [\begin{matrix} h_{11} & h_{12} \\ h_{12}^{*} & - h_{11}^{*} \end{matrix}], G_{1} = [\begin{matrix} g_{11} & g_{12} \\ g_{12}^{*} & - g_{11}^{*} \end{matrix}], H_{2} = [\begin{matrix} h_{21} & h_{22} \\ h_{22}^{*} & - h_{21}^{*} \end{matrix}], G_{2} = [\begin{matrix} g_{21} & g_{22} \\ g_{22}^{*} & - g_{21}^{*} \end{matrix}],

I is the diagonal matrix, h_ijis a channel coefficient estimate between a transmitting antenna i of a first transmitting unit and a receive antenna j of said M receiving antennas, and g_ijis a channel coefficient estimate between a transmitting antenna i of a second transmitting unit and a receive antenna j of said M receiving antennas.

8. The method of claim 7 where L=2.

9. The method of claim 7 where said maximum likelihood detection minimizes a metric ∥{tilde over (r)}₁−{tilde over (H)}·ê∥²over all potential code vectors ĉ, where {tilde over (r)}₁is one of said plurality of signals formed by said step of multiplying, and {tilde over (H)}=H₁−G₁G₂⁻¹H₂.

10. The method of claim 9 where said maximum likelihood detection develops an uncertainty measure.

11. The method of claim 10 where the uncertainty measure is given by Δ_c∥{tilde over (r)}₁−{tilde over (H)}·ĉ∥².

12. The method of claim 6 where said processing executes a subroutine ZF.DECODE(r_x,r_y,H_x,H_y,G_x,G_y) which returns an output ĉ by executing the following calculations

{tilde over (r)}=r_x−G_xG_y⁻¹r_y,

{tilde over (H)}=H_x−G_xG_y⁻¹H_y,

\hat{c} = \underset{c \in C}{argmin} {‖ \tilde{r} - \tilde{H} \cdot c_{x} ‖}^{2},

where

r_xis the first argument of the ZF.DECODE subroutine and corresponds to a vector of signals received over said L time intervals at receiving antenna x,

r_yis the second argument of the ZF.DECODE subroutine and corresponds to a vector of signals received over said L time intervals at receiving antenna y,

H_xis the third argument of the ZF.DECODE subroutine and corresponds to a matrix comprising channel coefficients between transmitting antenna x of a first terminal unit and said M receiving antennas,

H_yis the fourth argument of the ZF.DECODE subroutine and corresponds to a matrix comprising channel coefficients between transmitting antenna y of said first terminal unit and said M receiving antennas,

G_xis the fifth argument of the ZF.DECODE subroutine and corresponds to a matrix comprising channel coefficients between transmitting antenna x of a second terminal unit and said M receiving antennas,

G_yis the sixth argument of the ZF.DECODE subroutine and corresponds to a matrix comprising channel coefficients between transmitting antenna y of said second terminal unit and said M receiving antennas, and

ĉ is a selected code vector c from a set of code vectors C as the code transmitted by the terminal related to the second and third arguments of the ZF.DECODE subroutine.

13. The method of claim 12 where said processing comprises executing the ZF.DECODE subroutine a first time with argument (r₁, r₂, H₁, H₂, G₁, G₂),and executing the ZF.DECODE subroutine a second time with argument (r₂, r₁, H₂, H₁, G₂, G₁).

14. The method of claim 12 where said subroutine ZF.DECODE also returns an uncertainty measure output Δ_cby carrying out the calculation Δ_c=∥{tilde over (r)}−{tilde over (H)}·ĉ∥².

15. The method of claim 14 where said processing comprises the steps of:

executing the ZF.DECODE subroutine a first time, with argument (r₁, r₂, H₁, H₂, G₁, G₂), to obtain a first version of a probable vector signal of a first terminal unit, ĉ₀, and a measure of uncertainty, Δ_c,0;

obtaining a first version of a probable vector signal of a second terminal unit, ŝ₀, and a measure of uncertainty Δ_s,0;

executing the ZF.DECODE subroutine a second time, with argument (r₂, r₁, G₂, G₁, H₂, H₁), to obtain a second version of a probable vector signal of the second terminal unit, ŝ₀, and a measure of uncertainty, Δ_s,1;

obtaining a second version of a probable vector signal of a first terminal unit, ĉ₀, and a measure of uncertainty Δ_c,1;

selecting said first version of said probable vector signals if (Δ_c,0+Δ_s,0)<(Δ_x,1+Δ_s,1), and selecting said second version of said probable vector signal based otherwise.

16. The method of claim 1 where M≧2 and K≧2, and where said processing comprises the steps of:

combining said M received signals for form a signal that is essentially independent of transmissions of all but one of the terminal units,

performing a maximum likelihood detection to identify a probably signal vector transmitted by said one of the terminal units,

modifying remaining ones of said M received signals to account for contribution of the probable signals vector, and

returning to said step of combining as long as there is a terminal unit whose transmission is to be detected, to select another one of said terminal units whose signal is to be detected with said step of maximum likelihood detection.

17. The method of claim 1 where M≧2 and K≧2, and where said processing comprises executing a subroutine

(\hat{c}, Δ) = G_ZDECODE ({r_{m}}_{1 \leq m \leq M}, {H_{km}}_{1 \leq k \leq K, 1 \leq m \leq M})

{r_{m}^{(0)} = r_{m}, 1 \leq m \leq M H_{k, m}^{(0)} = H_{k, m}, 1 \leq m \leq M, 1 \leq k \leq K for j = 1 \to K - 1 M_{j} = M - j, K_{j} = K - j for i = 1 \to M_{j} r_{i}^{(j)} = r_{i}^{(j - 1)} - {H_{j^{, i}}^{(j - 1)} (H_{K_{j}, M_{j}}^{(j - 1)})}^{+} r_{M_{j}}^{(j - 1)} H_{k, i}^{(j)} = H_{k, i}^{(j - 1)} - {H_{K_{j}, m}^{(j - 1)} (H_{K_{j}, M_{j}}^{(j - 1)})}^{+} r_{k, M_{j}}^{(j - 1)}, 1 \leq k \leq K_{j} end end \hat{c} = \underset{{\hat{c}}_{1} \in C}{argmin} \sum_{m = 1}^{M - K + 1} {‖ r_{m}^{(K - 1)} - H_{1, m}^{(K - 1)} \cdot \hat{c} ‖}^{2} Δ = \sum_{m = 1}^{M - K + 1} {‖ r_{m}^{(K - 1)} - H_{1, m}^{(K - 1)} \cdot \hat{c} ‖}^{2}}

where

r_m^(a)is a received signal at antenna m, as modified in step (a);

H_k,m^(a)is a matrix of channel coefficients between terminal unit k and receive antenna m, at step a; and

(H_k,m^(a))⁺ is the generalized inverse of H_k,m^(a).

18. The method of claim 1 where said processing form a vector of signals from said M signals and said processing includes pre-multiplying said vector by a preconditioning matrix having components related to said channel coefficients.

19. The method of claim 18 where said components of the conditioning matrix are related to estimates of said channel coefficients.

20. The method of claim 1, further comprising the step of: post-processing signals developed in said of identifying, to create detected signals adapted for delivery to users.

21. A detector responsive to M received signals, where M≧2 from a plurality of terminal units, for canceling interference from K of said terminal units, where K N, comprising:

a processor for transforming said M received signals through algebraic operations involving data related to channel coefficients between transmitting antennas of said terminal units and said M received signals, to form a plurality of signals, to each of which is substantially unrelated to signals sent by all but one of the terminal units, where M, N, and K are positive integers, and to perform maximum likelihood detection calculations on said plurality of signals;

a memory coupled to said processor for storing information for controlling said processor; and

a subroutine stored in said memory and executed by said processor, which subroutine performs said transforming of said M received signals, where each execution of said subroutine develops one of said signals that is substantially unrelated to signals sent by all but one of the terminal units.

22. A The detector of claim 21 where said subroutine is executed not less than K times.

23. The detector of claim 21 where said data is

[\begin{matrix} I & - G_{1} G_{2}^{- 1} \\ - H_{2} H_{1}^{- 1} & I \end{matrix}], where

H_{1} = [\begin{matrix} h_{11} & h_{12} \\ h_{12}^{*} & - h_{11}^{*} \end{matrix}], G_{1} = [\begin{matrix} g_{11} & g_{12} \\ g_{12}^{*} & - g_{11}^{*} \end{matrix}], H_{2} = [\begin{matrix} h_{21} & h_{22} \\ h_{22}^{*} & - h_{21}^{*} \end{matrix}], G_{2} = [\begin{matrix} g_{21} & g_{22} \\ g_{22}^{*} & - g_{21}^{*} \end{matrix}],

I is the diagonal matrix, h_ijis a channel coefficient estimate between a transmitting antenna i of a first transmitting unit and a receive antenna j of said M receiving antennas developing said M received signals, and g_ijis a channel coefficient estimate between a transmitting antenna i of a second transmitting unit and a receive antenna j of said M receiving antennas.

24. The detector of claim 23 where said processor also performs maximum likelihood detection on said plurality of signals.

Patent Citations (4)

Patent	Date	Inventor	Cited By
US5790606(A)	1998-08-01	Dent	Examiner
US5822380(A)	1998-10-01	Bottomley	Examiner
US6178196(B1)	2001-01-01	Naguib et al.	Examiner
US6188736(B1)	2001-02-01	Lo et al.	Examiner

Abstract

Description

Background of the Invention

Summary of the Invention

[0012] Both zero-forcing (ZF) and minimum mean-squared error (MMSE) interference cancellation (IC) and maximum likelihood (ML) techniques are disclosed.

Brief Description of the Drawing

[0013] FIG. 1 depicts an arrangement that, illustratively, includes a receiving base station (20) and two transmitting terminal units (10 and 30).

Detailed Description

[0015] Single Transmitting Unit

[0018] In receiver 20, signals are received by antennas 21 and 22 and are applied to detector 25.

[0020] h_i(nT)=h_i((n+1)T), i=1, 2. (1)

[0022] When only transmitter 10 is considered, the received signals at antenna 21 can be expressed as

[0023] r₁₁=h₁₁c₁+h₁₂c₂+₁ (2)

[0024] r₁₂=−h₁₁c₂*+h₁₂c₁*+η₂ (3)

[0026] R=H·c+η (4)

[0027] where the channel matrix H is defined as

\begin{matrix} H = [\begin{matrix} h_{11} & h_{12} \\ h_{12}^{*} & - h_{11}^{*} \end{matrix}] . & (5) \end{matrix}

\begin{matrix} \hat{c} = \arg \min_{\hat{c} \in C} {‖ r - H \cdot \hat{c} ‖}^{2} . & (6) \end{matrix}

[0030] Δ_c=∥r−H₁·ĉ∥². (7)

[0032] r₁=H₁·c+η₁ (8)

[0033] r₂=H₂·c+η₂ (9)

[0034] where H₁and H₂are the channel matrices to receive antennas 21 and 22, respectively, and η₁and η₂are the corresponding noise vectors. That is,

\begin{matrix} H_{1} = [\begin{matrix} h_{11} & h_{12} \\ h_{12}^{*} & - h_{11}^{*} \end{matrix}], and H_{2} = [\begin{matrix} h_{21} & h_{22} \\ h_{22}^{*} & - h_{21}^{*} \end{matrix}], & (9a) \end{matrix}

\begin{matrix} \hat{c} = \arg \min_{\hat{c} = C} ({‖ r_{1} - H_{1} \cdot \hat{c} ‖}^{2} + {‖ r_{2} - H_{2} \cdot \hat{c} ‖}^{2}), & (10) \end{matrix}

[0036] and the uncertainty of the decoded symbols is defined as

[0037] Δ_c=∥r₁−H₁·ĉ∥+∥r₂−H₂·ĉ∥². (11)

[0039] Interference Cancellation and ML Decoding: BASIC CASE

[0042] At receiver 20, the received signals over two consecutive symbol periods at receive antenna 21, r₁₁and r₁₂, are

[0043] r₁₁=h₁₁c₁+h₁₂c₂+g₁₁s₁+g₁₂s₂+η₁₁ (12)

[0044] r₁₂=−h₁₁c₂*+h₁₂c₁*−g₁₁s₂*+g₁₂s₁*+η₁₂ (13)

[0045] Defining r₁=[r₁₁r₁₂*]^T, c=[c₁c₂]^T, s=[s₁s₂]^T, and n₁=[η₁₁η₁₂*]^Tequations (12) and (13) can be rewritten in matrix form as

[0046] r₁=H₁·c+G₁·s+n₁ (14)

[0047] where the channel matrices H₁and G₁between the transmitter units 10 and 30 and receive antenna 21 are given by

\begin{matrix} H_{1} = [\begin{matrix} h_{11} & h_{12} \\ h_{12}^{*} & - h_{11}^{*} \end{matrix}], and G_{1} = [\begin{matrix} g_{11} & g_{12} \\ g_{12}^{*} & - g_{11}^{*} \end{matrix}], & (15) \end{matrix}

[0049] r₂₁=h₂₁c₁+h₂₂c₂+g₂₁s₁+g₂₂s₂+η₂₁ (16)

[0050] r₂₂=−h₂₁c₂*+h₂₂c₁*−g₂₁s₂*+g₂₂s₁*+η₂₂ (17)

[0051] In a similar fashion, defining r₂=[r₂₁r₂₂*]^Tand n₂=[η₂₁η₂₂*]^Tequations (16) and (17) can be rewritten as

[0052] r₂=H₂·c+G₂·s+η₂ (18)

[0053] where the channel matrices H₂and G₂between transmitter units 10 and 30 and antenna 22 are given by

\begin{matrix} H_{2} = [\begin{matrix} h_{21} & h_{22} \\ h_{22}^{*} & - h_{21}^{*} \end{matrix}], and G_{2} = [\begin{matrix} g_{21} & g_{22} \\ g_{22}^{*} & - g_{21}^{*} \end{matrix}] & (19) \end{matrix}

[0054] Equations (14) and (18) can be combined to yield the matrix form

\begin{matrix} r = [\begin{matrix} r_{1} \\ r_{2} \end{matrix}] = [\begin{matrix} H_{1} & G_{1} \\ H_{2} & G_{2} \end{matrix}] [\begin{matrix} c \\ s \end{matrix}] + [\begin{matrix} n_{1} \\ n_{2} \end{matrix}] . & (20) \end{matrix}

[0055] Zero-Forcing IC and ML Decoding Scheme: In this case, a matrix W is chosen such that

\begin{matrix} W \cdot [\begin{matrix} r_{1} \\ r_{2} \end{matrix}] = [\begin{matrix} {\tilde{r}}_{1} \\ {\tilde{r}}_{2} \end{matrix}] = [\begin{matrix} \tilde{H} & 0 \\ 0 & \tilde{G} \end{matrix}] \cdot [\begin{matrix} c \\ s \end{matrix}] + [\begin{matrix} {\tilde{n}}_{1} \\ {\tilde{n}}_{2} \end{matrix}] & (21) \end{matrix}

[0056] We can find an appropriate matrix W by realizing that

\begin{matrix} {[\begin{matrix} H_{1} & G_{1} \\ H_{2} & G_{2} \end{matrix}]}^{- 1} = [\begin{matrix} A_{1}^{- 1} & 0 \\ 0 & A_{2}^{- 1} \end{matrix}] [\begin{matrix} I & - G_{1} G_{2}^{- 1} \\ - H_{2} H_{1}^{- 1} & I \end{matrix}] & (22) \end{matrix}

[0057] where

[0058] A₁=H₁−G₁G₂⁻¹H₂and A₂=G₂−H₂H₁⁻¹G₁ (23)

[0059] Hence, if we select W as

\begin{matrix} W = [\begin{matrix} I & - G_{1} G_{2}^{- 1} \\ - H_{2} H_{1}^{- 1} & I \end{matrix}] & (24) \end{matrix}

[0060] we will have

\begin{matrix} W \cdot [\begin{matrix} r_{1} \\ r_{2} \end{matrix}] = [\begin{matrix} {\tilde{r}}_{1} \\ {\tilde{r}}_{2} \end{matrix}] = [\begin{matrix} \tilde{H} & 0 \\ 0 & \tilde{G} \end{matrix}] \cdot [\begin{matrix} c \\ s \end{matrix}] + [\begin{matrix} {\tilde{n}}_{1} \\ {\tilde{n}}_{2} \end{matrix}] & (25) \end{matrix}

[0061] where

\begin{matrix} \tilde{H} = H_{1} - G_{1} G_{2}^{- 1} H_{2} \tilde{G} = G_{2} - H_{2} H_{1}^{- 1} G_{1} . {\tilde{n}}_{1} = n_{1} - G_{1} G_{2}^{- 1} n_{2} {\tilde{n}}_{2} = n_{2} - H_{2} H_{1}^{- 1} n_{1} & (26) \end{matrix}

[0063] 1) the modified noise vector ñ₁is a zero mean complex Gaussian random vector with covariance

\begin{matrix} R_{n1} = N_{o} \cdot (1 + \frac{D_{g1}}{D_{g2}}) \cdot I & (27) \end{matrix}

[0064] where D_g1=|g₁₁|²+|g₁₂|²and D_g2=|g₂₁|²+|g₂₂|². Hence, the modified noise vector ñ₁is also white.

[0065] 2) the modified noise vector ñ₂is also a zero mean Gaussian random vector with covriance

\begin{matrix} R_{n2} = N_{o} \cdot (1 + \frac{D_{h2}}{D_{h1}}) \cdot I & (28) \end{matrix}

[0066] where D_h1=|h₁₁|²+|h₁₂|²and D_h2=|h₂₁|²+|h₂₂|², and hence it is also white.

[0067] 3) the matrices {tilde over (H)} and {tilde over (G)} the form

\begin{matrix} \tilde{H} = [\begin{matrix} {\tilde{h}}_{1} & {\tilde{h}}_{2} \\ {\tilde{h}}_{2}^{*} & - {\tilde{h}}_{1}^{*} \end{matrix}], and \tilde{G} = [\begin{matrix} {\tilde{g}}_{1} & {\tilde{g}}_{2} \\ {\tilde{g}}_{2}^{*} & - {\tilde{g}}_{1}^{*} \end{matrix}] . & (29) \end{matrix}

[0068] 4) Conditioned on G₁and G₂, the random variables {tilde over (h)}₁and {tilde over (h)}₂are both zero mean complex Gaussian random variables with variance σ_g²=1+D_g1/D_g2.

[0069] 5) Conditioned on H₁and H₂, the random variables {tilde over (g)}₁and {tilde over (g)}₂are both zero mean complex Gaussian random variables with variance σ_g²=1+D_h2/D_h1.

[0071] Considering the modified received signal vector {tilde over (r)}₁, which contains signals only from transmitter 10, i.e.,

[0072] {tilde over (r)}₁={tilde over (H)}·c+ñ₁, (30)

\begin{matrix} \hat{c} = \underset{\hat{c} \in C}{\arg \min} {‖ {\tilde{r}}_{1} - \tilde{H} \cdot \hat{c} ‖}^{2} . & (31) \end{matrix}

[0074] The corresponding decoder uncertainty is given by

[0075] Δ_c=∥{tilde over (r)}₁−{tilde over (H)}·{tilde over (c)}∥₂. (32)

[0076] Moreover, since the channel Matrix {tilde over (H)} is orthogonal, the ML decoder will also decompose into two separate rules for c₁and c₂.

[0077] In a similar fashion, considering the modified received signal vector {tilde over (r)}₂, which contains signals only from transmitter 10, i.e.,

[0078] {tilde over (r)}₂={tilde over (H)}·c+ñ₂, (33)

\begin{matrix} \hat{s} = \underset{\hat{s} \in S}{\arg \min} {‖ {\tilde{r}}_{2} - \tilde{G} \cdot \hat{s} ‖}^{2} . & (34) \end{matrix}

[0080] The corresponding decoder uncertainty is given by

[0081] Δ_s=∥{tilde over (r)}₂−{tilde over (G)}·{tilde over (s)}∥₂. (35)

[0082] Moreover, since the channel Matrix {tilde over (G)} is also orthogonal, the ML decoder will also decompose into two separate rules for s₁and s₂.

(\hat{c}, Δ_{c}) = ZF . DECODE (r_{1}, r_{2}, H_{1}, H_{2}, G_{1}, G_{2})

{\tilde{r} = r_{1} - G_{1} G_{2}^{- 1} r_{2} \tilde{H} = H_{1} - G_{1} G_{2}^{- 1} H_{2} \hat{c} = \underset{\hat{c} \in C}{\arg \min} {‖ \tilde{r} - \tilde{H} \cdot c ‖}^{2} Δ_{c} = {‖ \tilde{r} - \tilde{H} \cdot \hat{c} ‖}^{2}}

[0084] With such a subroutine, both ŝ and ĉ can be estimated, as follows:

[0085] (ĉ,Δ)=ZF.DECODE(r₁, r₂, H₁, H₂, G₁, G₂) (36)

[0086] (ŝ,Δ)=ZF.DECODE(r₂, r₁, G₂, G₁, H₂, H₁) (37)

[0090] Recalling equation (20), the vector r can also be written as

\begin{matrix} r = H \cdot \tilde{c} + n where \tilde{c} = [\begin{matrix} {\tilde{c}}^{T} & {{\tilde{s}}^{T}]}^{T}, \end{matrix} r = [\begin{matrix} {\tilde{r}}_{1}^{T} & {\tilde{r}}_{2}^{T}] = {[\begin{matrix} r_{11} & r_{21}^{*} & r_{12} & r_{22}^{*} \end{matrix}]}^{T}, and \end{matrix} & (38) \\ H = [\begin{matrix} H_{1} & G_{1} \\ H_{2} & G_{2} \end{matrix}] & (39) \end{matrix}

[0091] To simplify notations, the vector r is also defined as r [r₁r₂r₃r₄]^T.

\begin{matrix} J (α, β) = {‖ \sum_{i = 1}^{4} α_{i} r_{i} - (β_{1} c_{1} + β_{2} c_{2}) ‖}^{2} = {‖ α \cdot r - β \cdot c ‖}^{2} & (40) \end{matrix}

\begin{matrix} J_{1} (α_{1}, β_{1}) = {‖ \sum_{i = 1}^{5} α_{1 i} r_{1 i} - c_{2} ‖}^{2} = {‖ {\tilde{α}}_{1} {\tilde{r}}_{r} - c_{2} ‖}^{2} & (41) \end{matrix}

[0094] where {tilde over (α)}₁=[α₁₁, α₁₂, α₁₃, α₁₄, −β₁]=[α₁−β₁] and {tilde over (r)}₁=[r^Tc₁]^T. From this it can be seen that

\begin{matrix} {\tilde{r}}_{1} = [\begin{matrix} H & 0^{T} \\ 0 & 1 \end{matrix}] [\begin{matrix} \tilde{c} \\ c_{1} \end{matrix}] + [\begin{matrix} n \\ 0 \end{matrix}] = R \cdot d + η & (42) \end{matrix}

[0095] where 0=[0 0 0 ].

[0096] What is needed is to select {tilde over (α)}₁so that the expected value of the expression in equation (41) is minimized. That is, select {tilde over (α)}×₁to minimize

[0097] E{J₁({tilde over (α)}₁)}=E{({tilde over (α)}₁{tilde over (r)}₁−c₂) ({tilde over (α)}₁{tilde over (r)}₁−c₂)*}. (43)

[0098] Taking the partial derivative with respect to {tilde over (α)}₁and setting it to zero, what results is

\begin{matrix} [\begin{matrix} M & h_{1} \\ h_{1}^{*} & 1 \end{matrix}] [\begin{matrix} α_{1}^{*} \\ - β_{1} \end{matrix}] = [\begin{matrix} h_{2} \\ 0 \end{matrix}] & (44) \end{matrix}

[0099] where M=HH* +{fraction (I/Γ)}I, Γ is the signal to noise ratio, I is the 4 by 4 identity matrix, h₁is the first column of H, and h₂is the second column of H. It follows that

[0100] α₁*=(M−h₁h₁*)⁻¹h₂and β₁*=h₁*(M−h₁h₁)⁻¹h₂. (45)

[0101] It can be shown that

\begin{matrix} {(M - h_{1} h_{1}^{*})}^{- 1} = M^{- 1} - \frac{M^{- 1} h_{1} h_{1}^{*} M^{- 1}}{1 - h_{1}^{*} M^{- 1} h_{1}}, & (46) \end{matrix}

[0102] which yields

\begin{matrix} β_{1}^{*} = \frac{h_{1}^{*} M^{- 1} h_{2}}{1 - h_{1}^{*} M^{- 1} h_{1}} & (47) \end{matrix}

[0103] From the structure of the matrix H we can easily verify that h₁and h₂are orthogonal. Using this fact and the structure of the matrix M, it can be shown that

[0104] β₁= (48)

[0105] α₁*=M⁻¹h₂ (49)

[0107] β₂= (47)

[0108] α₂*=M⁻¹h₁ (48)

(c, Δ_{c}) = MMSE . DECODE (r_{1}, r_{2}, H_{1}, H_{2}, G_{1}, G_{2}, Γ)

{\tilde{r} = {[\begin{matrix} r_{1}^{T} & r_{1}^{T} \end{matrix}]}^{T} \tilde{H} = [\begin{matrix} H_{1} & G_{1} \\ H_{2} & G_{2} \end{matrix}] M = {HH}^{*} + \frac{1}{Γ} I h_{1} = {[\begin{matrix} h_{11}^{T} & h_{21}^{T} \end{matrix}]}^{T} = first column of H h_{2} = {[\begin{matrix} h_{12}^{T} & h_{22}^{T} \end{matrix}]}^{T} = second column of H α_{1}^{*} = M^{- 1} h_{1}, α_{2}^{*} = M^{- 1} h_{2} c_{1} = \underset{{\hat{c}}_{1} \in C}{\arg \min} {‖ α_{1}^{*} \tilde{r} - {\hat{c}}_{1} ‖}^{2}, c_{2} = \underset{{\hat{c}}_{1} \in C}{\arg \min} {‖ α_{1}^{*} \tilde{r} - {\hat{c}}_{2} ‖}^{2} Δ_{c} = {‖ α_{1}^{*} \tilde{r} - {\hat{c}}_{1} ‖}^{2} + {‖ α_{1}^{*} \tilde{r} - {\hat{c}}_{2} ‖}^{2}}

[0110] With such a subroutine, both s and c can be estimated, as follows:

[0111] (ĉ,Δ)=MMSE.DECODE(r₁,r₂,H₁,H₂,G₁,G₂,Γ) (49)

[0112] (ŝ,Δ)=MMSE.DECODE(r₂,r₁,G₁,G₂,H₁,H₂,Γ) (50)

(\hat{c}, \hat{s}) = II . MMSE . DECODE (r_{1}, r_{2}, H_{1}, H_{2}, G_{1}, G_{2}, Γ)

({\hat{c}}_{0}, Δ_{c, o}) = MMSE . DECODE (r_{1}, r_{2}, H_{1}, H_{2}, G_{1}, G_{2}, Γ)

x_{1} = r_{1} - H_{1} \cdot {\hat{c}}_{o}, x_{2} = r_{2} - H_{2} \cdot {\hat{c}}_{o}

{\hat{s}}_{o} = \underset{s \in S}{\arg \min} ({‖ x_{1} - G_{1} \cdot s ‖}^{2} + {‖ x_{2} - G_{2} \cdot s ‖}^{2})

Δ_{s, o} = {‖ x_{1} - G_{1} \cdot s ‖}^{2} + {‖ x_{2} - G_{2} \cdot s ‖}^{2}

({\hat{s}}_{1}, Δ_{s, 1}) = MMSE . DECODE (r_{1}, r_{2}, G_{1}, G_{2}, H_{1}, H_{2}, Γ)

y_{1} = r_{1} - G_{1} \cdot {\hat{s}}_{1}, {\hat{y}}_{2} = r_{2} - G_{2} \cdot {\hat{s}}_{1}

{\hat{c}}_{1} = \underset{c \in C}{\arg \min} ({‖ y_{1} - H_{1} \cdot c ‖}^{2} + {‖ y_{2} - H_{2} \cdot c ‖}^{2})

Δ_{c, 1} = {‖ y_{1} - H_{1} \cdot c ‖}^{2} + {‖ y_{2} - H_{2} \cdot c ‖}^{2}

If (Δ_{c, o} + Δ_{s, o}) < (Δ_{c, 1} + Δ_{s, 1}) \Rightarrow (\hat{c}, \hat{s}) = ({\hat{c}}_{o}, {\hat{s}}_{o})

Else (\hat{c}, \hat{s}) = ({\hat{c}}_{1}, {\hat{s}}_{1})

(\hat{c}, \hat{s}) = II . ZF . DECODE (r_{1}, r_{2}, H_{1}, H_{2}, G_{1}, G_{2})

({\hat{c}}_{0}, Δ_{c, o}) = ZF . DECODE (r_{1}, r_{2}, H_{1}, H_{2}, G_{1}, G_{2},)

x_{1} = r_{1} - H_{1} \cdot {\hat{c}}_{o}, x_{2} = r_{2} - H_{2} \cdot {\hat{c}}_{o}

{\hat{s}}_{o} = \underset{s \in S}{\arg \min} ({‖ x_{1} - G_{1} \cdot s ‖}^{2} + {‖ x_{2} - G_{2} \cdot s ‖}^{2})

Δ_{s, o} = {‖ x_{1} - G_{1} \cdot s ‖}^{2} + {‖ x_{2} - G_{2} \cdot s ‖}^{2}

({\hat{s}}_{1}, Δ_{s, 1}) = ZF . DECODE (r_{1}, r_{2}, G_{2}, G_{1}, H_{2}, H_{1},)

y_{1} = r_{1} - G_{1} \cdot {\hat{s}}_{1}, {\hat{y}}_{2} = r_{2} - G_{2} \cdot {\hat{s}}_{1}

{\hat{c}}_{1} = \underset{c \in C}{\arg \min} ({‖ y_{1} - H_{1} \cdot c ‖}^{2} + {‖ y_{2} - H_{2} \cdot c ‖}^{2})

Δ_{c, 1} = {‖ y_{1} - H_{1} \cdot c ‖}^{2} + {‖ y_{2} - H_{2} \cdot c ‖}^{2}

If (Δ_{c, o} + Δ_{s, o}) < (Δ_{c, 1} + Δ_{s, 1}) \Rightarrow (\hat{c}, \hat{s}) = ({\hat{c}}_{o}, {\hat{s}}_{o})

Else (\hat{c}, \hat{s}) = ({\hat{c}}_{1}, {\hat{s}}_{1})

[0115] Interference Cancellation and ML Decoding: GENERAL CASE

[0119] r=[r₁r₂. . . r_L/2r*_L/2+1r*_L/2+2. . . r*_L]^T (51)

[0120] where the L×1 vector r can be written as

[0121] r=H·c+η (52)

\begin{matrix} \hat{c} = \underset{\hat{c} \in C}{\arg \min} {‖ r - H \cdot \hat{c} ‖}^{2} & (53) \end{matrix}

[0123] and the uncertainty, Δ_c, of the decoded symbols ĉ is given by

[0124] Δ_c=∥r−H·{tilde over (c)}∥₂. (54)

[\begin{matrix} c_{1} \\ c_{2} \\ c_{3} \\ c_{4} \end{matrix}] \to [\begin{matrix} c_{1} & - c_{2} & - c_{3} & - c_{4} & c_{1}^{*} & - c_{2}^{*} & - c_{3}^{*} & - c_{4}^{*} \\ c_{2} & c_{1} & c_{4} & - c_{3} & c_{2}^{*} & c_{1}^{*} & c_{4}^{*} & - c_{3}^{*} \\ c_{3} & - c_{4} & c_{1} & c_{2} & c_{3}^{*} & - c_{4}^{*} & c_{1}^{*} & c_{2}^{*} \\ c_{4} & c_{3} & - c_{2} & c_{1} & c_{4}^{*} & c_{3}^{*} & - c_{2}^{*} & c_{1}^{*} \end{matrix}]

[0127] r=H·c+η (55)

[0128] where η is the 8×1 AWGN noise vector and H is the 8×4 channel matrix given by:

\begin{matrix} H = [\begin{matrix} h_{1} & h_{2} & h_{3} & h_{4} \\ h_{2} & - h_{1} & h_{4} & - h_{3} \\ h_{3} & - h_{4} & - h_{1} & h_{2} \\ h_{4} & h_{3} & - h_{2} & - h_{1} \\ h_{1}^{*} & h_{2}^{*} & h_{3}^{*} & h_{4}^{*} \\ h_{2}^{*} & - h_{1}^{*} & h_{4}^{*} & - h_{3}^{*} \\ h_{3}^{*} & - h_{4}^{*} & - h_{1}^{*} & h_{2}^{*} \\ h_{4}^{*} & h_{3}^{*} & - h_{2}^{*} & - h_{1}^{*} \end{matrix}] & (56) \end{matrix}

[0129] We can immediately notice that the matrix H is orthogonal, that is H*H=D_h·I where

D_{h} = \sum_{i = 1}^{4} {| h_{i} |}^{2}

[0130] and I is a 4×4 identity matrix.

\begin{matrix} r_{m} = \sum_{k = 1}^{K} H_{km} \cdot c_{k} + η_{m}, m = 1, 2, \dots, M & (57) \end{matrix}

[0135] r_m⁽¹⁾=r_m⁽⁰⁾−H_KmH⁺_KMr_M_M⁽⁰⁾m=1, 2, . . . M−1 (58)

[0136] where H_km⁺ is the generalized inverse of the channel matrix H_kmand is given by

[0137] H_km⁺=(H*^kmH_km)⁻¹H*_km (59)

\begin{matrix} r_{m}^{(1)} = \sum_{k = 1}^{K - 1} H_{km}^{(1)} \cdot c_{k} + η_{m}^{(1)}, m = 1, 2, \dots, M - 1 & (60) \end{matrix}

[0139] where H⁽¹⁾_kmand η⁽¹⁾_mare given by

[0140] H⁽¹⁾_km=H⁽⁰⁾_km−H⁽⁰⁾_Km(H⁽⁰⁾_KM)⁺H⁽⁰⁾_kM, m=1, 2, . . . , M−1 (61)

[0141] η⁽¹⁾_m=η⁽⁰⁾_m−H⁽⁰⁾_km(H⁽⁰⁾_KM)⁺η⁽⁰⁾_M, m=1, 2, . . . , M− (62)

\begin{matrix} r_{m}^{(j)} = \sum_{k = 1}^{K - 1} H_{k, m}^{(j)} \cdot c_{k} + η_{m}^{(j)}, m = 1, 2, \dots, M - j & (63) \end{matrix}

[0144] where H^J_kmand η^(J)_mare given by

[0145] r^(j)_m=r^(J−1)_m−H^(J-1)_K−j,m(H^(J−1)_K−j,M−J)⁺r^(j−1)_M−jm=1, 2, . . . , M−j (64)

[0146] H^J)_k,m=H^(J−1)_k,m−H^(J−1)_K−j,m(H^J−1)_K−j,M−j)⁺H^(j−1)_k,M−1, 1≦m≦M−j,1≦k≦K−j (65)

[0147] η^(J)_m=η^(J−1)_m−H^(J−1)_K−j,m(H^(J−1)_K−j,M−J)⁺η^J−1)_M−j, m=1, 2, . . . , M−j (66)

[0149] r^(K−1)_m=H^(K−1)_1,m·c₁+η^(K−1)_m, m=1, 2, . . . , M−K+1 (67)

\begin{matrix} {\hat{c}}_{1} = \underset{{\hat{c}}_{1} \in C}{\arg \min} \sum_{m = 1}^{M - K + 1} {‖ r_{m}^{(K - 1)} - H_{1, m}^{(K - 1)} \cdot {\hat{c}}_{1} ‖}^{2} & (68) \end{matrix}

[0151] and the corresponding uncertainty will be given by

\begin{matrix} Δ_{1} = \sum_{m = 1}^{M - K + 1} {‖ r_{m}^{(K - 1)} - H_{1, m}^{(K - 1)} \cdot {\hat{c}}_{1} ‖}^{2} & (69) \end{matrix}

({\hat{c}}_{1}, Δ_{1}) = G_ZFDECODE ({r_{m}}_{1 \leq m \leq M}, {H_{km}}_{1 \leq k \leq K, 1 \leq m \leq M})

{r_{m}^{(0)} = r_{m}, 1 \leq m \leq M H_{k, m}^{(0)} = H_{k, m}, 1 \leq m \leq M, 1 \leq k \leq K for j = 1 \to K - 1 M_{J} = M - j, K_{j} = K - j for i = 1 \to M_{j} r_{i}^{(J)} = r_{i}^{(j - 1)} - {H_{K_{j} i}^{(j - 1)} (H_{K_{j}, M_{j}}^{(j - 1)})}^{+} r_{M_{J}}^{(j - 1)} H_{k, i}^{(j)} = H_{k, i}^{(j - 1)} - {H_{K_{j}, m}^{(j - 1)} (H_{K_{j}, M_{j}}^{(j - 1)})}^{+} H_{k, M_{j}}^{(j - 1)}, 1 \leq k \leq K_{j} end end {\hat{c}}_{1} = \underset{{\hat{c}}_{1} \in C}{\arg \min} \sum_{m = 1}^{M - K + 1} {‖ r_{m}^{(K - 1)} - H_{1, m}^{(K - 1)} \cdot {\hat{c}}_{1} ‖}^{2} Δ_{1} = \sum_{m = 1}^{M - K + 1} {‖ r_{m}^{(K - 1)} - H_{1, m}^{(K - 1)} \cdot {\hat{c}}_{1} ‖}^{2}}

\begin{matrix} r_{m} = \sum_{k = 1}^{K} H_{km} \cdot c_{k} + η_{m}, m = 1, 2, \dots, M & (70) \end{matrix}

[0155] This can be written in a matrix form as in equation (38)

[0156] r=H·{tilde over (c)}+n (71)

[0157] where r=[r^T₁r^T₂. . . r^T_M]^Tis a ML×1 vector, {tilde over (c)}=[c₁^Tc₂^T. . . c_K^T]^Tis QK×1 a vector, n=[η₁^Tη₂^T. . . η_M^T]^Tis a ML×1 vector, and

\begin{matrix} H = [\begin{matrix} H_{11} & H_{21} & \dots & H_{K1} \\ H_{12} & H_{22} & \dots & H_{K2} \\ ⋮ & ⋮ & ⋰ & ⋮ \\ H_{1 M} & H_{2 M} & \dots & H_{KM} \end{matrix}] & (72) \end{matrix}

\begin{matrix} J (α, β) = {‖ \sum_{i = 1}^{LM} α_{l} r_{l} - \sum_{j = 1}^{Q} β_{j} c_{lj} ‖}^{2} = {‖ α \cdot r - β \cdot c_{1} ‖}^{2} & (73) \end{matrix}

\begin{matrix} \begin{matrix} J_{j} (α_{j}, β_{j}) = {‖ \sum_{i = 1}^{LM + Q - 1} α_{ji} r_{ji} - c_{1 j} ‖}^{2} \\ = {‖ {\tilde{α}}_{j} {\tilde{r}}_{j} - c_{1 j} ‖}^{2}, 1 \leq j \leq Q \end{matrix} & (74) \end{matrix}

[0160] where

[0161] {tilde over (α)}_j=[α_J1, α_J12. . . , α_JLM,−β₁, . . . , −β_J−1,−β_J+1, . . . ,−β_Q]=[α_j−β_j] (75)

[0162] {tilde over (r)}_J=[r^T_Jc₁₁. . . c_1j−1c_1j+1. . . c_1Q]^T (76)

[0163] If we follow the same steps as in the basic case, we arrive at the conclusion that

[0164] α_j*=M⁻1h_j,

[0165] 1≦j≦Q (78)

[0166] α_j*=M⁻¹h_j, 1≦j≦Q (78)

[0167] where h is the j-th column of the channel matrix H, and M=HH*+{fraction (1/Γ)}I, is an ML×ML matrix, Γ is the signal to noise ratio, and I is the ML×ML identity matrix.

(\hat{c}, Δ) = G_ZDECODE ({r_{m}}_{1 \leq m \leq M}, {H_{km}}_{1 \leq k \leq K, 1 \leq m \leq M})

{r = {[r_{1}^{T} r_{2}^{T} ⃛ r_{M}^{T}]}^{T} H = [\begin{matrix} H_{11} & H_{21} & \dots & H_{K1} \\ H_{12} & H_{22} & \dots & H_{K2} \\ ⋮ & ⋮ & ⋰ & ⋮ \\ H_{1 M} & H_{2 M} & \dots & H_{KM} \end{matrix}] M = {HH}^{*} + \frac{1}{Γ} I for j = 1 \to Q h_{j} = j - th column of H α_{j}^{*} = M^{- 1} h_{j} c_{j} = \underset{{\hat{c}}_{j} \in C}{\arg \min} {‖ α_{j}^{*} r - {\hat{c}}_{j} ‖}^{2} Δ_{j} = {‖ α_{j}^{*} \tilde{r} - c_{j} ‖}^{2} end \hat{c} = {[c_{1} c_{2} ⃛ c_{Q}]}^{T} Δ = \sum_{j = 1}^{Q} Δ_{j}}

Claims

We claim:

1. A method for decoding signals received at M receiving antennas, where M≧2, from K synchronized terminal units, each of which employs N transmitting antennas, comprising the steps of:

identifying probable signals transmitted by said terminal units through maximum likelihood detection.

2. The method of claim 1 where transmissions of said terminal units are synchronized.

3. The method of claim 1 where each of the K terminal units employs at least two transmitting antennas, and said input interface includes M receiving antennas.

4. The method of claim 1 where each of the K terminal units employs N transmitting antennas, where N is greater than 1, and said input interface includes M receiving antennas.

6. The method of claim 4 where said K=M=2 and 1≦N≦2.

7. The method of claim 6 where said processing comprises multiplying said vector of received signals by conditioning matrix

[\begin{matrix} I & - G_{1} G_{2}^{- 1} \\ - H_{2} H_{1}^{- 1} & I \end{matrix}], where

H_{1} = [\begin{matrix} h_{11} & h_{12} \\ h_{12}^{*} & - h_{11}^{*} \end{matrix}], G_{1} = [\begin{matrix} g_{11} & g_{12} \\ g_{12}^{*} & - g_{11}^{*} \end{matrix}], H_{2} = [\begin{matrix} h_{21} & h_{22} \\ h_{22}^{*} & - h_{21}^{*} \end{matrix}], G_{2} = [\begin{matrix} g_{21} & g_{22} \\ g_{22}^{*} & - g_{21}^{*} \end{matrix}],

8. The method of claim 7 where L=2.

10. The method of claim 9 where said maximum likelihood detection develops an uncertainty measure.

11. The method of claim 10 where the uncertainty measure is given by Δ_c∥{tilde over (r)}₁−{tilde over (H)}·ĉ∥².

12. The method of claim 6 where said processing executes a subroutine ZF.DECODE(r_x,r_y,H_x,H_y,G_x,G_y) which returns an output ĉ by executing the following calculations

{tilde over (r)}=r_x−G_xG_y⁻¹r_y,

{tilde over (H)}=H_x−G_xG_y⁻¹H_y,

\hat{c} = \underset{c \in C}{argmin} {‖ \tilde{r} - \tilde{H} \cdot c_{x} ‖}^{2},