Optimized Energy and Information Relaying in Self-Sustainable IRS-Empowered WPCN
Abstract
This paper proposes a hybrid-relaying scheme empowered by a self-sustainable intelligent reflecting surface (IRS) in a wireless powered communication network (WPCN), to simultaneously improve the performance of downlink energy transfer (ET) from a hybrid access point (HAP) to multiple users and uplink information transmission (IT) from users to the HAP. We propose time-switching (TS) and power-splitting (PS) schemes for the IRS, where the IRS can harvest energy from the HAP’s signals by switching between energy harvesting and signal reflection in the TS scheme or adjusting its reflection amplitude in the PS scheme. For both the TS and PS schemes, we formulate the sum-rate maximization problems by jointly optimizing the IRS’s phase shifts for both ET and IT and network resource allocation. To address each problem’s non-convexity, we propose a two-step algorithm to obtain the near-optimal solution with high accuracy. To show the structure of resource allocation, we also investigate the optimal solutions for the schemes with random phase shifts. Through numerical results, we show that our proposed schemes can achieve significant system sum-rate gain compared to the baseline scheme without IRS.
Index Terms:
Wireless powered communication network, intelligent reflecting surface, time scheduling, phase shift optimization.I Introduction
With nearly 50 billion Internet of Things (IoT) devices by 2020 and even 500 billion by 2030 [1], we have already stepped into the new era of IoT. Having the vision of being self-sustainable, IoT has observed the energy limitation as a major issue for its widespread development. Recent advances in energy harvesting (EH) technologies, especially radio frequency (RF) EH [2], opened a new approach for self-sustainable IoT devices to harvest energy from dedicated or ambient RF sources. This has led to the emergence of wireless powered communication networks (WPCNs), in which low-cost IoT devices can harvest energy from a dedicated hybrid access point (HAP) and then use the harvested energy to transmit data to the HAP [3]. The development of WPCNs has been a promising step toward the future self-sustainable IoT networks [4].
Although possessing significant benefits and attractive features for low-cost IoT networks, WPCNs are facing some challenges which need to be addressed before they can be widely deployed in practice. In particular, the uplink information transmission (IT) of IoT devices in WPCNs relies on their harvested energy from downlink energy transfer (ET) of the HAP. However, the IoT devices typically suffer from doubly attenuations of RF signal power over distance [3], which severely limits the network performance. Reducing the distance between the HAP and IoT devices is one solution to enhance EH efficiency and achieve greater transmission rates. However, this is not a viable option because IoT devices are randomly deployed in practice, and thus we may not be able to control all of them over their locations. Hence, more efficient and cost-effective solutions are required to enhance the downlink ET efficiency and improve the uplink transmission rate for WPCNs in order to guarantee that WPCNs can be seamlessly fitted into the IoT environment with satisfying performance.
Relay cooperation is an efficient way to enhance the performance of WPCNs, which can be classified into two categories of active relaying and passive relaying. Active relaying refers to scenarios in which the communication between a transmitter and its destined receiver is assisted by a relay which forwards the user’s information to the destination via active RF transmission [5]-[7]. However, active relaying schemes have several limitations. Particularly, EH relays need to harvest sufficient energy from the RF sources and use the harvested energy to actively forward information to the receiver. Due to the high power consumption of active relays, it may take a long time for the relays to harvest enough energy. This thus reduces the IT time of the network. Moreover, most active relays operate in the half-duplex mode, which further shortens the effective IT time, resulting in a network performance degradation. Full-duplex (FD) relays can relax this issue; however, complex self-interference (SI) cancellation techniques are needed at the FD relays to ensure that the SI is effectively mitigated [8]. In addition, the number of antennas at EH relays is usually limited due to hardware constraints, which also leads to a limited performance enhancement. Passive relaying exploits the idea of backscatter communication (BackCom) for assisting in the source-destination communication [9]-[11]. Specifically, BackCom relay nodes do not need any RF components as they passively backscatter the source’s signals to strengthen the received signals at the receiver. Accordingly, the power consumption of BackCom relay nodes is extremely low and no dedicated time is needed for the relays’ EH [12]. Nonetheless, as no active signal generation is involved and the passive relays simply reflect the received signal from the source, passive relaying schemes suffer from a poor performance.
Intelligent reflecting surface (IRS), consisting of a large number of low-cost reflecting elements, has recently emerged as a promising solution to improve the performance of wireless communication networks [13, 14]. This technology enables transmitting information without any need for generating new signals but recycling the existing ones [15]. In this way, IRS can adjust the communication environment and create favorable conditions for energy and information transmission without using energy-hungry RF chains. Having the capability of cooperating in downlink ET and uplink IT, IRS has several advantages over the conventional active and passive relaying techniques [13]. First of all, IRS is a cost-effective technology and it can be readily integrated into existing wireless communication networks without incurring high implementation costs. Furthermore, IRS is more energy- and spectrum-efficient as compared to conventional relaying methods because it consumes very low power and uses the limited spectrum resources more efficiently. IRS essentially works in the full-duplex mode without causing any interference and adding thermal noise, which further improves the spectral efficiency. Moreover, it is easy to increase the number of IRS elements to achieve higher performance gains.
I-A Background
IRS has recently been applied to various wireless communication networks and demonstrated promising results for improving the performance in terms of spectrum efficiency (SE) and energy efficiency (EE). References [16] and [17] consider the EE maximization problem in an IRS-assisted multi-user downlink communication network. The authors jointly optimize the power allocation at the AP and the phase shifts at the IRS and show that the proposed IRS-assisted communication remarkably outperforms the conventional relay-assisted communication in terms of EE. The authors in [18] propose an architecture where a few IRS elements are assumed to be active. Based on the proposed architecture, the authors develop two solutions using compressive sensing and deep learning for designing IRS’s reflection matrices. The authors in [19] exploit deep reinforcement learning based algorithms to jointly design the transmit beamforming at the base station and phase shifts at the IRS to maximize the sum-rate of downlink multiple-input single-output (MISO) systems. A low-complexity channel estimation protocol is proposed in [20], which does not require any prior knowledge on channel state information (CSI) or any active participation from IRS. The authors design the near-optimal active beamforming at the power beacon and passive beamforming at the IRS in order to maximize the received power at an EH user. In [21], the authors propose two efficient algorithms for finding the optimal beamformer at the AP and phase shifts at the IRS in an IRS-assisted MISO communication system. The authors in [22] study the problem of transmit power minimization in a multi-user downlink communication network by jointly optimizing the active transmit beamforming at the AP and passive reflect beamforming at the IRS subject to the users’ individual signal-to-interference-plus-noise ratio (SINR) constraints. Compared to the conventional massive MIMO system, the proposed IRS-enhanced model in [22] can considerably reduce the required transmit power. [23] and [24] study the integration of IRS with simultaneous wireless information and power transfer (SWIPT) technology, where the transmit precoders at the AP and the passive beamforming at the IRS are jointly optimized for maximizing the weighted sum-power at EH users [23] and weighted sum-rate at information receiving users [24]. Physical layer security and outage probability analysis in IRS-assisted MISO networks are investigated in [25] and [26], respectively. Signal-to-noise-ratio (SNR) maximization problem in a self-sustainable single-user IRS-assisted MISO communication system is studied in [27], where IRS elements use part of the downlink information signal for harvesting their required energy.
A survey on recent research efforts in the area of IRS can be found in [28].
I-B Motivations
Although IRS has lately received significant interests from the research community, it is still at the very early stage of development and more investigations are needed to fully capture the potentials of IRS and make it applicable to practical scenarios. Specifically, the integration of IRS technology with WPCN is a great step toward the realization of efficient and self-sustainable IoT networks, which has not been well investigated in the literature. Recently, a few research works have investigated the application of IRS for improving the performance of WPCNs [29, 30]. In [29], the authors study the application of IRS for WPCN performance enhancement, where IRS elements assist in downlink ET from the HAP to the users and uplink IT from users to the HAP. The authors in [30] propose a similar idea to use the IRS as a hybrid energy and information relay, where the user cooperation is also investigated for a two-user WPCN scenario. These preliminary works on the integration of IRS with WPCN provide some insights on the performance enhancements offered by using IRS in WPCNs. However, this integration needs to be studied more deeply with practical considerations for the network setup and network elements.
One of the most important points that is often overlooked in the studies on IRS is the IRS’s power consumption. Although IRS elements passively reflect the incident signals, the power consumption of the IRS cannot be neglected [16, 17, 27]. However, the majority of the works in this area (e.g., [18]-[26]) assume that the IRS’s power consumption is negligible because it does not perform complex signal processing tasks. In practice, the power consumption of IRS depends on the type and characteristics of its reflecting elements [16, 17]. For example, the values of each reflecting element’s circuit power consumption are 1.5 and 6 mW for 3- and 5-bit resolution phase shifting, respectively [17]. As the number of IRS elements is typically large, the circuit power consumption of the IRS can be even comparable to its power supply and cannot be neglected.
In the self-sustainable IoT networks, devices are expected to operate in an uninterrupted manner and have theoretically perpetual lifespans. Considering the non-negligible power consumption of IRS elements, it is important to propose efficient strategies which can keep the IRS operational for very long periods. Although embedded batteries can power the IRS temporarily, they cannot be relied on for the long-term functionality and uninterrupted operation of the IRS. Wired charging may also be unavailable if the IRS is deployed in inaccessible places. Thus, equipping IRS elements with EH modules can resolve these issues and make the IRS energy-neutral [15, 31]. This is our main motivation for studying a self-sustainable IRS-empowered WPCN, where the EH-enabled IRS, powered by energy transmission of the HAP, can act as a hybrid energy and information relay assisting in both downlink ET and uplink IT.
I-C Contributions
We study a self-sustainable IRS-empowered multi-user WPCN, where the IRS is equipped with an EH circuit to harvest RF energy from the HAP to power its operations. Inspired by the conventional wireless-powered active relays [32], time-switching (TS) and power-splitting (PS) schemes are proposed to enable the IRS to harvest energy from the RF signals transmitted by the HAP. In the TS scheme, the ET phase is split into two sub-slots, where the IRS harvests energy in the first sub-slot and assists in the downlink ET to the users in the second sub-slot. Compared to the conventional TS scheme [32], the proposed TS scheme can efficiently improve the amount of harvested energy at the users. In the PS scheme, the IRS harvests energy from the HAP’s signal and assists in the downlink ET to the users by adjusting its amplitude reflection coefficients in the ET phase. Compared to the conventional PS scheme [32], the proposed PS scheme can enhance both ET and IT efficiency and is more spectrum-efficient. To make our study applicable to practical systems, we consider a piece-wise linear EH model for the IRS and the users to account for the saturation behavior of practical EH systems [33]-[36]. We investigate the problem of sum-rate maximization for both TS and PS schemes and optimize the IRS phase shift design and network resource allocation jointly with EH time and amplitude reflection coefficients of the IRS.
The main contributions of this paper are summarized as follows:
-
•
We propose a self-sustainable IRS-empowered WPCN, where a wireless-powered IRS acts as a hybrid relay to improve the performance of WPCN in both downlink ET from the HAP to the users and uplink IT from users to the HAP.
-
•
To enable energy collection and hybrid relaying functionalities at the IRS, we propose more efficient TS and PS schemes, which can enhance the ET efficiency from the HAP to the users and assist in the uplink information transmission. We consider a piece-wise linear EH model for the IRS and the users, which is mathematically tractable and is able to capture the saturation effect of practical energy harvesters.
-
•
We study the system sum-rate maximization problem for the TS scheme by jointly optimizing the IRS’s phase shift designs in both ET and IT phases, time allocation for the IRS and users’ EH, time allocation for each user’s IT, and the users’ power allocation. To deal with the non-convexity of the formulated problem, we propose a two-step algorithm to achieve the near-optimal solution: in the first step the phase shifts for the IT are obtained in closed-form, while an efficient method by using one-dimensional search, semidefinite relaxation (SDR) and Gaussian randomization is designed for optimizing the IRS phase shifts in the ET phase, time allocation and power allocation in the second step. In particular, we obtain a closed-form solution for the optimal IRS’s EH time and discuss its implications.
-
•
We then investigate the sum-rate maximization problem for the PS scheme and jointly optimize the IRS’s phase shift design in both ET and IT phases, time allocation for the EH and IT phases, power allocation at the users, and the amplitude reflection coefficient in the EH phase, using a similar two-step algorithm as for the TS scheme. In particular, we analyze the condition for activating the IRS in the PS scheme and obtain the optimal amplitude reflection coefficient as a function of the EH time, from which some interesting observations are revealed.
-
•
Finally, we evaluate the performance of our proposed schemes via numerical simulations which show that our proposed schemes can achieve significant system sum-rate gain compared to the baseline WPCN protocol.
I-D Organization
This paper is organized as follows. Section II describes the system model of the proposed IRS-empowered WPCN for both TS and PS schemes. Sections III and IV investigate the sum-rate maximization problems for TS and PS schemes, respectively. Section V evaluates the performance of the presented algorithms by conducting numerical simulations and Section VI concludes the paper.
II System Model
As illustrated in Fig. 1, we consider an IRS-assisted WPCN, consisting of an HAP with stable power supply, energy-constrained users (denoted by ), and an energy-constrained IRS. The IRS and users are each equipped with an EH circuit (rectenna) to harvest energy and an energy storage to store the harvested energy. The HAP serves as a central control point for the network, which coordinates the transmissions among all devices and also has the capability and constant energy supply for performing computational tasks. The HAP and users have single antenna each.111The model can be straightforwardly extended to the scenario that the HAP is with multiple antennas, which will be briefly discussed in Remark 3. The IRS is composed of passive reflecting elements, which can be configured to direct the incident signals to desired directions. The IRS assists in both downlink ET from the HAP to the users and uplink IT from the users to the HAP. The EH and energy/information relaying at the IRS are controlled by an attached micro-controller.
The downlink channels from the HAP to , from the HAP to the IRS, and from the IRS to are denoted by , , and , respectively. The counterpart uplink channels are denoted by , , and , respectively. All channels are assumed to be quasi-static flat fading, which remain constant during one block but may change from one block to another [23]. We assume that the channel state information (CSI) of all links is perfectly known.222The CSI of all links can be precisely obtained by existing channel estimation techniques [20, 37]. In the future work, the effect of channel estimation errors on system performance will be investigated.



