A Fast Generalized Sphere Decoder for Optimum Decoding of Under-determined MIMO systems

Transcription

1 A Fast Generalized Sphere Decoder for Optimum Decoding of Under-determined MIMO systems Pranav Dayal and Mahesh K. Varanasi fdayalp, University of Colorado, Boulder, CO Λ Abstract A new Generalized Sphere Decoder (GSD) is proposed for optimum decoding of lattice points in a system with fewer equations than unknowns. The proposed algorithm is signicantly faster than the one proposed in [1] and is a natural extension of the sphere decoding philosophy introduced in [, 3]. The complexity reduction is achieved by partitioning the entire range of maximum likelihood hypotheses for the GSD into ordered subsets. Such a partition enables the rejection of several hypothesis at once. A parameter called the depth of the GSD is shown to offer additional complexity reduction as the gap between the number of equations and the number of unknowns grows. Numerical results are presented to show the application of the proposed GSD for various Multiple-Input Multiple-Output (MIMO) systems. 1 Introduction The advent of smart closest lattice point search algorithms such as the Pohst strategy [] or the Schnorr-Euchner strategy [4] have enabled lattice encoding as an effective means for achieving good performance on both Gaussian and Rayleigh fading channels. The sphere decoding algorithm [3] based on the Pohst strategy has found application in several communication schemes such as the ones in [5 8]. Even though the general closest lattice point search problem is NP-hard, the expected complexity of the sphere decoding algorithm has been observed to be polynomial in the number of unknowns [6, 9] for systems designed for Rayleigh fading channels with additive white Gaussian noise. However, this polynomial nature is not preserved when the system is under-determined, i. e., the vector of observed statistics lies in a space of dimension smaller than the number of unknowns. The generalized sphere decoder presented in [1] solves the closest lattice point problem with a complexity exponential in the difference of the number of equations and the number of unknowns. In this paper, we present a new generalized sphere decoding algorithm that is signicantly faster than the one in [1] and possibly a step towards achieving non-exponential complexity for decoding under-determined systems. This paper is organized as follows. The description of the existing GSD is presented in Section. The new GSD algorithm is explained in Section 3. Further improvement of the proposed strategy is presented in Section 4. Numerical results are presented in Section 5 and the conclusions are given in Section 6. Λ This work was supported in part by NSF Grant CCR and by ARO grant DADD

2 Preliminaries Since any complex linear model can be reduced to an equivalent real model, the description of the generalized sphere decoding algorithm is presented for the real model only. Consider the q-pam constellation given by I q = f odd integer jj q +1» j» q 1g, for some even integer q. Let x I M q be the transmitted vector of information symbols. Let H IR N M be the equivalent channel transfer matrix with N < M and rank(h) = N. Let the received statistics be given by y = Hx + n; (1) where n IR N is the random vector of additive white Gaussian noise. Any linear system that can be reduced to the form in (1) will be referred to as under-determined because the number of equations N is less than the number of unknowns M. The matrix H can be interpreted as the generator matrix of the projection of a random lattice to IR N and y as the vector perturbed from the transmitted lattice point Hx. Assume that H is known exactly at the receiver. The maximum-likelihood (ML) estimate ^x for x is obtained by minimizing the Euclidean distance of y from the valid lattice points, ^x = arg min ky Hxk M xiq = arg min xiq M kf(ρ x)k ; () where ρ = H T (HH T ) 1 y and F IR M M is the upper triangular matrix corresponding to the generalized Cholesky decomposition [1] of H T H. Thus, F T F = H T H and» all entries of the F1 F last M N rows of F are zero. The matrix F can be represented as F =, where 0 0 F 1 IR N N is upper triangular and F IR N (M N ). Let the subscripts G and G μ refer to the indices corresponding to the rst N and the last M N elements of a vector, respectively. The minimum distance corresponding to the optimum rule in () is now given by ky H^xk = min x μg IM N q = min x μg IM N q min x G I N q min x G I N q k[f 1 ; F ]ρ F x μ G F 1 x G k k ~ρ F 1 x G k ; (3) where ~ρ =[F 1 ; F ]ρ F x μ G in the second equation. The Generalized Sphere Decoder (GSD) in [1] checks all valid constellation points that lead to the squared Euclidean distance in () smaller than some positive number C. This is done by exhaustively searching through all hypotheses corresponding to x μ G in (3) and employing the regular sphere decoder presented in [3] to compute the bracketed expression in (3). The full rank and upper triangular nature of F 1 is used by the regular sphere decoder to place constraints on the entries of x G successively so that k ~ρ F 1 x G k» C. Such successive constraints also allowcomputingonlyas manyentries of ~ρ as required. The regular sphere decoder will nd the valid hypothesis for x G that minimizes the squared Euclidean distance between ~ρ and F 1 x G if the squared minimumdistanceis less thanc. Otherwise, a failure of the regular sphere decoder for the given hypothesis of x μ G is declared. When a hypothesis pair (x μ G; x G ) is found within the sphere, the value of C is updated and the algorithm continues to search through the remaining hypotheses for x μ G. If the regular sphere decoder fails for every possible hypothesis of x μ G, then the entire algorithm is repeated

