Theoretical Derivation and Optimization Verification of BER for Indoor SWIPT Environments

Symmetrical antenna array is useful for omni bearing beamforming adjustment with multiple receivers. Beam-forming techniques using evolution algorithms have been studied for multi-user resource allocation in simultaneous wireless information and power transfer (SWIPT) systems. In a high-capacity broadband communication system there are many users with wearable devices. A transmitter provides simultaneous wireless information and power to a particular receiver, and the other receivers harvest energy from the radio frequency while being idle. In addition, the ray bounce tracking method is used to estimate the multi-path channel, and the Fourier method is used to perform the time domain conversion. A simple method for reducing the frequency selective effort of the multiple channels using the feed line length instead of the digital phase shifts is proposed. The feed line length and excitation current of the transmitting antennas are adjusted to maximize the energy harvest efficiency under the bit error rate (BER) constraint. We use the time-domain multipath signal to calculate the BER, which includes the inter symbol interference for the wideband system. In addition, we use multi-objective function for optimization. To the best of our knowledge, resource allocation algorithms for this problem have not been reported in the literature. The optimal radiation patterns are synthesized by the asynchronous particle swarm optimization (APSO) and self-adaptive dynamic differential evolution (SADDE) algorithms. Both APSO and SADDE can form good patterns for the receiver for energy harvesting. However, APSO has a faster convergence speed than SADDE.


Introduction
Simultaneous wireless information and power transfer (SWIPT) has been widely investigated recently [1][2][3][4][5][6][7][8][9][10][11][12]. SWIPT does not only provide a high spectral efficiency of wireless systems, but it also provides energy via ambient radio frequency electromagnetic waves. SWIPT is able to support various sensor networks, the Internet of Things (IoT), and cellular systems [13]. In SWIPT systems, signals with information are used for communication and for energy extraction. Energy harvesting can allow for a longer lifetime of energy-constrained wireless devices.
A trade-off between the data rate and energy transfer for a narrow band signal was proposed in the literature [1]. Then, the authors of [2] used a frequency-selective channel with additive Gaussian noise. Another piece of research proposed a dynamic power splitting receiver to spilt the signal and decoded information, with a variable ratio of power to harvest energy [3]. A broadband system with multiplexing of the orthogonal frequency division and beamforming of the transmission was used for SWIPT. [4] However, the direct transmission of the wide band was not considered in the paper (without carrier). The authors of [5] studied a multiple-input multiple-output (MIMO) SWIPT system. The beamforming technique transmitted to the energy and the information receiver of the narrow band signal. Robust beamforming was studied on the assumption of imperfect channel state information at the transmitter in the literature [6]. In the work of [7], the analytic solution was used to optimize the beamforming technology of the SWIPT system. In the literature [10], the allocation of a multi-object resource for secured communication in recognizable radio networks was studied. However, the study only used an analytic solution to calculate the bit error rate (BER) with a statistical channel. In the work of [14], precoding optimization research reduced the resources required for computing, but energy was not taken into consideration.
Most papers deal with the optimal analytic solution with the signal-to-noise ratio constraint for narrow-band signals. However, there are no analytic solutions for BER in the wideband system. For wide-band signals, the weighting vector is assumed to be adjusted according to the different frequencies, and BER should consider the inter symbol interference. In this paper, for design simplicity, we controlled the feed line length and current of the antenna element by beamforming techniques, based on the frequency selective effect of the multipath channel instead of the digital phase shifters. The circular antenna array is useful for adjustment with multiple receivers, as the deployment is symmetrical. Moreover, the bit error rate performance was considered in the wireless wideband communication system. In this paper, a wideband circular array was used to synthesize the array pattern to maximize the energy harvesting efficiency with the BER constraint using evolution algorithms. In Section 2, we explain the model for our wide band system. In Section 3, we introduce the criterion for different objective function optimization. In Section 4, the mechanisms of the algorithms are introduced for optimization. Section 5 includes the discussion on the simulation results, and Section 6 concludes with the key point. The definition of parameters are conclude in Appendix A.

