This paper proposes a new distributed spacetime block code (DSTBC) over frequencyselective fading channels for twohop amplify and forward relay networks, consisting of a source node (), two relay nodes ( and ), and a destination node (). The proposed DSTBC is designed to achieve maximal spatial diversity gain and decoupling detection of data blocks with a lowcomplexity receiver. To achieve these two goals, uses zerosequence padding, and relay nodes precode the received signals with a proper precoding matrix. The pairwise error probability (PEP) analysis is provided to investigate the achievable diversity gain of the proposed DSTBC for a general channel model in which one hop is modeled by Rayleigh fading and the other by Rician fading. This mixed RayleighRician channel model allows us to analyze two typical scenarios where are in the neighborhood of either or .
1. Introduction
The reliability of wireless communications over fading channels can be greatly improved by the use of diversity schemes. For multipleinput multipleoutput (MIMO) systems, transmitted diversity can be realized in the form of spacetime block codes (STBCs) [1, 2]. With lowcomplexity maximumlikelihood (ML) decodability and high achievable diversity gain, STBCs are widely used for wireless communications. Generally, the conventional STBCs were designed for the colocated antennas, and thus are easily deployed at the base station to improve the performance of the downlink transmission. Nevertheless, the realization of STBCs is impractical in the uplink transmission due to the constraints on size and hardware complexity in mobile handsets. Fortunately, mobile users can cooperate to form a virtual multipleantenna system, which is now known as cooperative diversity [3, 4].
The distributed spacetime block codes can be viewed as the distributed implementation of conventional STBCs for cooperative communications. Originally, the DSTBCs were proposed for flatfading channels [5–7]. The problem of DSTBC in frequencyselective fading channels was investigated in [8] with decodeandforward (DF) relaying, and in [9] with amplifyandforward (AF) relaying. However, these DSTBCs were devised for relay networks where there exists one active relay node and a direct communication link between the source and the final destination.
In this paper, we design a new DSTBC for twohop relay networks [6, 7] over frequencyselective fading channels with AF protocol, where there are two active relay nodes. The proposed DSTBC operates as follows: in the first time slot, the source () sends two blocks of information data to two relays ( and ). What is remarkable in our proposed DSTBC is that one of the two relays precodes its received signals that will be sent to the destination () in the next time slot. The precoding matrix is designed such that each relay conveys a distinct column of the block Alamouti scheme (see, e.g., [1, 10, 11]). Our main contributions in this paper are summarized in brief as follows.
(i)With our proposed DSTBC, the data rate of 1/2 is achieved, which is proved to be the maximum data rate for twohop relaying networks. As we can see later, the extension of [9] to two relays results in a rate of 1/3.
(ii)We propose the precoding matrix at the relays such that the decoupling detection of two data blocks in both time and frequency domains is possible at . The PEP analysis is carried out with ML detection in time domain, and numerical results are obtained with minimum mean square error (MMSE) receiver in frequency domain.
(iii)We study the achievable diversity gain of the proposed DSTBC for the general scenario where the relays are located near the source or the destination, that is, one of the two hops ( to or to ) is lineofsight (LOS) transmission, while the other is nonlineofsight (NLOS) transmission. Accordingly, the considered channel model is a mix of Rayleigh and Rician fading.
The theoretical results prove that our proposed scheme achieves the spatial diversity order of , where and are the channel memory lengths for the links from to and from to , respectively. The analysis also shows that the factor of Rician fading in the LOS component provides a coding gain to the PEP performance. It means that as the factor increases, a better performance is observed.
The rest of this paper is organized as follows In Section 2, we describe the system model of the proposed DSTBC and the proof of decoupling capability in time domain and frequency domain. Performance analysis is presented in Section 3. We present the numerical results in Section 4, and Section 5 concludes this paper.
Notation 1.
Bold lower and upper case letters represent vectors and matrices, respectively; , , and denote transpose, complex conjugate, and Hermitian transpose operations, respectively; and denote an identity matrix and an allzero matrix of size ; denotes the expectation; denotes the Euclidean norm of a vector; stands for a Fast Fourier Transform (FFT) matrix of size ; and represent the links from the source () to the th relay () and from to the destination (), respectively.
2. System Model and the Proposed DSTBC
We consider a fournode wireless relay network shown in Figure 1, where the source terminal cannot communicate directly with the intended destination. The data transmission from to is completed via twohop protocol [6, 7] with the assistance of two relays and . The frequencyselective channel from to is characterized by , where is the channel memory order. Two transmitted data blocks and of length , shown in Figure 1, are created by padding a zero sequence of length to two information data blocks , of length . To achieve the decoupling property of data detection, the length of the zero sequence must satisfy [9]. This condition makes circulant channel matrices from source to relays and relays to destination.
Figure 1. System model and data structures.
In the first time slot, the source serially transmits two data blocks to two relays. In the next time slot, one relay only amplifies and forwards its received signals, while the other precodes its received data blocks by a precoding matrix before transmitting to the destination as illustrated in Figure 1. The idea behind our design is that precoding in is designed to send the second column of the block Alamouti's scheme (see, e.g., [1, 11]) to . This enables the decoupling detection of two data blocks at the destination and achieves a rate of 1/2. To achieve the same goal in the considered scenario, the source with repetition code in [9], which is devised for onerelay system, must send two columns during two time slots. Thus, the rate of this scheme is reduced to 1/3. Recall that the maximum achievable data rate of an relay repetitioncoding network is [12]. Clearly, our design can achieve a higher data rate transmission.
We now proceed to prove that our proposed DSTBC can decouple the detection of two data blocks. Throughout this paper, the superscript denotes the relay index, while the subscript refers to data block index. The received signal at the relay is given by
where is the average energy of link; is the circulant channel matrix of link; is the white Gaussian noise vector at the th relay with each entry having zeromean and variance of per dimension. For any circulant matrix , its entry is written as .
At , the received data is conjugated, followed by the precoding operation which is denoted by precoding matrix
where the matrix is designed as
In (3), the matrix of size and the matrix of size have element given, respectively, by
We choose to ensure that, after precoding, at least last samples of and are all zeros to make the circulant channel matrix Before transmitting the signals to the destination, the relay normalizes its signals by to have unit average energy. The received signals at the destination are written by
where is the average energy of link; is the normalized received signal; is the circulant matrix, denoting the channel ; is white Gaussian noise vector at the destination with each entry having zeromean and variance of .
Using (1), we can rewrite (5) as
where and include the Gaussian noise of relays and destination. It is common to normalize the noise variance in (6) to be , which results in
The values of normalization factors in (7) and (8) are defined by
where and for .
By conjugating, and multiplying both sides of (8) with , and noting that for any circulant matrix , we can rewrite (8) as
For mathematical convenience, we group (7) and (10) in vectormatrix form as
where
Let us denote = + , where = for any circulant matrix . Then, = is a blockdiagonal matrix. By multiplying both sides of (11) with the unitary matrix , we can decouple the detection of two data blocks. That means two data blocks can be detected independently, rather than joint detection, without any loss of gain, which is based on the following model
where is the Gaussian noise vector resulting from the spacetime decoupling.
Since the general ML detection or standard equalization techniques based on (13) require complex receiver, we introduce the decoupling in the frequency domain that allows for lowcomplexity equalizers [13].
We notice that any circulant matrix can be diagonalized as where is the diagonal matrix whose diagonal elements are the DFT of the first column of . Taking the DFT of both sides of (11) results in
where , ,
Let us denote . Then, is a blockdiagonal matrix. By multiplying both sides of (14) with the unitary matrix , we can decouple the detection of each data block in frequency domain as
Since the matrix is diagonal, (16) can be decomposed into a set of scalar equations
In (17), we omit the dependency of subscript because the detection of two data blocks is based on the same model. Typical frequency domain equalizers in [13] can be applied to the outputs of decoupling process.
3. Performance Analysis
In this section, we derive the PEP expression of the proposed DSTBC model over frequencyselective fading channels based on the joint channel model (11), where the links and experience Rician fading and Rayleigh fading, respectively. For the case when and are Rayleigh fading and Rician fading, the PEP expression can be similarly obtained with some interchanged parameters.
Defining the decoded codeword vector as and the Euclidean distance between and as , the conditional PEP under fading channels is given by
where is the function. By applying Chernoff bound to function, this PEP is upper bounded by
The Euclidean distance in (19) is calculated by
Using the fact in [9], we can approximate (20) to
where and . We note that and ; . denotes the eigenvalue of codeword difference matrix, and is zeromean complex Gaussian vectors with unit variance. In (21), we mean that the distance can be approximated by one of four possible forms. Each component of the summations of the right hand side of (21) can be expressed by one of the two following factors
or
In (22), (23), and for the rest of the paper, the subscript stands for or for sake of generality because the links and can be treated in the same way. To derive the PEP, we differentiate three cases based on the relation of and because of the different characteristics of fading and .
Case 1 ().
We consider (22) and define where and . Applying Chernoff bound, the PEP corresponding to is upper bounded by , where denotes the momentgenerating function. If we consider (), will be corresponding to (). can be evaluated as [14]
where is the probability density function. Since the fading channels and are frequencyselective Rician and Rayleigh fading, respectively,
where is the Nakagami or Rician fading parameter, and represents the Gamma function defined by for any positive integer .
Substituting (26) and (25) into (24), we have
We can rewrite (27) as
Assuming high signaltonoise ratio (SNR), that is, and , (28) can be evaluated as
The integral in (29) when is given by [15]
Substituting (30) into (29), we obtain
Case 2 ().
We examine (23) and similarly define where and . Applying Chernoff bound, the PEP corresponding to is upper bounded by which is given by
where
The modified Bessel function of the first kind, in (33), is defined as
Substituting (33) and (34) into (32), we have
Under the assumption of high SNR (), (36) can be computed as
Let . After some mathematical manipulations [15], we obtain
where is hypergeometric regularized function and is defined as
Substituting (38) into (37), we obtain
Case 3 ().
Considering (22), can be calculated as
Let , (41) can be rewritten as
Let us define
By using some mathematical expansions [15], we can express as
where . Calculating the integral in (44), the value of is evaluated as [15]
where the incomplete Gamma function is defined by .
Substituting into (42), we have
Generalizing the three above cases with (31), (40) and (46), we can conclude that the diversity gain of our proposed DSTBC is by extracting the exponential terms of . Although the Rician fading parameter does not produce any diversity gain, it acts as a coding gain, and thus can improve the PEP.
We provide here one example of PEP calculation for the case and . Suppose that and , the normalization factor in (9) can then be approximated as . From the result in (31), the PEP is upper bounded by
We observe that the diversity order of this case is + . Additionally, PEP reduces to zero when the Rician parameter goes to infinity. That explains the positive effect of LOS component on the quality of received signals.
4. Numerical Results
As aforementioned, the PEP analysis is carried out based on (11). Moreover, the derivation of the PEP is accomplished under the assumption of an ML detection scheme at the receiver. However, such an ML scheme for the system model in our paper requires high complexity. Instead, we use MMSE frequencydomain linear equalizer to verify the analytical results. Since MMSE receiver is a suboptimal solution to the data detection problem, there is always a gap between the performance of ML and MMSE schemes. Note that the slopes of the performance curve of ML and MMSE receivers are similar at high SNR [16]. Consequently, the diversity gain can also be verified based on the performance of MMSE receiver.
In this section, we evaluate the BER performance of the proposed DSTBC via Monte Carlo simulation to justify the analysis of the achievable diversity gain. By observing the slope of these BER curves, we can confirm the validity of analysis since BER is proportional to PEP [16]. Each data block consists of 64 symbols including the zero sequence and informationcarrying data which is modulated by QPSK with Gray mapping. We assume that the receiver has perfect channel state information.
Figure 2 shows the BER performances of the proposed DSTBC for various combinations of channel lengths using frequency domain linear equalizer (FDLE), such as MMSE receiver. The fading and are assumed to be Rayleigh fading. We assume that the value of is fixed at 25 dB, dB, and plot the BER curves as a function of . Let us show the comparison of three cases.
Figure 2. BER performance of DSTBC with FDLE for various combinations of channel memory orders.
Case 1.
(black curve and square marker); the diversity order is calculated by .
Case 2.
(red curve and round marker); the diversity order is .
Case 3.
(green curve and uppertriangular marker); the diversity order is .
The simulation results indicate that the slopes of the curves for Case 2 and Case 3 are steeper than that for Case 1 since Case 1 achieves smaller diversity gain. We can also see that the curves for Case 2 and Case 3 have the same slope, which can also be shown using the diversity order expressions we derived. These facts confirm our conclusion of the achievable diversity gain in the analysis. Another observation is for the same value of and (e.g., red curve and green curve), the better performance is achieved as the number of paths from the relays to destination increases. However the slopes of these BER curves are identical at high SNR since they have the same diversity gain of .
The performance of DSTBC over wireless fading channels where the fading and are characterized by Rician and Rayleigh distributions, respectively, is drawn in Figure 3 with different values of factor and the assumption of channel memory order of . We remark that the better performance is achieved as the factor of Rician fading increases. However, the slopes of the performance curve are almost the same at high SNR, that is, they achieve the same diversity gain. This is because the term only produces the coding gain to the PEP as we can see in the analysis. This confirms our conclusion about the achievable diversity gain from the PEP analysis mentioned in the previous section.
Figure 3. BER performance of DSTBC with different values of factor over Rician fading channels.
The performances of two cases and for various combinations of channel memory orders are shown in Figure 4 to study the effects of LOS component and channel memory orders on the BER. For each value of , the BER curves corresponding to different channel memory orders are plotted to show the variations of slope due to different achievable diversity gains. We observe that the BER curves are shifted down if we increase . That reflects the effect of LOS component on the coding gain through all range of SNR from low to high SNR. Meanwhile, the channel memory orders only have noticeable effects on the BER curves at high SNR caused by the change of diversity gain.
Figure 4. BER performance of DSTBC with different values of channel memory orders and factor over Rician fading channels.
5. Conclusion
A DSTBC scheme for twohop cooperative systems over frequencyselective fading channels has been proposed. The proposed DSTBC can achieve halfdata rate transmission, spatial diversity order of . Furthermore, the data detection of the proposed DSTBC can be decoupled in frequency domain for a lowcomplexity receiver. The main idea in our design is that each relay conveys a distinct column of the conventional STBC developed for colocated antennas. We prove that the data rate of the proposed DSTBC is higher than that of repetition code which has the maximum data rate of 1/3 for a tworelay system. The design of DSTBC in our paper is inspired by the idea of linear dispersion codes where each column of the code matrix is a linear summation of transmit data and their conjugations [17]. In this paper, we analyze the PEP with a mixed Rayleigh and Rician frequencyselective fading channel model. The analysis shows that the LOS component effectively improves the error rate performance at the destination. Simulation results over different channel parameters and different factors were provided to verify the theoretical analysis.
Acknowledgments
This paper was partly supported by the IT R&D program of MKE/KEIT (KI001814, Game Theoretic Approach for Crosslayer Design in Wireless Communications) and MKE (The Ministry of Knowledge Economy), Korea, under the ITRC (Information Technology Research Center) support program supervised by the NIPA (National IT Industry Promotion Agency) (NIPA2010(C109010110001)).
References

