I Introduction
Recently, wireless backscatter has been proposed as a new promising technology to sustain wireless communications for Internet of Things (IoT) [21]. It is featured with extremely low power consumption by transmitting information in the passive mode via the modulation and reflection of the ambient RF signals [9]. In contrast, the conventional RF radios operate in the active mode, which transmits information by selfgenerated carrier signals relying on powerconsuming active components such as power amplifier and oscillator. The architectural differences between the active and passive radios lead to complement transmission capabilities and power demands in two modes, e.g. [5] and [4]. In particular, the active radios may transmit with a higher data rate and better reliability, however requiring more power consumption. While the passive radios generally achieve a lower data rate with significantly reduced power consumption. This implies that a hybrid wireless system with both active and passive radios can flexible schedule the data transmissions in different modes according to the channel conditions, energy status, and traffic demands, and thus achieve a higher network performance, especially for energy constrained IoT networks.
Due to the extremely low power consumption of the passive radios, it becomes promising to use the passive relays to assist the active RF communications, especially for wirelessly powered IoT networks with energy harvesting constraints, e.g., [13, 2, 19]. However, most of the existing works have considered a simple passive relay model in which one or more backscatter radios are employed to assist the active RF communications. In this paper, we propose a novel hybrid relay communications model in which both the passive and active relays are employed simultaneously to assist the active RF communications. Based on the dualmode radio architecture [4], each relay can switch between the active and passive modes independently to maximize the overall throughput, according to its channel conditions and energy status. With more relays in the passive mode, less active relays can be used for amplifying and forwarding the source signal. On the other hand, the active relays’ channels can be enhanced by the passive relays. It is obvious that the coupling between relays in two modes complicates the optimization of individual relays’ mode selection. As the conventional relay models focus on either the active or passive relays, the current relay strategies are not applicable to hybrid relay communications, which motivates our novel algorithm design in this paper.
Specifically, we focus on a twohop hybrid relay communication model in this paper. In the first hop, the multiantenna transmitter beamforms information to the relays and the receiver. The set of passive relays instantly backscatter the RF signals to enhance signal reception at both of the active relays and the receiver. In the second hop, the set of active relays jointly beamform the received signals to the receiver. Meanwhile, the passive relays can adapt their reflection coefficients to enhance the active relays’ forwarding channels. We aim to maximize the overall network throughput by jointly optimizing the transmit beamforming, the relays’ mode selection, and their operating parameters. It is clear that the throughput maximization problem is combinatorial and difficult to solve optimally. To overcome this difficulty, we propose a twostep solution to optimize the relay strategy. A preliminary study on hybrid relay communications has been reported in our previous work [24]. In this paper, we provide different evaluations of the relay performance, which motivate us to design a set of iterative algorithms that can be more effective than that in [24]
. Specifically, with a fixed relay mode, we firstly find two feasible lower bounds on the signaltonoise ratio (SNR) at the receiver under different channel conditions, based on which we devise a set of performance metrics to evaluate individual relay’s performance gain. Then, we provide an iterative procedure to update the relays’ mode selection that can improve the overall relay performance in each iteration. Simulation results verify that the proposed hybrid relay strategy can significantly improve the throughput performance compared to the conventional relay strategy with all active relays.
To be specific, our main contributions in this paper are summarized as follows:

A novel hybrid relay communications model: Different from the conventional relay communications, we allow multiple energy harvesting relays in both the active and passive modes to collaborate in relay communications. The active relays follow the amplifyandforward protocol while the passive relays can adapt their reflection coefficients to enhance the relay channels.

Lower bounds on relay performance: We propose two lower bounds to evaluate the receiver’s SNR under different channel conditions. The SNR evaluation serves as a performance metric for relay mode selection. Each lower bound requires to solve an optimization problem involving the transmit beamforming, the active relays’ power control, and the passive relays’ phase control.

