TwoWay RelayingBased TwoHop TwoUser MultipleInput MultipleOutput System
 Author: Cho Young Seek, Choi Seyeong
 Organization: Cho Young Seek; Choi Seyeong
 Publish: Journal of information and communication convergence engineering Volume 12, Issue2, p67~74, 30 June 2014

ABSTRACT
In multihop communication systems, twoway relaying is one of the solutions to mitigate the spectral efficiency loss caused by a halfduplex transmission. In this paper, a simple twoway relaying scheme is proposed for twohop twouser multiple input multiple output (MIMO) systems. In the proposed system, a base station and a relay station (RS), both equipped with two antennas, form a pointtopoint MIMO channel, while the RS and two singleantenna mobile users form a pointtomultipoint multiuser (MU)MIMO channel. Numerical examples show that the proposed system achieves a significant sum rate gain as compared to a oneway relaying system as the distance between a relay and the two users decreases. We also show that although we can expand the proposed scheme to more than two users, its performance gain as compared to that of oneway relaying decreases with an increase in the number of users.

KEYWORD
Halfduplex transmission , MUMIMO , Multihop communication , Twoway relaying

I. INTRODUCTION
Nextgeneration wireless communication systems will demand enhanced link capacity and larger coverage. One of the solutions to meet these demands is to use spacedivision multiple access in a multiuser multiple input multiple output (MUMIMO) system, which can achieve a multiplexing gain even when each user has a small number of antennas [13]. However, the problem arises when the users are far away from the base station (BS) and there are many obstructions between the BS and the users, such as in an urban cell edge area, where the users experience significant path loss and shadowing. The conventional solution to overcome this coverage shrinking problem is to increase the density of the BS or the repeater. Unfortunately, additional BS deployment is impractical due to its considerably high cost. Moreover, the repeater has the problem of amplifying interference and finds it difficult to perform complicated signal processing.
One of the practical schemes to enhance the coverage area in a cellular network is fixed relay station (RS)based multihop communication [4, 5]. In this system, the BS and the RS form a pointtopoint MIMO channel, while the RS and the users form a pointtomultipoint MUMIMO channel. Deploying the RS with multiple antennas in a cell edge area, where the nearby users experience significant shadowing, enables the system to achieve a high data transmission rate at the BS and the users [6]. In addition, since the fixed RS can perform more complex operations than a mobile relay, the RS can take over a portion of the operations in the BS [7]. Thus, it is possible to reduce the complexity of the BS by using the RS as a subBS. However, the problem of a relay system is that all the terminals in this system operate in the halfduplex mode; this leads to a loss in the spectral efficiency [8]. Conventional techniques to solve this problem include opportunistic relaying, twopath relaying, and twoway/bidirectional relaying [817]. Among them, in this study, we consider a twoway relaying system where a relay simultaneously transmits the uplink data symbols to the BS and the downlink data symbols to the users [1015].
The coded twoway relaying schemes have been shown in [10, 11] by using bitwise and symbolwise XOR operations at the RS, respectively. A combination of twoway relaying and cooperative relaying has been proposed in [12]. In [13, 14], an informationtheoretic approach for various twohop twoway relaying schemes has been introduced by showing that the sum rate can be increased significantly as compared to the oneway relaying scheme. In [15], the twohop twoway relayingbased MUMIMO system has been proposed, where the decodeandforward operation is performed at the RS, while the successive interference cancelation is performed by the BS and the users. However, the author in [15] did not describe the transmitting signal construction at the RS in detail. In addition, the feedback load and the computational complexity of the proposed scheme in [15] are very high.
In this paper, we propose a lowfeedbackload and lowcomplexity twoway relaying scheme for a twohop twouser MIMO system. In the proposed scheme, one needs three phases to transmit the uplink and downlink data symbols and the operation of the first two phases is the same as in [15]. In the third time slot, the RS decodes the received signal and constructs the transmitting signal by using singular value decomposition (SVD) and QR decomposition and simultaneously transmits this signal to the BS and the users. We also propose an optimal power allocation algorithm at the RS for the uplink and downlink data symbols.
The rest of this paper is organized as follows: Section II introduces the system model and notations, and Section III describes the process of the proposed scheme, computes its sum rate, and proposes an optimal power allocation. Section IV shows the simulation results of the proposed scheme in comparison with the conventional schemes. Finally, Section V presents the conclusion.
II. SYSTEM MODEL
The considered twohop twouser MIMO system is as illustrated in Fig. 1. There are only two users in this system, and the direct path between the BS and both the users is blocked due to the path loss and shadowing. Also, the time division duplex (TDD) mode is assumed, where the channels are assumed to be constant during one TDD frame but change independently between different frames. Both the BS and the RS are assumed to have two antennas, and each user has a single antenna. The matrix
γ_{b} H _{b} and the vectorγ_{ui} h _{ui} (i = 1, 2) model the downlink channels from the BS to the RS and the uplink channel from thei th user to the RS, respectively. All the elements ofH _{b} andH _{u} ( = [h _{u1},h _{u2}]) are independent, zero mean complex Gaussian random variables with unit variance. The distance between the RS and thei th user and the distance between the BS and the RS are denoted asd_{ui} andd_{b} km, respectively. A simple path loss model is considered, where the power of the transmitted signal decays as a function of ( j ∈ {b ,u _{1},u _{2}}) ; i.e., the signal magnitude decays likeγ_{j} = . In addition, the effect of path loss and shadowing at the direct path between the users is also assumed to be negligible. The uplink channel matrix from the RS to the BS and the downlink channel matrix from the RS to the users are denoted as and , respectively, where (⋅)^{H} denotes the Hermitian operator. The downlink and uplink transmitting data symbols between the BS and thei th user are denoted ass_{i} andt_{i} , respectively, where s_{i} =t_{i} =1. The total power for one TDD frame of the BS, the RS, and the two users is constrained asP_{b} ,P_{r} , andP_{u} , respectively.III. PROPOSED SCHEME
> A. Previous Work
In the proposed scheme in [15], a TDD frame is equally divided into three time slots. In the first and the second time slots, the BS transmits two data symbols to the RS by using the space time block code [18]. At the same time, user 1 and user 2 transmit the uplink data symbols to the RS in the first and the second time slots, respectively. While one of the two users, say user 1, transmits the uplink data symbol, the other user listens to the transmission of user 1. Then, the received signal vector at the RS in the
i th time slot, , can be written asand
The noise terms, (
i = 1, 2), are the additive white Gaussian noise (AWGN) vector with zero mean and covariance matrix , whereI_{n} denotes ann ×n identity matrix.After the transmissions of the first and the second time slots, the RS decodes the downlink and uplink data symbols by using minimum mean square estimation and successive interference cancelation (MMSESIC) [19] depending on the two decoding ordering methods. For example, at first, the RS decodes the uplink data symbols by treating the downlink data symbols as interference and subtracts the uplink data symbols from the received signal vector and then, decodes the downlink data symbols without interference. In this case, the achievable sum rates of the downlink and uplink transmission, and , are given as [20]
and
where tr {·} denotes the trace operator.
Similarly, if the RS first decodes the downlink data symbols and then, decodes the uplink data symbols, and are given as
and
In the third time slot, the RS executes the waterfilling algorithm [20] to transform uplink data symbols and the dirty paper code (DPC) [21] into the downlink data symbols, which are denoted as
b _{s} andb _{t} , respectively. Then, the RS simultaneously transmitsb =b _{s} +b _{t} to the BS and the users.Note that the BS and the users need to know the perfect channel state information (CSI) of the other side to cancel the interference, which leads to a high feedback load. In addition, although executing the DPC at the RS achieves the optimal sum rate for the MIMO broadcast channel, it may be difficult to implement this in the practical system and its complexity may be very high.
> B. Proposed TwoWay Relaying Strategy
For the first two time slots, the operation of the proposed scheme is the same as that in [15]. Let
ŝ_{i} and (i = 1, 2) denote the decoded downlink and uplink data symbols at the RS, respectively. Then, the RS constructs the transmitting signal vector as follows:The application of the SVD to the channel matrices
H _{b} and QR decomposition to the channel matrixH _{u} yields [22]where
U ,V , andQ are the unitary matrices, Σ has nonzero diagonal elements placed in the descending order (λ_{1} > λ_{2} > 0), and is the lower triangular matrix whose diagonal terms,g _{1} andg _{2}, are positive.Let and where
p_{si} andp_{ti} (i = 1, 2) denote the transmit powers for the downlink and uplink data symbols at the RS,p_{t1} +p_{t2} +p_{s1} +p_{s2} ≤P_{r} , and (⋅)^{T} represents the transpose operator. The RS constructs the transmitting signal vector,x _{r} , aswhere
D (g _{1} /α ,g _{2}) denotes the diagonal matrix whose diagonal terms areg _{1} /α andg _{2}, and .However, if
α is larger than a certain value, it would be better in terms of the sum rate to load all the power of the downlink data symbols to the user of the largest channel gain, which will be described in more detail at the end of the following subsection.At the BS, the received signal is multiplied by the unitary matrix
V ^{H} , and then, the output of the multiplier,n _{b} , can be expressed aswhere
n _{b} denotes the AWGN vector with zero mean and covariance matrix . Note that, sinceV ^{H} is a unitary matrix, the covariance matrix ofV ^{H} n _{b} is the same as that ofn _{b} . On the other hand, the received signal vector at the users,r _{u} , can be expressed aswhere
n _{u} denotes the AWGN vector with zero mean and covariance matrixD (, ).Before transmission, the RS feeds back the information of the correlation matrix
U ^{H} ,g _{3} /g _{2}, and power loading parameters to the BS and the users.Using this feedback information, the BS and the users execute the SIC to obtain
r _{b} andr _{u} , respectively. Then, the output signal vectors of the BS and the users after the execution of the SIC,y _{b} andy _{u} can be written asand
Note that the feedback load of the proposed protocol is considerably lower than that obtained in the previous work. Perfect channel knowledge of the other path
H _{u} andH _{b} is required to execute the SIC at the BS and the users, respectively. On the other hand, in the proposed protocol, the BS and the users require the same information ofU ^{H} Q andg _{3} /g _{2} (since is the multiplication of the unitary matrix, its feedback load is considerably less than that of or .).The results of the electromagnetic simulations for the lumped element and the distributed design show that the attenuation characteristics of the lumped design are superior to those of the distributed one.
> C. Sum Rate of the Proposed Scheme
After transmission of the third time slot, the sum rate of the uplink data symbols at the BS can be computed from (11), as follows:
In addition, the sum rate of the downlink data symbols at the both users can be computed from (12), as follows:
However, since additional power is used for presubtraction as in (8), there may be a situation in which the rate of the downlink data symbol is higher than the sum rate of the downlink data symbols of both users when the RS only transmits the downlink data symbol of the user with the largest channel gain . This situation occurs when the correlation of the two uplink channel vectors, i.e., is high, and thus,
α ^{2} >ρ_{th} (ifα ^{2} >ρ_{th} , then whereρ_{th} is given asTherefore, we need check whether
α ^{2} >ρ_{th} . If yes, then the RS transmits the downlink data symbol of the userm rather than the downlink data symbols of both the users. In this case, the transmitting signal vector at the RS can be constructed as (note that, in this case, the RS feeds back the information of instead of .)Then, the corresponding sum rate of the downlink data symbol can be calculated as
Consequently, the total system sum rate can be achieved in a TDD frame as
> D. Optimal Power Allocation at RS
In this paper, only the power allocation (PA) at the RS is considered since the BS and the users need to know the perfect CSI of two channel matrices,
H _{u} andH _{b} , to perform the PA, which is not available in the proposed system model.Based on the sum rate of the downlink and uplink data symbols and , which is achieved by using a certain decoding order, the optimal PA for
α ^{2} <ρ_{th} can be computed by solving the following optimization problem:where , and This optimization problem can be solved by the following algorithm:
• Step 1: Using the waterfilling algorithm, solve the sum rate optimization problem in (19) by considering only the power constraints as
where
µ denotes the Lagrange multiplier, which is chosen to satisfy the constraint and (x )^{+} max(x , 0). If the solution in (20) satisfies the other constraints in (19), the algorithm will be completed. If the solution does not satisfy both the downlink and the uplink sum rate constraints in (19), go to Step 2; else, go to Step 3.• Step 2: Since the power of the RS is enough to meet the sum rate of the downlink and the uplink data symbols at the RS, and , respectively, find the optimal PA by solving the following problem:
where
k ∈ {s ,t }. The solution of (21) can be found asThen, the algorithm is completed.
• Step 3: Let
k =s , if the solution in (20) does not satisfy the downlink sum rate constraint in (19); other,k =t . Find the optimal PA forp_{k} _{,1} andp_{k} _{,2}, by solving the following problem:The solution can be found as
Then, solve the following problem to calculate and , where
l ∈ {s ,t } andl ≠k ,The solution of (25) can be found as
Then, the algorithm is completed.
On the other hand, the optimal PA for
α ^{2} >ρ_{th} can be computed by solving the following optimization problem:where . This optimization problem can be solved by the similar processing for the proposed algorithm.
It should also be noted that the RS determines the decoding order to maximize the total sum rate if the RS can estimate the channels perfectly. For this reason, the RS needs to execute the power loading algorithm twice since there are two decoding ordering methods.
IV. NUMERICAL EXAMPLES
In this section, various MonteCarlo simulation results are illustrated by averaging more than 10,000 channel realizations. It is assumed that the noise variances are the distance between the BS and the users, d, is kept constant as 2 km; the distance between the RS and each user is
d_{u} ; andd =d_{u} +d_{b} . It is also assumed that the powers of the BS, the RS, and the users areP_{b} =P_{r} =P_{u} = 10 W.Fig. 2 shows the ergodic sum rate of the proposed scheme calculated on the basis of the optimal PA for a low and a high signaltonoise ratio (SNR). To show the optimality of the proposed PA, the suboptimal PA is illustrated, which achieves a nearoptimal sum rate at a high SNR. In the suboptimal PA, at first, the lower bound of the cost function
R_{LB} is found, which approximates + in the highSNR region, as given in [23].Then, the following optimization problem is solved by geometrical programming.
Fig. 2 shows that the proposed PA achieves the optimal sum rate of the system since suboptimal PA achieves the optimal sum rate as SNR reaches infinity. Note that the gap between the two PA schemes at a low SNR is caused by the lower bound optimization of suboptimal PA. Besides, Fig. 2 shows that fixing the decoding order to decode the uplink data symbol at the first time and the downlink data symbol at the next time can provide the nearoptimal sum rate as the distance between the RS and the user decreases.
In Figs. 3 and 4, the proposed scheme is compared to the previous work presented in [15] and the optimal oneway relaying scheme in order to demonstrate the performance gain of the proposed scheme. In the optimal oneway relaying scheme, a TDD frame is equally divided into four time slots. In the first and the second time slots, the BS transmits the downlink data symbols to the RS, and the users transmit the uplink data symbols to the RS. The RS decodes the downlink and uplink data symbols by using MMSESIC. In the third time slot, the RS transmits the decoded uplink data symbol to the BS by using the waterfilling algorithm. Finally, the RS transmits the decoded downlink data symbol to the users by using DPC in the fourth time slot. It is assumed that the optimal PA is performed at the RS for all the schemes.
Fig. 3 illustrates the sum rates versus of the three different schemes for 1 /
σ ^{2} = 10 dB. It shows that the proposed protocol achieves a significant gain over the one way relaying scheme and the gain increases as the distance between the RS and both the users becomes close tod_{u} = 0.7 km due to the decreased prelog factor. However, the performance of the twoway relaying schemes decreases asd_{u} becomes smaller or larger than 0.7 km. This is attributed to the fact that although the downlink sum rate increases asd_{u} decreases, the uplink sum rate decreases. Therefore, the maximum sum rate is achieved for an asymmetric channel situation [13]. It can also be seen in Fig. 3 that the performance gap between the previous work and the proposed scheme, which needs additional power to presubtract the other user interference as compared to the previous work, becomes negligible asd_{u} decreases.This is because as
d_{u} decreases, the allocated power to satisfy the downlink sum rate at the RS decreases and the portion of the uplink sum rate in the total system sum rate decreases. Fig. 4 illustrates the sum rate versus 1 /σ ^{2} of the three different schemes ford_{u} = 0.7 km. It can be observed that the proposed protocol still achieves the sum rate gain over the oneway relaying scheme irrespective of the SNR.V. DISCUSSION AND CONCLUSIONS
Twoway communication at the RS is one of the solutions to overcome the disadvantage of the time resource consumption due to the halfduplex constraint. A problem of this scheme is the degradation of the system sum rate mainly due to the interference. In the proposed scheme, the RS constructs the transmitting signal vector by employing SVD and QR decomposition, where the feedback load at the BS and the users is considerably lower than that in the previous work described in [15]. Simulation results show that the proposed protocol can achieve a significant sum rate gain as compared to the oneway relaying scheme. In addition, the DPC used in [15] has some difficulties for practical implementation. In this aspect, the proposed scheme can be easily implemented in various circumstances. Although we can expand the proposed scheme to more than two users, its performance gain as compared to oneway relaying decreases with an increase in the number of users. This is because each user needs one time slot for transmitting the uplink data symbol to the RS and the power consumption of the interference presubtraction for the downlink transmitting signal at the RS increases with an increase in the number of user.

[Fig. 1.] System model of twohop twouser multipleinput multipleoutput system. BS: base station, RS: relay station.

[Fig. 2.] Performance comparison between the proposed power allocation (PA) and the suboptimal PA for 1 / σ2 = 2 dB and 1 / σ2 = 10 dB.

[Fig. 3.] The sum rates versus the distance between the relay station and the users of three different schemes for 1 / σ2 = 10 dB. PA: power allocation, SNR: signaltonoise ratio.

[Fig. 4.] The sum rates versus the inverse of noise variance, 1 / σ2, of three different schemes for du = 0.7 km.