System Description
We considered an indoor SWIPT network for short-distance downlink communication, as shown in Figure 1. We had a base station with antennas of Nt > 1 and receivers of K equipped with one antenna. Both the transmitter and receivers used dipole antennas sharing the same spectrum. The transmitter provided SWIPT services to the receivers. In this system, the transmitter transfers information to one receiver, in which the information receiver divides 50% of the received signal for the harvested power, and the other receivers use all of the received signal as the harvested power. We formulated the object function for energy efficiency (EE) with the BER constraint for data communication. We adjusted the excited current and phase for each transmitter antenna, and maximized the wireless power transfer efficiency using the self-adaptive dynamic differential evolution (SADDE) algorithm and asynchronous particle swarm optimization (APSO).

Transmitter Architecture Design
We exploited eight dipole antennas with ultra-wideband (UWB) to form an array in a circle for the transmitter, as seen in Figure 2.

Transmitter Architecture Design
We exploited eight dipole antennas with ultra-wideband (UWB) to form an array in a circle for the transmitter, as seen in Figure 2. A UWB printed dipole antenna is the individual element for the array. The radiation pattern is omnidirectional in the azimuth plane. The pattern of the array antenna is closely related to the arrangement of the antenna, as well as to the amplitude and the phase of the element of the antenna. With N antennas arranged uniformly in a circular array, the antenna pattern is assumed to be isotropic, and the total field can be expressed as where and represent the spherical coordinate system, a is the distance from original position to antenna position, and and are the spherical coordinates of the transmitting antenna ( = + ). For the far field, the electrical field can be expressed as

Transmitter Architecture Design
We exploited eight dipole antennas with ultra-wideband (UWB) to form an array in a circle for the transmitter, as seen in Figure 2. A UWB printed dipole antenna is the individual element for the array. The radiation pattern is omnidirectional in the azimuth plane. The pattern of the array antenna is closely related to the arrangement of the antenna, as well as to the amplitude and the phase of the element of the antenna. With N antennas arranged uniformly in a circular array, the antenna pattern is assumed to be isotropic, and the total field can be expressed as where and represent the spherical coordinate system, a is the distance from original position to antenna position, and and are the spherical coordinates of the transmitting antenna ( = + ). For the far field, the electrical field can be expressed as (2) Figure 2. Circular ultra-wideband (UWB) dipole antennas.

BS
A UWB printed dipole antenna is the individual element for the array. The radiation pattern is omnidirectional in the azimuth plane. The pattern of the array antenna is closely related to the arrangement of the antenna, as well as to the amplitude and the phase of the element of the antenna. With N antennas arranged uniformly in a circular array, the antenna pattern is assumed to be isotropic, and the total field can be expressed as where θ and φ represent the spherical coordinate system, a is the distance from original position to antenna position, and θ n and φ n are the spherical coordinates of the transmitting antenna (θ n = θ + ψ n ). For the far field, the electrical field can be expressed as A n e − jkr r , where A n is the phase term.
where I n is the excitation current voltage, k = 2π/λ is the wave-number (λ is the wavelength of the sinusoidal wave), and R is the radius of the circle. α n ( f ) is the phase delay in different frequencies, which corresponds to the feed line length, l n . α n ( f ) = − 2π f c √ ε r ·l n , where C is the light speed and ε r is the permittivity of the feed line. We exploited eight antennas at the transmitter, positioned at ϕ n = 0, 45, 90, 135, 180, 225, 270, and 315 degrees in the circle. Equation (4) derives the array factor of the transmitter as the following equation: The total radiation field of the array antenna is the product of the element radiation pattern and the array factor. The vector of the total radiation at the transmitting site is described as follows: Here, the three-dimensional radiation field vector N e (θ, ϕ, f) can be calculated by the structure simulation of the high-frequency.

