Precoding Design and Power Allocation in Two-User MU-MIMO Wireless Ad Hoc Networks

: In this paper, we consider the precoding design and power allocation problem for multi-user multiple-input multiple-output (MU-MIMO) wireless ad hoc networks. In the first timeslot, the source node (SN) transmits energy and information to a relay node (RN) simultaneously within the simultaneous wireless information and power transfer (SWIPT) framework. Then, in the second timeslot, based on the decoder and the forwarding (DF) protocol, after reassembling the received signal and its own signal, the RN forwards the information to the main user (U1) and simultaneously sends its own information to the secondary user (U2). In this paper, when the transmission rate of the U1 is restricted, the precoding, beamforming, and power splitting (PS) transmission ratio are jointly considered to maximize the transmission rate of U2. To maximize the system rate, we design an optimal beamforming matrix and solve the optimization problem by semi-definite relaxation (SDR), considering the high complexity of implementing the optimal solution. Two sub-optimal precoding programs are also discussed: singular value decomposition and block diagonalization. Finally, the performance of the optimization and sub-optimization schemes are compared using a simulation.


Introduction
Recently, energy harvesting (EH) has received widespread attention as a promising technology capable of overcoming the problem of energy shortages in wireless communication networks [1].Many energy harvesting models have been proposed in [2,3].In particular, [4] studies a time allocation scheme to strike a balance between wireless power transfer and information transfer.These models demonstrate that wireless communication systems can effectively support EH [5,6].Since radio frequency (RF) signals can transmit information and energy simultaneously, a joint investigation of SWIPT has recently drawn significant attention [7], and various models for different SWIPT systems have been considered to characterize the performance trade-off between wireless energy transfer (WET) and wireless information transfer (WIT) in the same transmitted signal [8][9][10].Furthermore, one key issue in implementing SWIPT is the practical limitation that existing energy harvesting circuits cannot be used to decode and harvest the RF signals concurrently [11], to solve this problem, PS designs have been proposed for SWIPT, with PS, the received signal is split into two parts, with one part used for ID and the other part for EH [12].
Although EH can solve the problem of energy shortages in wireless networks, when numerous devices are connected to the network, the available spectrum will be seriously inadequate, Furthermore, signal coverage is also a problem in remote areas.These problems are expected to be solved using MIMO technology [13], which can effectively increase data throughput and extend system transmission distance without requiring additional bandwidth or increasing the total transmission power consumption [14].To significantly improve spectrum utilization in MIMO networks, cognitive radio has been widely considered by previous studies [15,16], in which [15] considered a linear precoder design for an underlay cognitive radio MIMO broadcast channel and [16] considered the power allocation for spectrum sharing in cognitive radio networks.
The two basic features of MIMO wireless propagation channels are broadcasting and multipleaccess [17].Therefore, signals transmitted on a certain link can be detected by other neighboring nodes in the wireless network.Those nodes can form an ad hoc network and serve as relays that help forward the message to the receiving terminal.In addition, although MU-MIMO has been used in ad hoc networks to improve system performance (e.g., in [18] for capacity enhancement and in [19] for stream control), using a precoding design and power allocation in MU-MIMO ad hoc networks to improve throughput has not been considered in previous works.Moreover, using a SWIPT MIMO strategy with relays can effectively improve the transmission rate, and the communication quality is determined by the RF energy receiver structure and the precoding performance.However, neither the singular value decomposition (SVD) algorithm nor the block diagonalization (BD) algorithms adopted in the precoding design of previous SWIPT MIMO cooperative systems have yet been considered.
Motivated by the aforementioned problem, to further improve the throughput of a SWIPT MIMO cooperative system, this paper studies a two-user MIMO system in a wireless ad hoc network (WAHN).First, by using a PS receiver under the DF protocol, relay transmission beamforming, and power splitting are analyzed.Then, under the premise of satisfying the quality of service of the main users, power allocation and an optimal precoding design are conducted jointly to maximize the transmission rate of the secondary users.Finally, to compare the performance, two sub-optimal precoding methods are designed.The optimal solution is found by applying semi-definite relaxation (SDR) techniques; it is proven that a rank-one optimal solution always exists.Furthermore, Karush-Kuhn-Tucker (KKT) conditions are also applied to solve the problem.
The remainder of this paper is organized as follows: Section 2 introduces the system model and presents the problem formulation.In Section 3, the SDR algorithm is adopted to derive the optimal solution, and we prove that an optimal solution is feasible.Based on the SVD and BD algorithms, Section 4 presents two sub-optimal solutions for the first and second timeslots, respectively.Section 5 presents numerical results to compare the performance of various approaches, and Section 6 concludes the paper.

