Transceiver Optimization for Multiuser Multiple-Input Multiple-Output Full-Duplex Amplify-and-Forward Relay Downlink Communications

: This paper considers the transceiver design in a multiuser multiple-input multiple-output (MIMO) full-duplex (FD) amplify-and-forward (AF) relay downlink communication system, where users simultaneously transmit data via an FD relay node. The design incorporates an imperfect loop interference (LI) cancellation which results in a residual LI. Linear precoders are employed at the sources and relay, and minimum mean-squared-error (MMSE) combiners are employed at the destinations to mitigate the effect of the residual LI. The corresponding design problem is highly nonconvex, so a closed-form solution is intractable. Thus, an iterative method is developed to solve this optimization problem. Simulation results are presented which show that the proposed iterative algorithm provides better performance than the corresponding half-duplex (HD) solution in terms of the achievable rate under residual LI.


Introduction
Multiple-input multiple-output (MIMO) relay communication systems have been the subject of considerable research due to their ability to improve both achievable rates and coverage [1][2][3][4][5][6].A relay is used between the source and destination, so the signal is transmitted from the source to the relay and then from the relay to the destination.The relay node can employ either the decode-and-forward (DF) or amplify-and-forward (AF) protocols.AF amplifies the received signal and forwards it to the destination, so it has a lower complexity than the DF protocol [7].
Precoding is a well-known technique for interference mitigation.A joint precoding optimization for a multiuser relay downlink system was investigated in [1].The sum capacity was maximized by using quadratic programming, but multiple antennas were employed only at the relay.The performance can be improved by using multiple antennas at the source and destination.A transceiver design for a multiuser non-regenerative MIMO relay system with multiple antennas at both the source and destination was investigated in [2][3][4][5][6].However, a half-duplex (HD) relay was employed, so transmission from the source to destination requires an extra time slot which limits the potential capacity.
In the past decade, full-duplex (FD) relay systems have attracted attention because data transmission from a source to a destination can be completed in one time slot [7][8][9][10][11].Thus, FD MIMO relaying can increase capacity compared to HD systems [9], but loop interference (LI) is a critical issue because the relay transmits and receives simultaneously.In general, the LI is much larger than the channel noise and so can significantly degrade performance.As a consequence, several LI cancellation techniques have been developed [10].Temporal cancellation methods such as antenna isolation and analog/digital precancellation have been shown to be effective.However, it is impossible to cancel the LI completely, and the residual LI can still be larger than the noise [12].
Estimating channel state information (CSI) is challenging due to the dynamic nature of wireless environments and the complexity of modern communication systems.However, current estimation algorithms have been shown to provide good accuracy, particularly in cellular systems [13].Thus, it can be assumed that CSI is available in the proposed system.In [14], deep neural networks were considered for channel estimation with multiuser precoding in the downlink of a frequency-division duplex massive MIMO system.In [15], an efficient transceiver design was examined for full-duplex (FD) multiuser massive MIMO systems operating in millimeter wave frequencies.An iterative solution was developed for the proposed digital MMSE transceiver beamforming.Secrecy in cell-free massive MIMO systems was proposed in [16] through the joint optimization of uplink power control, downlink beamforming, and duplex mode selection.Artificial noise was incorporated, and an iterative method was used to solve the associated nonconcave-convex approximation.In [17], beamforming was investigated for in-band FD multi-cell multiuser MIMO networks.Hardware impairments, channel uncertainty, and limited channel state information (CSI) were considered in the MMSE beamforming design to improve system performance under global and local CSI.However, none of the approaches in the literature considered the practical limitation of residual loop interference due to imperfect CSI cancellation.This interference should be considered in developing robust and effective wireless communication systems.The impact of residual self-interference (SI) and loop interference (LI) in a full-duplex space shift keying communication system was investigated in [18].Thus, the residual LI due to imperfect interference cancellation is considered in this paper.
Unlike the uplink communications in [19], this paper considers a MIMO FD relaying downlink communication system with multiple source-destination pairs.A transceiver design for the sources and relay, and a linear combiner design at the destinations, are investigated considering the residual LI with the goal of minimizing the MSE of the received signals at the destinations.An iterative algorithm is developed to update the source and relay transceiver matrices and the linear combining matrix at the destinations.
The remainder of this paper is organized as follows.Section 2 presents the system model, and the optimization problem is solved in Section 3. Results are presented in Section 4, which demonstrate the effectiveness of this solution, and some conclusions are given in Section 5.
Notation: Bold uppercase, bold lowercase, and normal letters denote matrices, vectors, and scalars, respectively.vec(•) denotes matrix vectorization, ⊗ denotes the matrix Kronecker product, tr{•} is the trace of a matrix, and I N is the N × N identity matrix.