Impulse Response
A ray-tracing technique can calculate the information on the channel state for different patterns of antenna radiation. It assumes that the transmitter sends out a pincushion of rays, and traces the environments until the rays either reach the receiver or lose enough power to the receiver. The weaker rays' contribution to the received signal is ignorable. Rays are transmitted, reflected, or diffracted to the receiver. The receiver accepts signals in narrow bands near a certain frequency, and transforms them to the similar and complex signals. By using shooting and bouncing rays (SBR)/image techniques [15,16], the following equation induces the channel frequency response where f is the frequency of the sinusoidal wave, i is the path index, θ i is the ith phase shift, b i is the ith receiving magnitude that depends on the radiation vector of the transmitter and receiver antenna, and N IM is the number of images. The frequency response of the ultra-wideband H UWB ( f ) form f 1 to f 2 can be expressed as In the calculation, the location of the image source is calculated by the method of ray tracing, while the channel frequency response is obtained using the height of the equivalent antenna. Then, applying an inverse fast Fourier transform converts the data in the frequency domain into time domain data in the wideband frequency system [17]. The channel frequency response is converted into a discrete time domain function, i.e., the ultra-wideband impulse response h UWB (t) with the Hermitian method and discrete fast inverse Fourier transform.

Receiver Architecture Design
We used power-splitting architecture in our system, as in Figure 3. The received radio frequency (RF) signals are divided for the information receiver and for the RF energy harvester with different power levels. When µ i is the fraction of RF signals for energy harvesting, µ i ∈ [0, 1] is defined as the power-splitting coefficient for receiver i. This means that η is the portion of RF signals utilized for power gathering. Similarly, the following equation calculates the power of RF energy that is harvested at a receiver j (this receiver splits the power from the source): Symmetry 2020, 12, x FOR PEER REVIEW 5 of 13 different power levels. When is the fraction of RF signals for energy harvesting, ∈ [0, 1] is defined as the power-splitting coefficient for receiver i. This means that is the portion of RF signals utilized for power gathering. Similarly, the following equation calculates the power of RF energy that is harvested at a receiver j (this receiver splits the power from the source):

Quality of Service (QoS) Factor for Indoor SWIPT System
The target for the indoor SWIPT network is providing a wireless power transfer to receivers while transferring data to the desired receiver. Thus, we can define the quality of service (QoS) measurement factor in two ways-data quality and energy efficiency.
Here, we set as the channel vectors between the desired receiver and the transmitter, as the channel vectors between the the kth idle receiver and transmitter, and as the corresponding beam-forming vector. The elements of w are as follows: The multi-path channel with an impulse response is where i is the index of the path, is the time delay for the ith path, and is the attenuation of the path. The pulse stream is where Td is the time span of the signal. The pulse-amplitude modulation symbols ∈ {±1} are considered to be independently identically distributed, and the Gaussian waveform can be expressed as where is the standard deviation. The received signal is expressed as The receiver samples the signal with suitably delayed references as

Quality of Service (QoS) Factor for Indoor SWIPT System
The target for the indoor SWIPT network is providing a wireless power transfer to receivers while transferring data to the desired receiver. Thus, we can define the quality of service (QoS) measurement factor in two ways-data quality and energy efficiency.
Here, we set H as the channel vectors between the desired receiver and the transmitter, G k as the channel vectors between the the kth idle receiver and transmitter, and w as the corresponding beam-forming vector. The elements of w are as follows: The multi-path channel with an impulse response is where i is the index of the path, τ i is the time delay for the ith path, and α i is the attenuation of the path. The pulse stream is where T d is the time span of the signal. The pulse-amplitude modulation symbols d ∈ {±1} are considered to be independently identically distributed, and the Gaussian waveform can be expressed as where σ is the standard deviation. The received signal is expressed as The receiver samples the signal with suitably delayed references as where τ 1 is the delayed time with the first revived wave. Here, the correlator output is Z(n) is compared with a threshold between 1 and −1. The variance of the output noise is σ 2 = N 0 2 E st . The probability of error on average for one bit is calculated as The data quality between the transmitter and the desired receiver for the bit error rate (BER) is as follows: where er f c(x) = 2 √ 2 ∞ x e −y 2 dy is the complementary error function and { d n } = {d 0 , d 1 , . . . d N } is the binary sequence. In the SWIPT system, energy is harvested by idle receivers harvested from the RF to extend their lifetime. The efficiency of energy harvesting is important in the system design of networks, and is defined as the ratio of total harvested power at the receivers and the total radiated power by the transmitter in the SWIPT system. The total harvested energy by the K idle receivers is calculated as where * denotes the conjugate transpose. The total energy harvested from all of the receivers is modeled as where we set η and η k as constants, 0 ≤ η, η k ≤ 1 for all k. This provides the conversion efficiency of the RF energy in converting the received radio signal to electrical energy by the desired receiver and K idle receiver. On the other hand, the energy radiated by the transmitter can be expressed as Thus, the efficiency of energy harvesting in the SWIPT system is expressed as

