Energy Efficiency Optimization for SWIPT-Enabled IoT Network with Energy Cooperation

As an advanced technology, simultaneous wireless information and power transfer (SWIPT), combined with the internet of things (IoT) devices, can effectively extend the online cycle of the terminal. To cope with the fluctuation of energy harvesting by the hybrid access points (H-AP), the energy cooperation base station is introduced to realize the sharing of renewable energy. In this paper, we study the SWIPT-enabled IoT networks with cooperation. Our goal is to maximize the energy efficiency of the system, and at the same time, we need to meet the energy harvesting constraints, user quality of service (QoS) constraints and transmission power constraints. We jointly solve the power allocation, time switching and energy cooperation problems. Because this problem is a nonlinear programming problem, it is difficult to solve directly, so we use the alternating variable method, the iterative algorithm is used to solve the power allocation and time switching problem, and the matching algorithm is used to solve the energy cooperation problem. Simulation results show that the proposed algorithm has obvious advantages in energy efficiency performance compared with the comparison algorithm. At the same time, it is also proved that the introduction of energy cooperation technology can effectively reduce system energy consumption and improve system energy efficiency.


Introduction
The development of 5G technology has contributed to the rapid spread of the internet of things (IoT). As an important application scenario for 5G, IoT connects a large number of physical objects such as wearable devices, smart home sensors, industrial sensors and agricultural sensors wirelessly to provide various services such as healthcare, smart industry and smart agriculture [1,2]. According to Cisco, the number of global IoT devices in IoT networks will reach 14.7 billion by 2023 [3], which means that the number of wireless devices and traffic demand will grow at an explosive rate, and it is clear that simply using macro base stations in cellular networks to connect this huge number of devices is not enough. To address this problem, dense hybrid access point (H-AP) deployment is seen as a promising technology to meet the quality of service (QoS) needs of devices [4].
The dense deployment of H-APs increases the energy consumption of communication systems, and in the context of a dual carbon strategy, green solutions need to be developed to reduce network-wide energy consumption. To address this challenge, a large number of scholars have explored efficient energy utilization and solutions, mainly in two dimensions: "open source" and "cost reduction". Specifically, the introduction of renewable energy, i.e., "open source", such as solar and wind energy, to reduce energy consumption on the grid. Efficient energy management strategies, i.e., "cost reduction", such as wireless resource management strategies or H-AP sleeping, control the transmitting power of base stations to reduce system energy consumption. split ratio and subcarrier allocation to enhance the rate of IoT networks. In [21], a penalty function algorithm is proposed for energy-carrying communication networks to jointly optimize transmit power and collect energy to reduce system energy consumption. The authors of [22] investigate the problem of maximizing the throughput of IoT devices and propose a Lagrangian-based algorithm for this gradient, which jointly allocates spectrum and transmit power to improve the total system throughput. The authors of [23] address the energy efficiency optimization problem of optimized networks and propose an iterative algorithm based on the Karush-Kuhn-Tucker (KKT) condition to combine the transmit power and power coefficient, and its algorithm improves the energy efficiency performance of the system. In [24], a multi-objective energy-carrying network optimization problem was investigated, and the multi-objective problem was converted into a single-objective problem by using the defined equivalent sum-rate method for solving the problem, and the scheme optimized both system throughput and system energy consumption. In [25], a Dinkelbach-based two-layer iterative algorithm was proposed to jointly optimize the time-switching and power allocation problems. A particle swarm algorithm was proposed in [26] to optimize the rate of the SWIPT-enabled network.
Driven by the carbon-peaking and carbon-neutrality strategic goals, communication systems will move towards a "greener" direction, where balancing the requirements of low power consumption and high speed is crucial. The authors of [7][8][9][10] focus on the optimization of resource allocation with energy harvesting, which addresses the optimization of base station energy consumption and does not consider terminal standby time. The authors of [14] and others demonstrate IoT networks with energy cooperation, low power consumption and good performance. In [20][21][22][23][24][25][26] studied SWIPT networks optimized the on-net operation of terminals without introducing energy cooperation to save the power consumption of the H-AP. The above literature shows that SWIPT and energy cooperation technologies can effectively reduce system power consumption, however, the energy efficiency performance of SWIP-enabled IoT with energy cooperation is still unknown, and as people attach importance to the green network, energy efficiency becomes more and more important. Therefore, this paper aims to maximize the energy efficiency of the system and realize the design of a green communication system. Inspired by the aforementioned literature, we consider IoT networks in using energy cooperation and SWIPT to optimize the energy efficiency of the system as a goal. The main contributions of this paper are summarized as follows: • We consider a downlink transmission model for SWIPT-enabled IoT with energy cooperation. A resource allocation problem is proposed that considers the quality of service (QoS) constraints for users, energy harvesting constraints, and jointly optimizes the power allocation, time switching coefficients and energy cooperation problems to maximize system energy efficiency as the optimization objective. Considering that the optimization problem is a mixed-integer non-linear programming problem that is difficult to solve directly, we consider decomposing the problem into three sub-problems of lower complexity, namely the power allocation, time switching coefficient and energy cooperation problems. We propose a two-stage algorithm for solving the problem • The first-stage algorithm is used to solve the power allocation and time-switching problems. The first-stage algorithm is a two-level iterative algorithm that the power allocation coefficient and the time switching coefficient are separated using the fixed variable method. In the outer layer, the power allocation solution is obtained using the Dinkelbach method iteratively. In the inner layer of the algorithm, the Dinkelbach method is used again to solve for the time switching coefficient under a fixed power allocation. The power allocation and time switching resolution are obtained by several iterations. Finally, in the second stage of the algorithm, the matching theory is used to obtain the resolution of energy cooperation.

