Energy-Efficient Analog Beamforming for RF-WET with Charging Time Constraint

Abstract

Internet of Things (IoT) sustainability may hinge on radio frequency wireless energy transfer (RF-WET). However, energy-efficient charging strategies are still needed, motivating our work. Specifically, this letter proposes a time division scheme to efficiently charge low-power devices in an IoT network. For this, a multi-antenna power beacon (PB) drives the devices’ energy harvesting circuit to the highest power conversion efficiency point via energy beamforming, thus achieving minimum energy consumption. Herein, we adopt the analog multi-antenna architecture due to its low complexity, cost, and energy consumption. The proposal includes a simple yet accurate model for the transfer characteristic of the energy harvesting circuit, enabling the optimization framework. The results evince the effectiveness of our RF-WET strategy over a benchmark scheme where the PB charges all the IoT devices simultaneously. Furthermore, the performance increases with the number of PB antennas.

I Introduction

RADIO frequency wireless energy transfer (RF-WET) is key for enabling the next generation of Internet of Things (IoT), including massive IoT, extreme edge devices, and zero-energy devices. As charging does not require physical contact when using RF-WET, one can enhance the devices’ mechanical durability, enable waterproof/dustproof designs, and reduce the risk of electric shocks. Moreover, RF-WET enables broadcasting energy to multiple users simultaneously in mobile scenarios, e.g., using unmanned aerial [1, 2] and terrestrial [3, 4] vehicles, and can potentially operate in non-line-of-sight (NLoS) channel conditions. More importantly, RF-WET mitigates the need for battery replacements, reducing network maintenance costs and the risk of environmental pollution [5].

The energy consumption of the power beacon (PB) is a key performance indicator in WET-enabled networks. It is desirable minimizing it to reduce operational costs and the network’s carbon footprint. Indeed, aiming for global decarbonization, the United Nations’ sustainable development goals [6] call for self-sufficient PBs powered by renewable energy [1, 7]. However, the availability of ambient energy sources is often erratic and may limit the PBs’ energy budget.

Fortunately, the PBs can leverage multi-antenna beamforming to dynamically steer the energy in specific spatial directions and increase the total harvested energy without necessarily requiring higher transmission power levels, thus improving the end-to-end energy conversion efficiency of the RF-WET system. In this regard, analog beamforming provides a more affordable solution regarding hardware complexity, cost, and energy consumption compared to other beamforming implementations [8]. This is because, as opposed to its multi-RF chain digital or hybrid counterparts, analog beamforming relies on a single RF chain and analog components for phase shifting the signal. Notice though that this hardware simplicity limits analog beamforming to single-beam transmissions with low spatial resolution at the time.

The most common analog beamforming implementation utilizes a network of phase shifters to adjust the phase of the equal-power antennas’ feeding signals. To improve spatial resolution, one can adjust the signal’s amplitude at the input of the antennas by using, for instance, variable gain amplifiers at the expense of increased hardware complexity and possibly higher power consumption [9]. Electronically steerable parasitic array radiators are also appealing for realizing analog beamforming. This architecture consists of an active antenna fed by a single RF chain and multiple parasitic antennas terminated in tunable analog loads [10]. Notice that this approach leverages the strong mutual coupling among antennas to adjust the radiation pattern; therefore, it is not suitable for arrays with widely spaced antennas. Load-modulated arrays overcome this limitation by feeding all the antennas with a centralized power amplifier while controlling the impedance of the coupling network to adjust the transmission pattern. Notice that no RF chain is required and the processing, including the array’s current distribution, is carried at baseband [11]. Unfortunately, the dynamic variations of the loads cause impedance mismatches to the effective array’s impedance resulting in unwanted power reflections that can deteriorate the system’s performance. Finally, one can resort to switches, instead of phase shifters, to control the set of radiating antennas connected at the power amplifier [12]. Since switches require less control/operation power than phase shifters, this architecture may become more appealing in low-cost/power applications. Despite its hardware simplicity, the switches-based implementation has the poorest spatial resolution compared to the aforementioned implementations.

Meanwhile, notice that the components’ non-linearities of the RF-WET chain hinder transmit beamforming optimization for minimum PB’s energy consumption. Indeed, the RF energy harvesting (EH) circuits have a non-linear transfer characteristic that is a function of the incident RF power, waveform, and modulation [13]. Unfortunately, accounting for all these factors may over-complicate the RF-EH circuit model, which hinders solving the specific WET-related problem. On the other hand, ignoring the non-linearities of RF-EH circuits to facilitate the solution process and providing interesting insights may lead to inaccurate/misleading conclusions.

In this letter, we propose an energy-saving RF-WET strategy in which a multi-antenna PB charges a low-power IoT network using a phase shifters-based analog beamforming. For this purpose, the PB schedules one device at a time, driving the input of the RF-EH circuits to their highest power conversion efficiency point. The main contributions are threefold: i) we propose a novel energy efficiency-oriented approach for modeling the transfer function of the RF-EH circuit in the region of interest, i.e., between the region of maximum conversion efficiency and saturation, which enables our proposed charging optimization procedure; ii) we propose a charging strategy that improves the end-to-end energy conversion efficiency, i.e., the ratio between the total harvested energy by the IoT devices and the energy consumed by the PB; and iii) our results show that, if time division is feasible, it outperforms the strategy of charging all devices simultaneously when they are separated by multiple times the wavelength of the transmit signal.