Optimization Problem Formulations
The problem formulation in our indoor SWIPT system aims to maximize the efficiency of energy harvesting and to provide reliable communication in the network. The problem formulation is defined as maximize : η e f f (w), constraint : BER ≤ W BER The system objective in Equation (23) is maximizing the efficiency of the energy harvesting of the system. W BER is the constraint for the system parameter, which denotes the required BER to maintain the data quality [18].

Evolution Algorithms
To get the synthetic radiation pattern, the best combination of feed line lengths needs to be found. This can be done by changing the length of signal feeding line of the element in the circular array. Then, the efficiency of energy harvesting can be maximized for multi-users. We search for the best solution of the radiation pattern using the evolution method that converts the pattern problem into an optimization problem. The evolution method defines the fitness function with the constraint in Equation (23): If the particle does not satisfy the constraint condition, then the particle will be discarded. However, for increasing the convergence speed, we can add some penalty weighting value to the objective function. For instance, the penalty weighting value is set as 2. In other words, the objective value in Equation (24) will multiply by 2 if the constraint is violated. We can check the final objective value to know whether the constraint is satisfied or not. The numerical results show that the final particle satisfies the constraint.

APSO
The PSO method selects the best solution through the calculation of the generation to search for the solution [18][19][20][21]. This is called a synchronous update. As the equation of the particle swarming for the velocity update depends on the group experience, synchronously updating does not allow for the particle to refer the experience of a group in time, causing a slower convergence. Thus, we changed the synchronous update mode to an asynchronous one. The advantage of the method is that the group experience updates with an optimal solution can be found. The search process of the other particles will immediately learn the experience of the last group, thereby accelerating convergence and saving computation time. Detailed steps of the asynchronous particle swarm optimization can be found in the literature [18][19][20][21].

Self-Adaptive Dynamic Differential Evolution
The evolution method of the self-adaptive dynamic differential evolution (SADDE) algorithm is shown in the literature [21][22][23]. The dynamic difference evolution method treats each solution as a parameter vector and selects the difference in the vector of the parameters of each group. Then, the method multiplies it by a weighting factor of inertia and adds it to the optimal solution to generate a perturbation vector. This vector is used for the perturbation of the best solution to avoid in order to be regional extreme. The mating process has a plausible solution. If the new parametrical vector performs better than the optimal vector, then the new vector becomes the optimal vector, so the population in the dynamic differential algorithm is updated in a dynamic form. Please refer to the literature [22][23][24] for details.