•
Our results show that our proposed algorithm has higher energy efficiency compared to the comparison algorithm. The system with SWIPT has good energy efficiency performance and can effectively extend the terminal on-grid operation cycle. In addition, the simulations show that energy cooperation can effectively reduce the energy consumption of the system. This demonstrates the performance benefits of integrating SWIPT and energy cooperation technologies in the IoT network.
The remainder of the paper is structured as follows. Section 2 presents the systembased model and the energy efficiency maximization problem modeling. In Section 3, we propose a two-layer algorithm to find the optimal joint power allocation and timeswitching allocation. Section 4 investigates the use of a matching algorithm to solve the energy cooperation problem. Numerical results are given in Section 5 to prove the theoretical results. Finally, we conclude the paper in Section 6.

Transmission Model
As shown in Figure 1, we consider a SWIPT-enabled IoT system with energy cooperation, consisting of hybrid access points (H-AP) and a smart grid with downlink transmission, where the H-AP is fed by a mix of smart grid and renewable energy sources. The system consists of m H-AP and j terminal. Let m ∈ {1, 2, 3, · · · , M} denote the set of H-APs. Let j ∈ {1, 2, 3, · · · , N} denote the set of terminals. Each terminal contains information decoding and RF energy harvesting circuits. Considering the low cost and easy implementation of the time-switching method circuit, we distinguish between information and energy signals through the time-switching method. It is assumed that the H-AP is capable of serving multiple terminals simultaneously, the terminals associate with the nearest H-AP, the H-AP improves the spectrum efficiency of the system by sharing the entire transmission band, and multiple terminals under a single H-AP service use orthogonal spectrum resources, only one H-AP can be associated with a terminal. It is also assumed that all H-APs and terminals have perfect channel state information (CSI). The signal S m transmitted by the m-th H-AP can be expressed as S m = P jm s m , with E |s m | H = 1, P jm indicates the transmitted power of the m-th H-AP. When the terminal is associated with the m-th H-AP, the received signal at the terminal can be expressed as = G m jm . The signal interference noise ratio (SINR) is expressed as According to Shannon's formula, the transmission rate of terminal j is given by where W is the bandwidth of the system. j  is expressed as the transmission time allocated to the information time slot.

Energy Model
Assuming that all terminals collect energy from the RF signal using time-switching According to Shannon's formula, the transmission rate of terminal j is given by where W is the bandwidth of the system. τ j is expressed as the transmission time allocated to the information time slot.

Energy Model
Assuming that all terminals collect energy from the RF signal using time-switching techniques and perform SWIPT techniques, the terminals split the received signal from the H-AP into two parts: in the first time slot for information transfer and in the second time slot for energy harvesting. Where τ j is expressed as the transmission time allocated to the information time slot and 1 − τ j is expressed as the time portion of the energy harvesting time slot. For the RF harvesting model of the IoT terminal, the paper uses the widely used linear energy harvesting model [25]. The energy harvesting circuit of the terminal is capable of converting the received power signal and the interfering power signal into DC power for use by the terminal. For the presence of multiple H-APs in an IoT system, the RF energy collected by the terminal consists of the superposition of the power emitted by multiple H-APs. The SWIPT-enabled terminals are used for energy harvesting within time slot 1 − τ j . The energy collected by a single terminal is expressed as where 0 < η j < 1, η j denotes the conversion efficiency of the energy harvesting. It is assumed that the conversion efficiency of all the terminals of the system energy harvesting is the same, satisfying η j = η(∀j).
In IoT systems with energy harvesting, each H-AP is equipped with an energy harvesting device. Due to the uneven distribution of renewable energy density and differences in transmitting power, some H-APs do not harvest enough energy to maintain their own standby power. Some H-APs harvest too much energy. To avoid wasting renewable energy, we have introduced energy cooperation technology, which is an important solution for effective energy dispatch. Through aggregators in the smart grid, the excess energy is transferred to the more power-consuming H-APs, effectively increasing the utilization of renewable energy. During the energy cooperation, the renewable energy received by H-AP where T mm represents energy transferred to other H-APs, T m m represents energy received from other H-APs and E RE m is the renewable energy collected by the H-AP. Typically, the power consumption of a conventional wireless communication system is defined as the following linear model [25], and the power consumption of an H-AP in a communication system is defined as the following model where ζ is the power amplification factor of the signal. P C m is the standby power Consumption of the H-AP, including power consumption such as baseband signal processing and cooling system.
In IoT network, the circuit power consumption of SWIPT-enabled terminals is not considered because of their low device power consumption. The total energy consumption of the system is expressed as