Notation: Here, we use boldface lowercase letters to denote column vectors and boldface uppercase letters to denote matrices. $\lVert\mathbf{x}\rVert_{p}$ denote the $\ell_{p}$-norm of $\mathbf{x}$, while $(\cdot)^{T}$, $(\cdot)^{H}$, $\mathrm{tr}(\cdot)$, and $|\cdot|$ denote the transpose, Hermitian transpose, trace, and absolute value operations, respectively. Also, $\mathcal{O}(\cdot)$ is the big-O notation, which specifies worst-case complexity, whereas $\mathbb{C}$ denotes the set of complex numbers with imaginary unit $j=\sqrt{-1}$. Finally, $\mathbb{E}[\cdot]$ is the expectation operator.

II System model and problem formulation

Consider the scenario where a multi-antenna PB charges a set $\mathcal{S}$ of single-antenna IoT devices. Herein, $\mathbf{h}_{s}\in\mathbb{C}^{N\times 1}$, where $N$ is the number of transmit antennas, captures the channel coefficients between the PB’s antennas and the $s^{\mathrm{th}}$ device. The system operates over quasi-static fading channels, for which the channel coefficients $\{\mathbf{h}_{s}\}$, assumed known, remain approximately constant over a transmission block with independent realizations between blocks. Moreover, as Fig. 1 illustrates, the PB uses a time division charging strategy in which only one device is scheduled, i.e., dedicatedly charged, at a time. Hence, during a timeslot of duration $\tau_{s}$τたう start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT, the PB transmits a normalized signal $x_{s}$, i.e., $\mathbb{E}[x_{s}^{H}x_{s}]=1$, to charge the corresponding $s^{\mathrm{th}}$ device, whose receive signal is $\mathbf{h}_{s}^{H}\mathbf{w}_{s}x_{s}$ for a given analog precoder $\mathbf{w}_{s}\in\mathbb{C}^{N\times 1}$. Thus, by ignoring the received noise power, the harvested energy is

$$E_{s}=\sum_{\forall s^{{}^{\prime}}\in\mathcal{S}}\tau_{s^{\prime}}g\left(\varpi_{s,s^{\prime}}\right),$$τたう start_POSTSUBSCRIPT italic_s start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT italic_g ( italic_ϖ start_POSTSUBSCRIPT italic_s , italic_s start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ) , (1)

where $\varpi_{s,s^{\prime}}\!=\!\mathbb{E}_{x}[|\mathbf{h}^{H}_{s}\mathbf{w}_{s^{\prime}}x_{s^{\prime}}|^{2}]=|\mathbf{h}_{s}^{H}\mathbf{w}_{s^{\prime}}|^{2}$ is the corresponding incident RF power at the $s^{\mathrm{th}}$ device when transmitting the signal vector $\mathbf{w}_{s^{\prime}}x_{s^{\prime}}$ for charging the device $s^{\prime}$. Moreover, $g:\mathbb{R}^{+}\rightarrow\mathbb{R}^{+}$ is a generic non-linear transfer function of the RF-EH circuit, which we assume only depends on the incident RF power. Several analytical models are commonly adopted in the literature for $g(\cdot)$, such as piece-wise linear, fractional, cubic fractional, and sigmoidal. We refer the reader to [5] and the references therein for additional information about these models.

Our objective is to minimize the PB’s energy consumption while meeting the energy demands $\{E_{s}^{\mathrm{th}}\}$ of the IoT devices per transmission block of duration $\tau_{\mathrm{max}}$τたう start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT. For that, we must find the precoders $\{\mathbf{w}_{s}\}$, where $w_{s,n}$ is the $n^{\mathrm{th}}$ entry of $\mathbf{w}_{s}$, and per-device charging time $\{\tau_{s}\}$τたう start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT } such that the required time to charge all devices does not exceed $\tau_{\mathrm{max}}$τたう start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT. The optimization problem can be formulated as

$$\mathbf{P1:} \underset{\{\mathbf{w}_{s}\},\{\tau_{s}\},\{p_{s}\}}{\mathrm{min.}}\ \sum_{s\in\mathcal{S}}\tau_{s}p_{s}$$ τたう start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT } , { italic_p start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT } end_UNDERACCENT start_ARG roman_min . end_ARG τたう start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT (2a)

s.t. $$E_{s}\geq E_{s}^{\mathrm{th}},\ \forall s\in\cal S,$$ (2b)

$$\sum_{s=1}^{S}\tau_{s}\leq\tau_{\mathrm{max}},$$τたう start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ≤ italic_τたう start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT , (2c)

$$|w_{s,n}|^{2}\!=\!\frac{p_{s}}{N},\ \!\forall s\!\in\!\mathcal{S},n\!=\!1,\!\ldots,N,$$ (2d)

where $p_{s}$ is the transmit power associated with the precoder $\mathbf{w}_{s}$. Notice that (2d) corresponds to the constant modulus constraint of the analog precoders. Notably, neither the objective function (2a) nor the constraints (2b) and (2d) are convex; therefore, $\mathbf{P1}$ is a non-convex optimization problem.

Figure 1: The multi-antenna PB employs a time division strategy to charge the IoT network via analog beamforming.

III Proposed solution

Solving $\mathbf{P1}$ requires an analytical expression for the RF-EH circuit’s transfer function $g(\cdot)$. In this section, we focus on constructing an accurate mapping that captures the necessary information for the purpose of devising an energy-efficient WET strategy. Without loss of generality, we focus on the point-to-point RF-WET scenario, where the PB charges a single and arbitrary device $s$. Hence, we remove the subindex $s^{\prime}$ from $\varpi_{s,s^{\prime}}$ to simplify the notation. Besides, we assume that the transfer characteristics of the EH circuit, as well as relevant key parameters exposed next, have been obtained through direct measurements for a certain operating frequency.