Environment Settings
The target for the indoor SWIPT network is to provide wireless power transfer to receivers while transferring data to the desired receiver. We set the frequency range from 4 GHz to 5 GHz with a spacing of ∆ f = 20 MHz for the simulation in this study, using the regulations of the ultra-wideband indoor communication system developed by the Federal Communications Commission (FCC), R = 0.125 m, which corresponds to the adjacent distance between the antennas for the half wavelength of 4 GHz. The QoS measurement factor can be defined in two ways-data quality and energy efficiency. We simulated the proposed resource allocation schemes to evaluate the system performance. We investigated a realistic room environment, as shown in Figure 4. The concrete room length was 10 m, the width was 10 m, and the height was 3 m. We had two wooden bookcases and one working desk in the room, with heights of 2.5 m and 1 m, respectively. The Tx and Rx1, Rx2, Rx3, and Rx4 antennas were 0.5 m in height. The transmitter was a UWB dipole antenna circular array with 16 antennas. The transmitter Tx position was at 7 m, 6 m, and 0.5 m. The receiver antenna was one UWB dipole antenna. The receivers Rx1, Rx2, Rx3, and Rx4 positions were at 3 m, 9.5 m, and 0.5 m; 3 m, 7.5 m, and 0.5 m; 3 m, 5.5 m, and 0.5 m; and 3 m, 3.5 m, and 0.5 m, respectively. The frequency for the UWB was from 4 GHz to 5 GHz. The transmitted power for each transmitter antenna to the noise power ratio at the receiver side was 25 dB. Note that the signal to noise ratio (SNR) at the receiver is dependent on the distance and multipath. The angles from the transmitter to receivers were Rx1: 138.8 • , Rx2: 159.4 • , Rx3: 187.1 • , and Rx4: 212.0 • . Transmitter Tx transferred signals to all of the receivers. Receiver Rx1 split 50% of the signal as power, and used the other 50% signal to decode the data. Rx2 and Rx3 received all of the signals as harvested power. The maximum required BER for the signal at R1 was set as 10 −5 . mmetry 2020, 12, x FOR PEER REVIEW 8 of her 50% signal to decode the data. Rx2 and Rx3 received all of the signals as harvested power. Th aximum required BER for the signal at R1 was set as 10 . With the impulse and the frequency responses of the channel, we calculated the received pow d BER, respectively. As a result of the omnidirectional characteristic of the array element, th nge of the frequency of the channel was from 4 to 5 GHz. The objective function was discussed f aximum harvested power and was searched by evolution algorithms. In our research, the objectiv nction was the inverse of the energy efficiency in Equation (14), and the number of generatio as set to 600. In other words, when the generation reached 600, the algorithm will be terminate e relative parameters for SADDE and APSO are as following: population size is 30 and crossov obability is 0.8. Note that the minimum objective function corresponds to the maximum harveste wer. Figure 5 shows the compared results of the objective values versus the generations for th DDE and APSO algorithms in our simulation. The horizontal axis is the number of generatio d the vertical axis is the fitness value that is equal to the inverse of the harvested power efficienc With the impulse and the frequency responses of the channel, we calculated the received power and BER, respectively. As a result of the omnidirectional characteristic of the array element, the range of the frequency of the channel was from 4 to 5 GHz. The objective function was discussed for maximum harvested power and was searched by evolution algorithms. In our research, the objective function was the inverse of the energy efficiency in Equation (14), and the number of generations was set to 600. In other words, when the generation reached 600, the algorithm will be terminated. The relative parameters for SADDE and APSO are as following: population size is 30 and crossover probability is 0.8. Note that the minimum objective function corresponds to the maximum harvested power. Figure 5 shows the compared results of the objective values versus the generations for the SADDE and APSO algorithms in our simulation. The horizontal axis is the number of generations and the vertical axis is the fitness value that is equal to the inverse of the harvested power efficiency. It is found that the objective value almost converges near 400 generations with SADDE, and near 100 generations with APSO. According to this result, APSO has a faster convergence speed than SADDE. However, the final harvesting power ratio compared to the harvesting power ratio at 100 generation had only 1% improvement with APSO. Figure 5. Objective value generations for the self-adaptive dynamic differential evolution (SADDE) and asynchronous particle swarm optimization (APSO) algorithms. Figures 6a-d and 7a-d show the radiation pattern for Tx to Rx in the x-y plane by SADDE and APSO at 100, 150, 250, and 600 generations. It is seen that the patterns of both algorithms point to all receivers at early generations, and they have similar radiation patterns as a result. In our case, we know that Rx1 uses half of the power received to decode information. To give the best power efficiency, Tx should send just enough power to Rx1 and try to maximize the power transfer to other receivers. Both SADDE and APSO adjusted the lowest amplitude gain for Rx1. For other receivers with SADDE, it was found that the amplitude gain point to Rx2, Rx3, and Rx4 reached near 15 to 13 at 250 generation. However, the distances of Rx2 and Rx3 were closer to the transmitter than Rx4. The amplitude gain point to Rx4 adjusted from near 15 to 12 at 600 generations. With APSO, the amplitude gain point to Rx2 and Rx3 reached 15 at 150 generations, and that to Rx4 adjusted from near 12 to 13 at 150 generations. The improvements of the total power harvesting ratio were about 2.45 by two algorithms.