Heuristic algorithms for mode selection: To bypass the complexity in SNR evaluation, we also propose a set of heuristic algorithms based on simple approximations of the SNR performance to evaluate each passive relay’s performance gain. Our simulation results verify that the SNRbased mode selection can achieve the optimal throughput performance, while the heuristic algorithms also perform well with significantly reduced complexity.
The rest of this paper is organized as follows. Section II summarizes recent works related to hybrid relay communications. Section III describes the system model for hybrid relay communications. Correspondingly, Section IV formulates the throughput maximization problem. In Section V, we devise iterative algorithms to search for the optimal passive relays and the operating parameters of the active relays. Simulation and conclusions are drawn in Sections VI and VII, respectively.
Ii Related Works
The concept of hybrid backscatter communications has been proposed in [7] to increase the transmission range of wirelesspowered communications networks. However, it restricts the user devices to switch between two typical configurations of backscatter communications, i.e., the bistatic and ambient configurations, e.g., [21] and [4]. The authors in [7] verify that an optimized time allocation between two backscatter configurations can achieve improved throughput performance and increased transmission range. Different from [7], in this work we allow each user device to operate in either the active or the passive mode. By optimizing the users’ radio modes, we envision to achieve enhanced network performance in terms of both energy consumption and throughput in a hybrid radio network, which have been corroborated by both analytical and simulation results in the literature. The authors in [5] allow the cognitive radios to switch among energy harvesting, backscatter and active communications. To maximize the performance gain, the time transitions are optimized under different channel conditions according to the presence of primary users. Considering pointtopoint communications between a pair of IoT devices, the authors in [11] present a theoretical analysis on the link capacities when the transmitter follows predesigned protocols to switch between the passive and active modes. The authors in [8, 17, 22] allow the transmitter to switch its radio mode during data frames and further optimize its time transitions to maximize the overall throughput performance. The authors in [26]
propose to integrate the passive backscatter communications with the cognitive radios, and verify a significant performance gain in terms of the coverage probability and ergodic capacity.
Besides mode switching of a single device, the user devices in different radio modes can further collaborate in data transmissions, exploiting both the radio diversity and user cooperation gains. The user cooperation can be envisioned by allowing the passive radios to serve as wireless relays for the active RF communications, e.g., [4] and [2]. Compared to the conventional relay communications, the backscatteraided passive relay communications can be more energy and spectrumefficient due to the extremely low power consumption of passive radios. The authors in [27] show that the backscattered signals can be used to enhance the signal reception performance of the active RF communications. Different from [27], the authors in [16] and [18] employ the active radios to assist information transmission of the passive radios. The rate maximization problem is formulated by jointly optimizing the time scheduling, power allocation, and energy beamforming strategies. The relay needs to decode the backscattered information in the first phase and then forward it to the receiver by active RF communications in the second phase. This implies that the relay’s operation requires significant power consumption in receiving and forwarding, which may prevent it from joining the cooperative transmissions. In [3], both the passive and active radios can assist each other in a twohop relay protocol. The source node firstly transmits information to the relay node via backscatter communications in the first hop. Then, the relay node will decode and forward the information by active RF communications, meanwhile a set of passive radios are selected and optimized to assist the active relay’s RF communications in the second hop. A similar twohop relay model is also studied in [12], where the relay firstly decodes the information and then uses backscatter communications to forward the information to the receiver.
The passive relay communications are also studied for multiuser wirelessly powered communications. In our previous work [2], we allow each passive relay to assist multiple active radios. Considering the couplings among different users, heuristic algorithms are devised to maximize the overall relay performance by jointly optimizing the transmit beamforming, mode switching, and the passive relays’ reflection coefficients, subject to the relays’ energy harvesting constraints. In [25], the couplings among multiple energy harvesting relays are further characterized in a distributed game model, in which each passive relay aims to optimize its reflection coefficients for the other receivers to maximize its utility function.
Iii System model
We consider a downlink communication system with a group of singleantenna user devices coordinated by a multiantenna hybrid access point (HAP), which constantly distributes data stream to individual users. For example, we can envision the HAP as the centralized controller in a wireless sensor network for control information distribution to spatially distributed wireless sensors, which can be lowpower and longlasting IoT devices. The HAP’s information transmissions to different user devices can be scheduled in orthogonal channels, e.g., via a time division multiple access (TDMA) protocol [2]. Without loss of generality, we focus on the simplest case with only one receiver. The user devices can serve as the wireless relays for each other via devicetodevice (D2D) communications. In particular, the data transmission from the HAP to the receiver can be assisted by a set of relays following the amplifyandforward (AF) protocol. The set of relays is denoted by . The HAP has a constant power supply, while the relays are wirelessly powered by RF signals emitted from the HAP. Via signal beamforming, the HAP can control the information rate and power transfer to the relays following the powersplitting (PS) protocol [20]. Our previous work in [14] reveals that the relay optimization with the typical timeswitching protocol has a similar structure with the PS protocol. Hence, we focus on the PS protocol in this work. Our goal is to maximize throughput from the HAP to the receiver by optimizing the HAP’s transmit beamforming and the relay strategy. Assuming that the HAP has antennas, let and denote the complex channels from the HAP to the receiver and from the HAP to the th relay, respectively. Let
denote the complex channels from the relays to the receiver. We assume that all the channels are block fading and can be estimated in a training period before data transmissions
[25].Iiia TwoHop Hybrid Relaying Scheme
The relayassisted information transmission follows a twohop halfduplex protocol. As shown in Fig. 1, the information transmission is divided into two phases, i.e., the relay receiving and forwarding phases, corresponding to the information transmission in two hops. Due to a short distance between transceivers in a dense D2D network, the direct links between the HAP and the receiver can exist in both hops and contribute significantly to the overall throughput. Moreover, leveraging the direct links, we allow the HAP to beamform the same information twice in two hops. This can significantly improve the data rate or transmission reliability, which may be more crucial for many applications such as industrial control process. Specifically, let denote the HAP’s signal beamforming strategies in two phases.
In the first hop, the HAP beamforms the information with a fixed transmit power
and the beamforming vector
. Conventionally, the beamforming information can be received by both the relays and the receiver directly, as shown in Fig. 1(a). Hence, the HAP’s beamforming design has to balance the transmission performance to the relays and to the receiver. A higher rate on the direct link potentially degrades the signal quality at the relays and reduces the data rate of relays’ transmission. Different from the conventional relay communications in Fig. 1(a), where all relays are operating in the AF protocol, in this paper we assume that each relay has a dualmode radio structure that can switch between the passive and active modes, similar to that in [4] and [25]. This arises the novel hybrid relay communications model. As illustrated in Fig. 1(b), when the HAP beamforms the information signal to the relays, the relay can turn into the passive mode and backscatter the RF signals from the HAP directly to the receiver. By setting a proper load impedance and thus changing the antenna’s reflection coefficient [21], the passive relay can backscatter a part of the incident RF signals, while the other part is harvested as the power to sustain its operations. Moreover, the backscattered signals from the passive relays can be coherently added with the active relays’ signals to enhance the signal strength at the receiver [27].The HAP’s beamforming in the first hop is also used for wireless power transfer to the relays. We consider a PS protocol for the energy harvesting relays, i.e., a part of the RF signal at the relays is harvested as power while the other part is received as information signal. Specifically, we allow each active relay to set the different PS ratio to match the HAP’s beamforming strategy and its energy demand. In the second hop, the active relays amplify and forward the received signals to the receiver. Meanwhile, the HAP also beamforms the same information symbol directly to the receiver with a new beamforming vector . Hence, the received signals at the receiver are a mixture of the signals forwarded by the relays and the direct beamforming from the HAP. With maximal ratio combining (MRC) at the receiver [10], the received signals in both hops can be combined together to enhance the data rate and reliability in transmission. Note that the passive relays in the second hop can optimize their reflection coefficients to enhance individual relays’ forwarding channels, as well as the direct channel from the HAP to the receiver.
IiiB Channel Enhancement via Passive Relays
The optimal selection of each relay’s radio mode is complicated by the relays’ couplings in transmission capabilities. The passive relays can enhance the active relays’ channel conditions, while more active relays can achieve user cooperation gain via network beamforming [6]. In the extreme case, the passive mode becomes the only afordable choice when the relay have low power supply, while the active mode can be preferred to provide more reliable transmissions if the relay has good channel conditions and sufficient power supply.
Let
denote the binary variable indicating the radio mode of the relay
for . Then the set of relays in Fig. 1(b) will be split into two subsets, i.e., and , denoting the sets of active and passive relays, respectively. Let and for denote the enhanced channels from the HAP to the receiver and to the active relays, respectively. Let denote the signal transmitted from the HAP with a constant transmit power . The received signal at the passive relay is given by , where is the normalized noise signal. The complex reflection coefficient of passive relay is given by , where denotes the magnitude of signal reflection and represents the phase offset incurred by backscattering. Hence, the signal received at the receiver is given bywhere denotes the aggregate noise signal. Hence, the enhanced channel from the HAP to the receiver can be rewritten as follows:
(1) 
For single relay case, the channel model in (1) is degenerated to a simple representation as that in [27] and [25]. In the same way, we can rewrite the enhanced channel from the HAP to the active relay as follows:
(2) 
where denotes the complex channel from the passive relay to the active relay. Note that the enhanced channels and depend not only on the binary indicator , but also the complex reflection coefficient of each passive relay in the set . We observe that the phase is a critical design variable for channel enhancement while can be simply set to its maximum to increase the reflected signal power.
The channel models in (1) and (2) are simplified approximations as we omitted the interactions among different passive relays. In fact, besides reflecting the active radios’ signals, each passive relay can also reflect the backscattered signals from the other passive relays, thus creating a feedback loop, which makes it very complicated to characterize the the signal models exactly. However, this simplified model can be reasonable as the double reflections can be significantly weakened. Most importantly, it not only provides the basis for analytical study, but also enables us to derive a lower bound on the throughput performance of hybrid relay communications.
Iv Problem Formulation
In the sequel, we focus on this simplified channel model and formulate the throughput maximization problem for hybrid relay communications, which involves the joint optimization of the HAP’s transmit beamforming, the relays’ mode selection, and their operating parameters. Assuming a normalized noise power, the receiver’s SNR in the first hop is given by
(3) 
where is the Hermitian transpose of . By controlling the beamformer in the first hop, the HAP can adjust its information and power transfer to different relays. Each active relay can choose the different PS ratio to balance its power supply and demand, taking into account the power budget constraint as follows:
(4) 
where denotes the relay’s transmit power and is the energy harvesting efficiency. The relay’s overall power consumption can also take into account a constant circuit energy consumption, which however will not bring new challenge in the following algorithm design. The PS ratio indicates the portion of RF power that is converted by the energy harvester. The other part is then used for signal detection and thus the received signal of the relay is given by
where we define for notational convenience and
is the complex Gaussian noise with zero mean and normalized unit variance.
In the second hop, each active relay forwards the information to the receiver with the transmit power . All the relays’ signals will be combined coherently at the receiver. The power amplifying coefficient of the relay can be denoted as , which has to be optimized to maximize the overall throughput [6]. Meanwhile, the HAP can transmit the same information directly to the receiver with a new beamformer . This can enhance the reliability and data rate of information transmission from the HAP to the receiver. Hence, the received signals at the receiver is a mixture of the HAP’s direct beamforming and the relays’ forwarding, i.e.,
(5) 
The first two terms in (5) correspond to the amplified signals by the active relays. The third term represents the direct beamforming from the HAP. Due to the passive relays’ backscattering, the channel is also an enhanced version of from the relay to the receiver. Till this point, we can formulate the SNR in the second hop as follows:
(6) 
where , and denotes the diagonal matrix with the diagonal element given by .
Hence, the overall SNR at the receiver can be evaluated as . We aim to maximize SNR in two hops by optimizing the HAP’s beamforming strategies , as well as the relays’ radio mode selection and operating parameters, including the power splitting factor and the complex reflection coefficient :
(7a)  
(7b)  
(7c)  
(7d)  
(7e)  
(7f) 
The constraints in (7b) denote the HAP’s feasible beamforming vectors in two hops. We assume that the HAP’s transmit power is fixed at while the beamforming vectors can be adapted to maximize the throughput performance. Generaly is not necessarily the same as as the optimization of has to take into account the data transmissions to both the relays and the receiver. The binary variable determines the division of the relays in different modes. The constraint in (7d) ensures that the phase offset of each passive relay in set is fully controllable via load modulation [4]. The constraints in (7e) and (7f) determine the active relays’ transmit power in the second hop, which is upper bounded by the energy harvested from the HAP’s signal beamforming in the first hop. Note that the power budget constraint (7f) is based on a linear energy harvesting model with a constant power conversion efficiency , which may require an ideal circuit design for energy harvester. In practice, due to the nonlinearity of diode devices, the measurement results show a nonlinear energy harvesting model as proposed in [1]. In this case, we can consider a successive linear approximation of the nonlinear power budget constraint. Then our focus is the solution to the subproblem in each iteration, which has a similar form as the problem in (7).
Besides power budget constraints for the active relays, the passive relays in set are also wirelessly powered via energy harvesting. Let denote the constant power consumption of the passive relays. Then, we can formulate the passive relays’ power budget as follows:
where the scaling factor denotes the part of RF power harvested as energy. For simplicity, we assume that each passive relay sets the same reflection coefficient in both hops. The tuning of the reflection coefficient is subject to a limited range due to the antenna’s structural scattering effect [4]. The maximum magnitude is typically less than one. With a small , the passive relays’ power budget constraints easily hold, and hence we omit it in the problem formulation.
V Performance Maximization with Hybrid Relay Communications
It is clear that the optimization of the relays’ mode selection is combinatorial and difficult to solve optimally. Even with fixed relay mode, the throughput maximization is still challenged by the couplings of multiple relays in different modes, which have very different tranmission capabilities and energy demands. In particular, the passive relays’ operating parameters, e.g., the reflecting phase and magnitude, affect the active relays’ channel conditions. A joint optimization is required to characterize their couplings and optimize all users’ operating parameters simultaneously. However, the optimal solution is generally unavavilable. In the sequel, we propose to solve the throughput maximization problem in a decomposed manner. Firstly, assuming a fixed relay mode, we evaluate the enhanced channels and formulate the throughput maximization with only active relays, similar to that in [25]. Secondly, with the fixed beamforming strategy, we evaluate individual relays’ energy status or performance gain. This motivates our algorithm design to update the relays’ mode selection in an iterative manner.
Va Relay Performance with Fixed Mode Selection
Given a set of the passive relays and their reflection coefficients , the enhanced channels for active RF communications are given as in (1) and (2). Then, we can formulate the throughput maximization problem with the set of active relays alone, which becomes a conventional twohop relay optimization problem similar to that in [6] and [25]. Our target is to maximize by optimizing the HAP’s beamforming in two hops and the active relays’ PS ratios , subject to the relays’ power budget constraints:
(8a)  
(8b)  
(8c)  
(8d) 
Note that problem (8) will achieve different performance when the passive relays set different phase offsets . The phase optimization can follow the alternating optimization method. In particular, each passive relay initially sets a random phase offset , based on which we can optimize and by solving the problem (8). Given the solution to (8), we then turn to phase optimization sequentially for each passive relay, which will be detailed in Section VB.
VA1 Lower Bounds on Relay Performance
The throughput maximization (8) is still challenging due to the nonconvex coupling between different active relays in the objective (8a). The HAP’s beamforming strategy is also coupled with the relays’ PS ratio in a nonconvex form via the power budget constraint (8d). In the sequel, we provide a feasible lower bound on (8), which is achievable by designing the HAP’s beamforming and relaying strategies.
Proposition 1 ():
A feasible lower bound on (8) can be found by the convex reformulation as follows:
(9a)  
(9b)  
(9c)  
(9d)  
(9e) 
where is a constant. At optimum, the PS ratio of the relay is given by for .
The proof of Proposition 1 follows a similar approach as that in [14], and thus we omit it here for conciseness. With the fixed relay mode, the channel information and can be estimated by a training process. The objective function in (9a) then becomes linear and the constraints (9b)(9d) define a set of linear matrix inequalities^{1}^{1}1The constraint (9d) can be rewritten into two linear matrix inequalities.. Hence, the resulting problem can be efficiently solved by semidefinite programming (SDP) [23]. Once we find the optimal matrix solution , we can retrieve the HAP’s beamforming vector by eigendecomposition or Gaussian randomization method [15].
Though exact solution to (8) is not available, problem (9) provides a lower bound on the SNR performance, which can serve as the performance metric for the relay’s mode selection. It is clear that (9c) will hold with equality at optimum and thus we can verify the following property.
Proposition 2 ():
The proof of Proposition 2 can be referred to [24]. It implies that the lower bound on SNR performance in (8) can be evaluated directly from (10). However, by an inspection on (10), the SNR evaluation can be much less than the optimum of (8) if the direct link is practically weak due to physical obstructions. To this end, we also provide another lower bound on problem (8) by ignoring the direct link in the problem formulation. This corresponds to the case when the direct link is blocked or there is a long distance between the transceivers, e.g., [10]. In this case, the objective in (8a) is lower bounded by . Define and then the lower bound on (8) can be evaluated as follows:
(11a)  
(11b)  
(11c)  
(11d)  
(11e) 
where and we define for notational convenience. The upper bounds on and are defined as and , respectively, which depend on the relays’ PS ratio and the HAP’s beamforming strategy . Note that and are auxiliary decision variables in problem (11), coupled with the PS ratio and the beamformer via constraints in (11b)(11c).
VA2 Alternating Optimization Solution to (11)
To the best of our knowledge, there is no exact solution to the nonconvex problem (11). Practically, it can be solved by the alternating optimization method that improves the objective (11a) in an iterative manner with guaranteed convergence. In particular, with fixed and , the auxiliary variables and are subject to the fixed upper bounds and , respectively. As such, problem (11) can be viewed as the conventional network beamforming optimization problem with perfect channel information [6], which can be solved optimally in a closed form. However, given the feasible and , the optimization of and is still very difficult due to their bilinear coupling in the constraints. This implies that we require further approximation within each iteration of the alternating optimization method. To proceed, we first exploit the structural properties of problem (11) that shred some insight on the algorithm design.
Proposition 3 ():
The proof of property (ii) in Proposition 3 is straightforward by inspection. To verify property (i), we can simply rewrite the upper bound on as follows:
which is obviously increasing in and . Then, we focus on the property (iii). Note that only appears in the objective function. If (11c) holds with strict inequality at the optimum, we can simply improve the objective (11a) by increasing properly while keeping the other variables unchanged, which brings a contradition. Therefore, we can guarantee that at the optimum of (11).
As both and are increasing functions of the variable , we consider a replacement of by its smallest value . As such, we can further derive a lower bound on (11) by the following problem:
(12a)  
(12b)  
(12c)  
(12d) 
Here is given by , which can be easily reformulated into an SDP as follows:
(13a)  
(13b)  
(13c) 
By solving problem (13) with the interiorpoint algorithms [23], we can easily determine the HAP’s beamformer from the matrix solution via Gaussian randomization.
Till this point, we can employ the alternating optimization method to solve (12), which provides the lower bound on (11) at the convergence. With a fixed PS ratio , the feasible solution to (12) can be easily obtained by the network beamforming optimization in [6], and then we turn to update for each active relay based on the following property:
Proposition 4 ():
The proof is straightforward by checking the properties in Proposition 3, which reveal that is increasing in , while is decreasing in . At the optimum of (12), if for some , we can choose such that . Meanwhile, we have , which implies that the inequality constraint on can be relaxed. With the relaxed upper bounds on and , the network beamforming optimization in [6] will produce a higher objective value.
All the above derivations lead to the iterative procedure in Algorithm 1. The algorithm starts from a random initialization of . The HAP’s beamforming strategy can be obtained by solving (13) and fixed during the algorithm iteration. With the fixed , the upper bounds are also fixed, and thus we can apply the network beamforming optimization to solve in problem (12). According to Proposition 4, the update the PS ratio can be based on the feasibility check of the inequality constraints in (12b)(12c). In particular, for , we search for the relay with the largest constraint gap, defined as , and then we reduce properly by a small amount such that
where is a constant parameter. By solving above equation, we can determine easily as follows:
(14) 
Once is updated as shown in line 9 of Algorithm 1, we turn to optimize by network beamforming optimization [6].
Moreover, we can show that the lower bound derived by Algorithm 1 is also applicable to a singleantenna case. In parituclar, when the HAP has one single antenna, the enhanced channels from the HAP to different relays now become complex variables instead of vectors. The power budget constraints in (4) are degenerated to . Each active relay’s power amplifying coefficient can be similarly defined as where . As such, the performance maximization problem is given as follows:
(15a)  
(15b)  
(15c)  
(15d) 
By simply setting , problem (15) has the same form as that in (12) and hence it is solvable by Algorithm 1.
VB Iterative Algorithm for Relay Mode Selection
The previous analysis provides the SNR evaluation with fixed relay mode. This can serve as the performance metric to update each relay’s mode selection. The basic idea of an iterative algorithm for the relays’ mode selection is to start from the special case with all active relays and then update the relay mode one by one depending on the relay’s performance gain. For simplicity, we allow the mode switch of a single relay in each iteration. Hence, the number of iterations will be linearly proportional to the number of relays and the main computational complexity lies in the SNR evaluation given the division within each iteration. Such an iterative process continues until no further improvement can be achieved by changing the relays’ mode.
VB1 Evaluation of SNR Performance
Considering different channel conditions in the direct link, the SNR evaluation can be performed by solving either the SDP in (9) or the nonconvex problem in (11) by Algorithm 1. With the fixed , we can evaluate the SNR performance of each relay when it is in the passive mode. After iterating over all relays, we can switch the relay with the maximum SNR performance to the passive mode. It is obvious that the SNR evaluations in (9) and (11) both rely on the passive relay’s complex reflection coefficient , which is critical for the channel enhancement in (1) and (2) and thus effects the SNR evaluation. To maximize the SNR performance, the passive relays can simply set the magnitude of reflection to its maximum . However the complex phase is more difficult to optimize due to its couplings cross different relays. The dependence of different channels makes it difficult to enhance all active relays’ channels simultaneously.
The optimization of phase can follow the alternating optimization method. In particular, we can firstly optimize the beamforming strategy and the relays’ PS ratios to maximize the SNR performance. After that, we fix and and then turn to optimize the passive relay’s reflecting phase to further improve the SNR performance. Considering the lower bound in (10), the maximum SNR via phase optimization can be approximated as follows:
(16) 
where we simply replace in (10) by a known approximation . This approximation stems from the observation that different active relays are in general spatially distributed with indepedent channel conditions. As such, the optimization of complex reflection coefficient of one single passive relay has very limited capability to enhance all the active relays’ forwarding channels simultaneously. For example, the forwarding channel of an active relay can be enhanced by the passive relay with the reflection coefficient (i.e., ), while another channel may become weakened (i.e., ) due to the indepedence among different relays. Hence, we simply view the term as a constant approximated by . By this approximation, problem (16) can be further rewritten as follows:
(17) 
Let be a known matrix coefficient were
denotes the identity matrix. It is clear that problem (
17) aims to maximize the compound channel gain given the HAP’s beamforming strategy , which can be easily solved by onedimension search algorithm.Let SNR denote the objective in (17) when the relay is selected as the passive relay with the reflecting phase . Assuming a fixed reflection magnitude , the enhanced channel in (1) can be simplified as . The optimal phase can be simply obtained by a onedimension search method. In particular, we can quantize the continuous feasible region into a finite discrete set , where denotes the size of . Then, we can devise a onedimension search algorithm to optimize the phase parameter that maximizes SNR. The detailed solution procedure is presented in Algorithm 2, which alternates between the HAP’s beamforming optimization and the relays’ phase optimization by solving two subproblems in (9) and (17), respectively. The iteration terminates when SNR can not be improved anymore.
VB2 Simplified Schemes for Relay Mode Selection
The convergent value of SNR in Algorithm 2 can be used as a performance metric to evaluate the relay performance when the relay is selected as the passive relay. We denote this SNRbased performance metric as the MaxSNR scheme. However, we note that the SNR evaluation in (9) or (11) is quite complicated as it requires to solve optimization problems in an iterative procedure. In this part, besides the MaxSNR scheme, we also devise a set of simplified performance metrics based on the problem structure and the corresponding iterative algorithms to update the relays’ mode selections.

MaxDirectRate (MaxDR): With a strong direct link from the HAP to the receiver, the data rate contributed by the direct links becomes significant. The MaxDR scheme aims to maximize the data rate via the direct links, which is given by . The phase optimization follows a similar approach as that in (17) while the beamforming optimization can be significantly simplified comparing to that of the SDP (9).

MaxRelayRate (MaxRR): In contrast to the MaxDR scheme, the MaxRR scheme picks the passive relay which acheives the maximum relay performance without direct links by solving problem (12), which avoids the beamforming optimization in every iteration and thus has significantly reduced complexity than that of (9).

MaxDirectGain (MaxDG): The MaxDG scheme further simplifies the optimization in the MaxDR scheme, by focusing on the maximization of the channel gain instead of the data rate, which can be separated from the beamforming optimization. Hence, the MaxDG scheme avoids the iterative procedure, making it more efficient than the alternating optimization method for solving (12).

MinRFEnergy (MinRF): This scheme stems from the intuition that the passive radios are more energy efficient than the active radios. Hence, the relay with low energy supply may prefer to operate in the passive mode. Given the beamforming in (9), we can sort the active relays by the RF power at their antennas, i.e.,