Let $\eta(\varpi_{s})=g(\varpi_{s})/\varpi_{s}$ηいーた ( italic_ϖ start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ) = italic_g ( italic_ϖ start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ) / italic_ϖ start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT denote the power conversion efficiency function of the rectenna. Besides, let $\varpi_{0}$ denote the incident RF power that maximizes the function $\eta(\cdot)$ηいーた ( ⋅ ). To drive the rectenna at such power level, we consider that the PB utilizes the precoder $\sqrt{p_{0}/N}\mathbf{u}_{s}$, where $\mathbf{u}_{s}$ with $|u_{s,n}|=1$, $n=1,\ldots,N$, is the optimal spatial direction for a given $\mathbf{h}_{s}$. Next, let us consider an arbitrary PB’s transmit power $p_{s}$ for which the resulting incident RF power satisfies $\frac{p_{s}}{N}|\mathbf{h}_{s}^{H}\mathbf{u}_{s}|^{2}\neq\varpi_{0}$. Finally, let $\tau_{0}$τたう start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT and $\tau_{s}$τたう start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT denote the required charging times to harvest $E_{s}$ energy units when the transmit powers are $p_{0}$ and $p_{s}$, respectively. Hence, we can state that

$$\tau_{0}\frac{p_{0}}{N}\left|\mathbf{h}_{s}^{H}\mathbf{u}_{s}\right|^{2}\ \eta(\varpi_{0}) =\tau_{s}\frac{p_{s}}{N}\left|\mathbf{h}_{s}^{H}\mathbf{u}_{s}\right|^{2}\eta\left(\frac{p_{s}}{N}\left|\mathbf{h}_{s}^{H}\mathbf{u}_{s}\right|^{2}\right)$$τたう start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT divide start_ARG italic_p start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG italic_N end_ARG | bold_h start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_u start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_ηいーた ( italic_ϖ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) τたう start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT divide start_ARG italic_p start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT end_ARG start_ARG italic_N end_ARG | bold_h start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_u start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_ηいーた ( divide start_ARG italic_p start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT end_ARG start_ARG italic_N end_ARG | bold_h start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_u start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT )

$$\tau_{0}p_{0} <\tau_{s}p_{s},$$τたう start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT τたう start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT , (3)

provided that $\eta(\varpi_{0})$ηいーた ( italic_ϖ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) is the maximum conversion efficiency. Thus, when $\tau_{0}$τたう start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT is feasible, the optimal strategy consists in driving the rectenna to its maximum conversion efficiency point. On the contrary, if $\tau_{0}$τたう start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT violates the maximum allowable charging time, the PB must increase its transmit power to deliver $E_{s}$ energy units within the given deadline. However, exceeding the saturation point $\varpi_{\mathrm{sat}}$ of the rectenna by increasing the incident RF power requires higher PB’s energy consumption and does not provide any additional harvesting power. To show this, let $p_{\mathrm{sat}}$ denote the transmit power that yields $\varpi_{\mathrm{sat}}$ units of incident power, and let $p_{s}$ denote any arbitrary transmit power satisfying $\frac{p_{s}}{N}|\mathbf{h}^{H}_{s}\mathbf{u}_{s}|^{2}>\varpi_{\mathrm{sat}}$. Besides, let $\tau_{\mathrm{sat}}$τたう start_POSTSUBSCRIPT roman_sat end_POSTSUBSCRIPT and $\tau_{s}$τたう start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT denote the required charging times to harvest $E_{s}$ energy units when the transmit powers are $p_{\mathrm{sat}}$ and $p_{s}$. Hence, we can state

$$\tau_{\mathrm{sat}}\frac{p_{\mathrm{sat}}}{N}\left|\mathbf{h}_{s}^{H}\mathbf{u}_{s}\right|^{2}\eta(\varpi_{\mathrm{sat}}) =\tau_{s}\frac{p_{s}}{N}\left|\mathbf{h}_{s}^{H}\mathbf{u}_{s}\right|^{2}\eta(\frac{p_{s}}{N}\left|\mathbf{h}_{s}^{H}\mathbf{u}_{s}\right|^{2}),$$τたう start_POSTSUBSCRIPT roman_sat end_POSTSUBSCRIPT divide start_ARG italic_p start_POSTSUBSCRIPT roman_sat end_POSTSUBSCRIPT end_ARG start_ARG italic_N end_ARG | bold_h start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_u start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_ηいーた ( italic_ϖ start_POSTSUBSCRIPT roman_sat end_POSTSUBSCRIPT ) τたう start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT divide start_ARG italic_p start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT end_ARG start_ARG italic_N end_ARG | bold_h start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_u start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_ηいーた ( divide start_ARG italic_p start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT end_ARG start_ARG italic_N end_ARG | bold_h start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_u start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) ,

$$\tau_{\mathrm{sat}}p_{\mathrm{sat}} <\tau_{s}p_{s},$$τたう start_POSTSUBSCRIPT roman_sat end_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT roman_sat end_POSTSUBSCRIPT τたう start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT , (4)

since $g(p_{s}\left|\mathbf{h}_{s}^{H}\mathbf{u}_{s}\right|^{2})\approx g(\varpi_{\mathrm{sat}})$ for $\frac{p_{s}}{N}|\mathbf{h}^{H}_{s}\mathbf{u}_{s}|^{2}>\varpi_{\mathrm{sat}}$. We summarize this result in the following remark.

Remark 1.

For the point-to-point scenario, driving the input of the RF-EH circuit to its highest power conversion efficiency point yields the optimum PB’s energy beamforming strategy