Simulation Result
(a) (b) Figure 5. Objective value generations for the self-adaptive dynamic differential evolution (SADDE) and asynchronous particle swarm optimization (APSO) algorithms.
Figures 6a-d and 7a-d show the radiation pattern for Tx to Rx in the x-y plane by SADDE and APSO at 100, 150, 250, and 600 generations. It is seen that the patterns of both algorithms point to all receivers at early generations, and they have similar radiation patterns as a result. In our case, we know that Rx1 uses half of the power received to decode information. To give the best power efficiency, Tx should send just enough power to Rx1 and try to maximize the power transfer to other receivers. Both SADDE and APSO adjusted the lowest amplitude gain for Rx1. For other receivers with SADDE, it was found that the amplitude gain point to Rx2, Rx3, and Rx4 reached near 15 to 13 at 250 generation. However, the distances of Rx2 and Rx3 were closer to the transmitter than Rx4. The amplitude gain point to Rx4 adjusted from near 15 to 12 at 600 generations. With APSO, the amplitude gain point to Rx2 and Rx3 reached 15 at 150 generations, and that to Rx4 adjusted from near 12 to 13 at 150 generations. The improvements of the total power harvesting ratio were about 2.45 by two algorithms.
Symmetry 2020, 12, x FOR PEER REVIEW 9 of 13 Figure 5. Objective value generations for the self-adaptive dynamic differential evolution (SADDE) and asynchronous particle swarm optimization (APSO) algorithms.
Figures 6a-d and 7a-d show the radiation pattern for Tx to Rx in the x-y plane by SADDE and APSO at 100, 150, 250, and 600 generations. It is seen that the patterns of both algorithms point to all receivers at early generations, and they have similar radiation patterns as a result. In our case, we know that Rx1 uses half of the power received to decode information. To give the best power efficiency, Tx should send just enough power to Rx1 and try to maximize the power transfer to other receivers. Both SADDE and APSO adjusted the lowest amplitude gain for Rx1. For other receivers with SADDE, it was found that the amplitude gain point to Rx2, Rx3, and Rx4 reached near 15 to 13 at 250 generation. However, the distances of Rx2 and Rx3 were closer to the transmitter than Rx4. The amplitude gain point to Rx4 adjusted from near 15 to 12 at 600 generations. With APSO, the amplitude gain point to Rx2 and Rx3 reached 15 at 150 generations, and that to Rx4 adjusted from near 12 to 13 at 150 generations. The improvements of the total power harvesting ratio were about 2.45 by two algorithms. (c) (d)  Both SADDE and APSO adjust a strong side-lobe at 247°. Figure 8 shows this side-lobe will reflect to Rx4 by the concrete wall at the point of 4.4736, 0.0, and 0.5, with reflection angle of 23°, making a constructive interference for power transfer.  Both SADDE and APSO adjust a strong side-lobe at 247°. Figure 8 shows this side-lobe will reflect to Rx4 by the concrete wall at the point of 4.4736, 0.0, and 0.5, with reflection angle of 23°, making a constructive interference for power transfer. Both SADDE and APSO adjust a strong side-lobe at 247 • . Figure 8 shows this side-lobe will reflect to Rx4 by the concrete wall at the point of 4.4736, 0.0, and 0.5, with reflection angle of 23 • , making a constructive interference for power transfer.