The transmission block with a duration of seconds, is divided into two phases, i.e., ET phase and IT phase. In the ET phase, the HAP transfers energy to the users and IRS in the downlink. The IRS uses the HAP’s signals for its own EH and energy relaying to the users. In the IT phase, the users use the harvested energy to transmit data to the HAP with the assistance of the IRS. Without loss of generality, we consider a normalized unit transmission block time in the sequel, i.e., second. The details of the ET and IT phases are shown in Fig. 2 and elaborated in the following subsections.
II-A Energy Transfer Phase
As mentioned earlier, the IRS is assumed to be energy-constrained, which needs to harvest energy from the HAP for powering its relaying operations. In this regard, we design efficient TS and PS schemes for the IRS.
II-A1 Time-switching scheme
For the TS scheme, the ET phase with the duration of 333The unit of all time coefficients is seconds. is divided into two sub-slots, having the duration of and , respectively, which satisfy . The users can harvest energy over the entire ET phase. For the IRS, it will spend the first sub-slot in the ET phase for its own EH and the second sub-slot for improving the EH efficiency at the users. In particular, in the first sub-slot, all incident signals at the IRS from the HAP are transferred to the EH harvester by setting the amplitude reflection coefficients to be zero, and thus no incident signals will be reflected by the IRS. While in the second sub-slot, the IRS cooperates with the HAP by adjusting its elements’ phase shifts to enhance the total received signal power at the users. The transmission block structure for the TS scheme is illustrated in Fig. 2 (a). Denote the transmit signal in the ET phase as , where is the transmit power and is the energy-carrying signal with .
The received signals at the IRS and in the first sub-slot are expressed as
(1) | |||
(2) |
where and denote the additive white Gaussian noises (AWGNs) at the IRS and , respectively. Note that the noise power is usually very small and ineffective for EH and can be thus neglected. Hence, the received power at the IRS, denoted by , is expressed as . Similarly, the received power at during is given by .
In the second sub-slot, the IRS assists in the downlink ET. The phase shift matrix of the IRS during is denoted by , where denotes the reflection efficiency and is typically set as a constant [38], and are the amplitude reflection coefficient and the phase shift of the -th element, respectively. Let , where . For the TS scheme, since the IRS only harvests energy during , all incident signals at the IRS during can be reflected to enhance the EH efficiency, i.e., [13]. Let . During , the received signal at for the TS scheme is given by
(3) |
The received power of during is then given by .
In practice, the EH circuits usually lead to a non-linear rectification efficiency, i.e., the RF power-to-direct current power conversion is a non-linear function with respect to the received RF power [39, 40]. In particular, the harvested power first improves with the increase of received power but finally becomes saturated when the received power is high [40]. To approximate the non-linear EH characteristics and account for the saturation region of practical energy harvesters, we employ a two-piece linear EH model,444There also exist other EH models, e.g., the logistic function based non-linear EH model [39] and the multi-piece linear EH model [41]. However, it is noted that the two-piece linear EH model is sufficiently accurate for modeling the behavior of practical EH circuits. Compared to the logistic function based non-linear EH model, the piece-wise linear EH model is mathematically appealing and easily tractable. In addition, the results obtained from the two-piece linear EH model can be straightforwardly extended to the multi-piece linear EH model. which is also widely used in the literature, e.g., [34]-[36]. According to this model, the harvested power is calculated as
(4) |
where is the EH efficiency in the linear regime,555In practice, the EH efficiency in this regime is not strictly linear. However, as mentioned in Footnote 4, assuming a constant is still sufficiently accurate for modeling the practical EH circuits. is the received power, and denotes the saturation power, beyond which there will be no increase in the amount of the harvested power. Therefore, the harvested energy at the IRS and can be obtained as
(5) | ||||
(6) |
where and represent the saturation power of the IRS and , respectively.
II-A2 Power-splitting scheme
Different from the TS scheme, the dedicated EH time is not required in the PS scheme and the IRS harvests energy from the HAP by adjusting the amplitude reflection coefficients ()666Adjusting the reflection coefficient can be achieved by using electronic devices such as positive-intrinsic-negative (PIN) diodes, field-effect transistors (FET), micro-electromechanical system (MEMS) switches, and variable resistor loads [13, 42]., as illustrated in Fig. 2 (b). To be specific, only a part of the HAP’s energy signals is fed into the IRS’s EH unit for harvesting and the remaining part is reflected by the IRS to enhance the amount of harvested energy at the users.
It is assumed that all the amplitude reflection coefficients of the IRS elements have the same value, i.e. .777In practice, the elements can have different amplitude reflection coefficients. However, the setting will greatly complicate the circuit design of the IRS as different circuits should be integrated to control the amplitude reflection coefficient and phase shift independently at each element [13, 42]. To guarantee the operations of the self-sustainable IRS, we should simplify its circuit design to reduce its circuit power consumption, which can be achieved by setting all amplitude reflection coefficients to be the same. The received signal at in the ET phase for the PS scheme is thus given by
(7) |
The harvested energy of the IRS and for the PS scheme is then given by
(8) | |||
(9) |
II-B Information Transmission Phase
In the IT phase, the users transmit information to the HAP via time division multiple access, using the harvested energy in the ET phase. Denote the duration of IT for as . Let be the information-carrying signal of with unit power. The transmit signal of during is then expressed as , where is ’s transmit power and satisfies
(10) |
with being the circuit power consumption of . As the amplitude reflection coefficients are set to be the same, the IRS’s circuit power consumption is mainly caused by performing each element’s phase shifting [16, 17]. The other power consumptions, such as powering the EH circuit and signaling overhead, can be considered to be negligible [27, 32, 43]. By denoting the power consumption of each element as , the circuit power consumption of the IRS is thus expressed as . To power its operations, IRS needs to harvest sufficient energy in the ET phase. We assume that all the harvested energy stored in the energy storage can be used to power the IRS’ circuits, the following constraints are thus held:
(11) | |||
(12) |
for TS and PS schemes, respectively. Note that the power consumption of the IRS in the first sub-slot of the TS scheme is neglected because the IRS’s power consumption is mainly determined by the reflection operation [16, 17], which do not take place during .
Denote the phase shift of the -th element for ’s IT as . Then, the phase shift matrix during is denoted by , where , , and . Note that we have set the amplitude reflection coefficients to be 1 to maximize the signal reflection in the IT phase [13]. The received signal at the HAP from , denoted by , is thus given by
(13) |
where is the AWGN at the HAP. The SNR at the HAP during , denoted by , is expressed as . The achievable rate from to the HAP in bits/second/Hz is then formulated as
(14) |
III Sum-rate maximization for the TS scheme
In this section, we aim to maximize the system sum-rate by jointly optimizing the phase shift design at the IRS in both ET and IT phases, time scheduling of the network, and power allocation at the users. The constraints for the TS scheme are given as follows: , , , , , , , , . The optimization problem is formulated as
(P1) | ||||
s.t. |
where , , and .
III-A Near-optimal solution to P1
It is obvious that P1 is a non-convex optimization problem due to the coupling of variables in the objective function and the constraints, and convex optimization techniques cannot be used to solve it directly. In the following, we propose a two-step algorithm to solve the sum-rate maximization problem in P1. Specifically, we first obtain the optimal phase shifts for the IT in closed-form and then propose an efficient algorithm to solve the simplified problem.
III-A1 Optimal phase shift design for IT
We first present a proposition for the optimal design of phase shifts of the IRS for the IT.
Proposition 1.
The optimal IRS phase shifts for the IT during () are given by
(15) |
where is the -th element of , is the -th element of , and represents the phase of .
Proof.
Refer to Appendix A. ∎
III-A2 Optimizing phase shift design for ET, time scheduling, and power allocation
According to Proposition 1, P1 can be simplified as
(P2) | ||||
s.t. |
where . Note that solving P2 is equivalent to solving P1. P2 is still non-convex because the variables are coupled in the objective function and the constraints. To make P2 tractable, we introduce , where and set . Let , and , where and . Based on these new variables, the constraint C2 is recast as follows:
(16) |
where
Then, P2 can be equivalently rewritten as
(P2.1) | ||||
s.t. | ||||
Due to the rank-one constraint in C14 and coupling of and in C10, P2.1 is still non-convex and difficult to be solved directly. However, it is straightforward to obtain the optimal duration of the first sub-slot in the ET phase, i.e., , as stated in the following proposition.
Proposition 2.
The optimal duration of the first sub-slot in the ET phase can be obtained as
(17) |
Proof.
Refer to Appendix B. ∎
Remark 2.
From Proposition 2, we can observe that the duration of the first sub-slot in the ET phase is mainly determined by the IRS’s setting, e.g., the number of reflecting elements, each element’s circuit power consumption, and the saturation power for the EH. With a higher circuit power consumption for each element, the IRS needs more time to harvest sufficient energy, which causes a shorter time for other network operations, e.g., users’ EH with the assistance of IRS and users’ IT. Furthermore, if the saturation power of IRS () is small such that , adding more reflecting elements will increase the IRS’s circuit power consumption, which subsequently increases the EH time. Otherwise if , increasing the number of elements provides additional transmission links between the HAP and IRS, and thus more energy from the HAP can be transferred to the IRS. Therefore, if the increase of the IRS’s circuit power consumption is smaller than that of its harvested power, the EH time of IRS can even be reduced by increasing the number of elements.
We now proceed to solve P2.1 with obtained in Proposition 2. For solving P2.1, we first fix and optimize time and energy allocation in the IT phase as well as the IRS phase shift design for the ET phase. We can then find the optimal by a one-dimensional search over . Denote . With fixed , P2.1 is reformulated as
(P2.2) | ||||
s.t. | ||||
P2.2 is still non-convex due to the rank-one constraint in C14, and its globally optimal solution is thus difficult to obtain. However, using the semidefinite relaxation (SDR) technique [44], we can relax the rank-one constraint to obtain a convex semidefinite programming (SDP) problem [45], which can be optimally solved using convex optimization tools, e.g., CVX [46]. However, the solution obtained for the relaxed version of P2.2 by CVX may not satisfy the rank-one constraint. The Gaussian randomization method is then employed to construct a rank-one solution to P2.2 from the solution obtained by CVX. Note that the constructed rank-one solution can be a near-optimal solution to P2.2 as it is constructed and searched with quite large times of randomization [47].
Denote the optimal solution to the relaxed problem as . The singular value decomposition (SVD) of is expressed as , where and are the unitary matrix and diagonal matrix, respectively. Then, the approximate solution to P2.2, denoted by , can be constructed as follows
(18) |
where is a random vector with . Note that as the objective function is an increasing function of , C10 must be an equality at the optimal solution. Therefore, based on the generated random vectors, the energy allocation of the users is computed as
(19) |
where means . We generate times of random vectors and compute the corresponding objective function values for P2.2. The near-optimal solution to P2.2, denoted by , is the one achieving the maximum objective function value. The near-optimal , is finally recovered by
(20) |
where represents that the first elements of are taken, denotes the -th element of . It has been numerically and mathematically proved in the literature that the SDR technique followed by Gaussian randomization can provide a good approximation of the optimal solution (see [44] and the references therein).
The procedure for solving the sum-rate maximization problem for the TS scheme is summarized in Algorithm 1, in which the two steps are implemented sequentially. According to [44], the worst-case computational complexity of Algorithm 1 is , where is the computational accuracy for using the interior-point method in CVX, and is the number of iterations for updating . By running Algorithm 1 at the HAP, we can obtain the globally-optimal solution for P1 approximately, the analysis of which is given as follows. First, we can obtain the optimal phase shifts for the IT and the optimal duration of the first sub-slot in the ET phase, which are the globally optimal solutions to P1. Second, the globally optimal duration of the second sub-slot in the ET phase can also be found by setting an appropriate step size. Third, the SDR technique followed by quite large times of randomization based on the Gaussian randomization scheme can guarantee at least approximation of the maximum objective function value achieved by solving P2.2 [47].
Remark 3.
For the multi-antenna HAP scenario, the sum-rate maximization can be achieved by jointly optimizing the phase shift designs, time scheduling, power allocation at the users, and transmit beamforming at the HAP, denoted by . The formulated problem can be solved by a similar two-step algorithm, where the problem is decoupled into two sub-problems. In particular, the first sub-problem optimizes the phase shifts for the IT, which can be solved by performing the SDR and Gaussian randomization. The second sub-problem can be solved by using the block coordinate descent (BCD) techniques. Specifically, we can divide the variables into two blocks, i.e., and , and iteratively optimize one of them with the other one fixed in an alternating manner. A sub-optimal solution can be finally obtained by sequentially solving the two sub-problems. However, the algorithm’s complexity for this scenario is much higher, i.e., , where is the number of antennas at the HAP, and is the designed accuracy for the BCD method. Hence, this scenario causes a much higher implementation cost and is not appropriate for the self-sustainable IoT network.
III-B Random phase shifts with optimized resource allocation for the TS scheme
To reduce the computational complexity and show more insights about resource allocation, we consider a special case with random design of phase shifts and focus on the time and power allocation optimization in the IRS-assisted WPCN. As will be shown in Section V, using IRS is beneficial for improving the performance of WPCN even with randomly designed phase shifts [22, 48]. Letting , we have
(21) |
and the sum-rate maximization problem with random phase shifts is formulated as
(P3) | ||||
s.t. |
where . The constraint C15 is an equality at the optimal solution as we discussed in Section III-A2. Hence, we have
(22) |
Substituting (III-B) into , we have
(23) |
where , and . Proposition 2 holds here as well. Hence, P3 is modified as
(P3.1) | ||||
s.t. | ||||
It can be verified that P3.1 is a convex optimization problem [45], which can be solved by standard convex optimization techniques, e.g., Lagrange duality method. The Lagrangian of P3.1 is given by
(24) |
where is the Lagrange multiplier.
Proposition 3.
With random design of phase shifts, the optimal time scheduling for the TS scheme is given by
(25) | ||||
(26) |
where is the unique solution of , and is the optimal dual variable.
Proof.
Refer to Appendix C. ∎
IV Sum-rate maximization for the PS scheme
In this section, we investigate the optimal solution to the sum-rate maximization problem for the PS scheme. The problem is formulated as
(P4) | ||||
s.t. | ||||
IV-A Near-optimal solution to P4
Similar to P1, P4 is a non-convex optimization problem due to the coupled variables in the objective function and the constraints. It is straightforward to observe that Proposition 1 also holds for P4. Accordingly, P4 can be equivalently reformulated as
(P4.1) | ||||
s.t. |
Lemma 1.
If the PS scheme is employed at the IRS, the following condition must be met in order to guarantee that IRS can be used for assisting in downlink ET and uplink IT:
(27) |
Proof.
Refer to Appendix D. ∎
Remark 4.
From Lemma 1, we can observe that the PS scheme cannot always be used, i.e., if (27) is not satisfied. That is to say, the applications of the PS scheme are restricted by the IRS’s setting (i.e., the number of IRS elements, the circuit power consumption, and the saturation power) and the network environment (i.e., the transmit power at the HAP and the channel power gain between the HAP and IRS). If , we can increase the transmit power at the HAP and/or reduce the distance between the HAP and IRS to enable the PS scheme. However, if , the maximum number of IRS elements for enabling the PS scheme is . Compared to the PS scheme, the TS scheme is free from the limitation and can be applied more widely.
In the following, we investigate P4.1 under the condition that (27) is satisfied, because otherwise the IRS is not able to improve the performance of WPCN. Following the same steps as in Section III-A2, the sum-rate maximization problem is formulated as
s.t. | |||
where
Similarly, solving P4.2 is equivalent to solving P4. From P4.2, we first obtain the following proposition about the optimal amplitude reflection coefficient.
Proposition 4.
The optimal value of the amplitude reflection coefficient is obtained as
(28) |
where .
Proof.
Refer to Appendix D. ∎
For solving P4.2, we first fix and optimize other variables. The optimal value of can then be obtained by a one-dimensional search over . Given , the optimal value of can be found from Proposition 4 and we will have the following optimization problem:
(P4.3) | ||||
s.t. |
After the relaxation of the rank-one constraint in C14, P4.3 is similar to P2.2 in Section III-A2 and can be solved following the same procedure. For brevity and to avoid repetition, we do not explain the details of solving P4.3 here.
Algorithm 2 describes the process of solving the sum-rate maximization problem for the PS scheme, which is implemented at the HAP. Similar to what has been mentioned for Algorithm 1, the computational complexity of Algorithm 2 is , where is the number of iterations for updating . Again, by setting the appropriate step size for updating and relatively large number of randomizations for the Gaussian randomization method, we can obtain the near-optimal solution to P4.
IV-B Random phase shifts with optimized resource allocation for the PS scheme
Similar to Section III-B, we consider the random design of phase shifts for the PS scheme and optimize the resource allocation in the network. With randomly generated phase shifts and after setting , we have the following resource allocation problem:
(P5) | ||||
s.t. | ||||
where . It can be observed that Proposition 4 also holds for P5. Due to the non-convexity of C20, it is still challenging to solve P5. Hence, we first fix and optimize the time and energy allocation in the IT phase. We then find the optimal value of by searching over .
We know from previous discussions that C20 must be met with equality at the optimal solution, i.e.,
(29) |
Consequently, given and , P6 is rewritten as
(P6.1) | ||||
s.t. |
where . The Lagrangian of the above convex problem is given by , where is the Lagrange multiplier.
Proposition 5.
With fixed and , the optimal time allocation in the IT phase for the PS scheme is given by
(30) |
where is the unique solution of , and is the optimal dual variable.
The proof of Proposition 5 is similar to that of Proposition 3 and is thus omitted for brevity. Updating by the one-dimensional search method, we can obtain its optimal solution. After that, the optimal energy allocation can be easily found via (29) and Proposition 5.