Hence, it is advisable to confine the input of the RF-EH circuit, and consequently the EH function, within the bounds of $[\varpi_{0},\varpi_{\mathrm{sat}}]$ to achieve this goal.

Herein, we use a second-degree polynomial function, i.e., $\tilde{g}(\varpi_{s})=a_{2}\varpi_{s}^{2}+a_{1}\varpi_{s}+a_{0}$ to fit the transfer function of the Powercast P2110B RF-EH circuit, as Fig. 2 depicts. The measurements depicted in Fig. 2 were extracted from a power conversion efficiency curve in [14] for a maximum input power of around $11~{}$dBm. Besides, we assume that from that point onward the output of the P2110B RF-EH circuit saturates. Notice that $\tilde{g}(\cdot)$ accurately models the region of interest and, as we will show in the next subsections, enables our proposed optimization framework. Hereinafter, we utilize $\tilde{g}(\cdot)$ as our proposed RF-EH model to distinguish from $g(\cdot)$ which models the supported operating range by the RF-EH circuit. In section V, we explain how the fitting of $g(\cdot)$ is done for performing numerical evaluations.

Figure 2: Transfer characteristic of the Powercast P2110B RF-EH circuit operating at $868~{}\text{MHz}$ frequency [14]. Herein, $a_{2}=-1.952$, $a_{1}=0.663$, and $a_{0}=-1.453\times 10^{-5}$ are obtained by standard curve fitting.

III-A Order-agnostic time division

Since the PB establishes a point-to-point connection at each charging timeslot, the optimal per-device analog precoder is $\mathbf{w}_{s}=\sqrt{\frac{p_{s}}{N}}\left[\frac{h_{s,1}}{|h_{s,1}|}\quad\frac{h_{s,2}}{|h_{s,2}|}\quad\cdots\quad\frac{h_{s,N}}{|h_{s,N}|}\right]^{T}$, which gives an incident RF power of $|\mathbf{h}_{s}^{H}\mathbf{w}_{s}|^{2}=\frac{p_{s}}{N}\lVert\mathbf{h}_{s}\rVert_{1}^{2}$. To ease the mathematical tractability of the problem, we assume that the contribution of the signals $\{x_{s^{\prime}}\},\forall s^{\prime}\in\mathcal{S},s^{\prime}\neq s$, on the $s^{\mathrm{th}}$ non-scheduled device is negligible, but we consider it in the next subsection as the criterion for establishing a charging order. Therefore, (1) reduces to

$$E_{s} \approx\tau_{s}g\left(\left|\mathbf{h}_{s}^{H}\mathbf{w}_{s}\right|^{2}\right)=\tau_{s}g\left(\frac{p_{s}}{N}\lVert\mathbf{h}_{s}\rVert^{2}_{1}\right).$$ τたう start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT italic_g ( | bold_h start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_w start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) = italic_τたう start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT italic_g ( divide start_ARG italic_p start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT end_ARG start_ARG italic_N end_ARG ∥ bold_h start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ∥ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ) . (5)

Then, we can use (5) and the quadratic EH function $\tilde{g}(\cdot)$ to modify the constraint (2b) as

$$\tau_{s}\tilde{g}\left(\frac{p_{s}}{N}\lVert\mathbf{h}_{s}\rVert^{2}_{1}\right) \geq E^{\mathrm{th}}_{s},$$τたう start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT over~ start_ARG italic_g end_ARG ( divide start_ARG italic_p start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT end_ARG start_ARG italic_N end_ARG ∥ bold_h start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ∥ start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT )

$$\frac{a_{2}p_{s}^{2}}{N^{2}}\lVert\mathbf{h}_{s}\rVert_{1}^{4}+\frac{a_{1}p_{s}}{N}\lVert\mathbf{h}_{s}\rVert_{1}^{2}+a_{0} \geq E_{s}^{\mathrm{th}}\tau^{-1}_{s},$$ τたう start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ,

$$c_{2}\tau_{s}^{-1}p_{s}^{-1}+c_{1}p_{s}+c_{0}p_{s}^{-1} \leq 1,$$τたう start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT italic_p start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT + italic_c start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT + italic_c start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT (6)

where $c_{2}=\frac{E_{s}^{\mathrm{th}}N}{a_{1}\lVert\mathbf{h}_{s}\rVert_{1}^{2}}$, $c_{1}=\frac{a_{2}\lVert\mathbf{h}_{s}\rVert_{1}^{2}}{a_{1}N}$, and $c_{0}=\frac{a_{0}N}{a_{1}\lVert\mathbf{h}\rVert_{1}^{2}}$ are constants. The following result generalizes this procedure for other PB/devices architectures.

Remark 2.

Notice that when relying on alternative analog beamforming implementations, e.g., [9, 10, 11, 12], and/or employing multi-antenna RF-EH circuits with RF combining, this procedure still holds and only the precoding/combiner optimization is different.

Now, we can reformulate $\mathbf{P1}$ as the following geometric programming problem

$$\mathbf{P2:}\ \underset{\{p_{s}\},\{\tau_{s}\}}{\mathrm{min.}}\ \sum_{s\in\mathcal{S}}\tau_{s}p_{s}$$ τたう start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT } end_UNDERACCENT start_ARG roman_min . end_ARG τたう start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT (7a)

s.t. $$c_{2}\tau_{s}^{-1}p_{s}^{-1}+c_{1}p_{s}+c_{0}p_{s}^{-1}\leq 1,\ \forall s\in\mathcal{S},$$τたう start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT italic_p start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT + italic_c start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT + italic_c start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT italic_p start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ≤ 1 , (7b)

