the scenario of @xmath0 users wanting to communicate to a single destination using cooperative strategies . for simplicity , we consider one of the @xmath0 users as the source and the remaining @xmath1 users as relays .
predominantly there are two types of relaying strategies discussed in the literature , i.e. , 1 ) amplify and forward and 2 ) decode and forward . in this work ,
we focus only on amplify and forward based relaying strategies wherein the relays are allowed to perform linear processing of the received signal . amplify and forward based relaying is simpler and moreover does not require the relays to have knowledge of the codebook used by the source or the knowledge of the source to relay channel fading gains . recently in @xcite , the idea of space - time coding for collocated mimo ( multiple input multiple output ) systems has been applied in the setting of cooperative wireless relay networks in the name of distributing space - time coding , wherein coding is performed across users and time .
this strategy provides each user a diversity order equal to the number of relays even though all the users are only equipped with a single antenna .
the diversity thus achieved is called as cooperative diversity .
however , such strategies require that the destination have complete knowledge of the fading coefficients from all the users to itself as well as that of the fading coefficients between users .
obtaining the knowledge of the fading coefficients between the users at the destination requires additional resources . to solve this problem , in @xcite , kiran _
et al _ have proposed a differential encoding / decoding setup for cooperative wireless networks that does not require the knowledge of fading coefficients between the users at the destination .
such codes were named as partially coherent distributed space - time codes in @xcite . in a recent work @xcite
, it has been shown that the same strategy of @xcite offers full diversity with a suboptimal receiver that does not require the knowledge of any of the fading coefficients . in an independent work @xcite , jing and jafarkhani have proposed a differential encoding / decoding setup for cooperative wireless relay networks which is more general than the setup proposed in @xcite and they have also provided few code constructions .
we call the class of differential space - time codes ( dstcs ) which can be used in a distributed manner for cooperative diversity as distributed dstcs ( and denote by ddstcs ) to differentiate them from dstcs for collocated mimo systems
. the problem of designing ddstcs is more challenging than that of dstcs for collocated mimo systems , since in this scenario we have additional constraints to be satisfied which are due to the cooperative diversity protocol . in @xcite ,
the setup of @xcite has been generalized to the case when the source and destination nodes have multiple antennas . in @xcite ,
the authors propose the use of a random diagonal unitary matrix codebook in general and for certain specific number of relays ( @xmath2 , @xmath3 and @xmath4 ) propose the use of cyclic codes .
however , a random code is not guaranteed to offer full diversity . moreover designing good cyclic codes
becomes difficult for large number of relays . in @xcite ,
the authors propose the use of alamouti code and @xmath5 code for the specific case of @xmath6 and @xmath7 relays .
for an arbitrary number of relays , the authors construct fully diverse ddstcs based on circulant matrices .
but , except for the alamouti code all other ddstc constructions in the literature @xcite have large decoding complexity .
thus , a general full diversity code construction targeting the requirements of low encoding complexity as well as low decoding complexity is not available .
this issue gains significant importance especially if the number of cooperating terminals is large , which is quite expected in applications like wireless sensor networks .    in this paper
we address this issue and present the following contributions :    * the notion of encoding complexity for space - time codes in terms of @xmath8-group encodability is introduced ( definition [ defn1 ] ) and its significance is highlighted .
* the differential encoding / decoding setup introduced by kiran _
et al _ @xcite , oggier _ et al _ @xcite and jing _ et al _ @xcite for wireless relay networks that use codebooks consisting of unitary matrices is generalized to allow codebooks consisting of scaled unitary matrices .
this generalization aids in reducing the encoding / decoding complexity of ddstcs .
when the codebook of scaled unitary matrices is used in the jing - hassibi protocol @xcite for differential encoding at the source , the extra conditions involving the interrelationship between relay matrices and the codeword matrices that need to be satisfied are identified ( eqn . ) . * using the algebraic framework of extended clifford algebras , a new class of ddstcs satisfying these extra conditions ( constructions [ cons_design ] and [ cons_signalset ] ) is proposed which achieves full diversity with a low complexity sub - optimal receiver ( theorem [ thm_relay ] ) .
explicit construction of multidimensional signal sets that lead to full diversity ( theorem [ thm_signalset ] ) are also provided . to the best of our knowledge ,
this is the first known low decoding complexity ddstc scheme for cooperative wireless networks which is available for all power of two number of relays .
section [ sec2 ] introduces the system model for the distributed differential encoding / decoding setup employing a differential scaled unitary matrix codebook at the source . in section [ sec3 ] , the notion of @xmath8-group encodable and @xmath8-group decodable space - time codes is discussed for the system model of section [ sec2 ] and the problem statement for ddstcs with low encoding and decoding complexity is given .
the extra conditions on the code structure imposed by the cooperative diversity protocol of @xcite are then described . in section [ sec4 ] , we explicitly construct a family of dstcs using extended clifford algebras and prove using algebraic techniques that they satisfy the extra conditions needed for ddstcs .
simulation results and discussion on further work comprise sections [ sec5 ] and [ sec6 ] respectively .      for two sets @xmath9 and @xmath10
, @xmath11 denotes the cartesian product of @xmath9 and @xmath10 . for a matrix @xmath12 ,
the matrices @xmath13 , @xmath14 and @xmath15 denote the conjugate , transpose and conjugate transpose of @xmath12 respectively . for a complex number @xmath16 , @xmath17 and @xmath18 denote its in - phase and quadrature - phase components respectively . for a vector @xmath19 ,
@xmath20 $ ] denotes the @xmath21 entry of @xmath19 .
an all zero matrix of appropriate size will be denoted by @xmath22 .
a complex gaussian vector with zero mean and covariance matrix @xmath23 will be denoted by @xmath24 .
in this section , we briefly explain the distributed differential encoding / decoding setup proposed in @xcite with a slight modification .
we consider a network consisting of a source node , a destination node and @xmath1 other relay nodes which aid the source in communicating information to the destination .
all the nodes are assumed to be equipped only with a single antenna and are half duplex constrained , i.e. , a node can not transmit and receive simultaneously at the same frequency .
the wireless channels between the terminals are assumed to quasi - static and flat fading .
the channel fading gains from the source to the @xmath25-th relay , @xmath26 and from the @xmath27-th relay to the destination @xmath28 are all assumed to be independent and identically distributed complex gaussian random variables with zero mean and unit variance .
moreover , we assume symbol synchronization among all the nodes .
every transmission cycle from the source to the destination comprises of two stages . in the first stage ,
the source transmits a @xmath29 length vector @xmath30 which the relays receive . here
, @xmath31 denotes the total average power spent by all the relays and the source . the fraction of total power @xmath31 spent by the source is denoted by @xmath32 .
the vector @xmath33 satisfies @xmath34=1 $ ] and is a function of the information that the source intends to communicate .
the exact relation of @xmath33 with the information sent by the source will be made precise in the sequel . the received vector at the @xmath27-th relay node is then given by @xmath35 . in the second half of the cycle , all the relay nodes are scheduled to transmit together .
the @xmath27-th relay node transmits a @xmath36 length vector @xmath37 which is a function of @xmath38 .
the relays are only allowed to linearly process the received vector @xmath38 or its conjugate @xmath39 . to be precise ,
the @xmath27-th relay node is equipped with a @xmath40 unitary matrix @xmath41 ( called relay matrix ) and it transmits either @xmath42 or @xmath43 . without loss of generality , we may assume that @xmath44 relays linearly process @xmath38 and the remaining @xmath45 relays linearly process @xmath39 .
here , @xmath46 denotes the fraction of total power @xmath31 spent by a relay .
if the quasi - static duration of the channel is much greater than @xmath47 time slots , then the received vector at the destination is given by @xmath48 where , @xmath49,\\ \mathbf{h}=&\left[\begin{array}{ccccccc}f_1g_1 & \dots f_mg_m & f_{m+1}^*g_{m+1 } & \dots & f_r^*g_r\end{array}\right]^t,\\ \mathbf{n}=&\sqrt{\frac{\pi_2p}{\pi_1p+1}}\left(\sum_{j=1}^{m}g_j\mathbf{a_jv_j}+\sum_{j = m+1}^{r}g_j\mathbf{a_jv_j}^*\right)+w \end{array}\ ] ] and @xmath50 represents the additive noise at the destination .
the power allocation factors @xmath32 and @xmath46 are chosen in such a way to satisfy @xmath51 . throughout this paper , we chose @xmath52 and @xmath53 as suggested in @xcite .
the differential encoding is performed at the source as follows .
an initial vector @xmath54 known to the destination which satisfies @xmath55 is transmitted by the source in the first cycle .
the initial vector @xmath56 should be such that the initial matrix @xmath57\ ] ] is unitary .
the transmitted vector at the @xmath58-th cycle is then given as follows @xmath59 where , @xmath60 is the codeword containing the information at the @xmath58-th cycle which satisfies @xmath61 .
the codebook @xmath62 at the source consisting of scaled unitary matrices should satisfy @xmath63=1 $ ] in order to meet the power constraint at the source and relays .
the originally proposed coding strategies in @xcite force @xmath64 for all codewords . the received vector at the destination in the @xmath58-th cycle can then be written as @xmath65 where , .
if @xmath66 then , substituting for @xmath67 from we get @xmath68 .
if the channel remains approximately constant for @xmath69 channel uses , then we can assume that @xmath70 .
then , @xmath71 where , @xmath72 .
now we propose to decode the codeword @xmath73 as follows @xmath74 where , @xmath75 can be estimated from the previous decision .
note that the above decoder does not require the knowledge of the channel gains and the relay matrices @xmath76 .
also , note that this decoder is not a maximum - likelihood ( ml ) decoder since the equivalent noise vector @xmath77 is dependent on the previous codeword @xmath78 through @xmath75 .
however , for the collocated mimo case , it has been shown in @xcite that there is no significant performance loss due to this . in this setup ,
we call @xmath62 a ddstc to distinguish it from dstcs for collocated mimo systems .
in this section , we describe the problem statement of designing ddstcs with low encoding complexity and low decoding complexity .    in
order to benefit from low encoding complexity and low decoding complexity , we propose to choose the ddstc @xmath62 to be a @xmath79group encodable ( @xmath80 ) linear stbc .
[ defn1]a linear design @xmath81 of size @xmath82 in @xmath83 real indeterminates or variables @xmath84 is a @xmath82 matrix with entries being a complex linear combination of the variables @xmath84 . to be precise
, it can be written as follows : @xmath85 where , @xmath86 are square complex matrices called the weight matrices .
a linear stbc @xmath62 is a finite set of @xmath82 complex matrices which can be obtained by taking a linear design @xmath81 and specifying a signal set @xmath87 from which the information vector @xmath88^t$ ] take values , with the additional condition that @xmath89 .
a linear stbc @xmath90 is said to be @xmath8-group encodable ( or @xmath91 real symbol encodable or @xmath92 complex symbol encodable ) if @xmath8 divides @xmath83 and if @xmath93 where each @xmath94 .
thus , given a @xmath8-group encodable stbc , it defines a natural partitioning of the real variables of its associated linear design into @xmath8-groups .
this partitioning defines a partitioning of the set of weight matrices of @xmath95 into @xmath8-groups , the @xmath96-th group containing @xmath97 matrices . for simplicity
we assume the simplest partitioning of the information symbol vector as @xmath98^t$ ] where , .
now @xmath95 can be written as , @xmath99 minimizing the decoding metric corresponding to , @xmath100 is in general not the same as minimizing @xmath101 for each @xmath102 individually .
however if it so happens , then the decoding complexity is reduced by a large extent .
[ defn2 ] a linear stbc @xmath90 is said to be @xmath8-group decodable ( or @xmath91 real symbol decodable or @xmath92 complex symbol decodable ) if it is @xmath8-group encodable and if its decoding metric in can be simplified as in .
@xcite the decoding metric in can be simplified as in if @xmath103 for all weight matrices @xmath104 and @xmath105 belonging to two different groups .
we illustrate the significance of the above definitions concerning encoding and decoding complexity by giving few examples .
consider the golden code for @xmath6 transmit antennas .
it has @xmath106 real variables . for the coherent mimo channel , the signal set used is qam for each complex variable . hence , this code is a @xmath106-group encodable ( since qam is a cartesian product of two pam signal sets ) and @xmath107-group decodable linear stbc . however , if we now impose the requirement for unitary codewords , then we will have to solve for signal sets which will yield unitary codewords inside the division algebra
. this may amount to entangling all the @xmath106-real variables which will make the code @xmath107-group encodable and @xmath107-group decodable .
this approach was recently attempted in @xcite .
[ eg_alamouti ] let us take the example of the alamouti code for @xmath6 transmit antennas .
it has @xmath7 real variables .
if the signal set is chosen to be psk for every complex variable , then all the codewords become unitary matrices .
hence the resulting code is @xmath6-group encodable and @xmath6-group decodable .
however , note that if we take square qam to be the signal set for each complex variable , then we get a @xmath7-group encodable ( since square qam is a cartesian product of two pam signal sets ) and @xmath7-group decodable code , but now the codewords are scaled unitary matrices as opposed to unitary matrices . thus relaxing the codewords to be scaled unitary matrices allows us to lower the encoding and decoding complexity .
the above two examples show that _ the choice of signal sets is crucial in obtaining low encoding and decoding complexity_. example [ eg_alamouti ] explicitly shows how allowing scaled unitary codebooks aids in reducing the encoding / decoding complexity .
the ddstc design problem is then to design a @xmath8-group decodable linear stbc @xmath108\right)|\mathbf{x}\in\mathscr{a}\right\}\ ] ] of size @xmath40 to be used at the source such that    1 .
all codewords are scaled unitary matrices satisfying the transmit power constraint .
2 .   the parameters @xmath83 and @xmath8 are maximized .
increasing @xmath83 is motivated by the need for high rate transmission and maximizing @xmath8 is motivated by the requirement for low decoding complexity .
[ cond_ddstc4 ] there exist @xmath1 unitary matrices @xmath109 of size @xmath40 such that the first @xmath44 of them satisfy and the remaining @xmath45 of them satisfy . 4 .
[ cond_ddstc5 ] there exists an initial vector @xmath54 such that the initial matrix @xmath110 is unitary .
@xmath111 is maximized , i.e. , the coding gain is maximized .
observe that the requirements for designing ddstcs are much more restrictive than that for dstcs .
note that condition and condition are not required for designing dstcs . as an additional requirement it would be nice to have a single design @xmath112 and a family of signal sets , one for each transmission rate such that all the required conditions are met .
this means that we need to be able to find @xmath1 relay matrices satisfying the required conditions irrespective of the size of the code @xmath113 .
in this section we give an explicit construction of full diversity @xmath7-group decodable ddstcs for all power of two number of relays . optimizing the coding gain for the proposed ddstcs is difficult in general and is beyond the scope of this paper .
constructing @xmath79group decodable ddstcs for arbitrary @xmath8 appears to be rather difficult and a brief discussion on the issues involved is given is section [ sec6 ] . in subsection
[ subsection_design ] , the construction of the proposed linear designs using tools from extended clifford algebras is given .
in subsection [ subsection_signalset ] , an explicit construction of signal sets @xmath114 leading to full diversity for arbitrary transmission rate is provided for the linear designs given in subsection [ subsection_design ] , thus completely describing the construction of the linear stbc @xmath62 . in subsection
[ subsection_relay ] , we use algebraic techniques to explicitly construct the @xmath1 relay matrices satisfying the required conditions for arbitrary size of the codebook @xmath62 .      in this subsection , we briefly describe a construction of a class of rate one , linear designs satisfying the conditions for four - group decodability which were first obtained using extended clifford algebras in @xcite . extended clifford algebras and linear designs from them are addressed in detail in @xcite .
however @xcite chooses signal sets as applicable for coherent communication in a distributed space - time coding setup .
a brief explanation of extended clifford algebras and linear designs using them is provided in appendix [ appendix1 ] . this algebraic framework given in appendix [ appendix1 ]
is needed to explicitly construct the relay matrices and the initial vector in subsection [ subsection_relay ] .    given a @xmath82 linear design @xmath115 in @xmath116 complex variables @xmath117 , one can construct a @xmath118 linear design @xmath119 as follows :    @xmath120 } \ ] ]    where , the linear design @xmath121 is identical to the linear design @xmath115 except that it is a design in a different set of complex variables @xmath122 .
we call this construction as the abba construction which was first introduced in @xcite .    [ cons_doubling ] given a @xmath82 linear design @xmath115
, one can construct a @xmath118 linear design @xmath123 as follows .
@xmath124 } \ ] ]    where , the linear design @xmath121 is identical to the linear design @xmath115 .
we call this construction as the doubling construction.    we are now ready to describe the construction of the proposed linear designs .
[ cons_design ] to obtain a linear design for @xmath125 relays , follow the steps given below .
+ step 1 : starting with the linear design @xmath126 $ ] in complex variable @xmath127 , keep applying abba construction iteratively on it till a @xmath128 linear design @xmath119 is obtained .
+ step 2 : apply doubling construction on @xmath119 and scale it by @xmath129 to obtain the linear design .
following the steps given above for @xmath130 relays , we get    @xmath131.\ ] ]      in this subsection , we construct signal sets for the linear designs constructed in the previous subsection that lead to full diversity stbcs . in general , the signal sets should be designed such that the resulting linear stbc @xmath62 meets the following requirements .    1
.   scaled unitary codeword matrices meeting power constraint .
four - group encodable and four - group decodable .
3 .   difference of any two different codeword matrices should be full rank .
we call such a code to be fully diverse. 4
.   the minimum determinant of the difference of any two codewords matrix should be maximized ( coding gain ) .
we shall first illustrate the signal set construction procedure for @xmath7 relays and derive important insights form it .
then , we generalize the ideas thus obtained for constructing signal sets for any @xmath125 relays .
the design for @xmath7 relays according to the construction in the previous subsection is given by @xmath132.\ ] ] let us look at @xmath133 shown in at the top of the next page to find out the conditions on the signal sets under which the codewords are scaled unitary matrices .
@xmath134\ ] ]    ' '' ''    from , we see that the signal set should be chosen such that the following condition is satisfied for all the signal points .
@xmath135 firstly , we identify the grouping of the real variables into groups such that the corresponding weight matrices in each group satisfy . according to the construction , the four groups or real variables are as follows : first group : @xmath136 , second group : @xmath137 , third group : @xmath138 , fourth group :
it is important to choose signal sets such that they do not enforce joint constraints on variables from different groups .
the requirement for scaled unitary codewords as in can be satisfied without disturbing @xmath7-group encodability if all the signal points satisfy the following equation : @xmath140 where , @xmath141 and @xmath142 are positive real constants .
then , the average power constraint requirement can be met by satisfying the conditions , @xmath143 a common set of solutions for and can be obtained by taking points on the intersection of circles and hyperbolas . to meet the third requirement of full diversity we make use of the structure of the constructed designs .
note that    @xmath144\right|\\ \\ & = \frac{1}{4}\left|\left[\begin{array}{cc}\delta \mathbf{a}^h\delta \mathbf{a}+\delta \mathbf{b}^h\delta \mathbf{b } & \mathbf{0}\\ \mathbf{0 } & \delta \mathbf{a}\delta \mathbf{a}^h+\delta\mathbf{b}\delta \mathbf{b}^h\end{array}\right]\right|\\ & \geq \frac{1}{4}\mathrm{max}(|\delta \mathbf{a}|^*|\delta \mathbf{a}|,|\delta \mathbf{b}|^*|\delta \mathbf{b}|)^2 \end{array } } \ ] ]    where , the notation @xmath145 has been used to denote the difference matrix which will have the same form as the associated linear design .
the second equality in follows because the linear designs @xmath146 and @xmath147 commute .
this can be proved as follows .
let @xmath148 $ ] and .
then , we have    @xmath149\left[\begin{array}{cc}\mathbf{w ' } & \mathbf{x'}\\\mathbf{x ' } & \mathbf{w ' } \end{array}\right]=\left[\begin{array}{cc}\mathbf{ww'}+\mathbf{xx ' } & \mathbf{wx'+xw'}\\\mathbf{xw'+wx ' } & \mathbf{xx'+ww ' } \end{array}\right],\\ \\ \left[\begin{array}{cc}\mathbf{w ' } & \mathbf{x'}\\\mathbf{x ' } & \mathbf{w ' } \end{array}\right]\left[\begin{array}{cc}\mathbf{w } & \mathbf{x}\\\mathbf{x } & \mathbf{w } \end{array}\right]=\left[\begin{array}{cc}\mathbf{w'w}+\mathbf{x'x } & \mathbf{w'x}+\mathbf{x'w}\\\mathbf{x'w}+\mathbf{w'x } & \mathbf{x'x}+\mathbf{w'w } \end{array}\right ] .
\end{array } } \ ] ]    thus we observe that linear designs @xmath146 and @xmath147 commute if their constituent sub - matrices commute .
applying this argument recursively , the claim follows .
thus @xmath123 is fully diverse , i.e. , @xmath150 if either of the submatrices @xmath146 or @xmath147 is fully diverse .
for the example of @xmath7 relays , we have .
hence we can guarantee full diversity by ensuring that @xmath151 and @xmath152 .
just like before , we should be careful not to disturb @xmath7-group encodability in the process .
we take care of that requirement also by satisfying the following conditions :    @xmath153        putting together all the conditions that need to be met , we have to choose constellation points satisfying , and .
the solution can be found simply by finding the intersection of points on the unit circle @xmath154 with a hyperbola @xmath155 , where @xmath156 on the two dimensional @xmath157 plane .
this is illustrated in fig .
[ fig_signal2dimstruct ] .    observe from fig .
[ fig_signal2dimstruct ] that the hyperbola intersects the circle at four different points .
but the full diversity criterion demands that @xmath158 . after enforcing this condition ,
only two points survive out of the four points
. they can be either the set of points marked a or the set of points marked b in fig .
[ fig_signal2dimstruct ] .
thus we have obtained a signal set containing @xmath6 points .
supposing that we need more than two points , we can draw more circles ( centered at origin ) with radii such that the average power constraint is met and then find those points intersecting with the hyperbola .
more precisely , to get @xmath159 points , we draw @xmath160 concentric circles with increasing radii @xmath161 such that @xmath162 .
then we find those points intersecting with the hyperbola @xmath155 where , @xmath163 is a positive number less than @xmath164 . in this manner
we can get the desired signal set for the variables @xmath165 and @xmath166 . the signal set for the variables @xmath167 and @xmath168 can be obtained by considering a different hyperbola .
this is illustrated in fig .
[ fig_signal2dim ] shown at the top of the next page .        based on empirical studies ,
we propose to choose points in the set * a * and set * b * alternatively on each circle . also , choosing the constants @xmath169 seems to give the largest coding gain .
similar observations have also been made in @xcite for another class of linear stbcs applicable in co - located mimo systems .
implementing these two suggestions makes the signal identical for all the @xmath7-groups and the resulting signal set for the case of @xmath106 points is shown in fig.[fig_signal2dim2 ] .
note that the proposed signal set construction is far from general and the most general solution would correspond to the following equations : @xmath170        in general , we need to choose the parameters @xmath171 and the solution points such that the coding gain is maximized .    based on the detailed illustration of signal set construction for @xmath7 relays antennas , a natural generalization of it to higher dimensions is given in construction [ cons_signalset ] as follows .
[ cons_signalset ] suppose that we want a @xmath172 points signal set@xmath173 for the constructed design for @xmath125 relays .
then the resulting signal set @xmath173 should be a cartesian product of @xmath7 signal sets in @xmath174 since we insist on @xmath7-group encodability . in our case , we choose all the four sets to be identical and each will contain @xmath175{q}$ ] points .
let the signal points in @xmath174 be labeled as @xmath176{q}$ ] .
if @xmath177 , for some integers @xmath178 and @xmath179 where @xmath180 , then @xmath181 is given by @xmath182=0\ \forall j\neq(q\ \mathrm{mod}\ 2^{\lambda-1})+1\\ \mathbf{p}_i[(q\ \mathrm{mod}\ 2^{\lambda-1})+1]=+r_{q+1},\ \mathrm{if}\ r=1\\ \mathbf{p}_i[(q\ \mathrm{mod}\ 2^{\lambda-1})+1]=-r_{q+1},\ \mathrm{if}\ r=2 \end{array}\ ] ] where , @xmath183{q}}{2}$ ] are positive real numbers such that @xmath184{q}}{2}-1 $ ] and @xmath185{q}}{2}}r_i^2=\frac{\sqrt[4]{q}}{2}$ ] .
let @xmath186 and @xmath187 .
thus the rate of transmission of this code will be @xmath188 bits per channel use .
the design is given in .
the corresponding four dimensional signal set is described below .
@xmath189^t,~ \mathbf{p}_2=\left[\begin{array}{cccc}-r_1 & 0 & 0 & 0\end{array}\right]^t,\\ \mathbf{p}_3=\left[\begin{array}{cccc}0 & r_2 & 0 & 0\end{array}\right]^t,~ \mathbf{p}_4=\left[\begin{array}{cccc}0 & -r_2 & 0 & 0\end{array}\right]^t,\\ \mathbf{p}_5=\left[\begin{array}{cccc}0 & 0 & r_3 & 0\end{array}\right]^t,~ \mathbf{p}_6=\left[\begin{array}{cccc}0 & 0 & -r_3 & 0\end{array}\right]^t,\\ \mathbf{p}_7=\left[\begin{array}{cccc}0 & 0 & 0 & r_4\end{array}\right]^t,~ \mathbf{p}_8=\left[\begin{array}{cccc}0 & 0 & 0 & -r_4\end{array}\right]^t,\\ \mathbf{p}_9=\left[\begin{array}{cccc}r_5 & 0 & 0 & 0\end{array}\right]^t,~ \mathbf{p}_{10}=\left[\begin{array}{cccc}-r_5 & 0 & 0 & 0\end{array}\right]^t,\\ \mathbf{p}_{11}=\left[\begin{array}{cccc}0 & r_6 & 0 & 0\end{array}\right]^t,~ \mathbf{p}_{12}=\left[\begin{array}{cccc}0 & -r_6 & 0 & 0\end{array}\right]^t,\\ \mathbf{p}_{13}=\left[\begin{array}{cccc}0 & 0 & r_7 & 0\end{array}\right]^t,~ \mathbf{p}_{14}=\left[\begin{array}{cccc}0 & 0 & -r_7 & 0\end{array}\right]^t,\\ \mathbf{p}_{15}=\left[\begin{array}{cccc}0 & 0 & 0 & r_8\end{array}\right]^t,~ \mathbf{p}_{16}=\left[\begin{array}{cccc}0 & 0 & 0 & -r_8\end{array}\right]^t \end{array}\ ] ]    where , @xmath190 , @xmath191 , @xmath192 , @xmath193 , @xmath194 , @xmath195 , @xmath196 , @xmath197 .
though the dimension of the signal set is @xmath7 , due to the special structure of the signal set , we can study the two dimensional projections of the signal points which is graphically shown in fig .
[ fig_signal4dim ] at the top of the next page .
the proposed multidimensional signal sets in construction [ cons_signalset ] are far from generality and obtaining a general solution targeting to maximize the coding gain appears to be difficult .
[ thm_signalset ] construction [ cons_signalset ] provides fully diverse signals sets for the designs given in construction [ cons_design ] .
the constructed designs for @xmath198 complex variables has the structure as shown in construction [ cons_doubling ] .
the partitioning of the real variables into four groups is as follows .    1 .
first group : @xmath199 2 .
second group : @xmath200 3 .
third group : @xmath201 4 .
fourth group : @xmath202    the signal set is a cartesian product of four smaller dimensional signal sets @xmath203 . since and because the designs @xmath146 and @xmath147 are identical , it is sufficient to show that the design @xmath115 gives rise to fully diverse scaled unitary codewords for the chosen signal set .
let the given signal set be denoted by @xmath204 . then the following properties are satisfied @xmath205 . @xmath206x[j]&=&0,~
1\leq i\neq j\leq l\\ \label{eqn_fulldiv1 } 1\leq|\left\{i|~(x - y)[i]\neq 0 \right\}|&\leq&2\\ \label{eqn_fulldiv2 } ( x - y)[i]&\neq&\pm(x - y)[j],~ i\neq j\end{aligned}\ ] ] the design @xmath115 can be obtained by repeatedly applying abba construction on the linear design @xmath126 $ ] .
the off - diagonal entries of the matrix @xmath207 are sum of terms like @xmath208 .
but we have , which is equal to zero for the signal set @xmath209 by virtue of . thus we are guaranteed of scaled unitary codewords .    due to the abba structure of the linear design @xmath115
, the determinant can be obtained in an iterative manner .
let @xmath210 $ ] .
then we have    @xmath211\right|\\ & = & \left(|\mathbf{a_1}+\mathbf{a_2}|\right)\left(|\mathbf{a_1}-\mathbf{a_2}|\right ) .
\end{array}\ ] ]    the above equation suggests that @xmath212 can be computed recursively . for @xmath213 , . using this
, we can easily show that @xmath212 is a product of @xmath116 terms .
a typical term in the product looks like @xmath214 .
thus to ensure @xmath215 , each term in the product should not equal zero .
thus we need .
let us look at the real part of this expression .
we get @xmath216 .
for the given signal set @xmath209 , from and , we have .
thus the theorem is proved .
[ thm_relay ] for the linear designs given by construction [ cons_design ] , there exist @xmath1 relay matrices satisfying condition for arbitrary signal sets
. moreover , if the initial vector , then the initial matrix @xmath110 given in becomes unitary .
see appendix [ appendix2 ] .
[ eg_4relay ] let @xmath217 .
then the ddstc @xmath62 is obtained using the design shown in and the signal set given in construction [ cons_signalset ] .
the signal set is a cartesian product of four @xmath6-dimensional signal sets .
the relay matrices are given as follows : @xmath218,\\ \mathbf{a_3}=\left[\begin{array}{ccrr}0 & 0 & -1 & 0\\ 0 & 0 & 0 & -1\\ 1 & 0 & 0 & 0\\ 0 & 1 & 0 & 0\end{array}\right],~\mathbf{a_4}=\left[\begin{array}{ccrr}0 & 0 & 0 & -1\\ 0 & 0 & -1 & 0\\ 0 & 1 & 0 & 0\\ 1 & 0 & 0 & 0\end{array}\right ] .
\end{array}\ ] ] the initial vector @xmath219^t$ ] and the initial matrix @xmath220 .
this ddstc is single complex symbol decodable ( or 2 real symbol decodable ) .
in this section , we compare the error performance of the proposed ddstc in example [ eg_4relay ] with that of the circulant codes in @xcite and the cyclic codes in @xcite for transmission rates of @xmath107 bit per channel use ( bpcu ) and @xmath221 bpcu . the signal set chosen for the proposed code is as given by construction @xmath222 with parameters @xmath223 , @xmath224 for a transmission rate of @xmath107 bpcu and for the case of rate=1.5 bpcu , the signal set parameters are @xmath225 , @xmath226 , @xmath227 and @xmath228 .
the circulant code of @xcite is given by @xmath229 .
the relay matrices @xmath230 are taken to be powers of the @xmath231 circulant matrix .
the signal set @xmath232 is chosen to be @xmath233 radians rotated version of @xmath234-psk for a transmission rate of @xmath107 bpcu and @xmath235-psk for a rate of @xmath221 bpcu .
the rotation angles are chosen as follows to guarantee full diversity : @xmath236 , @xmath237 , @xmath238 , @xmath239 .
another reference for comparison is the cyclic codes proposed in @xcite whose codewords for rate=@xmath107 bpcu and rate=@xmath221 bpcu are given as follows : + @xmath240^i , i=0,\dots,255\right\}$ ] and @xmath241^i , i=0,\dots , 4095\right\}$ ] where , and @xmath242 . for simulations ,
we have assumed a block fading channel which is quasi - static for @xmath243 channel uses and varies independently from one block to another .
shows the error performance curves of the proposed codes in comparison with those of @xcite .
it can be observed from fig .
[ fig_simulation ] that for a transmission rate of @xmath107 bpcu , the proposed code outperforms the cyclic code by about @xmath244 db and the circulant code by more than @xmath245 db . similarly for a transmission rate of @xmath221 bpcu , it can be observed that the proposed code outperforms the cyclic code by about @xmath244 db and the circulant code by more than @xmath246 db . it is worthwhile to note that for rate=@xmath107 bpcu , the decoding search space for the proposed code is only @xmath7 whereas it is @xmath247 for the other two codes .
similarly for rate=@xmath221 bpcu , the decoding search space for the proposed code is only @xmath106 whereas it is @xmath248 for the other two codes .
because of the exponential growth of decoding search space with increasing transmission rate for the codes in @xcite , it takes a prohibitively large time to complete a error performance simulation for higher rates .
we have thus constructed a class of four group decodable ddstcs for any power of two number of relays using algebraic techniques .
it is important to note that relaxing the unitary matrix codebook to scaled unitary matrix codebook has aided us in obtaining decoding complexity benefits .
a limitation of this algebraic method is that it is available only for power of two number of relays . extending this method to obtain @xmath8-group
decodable ddstcs for @xmath249 appears to be difficult for two reasons :    1 .
the rate of the known linear designs in the literature for @xmath249 fall below @xmath107 and obtaining @xmath1 unitary relay matrices satisfying condition is difficult and it is not clear to the authors whether such matrices exist at all . for example , constructing the relay matrices satisfying condition for orthogonal designs is difficult .
2 .   even if the relay matrices are obtained , an initial vector satisfying condition may not exist for the chosen relay matrices .    for the specific case of @xmath250 ,
things were possible because the linear designs correspond to left regular matrix representation of an algebra over @xmath251 and hence using algebraic techniques , it was possible to construct the relay matrices and the initial vector satisfying all the required conditions .    optimizing the signal sets for coding gain might be possible for small number of relays and for small number of points .
however it becomes difficult for larger number of relays and/or larger number of points because it involves many parameters and a general closed form solution is difficult .
hence optimizing the signal sets for coding gain is an important direction for further work .
let @xmath252 . an extended clifford algebra denoted by @xmath253 is the associative algebra over @xmath254   generated by @xmath255 objects @xmath256 and @xmath257 which satisfy the following relations :    * @xmath258 ;  @xmath259 * @xmath260 ;  @xmath261 * @xmath262 .
the classical clifford algebra , denoted by @xmath263 , is obtained when only the first two relations are satisfied and there are no @xmath264 . thus @xmath263 is a sub - algebra of @xmath253 .
let @xmath265 be the natural @xmath254   basis for this sub - algebra .
@xmath266 then a natural @xmath254   basis for @xmath253 is : + @xmath267 we need a unitary matrix representation for the symbols @xmath107 , @xmath268 , @xmath269 , @xmath270 , , in the algebra @xmath271 .
such matrices are naturally provided by the left regular representation .
we first view @xmath271 as a vector space over @xmath251   by thinking of @xmath268 as the complex number @xmath272 . a natural @xmath251   basis for @xmath271
is given by : @xmath273 thus the dimension of @xmath271 seen as a vector space over @xmath251   is @xmath274 .
we have a natural embedding of @xmath271 into @xmath275 ( the set of all @xmath251-linear maps from @xmath271 to itself ) given by left multiplication as shown below . @xmath276
since @xmath277 is @xmath251-linear , we can get a matrix representation of @xmath277 with respect to the natural @xmath251   basis @xmath278 . left regular representation yields unitary matrix representations for the required symbols in the algebra .
the resulting linear designs are precisely those given by construction [ cons_design ] .
this has been explicitly shown in more detail in @xcite .
let us prove the theorem assuming that all the complex variables of the linear design take values from the entire complex field .
we use the fact that the linear design for @xmath125 relays was obtained as a matrix representation of the extended clifford algebra @xmath279 .
thus @xmath280 .
we choose @xmath281 .
the @xmath44 relay matrices are explicitly given by the union of the elements of the sets @xmath282 and @xmath283 . by virtue of the property
that @xmath284 is a ring homomorphism , these matrices are guaranteed to commute with all the codewords because they are matrix representations of elements belonging to the center of the algebra @xmath279 . to obtain the remaining @xmath45 relay matrices , we need to find unitary matrices which satisfy @xmath285 where , @xmath12 is any codeword .
but the codeword @xmath12 is simply a matrix representation of some element belonging to the extended clifford algebra .
one method to get these relay matrices is to take them from within the extended clifford algebra itself . by doing so , we can translate the condition in into a condition on elements of the algebra which will then provide us a handle on the problem . towards that end ,
we first identify a map in the algebra which is the analogue of taking the conjugate of the matrix representation of an element .
note from appendix [ appendix1 ] that we used the fact that @xmath268 can be thought of as the complex number @xmath272 .
when we take the conjugate of a matrix , we simply replace @xmath25 by @xmath286 .
hence the analogue of this action in the algebra is to replace @xmath268 by @xmath287 .
thus , the analogous map @xmath288 in the algebra is defined as follows : @xmath289 where , the element @xmath290 is obtained from @xmath16 by simply replacing @xmath268 by @xmath287 in the expression of @xmath16 in terms of the natural @xmath291-basis of extended clifford algebra .
now the problem is to find @xmath45 distinct elements denoted by @xmath292 of the algebra @xmath279 which satisfy @xmath293 .
the elements of the union of the following two sets satisfy the above required condition .
@xmath294 this can be proved by using the anti - commuting property , i.e. , @xmath295 .
hence the matrix representation of these specific elements gives the unitary relay matrices @xmath296 .
suppose we plug in these relay matrices to form a linear design where , @xmath297^t$ ] , it turns out that we get exactly the same linear design which is used at the source . thus the initial vector choice of @xmath219^t$ ]
guarantees that the initial matrix @xmath110 is an identity matrix .    however , we would like to point out that there are also other elements of the algebra which satisfy .
for example , consider the union of the elements of the sets @xmath298 and @xmath299 .
the authors thank prof .
hamid jafarkhani , dr . yindi jing and dr .
frdrique oggier for providing us with preprints of their recent works @xcite .
the authors are grateful to the anonymous reviewers for providing constructive comments which helped in improving the presentation of this paper .
f. oggier , b. hassibi , " cyclic distributed space - time codes for wireless relay networks with no channel information ,  submitted for publication .
m. tao and r. s. cheng ,  differential space - time block codes , " proceedings of _ ieee globecom 2001 _ , vol .
2 , pp . 1098 - 1102 , san antonio , usa , nov 25 - 29 , 2001 .
o. tirkkonen , a. boariu and a. hottinen , `` minimal non - orthogonality rate 1 space - time block code for @xmath300 tx antennas , '' in _ ieee int .
symp . on spread - spectrum tech . &
_ , new jersey , sep 6 - 8 , 2000 , pp .
429 - 432 .
g. susinder rajan and b. sundar rajan , `` algebraic distributed space - time codes with low ml decoding complexity , '' proceedings of _ ieee intl .
inform . theory _
, nice , france , june 24 - 29 , 2007 , pp .
1516 - 1520 .
g. susinder rajan ( s2006 ) was born in chennai , india in 1983 .
he completed his b.e .
degree at the college of engineering guindy , anna university , india in 2005 .
he is currently a ph.d .
student in the department of electrical communication engineering , indian institute of science , bangalore , india .
his primary research interests include space - time coding for mimo channels and wireless relay channels with an emphasis on algebraic code construction techniques .
b. sundar rajan ( s84-m91-sm98 ) was born in tamil nadu , india .
he received the b.sc .
degree in mathematics from madras university , madras , india , the b.tech degree in electronics from madras institute of technology , madras , and the m.tech and ph.d .
degrees in electrical engineering from the indian institute of technology , kanpur , india , in 1979 , 1982 , 1984 , and 1989 respectively .
he was a faculty member with the department of electrical engineering at the indian institute of technology in delhi , india , from 1990 to 1997 .
since 1998 , he has been a professor in the department of electrical communication engineering at the indian institute of science , bangalore , india .
his primary research interests are in algebraic coding , coded modulation and space - time coding .
rajan is an editor of ieee transactions on wireless communications from 2007 and also a editorial board member of international journal of information and coding theory .
he is a fellow of indian national academy of engineering and is a member of the american mathematical society .