System Model and Problem Formulation
As shown in Figure 1, this paper considers a two-user MIMO wireless ad hoc network that consists of one source node, one relay node, and two user terminal nodes (denoted by U1 and U2, respectively).We assume that the source node (SN) and the relay node (RN) are both equipped with M ≥ 1 transmitting antennas and that the SN has a stable energy supply, whereas the RN must replenish its energy from the received signals broadcast by the SN in the first timeslot.This energy is then stored and used to forward the information from the SN to the main user, U1, and to transmit its own information to the secondary user, U2, in the second timeslot.The users' terminals are both equipped with N ≥ 1 receiving antennas.
In this paper, the SN is considered to improve the efficiency of information and energy transmission through a precoding design, and the precoding design at the RN is used to ensure that the receiving user can efficiently decode the received signal.In the first phase (timeslot 1), assume that the EH and information decoding (ID) receivers are co-located in the RN, the channel state information between the SN and the RN is perfectly known, and the channel matrix is denoted as . The energy receiver of the relay node does not need to collect energy by switching the signal from the RF band to the baseband.Furthermore, based on the principle of energy conservation, P denotes the total RF-band power (energy normalized by the baseband symbol period) harvested from the receiver's receiving antennas.It is worth noting that P is proportional to the received baseband signal: where V denotes the energy conversion efficiency.For convenience, in this paper, unless stated otherwise, it is assumed that 1 V = .
Then, the received signals at the RN can be expressed as follows: where denotes the channel coefficients from the SN to the RN, and is the signal power splitting ratio for the ID phase, and the ratio of the EH phase can be expressed as 1 r -.Thus, the transmission rate can be expressed as follows: The second phase requires harvesting more energy to ensure continued information transmission.If the threshold of the rate R g ³ is satisfied, the remaining power will all be used to harvest energy.The energy harvested by the RN can be expressed as: ( ) where L = -L , and g is the minimum transmission rate that the SN requires.
In the first timeslot, under the premise of satisfying the SN rate requirement, maximizing the energy received by the RN can allow it to transmit more information in the second timeslot, thus achieving the goal of maximizing system throughput.Then, the optimal problem in the first phase of the PS scheme based on a cooperative network can be formulated as follows: In the first phase (timeslot 1), assume that the EH and information decoding (ID) receivers are co-located in the RN, the channel state information between the SN and the RN is perfectly known, and the channel matrix is denoted as h ∈ C M×M .The energy receiver of the relay node does not need to collect energy by switching the signal from the RF band to the baseband.Furthermore, based on the principle of energy conservation, P denotes the total RF-band power (energy normalized by the baseband symbol period) harvested from the receiver's receiving antennas.It is worth noting that P is proportional to the received baseband signal: where ς denotes the energy conversion efficiency.For convenience, in this paper, unless stated otherwise, it is assumed that ς = 1.Then, the received signals at the RN can be expressed as follows: where h ∈ C M×1 denotes the channel coefficients from the SN to the RN, and z ∈ C M×1 denotes the relay noise vector.It is assumed that z are independent overall and that z ∼ CN (0, I).The covariance matrix of x is denoted by S S = E xx H , the power constraint in the source terminal is denoted by E[ x 2 ] = tr(S S ) ≤ P S , ρ(0 < ρ < 1) is the signal power splitting ratio for the ID phase, and the ratio of the EH phase can be expressed as 1 − ρ.Thus, the transmission rate can be expressed as follows: The second phase requires harvesting more energy to ensure continued information transmission.If the threshold of the rate R ≥ γ is satisfied, the remaining power will all be used to harvest energy.The energy harvested by the RN can be expressed as: where and γ is the minimum transmission rate that the SN requires.
In the first timeslot, under the premise of satisfying the SN rate requirement, maximizing the energy received by the RN can allow it to transmit more information in the second timeslot, thus achieving the goal of maximizing system throughput.Then, the optimal problem in the first phase of the PS scheme based on a cooperative network can be formulated as follows: It is assumed that the relay is capable of decoding the source signal and forwarding it appropriately.In the second phase, the relay simultaneously transmits signals to U1 and U2.The relay node will decode the signal and reconstruct the message.The RN transmission signal consists of two parts: the signal revived from the SN is S S ∈ C N×1 , and the RN's independent signal is S R ∈ C N×1 .Then, the signal vector received at U1 is: and the signal vector received at U2 is: where h 1 ∈ C N×M and h 2 ∈ C N×M denote the channel coefficients between the RN and U1 and U2, respectively.The RN uses the pre-decoding vector W S to forward the signal from the SN and the pre-decoding vector W R to transmit its own signal.The additional Gaussian noise at the destination terminal is ). Accordingly, the signal-to-interference plus the noise ratio (SINR) of U1 can be expressed as follows: and the SINR of U2 can be expressed as: In the next section, this paper will investigate the optimal precoding vectors W S and W R to eliminate the interference in the signal between two users and maximize the transmission rate.
Since the RN forwards the SN's signal based on the DF protocol, as was concluded in [14], the most effective condition is achieved when the transmission rate from the SN to the RN in the first timeslot is equal to the transmission rate of the RN to U1 in the second timeslot.Thus, in the second timeslot, the optimal target is to maximize the rate of U2 while still ensuring the required minimum rate of U1.Consequently, the optimization problem can be formulated as follows: where ψ denotes the required minimum rate of U1.