3 with a larger value of the squared search radius C. We will refer to the GSD based on the exhaustive search over xg μ and the regular sphere decoder for every value of xg μ as the conventional GSD. The shortcoming of the conventional GSD algorithm described above can be explained as follows. The under-determined system in consideration is interpreted as consisting of N information symbols drawn from I q and one super-symbol drawn from I M N q. The conventional GSD requires one to go through each of the q M N possible hypotheses for the super-symbol at least once, irrespective of whether or not the regular sphere decoder fails subsequently for the given search radius C. Hence, the philosophy of the sphere decoding algorithm is not effectively applied to the larger super-symbol to predict a-priori the values of xg μ for which no part of the effective statistic ~ρ need be computed. In other words, the knowledge of the failure of the regular sphere decoder for a particular hypothesis of xg μ is not used to predict the failure for any other hypotheses. The new algorithm we propose next is motivated by the key observation stated above and leads to signicantly faster decoding using the ML rule than with the conventional GSD. The algorithm presented here is also a more natural generalization of the sphere decoder to systems with fewer equations than the number of unknowns. 3 A New Generalized Sphere Decoder To simplify the description of the proposed algorithm, we shall restrict to q = so that the signal constellation I q = f±1g corresponds to the BPSK set. The extension for q > is provided at the end of this section. The Euclidean distance in the ML decoding rule (), for x I M q,isgivenby kf(ρ x)k = = NX MX i=1 j=i N X 1 MX i=1 F ij (ρ j x j ) j=i F ij (ρ j x j ) 1 X M A + F Nj (ρ j x j ) (4) : (5) As in the ordinary sphere decoding algorithm, the requirement of kf(ρ x)k» C translates into a constraint on the value of each component in the summation given above. Thus, kf(ρ x)k» C ) =) p C» MX We rst make a bijective transformation given by MX F Nj (ρ j x j )» C (6) F Nj (ρ j x j )» p C: (7) b 1+x j j = ; N» j» M; (8) so that b j f0; 1g;N» j» M. Dening a j = F Nj ;N» j» M, the constraint in (7) becomes ψ MX F Nj (1 + ρ j )! p C» MX a j b j» ψ MX F Nj (1 + ρ j )! + p C: (9)