$$w_{0}\leq p_{s}\frac{\lVert\mathbf{h}_{s}\rVert_{1}^{2}}{N}\leq w_{\mathrm{sat}}, \forall s\in\mathcal{S},$$ (7c)

$$\sum_{s=1}^{S}\tau_{s}\leq\tau_{\mathrm{max}},$$τたう start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT ≤ italic_τたう start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT , (7d)

where the additional constraint (7c) guarantees that the incident power remains within the interval $[\varpi_{0},\varpi_{\mathrm{sat}}]$. Although $\mathbf{P2}$ is still non-convex, it can be recast as a convex problem by making the change of variables $\tilde{\tau}_{s}=\log{\tau_{s}}$τたう end_ARG start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT = roman_log italic_τたう start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT and $\tilde{p_{s}}=\log{p_{s}}$ and then apply the function $\log{(\cdot)}$ to the resulting objective and constraints [15]. The resulting optimization problem can be then efficiently solved with accuracy error $\epsilon$ using, e.g., interior-point methods, for which the required number of iterations is in the order of $\mathcal{O}(\sqrt{9S+1}\log{\frac{9S+1}{\epsilon}})$, each requiring $\mathcal{O}\big{(}(40S^{2}+8S)\sqrt{9S+1}\log{\frac{9S+1}{\epsilon}}\big{)}$ arithmetic operations [16].

The following result establishes the feasibility of $\mathbf{P2}$.

Remark 3.

Since the maximum output power of the RF-EH circuit is $g(\varpi_{\mathrm{sat}})$, $\mathbf{P2}$ is feasible when $\sum_{\forall s\in\mathcal{S}}\frac{E_{s}^{\mathrm{th}}}{g(\varpi_{\mathrm{sat}})}\leq\tau_{\mathrm{max}}$τたう start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT. Notice that this result is stricter than the feasibility condition of $\mathbf{P1}$, which is feasible when $\tau_{\mathrm{max}}g(\varpi_{\mathrm{sat}})\geq\mathrm{max}_{s}\ E_{s}^{\mathrm{th}}$τたう start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT italic_g ( italic_ϖ start_POSTSUBSCRIPT roman_sat end_POSTSUBSCRIPT ) ≥ roman_max start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT italic_E start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT roman_th end_POSTSUPERSCRIPT, as it potentially allows driving the EH circuits at the saturation level during the transmission block.

III-B Order-aware time division

The previous time-beamforming allocation strategy may unintentionally power non-scheduled devices, in which case the approximation (5) will no longer hold. Motivated by this, we propose an order-aware time division algorithm in which the power/time allocation corresponding to each device is recomputed after $\mathbf{P2}$ is solved based on the energy they can harvest from others’ time-beamforming allocations. Specifically, the PB schedules at each time slot the device whose timeslot-beamforming allocation is the most polluting for other devices. For instance, during the first timeslot, the PB schedules the $s^{\mathrm{th}}$ device that satisfies

$$s=\underset{s^{\prime\prime}}{\operatorname{arg\,max}}\sum_{\forall s^{\prime}\in\mathcal{S}\backslash s^{\prime\prime}}\tau_{s^{\prime\prime}}g(\varpi_{s^{\prime},s^{\prime\prime}}),$$τたう start_POSTSUBSCRIPT italic_s start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT italic_g ( italic_ϖ start_POSTSUBSCRIPT italic_s start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , italic_s start_POSTSUPERSCRIPT ′ ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ) , (8)

where $\varpi_{s^{\prime},s^{\prime\prime}}$ is the received power at device $s^{\prime}$ when transmitting the signal vector $\mathbf{w}_{s^{\prime\prime}}x_{s^{\prime\prime}}$, and then, we recompute the power allocation for the non-scheduled devices, i.e., $\forall s^{\prime}\in\mathcal{S}\backslash s$, according to

$$p_{s^{\prime}}=\frac{N}{\lVert\mathbf{h}_{s^{\prime}}\rVert_{1}^{2}}g^{-1}\left(\frac{\tilde{E}_{s^{\prime}}^{\mathrm{th}}}{\tau_{s^{\prime}}}\right),$$τたう start_POSTSUBSCRIPT italic_s start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT end_ARG ) , (9)

unless $\frac{p_{s^{\prime}}}{N}\lVert\mathbf{h}_{s^{\prime}}\rVert_{1}^{2}<\varpi_{0}$, in which case

$$\Big{\{}p_{s^{\prime}}=\frac{\varpi_{0}N}{\lVert\mathbf{h}_{s^{\prime}}\rVert_{1}^{2}},\qquad\tau_{s^{\prime}}=\frac{\tilde{E}_{s^{\prime}}^{\mathrm{th}}}{g(\varpi_{0})}\Big{\}},$$τたう start_POSTSUBSCRIPT italic_s start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT = divide start_ARG over~ start_ARG italic_E end_ARG start_POSTSUBSCRIPT italic_s start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT roman_th end_POSTSUPERSCRIPT end_ARG start_ARG italic_g ( italic_ϖ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) end_ARG } , (10)

where $\tilde{E}_{s^{\prime}}^{\mathrm{th}}=E_{s^{\prime}}^{\mathrm{th}}-\tau_{s}g(\varpi_{s^{\prime},s})$τたう start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT italic_g ( italic_ϖ start_POSTSUBSCRIPT italic_s start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , italic_s end_POSTSUBSCRIPT ) is the updated energy requirement for the non-scheduled devices. Similarly, we schedule the remaining devices, discarding previously charged ones from the updates. Finally, Algorithm 1 sketches the full procedure.