Feasibility Proof and Optimal Solution
The optimization problem in the first timeslot is similar to that in traditional SWIPT and has been widely shown to be feasible [20]; therefore, this section focuses on the feasibility of the optimization problem P2 in the second timeslot.

A Feasibility Analysis of Problem P2
To prove the feasibility of problem P2, first, the SINR of U2 is set to zero (SI NR 2 = 0): the relay forwards only the SN's information and does not transmit its own signal.In this case, problem P2 is transformed into a typical single-user MIMO system relay node problem, which has been solved in many studies [21].Based on the DF protocol, the maximum transmission rate can be achieved when the SINR is met in the first timeslot (SI NR = SI NR 1 ).Thus, P2 can be reformulated as follows: where W S is the only variable in Equation (11).By setting W = W S W S H , H 1 = h 1 h 1 H and applying semi-definite planning (SDP) in P2-A, P2-A can be expressed as: If for P2-A-SDP, there exists an optimal solution W * , then the optimal solution W * S for problem P2-A can be obtained by eigenvalue decomposition (EVD).Since problem P2-A-SDP is a convex problem, it can be solved via existing software (e.g., Matlab CVX developed by MathWorks, Natick, MA, USA [22,23]).
Since problem P2-A is feasible, and problem P2 is similar to problem P2-A, P2 is also a feasible problem.
Next, the optimization algorithms for problems P1 and P2 are described.

The Optimization Algorithm for Problem P1
The optimization goal of problem P1 is to maximize a concave function, and the function in P1's constraint is a convex function; therefore, P1 is a typical convex optimization problem [11].In the first phase, P1 is non-convex in its initial form because both the object function and the subject function involve coupled S S and Λ ρ .This problem can be solved via a two-step procedure.First, fix ρ = ρ through a simple one-dimension search over 0 < ρ < 1.Then, P1 can be re-described as follows: max After transforming P1 into a convex problem with only one unknown dimension, it can be solved by existing solvers, such as Matlab CVX [22].An one-dimensional search is performed on ρ, and the optimal precoding matrix S S is obtained using the CVX toolbox.At that point, the ρ corresponding to the optimal precoding matrix S S is the optimal power allocation factor.