SM Alamouti, A simple transmit diversity technique for wireless communications. IEEE Journal on Selected Areas in Communications 16(8), 1451–1458 (1998). Publisher Full Text

V Tarokh, H Jafarkhani, AR Calderbank, Spacetime block codes from orthogonal designs. IEEE Transactions on Information Theory 45(5), 1456–1467 (1999). Publisher Full Text

A Sendonaris, E Erkip, B Aazhang, User cooperation diversity—part I. System description. IEEE Transactions on Communications 51(11), 1927–1938 (2003). Publisher Full Text

JN Laneman, DNC Tse, GW Wornell, Cooperative diversity in wireless networks: efficient protocols and outage behavior. IEEE Transactions on Information Theory 50(12), 3062–3080 (2004). Publisher Full Text

JN Laneman, GW Wornell, Distributed spacetimecoded protocols for exploiting cooperative diversity in wireless networks. IEEE Transactions on Information Theory 49(10), 2415–2425 (2003). Publisher Full Text

Y Jing, B Hassibi, Distributed spacetime coding in wireless relay networks. IEEE Transactions on Wireless Communications 5(12), 3524–3536 (2006)

Z Yi, IM Kim, Singlesymbol ML decodable distributed STBCs for cooperative networks. IEEE Transactions on Information Theory 53(8), 2977–2985 (2007)