V Performance Evaluation
In this section, we present numerical results to evaluate the performance of the proposed solutions for the IRS-assisted WPCN. The simulated network topology is a 2-D coordinate system as shown in Fig. 3, where the coordinates of the HAP and the IRS are given as (0,0) and , the users are randomly deployed within a circular area centered at with radius 1 m. We follow the channel model considered in [23]. In particular, the large-scale path-loss is modeled as , where is the path-loss at the reference distance m and set at dB [25], denotes the distance between two nodes, and is the path-loss exponent. The path-loss exponents of the links between the HAP and users are assumed to be since the users are randomly deployed, while the path-loss exponents of the links between the HAP and IRS and between the IRS and users are set at because the IRS can be carefully deployed to avoid the severe signal blockage. The small-scale fading coefficients are modeled to be Rician fading. In particular, the small-scale channel from the HAP to the IRS can be expressed as , where is the Rician factor for the HAP-IRS link, is the deterministic line of sight (LoS) component, and is the Rayleigh fading component with circularly symmetric complex Gaussian random variables with zero mean and unit variance. The small-scale channels for the other links are similarly defined. The Rician factors for the HAP-IRS link, the HAP- link, and the IRS- link are set at , , and , respectively. Unless otherwise stated, other parameters are given as follows: , , dBm, mW, W, mW [17], mW, , , dBm, m, m, and m. The scheme with random design of phase shifts, the scheme with random EH time, and the scheme without IRS are used as benchmarks for performance comparisons.
Fig. 4 shows the influence of the HAP’s transmit power on the average system sum-rate. As expected, the average sum-rate is improved with the increase of the HAP’s transmit power because the users can harvest more energy when the HAP’s transmit power is higher. Further, according to Proposition 2, the time needed for the IRS’s EH in the TS scheme can be reduced when the transmit power of the HAP is increased. This gives more time for the IRS to assist in downlink ET from the HAP to the users, which boosts the harvested energy at the user and consequently improves the sum-rate. As for the PS scheme, increasing the HAP’s transmit power results in higher amplitude reflection coefficient according to Proposition 4, which enhances the users’ harvested energy. It can be seen that our proposed schemes with optimized phase shift design outperform the benchmark ones for both the TS and PS schemes. The figure also shows that when dBm, there is no gain in using the IRS for improving the performance of WPCN for the PS scheme, which is consistent with what has been noted in Lemma 1. When dBm, the average sum-rate achieved by the PS scheme becomes stable. It is because the maximum harvested power by the IRS and users is limited by the values of their saturation power.