The Optimization Algorithm for Problem P2
Applying beamforming (BF) to the MU-MIMO transmit signals can effectively eliminate inter-signal interference and improve system throughput.In this section, the optimal precoding matrices W S and W R in the optimization problem P2 are solved by BF.
First, H = h n h H n , n = 0, 1, 2 denotes the channel matrixes of h, h 1 , h 2 , S = W s W s H and H by dropping the rank-one constraint.Then, the SDR of problem (P2) can be expressed as follows: P2 − SDR : max If the optimal solution of P2-SDR results in an S * and Q * that satisfy the conditions rank(S * ) = 1 and rank(Q * ) = 1, respectively, then the optimal W * s and W * R can be obtained by performing EVD on S * and Q * , respectively.Furthermore, if rank(S * ) > 1 or rank(Q * ) > 1, then the optimal values S * and Q * of P2-SDR can be used only as the upper limit of the optimal solution of problem P2, and W * S and W * R must be obtained by other algorithms, such as rank-approximation.
Note that problem P2 is a non-convex problem, not only because the SINR expression contains two matrix vectors, but also because the constraints are quadratic.Therefore, because the numerator and denominator of the objective function in problem P2-SDR contain S and Q, respectively, P2-SDR is also a non-convex problem.Thus, the Charnes-Cooper transformation can be used to reformulate P2-SDR as an equivalent convex problem.In addition, both have the same optimal value (for a proof, see [5]).After introducing an equivalent substitution factor t in P2-SRD, the problem can be expressed as: Obviously, P2-SDR-Equ is convex and satisfies Slater's condition-the duality gap between P2 and P2-SDR-Equ is zero.Therefore, P2-SDR-Equ can be efficiently solved by existing solvers, such as Matlab CVX developed by MathWorks Natick, MA, USA [23].Since there are two design variables, based on the result in [2], an optimal solution (S * , Q * ) exists for P2-SDR-Equ such that rank 2 (S * ) + rank 2 (Q * ) ≤ 2K + 2. Under our system model, there is only one secondary user, U2, where k = 1.Moreover, the rank of the precoding matrix must be an integer.Consequently, there is only one solution in which rank(S * ) = 1 , rank(Q * ) = 1.Thus far, the relaxation rank-one condition is satisfied.

Two Suboptimal Solutions
Since solving the optimization algorithm is complex and has high computational complexity, to facilitate practical applications, this section proposes two sub-optimal algorithms for slots 1 and 2, respectively.