Algorithm 1 Order-aware time division

1:  Input: $\{E_{s}^{\mathrm{th}},\mathbf{h}_{s}\}_{\forall s\in\mathcal{S}},\tau_{\mathrm{max}}$τたう start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT

2:  Solve $\mathbf{P2}$

3:  repeat

4:     Schedule the $s^{\mathrm{th}}$ device satisfying (8)

5:     Update: $\{p_{s^{\prime}},\tau_{s^{\prime}},\tilde{E}_{s^{\prime}}^{\mathrm{th}}\}_{\forall s^{\prime}\in\mathcal{S}\backslash s}$τたう start_POSTSUBSCRIPT italic_s start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT , over~ start_ARG italic_E end_ARG start_POSTSUBSCRIPT italic_s start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT roman_th end_POSTSUPERSCRIPT } start_POSTSUBSCRIPT ∀ italic_s start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ∈ caligraphic_S \ italic_s end_POSTSUBSCRIPT using (9)-(10)

6:     Update: $\mathcal{S}\leftarrow\mathcal{S}\backslash\{s\}$

7:  until $|\mathcal{S}|=1$

8:  Schedule the remaining device in $\mathcal{S}$

9:  Output: $\{p_{s},\tau_{s}\}$τたう start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT } for all devices

III-C Charging all-at-once

In this subsection, we consider the benchmark technique allowing the PB to simultaneously charge all the devices. The PB uses a single analog precoder over the charging time. Hence, we can state the resulting optimization problem as

$$\mathbf{P3:}\quad \underset{p,\mathbf{w},\tau}{\mathrm{min.}}\quad p\tau$$ τたう end_UNDERACCENT start_ARG roman_min . end_ARG τたう (11a)

s.t. $$\tau g\left(|\mathbf{h}_{s}^{H}\mathbf{w}|^{2}\right)\geq E_{s}^{\mathrm{th}},\quad \forall s\in\cal S,$$τたう italic_g ( | bold_h start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT bold_w | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) ≥ italic_E start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT roman_th end_POSTSUPERSCRIPT , (11b)

$$\tau\leq\tau_{\mathrm{max}},$$τたう ≤ italic_τたう start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT , (11c)

$$|w_{n}|^{2}=\frac{p}{N},\ n=1,\ldots,N.$$ (11d)

Since the objective (11a) is quasiconcave and the constraints (11b) and (11d) are non-convex, $\mathbf{P3}$ is a non-convex problem. To circumvent this, we resort to a heuristic method in which the optimal spatial direction is obtained from the following max-min weighted incident RF power (non-convex) problem

$$\mathbf{P4:}\quad \underset{\mathbf{w}}{\mathrm{max.}}\quad\underset{s}{\mathrm{min}}\ \frac{|\mathbf{h}_{s}^{H}\mathbf{w}|^{2}}{E^{\mathrm{th}}_{s}}$$ (12a)

such that the coefficients $1/E^{\mathrm{th}}_{s}$ prioritize the devices with higher energy demands. This guarantees near-optimal results for $\mathbf{P3}$. Then, $\mathbf{P4}$ can be recast into its epigraph form as

$$\underset{\mathbf{w}}{\mathrm{max.}}\quad \xi$$ ξくしー (13a)

s.t. $$|\mathbf{h}^{H}_{s}\mathbf{w}|^{2}\geq\xi E^{\mathrm{th}}_{s},\quad\forall s\in\mathcal{S}.$$ξくしー italic_E start_POSTSUPERSCRIPT roman_th end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT , ∀ italic_s ∈ caligraphic_S . (13b)

Next, let us rewrite $|\mathbf{h}^{H}_{s}\mathbf{w}|^{2}=\mathbf{h}^{H}_{s}\mathbf{w}\mathbf{w}^{H}\mathbf{h}_{s}=\mathrm{tr}(\mathbf{W}\mathbf{H}_{s})$, with $\mathbf{W}=\mathbf{w}\mathbf{w}^{H}$ and $\mathbf{H}_{s}=\mathbf{h}_{s}\mathbf{h}^{H}_{s}$ being rank-$1$ matrices, and recast $\mathbf{P4}$ into a semidefinite programming problem

$$\mathbf{P5:}\quad \underset{\xi,\mathbf{W}}{\mathrm{max.}}\quad \xi$$ ξくしー , bold_W end_UNDERACCENT start_ARG roman_max . end_ARG ξくしー (14a)

s.t. $$\mathrm{tr}(\mathbf{W}\mathbf{H}_{s})\geq\xi E^{\mathrm{th}}_{s},\quad\forall s\in\mathcal{S},$$ξくしー italic_E start_POSTSUPERSCRIPT roman_th end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT , ∀ italic_s ∈ caligraphic_S , (14b)

$$\mathbf{W}\succeq 0.$$ (14c)

Here, we relaxed the rank-$1$ requirement of $\mathbf{W}$ to make the problem convex. Then, we take $\bm{\theta}$θしーた, where $\mathbf{w}=\sqrt{\frac{1}{N}}e^{j\bm{\theta}}$θしーた end_POSTSUPERSCRIPT, as the phase of the dominant eigenvector of $\mathbf{W}$, which is optimal when $\mathbf{P5}$ yields a rank-$1$ solution.

Then, from the constraint (11b), we derive the optimal power allocation for a given $\tau$τたう as