In Fig. 5, we study the impact of the number of IRS reflecting elements on the average sum-rate. It can be clearly observed that our proposed schemes can achieve a significant gain in terms of the average sum-rate compared with other schemes. As the number of IRS elements increases, the sum-rate achieved by our proposed schemes first increases and then reduces. It is because increasing the number of elements can provide additional transmission links for the ET and IT but also increases the circuit power consumption of the IRS, which thus reduces the users’ IT time . If the improved channel power gains can compensate for the reduction of IT time, the sum-rate can be improved; otherwise, the sum-rate reduces. This observation indicates that setting an appropriate number of IRS elements is important for performance enhancement. For the scheme without IRS, the sum-rate is smallest. It is because the received power at each user from the HAP through the direct link only is limited, thus more time is required to harvest energy to power its circuit, and the remaining energy and time for the IT is relatively small. It is also worth mentioning that even the schemes with random phase shifts and the schemes with random EH time can bring performance gains to the WPCN. That is because the RF energy can still be transferred from the HAP to the users through the reflecting links [48]. It endorses the effectiveness of using the IRS for performance enhancement.


In Fig. 6, we study the effect of the number of network users on the average sum-rate. Again, the proposed IRS-assisted WPCN with optimal phase shift design notably outperforms the other schemes. It can be observed that the average sum-rate is increasing with the number of users because more energy can be harvested with the increase of the number of users. Nevertheless, the average sum-rate does not increase when the number of users reaches a high number, e.g., over 10 users. The reason for this observation is that adding new users implies that more time is needed for the IT phase, which in consequence decreases the ET phase duration. Shorter ET duration in the TS scheme means that less time will be left for the IRS to assist in the downlink ET. In the PS scheme, the IRS needs to decrease its amplitude reflection coefficient to compensate for the loss of energy incurred by shortening the ET duration. Therefore, the gain brought by incrementing the number of users is neutralized by shortened ET time and the average sum-rate converges to an upper bound.