Problem Formulation
The constraints in this paper include QoS constraints for terminals, maximum transmit power constraints for H-APs, and maximizing the total energy efficiency of the IoT system under energy harvesting constraints. According to [8], the energy efficiency (EE) of a system is defined as the ratio of the total achievable rate to the total power consumption. The optimization problem can then be expressed as We can observe that the P1 problem is a complex fractional form, and the energy cooperation problem is an integer programming problem, then the P1 problem is a mixedinteger non-linear optimization problem that is difficult to solve directly. We consider reducing the problem to three less complex sub-problems, namely power allocation, time switching and energy cooperation. In addition the two sets of variables for the power allocation P and time switching coefficients τ, which are mutually coupled, as in the scheme in the literature [12], can be solved for the multi-variable problem by considering the remaining variables as constants and how to solve the remaining variables, and finally by an iterative scheme of alternating variables to obtain the resolution of the multi-variables. Because the power allocation P and time switching coefficients τ are coupled, we propose to solve the problem in the first stage using an iterative algorithm. In contrast, the subproblem of energy cooperation is uncoupled from the two sets of variables mentioned above, and for this reason we consider proposing a matching algorithm in the second stage to find the resolution of energy cooperation.

Joint Power Distribution and Time-Switching Control Algorithms
This section may be divided by subheadings. It should provide a concise and precise description of the experimental results, their interpretation, as well as the experimental conclusions that can be drawn.

Power Distribution Problem
In this section, we propose control algorithms for joint power allocation and time switching. From the P1 solution objective formulation, it is clear that optimizing the transmit power and time switching problems is the focus of ensuring terminal QoS and terminal standby time and improving system energy efficiency. Under the subproblem of ground-given energy collaboration. We have adopted the resource allocation scheme in [25], which is applied to the case of a single H-AP. We have expanded the scheme and applied it to multiple H-APs. we develop a two-layer iterative algorithm to solve for the power allocation and time-switching coefficients. Firstly, in the outer iteration, the iterative algorithm is used to solve for the transmit power given a time factor. In the inner layer of the algorithm, the transmit power is fixed and then the time switching coefficient is solved. The specific algorithmic analysis of the solution process is shown below.

A.
Power allocation method under timed switching allocation Given the two sets of variables for the time switching factor τ and energy cooperation T, only one set of variables for the transmit power needs to be solved. The original problem P1 is downscaled to solve a one-dimensional power distribution problem, then the P1 problem is reformulated as As can be seen from the P2 problem, the problem is a fractional objective function, making the problem non-linear and difficult to solve directly. First, we need to determine the non-concave nature of the optimization problem and then adopt the appropriate solution. The first step is to determine the non-concave nature of the constraints. Condition C1 can be converted to As can be seen from the inequality of the transformation deformation, the constraint is that the feasible domain of the resolution is a convex set. Similarly, the analytic feasible region of the constraint C2, C3, C5 is also a convex set. Since G total is fractional, we first prove the concavity of the numerator N ∑ n=1 R j of the objective function. Then, the first order derivative of N ∑ n=1 R j with respect to P jm is expressed as follows: Then, the second order derivative with respect to P jm is expressed as follows: Let 2 , according to the above formula, we get Then, the Hessian matrix with respect to variables expressed as Then, the opposite matrix of the Hessian matrix is Q = −H, then the j-th order principal subformula of the matrix is expressed as According to the formula, we can see that the power and channel is constantly greater than 0, and from the matrix properties can be obtained H j−1 − H j ≥ 0. Any j-order sequential principal subformula of matrix Q, Q j ≥ 0. It is possible to obtain Q ≥ 0, H ≤ 0. R j is a non-negative concave function with respect to the variable P jm . Similarly, the denominator of the objective function is non-negative. In summary, it is shown that the optimization objective function EE(P) = G total is a concave fractional programming problem with respect to the transmit power P.
It can be seen that the objective function is a fractional programming problem, which is difficult to solve directly, and the Dinkelbach algorithm [26] has been widely used with solving non-linear fractional optimization problems. According to the nature of the Dinkelbach algorithm [27], the objective function needs to be transformed into the form of subtracting the numerator from the denominator. According to the Dinkelbach method, we need to introduce a parameter e = EE(P), and the P2 problem is converted into the following form Proposition 1. Assuming that the optimal transmit power of the H-AP is P * , the e * is an optimal resolution of problem P2.1 for which the sufficient conditions are Max: Proof. The proposition is a classical conclusion in generalized dispersion planning, the proof of which has been proved in [28], and the proof process is not described in this paper. The above proposition provides the theoretical support for the transformation of the optimization problem P2.1. Thus, an approximate equivalent problem for problem P2.1 can be obtained by iteration. the Dinkelbach algorithm requires several iterations to obtain the resolution of the problem, where the power allocation and energy efficiency resolution at the t-th iteration are P(t) and e(t), respectively. Needs to be satisfied in the t-th iteration G total (P (t) ) ∼ = 0. The Dinkelbach parameter e is updated by iterations until the convergence condition is met and the iteration is exited. From the above proof it can be seen that G total is a convex optimization problem and its solution function and constraints satisfy the scope of application of the Lagrangian dual method, so we can use the Lagrangian dual method to solve its optimization problem where denote the Lagrange multipliers of the constraints, respectively, and the pairwise function expressions are g(µ, ν, ψ) = maxL(P, µ, ν, ψ) The pairwise optimization problem expression for the problem is In this section, we use a gradient descent-based algorithm to obtain a power allocation solution to the pairwise optimization problem by multiple iterations, and we need to find the first-order derivative of the Lagrangian function [29], whose derivative is Based on the Lagrangian derivative of the dual function, we give an updated formula for the power distribution, expressed as follows: To ensure that the iteration values converge, we update the step size to satisfy We use the subgradient method to update the Lagrange multipliers of the constraints [25,30]. The subgradient formulation of the dual function g(µ, ν, ψ) is given by The Lagrange multiplier update equation is as follows: where δ is the iteration step.