$$p(\tau)=\underset{s}{\mathrm{max}}\ g^{-1}\Big{(}\frac{E^{\mathrm{th}}_{s}}{\tau}\Big{)}\frac{1}{|\mathbf{h}^{H}_{s}\mathbf{w}|^{2}},$$τたう ) = underitalic_s start_ARG roman_max end_ARG italic_g start_POSTSUPERSCRIPT - 1 end_POSTSUPERSCRIPT ( divide start_ARG italic_E start_POSTSUPERSCRIPT roman_th end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT end_ARG start_ARG italic_τたう end_ARG ) divide start_ARG 1 end_ARG start_ARG | bold_h start_POSTSUPERSCRIPT italic_H end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT bold_w | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG , (15)

which is the transmit power required for charging the device with the highest ratio between the needed harvesting power and the input RF power when using the power-normalized precoder $\mathbf{w}$. Finally, we perform an exhaustive search to find the pair $(p,\tau)$τたう ) that minimizes the objective function in (11a) by evaluating (15) over the interval $\tau\in[\frac{\mathrm{max}_{s}\ E^{\mathrm{th}}_{s}}{g(\varpi_{\mathrm{sat}})},\tau_{\mathrm{max}}]$τたう ∈ [ divide start_ARG roman_max start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT italic_E start_POSTSUPERSCRIPT roman_th end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT end_ARG start_ARG italic_g ( italic_ϖ start_POSTSUBSCRIPT roman_sat end_POSTSUBSCRIPT ) end_ARG , italic_τたう start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT ]. After computing the optimal power allocation, we update the optimal (analog) precoder as $\mathbf{w}^{*}=\sqrt{p}\mathbf{w}$.

IV Practical considerations

In practice, the resolution of the phase shifters constrains the spatial flexibility of analog beamforming. Hence, our assumption of infinite resolution phase shifters upper bounds the actual performance of both time division and all-at-once schemes when using practical limited-resolution shifters [17]. This limitation motivates exploring alternative analog beamforming implementations in which the array pattern is controlled by switches, variable gain amplifiers, or even the induction of a nearby antenna. Notice that a formal comparison of those architectures requires careful consideration of power consumption and costs of the corresponding components. Besides, notice that scheduling groups of devices instead of using our proposed time division strategy becomes more appealing in massive IoT deployments. Therefore, charging one device at a time as the network grows may not suffice the energy demands due to excessive delays. Moreover, notice that our strategy extends to the multi-antenna RF-EH case when the signal is combined in the RF domain. In other cases, such as when the receiver comprises multiple RF-EH circuits, the extension is not straightforward. Finally, notice that our proposed strategy relies on estimating the channel state information at the beginning of the charging block for computing the optimal power and time allocation. However, depending on the coherence time of the channel, the channel state information of the latter scheduled devices may be inaccurate. In this regard, we can rely on a robust analog beamforming formulation [18] of our problem or, instead, exploit the channel’s statistics [19], which are easier to estimate and vary more slowly. We left the aforementioned concerns as topics for future work.

V Numerical results

Here, we evaluate the PB energy consumption under the proposed time division schemes, referred to as order-agnostic and order-aware for short, against the all-at-once strategy.

Figure 3: Geometry of the simulation setup. The antenna elements are separated by half the wavelength of the transmitted signal in both $z$ and $y$ directions.

As Fig. 3 illustrates, the PB, deployed at the origin, is equipped with a uniform rectangular array with the same number of antennas in both horizontal and vertical dimensions, i.e., $\sqrt{N}$. Meanwhile, the path-loss for the $s^{\mathrm{th}}$ user located at a distance $d_{s}$ from the PB is modeled as $\beta_{s}=10^{-1.6}d_{s}^{-2.7}$βべーた start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT = 10 start_POSTSUPERSCRIPT - 1.6 end_POSTSUPERSCRIPT italic_d start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - 2.7 end_POSTSUPERSCRIPT [20], while the channel is subject to Rician fading, thus

$$\mathbf{h}_{s}=\sqrt{\frac{\kappa\beta_{s}}{1+\kappa}}\mathbf{h}_{s}^{\mathrm{los}}+\sqrt{\frac{\beta_{s}}{2(1+\kappa)}}\mathbf{h}_{s}^{\mathrm{nlos}},$$κかっぱ italic_βべーた start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT end_ARG start_ARG 1 + italic_κかっぱ end_ARG end_ARG bold_h start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT roman_los end_POSTSUPERSCRIPT + square-root start_ARG divide start_ARG italic_βべーた start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT end_ARG start_ARG 2 ( 1 + italic_κかっぱ ) end_ARG end_ARG bold_h start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT start_POSTSUPERSCRIPT roman_nlos end_POSTSUPERSCRIPT , (16)