Next, we investigate the effect of users’ locations on the sum-rate performance. As shown in Fig. 7, increasing results in sum-rate reduction because as increases, the users move further from both the HAP and IRS. Therefore, the signals received by the users in the ET phase from both the HAP and the IRS become weaker. Similarly, the signals received by the HAP in the uplink IT also get weaker. Once again, the proposed schemes can significantly outperform the other benchmark schemes.
Finally, we investigate the impact of the IRS’s location on the sum-rate performance in Fig. 8. It can be observed that increasing the distance between the HAP and IRS reduces the sum-rate for the proposed schemes. It is because as the distance increases, the IRS has to spend more time to harvest sufficient energy to power its operations, which thus reduces the IT time for users. Compared to the TS scheme, the PS scheme is more susceptible to the IRS’s location. The reason is that for the PS scheme, as the distance increases, not only the users’ IT time but also the amplitude reflection coefficient at the IRS will be reduced.
VI Conclusions
This paper has proposed the hybrid-relaying scheme empowered by a self-sustainable IRS to enhance the performance of WPCN, where the IRS is deployed to improve the efficiency of downlink ET from the HAP to a number of users and uplink IT from the users to the HAP. In addition, we have proposed the TS and PS schemes for the IRS to harvest sufficient energy from the HAP to power its operations and investigated system sum-rate maximization problems for both schemes. To address the non-convexity of each formulated problem, we have developed the two-step algorithms to efficiently obtain the near-optimal solution with satisfying accuracy. The special problems with random phase shifts have also been investigated to revel the structure of time and energy allocation. Then, we have performed simulations to evaluate the superiority of our proposed schemes, which have shown that our proposed schemes can achieve remarkable sum-rate gain compared to the baseline WPCN without IRS. From simulation results, we have also observed that the PS scheme can achieve a better performance than the TS scheme if the transmit power at the HAP is large enough or the channel between the HAP and IRS is strong. However, compared to the PS scheme, the TS scheme can be more widely applied because it is free from the constraint defined in Lemma 1 for the PS scheme.
Appendix A Proof of Proposition 1
It is straightforward that is an increasing function with respect to for . Therefore, the optimal solution of P1 is found when is maximized. In addition, only depends on . As a result, for any given and feasible and , maximizing the objective function of P1 with respect is equivalent to solving the following problem for
(P-A) | ||||
s.t. |
The objective function can be rewritten as , where . It is obvious that the maximum of is achieved if , i.e., = . Denoting and , we have . Then, P-A can be rewritten as
(P-B) | ||||
s.t. | ||||
Appendix B Proof of Proposition 2
It can be verified that the objective function of P2.1 is an increasing function with respect to and for . Therefore, at the optimal solution, C10 must be met with equality. The constraint C1 must also be satisfied with equality, because otherwise we can decrease and increase , which results in more harvested energy at the users and larger transmit energy . We can also observe that the right hand side of C1 is increasing with respect to . Thus, the constraint C3 must be met with equality at the optimal solution because otherwise we can always increase as a result of which and users’ harvested energy can also be increased. Similarly, the constraint C4 must also be an equality at the optimal solution as otherwise we can increase , leading to the increase of and . Based on the three equalities from the constraints C1, C3 and C4, we can straightforwardly obtain the optimal value of as given by (17).
Appendix C Proof of Proposition 3
The dual function of P3.1 is given by . Karush-Kuhn-Tucker (KKT) conditions are both necessary and sufficient for the optimality of P3.1 [45], which are given by
(31) | ||||
(32) | ||||
(33) |
Setting and substituting it into (C) and (32), we have
(34) | ||||
(35) |
It is straightforward to verify that the left hand side of (34) is a strictly increasing function with respect to . Hence, there exists a unique solution, denoted by , satisfying (34). From (35), we can observe that is upper-bounded by and can be thus found by the bisection method. Also, (35) indicates that . Having from (33) and , (25) and (26) are obtained with some simple mathematical calculations. This thus proves Proposition 3.
Appendix D Proof of Lemma 1
First of all, for the IRS to be able to assist in downlink ET and uplink IT, we must have
(36) |
according to C16. Otherwise if , the IRS cannot harvest enough energy to power its circuit operations even if the harvested power reaches its maximum value (i.e., saturation power). Furthermore, at optimality, the received power at the energy harvester of the IRS must not be greater than the saturation power, because otherwise, the amount of the reflected power by the IRS can be increased by increasing the amplitude reflection coefficient, without affecting the amount of harvested power at the IRS. Therefore, we must have . Therefore, . Now, according to the energy causality constraint of the IRS in C16, we have . Thus, is upper-bounded by
(37) |
To ensure a feasible value for , the following condition must be satisfied:
(38) |
From (36) and (38) and the fact that , we obtain . Lemma 1 is thus proved. At the optimal solution, the amplitude reflection coefficient must be set to its upper-bound to maximize the amount of reflected power from the IRS. Therefore, according to (37), is calculated as , where according to (36) and (38). This thus proves Proposition 4.
References
- [1] Cisco edge-to-enterprise IoT analytics for electric utilities. Available Online: https://www.cisco.com/c/en/us/solutions/collateral %****␣Final_Paper-accepted.tex␣Line␣825␣****/data-center-virtualization/big-data/solution-overview-c22-740248.html, Feb. 2018.
- [2] X. Lu, P. Wang, D. Niyato, D. I. Kim, and Z. Han, “Wireless networks with RF energy harvesting: A contemporary survey,” IEEE Commun. Surv. Tut., vol. 17, no. 2, pp. 757-789, Secondquarter 2015.
- [3] H. Ju and R. Zhang, “Throughput maximization in wireless powered communication networks,” IEEE Trans. Wireless Commun., vol. 13, no. 1, pp. 418-428, Jan. 2014.
- [4] P. Ramezani and A. Jamalipour, “Toward the evolution of wireless powered communication networks for the future Internet of Things,” IEEE Network, vol. 31, no. 6, pp. 62-69, Nov./Dec. 2017.
- [5] H. Chen, Y. Li, J. L. Rebelatto, B. F. Uchoa-Filho, and B. Vucetic, “Harvest-then-cooperate: Wireless-powered cooperative communications,” IEEE Trans. Signal Process., vol. 63, no. 7, pp. 1700-1711, Apr., 2015.
- [6] H. Ju and R. Zhang, “Uer cooperation in wireless powered communication networks,” in Proc. IEEE GLOBECOM, Austin, TX, USA, Dec. 2014, pp. 1430-1435.
- [7] Y. Zeng, H. Chen, and R. Zhang, “Bidirectional wireless information and power transfer with a helping relay,“ IEEE Commun. Letters, vol. 20, no. 5, pp. 862-865, May 2016.
- [8] C. Zhong, H. A. Suraweera, G. Zheng, I. Krikidis, and Z. Zhang, “Wireless information and power transfer with full duplex relaying,“ IEEE Trans. Commun., vol. 62, no. 10, pp. 3447-3461, Oct. 2014.
- [9] B. Lyu, D. T. Hoang, and Z. Yang, “User Cooperation in wireless-powered backscatter communication networks,” IEEE Wireless Commun. Lett., vol. 8, no. 2, pp. 632-635, Apr. 2019.
- [10] S. H. Kim and D. I. Kim, “Hybrid backscatter communication for wireless-powered heterogeneous networks,” IEEE Trans. Wireless Commmun., vol. 16, no. 10, pp. 6557-6570, Oct. 2017.
- [11] S. Gong, X. Huang, J. Xu, W. Liu, P. Wang, and D. Niyato, “Backscatter relay communications powered by wireless energy beamforming,” IEEE Trans. Commun., vol. 66, no. 7, pp. 3187-3200, July 2018.
- [12] V. Liu, A. Parks, V. Talla, S. Gollakota, D. Wetherall, and J. R. Smith, “Ambient backscatter: Wireless communication out of thin air,” in Proc. SIGCOMM, pp. 39-50, Hong Kong, Aug. 2013.
- [13] Q. Wu and R. Zhang, “Towards smart and reconfigurable environment: Intelligent reflecting surface aided wireless network,” IEEE Commun. Mag., vol. 58, no. 1. pp 106-112, Jan. 2020.
- [14] S. Gong et al. “Towards smart radio environment for wireless communications via intelligent reflecting surfaces: A comprehensive survey,” IEEE Commun. Surveys Tuts.,, doi: 10.1109/COMST.2020.3004197, 2020.
- [15] M. D. Renzo et al., “Smart radio environments empowered by reconfigurable AI meta-surfaces: An idea whosetime has come,” EURASIP J. Wireless Commun. Netw., vol. 129, pp. 1-20, May 2019.
- [16] C. Huang, G. C. Alexandropoulos, A. Zappone, M. Debbah, and C. Yuen, “Energy efficient multi-user MISO communication using low resolution large intelligent surfaces,” in Proc. IEEE GLOBECOM Workshops, Abu Dhabi, United Arab Emirates, 2018, pp. 1-6.
- [17] C. Huang, A. Zappone, G. C. Alexandropoulos, M. Debbah, and C. Yuen, “Reconfigurable intelligent surfaces for energy efficiency in wireless communication,” IEEE Trans. Wireless Commun., vol. 18, no. 8, pp. 4157-4170, Jun. 2019.
- [18] A. Taha, M. Alrabeiah, and A. Alkhateeb, “Enabling large intelligent surfaces with compressive sensing and deep learning,” Available Online: https://arxiv.org/pdf/1904.10136.pdf, Apr. 2019.
- [19] C. Huang, R. Mo and C. Yuen, “Reconfigurable intelligent surface assisted multiuser MISO systems exploiting deep reinforcement learning,” IEEE J. Sel. Areas Commun., vol. 38, no. 8, pp. 1839-1850, Aug. 2020.
- [20] D. Mishra and H. Johansson, “Channel estimation and low-complexity beamforming design for passive intelligent surface assisted MISO wireless energy transfer,” in Proc. IEEE ICASSP, Brighton, UK, May 2019, pp. 4659–4663.
- [21] X. Yu, D. Xu, and R. Schober, “MISO wireless communication systems via intelligent reflecting surfaces,” in Proc. IEEE/CIC ICCC, Changchun, China, Aug. 2019, pp. 735–740.
- [22] Q. Wu and R. Zhang, “Intelligent reflecting surface enhanced wireless network via joint active and passive beamforming,” IEEE Trans. Wireless Commun., vol. 18, no. 11, pp. 5394-5409, Nov. 2019.
- [23] Q. Wu and R. Zhang, “Weighted sum power maximization for intelligent reflecting surface aided SWIPT,” IEEE Wireless Commun. Lett., vol. 9, no. 5, pp. 586-590, May 2020.
- [24] C. Pan, H. Ren, K. Wang, M. Elkashlan, A. Nallanathan, J. Wang, and L. Hanzo “Intelligent reflecting surface aided MIMO broadcasting for simultaneous wireless information and power transfer,” IEEE J. Sel. Area. Commun., vol. 38, no. 8, pp. 1719-1734, Aug. 2020.
- [25] Z. Chu, W. Hao, P. Xiao, and J. Shi, “Intelligent reflect surface aided multi-antenna secure transmission,” IEEE Wireless Commun. Lett., vol. 9, no. 1, pp. 108-112, Jan. 2020.
- [26] C. Guo, Y. Cui, F. Yang, and L. Ding, “Outage probability analysis and minimization in intelligent reflecting surface-assisted MISO systems,” IEEE Commun. Lett., vol. 24, no. 7, 1563-1567, July 2020.
- [27] Y. Zou, Y. Liu, S. Gong, W. Cheng, D. T. Hoang, and D. Niyato, “Joint energy beamforming and optimization for intelligent reflecting surface enhanced communications,” in Proc. WCNC Workshops, Seoul, South Korea, May 2020, pp. 1-6.
- [28] J. Zhao and Y. Liu, “A survey of intelligent reflecting surfaces (IRSs): Towards 6G wireless communication networks,” Available Online: https://arxiv.org/pdf/1907.04789.pdf, Nov. 2019.
- [29] B. Lyu, D. T. Hoang, S. Gong, and Z. Yang, “Intelligent reflecting surface assisted wireless powered communication networks,” in Proc. WCNC Workshops, Seoul, South Korea, May 2020, pp. 1-6.
- [30] Y. Zheng, S. Bi, Y. J. Zhang, Z. Quan, and H. Wang, “Intelligent reflecting surface enhanced user cooperation in wireless powered communication networks,” IEEE Wireless Commun. Lett., vol. 9, no. 6, pp. 901-905, June 2020.
- [31] C. Huang et al., “Holographic MIMO surfaces for 6G wireless networks: Opportunities, challenges, and trends,” IEEE Wireless Commun. Mag., doi:10.1109/MWC.001.1900534, 2020.
- [32] A. A. Nasir, X. Zhou, S. Durrani, and R. A. Kennedy, “Relaying protocols for wireless energy harvesting and information processing,” IEEE Trans. Wireless Commun., vol. 12, no. 7, pp. 3622-3636, Jul. 2013.
- [33] P. N. Alevizos and A. Bletsas, “Sensitive and Nonlinear Far-Field RF Energy Harvesting in Wireless Communications,” IEEE Trans. Wireless Commun., vol. 17, no. 6, pp. 3670-3685, Jun. 2018.
- [34] Y. Dong, M. J. Hossain, and J. Cheng, “Performance of wireless powered amplify and forward relaying over Nakagami-m fading channels with nonlinear energy harvester,” IEEE Commun. Lett., vol. 20, no. 4, pp. 672-675, Apr. 2016.
- [35] A. El Shafie, D. Niyato, and N. Al-Dhahir, “Security of an ordered-based distributive jamming scheme,” IEEE Commun. Lett. vol. 21 no. 1 pp. 72-75 Jan. 2017.
- [36] S. Pejoski, Z. Hadzi-Velkov, and R. Schober, “Optimal power and time allocation for WPCNs with piece-wise linear EH model,” IEEE Wireless Commun. Lett., vol. 7, no. 3, pp. 364-367, June 2018.
- [37] Z. Wang, L. Liu, and S. Cui, “Channel estimation for intelligent reflecting surface assisted multiuser communications,” Proc. WCNC, Seoul, South Korea, May 2020, pp. 1-6.
- [38] Y. Gao et al., “Reconfigurable intelligent surface for MISO systems with proportional rate constraints,” Proc. ICC, Dublin, Ireland, June 2020, pp. 1-7.
- [39] E. Boshkovska et al., “Practical non-linear energy harvesting model and resource allocation for SWIPT systems,” IEEE Commun. Lett., vol. 19, no. 12, pp. 2082-2085, Dec. 2015.
- [40] D. Mishra and G. C. Alexandropoulos, “Transmit precoding and receive power splitting for harvested power maximization in MIMO SWIPT systems,” IEEE Trans. Green Commun. Netw., vol. 2, no. 3, pp. 774-786, Sept. 2018.
- [41] C. Lu, L. Shi, and Y. Ye, “Maximum throughput of TS/PS scheme in an AF relaying network with non-linear energy harvester,” IEEE Access, vol. 6, pp. 26617-26625, 2018.
- [42] H. Yang et al., “Design of resistor-loaded reflectarray elements for both amplitude and phase control,” IEEE Antennas Wireless Propag. Lett., vol. 16, pp. 1159–1162, 2017.
- [43] S. Hu, Z. Wei, Y. Cai, D. W. K. Ng, and J. Yuan, “Sum-rate maximization for multiuser MISO downlink systems with self-sustainable IRS,” Available Online: https://arxiv.org/pdf/2005.11663.pdf, May 2020.
- [44] Z. Q. Luo, W.-K. Ma, A. M.-C. So, Y. Ye, and S. Zhang, “Semidefinite relaxation of quadratic optimization problems,” IEEE Signal Process., vol. 27, no. 3, pp. 20-34, May 2010.
- [45] S. Boyd and L. Vandenberghe, Convex Optimization. Cambridge University Press, 2004.
- [46] Michael Grant et al., “CVX: Matlab software for disciplined convex programming,” Available Online: http://cvxr.com/cvx, September 2013.
- [47] A. M.-C. So, J. Zhang, and Y. Ye, “On approximating complex quadratic optimization problems via semidefinite programming relaxations,” Mathematical Programming, vol. 110, no. 1, pp. 93–110, Jun. 2007.
- [48] V. Arun et al., “RFocus: Practical beamforming for small devices”. Available Online: https://arxiv.org/pdf/1912.07794.pdf, May 2019.