Time Switching Problem
This section analyzes the time-switching scheme in detail. In the inner layer of the algorithm, the transmit power is fixed and then the time switching coefficients are solved by iteration. After fixing the transmit power P, the time switching coefficients τ are a set of unknown variables to be solved, then the problem is reformulated as To simplify the expression of the formula we make It can be seen that the problem is a fractional programming problem with time coefficients, and as above, we also use the Dinkelbach method for the fractional problem, for which the problem is reformulated after the Dinkelbach treatment as It is clear that the problem is a one-time programming problem with a one-time switching coefficient τ, and it is only necessary to prove whether the function is monotonic to find the resolution of the time switching factor. The derivative of the problem P2. 3 is If the derivative is positive, then the function is monotonically increasing, then the time switching coefficient τ for its positive direction of the boundary value, if the negative, the function is monotonically decreasing, then the minimum value of the boundary value for the negative direction. The time switching coefficient discriminant is expressed as follows: We also consider a special scenario. When the terminal is close to the H-AP, the throughput is almost similar. In this case, we can think that the time switching coefficients are the same, and the computational complexity can be reduced. Finally, we obtain the resolution of the power distribution and the time switching coefficients by iteration of the alternating variables method. The two-level iterative algorithm for power allocation and time switching is shown in Algorithm 1. Algorithm 1. Two-layer iteration for power allocation and time switching 1: Input: transmit power P, time switching coefficient τ, Lagrange multiplier µ, ν, ψ, Maximum number of iterations T out , Update step β(t),δ(t), QoS threshold R min , threshold E min number of iterations t = 1, and ε; 2: Output: energy efficiency EE(P, τ) 3: for j = 1 : 1 : N 4: for 1 < t < T out 5: update transmit power P jm (t + 1) according to (21) 6: update Lagrange multiplier µ, ν, ψ according to (23) 7: for 1 < t < T out 8: update time switching coefficient τ according to (28) 9: end for 10: if EE(t + 1)−EE(t + 1) ≤ ε exit loop; 11: break 12: end if 13: end for 14: end for

