Decoupled beamforming techniques for distributed MIMO two-way multi-relay networks with imperfect CSI
Wei Duan1,2,
Xiaojun Zhu1,4,
Li Jin1,
Wei Wang1,4,
Guoan Zhang1 &
Jeaho Choi3
This paper investigates decoupled beamforming techniques for the distributed multi-input multi-output (MIMO) two-way relay networks (TWRN) with imperfect channel state informations (CSIs). The objective of this paper is to maximize the weighted sum rate (SR) with semi-infinite relay power constraints. Since the objective problem is difficult to be solved directly, considering the high signal-to-residual-interference-plus-noise ratio (SRINR) and employing the Cauchy-Schwarz inequality and S-lemma, the problem can be approximately converted into a source beamforming decoupled one. In addition, with the optimal relay beamforming design and maximum ratio combining (MRC) at the receiver, the MIMO channels are decoupled into parallel single-input single-output (SISO) channels. By this way, the suboptimal relay beamforming matrix can be efficiently obtained with minimal relay power constraint. Specifically, the semi-infinite constraints can be reformulated into a linear matrix inequality (LMI), which can be efficiently solved by using an alternating optimization algorithm. Numerical results demonstrate that our proposed decoupled beamforming scheme outperforms the existing works in terms of the SR and the computational complexity with satisfactory convergence.
Recently, cooperative and multiple-input multiple-output (MIMO) systems [1–4] have been widely considered as a candidate for the fifth generation (5G) wireless communication due to their transmission reliability, where multiple users (MU) are simultaneously served at the same frequency band by the base station (BS) equipped with number of antennas. Cooperative two-way relay network (TWRN) technology has attracted significant interests due to the superior spectral efficiency. Several relaying schemes haven proposed, i.e., amplify-and-forward (AF) [5–8], decode-and-forward (DF) [9, 10], and denoise-and-forward (DNF) [11, 12]. Specifically, the cooperative relaying beamforming design has been studied in [13–16], where the authors showed the optimal solutions of the source and/or relay beamforming matrices. Moreover, the spatially correlated fading channels are considered in [17, 18], which is a more practical assumption but makes the training problem more challenging. The authors of [17] presented training designs for estimation of spatially correlated MIMO AF two-way multi-relay channels, where an optimal training structure is initially derived to minimize total mean-square-error (MSE) of the channel estimation. Based on [17], an optimal training scheme is efficiently designed to minimize the total MSE of the channel estimation under the transmit power constraints at the source nodes and at the relay in [18]. In addition, the authors considered a TWRN with an AF protocol over either two, three, or four time slots [19].
Considering inaccurate channel estimation and feedback delay, the perfect channel state information (CSI), which is proposed in the above works [5–18], is usually hard to obtain in practice. By taking account into the channel uncertainties, the imperfect CSI scenario has been studied in [20–25]. In [20], joint relay and jammer selection and power control for physical layer security issues in two-way relay networks are studied to maximize the secrecy capacity of the network. In [21], the authors addressed the robust multiple-antenna relay design problem in TWRN and provided the robust multiple-antenna relay design based on the channel estimates. A low complexity processing matrices are presented for a multi-pair two-way massive MIMO AF full-duplex relay system [22]. Specifically, in [23], the problem of optimal beamforming and power allocation for an AF-based TWRN is studied in the presence of interference and CSI uncertainty, where two different approaches, namely the total power minimization method and the signal-to-interference-plus-noise-ratio (SINR) balancing technique, are proposed. Particularly, the robust two-way relay precoder design for a cognitive radio network is investigated in [24], where two different types of CSI errors with corresponding robust designs are proposed. It is worth noting, in [25], the joint optimal robust beamforming designs in multi-pair two-way nonregenerative relaying systems, where an insightful closed form solution is obtained that is not only robust to the imperfect CSI but also adjustable to various CSI circumstances at the users and the relay.
However, most existing works for the robust two-way relaying networks were unavoidable to focus on the joint design of source and relay beamforming matrices that may lead to a higher computational complexity. In addition, as shown in [26], for the multiple relays scheme, the performance of the capacity outperforms that of the single relay one. On the other hand, multiple relays scheme is more practical and challenging for the wireless communication scenarios. Moreover, the weighed scheme can be regarded as one kind of the resource allocations, which is general. Therefore, in order to reduce the complexity, our study aims to design a decoupled beamforming scheme for the distributed MIMO two-way relaying scheme with imperfect and reciprocal CSI. By converting the weighted sumrate (SR) maximization problem into a source beamforming decoupled one whose target is to maximize the signal-to-residual-interference-plus-noise ratio (SRINR) by means of approximations, the S-lemma, and the Cauchy-Schwarz inequality, the objective problem can be efficiently solved by an alternating algorithm. With the decoupled source beamforming design and the maximum ratio combining (MRC) at the receiver, the optimal relay beamforming can be efficiently obtained which leads to that the MIMO channels can be decoupled into parallel single-input single-output (SISO) channels. Numerical results are presented to corroborate that the performance of the decoupled scheme is improved and to compare it to the existing works with satisfactory convergence.
The rest of this paper is organized as follows. Section 2 describes the system model of the TWRN and objective problems. In Section 3, the proposed decoupled beamforming and optimal relay beamforming designs are investigated. Numerical results are presented to show the excellent performance of our proposed scheme in Section 4. Section 5 concludes this paper.
Notations: For an M×N matrix A, E(A), ∥A∥, AT, tr(A), vec(A), and A† denote the statistical expectation, Frobenius norm, transpose, trace, vectorization, and Hermitian transpose of A, respectively. I N represents an N×N identity matrix.
System model and optimal relay beamforming design
System model
Consider a simple two-way relay channels (TWRC) consisting of two user equipments, i.e., UE1 and UE2, and L relay nodes {R1,R2,…,R L } as shown in Fig. 1. Each source and relay nodes are equipped with M and N antennas, respectively. Assume that all nodes operate in half-duplex mode and the direct link between two UEs does not exist. By taking the estimation error into account, the CSI is assumed to be reciprocal and partially known at each node. Denote the estimated channels from the UE t to the relay R i , from the relay R i to the UE t by $\widetilde {\mathbf {F}}_{i,t}\in {\mathbb {C}^{N\times {M}}}$ and $\widetilde {\mathbf {G}}_{t,i}\in {\mathbb {C}^{M\times {N}}}$, for t∈{1,2}, respectively, as follows:
$$\begin{array}{@{}rcl@{}} \mathbf{F}_{i,t}\triangleq\left\{\widetilde{\mathbf{F}}_{i,t}+\Delta_{\mathbf{F}_{i,t}}\right\} \end{array} $$
Two-way MIMO multi-relay networks
$$\begin{array}{@{}rcl@{}} \mathbf{G}_{t,i}\triangleq\left\{\widetilde{\mathbf{G}}_{t,i}+\Delta_{\mathbf{G}_{t,i}}\right\} \end{array} $$
with $\left \{\Delta _{\mathbf {F}_{i,t}},~\Delta _{\mathbf {G}_{t,i}}\right \}$ as the channel uncertainties. For a better readability and distinction, in this paper, we simply consider two channel matrices (F i,t , G t,i ) with $\mathbf {F}_{i,t}=\mathbf {G}^{T}_{t,i}$ and $\mathbf {F}_{i,\overline {t}}=\mathbf {G}^{T}_{\overline {t},i}$. Since that, the channel uncertainty is provided by the inaccurate channel estimation, feedback delay, and so on. Therefore, the norm bounded error (NBE) of the channel uncertainty is normally considered as a small slack value [27], which is a reasonable assumption in a practical system. Therefore, in this paper, the channel uncertainties are considered as a NBE model for its simplicity:
$$\begin{array}{@{}rcl@{}} \left\|\Delta_{\mathbf{F}_{i,1}}\right\|=\left\|\Delta_{\mathbf{G}_{1,i}}\right\|=\alpha_{i},~ \left\|\Delta_{\mathbf{F}_{i,2}}\right\|=\left\|\Delta_{\mathbf{G}_{2,i}}\right\|=\beta_{i}, \end{array} $$
where 0≤{α i ,β i }≪1. By this way, a worst-case design methodology can be adopted, resulting in the proposed system design which maximizes performance for the worst possible CSI realization as defined by the NBE.
It is assumed that two UEs exchange the message using two consecutive time slots. At the first time slot, the information data x t =[x1,t,x2,t,…,x m,t ]T, with $\mathrm {E}\left (\mathbf {x}_{t}\mathbf {x}_{t}^{\dag }\right)=\mathbf {I}_{M}$, for i∈{1,…,M}, is linearly processed by a precoding matrix $\mathbf {B}_{t}\in {\mathbb {C}^{M\times {M}}}$ with the power constraint ∥B t ∥2≤P t and then transmitted to the relay nodes. Thus, the received signal at R i can be expressed as
$$\begin{array}{@{}rcl@{}} \mathbf{y}_{R_{i}}=\mathbf{F}_{i,1}\mathbf{B}_{1}\mathbf{x}_{1}+\mathbf{F}_{i,2}\mathbf{B}_{2}\mathbf{x}_{2}+\mathbf{n}_{R_{i}}, \end{array} $$
where $\mathbf {n}_{R_{i}}\sim {\mathcal {CN}\left (0,\sigma ^{2}_{R_{i}}\mathbf {I}_{N}\right)}$ denotes the additive white Gaussian noise (AWGN) vector with zero mean and variance $\sigma ^{2}_{R_{i}}\mathbf {I}_{N}$.
At the second time slot, the relay node R i generates the transmit signal $\phantom {\dot {i}\!}\mathbf {W}_{i}\mathbf {y}_{R_{i}}$, for ∥W i ∥2≤ω i , with the N×N beamforming matrix W i and forwards to the UE1 and UE2 which results in $\phantom {\dot {i}\!}\mathbf {x}_{R_{i}}=\mathbf {W}_{i}\mathbf {y}_{R_{i}}$ with the average transmit power as
$$\begin{array}{@{}rcl@{}} \mathrm{E}\left\{\left\|\mathbf{x}_{R_{i}}\right\|^{2}_{2}\right\} \,=\, \left\|\mathbf{W}_{i}\mathbf{F}_{i,1}\mathbf{B}_{1}\right\|^{2} \,+\, \left\|\mathbf{W}_{i}\mathbf{F}_{i,2}\mathbf{B}_{2}\right\|^{2}+\sigma^{2}_{R_{i}}\left\|\mathbf{W}_{i}\right\|^{2},\\ \end{array} $$
where we have assumed that the elements of the data stream vector x t transmitted by each user are independent, i.e., $\mathbf {E}\left \{\mathbf {x}_{t}\mathbf {x}_{t}^{\dag }\right \}=\mathbf {I}_{M}$. Thus, the received signal at source node S t , for t∈{1,2}, becomes
$$\begin{array}{@{}rcl@{}} \mathbf{y}_{t}=\sum^{L}_{i=1}\mathbf{G}_{t,i}\mathbf{W}_{i}\mathbf{y}_{R_{i}}+\mathbf{n}_{t}, \end{array} $$
where n t denotes the noise vector at the source node S t with zero mean and variance $\sigma _{S_{t}}^{2}\mathbf {I}_{N}$.
In practical systems, however, achieving perfect SIC is a hard problem due to channel estimation errors. Thus, the self-interference cannot be completely eliminated which leads to the residual self-interference cannot be ignorable. In order to subside the adverse effect of the interference, with the estimated channel coefficients, the final observation of the residual self-interference at the $\text {UE}_{\overline {t}}$ can be obtained as
$$\begin{array}{@{}rcl@{}} \chi_{\overline{t}}&=&\sum_{i=1}^{L}\left(\mathbf{G}_{\overline{t},i} \mathbf{W}_{i}\mathbf{F}_{i,\overline{t}}\mathbf{B}_{\overline{t}}\mathbf{x}_{\overline{t}} - \widetilde{\mathbf{G}}_{\overline{t},i}\mathbf{W}_{i}\widetilde{\mathbf{F}}_{i,\overline{t}} \mathbf{B}_{\overline{t}}\mathbf{x}_{\overline{t}}\right)\\ &=&\!\sum_{i=1}^{L}\!\left(\widetilde{\mathbf{G}}_{\overline{t},i} \mathbf{W}_{i}\Delta_{\mathbf{F}_{i,\overline{t}}}\,+\,\Delta_{\mathbf{G}_{\overline{t},i}} \mathbf{W}_{i}\widetilde{\mathbf{F}}_{i,\overline{t}}\,+\,\underbrace{\Delta_{\mathbf{G}_{\overline{t},i}} \mathbf{W}_{i}\Delta_{\mathbf{F}_{i,\overline{t}}}}_{\mathcal{K}}\right)\!\mathbf{B}_{\overline{t}}\mathbf{x}_{\overline{t}}.\\ \end{array} $$
It is clear that the covariance of the term $\mathcal {K}$ can be obtained as $\left \|\Delta _{\mathbf {G}_{\overline {t},i}}\mathbf {W}_{i}\Delta _{\mathbf {F}_{i,\overline {t}}}\right \|^{2}$. Since that, the norm bounded errors of $\left \{\Delta _{\mathbf {G}_{\overline {t},i}}, \Delta _{\mathbf {F}_{i,\overline {t}}}\right \}$ are assumed to be small slack values, which is a reasonable assumption in a practical system, i.e., 0≤{α i ,β i }≪1. Using ∥AB∥≤∥A∥∥B∥, we have $\left \|\Delta _{\mathbf {G}_{\overline {t},i}}\mathbf {W}_{i}\Delta _{\mathbf {F}_{i,\overline {t}}}\right \|^{2}\leq \left \{\alpha ^{4}_{i}\left \|\mathbf {W}_{i}\right \|^{2},\beta ^{4}_{i}\left \|\mathbf {W}_{i}\right \|^{2}\right \}$, which is very close to 0. Therefore, in this paper, the negligible term involving only CSI uncertainty is omitted, which is practical and relatively easy to achieve.
Since the sum of the individual rate in a MIMO system is difficult to be obtained, therefore, we turn to design the relay beamforming to recast the optimization problem into the one of achievable data rate for each data stream.
Optimal relay beamforming design
Suppose the singular value decomposition (SVD) of $\widetilde {\mathbf {G}}_{\overline {t},i}$ and $\widetilde {\mathbf {F}}_{i,t}$ as follows
$$\begin{array}{@{}rcl@{}} \widetilde{\mathbf{G}}_{\overline{t},i}=\mathbf{\Pi}_{\overline{t},i}{\left[\begin{array}{cc} \left[\mathbf{G}^{\diamondsuit}_{\overline{t},i}\right]_{M\times M} &\mathbf{0}_{M\times (N-M)}\end{array}\right]}\mathbf{V}_{\overline{t},i}^{\dag}\triangleq \mathbf{\Pi}_{\overline{t},i}\mathbf{\Sigma}_{\overline{t},i}\mathbf{V}_{\overline{t},i}^{\dag},\\ \end{array} $$
$$\begin{array}{@{}rcl@{}} \widetilde{\mathbf{F}}_{i,t}=\mathbf{U}_{i,t}{\left[\begin{array}{c} \left[\mathbf{F}^{\diamondsuit}_{i,t}\right]_{M\times M}\\ \mathbf{0}_{(N-M)\times M}\end{array}\right]}\mathbf{\Omega}^{\dag}_{i,t}\triangleq \mathbf{U}_{i,t}\mathbf{\Gamma}_{i,t}\mathbf{\Omega}^{\dag}_{i,t}, \end{array} $$
where $\mathbf {V}_{\overline {t},i}\in {\mathbb {C}^{N\times {N}}}$, $\mathbf {U}_{i,t}\in {\mathbb {C}^{N\times {N}}}$, $\mathbf {\Pi }_{\overline {t},i}\in {\mathbb {C}^{M\times {M}}}$, and $\mathbf {\Omega }_{i,t}\in {\mathbb {C}^{M\times {M}}}$ are unitary matrices as well as $\mathbf {\Sigma }_{\overline {t},i}={\left [\begin {array}{cc} \left [\mathbf {G}^{\diamondsuit }_{\overline {t},i}\right ]_{M\times M} &\mathbf {0}_{M\times (N-M)}\end {array}\right ]} $ and $\mathbf {\Gamma }_{i,t}={\left [\begin {array}{c} \left [\mathbf {F}^{\diamondsuit }_{i,t}\right ]_{M\times M}\\ \mathbf {0}_{(N-M)\times M}\end {array}\right ]}. $
Using the SVDs in (7) and (8), the optimal relay beamforming matrix as the solution to the objective problem T3 can be obtained as
$$\begin{array}{@{}rcl@{}} \mathbf{W}_{i}=\mathbf{V}_{\overline{t},i}\mathbf{C}_{i}\mathbf{U}_{i,t}^{\dag}, \end{array} $$
where $\mathbf {C}_{i}\in {\mathbb {C}^{N\times {N}}}$ is a matrix to be determined.
The proof is similar to [8, 28]. □
Without loss of generality, C i can be further partitioned as follows:
$$\begin{array}{@{}rcl@{}} \mathbf{C}_{i}={\left[\begin{array}{cc} \left[\mathbf{W}^{\diamondsuit}_{i}\right]_{M\times M} & \mathbf{X}_{M\times(N-M)}\\ \mathbf{Y}_{(N-M)\times M} & \mathbf{Z}_{(N-M)\times(N-M)} \end{array}\right]}. \end{array} $$
On the one hand, since the channels from different relays, the estimated channel, and its corresponding channel uncertainties are independent, denoting $\widetilde {\chi }_{\overline {t}}=\sum _{i=1}^{L} \left (\widetilde {\mathbf {G}}_{\overline {t},i}\mathbf {W}_{i}\Delta _{\mathbf {F}_{i,\overline {t}}} + \Delta _{\mathbf {G}_{\overline {t},i}}\mathbf {W}_{i}\widetilde {\mathbf {F}}_{i,\overline {t}}\right)$, we have
$$ \begin{aligned} \left\|\widetilde{\chi}_{\overline{t}}\right\|^{2}&\leq\sum_{i=1}^{L}\varrho^{2}_{\overline{t},i}\left(\left\|\mathbf{G}^{\diamondsuit}_{\overline{t},i}\right\|^{2}+\left\|\mathbf{F}^{\diamondsuit}_{i,\overline{t}}\right\|^{2} \right)\\ &\quad\times\left\|{\left[\begin{array}{cc} \left[\mathbf{W}^{\diamondsuit}_{i}\right]_{M\times M} & \mathbf{X}_{M\times(N-M)}\\ \mathbf{Y}_{(N-M)\times M} & \mathbf{Z}_{(N-M)\times(N-M)} \end{array}\right]}\right\|^{2}, \end{aligned} $$
where ϱ t,i =α i for t=1 as long as ϱ t,i =β i for t=2. On the other hand, since the channel uncertainty $\phantom {\dot {i}\!}\Delta _{\mathbf {F}_{i,t}}$ can be further partitioned following
$$\begin{array}{@{}rcl@{}} \Delta_{\mathbf{F}_{i,t}}&=&{\left[\begin{array}{cc} \left[\Delta_{\mathbf{F}_{i,t}}^{\diamondsuit}\right]_{M\times M}\\ \left[\Delta_{\mathbf{F}_{i,t}}^{\flat}\right]_{(N-M)\times M} \end{array}\right]}, \end{array} $$
upon substituting (9), (10), and (12) into the relay power constraint (4), we have
$$ \begin{aligned} &\mathbf{E}\left\{\left\|\mathbf{x}_{R_{i}}\right\|^{2}_{2}\right\}\\ &=\text{tr} \left\{\mathbf{W}_{i}\mathbf{y}_{R_{i}}\mathbf{y}^{\dag}_{R_{i}}\mathbf{W}_{i}^{\dag}\right\}\\ &=\text{tr}\left\{\mathbf{W}_{i}\left(\mathbf{F}_{i,1}\mathbf{B}_{1}\mathbf{B}_{1}^{\dag}\mathbf{F}_{i,1}^{\dag}+ \mathbf{F}_{i,2}\mathbf{B}_{2}\mathbf{B}_{2}^{\dag}\mathbf{F}_{i,2}^{\dag}+\sigma^{2}_{R_{i}}\mathbf{I}_{N}\right) \mathbf{W}_{i}^{\dag}\right\}\\ &=\text{tr}\!\left\{\mathbf{W}^{\diamondsuit}_{i}\!\!\left(\sum_{t=1}^{2}\mathbf{F}^{\sharp}_{i,t}\mathbf{B}_{t} \!\left(\mathbf{F}^{\sharp}_{i,t}\mathbf{B}_{t}\right)^{\dag}\,+\,\sigma^{2}_{R_{i}}\mathbf{I}_{M}\right)\!\! \left(\mathbf{W}^{\diamondsuit}_{i}\right)^{\dag}\,+\,\mathbf{X}\mathbf{\Phi}_{i}\mathbf{X}^{\dag}\right\} \\ &\quad+\text{tr}\left\{\mathbf{Z}\mathbf{\Phi}_{i}\mathbf{Z}^{\dag}+\mathbf{Y}\left(\sum_{t=1}^{2}\mathbf{F}^{\sharp}_{i,t}\mathbf{B}_{t} \left(\mathbf{F}^{\sharp}_{i,t}\mathbf{B}_{t}\right)^{\dag} +\sigma^{2}_{R_{i}}\mathbf{I}_{M}\right)\mathbf{Y}^{\dag}\right\}, \end{aligned} $$
where $\mathbf {F}^{\sharp }_{i,t}\,=\,\mathbf {F}^{\diamondsuit }_{i,t}\,+\,\Delta ^{\diamondsuit }_{\mathbf {F}_{i,t}} \text {and}\,\mathbf {\Phi }_{i}=\sum _{t=1}^{2}\Delta _{\mathbf {F}_{i,t}}^{\flat }\mathbf {B}_{t}\left (\Delta _{\mathbf {F}_{i,t}}^{\flat }\mathbf {B}_{t}\right)^{\dag } +\sigma ^{2}_{R_{i}}\mathbf {I}_{N-M}$. From (11) and (13), it is clear that, for any feasible C i with {X,Y,Z}≠0, one can always find
$$\begin{array}{@{}rcl@{}} \mathbf{C}^{\prime}_{i}={\left[\begin{array}{cc} \left[\mathbf{W}^{\diamondsuit}_{i}\right]_{M\times M} & \mathbf{0}_{M\times(N-M)}\\ \mathbf{0}_{(N-M)\times M} & \mathbf{0}_{(N-M)\times(N-M)} \end{array}\right]}, \end{array} $$
which can achieve the smaller relay and interference power constraints. By this way, the optimal expression of W i with minimal relay and interference power constraints can be represented as
$$\begin{array}{@{}rcl@{}} \mathbf{W}^{\sharp}_{i}=\mathbf{V}_{\overline{t},i}{\left[\begin{array}{cc} \left[\mathbf{W}^{\diamondsuit}_{i}\right]_{M\times M} & \mathbf{0}_{M\times(N-M)}\\ \mathbf{0}_{(N-M)\times M} & \mathbf{0}_{(N-M)\times(N-M)} \end{array}\right]}\mathbf{U}_{i,t}^{\dag}. \end{array} $$
Objective problem
It is easy to see that the actual equivalent channel from UE t to $\text {UE}_{\overline {t}}$ is with dimension M×M. Therefore, by employing the MRC at $\text {UE}_{\overline {t}}$, the equivalent received signals from each channel can be linearly added together. With the condition $\mathrm {E}\left (\mathbf {x}_{t}\mathbf {x}_{t}^{\dag }\right)=\mathbf {I}_{M}$, the objective weighted sum-rate (SR) problem for the total system can be recast into the problem achievable data rate for each data stream, i.e., the MIMO channels are decoupled into M parallel SISO channels. With these observations, from (4), (6), and (15) similar to [29–31], the received SRINR over the mth channel (1≤m≤M) for $\text {UE}_{\overline {t}}$ can be simply written as
$$\begin{array}{@{}rcl@{}} \zeta_{\overline{t}}^{(m)}=\frac{\left\|\sum_{i=1}^{L}\mathbf{G}_{\overline{t},i}\mathbf{W}_{i} \mathbf{F}_{i,t}\mathbf{B}_{t}\right\|{~}^{2}}{\left\|\widetilde{\chi}_{\overline{t}}\mathbf{B}_{\overline{t}}\right\|{~}^{2}+\sigma^{2}_{R_{i}}\left\|\sum_{i=1}^{L}\mathbf{G}_{\overline{t},i}\mathbf{W}_{i}\right\|{~}^{2}+\sigma_{S_{t}}^{2}}, \end{array} $$
where $\overline {t}=2$ for t=1 while $\overline {t}=1$ for t=2. Based on which, we have the weighted SR maximization problem as
$$ \begin{array}{ccl} T_{1}: & \max\limits_{\mathbf{B}_{t},\mathbf{W}_{i}} & \sum\limits_{t=1}^{2}\lambda_{\overline{t}}\text{log}\left(1+\frac{\left\|\mathbf{H}_{t}\mathbf{B}_{t}\right\|{~}^{2}} {\left\|\widetilde{\chi}_{\overline{t}}\mathbf{B}_{\overline{t}}\right\|{~}^{2}+N_{t}}\right)\\ & s.t. & \text{(4)} \leq P_{R_{i}},~\left\|\mathbf{B}_{t}\right\|{~}^{2}\leq P_{t}, \forall t=\left\{1,2\right\}, \end{array} $$
where $\mathbf {H}_{t}=\sum _{i=1}^{L}\mathbf {G}_{\overline {t},i}\mathbf {W}_{i}\mathbf {F}_{i,t}$, $N_{t}=\sigma ^{2}_{R_{i}}\left \|\sum _{i=1}^{L}\mathbf {G}_{\overline {t},i}\mathbf {W}_{i}\right \|^{2}+\sigma _{S_{t}}^{2}$, and $P_{R_{i}}$ is the maximum allocated power to the relay. In addition, the weight term λ t is determined depending on the required quality of service (QoS).
Decoupled source beamforming designs
Since the problem T1 is neither convex nor concave and the optimal B t is intractable, it is difficult to obtain the globally optimal solution. In this section, we propose a decoupled way to obtain the sub-optimal solution of the worst-case weighted SR.
Considering the high signal-to-interference-plus-noise ratio (SINR), the objective weighted SR Rsum in (17) can be approximately reexpressed as
$$\begin{array}{@{}rcl@{}} R_{\text{sum}}&\approx&\lambda_{1}\text{log}\frac{\left\|\mathbf{H}_{2}\mathbf{B}_{2}\right\|{~}^{2}}{\left\|\widetilde{\chi}_{1} \mathbf{B}_{1}\right\|{~}^{2}+N_{1}}+\lambda_{2}\text{log}\frac{\left\|\mathbf{H}_{1}\mathbf{B}_{1}\right\|{~}^{2}}{\left\|\widetilde{\chi}_{2} \mathbf{B}_{2}\right\|{~}^{2}+N_{2}}\\ &=&\text{log}\frac{\left(\left\|\mathbf{H}_{2}\mathbf{B}_{2}\right\|{~}^{2}\right){~}^{\lambda_{1}}}{\left(\left\|\widetilde{\chi}_{2}\mathbf{B}_{2} \right\|{~}^{2}+N_{2}\right)^{\lambda_{2}}}+\text{log}\frac{\left(\left\|\mathbf{H}_{1}\mathbf{B}_{1}\right\|{~}^{2}\right){~}^{\lambda_{2}}}{\left(\left\|\widetilde{\chi}_{1}\mathbf{B}_{1}\right\|{~}^{2}+N_{1}\right)^{\lambda_{1}}}.\\ \end{array} $$
It is interesting to see that, since log(·) is a monotonic function and assuming the nominal CSI and the channel errors are independent, with fixed λ t , from (31), the objective problem T1 can be equivalently converted into the following problem as
$$\begin{array}{@{}rcl@{}} T_{2}: & \max\limits_{\mathbf{B}_{t},\mathbf{W}_{i}} & \prod_{t=1}^{2}\gamma_{t}\ \ \text{for}\ t=1,2,\\ & s.t. & q_{t}(\mathbf{\mathbf{B}}_{t})\geq\gamma_{t}, \text{(4)}\leq P_{R_{i}},\left\|\mathbf{B}_{t}\right\|^{2}\leq P_{t}, \end{array} $$
where $q_{t}(\mathbf {B}_{t})=\frac {\left (\left \|\mathbf {H}_{t}\mathbf {B}_{t}\right \|^{2}\right)^{\lambda _{\overline {t}}}}{\left (\|\widetilde {\chi }_{t}\mathbf {B}_{t}\|^{2}+N_{t}\right)^{\lambda _{t}}}$ and γ t is a slack value, where, similar to [29, 32–34], by introducing a slack value a, the optimization problem $\left \{\text {max}~\mathcal {A}\right \}$ can be equivalently rewrite as $\left \{\text {max}~a\right \}~s.t. \mathcal {A}\geq a$. In order to relax N t in q t (B t ), after introducing the slack value ς t , we have
$$\begin{array}{@{}rcl@{}} \sum_{i=1}^{L}\left(\varrho^{2}_{t,i}+\left\|\widetilde{\mathbf{G}}_{\overline{t},i}\mathbf{W}_{i}\right\|^{2}\right) \geq\frac{\varsigma_{t}-\sigma_{1}^{2}}{\sigma^{2}_{R_{i}}} \end{array} $$
which leads to N t ≥ς t , where ϱ t,i =α i for t=1 as long as ϱ t,i =β i for t=2, and τ t is a slack value. With this observation and employing the Cauchy-Schwarz inequality, we have
$$\begin{array}{@{}rcl@{}} \frac{\left(\left\|\mathbf{H}_{t}\mathbf{B}_{t}\right\|^{2}\right)^{\lambda_{\overline{t}}}}{\left(\|\widetilde{\chi}_{t} \mathbf{B}_{t}\|^{2}+N_{t}\right)^{\lambda_{t}}}&=&\frac{\text{tr}\left(\mathbf{B}^{\dag}_{t}\mathbf{H}^{\dag}_{t}\mathbf{H}_{t} \mathbf{B}_{t}\right)^{\lambda_{\overline{t}}}}{\text{tr}\left[\mathbf{B}^{\dag}_{t}\left(\widetilde{\chi}^{\dag}_{t}\widetilde{\chi}_{t}+ \frac{\varsigma_{t}}{mP_{t}}\mathbf{I}_{M}\right)\mathbf{B}_{t}\right]^{\lambda_{t}}}\\ &\leq&\frac{\text{tr}\left(\mathbf{H}^{\dag}_{t}\mathbf{H}_{t}\right)^{\lambda_{\overline{t}}}}{\text{tr}\left[\left(\widetilde{\chi}^{\dag}_{t} \widetilde{\chi}_{t}+\frac{\varsigma_{t}}{mP_{t}}\mathbf{I}_{M}\right)\right]^{\lambda_{t}}P^{2(\lambda_{t}-\lambda_{\overline{t}})}_{t}}\\ &=&\frac{\left(\left\|\mathbf{H}_{t}\right\|^{2}\right)^{\lambda_{\overline{t}}}}{\left(\|\widetilde{\chi}_{t}\|^{2}+\frac{\varsigma_{t}}{P_{t}}\right)^{\lambda_{t}}P^{2(\lambda_{t}-\lambda_{\overline{t}})}_{t}}, \end{array} $$
which decouples the source beamforming matrix B t . Denoting that $\frac {\left (\left \|\mathbf {H}_{t}\right \|^{2}\right)^{\lambda _{\overline {t}}}}{\left (\|\widetilde {\chi }_{t}\|^{2}+\frac {\varsigma _{t}}{P_{t}}\right)^{\lambda _{t}}} $ and $\widehat {\gamma }_{t}=\gamma _{t}P^{2(\lambda _{t}-\lambda _{\overline {t}})}_{t}$, for the objective problem T2, the SRINR constraint can be equivalently expressed as $\widehat {q}_{t}(\mathbf {B}_{t})\geq \widehat {\gamma }_{t}$. It is clear that, after decoupling the source beamdforming B t in the objective problem T2, B t only exists at the relay power constraint (4). By considering $\mathbf {\widehat {B}}_{t}=\text {vec}(\widetilde {\mathbf {B}}_{t})\text {vec}(\widetilde {\mathbf {B}}_{t})^{\dag }$ with $\text {vec}(\widetilde {\mathbf {B}}_{t})=\left [\text {vec}(\mathbf {B}_{t})^{T}~~ \mathbf {0}_{1\times M(N-M)}\right ]^{T}$, similar to Appendix 2 in [8], we have the following linear matrix inequality (LMI) of the individual relay power constraints:
$$\begin{array}{@{}rcl@{}} {\left[\begin{array}{ccc} \!\!\mathbf{\widehat{B}}_{1}+\lambda_{1}\mathbf{I} & \!\!\mathbf{\widehat{B}}_{1}\text{vec}\left(\widetilde{\mathbf{Q}}_{1}\right) & \!\!\mathbf{0}\\ \!\! \text{vec}\left(\widetilde{\mathbf{Q}}_{1}\right)^{\dag}\mathbf{\widehat{B}}_{1} & \!\!\aleph_{t} & \!\! \text{vec}\left(\widetilde{\mathbf{Q}}_{2}\right)^{\dag}\mathbf{\widehat{B}}_{2}\\ \!\! \mathbf{0} & \!\!\mathbf{\widehat{B}}_{2}\text{vec}\left(\widetilde{\mathbf{Q}}_{2}\right) & \!\!\mathbf{\widehat{B}}_{2}+\lambda_{2}\mathbf{I} \end{array}\right]} \!\!\succeq0, \end{array} $$
where $\aleph _{t}=\sum _{t=1}^{2}\text {vec}{\left (\widetilde {\mathbf {F}}_{i,t}\right)}^{\dag }\mathbf {\widehat {B}}_{t}\text {vec}\left (\widetilde {\mathbf {F}}_{i,t}\right)-\lambda _{1}\omega _{i}\alpha _{i}^{2}-\lambda _{2}\omega _{i}\beta _{i}^{2}-\widehat {P}_{R_{i}}$ and $\widetilde {\mathbf {Q}}_{t}=\mathbf {W}_{i}\widetilde {\mathbf {F}}_{i,t}$. Based on which, the optimization problem T2 can be approximately converted into the following version as
$$\begin{array}{@{}rcl@{}} T_{3}: & \max\limits_{\mathbf{B}_{t},\mathbf{W}_{i}} & \qquad \prod_{t=1}^{2}\widehat{\gamma}_{t}\ \text{for}~~t=1,2,\\ &s.t.&\ \ \widehat{q}_{t}(\mathbf{\mathbf{B}}_{t})\geq\widehat{\gamma}_{t}, \text{(21)},\left\|\mathbf{B}_{t}\right\|^{2}\leq P_{t}. \end{array} $$
It is easy to see that the problem T3 is still hard to obtain the optimal solution straightforwardly, therefore, we try to convert it into LMI version. It is clear that the denominator of $\widehat {q}_{t}(\mathbf {B}_{t})$ can be approximately reexpressed as
$$\begin{array}{@{}rcl@{}} \tau_{t}+\frac{\varsigma_{t}}{P_{t}}, \end{array} $$
where τ t is obtained from
$$\begin{array}{@{}rcl@{}} \left\|\widetilde{\chi}_{\overline{t}}\right\|^{2} &\leq&\left\|\sum_{i=1}^{L}\widetilde{\mathbf{G}}_{\overline{t},i}\mathbf{W}^{\sharp}_{i}\Delta_{\mathbf{F}_{i,\overline{t}}}\right\|^{2}+\left\|\sum_{i=1}^{L}\Delta_{\mathbf{G}_{\overline{t},i}}\mathbf{W}^{\sharp}_{i}\widetilde{\mathbf{F}}_{i,\overline{t}}\right\|^{2}\\ &\leq&2\sum_{i=1}^{L}\varrho_{t,i}^{2}\omega_{i}\left\|\widetilde{\mathbf{G}}_{\overline{t},i}\right\|^{2}=\tau_{t}. \end{array} $$
Thus, letting $\varepsilon _{t}=\sqrt [\lambda _{\overline {t}}]{\widehat {\gamma }_{t}\left (\tau _{t}+\frac {\varsigma _{t}}{P_{t}}\right)^{\lambda _{t}}}$, the SRINR constraint becomes ∥H t ∥2≥ε t which has the LMI form as follows:
$$\begin{array}{@{}rcl@{}} {\left[\begin{array}{cc} -\varepsilon_{t} & -\mathbf{H}_{t} \\ \mathbf{H}^{\dag}_{t} & \mathbf{I}_{M} \end{array}\right]\succeq0}. \end{array} $$
$$\begin{array}{@{}rcl@{}} \mathbf{H}_{t}=\mathbf{\Lambda}_{\mathbf{G},\mathbf{F}}\!+\sum_{i=1}^{L}\mathbf{M}_{\mathbf{G}_{\overline{t},i}}\text{vec} \left(\Delta_{\mathbf{G}_{\overline{t},i}}\right)\,+\,\sum_{i=1}^{L}\mathbf{M}_{\mathbf{F}_{i,t}}\text{vec}\left(\Delta_{\mathbf{F}_{i,t}}\right),\\ \end{array} $$
$$\begin{array}{@{}rcl@{}} \mathbf{M}_{\mathbf{G}_{\overline{t},i}}&=&\sum_{i=1}^{L}\left(\mathbf{W}^{\sharp}_{i}\widetilde{\mathbf{F}}_{i,t}\right)^{T}\otimes \mathbf{I}, \end{array} $$
$$\begin{array}{@{}rcl@{}} \mathbf{\Lambda}_{\mathbf{G},\mathbf{F}}&=&\sum_{i=1}^{L}\text{vec}\left(\widetilde{\mathbf{G}}_{\overline{t},i}\mathbf{W}_{i}\widetilde{\mathbf{F}}_{i,t}\right) \end{array} $$
$$\begin{array}{@{}rcl@{}} \mathbf{M}_{\mathbf{F}_{i,t}}=\sum_{i=1}^{L}\mathbf{I}\otimes\left(\widetilde{\mathbf{G}}_{\overline{t},i}\mathbf{W}_{i}\right), \end{array} $$
substituting (25) back into (24) and using the S-Lemma [35], the LMI version of the SRINR constraint can be finally denoted as
$$\begin{array}{@{}rcl@{}} {\left[\begin{array}{cccc} \mathbf{\Theta}_{t}-\sum_{j=1}^{2l}\phi_{j}\mathbf{q}_{2}^{\dag}\mathbf{q}_{2} & \xi_{1,t}\mathbf{\Xi}_{1} & \cdots & \xi_{2l,t}\mathbf{\Xi}_{2l}\\ -\xi_{1,t}\mathbf{\Xi}^{\dag}_{1} & \phi_{1}\mathbf{I}_{MN} & \cdots & \mathbf{0}\\ \vdots & \vdots & \ddots & \vdots\\ -\xi_{2l,t}\mathbf{\Xi}^{\dag}_{2l} & \mathbf{0} & \cdots & \phi_{2l}\mathbf{I}_{MN}\\ \end{array}\right]}\succeq0, \end{array} $$
where $\mathbf {q}_{2}\,=\,\left [-1,\mathbf {0}_{1\times M^{2}}\right ]$ and ϕ j ≥0 are the slack variables for j=1,…,2l. $\mathbf {\Xi }_{i}\,=\,\left [\mathbf {0}_{M^{2}\times 1},\mathbf {M}_{\mathbf {G}_{\overline {t},i}}^{H}\right ]$, and $\xi _{i,t}=\left \|\Delta _{\mathbf {G}_{\overline {t},i}}\right \|$ for i=1,…,l, as well as $\mathbf {\Xi }_{i}\,=\,\left [\mathbf {0}_{M^{2}\times 1},\mathbf {M}_{\mathbf {F}_{i,t}}^{H}\right ]$, and $\xi _{i,t}=\left \|\Delta _{\mathbf {F}_{i,\overline {t}}}\right \|$ for i=l+1,…,2l, with
$$\begin{array}{@{}rcl@{}} \mathbf{\Theta}_{t}={\left[\begin{array}{cc} -\varepsilon_{t} & -\mathbf{\Lambda}_{\mathbf{G},\mathbf{F}} \\ \mathbf{\Lambda}_{\mathbf{G},\mathbf{F}}^{\dag} & \mathbf{I}_{M^{2}} \end{array}\right]}. \end{array} $$
Thus, optimization problem T3 is recast as
$$ \begin{aligned} &T_{4}:~\max\limits_{\mathbf{B}_{t},\mathbf{W}_{i}}\qquad\qquad\qquad\! \prod_{t=1}^{2}\widehat{\gamma}_{t}\\ &s.t.\quad \text{(21), (29)},\left\|\mathbf{B}_{t}\right\|^{2}\leq P_{t},\left\|\mathbf{W}_{i}\right\|^{2}\leq \omega_{i},\tau_{t}\geq0,\varsigma_{t}\geq0,\\ &\lambda_{t}\geq\!0,\nu_{i}\geq0,\phi_{j}\geq\!0,\forall t=1,2,\forall i=1,\ldots,l,\forall j\,=\,1,\ldots,2l. \end{aligned} $$
Since the objective problem T3 is a biconvex semi-definite program (SDP) which can be efficiently solved by the following alternating algorithm by using CVX [34]:
1. Initialize: ξ, N max , $\mathbf {B}_{t}^{(n)}$ and $\mathbf {W}^{(n)}_{i}$, set n=0;
2. Repeat: For n=0 to N max
1: with fixed $\mathbf {B}_{t}^{(n-1)}$ update $\mathbf {W}^{(n)}_{i}$ by solving T4;
2: for given $\mathbf {W}^{(n)}_{i}$ update $\mathbf {B}^{(n)}_{t}$ by solving T4;
3: If $\prod _{t=1}^{2}\widehat {\gamma }^{(n)}_{t}-\prod _{t=1}^{2}\widehat {\gamma }_{t}^{(n-1)}\leq \xi $, break;
By initializing the small ξ and setting the limitation of the number of iterations N max , the suboptimal solution can be obtained when $\prod _{t=1}^{2}\widehat {\gamma }^{(n)}_{t}-\prod _{t=1}^{2}\widehat {\gamma }_{t}^{(n-1)}\leq \xi $.
Numerical results
In this section, we examine the performance of the proposed decoupled scheme in terms of the average weighted SR and the convergent performance compared with the case with perfect CSI, the maximal-ratio transmission (MRT) and the SI-SRINR scheme in [8]. In all simulations, spatially uncorrelated Rayleigh fading channels with unit variance are assumed. The number of relay nodes is set to be L=3 with M=N=2, and the noise variance $\sigma _{R_{i}}^{2}$ and $\sigma _{S_{t}}^{2}$ are equally given as σ2=1. We further set up that, in the proposed alternating algorithm, Nmax=400, ξ=10−3, and the relay beamforming matrices are initialed at random, respectively. In the following figures, we use "Optimal," "Proposed perfect," "Proposed α i / β i ," and "SNR" to denote the performance upper bound, proposed scheme with perfect SCI, proposed scheme with NBEs of α i / β i , and transmit SNR, respectively. Moreover, the performance upper bound "Optimal" is obtained by using the exhaustive search as a benchmark.
Figure 2 depicts the weighted SR performance versus transmit SNR in two considered cases: (1)λ1=6,λ2=1; (2)λ1=3,λ2=1, which compares to the optimal, perfect CSI and MRT schemes, respectively, with α i =β i =0.01 and α i =β i =0.1. In addition, it is worth noting that the perfect one is obtained with the NBE as α i =β i =0 while the robust solutions are obtained by using the proposed alternating algorithm. It is clear that our proposed scheme has a relative advantage over the MRT maximization scheme and close to the optimal and perfect ones, especially for the smaller NBE case. In addition, the performance gain grows as the ratio between λ1 and λ2 increases. This is reasonable, since considering the high transmit SNR, (31) is very close to the original one which has the the negligible performance lose.
The average weighted SR versus the transmit SNR
Figure 3 compares the convergence performance of the weighted SR for the proposed method with the SI-SRINR in [8], where we have fixed $\alpha ^{2}_{i}=\beta ^{2}_{i}=0.01$. Two cases are considered: (1) λ1=6,λ2=1 with SNR=30 dB and (2) λ1=3,λ2=1 with SNR=20 dB. Results reveal that our proposal has an obvious satisfactory convergence over the SI-SINR scheme. This is because, for our proposed method, the weighted SR problem is approximately converted into a source beamforming decoupled one which significantly reduces the computational complexity. Specifically, by employing the Cauchy-Schwarz inequality and S-lemma, the optimal relay beamforming W i is efficiently relaxed in q t (B t ) which also simplifies the corresponding algorithm. In addition, it is observed that, for the higher transmit SNR, the convergence rate of our proposed scheme overwhelms the SI-SRINR one and both these two schemes provide the optimal solution with more iterations.
The average weighted SR versus the number of iterations
Figure 4 exhibits the average SR performance for the proposed method with the optimal scheme versus the transmission number L with fixed α i =β i =0.01 and λ1=λ2=1. Two cases are considered: (1) SNR=30 dB and (2) SNR=25 dB. It is clear that, from Fig. 3, it is easy to see that the performances of our proposed scheme have a negligible performance loss compared to the optimal scheme, especially for the smaller number of the relay, which supports the practical utility of our design. Remarkably, the performance gap between our proposed scheme and the optimal one is smaller for a lower SNR case.
The average weighted SR versus the number of relay L
Denoting d1,i and d2,i as the distance between {UE1, UE2} and the relay R i , respectively, the average transmit SNRs, which take into account the relay positions, can be written as
$$\begin{array}{@{}rcl@{}} \zeta^{*}_{{1,i}}=\frac{\overline{\zeta}_{{1,i}}}{d^{\delta}_{1,i}},~\zeta^{*}_{{2,i}}=\frac{\overline{\zeta}_{{2,i}}}{d^{\delta}_{2,i}}, \end{array} $$
where $\overline {\zeta }_{{1,i}}=\frac {P_{1}}{\sigma ^{2}_{R_{i}}}$, $\overline {\zeta }_{{2,i}}=\frac {P_{2}}{\sigma ^{2}_{R_{i}}}$ and δ is the path loss exponent. Figure 5 depicts the maximum weighted SR of our proposed decoupled beamforming design for different values of d1,i with fixed λ1=λ2=1, δ=3 in two considered system setups: (1) case 1: $\overline {\zeta }_{1}=20~\text {dB},~ \overline {\zeta }_{2}=15~\text {dB}$; (2) case 2: $\overline {\zeta }_{1}=20~\text {dB},~ \overline {\zeta }_{2}=25~\text {dB}$. For simplicity, we further assume that d1,1=d1,2=d1,3=d1 and d2,1=d2,2=d2,3=d2 with d1+d2=1. One can see that the achieved SR for case 2 shows a higher performance than case 1 due to a lager source power considered for UE2. Moreover, the optimal relay positions occur when d1 is close to 0.1 for case 1 and d1→0.95 for case 2. It is reasonable since that the power at UE1 is greater than the power at UE2 for case 1, while an opposite consideration is performed in case 2.
The achievable weighted SR versus the distance between UE t and R i
In this paper, we investigated the decoupled beamforming techniques in the distributed two-way multi-relay networks with imperfect and reciprocal CSI. In order to maximize the weighted SR, the objective problem is first converted into a decoupled one by employing the Cauchy-Schwarz inequality and S-lemma, and then, the optimal relay beamforming design is also investigated with the LMI versions of the semi-infinite constraints. Based on which, the objective problem can be efficiently solved by the proposed alternating algorithm. By means of the numerical results, it has been shown that our proposed scheme not only has the advantage in the term of the SR but also is with the satisfactory convergence compared to the existing works.
S Han, C Yang, M Bengtsson, User scheduling for cooperative base station transmission exploiting channel asymmetry. IEEE Trans. Commun. 61(4), 1426–1435 (2013).
H Lin, F Gao, S Jin, G Ye Li, A new view of multi-user hybrid massive MIMO: non-orthogonal angle division multiple access. IEEE J. Sel. Areas Commun. 35(10), 2268–2280 (2017).
H Xie, F Gao, S Jin. An overview of low-rank channel estimation for massive MIMO systems.IEEE Access.4(99), 7313–7321 (2016).
H Xie, F Gao, S Zhang, S Jin, A unified transmission strategy for TDD/FDD massive MIMO systems with spatial basis expansion model. IEEE Trans. Veh. Technol. 66(4), 3170–3184 (2017).
J Joung, J Choi, Linear precoder design for an AF two-way MIMO relay node with no source node precoding. IEEE Trans. Veh. Technol. 66(11), 10526–10531 (2017). to appear.
R Budhiraja, B Ramamurthi, Transceiver design for nonconcurrent two-way MIMO AF relaying with QoS guarantees. IEEE Trans. Veh. Technol. 65(12), 9651–9661 (2016).
Y Dai, X Dong, Power allocation for multi-pair massive MIMO two-way AF relaying with linear processing. IEEE Trans. Wireless Commun. 15(9), 5932–5946 (2016).
W Duan, M Wen, X Jiang, Y Yan, M Lee, Sum-rate maximization and robust beamforming design for MIMO two-way relay networks with reciprocal and imperfect CSI. EURASIP J. Commun. and Network. (2016). Atical ID: 157.
Q Cui, T Yuan, X Tao, AA Dowhuszko, R Jantti, Energy efficiency analysis of two-way DF relay system with non-ideal power amplifiers. IEEE Commun. Lett. 18(7), 1254–1257 (2014).
J Gao, SA Vorobyov, H Jiang, J Zhang, M Haardt, Sum-rate maximization with minimum power consumption for MIMO DF two-way relaying—part II: network optimization. IEEE Trans. Signal Process. 61(14), 3578–3591 (2013).
Z Zhao, M Peng, Z Ding, W Wang, H Chen, Denoise-and-forward network coding for two-way relay MIMO systems. IEEE Trans. Veh. Technol. 63(2), 775–788 (2014).
M Masjedi, A Hoseini, S Gazor, Non-coherent detection and denoise-and-forward two-way relay networks. IEEE Trans. Commun. 64(11), 4497–4505 (2016).
C Wang, H Chen, Q Yin, A Feng, AF Molisch, Multi-user two-way relay networks with distributed beamforming. IEEE Trans. Wireless Commun. 10(10), 3460–3471 (2011).
V Nassab, S Shahbazpanahi, A Grami, Optimal distributed beamforming for two-way relay networks. IEEE Trans. Signal Process. 58(3), 1238–1250 (2010).
Y Rong, Joint source and relay optimization for two-way linear nonregenerative MIMO relay communications. IEEE Trans. Signal Process. 60(12), 6533–6546 (2012).
W Cheng, M Ghogho, Q Huang, D Ma, J Wei, Maximizing the sum-rate of amplify-and-forward two-way relaying networks. IEEE Commun. Lett. 18(11), 635–638 (2011).
J-M Kang, H-M Kim, Training designs for estimation of spatially correlated fading channels in MIMO amplify-and-forward two-way multi-relay networks. IEEE Commun. Lett. 20(4), 772–775 (2016).
J-M Kang, I-M Kim, H-M Kim, Optimal training design for MIMO-OFDM two-way relay networks. IEEE Trans. Commun. 65(9), 3675–3690 (2017).
R Louie, Y Li, B Vucetic, Practical physical layer network coding for two-way relay channels: performance analysis and comparison. IEEE Trans. Wireless Commun. 9(2), 764–777 (2010).
F Jiang, C Zhu, J Peng, W Liu, Z Zhu, Y He, Joint relay and jammer selection and power control for physical layer security in two-way relay networks with imperfect CSI. Wireless Pers. Commun. 85:, 841–862 (2015).
C Wang, X Dong, Y Shi, Robust relay design for two-way multi-antenna relay systems with imperfect CSI. J. Commun. and Networks.16(1), 45–55 (2014).
F Jiang, C Zhu, J Peng, W Liu, Z Zhu, Y He, Multi-pair two-way massive MIMO AF full-duplex relaying with imperfect CSI over Ricean fading channels. IEEE Access. 4:, 4933–4945 (2016).
S Salari, MZ Amirani, I Kim, D Kim, J Yang, Distributed beamforming in two-way relay networks with interference and imperfect CSI. IEEE Trans. Wireless Commun. 15(6), 4455–4469 (2016).
P Ubaidulla, S Aissa, Robust two-way cognitive relaying: precoder designs under interference constraints and imperfect CSI. IEEE Trans. Wireless Commun. 13(5), 2478–2489 (2014).
C Song, H Park, H Lee, I Lee, Robust beamforming designs for nonregenerative multipair two-way relaying systems. IEEE Trans. Wireless Commun. 65(9), 7802–7808 (2016).
M Gastpar, M Vetterli, On the capacity of large Gaussian relay networks. IEEE Trans. Inf. Theory. 51(3), 765–779 (2005).
Y Huang, D Palomar, S Zhang, Lorentz-positive maps and quadratic matrix inequalities with applications to robust MISO transmit beamforming. IEEE Trans. Signal Process. 61(5), 1121–1130 (2013).
Y Rong, Joint source and relay optimization for two-way linear non-regenerative MIMO relay communications. IEEE Trans. Signal Process. 60(12), 6533–6546 (2003).
Q Zhang, Q Li, J Qin, Beamforming design for OSTBC-based AF-MIMO two-way relay networks with simultaneous wireless information and power transfer. IEEE Trans. Veh. Technol. 65(9), 7285–7296 (2016).
X Xie, H Yang, AV Vasilakos, Robust transceiver design based on interference alignment for multi-user multi-cell MIMO networks with channel uncertainty. IEEE Access. 5:, 5121–5134 (2017).
Y Tang, J Xiong, D Ma, X Zhang, Robust artificial noise aided transmit design for MISO wiretap channels with channel uncertainty. IEEE Commun. Lett. 17(1), 2096–2099 (2013).
Q Li, Q Zhang, J Qin, Robust beamforming for cognitive multi-antenna relay networks with bounded channel uncertainties. IEEE Trans. Commun. 62(2), 478–487 (2014).
S Boyd, L Vandenberghe, Convex optimization (Cambridge University Press, Cambridge, 2004).
J Liu, F Gao, Z Qiu, Robust transceiver design for downlink multiuser MIMO AF relay systems. IEEE Trans. Wireless Commun. 14(4), 2218–2231 (2015).
The authors would gratefully acknowledge the grants from the National Natural Science Foundation of China (61371113 and 61401241), Nantong University-Nantong Joint Research Center for Intelligent Information Technology (KFKT2016B04 and KFKT2017B02), the Science and Technology Program of Nantong (GY22017013), and the Brain Korea BK21 plus.
Finally, we would like to thank the Editor for your constructive remarks and careful reading of our paper, which were essential in improving the overall presentation of the paper.
School of Electronics and Information, Nantong University, Nantong, China
Wei Duan
, Xiaojun Zhu
, Li Jin
& Guoan Zhang
Division of Electronic and Information Engineering, Chonbuk National University, Chonju, Korea
Department of ECE, CAIIT, Chonbuk National University, Chonju, Korea
Jeaho Choi
Nantong Research Institute for Advanced Communication Technologies, Nantong, 226019, Jiangsu, China
Xiaojun Zhu
& Wei Wang
Search for Wei Duan in:
Search for Xiaojun Zhu in:
Search for Li Jin in:
Search for Wei Wang in:
Search for Guoan Zhang in:
Search for Jeaho Choi in:
WD, JC, and GZ conceived and designed the study. WD and XZ performed the simulations. WD and GZ wrote the paper. LJ, WW, JC, and GZ reviewed and edited the manuscript. All authors read and approved the final manuscript.
Correspondence to Guoan Zhang.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License(http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
Two-way relay network
Amplify-and-forward
Beamforming design
Imperfect CSI