System Model
Consider a MIMO full-duplex (FD) relay system with K source-destination pairs communicating simultaneously with the aid of a relay.The source-to-destination links are assumed to be inconsequential due to the large-scale fading and the large distances between them.
Figure 1 presents the system model with the kth source-destination pairs equipped with N s k and N d k antennas, respectively.The FD relay has N r and N t antennas to simultaneously receive and transmit signals, respectively.Therefore, communication between the source-destination pairs is accomplished in one time slot compared to an HD system that requires two time slots. where is an independent and identically distributed (i.i.d.) additive white Gaussian noise (AWGN) matrix.As the received signal y r [n] is corrupted by residual LI, from [12], it can be written as where is the residual LI.At time n + 1, the relay amplifies the received signal with a relay transceiver matrix Γ[n + 1] ∈ C N t ×N r and forwards the resulting signal to the destinations.The received signal at the kth destination is then where is the channel matrix between the relay and destination k, and AWGN vector with a zero mean and unit variance.As in [6, 7,12], channel state information (CSI) is assumed to be available at all nodes.In addition, channel variations during the transceiver matrix update interval are assumed to be relatively small [19] and so do not influence the transceiver design [12].Thus, the received signal at the kth destination can be expressed as where As in [18,20], ∆t is modeled as white Gaussian noise which is independent of n r and n d k , and has variance σ 2 t .Consequently, the covariance matrix of n k is A linear combiner M k ∈ C N k ×N b is used at the kth destination to obtain the signal estimate The mean square error (MSE) at the kth user is defined as where Optimization is employed to obtain the linear transceiver matrices Θ k and Γ at the source and relay, and to obtain the linear combiners M k at the destinations.The minimization of the sum mean squared error (SMSE) of the received signals at the destinations can be expressed as min where SMSE = ∑ K k=1 tr(MSE k ), and P s k > 0 and P r > 0 are the power constraints at the kth source and relay, respectively.

Optimization Problem Solution
The problem in ( 9) is highly nonconvex which makes obtaining a global optimal solution intractable, so an iterative algorithm is presented.This solution is based on an alternating optimization that updates Θ k , Γ and M k individually with the other parameters kept fixed to solve the three convex subproblems.
First, given Θ k and Γ, the optimal combiner at the kth destination M k is obtained by solving the unconstrained convex problem for M k , which is independent of the constraints in (9b) and (9c).The optimal solution can then be obtained by taking the derivative of (8) with respect to M k and setting it to zero.Solving ∂ ∂M k tr(MSE k ) = 0 gives where Second, with M k from (10) and given Θ k , Γ can be obtained by solving the following problem.The objective function in (9) can be expressed as Note that since Θ k is known, the constraint in (9c) is eliminated.The optimization problem then becomes min Consider the singular-value decomposition (SVD) of the equivalent source-to-relay and relay-to-destination channels: and The optimal relay precoding matrix is similar to that in [7]: and from ( 13) and ( 14), we have and Note that T , and U g = [U T g1 , . . ., U T g K ] T , which have dimensions d × L 1 and N i × L 2 , respectively.Substituting (13)-( 15) in (11) gives Using the matrix identities where vec(•) concatenates the columns of a matrix into a single vector, the SMSE in (11) can be expressed as a function of a = vec(A): where does not depend on a and The relay transmit power constraint in ( 12) can be written as 24) can be written as The original relay optimization problem is then given by min This is a quadratically constrained quadratic programming (QCQP) problem which is convex and so can be efficiently solved using the interior point method.For example, the CVX toolbox for disciplined convex programming [21] can be employed.Third, the kth source precoding matrix Θ k can be obtained using M k and Γ given above.The corresponding optimization problem can be formulated as the following QCQP problem: where ) can be ignored in the optimization as they do not depend on Θ k .Using the matrix identities, (27) can be written as a function of where and which is independent of Θ k and so can be ignored.Introducing T = bd(T 1 , . . ., T K ), where bd(•) denotes a block-diagonal matrix, and Sk = [S k1 , . . ., S kK ], where S kk = S k and Now, defining and Ē = bd( Ēi1 , Ēi2 , . . ., ĒiK ), where ĒiK = I dN s and Ēij = 0, i ̸ = j, the optimal θ can be obtained by solving the following problem This is a QCQP problem which can be solved using the CVX toolbox for disciplined convex programming.The symbols employed in this work are summarized in Table 1.The proposed iterative algorithm is given in Algorithm 1.The convergence of this algorithm is guaranteed as follows.As the three subproblems are convex, each update of Θ k , Γ, and M k may decrease or at least not increase the value of the objective function.Thus, the iterative algorithm will converge to an optimum solution.
Algorithm 1: Iterative Design of Θ k , Γ, and M k .
1: Initialize the algorithm with Θ (0) = P s L I L and I N r , and set n = 0.