Energy Cooperation Programmer
In this section, we investigate IoT systems with energy cooperation scenarios. We develop a many-to-many matching algorithm to solve the renewable energy dispatch problem. As can be seen from the system diagram, energy cooperation is done through aggregators in the smart grid, which act as intermediaries between the base stations and the grid, so that the grid operator charges a fee for the energy exchange through the aggregators, but it is lower than the actual fee, due to the fact that the renewable energy is acquired by the base stations [31]. The issue of cost is beyond the scope of this paper. In this paper, we only consider the scheduling of incoming energy from the perspective of energy consumption and quality of service optimization.
When given two sets of variables using the transmit power P and time switching coefficients τ, only one set of variables for the energy cooperation T needs to be solved. The original optimization problem P1 can be reformulated as optimization problem P4, which is formulated as follows: The energy cooperation subproblem is an integer combinatorial optimization problem. In this paper, we consider the use of many-to-many matching theory to solve it. Matching theory is one of the effective solution tools for studying decentralized resource allocation and can transform the resource allocation problem into a simple distributed problem. In existing studies, the many-to-many matching theory has been used to solve optimization problems related to wireless networks [15,32,33], demonstrating that matching theory has the characteristics of fast convergence and stable configuration results.
We adopt the allocation scheme based on the matching theory in [15,32]. We extend the scheme to the IoT network and adjust the preference according to the characteristics of the IoT network. According to the matching theory allocation scheme, the H-APs are first divided into two sets of categories, where the set M + = m + ∈ M E RE m − P W > 0 denotes the set of H-APs that have excess energy while satisfying their own power consumption, both in terms of energy output. M − = m − ∈ M E RE m − P W < 0 indicates a collection of H-APs that do not collect enough energy to sustain their own power consumption. Based on the principle of matching bilateral benefits [15], the utility functions (preference degrees) corresponding to each other between the two types of H-APs are established, and setting a suitable preference degree function can effectively reduce the loss of collected energy and improve the energy utilization of the system. The matching diagram of its two types of collections. Each H-AP in set M + corresponds to all H-APs in set M − and has a preference list corresponding to it. Similarly, the base stations in set M − have their own preference for the H-APs in set M + . The premise of matching is to first match the base station with the corresponding base station according to its preference to complete the energy cooperation.
In the process of energy cooperation, the transmission efficiency is mainly related to the resistance value of the power line, the greater the resistance value, the greater the energy loss, the loss of energy is expressed as follows: where I is the current in the transmission line, R(l) is the total resistance of the power line, R(l) = ρl, ρ is the resistance factor and l is the length of the power line [34]. It can be seen that the lost energy is positively related to the length of the power line. The transmission efficiency α mm from H-AP m + to H-AP m − is expressed as where T mm indicates the renewable energy allocated by H-AP m + to H-AP m − . There is a difference in transmission efficiency α mm between the two H-APs due to the different lengths of the power lines. Set M + prefers to transfer energy to a base station with higher transmission efficiency to reduce energy losses, so the preference of H-APs in Set M + over those in Set M − is expressed as When H-AP m − sends an energy request to the H-APs in Set M + , the H-APs in set M + will select the H-APs with the highest ranking according to the preference ranking in (32) and accept its request, passing the energy to H-APs m − .
H-Aps M − within the set prefer H-APs with more energy remaining in set M + , as this reduces the number of passes and responses from the H-AP and the H-AP is able to obtain energy faster. The preference of the H-APs in set M − for the H-APs in set M + is expressed as H-AP m − is ranked according to the preference of (33) and H-AP m − selects the H-AP with the highest preference in set m + .
Based on the previous analysis, we propose a solution of joint power allocation, time switching and energy cooperation. The specific solution process is summarized in Algorithm 2.