4 Let n = M N+1 and dene the vectors a =[a N ;a N +1 ;:::;a M ] T and b =[b N ;b N +1 ;:::;b M ] T. Let LB and UB represent the lower and upper bounds given in (9). Among all possible n binary sequences for the vector b, the ones that do not satisfy (9) will not lead to constellation points within the sphere of radius C. The conventional GSD effectively veries the validity of the inequality in (9) by an exhaustive search. We now propose a systematic method of checking (9) so that not all n possible binary sequences need be searched. The new method also does not sacrice the optimality inherent in the ML decoding rule. Let S be the set of all possible n binary sequences for b. Let S denote the power set of S, i. e., the set of all subsets of S. For every set A S, we associate two quantities, lb(a) and ub(a), that depend only on a and provide lower and upper bounds for the quantity a T b so that lb(a)» a T b» ub(a); 8b A: (10) Dene a relation μ H on the subsets of S based on the rule that if A; B S,thenA μ H B implies that lb(a)» lb(b). For any given value of y and C, iflb(a) is strictly greater than the upper bound UB in (9), then no sequence in A satises the constraints in (9). Additionally, if A μ H B, then there is no sequence in B that satises (9) either. The subscript H will be provided to explicitly represent the dependence of a generalized inequality or a mapping on the channel matrix H. We also introduce the symbol v H based on the upper bound ub(:) so that if A; B S,thenA v H B implies that ub(b)» ub(a). Thus, if ub(a) is strictly less than LB and A v H B, then none of the elements of either A or B will satisfy the constraint in (9). Consider the construction of a certain disjoint partition of S that leads to the same ordering with respect to μ H when all the elements of every set in the partition are permuted by a mapping ff H. Such a partition of S is required to be independent of the received statistics y, the channel matrix H and even the search radius C. Let the result of the possible orderings be summarized into the following array of D +1set inequalities, ff H (S 0;1 ) μ H ff H (S 0; ) μ H :::μ H ff H (S 0;L0 ) ff H (S 1;1 ) μ H ::: ::: μ H ff H (S 1;L1 ) ::: ::: ::: ff H (S D;1 ) μ H ff H (S D; ) :::μ H ff H (S D;L D ); (11) where [ D d=0 [L d l=1 S d;l = S and S d1 ;l 1 S d ;l = f if (d 1 ;l 1 ) 6= (d ;l ). The elements of S are now tested for the constraint in (9) starting with the zeroth row in the array, then the next until the D-th row. (ff H (S 0;1 );ff H (S 0; );:::;ff H (S 0;L0 );:::;ff H (S D;1 );:::;ff H (S D;L D )): Thus, one goes through all the row inequalities in (11) sequentially. In any row, if the lower bound lb(ff H (S d;l )) becomes greater than the upper bound in (9), then ff H (S d;l ) and all the subsequent sets to the right of ff H (S d;l ) in the same row, namely ff H (S d;l 0);l 0 >l, need not be checked and the search advances to the next row. This way we may possibly avoid searching through all the n hypotheses. The immediate question that arises is whether such a xed partition of S and a channel dependent permutation ff H is indeed possible that leads to the same ordering in (11) for all H. One solution for this problem is provided next. Since the matrix H is known at the receiver, it is possible to sort the entries of a in increasing order to form the vector a 0,sothat a 0 1» a 0» :::» a 0 n : (1)

5 Set ff H to be the permutation that sends a 0 back to a. Then, ff H (a 0 ) = a and for any binary sequence b S,wehavethat a 0T b = a T ff H (b): (13) Let w(:) denote the Hamming weight of a binary sequence. Consider the following subsets of S, S 0;1 = f0 n g; S 1;1 = fb Sjw(b) =1g; S d;l = fb Sjw(b) =d; b l =1and b j =0; 1» j» l 1g; (14) for» d» n; 1» l» L d and L d = n d +1for d. It is evident that the sets S d;l form a n l disjoint partition of S and that js d;l j = d 1 for d. The following proposition establishes the ordering that we seek in (11). Proposition 1 For the sets in (14), lb(ff H (S 1;1 )) = a 0 1 ;ub(ff H(S 1;1 )) = a 0 n and lb(ff H (S d;l )) = d+l 1 X j=l ub(ff H (S d;l )) = a 0 l + nx a 0 j ; (15) j=n d+ a 0 j ; (16) for» d» n; 1» l» L d so that ff H (S d;l1 ) μ H ff H (S d;l ) whenever d ;l 1» l. Proof: Consider d. The Hamming weight of any element b S d;l is d and thus a T ff H (b) = a 0T b is a sum of d components of a 0. This sum includes a 0 l and precludes all entries of a 0 with index strictly smaller than l. Since the vector a 0 is ordered, the smallest and largest values of such a sum are given by (15) and (16), respectively. Similarly, a 0 1 and a 0 n are lower and upper bounds for a T ff H (b); b S 1;1.Now, lb(ff H (S d;l1 ))» lb(ff H (S d;l )); 8d ;l 1» l ; (17) as can be seen by a term by term comparison of the two sides in (17) based on (15). Hence, by denition of μ H, ff H (S d;l1 ) μ H ff H (S d;l ) whenever d ;l 1» l. Before searching through all the elements of a subset ff H (S d;l ) in the order specied in (11), we rst check that both of the following conditions are satised : lb(ff H (S d;l ))» UB and ub(ff H (S d;l )) LB: (18) If the rst condition in (18) is not met then all the subsets ff H (S d;l 0), l 0 l, can be discarded. If the the second condition in (18) is not met, then the set ff H (S d;l ) can be discarded but not the ones given by ff H (S d;l 0), l 0 > l. This is because the ordering for every d is with respect to μ H and not v H. To make the sets ff H (S d;l ); 1» l» L d ; ordered with respect to v H also, the assignment of the upper bound ub(ff H (S d;l )) must be weakened to the following quantity ub(ff H (S d;l )) = nx j=n d+1 a 0 j ;» d» n; 1» l» L d: (19) The upper bound in (19) is, in fact, the maximum possible value of a T b among all binary sequences b with Hammingweightd. Theupperbound in (19) is quiteloose, but theoretically, dening lb(ff H (S d;l )) and ub(ff H (S d;l )) by (15) and (19), respectively, allows expurgation of all the sets ff H (S d;l 0) for l 0 l if any of the two conditions in (18) is not met.

6 If both the conditions in (18) are satised for the current set ff H (S d;l ), then we go through each element of ff H (S d;l ) and execute a regular sphere decoder as explained next. Let x b denote the n-length vector of BPSK symbols corresponding to the binary sequence b ff H (S d;l ). Dene the following sub-blocks of the matrix F, F 0 1 IR N 1 N 1 ; F 0 1 (i; j) =F(i; j); 1» i; j» N 1 F 0 IR N 1 M N +1 ; F 0 (i; j) =F(i; j + N 1); 1» i» N 1; 1» j» M N +1: The regular sphere decoder is executed with the effective received statistics ~ρ = [F 0 ; 1 F0 ]ρ F 0 x b, the effective generator matrix F 0 and the squared search radius set to C jp M F 1 Nj(ρ j x j )j. If the sphere decoder succeeds in nding a valid solution for x j ; 1» j» N 1, the value of C is updated to the minimum distance thus found and the corresponding lattice point is stored as the current solution. The bounds LB and UB in (9) are then updated. Since we change the quantities LB and UB while searching through the elements of ff H (S d;l ), we must verify that the conditions in (9) are met before executing the regular sphere decoder each time. After all the sequences in the set S d;l have been tested, the algorithm proceeds to verify the validity of (18) for the remaining groups in the order specied by (11). If no valid lattice point is found after either checking or discarding every set ff H (S d;l ),the value of C is increased by some factor and all the sets are searched again in the same order specied by (11). If all the sets have been checked or discarded, and the current solution is non-empty, the algorithm terminates, declaring the current solution as optimal. The proposed algorithm is now extended for q = m for some positive integer m. Replace the bijective transformation in (8) by b j = q 1+x j ; N» j» M; so that 0» b j» q 1 for N P» j» M. Consider P P m 1 the binary representation b j = i=0 i b i;j, M b i;j f0; 1g. We get that a M P m 1 jb j = i=0 (i a j )b ij. Let n = (M N +1)m. Dene the n length vectors a = [a N ; a N ;:::; m 1 a N ;:::;a M ; a M ;:::; m 1 a M ] and b = [b 0;N ;:::;b m 1;N ;:::;b 0;M ;:::;b m 1;M ] and proceed as before from (1). 4 Multiple Depth GSD In the previous section, we formulated a method for grouping n-length binary sequences into disjoint groups ff H (S d;l ) and using the information about a group, namely ub(ff H (S d;l )) and lb(ff H (S d;l )), to facilitate rejection of several sequences. For large n, the cardinality of some sets S d;l may be quite large, such as for l =1;dß n. For such sets, the upper and lower bounds [lb(ff H (S d;l ));ub(ff H (S d;l ))] may be very weak. Thus, we may not reject many sequences and the number of hypotheses tested may even be comparable to the total number of possibilities n. We now propose an extension of the method in Section 3 to accommodate for the increase in the average group size with increasing n. The multiple depth GSD is obtained by a simple recursive application of the same key idea of dividing a group of sequences into ordered subgroups. For illustration, consider the equations (11) and the partition of S given by (14). Suppose, for d 3 and 1» l» L d, both the relations in (18) are true so that all elements of the set ff H (S d;l ) need to be searched. We now further partition the set S d;l into disjoint subsets that can be ordered with respect to μ H after the permutation by ff H. Consider the following subdivision of S d;l for d 3, S d;l = [ n l d+ v=1 S d;l;v ; 1» l» n d +1 S d;l;v1 S d;l;v = f; v 1 6= v ;

7 where the subsets S d;l;v are obtained as S d;l;v = fb S d;l jb l+v =1; b l+j =0; 1» j» v 1g; (0) for 1» v» n l d +. The quantities lb and ub are dened as l+v+d X lb(ff H (S d;l;v )) = a 0 l + j=l+v ub(ff H (S d;l;v )) = a 0 l + a0 l+v + nx a 0 j ; (1) n d+3 a 0 j : () A term-by-term comparison of (1) for v = v 1 and v = v > v 1 shows that lb(ff H (S d;l;v1 ))» lb(ff H (S d;l;v )) and so the sets ff H (S d;l;v ) are indeed ordered according to μ H,i.e., ff H (S d;l;1 ) μ H ff H (S d;l; ) :::μ H ff H (S d;l;n l d+ ): All the elements of the set ff H (S d;l;v ) are checked only if the following two conditions are satised lb(ff H (S d;l;v ))» UB and ub(ff H (S d;l;v )) LB: (3) Note that lb(ff H (S d;l;1 )) = lb(ff H (S d;l )) and so lb(ff H (S d;l;1 ))» UB because (18) is true. If lb(ff H (S d;l;v )) >UBfor some v > 1, then all the sets ff H (S d;l;v 0), v 0 v can be discarded. If ub(ff H (S d;l;v )) < LB, then the set ff H (S d;l;v ) can be discarded but not ff H (S d;l;v 0) for v 0 > v because the sets ff H (S d;l;v ) foraxedd and l are not ordered according to v H. Once again we can weaken the upper bound ub to ub(ff H (S d;l;v )) = a 0 l + nx j=n d+ a 0 j ; (4) so that ub(ff H (S d;l;v )) <LBwould imply that all the sets ff H (S d;l;v ), v 0 v can be discarded. We refer to the proposed algorithm involving the partition of S into S d;l only as a Depth 1 GSD and the algorithm that additionally involves the sets S d;l;v as a Depth GSD. Further splitting of S d;l;v into disjoint ordered groups leads to higher Depth GSD schemes. In the particular partitioning scheme described above, the subsets of S consisting of sequences of Hamming weight up to X > 1 are chosen to be the same for both Depth X and Depth X 1 GSD. As the depth of the algorithm increases, the maximum cardinality of the subsets in the partition of S decreases and the upper and lower bounds lb and ub become more accurate. 5 Numerical Results In this section, we present numerical comparison of the proposed algorithm with respect to the conventional GSD. The purpose of this section is to show the potential advantage available from the best partitioning scheme for S by presenting the results for the specic partitions provided in this paper. In the rst subsection, we discuss the uncoded BLAST system with fewer receive than transmit antennas. In the second subsection, we discuss the application of the improved GSD algorithm to the Threaded Algebraic Space-Time (TAST) [7] codes. The proposed GSDs are implemented for a Depth up to 4. The Depth 1 GSD is obtained by the ordering rule in (11) and the partitioning scheme in (14). The values of lb(ff H (S d;l )) and ub(ff H (S d;l )) are taken from (15) and (16), respectively. The Depth GSD is obtained

8 from the additional partitioning in (0) and the values of lb(ff H (S d;l;v )) and ub(ff H (S d;l;v )) are taken from (1) and (), respectively. The Depths 3 and 4 GSD are obtained by the further partitioning of S d;l;v as described in Section 4. The conventional GSD is obtained from any of the above mentioned GSD by not exploiting the bounds lb and ub and exhaustively searching through all binary sequences in S. The complexity of the decoders is given in terms of the total number of CPU cycles required to perform the arithmetic operations in the algorithm. We have counted 1; 4; 7 and 0 cycles for each instance of integer addition/subtraction, floating point addition/subtraction, floating point multiplication and the square root operation, respectively. Since the sets S d;l and S d;l;v are xed irrespective of the channel realization, they are generated off-line and provided as a multi-dimensional array to the function that implements the GSD algorithm. 5.1 BLAST Scheme The system considered in this subsection is comprised of an N t transmit and N r receive antennas fading channel. Independent symbols from the 4-QAM constellation are transmitted from each antenna in every symbol interval. The received statistics in a given symbol interval for the BLAST scheme can be written as r ρ r = Ms + w; (5) N t where M IC Nr Nt is the matrix of fading coefcients, w IC Nr 1 is the noise vector and s 4 QAM Nt is the vector of transmitted symbols. The entries of M and w are independent and identically distributed zero mean complex Gaussian random variables with unit variance. In the quasi-static assumption, the fading process remains constant for T symbols periods before the next realization. The quantity ρ is the SNR per bit. The equivalent real model of Nt (5) consists of N r equations and N t integer unknowns that assume the values from f±1g. When N r N t, the regular sphere decoding algorithm is applicable. However, large number of receive antennas may not always be practical, for instance, on a space-limited mobile receiver. In such a case, ML decoding can be performed using the GSD. We compute the total number of CPU cycles taken by the conventional GSD and the proposed GSDs to decode a common sequence of channel realizations in (5). The operations involved in Cholesky factorization and other channel operations that depend only on M are not counted because these can be neglected as T!1. The initial value of the squared search radius C is chosen to be N r. If no constellation point is found within the initial sphere, the algorithm is repeated with the squared search radius increased by a factor of. An interesting trend is observed when we plot the ratio of the total CPU cycles for the proposed GSDs with respect to the total CPU cycles for the conventional GSD. In Figure 1, such a comparison is made for N r =;T =10; 4» N t» 9 and 3000 fading realizations. Not only do the proposed GSDs require fewer CPU cycles on average, but the ratio corresponding to the best of the proposed GSDs is decreasing with increasing N t. For a xed depth GSD, the ratio eventually tends to 1 for large N t as predicted in Section 4. We therefore propose to increase the depth of the GSD with increasing N t so that the rate of increase in complexity of the multiple depth GSD is strictly less than that of the conventional GSD. 5. Full Layer TAST codes Space-time codes are a means of achieving higher transmit diversity than the BLAST scheme on a quasi-static fading channel. The family of TAST codes presented in [7] are designed to

9 achieve the full transmit diversity N t and a rate of L (1» L» N t ) information symbols per channel use on a MIMO Rayleigh fading channel. The equivalent linear model for decoding each TAST codeword consists of LN t unknowns and N t N r equations. It was recommended in [7] that the TAST codes should transmit at the maximum rate of L = min(n t ;N r ) symbols per channel use to prevent the system from becoming under-determined. However, for L = N r <N t, the capacity of the equivalent channel provided by the L-layer TAST code decreases to the capacity of an N r transmit and N r receive antenna system at the same SNR. The full capacity of the N t transmit antenna and N r receive antenna system can be achieved only with the full layer (L = N t ) TAST code that unfortunately makes the system under-determined for N r < N t. Hence, the GSD algorithm is required to decode the full layer TAST code in this scenario and the size of the ML part of the conventional GSD in the equivalent real model is N t (N t N r ). The well known result of the capacity of the Rayleigh fading channel growing as min(n t ;N r ) log (SNR) [10,11] suggests that the difference in the capacities of the N t N r and N r N r systems may not be signicant for N r <N t. However, for reasonably small values of N t and N r, the difference is not negligible. For instance, the capacity at an SNR of 30 db for the (N t = 1;N r = 1) system is 9:14 bits/channel use and that for the (N t = 3;N r = 1) system is 9:71 bits/channel use. This is also reflected in the better performance of the 3-layer TAST code compared to the 1-layer TAST code for the 3 transmit antenna 1 receive antenna channel as shown in Figure. Thus, to enjoy the improved performance of the full-layer TAST code with N r <N t, one must deal with the high complexity of the generalized sphere decoder. In Figure 3, the ratio of the total numberof CPU cycles to decode 1000 codeword matrices from a 3-layer TAST code for (N t =3;N r =1)at 30 db SNR and 4-QAM symbols is shown for the Depth X GSD, 1» X» 4, with respect to the conventional GSD. The Depth 4 GSD decreases the total CPU cycles by a factor of 1:8 compared to the conventional GSD. While the Depth 4 GSD still requires about 460 times more CPU cycles than the sphere decoder for the 1-layer TAST code, it would be interesting to see how this factor for higher depth GSD or a different set partitioning scheme compares to the factor of 7 which corresponds to a cubic increase with 3 times the number of unknowns. 6 Conclusions A new generalized sphere decoding algorithm was proposed and shown to be much faster than the existing algorithm presented in [1]. The key idea of partitioning the set of all maximum likelihood hypotheses into disjoint ordered groups was applied in conjunction with the sphere decoding algorithm. It was recognized that the complexity of the proposed algorithm depends upon the largest size of the subgroup in the partition and, therefore, further improvement was achieved with the concept of multiple-depth GSD algorithm. Higher depth for the GSD algorithm are necessary to effectively reduce the complexity for severely under-determined systems. References [1] M. O. Damen, K. Abed-Meraim, and J.-C. Belore, Generalised sphere decoder for asymmetrical space time communication architecture, IEE Electronics Letters, vol. 36, no., pp , Jan. 000.

10 Ratio of Total CPU cycles N r =, 4 QAM, 8 db SNR/bit Conventional GSD Depth 1 GSD Depth GSD Depth 3 GSD Depth 4 GSD N t Figure 1: Ratio of CPU cycles with respect to the conventional GSD for BLAST. Frame Error Probability Tx, 1 Rx, 6 bits/channel use, TAST codes 1 Layer TAST code 3 Layer TAST code Average Received SNR (db) Figure : Performance of TAST codes with 3 transmit and 1 receive antenna. GSD Ratio Conv. 1 X =1 0:89 X = 0:77 X =3 0:64 X =4 0:55 Figure 3: CPU cycle ratio for full layer TAST code, N t = 3, N r = 1, 30 db SNR, Depth X vs conventional GSD. [] U. Fincke and M. Pohst, Improved methods for calculating vectors of short length in a lattice, including a complexity analysis, Math. Comput., vol. 44, pp , Apr [3] E. Viterbo and J. Boutros, A universal lattice code decoder for fading channels, IEEE Trans. Inform. Theory, vol. 45, no. 5, pp , July [4] C. P. Schnorr and M. Euchner, Lattice basis reduction: Improved practical algorithms and solving subset sum problems, Math. Programming, vol. 66, pp , [5] B. Hassibi and B. M. Hochwald, High-rate codes that are linear in space and time, IEEE Trans. Inform. Theory, vol. 48, no. 7, pp , July 00. [6] M. O. Damen, K. Abed-Meraim, and J.-C. Belore, Diagonal algebraic space time block codes, IEEE Trans. Inform. Theory, vol. 4, no. 3, pp , Mar. 00. [7] H. El Gamal and M. O. Damen, Universal space time coding, IEEE Trans. Inform. Theory, vol. 49, no. 5, pp , May 003. [8] L. Brunel and J. Boutros, Euclidean Space Lattice Decoding for Joint Detection in CDMA systems, in ITW 99, June [9] M. O. Damen, K. Abed-Meriam, and M. S. Lemdani, Further results on the sphere decoder, in Proc. IEEE Intl. Symposium on Information Theory, June 001. [10] G. J. Foschini, Layered space time architecture for wireless communication in fading environments when using multiple antennas, Bell Labs Tech. J., vol. 1, no., pp , Autumn [11] L. Zheng and D. N. C. Tse, Diversity and multiplexing: A fundamental tradeoff in multiple antenna channels, IEEE Trans. Inform. Theory, vol. 49, pp , May 003.