Energy Harvesting and Information Transmission Mode Design for Cooperative EH-Abled IoT Applications in beyond 5G Networks

Energy harvesting (EH) technology is considered to be a promising approach to provide enough energy for energy-constrained Internet of (ings (IoT). In this paper, we propose an energy harvesting and information transmission mode for the spectrum sharing system with cooperative EH-abled IoT applications in beyond 5G networks. Different from most existing IoT spectrumsharing research studies, in our system, both primary user (PU) and IoT devices (IDs) collect energy for their information transmission. In addition, for all IDs, they should realize two communication functions: working as relays to help the information transfer process of PU and completing their own information transmission. We analytically derive exact expressions for the throughput of the primary system and IoTsystem and then formulate two objective functions. It is easy to see that power splitting ratio, dynamic EH ratio, power sharing ratio, and relay selection should be optimized to get the best performance for different communication circumstances. Actually, it is a hybrid NP-hard problem to optimize these parameters and traditional algorithms cannot solve it well. (erefore, a novel algorithm-quantum whale optimization algorithm (QWOA) is proposed to obtain the best performance. Simulation results show the good performance of QWOA in different simulation scenarios.


Introduction
Internet of ings (IoT) is an emerging technique that can provide intelligent information exchange opportunity for smart devices in the future communication scenarios and the beyond 5G networks [1,2]. With the development of IoT, energy efficiency and spectrum efficiency have been two essential concerns due to the rapidly growing number of IoT devices [3,4]. Since the technique of IoT spectrum sharing is able to make better use of the spectrum resource, many experts and scholars have a deep study on it then apply to IoT systems [5][6][7][8]. Some IoT spectrum-sharing protocols are proposed in [9][10][11]. In the investigated IoT system, the licensed primary user (PU) can use its own spectrum if it needs at any time, and IoT devices (IDs) can share the PU's spectrum in the case that IDs have no worse effect on the PU. In [10], Khan et al. studied the spectrum-sharing IoT structure. Some future application scenarios, such as smart grid, smart city, and wireless sensor are investigated. An arbitrarily-shaped underlay cognitive network was investigated by Guo et al. [11]. To satisfy interference constraint and meet outage probability of the PU, a cooperative protocol was proposed. In [12], Aslam et al. designed opportunistic spectrum sharing protocol, frame structure, and energy consumption function.
In order to prolong the life of the network, EH is seen as an effective technique which is widely applied in energylimited environment, where the batteries of communication devices are hard to recharge or replace [13,14]. EH technique can supply enough energy for cooperative communication, information transfer, and many other telecommunication works [15][16][17]. To explore potential of EH, many study works have been conducted in recent years [18,19]. To satisfy energy demand for the communication system, a fixed harvesting energy ratio was considered, and the communication system could harvest energy in one time slot for information transfer in [20][21][22]. In [21], Gao et al. investigated an energy harvesting cooperative communication protocol. Network users can collect energy from ambient environment and complete their own information transmission. In [22], relay nodes can harvest energy from wireless signals. Maji et al. presented the mathematical expressions for the system outage probability. e combination of IoT and EH, which is aimed to exert advantage of them, has been a hot topic in recent years [23,24]. It is necessary to investigate the EH-IoT applications because some IoT scenarios are limited by energy constraint and cannot have a good performance [25][26][27]. An EH-abled IoT system which consists of IoT devices and an access point were considered by Lee and Lee [26]. Hybrid energy resources, i.e., wireless signals and on-grid energy resource were combined to realize information transmission. Gurjar et al. [27] proposed an energy transfer approach in an overlay cognitive IoT system. e system throughput, outage probability, and energy efficiency were investigated. Human kinetic energy is also seen as an energy resource by Ju et al. [28]. e power allocation structure, sink selection, and power control blocks were discussed. Research studies [26][27][28] illustrated the effectiveness of EH technology in IoT applications.
For improving quality of service (QoS) and extend network communication range, relay cooperative communication is an essential technique in modern communication which can be combined with spectrum-sharing networks to further improve the performance [29][30][31][32]. Some works technically studied the relay selection problem for cooperative networks [33][34][35]. In [33], to measure performance of the system more comprehensively, Atapattu et al. considered the minimum achievable rate of the system. To further exploit benefits of relay cooperative communication, there are more and more studies on relay selection problem in recent years. For maximizing sum rate for the whole network in a multiuser network, Sharma et al. [35] presented a relay selection protocol for the cooperative ad hoc network. e optimization problem was assigning available relay nodes for different communication nodes. However, common channel interference (CCI) was not considered in [33,35]. To find a relay selection scheme with CCI, Singh et al. [34] presented a new relay selection scheme considering CCI among multiple users. However, the performance of the proposed protocol is not optimal solution and only a suboptimal solution is obtained.
Most existing works as mentioned above mainly discussed about designing a better system of spectrumsharing IoT [8,10,12] or EH-abled IoT [19,27,28]. e research studies about cooperative information exchange in EH-abled IoT with spectrum sharing are few. However, for some IoT distance transmission situations, cooperative relay transmission is necessary and can greatly improve network performance [29,36,37]. Developing the cooperative information transmission mode for the EHabled IoT with spectrum sharing is essential. Besides, considering the problem of existing wireless EH instability (susceptible to the channel gain), it is necessary to design a dynamic EH time slot structure to collect energy (ensuring energy supply). A dynamic EH time slot structure can maintain the stability of system information transmission and EH.
In this paper, we consider a cooperative EH-abled IoT with spectrum sharing (CEH-IoT-SS), where each ID can not only share the licensed spectrum of the PU to transmit its own information but also can help PU's information transmission. For different practical requirements, two objective functions are formulated and the main contribution of this paper can be summarized as follows: (i) We propose a novel CEH-IoT-SS model for existing IoT applications. Different from models in existence, there are multiple IDs and single PU in CEH-IoT-SS and some IDs can work as relays to help PU's information transfer. At the same time, each ID transmits its own information during one time slot. e PU is composed of a primary transmitter (PT) and a primary receiver (PR). e PT collects energy from ambient energy source while each IoT device transmitter (IDT) harvests energy from ambient energy source and a portion of received PT's signal. e energy harvested from PT's signal for each ID can be seen as a transmission energy supplement. (ii) We derive analytical expressions for the throughput of primary network and IDs. Two objective functions are given for different communication situations. Analytical expression indicates that the scheme of power splitting ratio, dynamic EH ratio, power sharing ratio, and relay selection should be optimized. (iii) A new algorithm-quantum whale optimization algorithm (QWOA) is designed to optimize problems which are proposed. Actually, finding the relay selection scheme is a discrete optimization problem. However, the process for obtaining the best power splitting ratio, power sharing ratio, and dynamic EH ratio is solving a continuous optimization problem.
Since it is a hybrid optimization problem, traditional continuous or discrete algorithms cannot get the best solution. erefore, QWOA is proposed to solve it. rough simulation results, we can see QWOA has advantages which traditional algorithms do not have.
Other parts for this paper are structured in the following way. In Section 2, we will give our proposed CEH-IoT-SS model and time slot structure. For Section 3, expressions of throughput for primary network and IDs are derived. For different communication scenarios, two objective functions are formulated. In Section 4, QWOA is proposed to find the best power splitting ratio, dynamic EH ratio, power sharing ratio, and relay selection scheme for different communication scenarios. Simulation results and conclusion will be provided in Sections 5 and 6.

System Model and Time Slot Structure
In this part, a CEH-IoT-SS is considered. As shown in Figure 1, the primary network consists of one PT and one PR. PT tries to deliver information to PR. We consider the situation where PT cannot transmit information to PR directly so PT needs some relays to complete information transfer [27]. M pairs of IDs, which consist of M IDTs and M IoT device receivers (IDRs), implement two communication functions: (1) some IDTs serve as relays to help PT transmit information; (2) each IDT has its own information for transmitting to its corresponding IDR. Both PT and IDTs can harvest energy from ambient energy sources to complete their own information transmission. Besides, IDTs can scavenge energy from a portion of received signal from PT as transmission energy supplement. In the proposed CEH-IoT-SS, energy storage equipment can be a super capacitor and there is no energy rested at the end of each time slot [21,38].
In this model, each information channel follows Rayleigh distribution, and parameters of different Rayleigh distributions are related with distance between two devices. e channel state information (CSI) between two devices is the same during the same time slot. Meanwhile, two channels are mutual independent from each other in one time slot and CSI of the same channel in two different time slots is independent. CSI from PT to the mth IDT (m � 1, 2, . . . , M) is represented as G PT,IDT m , CSI from PT to the mth IDR is represented as G PT,IDR m , CSI from the mth IDT to the mth IDR is represented as G IDT m ,IDR m , and CSI from the mth IDT to PR is represented as G IDT m ,PR .
To simplify representation of the distance between two communication devices, we make definition as follows: the distance between PT and the mth IDT is denoted as d PT,IDT m . Similarly, the distance between the PT and the mth IDR is denoted as d PT,IDR m , the distance between the mth IDT to the mth IDR is denoted as d IDT m ,IDR m , and the distance between the mth IDT and PR is denoted as d IDT m ,PR . e CSI G PT,IDT m , G PT,IDR m , G IDT m ,IDR m , and G IDT m ,PR follow exponentially distribution with parameters IDT m ,PR , respectively. ℓ is the factor of path loss. For any two different exponentially distribution, it is independent from each other. In proposed CEH-IoT-SS, each PT transmits its information to PR by cooperative communication with IDTs, and IDTs also transmit their own information to corresponding IDRs. In addition, both PT and IDTs collect energy for information transmission. Figure 2 shows one time slot of EH and information transmission mode for the proposed CEH-IoT-SS. Each time slot for the CEH-IoT-SS is composed of three phases: (1) in the first phase (0, αT], both PT and IDTs harvest energy from ambient energy source, α is dynamic EH ratio, 0 < α < 1, and T is the duration of one time slot. (2) In the second phase (αT, (1 + α)T/2], PT transmits its information to IDTs and IDRs. One portion of each IDT's received signal is used for cooperative communication. e other portion of received signal is the collected energy by each IDT as the transmission energy supplement. (3) For the third phase ((1 + α)T/2, T], each IDT which is selected as a relay will transmit its own signal and received signal to IDR and PR simultaneously. Other IDTs transmit their own signals to their corresponding IDRs.

Throughput Analysis and Objective
Function Design is part we will derive expressions for the throughput of primary network and IDs in CEH-IoT-SS. For different communication needs, we formulate two objective functions.

EH of PT and IDT in the First Phase.
For the first phase, both PT and IDTs collect energy from ambient energy source. e length of the first phase is αT, so energy stored by PT and IDTs can be shown as where Y PT is EH rate of PT [38,39]. Similarly, the energy stored for the mth IDT can be calculated as E where Y IDT m is EH rate of the mth IDT. e EH process of PT is ended and PT no longer collects energy in one time slot. However, IDTs will still gather energy in the second phase. Since all the harvested energy is used up in one time slot, so transmit power provided by PT can be shown as

Information Transmission of PT in the Second Phase.
In the second phase for one time slot, PT transmits its signal x to IDTs and IDRs, where x is normalized as E|x| 2 � 1. e received signals at the mth(m � 1, 2, . . . , M) IDT and the mth IDR are shown as follows: where n 1 is additive white Gaussian noise (AWGN), in which power is w 1 , and n 2 is AWGN with power w 2 .
On the basis of power splitting approach, the power splitting circuit of the mth IDT will divide the received signal into two parts, � � βy IDT m and ���� � 1 − βy IDT m , which is used for EH and information transfer, respectively. β(0 < β < 1) is the power splitting ratio. Hence, the portion of the received signal for EH can be given by

Wireless Communications and Mobile Computing
Energy stored by the mth IDT at the second phase is calculated as where ρ is the efficiency of EH. Since collected energy from AWGN is small, it can be neglected. erefore, total energy stored by the mth IDT in the first two phases can be shown as e mth IDT utilizes all the harvested energy in the third phase to complete cooperative communication and its own information transfer, so the power at the mth IDT can be given by On the other hand, the portion of the received signal at the mth IDT for information transfer is as follows: where n 3 is AWGN which is caused by the signal conversion from the passband to baseband [40] at the mth IDT and the power of n 3 is w 3 . e received signal at the mth IDT is needed to be normalized, and the power normalization factor ψ m is in the following: e received normalization signal of the mth IDT is

Information Transmission of IDT in the ird Phase.
In the third phase, some IDTs will be selected as relays. If the mth IDT is selected as a relay, the mth IDT transmits the received normalized-signal y IT− Norm IDT m and its own signal z IDT m by superposition coding in the following: where σ(0 < σ < 1) is power sharing ratio, z m is the combined signal which is divided into two parts. e first is used to transmit PT's signal, and the second part z IDT m is used to transmit IDT's own information. If the mth IDT is not selected as a relay, it just needs to transmit its own information to its corresponding IDR: According to different communication functions of IDT, the signal transferred by the mth IDT can be shown as

Wireless Communications and Mobile Computing
erefore, the received signal at PR is given by where b m (m � 1, 2, . . . , M) shows whether the mth IDT is selected as a relay for information transfer of PT. If the mth IDT is selected as a relay, b m � 1. e relay selection scheme of IDTs can be shown as ����� P IDT m z IDT m , and n 4 are IDT signal interference. erefore, the signal-to-interference plus noise ratio (SINR) of the PR is given by According to (9) and (17), the throughput of PR, i.e., the throughput of the primary network can be shown as follows: Similarly, the average throughput of M IDT-IDR transmission pairs in CEH-IoT-SS is expressed as follows: if IDT m is not selected as relay, if IDT m is selected as relay, Wireless Communications and Mobile Computing 7 where w 5 is the AWGN power.
Proof. See Appendix.

Objective Function Design for Different Communication
Requirements. According to the throughput results of the primary network and IDT-IDR transmission pairs, we design two objective functions for improving the data transmission rate in different communication scenarios.

Objective Function of IDT-IDR Transmission Pairs'
Average roughput with Primary Network Constraint. For this communication scenario, we aim to maximize the average throughput of IDT-IDR transmission pairs. At the same time, the performance of the primary network should satisfy requirement. at is to say, throughput of the primary network should be higher than a certain value. According to (18), we can find that the IDT-IDR transmission pairs' throughput is impacted by power splitting ratio, dynamic EH ratio, power sharing ratio, and relay selection scheme. erefore, we focus on optimizing these parameters to maximize the primary network. Under this constraint of the primary network, we aim to maximize the average throughput of IDT-IDR transmission pairs. Combine (18) and (19), objective function can be given by where R require PR is threshold that the throughput of the primary network should achieve.

Objective Function of Sum roughput with Primary
Network Constraint. In this situation, we aim to maximize sum throughput of the whole network with constraint of the primary network. Optimal power splitting ratio, dynamic EH ratio, power sharing ratio, and relay selection scheme should be found and objective function of sum throughput is formulated as follows: To obtain the best throughput of two communication circumstances proposed above is the process for solving the hybrid NP-hard optimization problem [41]. However, traditional intelligent algorithms cannot solve it well because of slow convergence speed and poor convergence accuracy. Hence, we proposed a novel algorithm-QWOA.

Optimization for Throughput Based on Quantum Whale Optimization Algorithm
For the problems which are proposed above, relay selection is a discrete optimization problem. However, getting the best power splitting ratio, dynamic EH ratio, and power sharing ratio is a progress of finding continuous solution. Since domains of them are different and dimensions are complex, traditional intelligent algorithms cannot be easily applied to get the optimal solution. Hence, a new intelligent algorithmquantum whale optimization algorithm (QWOA) is proposed to solve it.

QWOA for Hybrid Optimization Problem.
QWOA combines the advantages of whale optimization algorithm (WOA) [42] and quantum evolution theory of quantum computing [43][44][45]. Since relay selection is a discrete optimization problem and obtaining other parameters is a continuous optimization problem, for variables in different domains, we will design different evolution approaches. For the optimization problems (20) and (21) Both discrete variables and continuous variables will evolve in different ways. e hth(h � 1, 2, . . . , H) quantum whale x t h in the tth iteration is given by where i � 1, 2, . . . , M + 3 and |φ t hi | 2 + |ϕ t hi | 2 � 1. To lower computational complexity and improve evolutionary efficiency of QWOA, we define 0 ≤ φ t hi ≤ 1 and 0 ≤ ϕ t hi ≤ 1. For each φ t hi , φ t hi � �������� 1 − (ϕ t hi ) 2 , each ϕ t hi can be computed by its corresponding φ t hi . Hence, x t h is simply expressed in the following: where i � 1, 2, . . . , M + 3 and 0 ≤ x t hi ≤ 1. Each x t hi is the ith quantum bit of the hth quantum whale x t h .

Wireless Communications and Mobile Computing
For the hth quantum whale x t h , x t hi (i � 1, 2, . . . , M + 3) should be measured to its own definition domain x t hi according to the following rules: where ϑ t hi is a uniform random number which is distributed from zero to one, l i is the ith dimensional variant's lower bound for continuous variables, and u i is the ith dimensional variant's upper bound for the continuous variables. After the process of mapping, we calculate fitness of each quantum whale according to the fitness function. Until the tth iteration, the global optimal quantum whale p t g of the whole population is shown as p t g � [p t g1 , p t g2 , . . . , p t gi , . . . , p t gM , p t g(M+1) , p t g(M+2) , p t g(M+3) ], and its mapping state is For the whole population, each quantum whale is evolved in different ways, which contains encircling prey, searching for prey, and bubble-net attacking. For each evolution mode, we will update their quantum rotation angle and quantum bit with different methods. During each iteration, we first calculate the coefficient vector ] of the hth quantum whale, and the calculation rules can be shown by where ξ 1 is a uniform random number that is distributed from zero to one. ψ � 1 − t/K and K is the maximum iteration number. For the hth quantum whale, we generate a uniform random number ξ 2 that is distributed from zero to one. When ξ 2 < 0.5 and |A t h | < 1 (|•| is the modulus function), the hth quantum whale will encircle the prey and the ith (i � 1, 2, . . . , M + 3) quantum rotation angle and the ith (i � 1, 2, . . . , M + 3) quantum bit of the hth quantum whale are updated as where ξ 3 is a uniform random number that is distributed from zero to one. abs(·) is the absolute value function. When ξ 2 < 0.5 and |A t h | ≥ 1, the hth quantum whale will enter searching for the prey mode. In this mode, the quantum rotation angle is given by where ξ 4 is a uniform random number, a ∈ 1, 2, . . . , H { } is a random integer, and the quantum bit of each quantum whale is updated as (28).
If ξ 2 is no less than 0.5, the quantum whale will take bubble-net attacking behaviour. In this situation, the quantum rotation angle and the quantum bit can be shown as where ξ 5 is a uniform random number that is distributed from zero to one, c 1 , c 2 , and b are constants which determine the radius of bubble-net, and ξ 6 is a uniform random number which is distributed from − 1 to 1.
e quantum whale μ t+1 h (h � 1, 2, . . . , H) is measured to its definition domain μ t+1 h according to (24) and (25). en, the fitness of μ t+1 h is computed by fitness function. If the fitness of μ t+1 h is better than that of x t h , x t+1 h � μ t+1 h ; else,

Wireless Communications and Mobile Computing
x t+1 h � x t h . At last, the best quantum whale of current iteration is used to update the global optimal quantum whale.

4.2.
roughput Optimization Process Based on QWOA. QWOA can optimize problem of maximizing IDT-IDR transmission pairs' average throughput or network sum throughput in CEH-IoT-SS. We convert the maximizing throughput problem into the process of finding global optimal quantum whale. e fitness function is set as , satisfy constraint condition 0, else . Hence, the throughput optimization process based on QWOA is presented in Algorithm 1:

Simulation Results and Analysis
In this part, the performance of proposed QWOA for different communication scenarios is investigated. e simulation results are divided into two parts. In the first part, we compare the optimization results of proposed QWOA and other intelligent algorithms for different communication needs. In fact, since existing algorithms cannot optimize both the continuous problems and discrete problems, we will give the comparison results of QWOA and the circumstances where the classical continuous algorithms are applied to the proposed problems directly. For the second simulation part, we aim to study the impact of different system parameters on the whole system according to different objective functions. In the simulation, PT and PR are located at (0, 1) and (2, 1) on X-Y plane. For the IDT-IDR transmission pairs, all IDTs are generated in a circle and the center of the circle is (1, 1) and the radius of the circle is 0.5, while IDRs are located in a circle in which center and radius are (1, 0) and 0.5, respectively. e factor of path loss ℓ is equal to 3 [40] and ρ � 1. For all nodes, the power of AWGN is the same, i.e., w 1 � w 2 � w 3 � w 4 � w 5 � 10 − 1 . In order to not lose generality, the EH rate of each IDT is the same. Unless otherwise stated, Y PT � 10, Y IDT � Y IDT m � 10 (m � 1, 2, . . . , M), R require PR � 0.3, and the number of IDT-IDR pairs is 10. All simulation results are average of 500 trials.

Comparison Simulation Result of QWOA and Other
Intelligent Algorithms. In this section, we show the performance of QWOA and other intelligent algorithms for the proposed two problems, which are sum throughput and IDT-IDR transmission pairs' average throughput. Figures 3-6 show the simulation results when the objective function is maximizing the average throughput of IDT-IDR transmission pairs, i.e., equation (20). Meanwhile, Figures 7-9 investigate comparison results for maximizing sum throughput, and the objective function is equation (21). In order to intuitively compare the performance of the traditional algorithm and the proposed algorithm, we apply some continuous intelligent algorithms to CEH-IoT-SS. e process for optimizing power splitting ratio, dynamic EH ratio, and power sharing ratio by continuous intelligent algorithms is the same as the course of optimizing continuous variables. However, for optimizing the relay selection scheme, we are going to round continuous variables up to integer for the compared continuous intelligent algorithms. Compared algorithms can be shown as follows: whale optimization algorithm (WOA), differential evolutionary algorithm (DEA), and the proposed QWOA. In all simulation results, the maximal iteration number is 500 and the number of population size H is set as 20. For QWOA, c 1 � 2, c 2 � 2, and b � 1.5. For WOA, the parameter settings can refer to [42]. For DEA, the parameter settings of DEA can refer to [46]. Figures 3 and 7 give the comparison results of QWOA and two classical continuous intelligent algorithms when considering equations (20) and (21) as optimization object, respectively. It is easy to see that QWOA has a faster convergence speed and higher convergence accuracy than WOA and DEA under the same simulation conditions. Since the proposed QWOA uses the thinking of quantum evolution theory and whale hunting characteristics, QWOA is more suitable to solve hybrid optimization problems. For continuous variables and discrete variables, we design different evolution strategies to obtain the solution. e resolving method of traditional WOA is updating individual in the population by some certain equations, and it is easy to get into a local optimum in the early searching stage. From the simulation results, we can see that WOA gets a better solution in the last dozens of iterations. is is because the traditional WOA will process an accurate search due to the searching step length which is generated in a fixed way. However, it is still a local optimum. Unlike WOA, since designed evolution strategies can make full use of the advantages of quantum theory and WOA, the best solution is obtained by QWOA. For DEA, it has disadvantages such as local convergence and randomness and it cannot have a good performance on the proposed problem. To summarize, QWOA can overcome the shortcomings of traditional algorithms and show its advantage in convergence speed and convergence accuracy.
In Figure 4, we investigate the impact of Y IDT on the IDT-IDR transmission pairs' average throughput. Y PT � 10 and Y IDT varies from 2 to 16 in simulation. It is easy to see that as Y IDT increases, IDT-IDR transmission pairs' average throughput becomes larger and larger for QWOA. is is because as Y IDT increases, each IDT can get more energy to store for information transfer. A larger Y IDT can make each IDT use less time to store enough energy. However, due to the fact that WOA and DEA cannot find the best solution, IDT-IDR transmission pairs' average throughput for WOA and DEA is not increasing all the time as Y IDT becomes larger. Simulation results show that QWOA can find the best solution for the whole network and have a better performance compared with WOA and DEA for different Y IDT . Figure 5 considers the circumstance where Y IDT � 10 and Y PT increases from 2 to 16. For different simulation situations, QWOA has the best performance all the time. e average throughput of IDT-IDR transmission pairs becomes larger and larger when Y PT obtains a higher value. It is easy to understand. As Y PT increases, PT can collect more energy for the same time and PT has more energy for transmitting its own information. In addition, each IDT can store more energy from PT as Y PT increases. More energy of IDTs can bring a bigger SINR for IDT-IDR transmission pairs and more easily satisfy the need of the primary network. increases. is performance matches the theory. e power need of the primary network is increasing as R require PR reaches a larger value. Since the power for transmission of the primary network becomes larger, the power that IDT can use becomes less. erefore, the average throughput of IDT-IDR transmission pairs becomes less. Figure 8 considers the sum throughput of the whole system as Y PT increases. Y PT increases from 2 to 16. It is clear that the sum throughput reaches a higher value when Y PT increases. When Y PT becomes larger, the energy collected by PT increases in one time slot. More energy of PT not only helps information transfer of the primary network but also provides more power for IDT-IDR transmission pairs and (1) Input system parameters of CEH-IoT-SS, the fitness function; (2) Initialize the initial population of H quantum whales; (3) t � 1 // the first iteration; (4) Measure each quantum whale to the definition domain according to (24) and (25); (5) Compute fitness of each quantum whale in accordance with the fitness function and find out the global optimal quantum whale p t g ; (6) while t ≤ K (7) Calculate the coefficient vector of each quantum whale (8) if ξ 2 < 0.5 and |A t h | < 1 (9) Update the quantum whale through (27) and (28) (10) else if ξ 2 < 0.5 and |A t h | ≥ 1 (11) Update the quantum whale through (29) and (28) (20) According to the global optimal quantum whale, obtain the power splitting ratio, dynamic EH ratio, power sharing ratio, and relay selection scheme after K iterations; (21) Output: Power splitting ratio, dynamic EH ratio, power sharing ratio, and relay selection scheme.   each IDTs can share more energy to process cooperative communication and its own information transmission. Besides, by comparing with WOA and DEA, we can find that QWOA can get the best simulation results all the time.
In Figure 9, the impact of R require PR on the whole system is presented. For a certain R require PR , QWOA can obtain the best performance compared with WOA and DEA. As R require PR increase, the sum throughput of the whole system becomes smaller and smaller. is is because a higher R require PR will bring more constraint for the primary network and IDT-IDR transmission pairs. e number of solutions for the best power splitting ratio, dynamic EH ratio, power sharing ratio, and the relay selection decreases with a higher R require PR . erefore, the sum throughput becomes smaller with a higher R require PR . From Figures 3 to 9, the proposed QWOA can get the best performance all the time. All the simulation results show the advantages of QWOA in convergence speed and convergence accuracy.

Simulations of Different Parameters Based on QWOA for
CEH-IoT-SS. In this part, the impact of partial different parameters on two proposed objective functions for  CEH-IoT-SS is presented. We show part of simulation results because of the similarity of parameter effects. Figure 10  Due to the fact that a higher Y PT can help PTstore more energy at the same time and IDT can share more energy from PT, the average throughput of IDT-IDR transmission pairs becomes higher and higher when Y PT increases. Meanwhile, the more need of the primary network will bring part of energy from IDT-IDR transmission pairs. Hence, IDT-IDR transmission pairs' average throughput will decrease with a higher R require PR . Figure 11 considers the impact of Y IDT and Y PT on IDT-IDR transmission pairs' average throughput. Y PT varies from 2 to 16. Y IDT is 2, 4, 6, and 8, respectively. From the simulation results, we can easily find that IDT-IDR transmission pairs' average throughput increases larger with a higher Y IDT and Y PT . For a certain period of time, the higher Y IDT and Y PT , the more energy will be stored. More energy will bring a higher power to the primary network and IDT-IDR transmission pairs for information transmission and it can increase the throughput for the certain communication problem. erefore, boosting Y IDT and Y PT can improve the performance of IDT-IDR transmission pairs.
In Figure 12, the impact of different R require PR and Y IDT on the IDT-IDR transmission pairs is investigated. R require PR is increasing from 0.1 to 0.4, while Y IDT is 2, 4, 6, and 8, respectively. For a certain Y IDT , the average throughput of IDT-IDR transmission pairs decreases with R require PR becomes higher and higher. Meanwhile, the average throughput of IDT-IDR transmission pairs increases with a higher Y IDT . Figure 13 presents the simulation results which consider different Y PT and R require PR for the sum throughput of the primary network and IDT-IDR transmission pairs. In the simulation, R require PR is set as 0.2, 0.4, 0.6, and 0.8, respectively. Y PT varies from 2 to 16. It is easy to see that a higher Y PT is conducive to the whole network and a higher R require PR will reduce the throughput of the network. To improve Y PT is a meaningful thing for the whole network. It is obvious that the proposed QWOA can satisfy different throughput requirements of the primary network.

Conclusion
In this paper, we design a novel EH and information transmission mode for CEH-IoT-SS. In the network, IDTs can work as relays for PU's information transmission and finish its own information transfer in one time slot. e expressions for the throughput of the primary network and IDT-IDR transmission pairs are mathematically derived. For different communication scenarios, we formulate two objective functions with different communication demands. Since the formulated objective functions are hybrid optimization problems, traditional intelligent algorithms cannot solve it well. A new algorithm-QWOA is proposed to solve it. Compared with the traditional algorithms, QWOA can get the best solution for hybrid optimization problems and various system parameters. Simulation results validate that QWOA has a wide application for different communication scenarios and hybrid problems. In our future work, we will apply our work to other specific IoT scenarios, such as A-IoT, fog computing, and blockchain.
where n 5 is the AWGN with power w 5 .
���� � P IDT m G IDT m ,IDR m z IDT m is a useful signal part for the mth IDR and the other parts are interference signals. According to (5), the mth IDR has received information x of PT in the