Numerical Results
This section presents the evaluation of the performance of the proposed multiuser MIMO FD relay downlink algorithm using numerical simulation.The CVX toolbox was employed, which is a Matlab-based software package for solving convex optimization problems [21].For simplicity, we considered a system with two sources and two destinations.The extension to more than two source-destination pairs is straightforward.We assumed that the sources and destinations were each equipped with two antennas, and the relay was equipped with four receive antennas and four transmit antennas.The simulation parameters are given in Table 2.
As in the related literature, flat-fading MIMO channels were considered.It was assumed that the elements of H k , G k , and R LI were i.i.d.complex Gaussian random variables with a zero mean and unit variance.Moreover, all noise terms were i.i.d.complex, circularly symmetric Gaussian random variables with a zero mean and unit variance.According to [9], the residual LI can be 0 dB to 20 dB larger than the channel noise.Therefore, LI levels of 0 dB, 5 dB and 10 dB were employed here.In all cases, the results are given for an average of 1000 independent channel realizations.Note that the HD system mentioned in this paper is the same as the proposed FD system except that the residual LI term is zero and the achievable rate of the HD system is half of the FD system since two time slots are required for transmission between a source and destination.Thus, the same optimization procedure is employed.
Figure 2 presents the SMSE for the HD and FD systems with different levels of residual LI.The signal-to-noise ration (SNR) between the sources and relay varies from 0 dB to 30 dB, and the SNR between the relay and destinations is fixed at 30 dB.These results indicate that the HD system provides the best SMSE as the FD performance is degraded as the LI increases.The achievable rate for the system in ( 5) can be obtained using an approach similar to that in [8], and can be written as where Figure 3 gives the achievable rate for the HD and FD systems.In this figure, the SNR between the sources and relay varies from 0 dB to 30 dB while the SNR between the relay and destinations is fixed at 30 dB.The HD system corresponds to the case when the residual LI is zero and two time slots are required for transmission from sources to destinations.Therefore, the FD achievable rate is twice the HD achievable rate if the LI is cancelled completely.Further, the achievable rate with the FD system is higher when the residual LI is 0 dB to 10 dB.The achievable rate of the FD system is degraded as the residual LI increases in the high-SNR region since the residual LI is greater than the multiuser interference.In the low-SNR region, the multiuser interference dominates the residual LI, so the effect of the residual LI is negligible.Figures 4 and 5 present the SMSE and achievable rate with a fixed SNR of 30 dB between the sources and relay and an SNR between the relay and destinations from 0 dB to 30 dB.The SMSE in Figure 4 is better than in Figure 2 in the low-SNR region because a higher transmit power results in more residual LI. Figure 5 shows that the achievable rate of the FD system is always higher than that of the HD system for residual LI levels up to 10 dB.These results confirm that the proposed approach is effective in addressing the effects of the residual loop interference due to imperfect LI cancellation in full-duplex systems.Further, this solution is superior to conventional half-duplex communications.Table 3 gives the percentage improvement in achievable rate with the proposed system over the conventional half-duplex system for different levels of residual LI and SNRs.These results demonstrate the robustness of the proposed design as the improvements exceed 75% when the residual LI is low.Further, they confirm the effectiveness in optimizing performance in practical scenarios.Achievable Rate (bps/Hz)

Conclusion
This paper considers a transceiver design for the downlink of a multiuser non-regenerative MIMO full-duplex (FD) relay system with residual loop interference (LI).The source and relay transceiver matrices and destination combining matrix are optimized to minimize the sum mean squared error (SMSE).The original nonconvex problem is converted to three convex subproblems, and an iterative algorithm is developed to optimize the three matrices.Results are presented which demonstrate that the proposed iterative method outperforms the corresponding half-duplex (HD) relay system in terms of the SMSE and achievable rate.The effectiveness of the proposed solution was illustrated via simulation.Future research can consider interference cancellation techniques, integrating machine learning for dynamic system optimization, and next-generation wireless technologies.Internet of things (IoT) and mobile communication applications can also be examined considering the adaptability and scalability of the proposed approach and to further confirm its impact.

Figure 1 .
Figure 1.The multiple-user MIMO full-duplex (FD) relay system.Let s k [n] represent the length d signal vector at time n for the kth source.A linear transceiver matrix Θ k [n] is applied to s k [n] before transmission.The received signal at the relay can then be expressed as

Table 3 .
Percentage improvement in achievable rate.