Performance Analysis
This section verifies the effectiveness of the algorithm through simulation. It is assumed that there are 5 H-APs in the network. The cell range of the H-APs is 100 × 100 m 2 , the terminals N = 10 is uniformly distributed within a 10m radius of the H-AP with a terminal RF energy conversion rate of 0.5 [25]. The communication rate threshold is 10 Mbit/s, terminal energy collection thresholds is 0.01 mW. The channel fading model contains Rayleigh fading and path loss, the channel gain is denoted as d −β g, where d −β is the path loss, β = 2, and g is the small-scale fading, generated by the Rayleigh distribution, with a mean difference of 0 and a variance of 1 [19]. The static power consumption is 6 W and the energy transfer efficiency is rand (0.7-0.9) [12]. The joint optimization problem solved by the proposed algorithm following algorithms to compare the performance of each aspect: the rate-maximization (Max-rate) algorithm, which optimizes the power allocation and time switching [25]. The wireless portable energy resource optimization algorithm, which optimizes power allocation and time switching using particle swarm optimization (PSO) algorithms [26]. The energy cooperation part adopts the matching algorithm in this paper. The simulation parameters are shown in Table 1.  Figure 2 shows the iterative convergence of the algorithm. It can be seen from the figure that the algorithm proposed in this paper is the fastest in terms of convergence speed and also obtains the highest energy efficiency performance. As can be seen from the algorithm flow diagram, the algorithm in this paper is given certain initial values in the setting of parameters, such as energy efficiency and firing power, such that the initial values will reduce the number of iterations to some extent. The particle swarm algorithm, on the other hand, searches from a global resolution, so the search resolution requires a large number of iterations to complete, and therefore iterative convergence becomes slower. This suggests that, for this system, our proposed algorithm, has some performance advantage in terms of convergence speed.  Figure 2 shows the iterative convergence of the algorithm. It can be seen from the figure that the algorithm proposed in this paper is the fastest in terms of convergence speed and also obtains the highest energy efficiency performance. As can be seen from the algorithm flow diagram, the algorithm in this paper is given certain initial values in the setting of parameters, such as energy efficiency and firing power, such that the initial values will reduce the number of iterations to some extent. The particle swarm algorithm, on the other hand, searches from a global resolution, so the search resolution requires a large number of iterations to complete, and therefore iterative convergence becomes slower. This suggests that, for this system, our proposed algorithm, has some performance advantage in terms of convergence speed.  Figure 3 shows the curve of the number of terminals versus energy efficiency. As can be seen from the figure, our proposed algorithm achieves higher energy efficiency compared to the PSO algorithm and the maximum rate algorithm. This is because our algorithm achieves a higher resolution of power allocation during the convergence iterations, which effectively suppresses the co-channel interference to its users, which results in a higher throughput for the users, and at the same time, the system consumes less transmit power. The particle swarm algorithm, on the other hand, tends to fall into localized resolution and does not have high search accuracy. The maximum rate algorithm, although able to obtain  Figure 3 shows the curve of the number of terminals versus energy efficiency. As can be seen from the figure, our proposed algorithm achieves higher energy efficiency compared to the PSO algorithm and the maximum rate algorithm. This is because our algorithm achieves a higher resolution of power allocation during the convergence iterations, which effectively suppresses the co-channel interference to its users, which results in a higher throughput for the users, and at the same time, the system consumes less transmit power. The particle swarm algorithm, on the other hand, tends to fall into localized resolution and does not have high search accuracy. The maximum rate algorithm, although able to obtain higher rates, consumes more transmit power, resulting in a less energy efficient system. On the other hand, it can be seen from the figure that the use of energy cooperation techniques can improve the energy efficiency of the system compared to scenarios where no energy cooperation techniques are used, because the excess renewable energy is fully utilized and the consumption of the grid is reduced. As can be seen from the figure, the algorithm of this paper can be effectively applied to a multi-terminal scenario. Figure 4 shows the curve of the effect of the number of H-APs on energy consumption. From the figure, it can be seen that as the number of H-APs increases, the energy consumption of the system also increases, which is due to the increase in static power consumption of the H-APs. The algorithm in this paper and the PSO algorithm both use energy cooperation technology, which makes full use of renewable energy to reduce the system energy consumption. H-APs that do not use energy cooperation will consume more energy. The maximum rate algorithm has the highest energy consumption because, in order to obtain a higher throughput, the transmit power is high, which leads to more energy consumption of the system on the grid. The algorithm proposed in this paper consumes less energy than the PSO algorithm because the algorithm achieves a resolution that is closer to the optimal solution through multiple iterations and consumes less transmit power. As a result, the algorithm is more suitable for multi-H-AP scenarios.
higher rates, consumes more transmit power, resulting in a less energy efficient system. On the other hand, it can be seen from the figure that the use of energy cooperation techniques can improve the energy efficiency of the system compared to scenarios where no energy cooperation techniques are used, because the excess renewable energy is fully utilized and the consumption of the grid is reduced. As can be seen from the figure, the algorithm of this paper can be effectively applied to a multi-terminal scenario.  it can be seen that as the number of H-APs increases, the energy consumption of the system also increases, which is due to the increase in static power consumption of the H-APs. The algorithm in this paper and the PSO algorithm both use energy cooperation technology, which makes full use of renewable energy to reduce the system energy consumption. H-APs that do not use energy cooperation will consume more energy. The maximum rate algorithm has the highest energy consumption because, in order to obtain a higher throughput, the transmit power is high, which leads to more energy consumption of the system on the grid. The algorithm proposed in this paper consumes less energy than the PSO algorithm because the algorithm achieves a resolution that is closer to the optimal solution through multiple iterations and consumes less transmit power. As a result, the algorithm is more suitable for multi-H-AP scenarios.    Figure 4 shows the curve of the effect of the number of H-APs on energy consumption. From the figure, it can be seen that as the number of H-APs increases, the energy consumption of the system also increases, which is due to the increase in static power consumption of the H-APs. The algorithm in this paper and the PSO algorithm both use energy cooperation technology, which makes full use of renewable energy to reduce the system energy consumption. H-APs that do not use energy cooperation will consume more energy. The maximum rate algorithm has the highest energy consumption because, in order to obtain a higher throughput, the transmit power is high, which leads to more energy consumption of the system on the grid. The algorithm proposed in this paper consumes less energy than the PSO algorithm because the algorithm achieves a resolution that is closer to the optimal solution through multiple iterations and consumes less transmit power. As a result, the algorithm is more suitable for multi-H-AP scenarios.   Figure 5 shows the relationship between the number of H-APs and energy efficiency. As can be seen from the figure, the energy efficiency decreases as the number of H-APs increases. This is due to the fact that the standby energy consumption of the system increases as the number of H-APs increases, resulting in a decrease in energy efficiency. On the other hand, the algorithm proposed in this paper outperforms other algorithms in terms of energy efficiency. This is because our proposed algorithm achieves higher throughput and lower energy consumption in the power allocation problem, and therefore higher energy efficiency performance, which indicates that our proposed algorithm can be applied to multi-H-APs scenarios. Figure 6 shows the effect of QoS on the energy efficiency of the system. The graph shows that as the QoS increases, the energy efficiency of the system decreases. This is because, in order to ensure the QoS constraint, the H-APs needs to increase the transmitting power to ensure the QoS demanded by terminals with poor channel conditions, which increases the energy consumption of the system and consequently decreases the energy efficiency. The graph shows that the energy efficiency of the system decreases more slowly when the QoS is increased compared to other algorithms, which indicates that this paper has a higher resolution accuracy in solving the power allocation and therefore achieves a higher energy efficiency than other algorithms.
increases. This is due to the fact that the standby energy consumption of the system in-creases as the number of H-APs increases, resulting in a decrease in energy efficiency. On the other hand, the algorithm proposed in this paper outperforms other algorithms in terms of energy efficiency. This is because our proposed algorithm achieves higher throughput and lower energy consumption in the power allocation problem, and therefore higher energy efficiency performance, which indicates that our proposed algorithm can be applied to multi-H-APs scenarios.  Figure 6 shows the effect of QoS on the energy efficiency of the system. The graph shows that as the QoS increases, the energy efficiency of the system decreases. This is because, in order to ensure the QoS constraint, the H-APs needs to increase the transmitting power to ensure the QoS demanded by terminals with poor channel conditions, which increases the energy consumption of the system and consequently decreases the energy efficiency. The graph shows that the energy efficiency of the system decreases more slowly when the QoS is increased compared to other algorithms, which indicates that this paper has a higher resolution accuracy in solving the power allocation and therefore achieves a higher energy efficiency than other algorithms.    Figure 6 shows the effect of QoS on the energy efficiency of the system. The graph shows that as the QoS increases, the energy efficiency of the system decreases. This is because, in order to ensure the QoS constraint, the H-APs needs to increase the transmitting power to ensure the QoS demanded by terminals with poor channel conditions, which increases the energy consumption of the system and consequently decreases the energy efficiency. The graph shows that the energy efficiency of the system decreases more slowly when the QoS is increased compared to other algorithms, which indicates that this paper has a higher resolution accuracy in solving the power allocation and therefore achieves a higher energy efficiency than other algorithms.   Figure 7 shows the curve of the effect of energy harvesting on the energy efficiency of the system. As can be seen from the figure, the energy efficiency of the system decreases as the energy harvesting constraint increases. This is because, in order to satisfy the energy harvesting constraint, on the one hand, the H-APs needs to transmit more power to satisfy the constraint, which increases the energy consumption of the system and leads to a decrease in energy efficiency. On the other hand, in order to satisfy the energy collection constraint, the time for the energy collection part will increase, while the time for the message decoding will decrease accordingly, which will increase the amount of energy collected by the terminal, but at the same time will decrease the throughput of the terminal. The introduction of the SWIPT technology, which increases the standby time of the terminal, the effect of this technology on the energy efficiency of the system is negligible, is very promising. There is a trade-off between throughput and terminal standby time, and in practical scenarios the values can be set according to the different terminal categories. Figure 8 gives the curve of the effect of the number of H-APs on the energy consumption of the system. As can be seen from the graph, as the number of H-APs increases, the energy consumption of the system also increases. This is because the energy consumption increases as the standby power consumption of the H-APs increases, and although energy harvesting techniques are introduced, they are not yet able to balance the standby energy consumption. In this section, we mainly show the performance comparison between our proposed energy cooperation algorithm and the comparison algorithm DES [15]. It can be seen from the figure that our algorithm performs better due to the impact of the transmission efficiency that we prioritize in setting the matching preference, which, to a certain extent, reduces the energy loss during transmission. On the other hand, the graph shows that the two sets of algorithms that introduce energy cooperation significantly outperform the scenario without energy cooperation in terms of energy consumption. The energy cooperation technique makes full use of renewable energy sources and avoids the waste of excess energy. Therefore, the introduction of energy cooperation techniques has a positive effect on the energy consumption of the communication system. Figure 7 shows the curve of the effect of energy harvesting on the energy efficiency of the system. As can be seen from the figure, the energy efficiency of the system decreases as the energy harvesting constraint increases. This is because, in order to satisfy the energy harvesting constraint, on the one hand, the H-APs needs to transmit more power to satisfy the constraint, which increases the energy consumption of the system and leads to a decrease in energy efficiency. On the other hand, in order to satisfy the energy collection constraint, the time for the energy collection part will increase, while the time for the message decoding will decrease accordingly, which will increase the amount of energy collected by the terminal, but at the same time will decrease the throughput of the terminal. The introduction of the SWIPT technology, which increases the standby time of the terminal, the effect of this technology on the energy efficiency of the system is negligible, is very promising. There is a trade-off between throughput and terminal standby time, and in practical scenarios the values can be set according to the different terminal categories.  Figure 8 gives the curve of the effect of the number of H-APs on the energy consumption of the system. As can be seen from the graph, as the number of H-APs increases, the energy consumption of the system also increases. This is because the energy consumption increases as the standby power consumption of the H-APs increases, and although energy harvesting techniques are introduced, they are not yet able to balance the standby energy consumption. In this section, we mainly show the performance comparison between our proposed energy cooperation algorithm and the comparison algorithm DES [15]. It can be seen from the figure that our algorithm performs better due to the impact of the transmission efficiency that we prioritize in setting the matching preference, which, to a certain extent, reduces the energy loss during transmission. On the other hand, the graph shows that the two sets of algorithms that introduce energy cooperation significantly outperform the scenario without energy cooperation in terms of energy consumption. The energy cooperation technique makes full use of renewable energy sources and avoids the waste of excess energy. Therefore, the introduction of energy cooperation techniques has a positive effect on the energy consumption of the communication system.   Figure 9 shows the effect of the number of terminals on the energy collected. It is clear from the figure that as the number of terminals increases, the amount of energy collected also increases. It is clear that based on the PSO algorithm its obtained energy collection is significantly better than the algorithm proposed in this paper, this is because the particle swarm algorithm has a larger value of power allocation, although the impact on the energy consumption and energy efficiency of the H-APs is negative, this increases the amount of energy collected, which is positive for the standby time of the terminals. The present algorithm, on the other hand, mainly optimizes the energy efficiency of the system and therefore obtains a smaller resolution of the power allocation, which is able to suppress co-channel interference and increase the throughput, but the amount of energy collected by the terminal is then reduced. Additionally, for different systems under the system, it is necessary to weigh the H-AP power consumption and terminal standby time and set the energy collection constraint of the terminal for different needs.  Figure 9 shows the effect of the number of terminals on the energy collected. It is clear from the figure that as the number of terminals increases, the amount of energy collected also increases. It is clear that based on the PSO algorithm its obtained energy collection is significantly better than the algorithm proposed in this paper, this is because the particle swarm algorithm has a larger value of power allocation, although the impact on the energy consumption and energy efficiency of the H-APs is negative, this increases the amount of energy collected, which is positive for the standby time of the terminals. The present algorithm, on the other hand, mainly optimizes the energy efficiency of the system and therefore obtains a smaller resolution of the power allocation, which is able to suppress co-channel interference and increase the throughput, but the amount of energy collected by the terminal is then reduced. Additionally, for different systems under the system, it is necessary to weigh the H-AP power consumption and terminal standby time and set the energy collection constraint of the terminal for different needs.
lection is significantly better than the algorithm proposed in this paper, this is because the particle swarm algorithm has a larger value of power allocation, although the impact on the energy consumption and energy efficiency of the H-APs is negative, this increases the amount of energy collected, which is positive for the standby time of the terminals. The present algorithm, on the other hand, mainly optimizes the energy efficiency of the system and therefore obtains a smaller resolution of the power allocation, which is able to suppress co-channel interference and increase the throughput, but the amount of energy collected by the terminal is then reduced. Additionally, for different systems under the system, it is necessary to weigh the H-AP power consumption and terminal standby time and set the energy collection constraint of the terminal for different needs.