where $\mathbf{h}_{s}^{\mathrm{los}}$ and $\mathbf{h}_{s}^{\mathrm{nlos}}$ are the line-of-sight (LoS) and the NLoS channel components. Moreover, the LoS factor $\kappa$κかっぱ is set to $10$ and we model the LoS component of the channel as $\mathbf{h}_{s}^{\mathrm{los}}=\mathbf{h}_{s}^{h}\otimes\mathbf{h}_{s}^{v}$, where $\mathbf{h}_{s}^{h}\!=\![1,e^{j\pi\sin{\theta_{s}}\sin{\phi_{s}}},\ldots,e^{j\pi(\sqrt{N}-1)\sin{\theta_{s}}\sin{\phi_{s}}}\!]^{T}\!$πぱい roman_sin italic_θしーた start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT roman_sin italic_ϕ start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT end_POSTSUPERSCRIPT , … , italic_e start_POSTSUPERSCRIPT italic_j italic_πぱい ( square-root start_ARG italic_N end_ARG - 1 ) roman_sin italic_θしーた start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT roman_sin italic_ϕ start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT end_POSTSUPERSCRIPT ] start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT, $\mathbf{h}_{s}^{v}\!=\!\![\!1,e^{j\pi\cos{\theta_{s}}},\ldots,e^{j\pi(\sqrt{N}-1)\cos{\theta_{s}}}\!]^{T}\!$πぱい roman_cos italic_θしーた start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT end_POSTSUPERSCRIPT , … , italic_e start_POSTSUPERSCRIPT italic_j italic_πぱい ( square-root start_ARG italic_N end_ARG - 1 ) roman_cos italic_θしーた start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT end_POSTSUPERSCRIPT ] start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT are the horizontal and vertical LoS components, respectively, for a given elevation angle $\theta_{s}$θしーた start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT and azimuthal angle $\phi_{s}$ [21]. Finally, unless stated otherwise, we consider a setup with three IoT devices respectively located at $\mathbf{v}_{1}=[0,\ -5,\ 0]^{T}$, $\mathbf{v}_{2}=[5,\ 0,\ -5]^{T}$, and $\mathbf{v}_{3}=[10,\ 10,\ 5]^{T}$ (all distances in meters), given Cartesian coordinates, and $E^{\mathrm{th}}_{s}=4~{}\text{mJ},\forall s\in\mathcal{S}$.

Fig. 4 shows the PB’s energy consumption (hereinafter denoted as $E$) vs the maximum charging time $\tau_{\mathrm{max}}$τたう start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT for all strategies. As expected, the all-at-once strategy provides the shortest charging time, which is approximately $\frac{\mathrm{max}_{s}E^{\mathrm{th}}_{s}}{g(\varpi_{\mathrm{sat}})}$, since all devices are charged simultaneously. However, once $\mathbf{P2}$ becomes feasible the order-agnostic outperforms the all-at-once strategy since scheduling one device at a time maximizes its achievable conversion efficiency. Further improvements are obtained by exploiting our order-aware scheduling (at the expense of a slight increase in computational complexity), specially when $N=16$, showing that non-scheduled devices can also harvest energy. However, as the number of antennas increases, the energy beam narrows which reduces the amount of interference at non-scheduled devices. Notice that the energy consumption curves flatten for large values of $\tau_{\mathrm{max}}$τたう start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT. In such a case, all the RF-EH circuits operate at their maximum efficiency input power $\varpi_{0}$ when using the time division strategies. Meanwhile, for the all-at-once strategy, only the worst-performing device according to (15) is guaranteed to operate nearly at maximum conversion efficiency.

Figure 4: PB’s energy consumption vs $\tau_{\mathrm{max}}$τたう start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT for the time division and all-at-once charging strategies and $N\in\{16,64\}$. The vertical dash-dotted lines delimit the feasible regions for both approaches.

Fig. 5 shows how the network deployment, captured by the parameter $\alpha$αあるふぁ, can change the performance gap among all strategies.

Figure 5: PB’s energy consumption vs $\alpha$αあるふぁ for $N=\{16,64\}$, and $\tau_{\mathrm{max}}=8~{}\text{s}$τたう start_POSTSUBSCRIPT roman_max end_POSTSUBSCRIPT = 8 s.

Mathematically, the devices’ coordinates are updated for each value of $\alpha\in[0,1]$αあるふぁ ∈ [ 0 , 1 ] as $(1-\alpha)\mathbf{v}_{s}+\alpha[2.5,\ 2.5,\ 2.5]^{T}$αあるふぁ ) bold_v start_POSTSUBSCRIPT italic_s end_POSTSUBSCRIPT + italic_αあるふぁ [ 2.5 , 2.5 , 2.5 ] start_POSTSUPERSCRIPT italic_T end_POSTSUPERSCRIPT such that all the devices approach the common spatial point $[2.5,\ 2.5,\ 2.5]^{T}$ as $\alpha$αあるふぁ increases. Notice that, increasing the number of antennas reduces the performance gap between the time division strategies as a result of a narrower energy beam. Moreover, when all devices are deployed near each other, around $\alpha=0.9$αあるふぁ = 0.9, the order-aware scheduling and the all-at-once strategies significantly outperform the order-agnostic scheduling approach. The poorer performance of the order-agnostic scheduling is due to the non-negligible harvested energy at the non-scheduled devices, which is not accurate when all devices are deployed close to each other, separated explicitly by sub-wavelength distances, and a single beam efficiently charges all devices simultaneously. Additionally, observe that under such conditions even the order-aware can become sub-optimal with respect to the all-at-once strategy, specially when $N=64$ for which harvesting from non-intentional charging becomes less feasible.

VI Conclusions

We proposed an energy-efficient time division RF-WET strategy to charge an IoT network using analog beamforming. To this end, we demonstrated that for the point-to-point RF-WET scenario fitting the transfer function of the RF-EH circuit in the interval between maximum conversion efficiency and saturation suffices for computing the optimal power/time allocation. Results evinced the outstanding performance of the proposed strategy, which improves as the number of transmit antennas increases. Nevertheless, we acknowledge the limitations of this scheme when the network has urgent energy demands, or the number of devices grows massive. Moreover, our findings show that all-at-once strategy is appealing when the devices are tightly spatially clustered. As a result, we leave the problem of developing a scheduling algorithm combining both strategies’ advantages as future work.

Reproducible research: The simulation results can be reproduced using the Matlab code available at: https://github.com/Osmel-dev/EE-analog-beamforming-WET