Conclusions
We studied a design of the algorithm for resource allocation in indoor networks with wireless power transfer on the basis of an objective optimization framework. We used wideband circular antenna arrays at the transmitter to maximize the power efficiency performance and keep the data quality. The BER should consider the inter symbol interference in a wideband system; instead, the BER only considered SNR in narrow band system. We employed SADDE and APSO algorithms to adjust the feed length and excitation voltages to synthesize the radiation pattern. Both the SADDE and APSO patterns pointed to all of the receivers at early generations, and had similar radiation patterns in the final result. The improvements in the total power harvesting ratio were about 2.45. Moreover, the results also reveal that the speed of convergence for APSO was faster than that for SADDE. The APSO uses only 100 generation for the final result.

Conflicts of Interest:
The authors declare no conflict of interest.

Appendix A
List of mathematical symbols and abbreviations Proper noun/mathematical symbols Abbreviations/Description Simultaneous wireless information and power transfer SWIPT Internet of Things IoT Multiple-input multiple-output MIMO Energy efficiency EE Self-adaptive dynamic differential evolution SADDE Asynchronous particle swarm optimization APSO Bit error rate BER Ultra-wideband UWB Quality of service QoS

Conclusions
We studied a design of the algorithm for resource allocation in indoor networks with wireless power transfer on the basis of an objective optimization framework. We used wideband circular antenna arrays at the transmitter to maximize the power efficiency performance and keep the data quality. The BER should consider the inter symbol interference in a wideband system; instead, the BER only considered SNR in narrow band system. We employed SADDE and APSO algorithms to adjust the feed length and excitation voltages to synthesize the radiation pattern. Both the SADDE and APSO patterns pointed to all of the receivers at early generations, and had similar radiation patterns in the final result. The improvements in the total power harvesting ratio were about 2.45. Moreover, the results also reveal that the speed of convergence for APSO was faster than that for SADDE. The APSO uses only 100 generation for the final result.

Conflicts of Interest:
The authors declare no conflict of interest.

Appendix A
List of mathematical symbols and abbreviations Proper noun/mathematical symbols Abbreviations/Description Simultaneous wireless information and power transfer SWIPT Internet of Things IoT Multiple-input multiple-output MIMO Energy efficiency EE Self-adaptive dynamic differential evolution SADDE Asynchronous particle swarm optimization APSO Bit error rate BER Ultra-wideband UWB Quality of service QoS Federal Communications Commission FCC E tot (r, θ, φ) Electrical field θ The spherical coordinate system φ The spherical coordinate system θ n The spherical coordinates of the transmitting antenna φ n The spherical coordinates of the transmitting antenna A n The phase term AF (θ, ϕ, f ) The array factor I n The excitation current voltage λ The wavelength k The wave-number R the radius of circle array ε r permittivity l n the feed line length N(θ, ϕ, f ) the three-dimensional radiation field vector H( f ) the channel frequency response b i the ith receiving magnitude f the frequency of sinusoidal wave θ i the ith phase shift f 1 Start frequency The frequency response of the ultra-wideband P i, j The received RF signals power µ i The fraction of RF signals η The portion of RF signals G k The channel vectors between the transmitter and the kth idle receiver W n The corresponding beam-forming vector h UWB (t) The multi-path channel with an impulse response x(t) The pulse stream τ i The time delay for the ith path p(t) The Gaussian waveform r(t) The received signal q(t) The receiver samples the signal with suitably delayed references Z(n) The output of the correlator σ The variance of the output noise α i The attenuation of the pat T d The duration of the signal d The pulse-amplitude modulation symbols E k The total amount of energy harvested by the K idle receivers HP(w) total The total energy harvested TP(w) The energy radiated by the transmitter η e f f (w) The energy harvesting efficiency W BER The constraint for the system parameter