Conclusions
In this paper, we study the SWIPT-enabled IoT network with energy cooperation. We develop a mathematical model with energy efficiency as the optimization objective, while needing to satisfy quality of service and minimum energy harvesting constraints. The problem is nonlinear and difficult to solve directly. We propose an iterative algorithm to solve the problem of power allocation, time switching and energy cooperation. Simulation results show that. Our proposed algorithm outperforms the comparison algorithm in terms of energy efficiency performance. Moreover, this algorithm has good performance for multi H-APs and multi terminal scenarios. In addition, simulation shows that SWIPT technology can effectively extend the operation cycle of the terminal, and the energy cooperation technology can effectively reduce the system energy consumption, which is positive for the development of green communication.
Our algorithm can be extended to other networks with energy collection, especially low-power terminal device networks, such as the current research focus on 5G networks with NOMA, or heterogeneous network systems, where the energy efficiency of the system can be effectively improved using our proposed algorithm for systems with multiple H-APs. In future work, there is still room for improvement. In this paper, we consider the case of having perfect channel conditions; according to the literature [34], the case of imperfect channel conditions information can cause interruptions and rate degradation, for such problems still need further analysis and proposed solutions. On the other hand, with the rapid development of smart grids and energy cooperation technology involving the trading of energy in order to weigh the interests of energy intermediaries and communication operators, further research is needed to solve it.