G Scutari, S Barbarossa, Distributed spacetime coding for regenerative relay networks. IEEE Transactions on Wireless Communications 4(5), 2387–2399 (2005)

H Mheidat, M Uysal, N AlDhahir, Equalization techniques for distributed spacetime block codes with amplifyandforward relaying. IEEE Transactions on Signal Processing 55(5), 1839–1852 (2007)

S Zhou, GB Giannakis, Spacetime coding with maximum diversity gains over frequencyselective fading channels. IEEE Signal Processing Letters 8(10), 269–272 (2001). Publisher Full Text

N AlDhahir, Singlecarrier frequencydomain equalization for spacetime blockcoded transmissions over frequencyselective fading channels. IEEE Communications Letters 5(7), 304–306 (2001). Publisher Full Text

T Wang, GB Giannakis, Complex field network coding for multiuser cooperative communications. IEEE Journal on Selected Areas in Communications 26(3), 561–571 (2008)

D Falconer, SL Ariyavisitakul, A BenyaminSeeyar, B Eidson, Frequency domain equalization for singlecarrier broadband wireless systems. IEEE Communications Magazine 40(4), 58–66 (2002). Publisher Full Text

A Stuart, JK Ord, Kendall's Advanced Theory of Statistics, h (C. Griffin, London, UK, 1987) PubMed Abstract  Publisher Full Text  PubMed Central Full Text

IS Gradshteyn, IM Ryzhik, Table of Integrals, Series, and Products, h (Academic Press, San Diego, Calif, USA, 2007)

JG Proakis, Digital Communications, h (McGraw Hill, New York, NY, USA, 2001)

B Hassibi, BM Hochwald, Highrate codes that are linear in space and time. IEEE Transactions on Information Theory 48(7), 1804–1824 (2002). Publisher Full Text