Sub-Optimal Algorithm with SVD in the First Timeslot
First, we consider the first phase's sub-optimal algorithm with SVD.In slot 1, the SN transmits the signal x to the RN with the precoding matrix F ∈ C M×M .When the channel state information is perfectly known, the received signals at the RN can be represented by: where P s represents the source transmission power and h ∈ C M×M denotes the channel coefficients between the SN and the RN.Since the SN follows the DF protocol, the received signal must be decoded at the receiving terminal; thus, SVD is used to code the transmit signal.When the channel matrix h = UDW H is decomposed, the precoding matrix becomes a unitary matrix F = W. Based on the conjugate transpose matrix, we can obtain the equalization matrix G = U H at the receiving terminal.Through the decoding process, the signal received at the terminal can be expressed as follows: Since the received signal is in the form of non-interference diagonalization, after signal processing, the decoded signal can be used for transmission during the second phase.The additional Gaussian noise at the RN is z ∈ CN (0, σ 2 R ).Since the RN is a battery-free terminal, it needs to scavenge energy from the received signal transmitted by the SN.Furthermore, using the harvested energy, the RN simultaneously decodes and forwards both the received signal and its own signal to U1 and U2, respectively.By applying the power splitting scheme, the signals at the RN are split into two parts: one is used for energy harvesting and the other to receive information.We denote ρ(0 < ρ < 1) as the power splitting ratio, where the ratio ρ of the available power is used for receiving information, and the remaining portion (1 − ρ) is used for energy harvesting.The harvested energy P (in joules) at the RN during the first timeslot is: where ζ (0 < ζ < 1) denotes the energy conversion efficiency at the RN.Then, the remaining power is split at the information receiver.As a result, the signal split at the information receiver is expressed as follows: The SINR between the SN and the RN is: In the first timeslot, the optimization problem is to determine the power allocation factor ρ, which can be expressed as follows: The objective of P1-SVD is to maximize a concave function, and all the functions in the constraint condition are linear functions that include the unknown variable ρ.Therefore, P1-SVD is a convex optimization problem, and its optimal solution can be obtained by the KKT condition.
as to solve the sub-optimal algorithm with BD in the second timeslot.Furthermore, we compared the performances between three different precoding and power allocation strategies.
First, consider that both the SN and the RN are equipped with four antennas (N = 4) and that both U1 and U2 are equipped with two antennas (M = 2).
Figure 2 shows a comparison of the system rates achieved by different pre-decoding schemes.
the RN.Two terminals are separated at a normalized distance.We adopt the simulation parameters that have been widely used in many classical models of MIMO systems [17,20]; to better justify the comparison results, we set the maximum power of SN to 1 S P = Watt (W) or 30 dBm, , and Note that the optimal solutions obtained through the BF algorithm in the two timeslots use the Matlab CVX toolbox to solve Equations ( 13) and (15), as well as to solve the sub-optimal algorithm with BD in the second timeslot.Furthermore, we compared the performances between three different precoding and power allocation strategies.First, consider that both the SN and the RN are equipped with four antennas ( 4 N = ) and that both U1 and U2 are equipped with two antennas ( 2 M = ).
Figure 2 shows a comparison of the system rates achieved by different pre-decoding schemes.Based on Figure 2, we can draw the following conclusions: (1) The higher the rate requirement is for the main user, U1, the smaller the achievable rate is for the secondary user, U2.This is because for U1 to achieve a high rate, the RN can harvest less energy in the first timeslot; consequently, less energy is allocated to U2 in the second timeslot.(2) The precoding design that uses the BF algorithm at both the first and second timeslots achieves the optimal performance; however, considering the high complexity of the optimization algorithm, we proposed two sub-optimal algorithms: SVD and BD.(3) When there are fewer receiving antennas than transmitting antennas, the performance of the sub-optimal precoding design with the BD algorithm achieves a performance close to that of the optimal precoding design with the BF algorithm.This occurs because the interference signals can effectively map to the null space of the useful signal with the BD algorithm; thus, it achieves an excellent performance in cancelling interference effects.In contrast, the performance with the SVD algorithm significantly impacts the performance of the sub-channel; its effect is relatively poor, but the receiver is simple to design.Based on Figure 2, we can draw the following conclusions: (1) The higher the rate requirement is for the main user, U1, the smaller the achievable rate is for the secondary user, U2.This is because for U1 to achieve a high rate, the RN can harvest less energy in the first timeslot; consequently, less energy is allocated to U2 in the second timeslot.(2) The precoding design that uses the BF algorithm at both the first and second timeslots achieves the optimal performance; however, considering the high complexity of the optimization algorithm, we proposed two sub-optimal algorithms: SVD and BD.(3) When there are fewer receiving antennas than transmitting antennas, the performance of the sub-optimal precoding design with the BD algorithm achieves a performance close to that of the optimal precoding design with the BF algorithm.This occurs because the interference signals can effectively map to the null space of the useful signal with the BD algorithm; thus, it achieves an excellent performance in cancelling interference effects.In contrast, the performance with the SVD algorithm significantly impacts the performance of the sub-channel; its effect is relatively poor, but the receiver is simple to design.
Finally, the condition that all nodes in the system are equipped with four antennas (N = M = 4) has been considered.Figure 3 shows a comparison of the system rate for different pre-decoding schemes.
Finally, the condition that all nodes in the system are equipped with four antennas ( has been considered.Figure 3 shows a comparison of the system rate for different pre-decoding schemes.From Figure 3, we can conclude the following: (1) Increasing the number of antennas in a user's receiving terminal effectively increases the system's capacity.The optimized beamforming scheme is not affected by the performance of the separated sub-channels nor by the number of receiver antennas; therefore, it achieves the maximum rate, and the optimization effect is remarkable.(2) When the number of receiver antennas is larger than the number of transmitter antennas, the channel dimensions also increase, and BD can no longer effectively eliminate the signal interference between users, therefore, its performance is the worst among the three tested precoding designs.

Conclusions
This paper presented optimal and sub-optimal precoding designs and power allocation for a multiuser MIMO ad hoc network.The system rate is maximized based on a specified U1 rate requirement and the power-harvesting constraints of the relay nodes.First, an optimized precoding scheme based on a BF algorithm is proposed, and the optimal solution is obtained using the SDR strategy.Then, the suboptimal precoding schemes for two timeslots are given by the SVD and BD algorithms, respectively, to reduce the computational complexity caused by implementing the optimized precoding scheme.Finally, the results of simulations show that the precoding scheme can significantly improve system speed and has practical significance.

Figure 2 .
Figure 2. When N = 4, M = 2, the second user's rates with different precoding designs based on different rate requirement for the main user.

Figure 2 .
Figure 2. When N = 4, M = 2, the second user's rates with different precoding designs based on different rate requirement for the main user.

Figure 3 .
Figure 3.When N = 4, M = 4, the second user's rates with different precoding designs based on different rate requirements for the main user.