#### Abstract

In an Energy Harvesting system (EHS) the gamma process is used to model the electromagnetic energy received from radiofrequency (RF) radiation. The stochastic characterization of the harvested energy as a continuous-time stochastic process, namely, gamma process, is obtained from the Nakagami-m fading model, which describes the signal reception in a large amount of types of radiofrequency channels. Using the gamma process, some performance measures of the EHS system are obtained. Also, a transmission policy subject to different fading conditions is considered.

#### 1. Introduction

The energy present in the environment, from sources as wind, solar energy, vibration, and electromagnetic radiofrequency (RF), can be harvested to power systems as low and ultralow power electronics and sensor networks [1]. Nowadays, Energy Harvesting Systems (EHS) are emerging as a promising technology for charging batteries used to supply ultra low power micro controllers, sensors, and wearable systems [2].

In an EHS, the energy is incoming from an exogenous source and is managed using two points of view(i)Harvest-store-use (HSU) where a battery stores the harvested energy before its future use. It requires an energy storage device or battery along with an appropriate management of the stored energy.(ii)Harvest-use (HU) where the harvested energy is not stored and it is directly consumed.

Many works in the literature deal with the HSU architectures due to their interest in commercial applications. Although HU architecture literature is not very common, it is a suitable technology for limited battery systems. For a HSU system, two approaches have been developed taking into account the capacity of storage of the battery: infinite and finite capacity. Assuming infinite capacity means that the harvested energy can always be stored in the battery. This infinity assumption is valid for devices with large storage capacity compared to the incoming energy [3]. On the other hand, for a device with finite capacity, the harvested energy can be stored up to a limit. Beyond this limit, incoming energy cannot be stored and would be lost [4].

From the system theoretical perspective, the harvested energy can be modeled as a deterministic or a stochastic process. In a deterministic model, the incoming energy and its fluctuation are known in advance. These models are suitable for energy sources with predictable or with low fluctuation power sources. For example, assuming no cloud coverage, the sun can be modeled as a deterministic source of energy.

In a RF energy harvesting system, the energy is captured from the electromagnetic spectrum, from sources as radio FM, wireless radio (for example, Wi-Fi communication systems, mobile communications systems), TV broadcasting, or cognitive wireless sensor systems [5, 6]. One of the advantages of the RF harvesting is that the energy is ubiquitous and permanently available; in contrast, energy levels are low, except in the surroundings of the transmitter, because the power of the electromagnetic wave is in inverse proportion to the square of the distance to the transmitter [7]. Other kinds of ambient sources for harvesting as mechanical, thermal, or acoustics are nor ubiquitous but with levels of energy greater than electromagnetic in many cases [8].

In practice, the amount of energy that can be harvested from the environment is generally a random quantity [9]. For example, when RF harvested energy is used, the harvested energy is randomly distributed since it depends on the channel fading [10]. Channel fading is a random adverse phenomenon that affects the signal and energy reception due to multipath and shadowing. Some authors have characterized the RF power received by an energy harvester node assuming a specific distribution for the fading. For example, if the fading follows a Rayleigh distribution and assuming Gamma shadowing, and the RF received power a the antenna can be modelled using a sum of gamma distributions [11].

In this paper, we focus on the distribution probability for the energy harvested. In the literature, the stochastic models used to describe the energy harvested by the system include Poisson processes [12], Markov models [13], or compound Poisson processes [14, 15].

For a compound Poisson process, energy arrivals occur at random times (following a Poisson process) and in random amounts. Using a compound Poisson process model, the total harvested energy at time can be expressed as where represents the number of energy arrivals at time and the harvested energy at the -th arrival. Under the HSU architecture, to characterize the storage model in the battery, some authors use the so-called compound Poisson dam storage process [14, 15]. So, the energy stored in the battery at time is given bywhere represents the total energy expenditure until time being the power delivered. In the case of a finite capacity battery (), expression (2) incorporates a process () called reflection process, that it is null at zero, nondecreasing, continuous almost everywhere as follows:This reflection process ensures that the storage process does not exceed the threshold for any energy arrival, These stochastic models assume a discrete arrival scheme for the incoming energy (more discrete models of energy arrival can be found in [3, 12]). However, energy harvesting devices receive their energy from external nondiscrete fluctuating sources (such as radio waves) and a continuous model would be a better model when these sources are used. Despite this fact, there are very limited results for the energy harvesting systems with continuous energy arrival except from a few papers such as [16, 17].

Considering a continuous scheme for the energy arrival, in this paper we propose a stochastic process to model the energy harvested by an RF harvesting system. This stochastic process is the gamma process. The reason to choose this process is based on electromagnetic propagation statistical models and its use is justified in the paper. When the electromagnetic RF channel is dominated by nonline-of-sight (LOS) propagation, the energy stored in the battery can be modeled accurately by a gamma process.

Before defining a gamma process, the definition of gamma distribution is first recalled. We say that a random variable follows a gamma distribution with parameters if its probability density function is given bywhere is the gamma function given byIn mathematical terms, an homogeneous gamma process with shape parameter and scale parameter is a continuous-time stochastic process with the following properties: (i).(ii)Increments are independent for all and .(iii)The increments follow a gamma distribution with parameters and .

A gamma process is a pure-jump increasing Lévy process mainly used to model the deterioration of a system over the time [18]. Other examples of the application of gamma processes are the theory of water storage by dams and the risk theory [19]. The gamma process can be thought as a compound Poisson process in which the Poisson rate tends to infinity while the sizes of the increments tend to zero in proportion [20]. Although the compound Poisson process has been used to model the harvested energy process, as far as the authors are concerned, the gamma process has not been used to model the harvested energy. An advantage of using this gamma process in the stochastic modeling is that the required mathematical calculations are relatively straightforward. We again emphasize that the harvested energy process is continuous rather than slotted [21].

The general framework of this paper is the following. The use of the gamma process to model the harvested energy is justified. Afterwards, this process models the harvested energy by an EHS. This EHS contains an antenna that gathers energy from the environment, a rectifier, and a battery with finite capacity to store the harvested energy. Furthermore, a device is powered by the battery. Due to the finite capacity of the battery, energy may overflow without being utilized and hence is wasted. To model the energy harvested by the system, a gamma process is used. The battery delivers to the device a deterministic quantity of energy at certain instants of time that we call feeding times. Denoting by the successive feeding times, we assume a uniform scheduling policy for the feeding times, that is, for each .

Some constraints are imposed to this system.(1)The energy stored in the battery at time , , is limited by the battery capacity (that is, ).(2)The energy drawn from the battery can not exceed the energy stored in the battery (*Energy Casualty*).(3)When the system detects that its stored energy is going to drop below a certain threshold, the feeding process to the device stops. This threshold corresponds to the minimum energy reserve of the battery that can not be consumed. Denoting by this threshold, the stored energy in the battery at time fulfills for all t.

Under this general framework, a transmission policy is implemented. Renewal process theory is used to analyze the optimal transmission policy.

The objectives of this paper and its main contributions are the following.(i)Firstly, the use of the gamma process as a probabilistic tool to model the energy harvested by a RF harvesting system is justified. Since the gamma process is a continuous stochastic process, the continuous scheme for the energy arrival in EHS system is extended. Up to our knowledge, there are not works in the literature that deal with the use of the gamma process to model the harvested energy. It is, therefore, the main novelty and contribution of this paper.(ii)Secondly, to establish an optimal transmission policy to optimize the functioning of the RF harvesting system. For that, probabilistic efficiency measures are obtained. Another novelty of the paper is the approach used to deal with the transmission policy. In this work, renewal techniques are used to find the optimal transmission policy, in particular, the renewal-reward theorem. Although the renewal techniques are not new in the harvesting system literature [22, 23], the use of the renewal reward theorem to achieve a transmission policy is not a common practice in energy harvesting system literature (except some particular cases as [24]).

To obtain these objectives, the paper is structured as follows. In Section 2, the stochastic process used to model the harvested energy is described and its use is justified. Section 3 describes the main assumptions of the system functioning. Sections 4 and 5 obtain probabilistic efficiency measures related to the system functioning to analyze the optimal transmission policy for this system. Section 6 shows some numerical examples associated with the system and finally Section 7 concludes.

#### 2. Harvested Energy Model as a Gamma Process

In this section, we shall analyze the gamma process as a feasible tool to model the energy harvested by an EHS located in a RF fading channel environment. To obtain a first practical confirmation of this hypothesis, the storage of the energy received by an antenna in a Rayleigh fading environment has been simulated using the Matlab* Communication System Toolbox*™. The channel is modeled with the command channel = rayleighchan (Ts, FD) with Ts the sampling period equals to and FD the maximum Doppler shift of 15 Hz that corresponds to a frequency-flat Rayleigh fading channel. We establish a mean received power at antenna output of . The instantaneous received power is harvested by a RF rectifier with no efficiency loss obtaining the energy harvested energy at time as where denotes the instantaneous received power at time .

In Figure 1, the simulated received power and the harvested energy by the system in the first are represented. This harvested energy can be easily interpreted as the realization of an increasing stochastic process. In particular, we shall show that the harvested energy can be modeled by a gamma process. To validate this assumption, we have performed goodness-of-fit tests to compare the distribution of the sampling increments for different values of with a gamma distribution of parameters and being and . These statistical goodness-of-fit tests characterize the discrepancy between the simulated harvested energy and the expected values under a theoretical gamma distribution [25]. The Kolmogorov-Smirnov test is used to validate the gamma distribution of the increments. The p-values of the test are represented in Figure 2. As we can see, Figure 2 shows how well the gamma distribution fits the observed data for (approximately 4 times the value of or 2 times the value of the inverse of the maximum Doppler shift). In general, if the increments of time are large compared to the length of the fading, we can consider that the harvested and stored energy can be modeled by a gamma process. In practice , the value of the increment of from which the simulated process fits as a gamma process is proportional to the correlation length of the fading. In any case, the length of the fading is tiny enough (usually of order of milliseconds) compared to the order of magnitude of the times in which the EHS supplies energy to other devices. Consequently, considering fine and coarse time increments, harvested and stored energy is well modeled by a gamma process except for very fine time requirements.

Nakagami-m fading model is a very general model which fits a lot of different empirical measurements for signal reception levels in different multipaths and multisignal reception. This distribution includes and generalizes Rayleigh fading (as the distribution used in the previous example) used for nonline-of-sight propagation models (NLOS). It also approximates well other set of channel variations dominated by line-of-sight (LOS) propagation as Rician fading. The density function of the Nakagami-m model with parameters and isfor , where is the gamma function given by (6). The Nakagami-m distribution has two parameters, a shape parameter and a second parameter controlling the spread. It is well known that if follows a Nakagami-m distribution with parameters and , follows a gamma distribution with parameters [26].

Since the received power is proportional to the square of the signal amplitude, assuming a signal amplitude modeled as a Nakagami-m distribution, the received power at time is gamma distributed with parameters with being the mean received power.

Assuming that there is not efficiency loss, the harvested energy by the EHS at time is given bywhere corresponds to the received power at time . In the sequel of this work, we denote by the harvested energy process.

##### 2.1. as a Gamma Process

Let be a partition of in as follows: where denotes the floor function. Considering this partition we approximate asSince follows a gamma distribution with parameters and , then the approximation of given by (11) also follows a gamma distribution with parameters and with since the sum of independent random and gamma distributed variables with the same scale parameter is also gamma distributed (*gamma additivity property*). Notice that Obviously, in this approximation, the parameter of the gamma process depends on the length of discretization of (11), which looks inconsistent. Nevertheless, neither small nor large can be used. Firstly, large values of might invalidate the Simpson rule in (11). Also, small values of might imply that and are correlated and the assumption of independence in the gamma additivity property is violated. For small intervals of time, the fading of the channel is highly correlated; otherwise, the independence assumption can be assumed as true. The exact computation of is out of the scope of this paper. In general, is inversely proportional to the shape factor of the Nakagami-m fading channel and directly proportional to the correlation length of the fading (because the length of the fading is strongly related to the best ). In the example shown in Figures 1 and 2, we obtained the values of , that is, equal to the power reception provided in the simulation and Hz which is approximately the double of the maximum Doppler frequency shift. Since the fading time variations are very fast compared to the orders of magnitude of time used in the functioning of the harvesting system, we can consider that the discrete approximation in (11) can be used in continuous mode without loss of accuracy. In fact, for the example given in Figure 2, we showed that, for intervals length bigger than , the hypothesis of gamma process fits well.

For the same received power at antenna output, in a Rayleigh Channel environment, simulations have been performed to obtain an empiric relation between the parameters and of the gamma process (used to model the harvested energy) and the fading length. In the case of the Rayleigh Channel simulation provided by Matlab, the Doppler frequency shift, FD, is the parameter which is inversely proportional to the fading length. Figure 3 shows the and parameters estimated for the gamma process versus the Doppler shift varying from 2 to 100 Hz. A total of 5000 simulations of of time length have been performed for each Doppler shift frequency. In Figure 3, we can see that is actually equal to the mean power received () independently of the Doppler shift and the parameter varies linearly with the Doppler frequency shift, approximately FD.

Therefore, for a Nakagami-m or a Rayleigh RF channel, the process is modeled using an homogeneous gamma process that fulfills the following properties:(i).(ii)Increments are independent for all and .(iii)The harvested energy in an interval of amplitude , , follows a gamma distribution with parameters and with density for where we recall that corresponds to the mean received power.

The mean level of harvested energy and the variance in an time interval of amplitude is given byIn the remainder of this work, the expectation operator is used to avoid confusion between energy and expectation. So, we have a simplified model for the harvesting process using a continuous cumulative stochastic process.

#### 3. System Description

We assume a harvesting system which contains an antenna, a RF rectifier and a battery to store the harvested energy. The battery powers a device. The system functioning is described as follows:(1)We suppose that no efficiency loss is present in the RF rectifier and the battery.(2)The fading of the channel is modeled using a Nakagami-m distribution with density given in (8).(3)Let be the harvested energy at the output of the rectifier at time , and follows an homogeneous gamma process with parameters and , where the ratio corresponds to the mean received power.(4)Let be the stored energy in the battery at time with , where denotes the initial stored energy. For operational reasons, the stored energy in the battery must always exceed a threshold .(5)We assume that the capacity of the battery is finite and equals to , hence The battery wear is assumed to be much slower than the time scale of the problem.(6)The device can be in two modes: ON and OFF depending on its* feeding process*. In ON mode, the device is operational and powered by the battery. In OFF mode the device is turned off and the feeding process of the device stops.(7)* Feeding times*. At certain instantaneous times named* feeding times* the device consumes units of energy performing, for example, an instantaneous measurement. Let be the feeding times. We consider a block-based policy for the feeding times where they are scheduled periodically during an ON state with for . We assume a large enough value of to consider that the gamma process model for the harvested energy is satisfied.(8)The standby consumption of the device is negligible compared to the periodic consumption.(9)* Device turn on.* At the time of the turn on, the device needs to perform transitory tasks that consume energy units. The turn on instant is also considered a feeding time. At short, there is an instantaneous requirement of energy of units in the turn on of the device. We assume that i.e., the device turn on is ensured.(10)Let be the stored energy in the battery right before the feeding time for (i)If , then the device is properly powered and it is called a* feasible feeding task*.(ii)If , the feeding of the device starts but, when the system detects that the stored energy is going to drop below , the feeding process is stopped. In this case, the device is not properly powered and it is called an* unfeasible feeding task*. Furthermore, there is a waste of energy units. Right after an* unfeasible feeding task*, the energy stored in the battery is equal to the threshold and the device enters in OFF mode until the energy stored in the battery exceeds again the level (see Figure 4). Denoting by the energy stored in the battery right after , we get that Notice that condition implies that as long as there is energy in the battery, the device remains in ON mode.

Figure 4 shows a realization of the process that represents the energy stored in the battery. Initially, the energy stored in the battery is equal to units and the device consumes energy units. Furthermore, the device consumes units of energy each units of time whenever the device is in ON mode. In red color, an* unfeasible feeding task* is shown since the energy stored in the battery was not enough to feed the device fulfilling the energy constraint given in (16). Right after an* unfeasible feeding task*, the OFF period starts and it remains in this mode until the energy stored in the battery exceeds the threshold .

At the activation of the device, its consumption is equal to . Assuming an activation at , we get thatTo fulfill the constraint given in (16), has to verify This lower limit is imposed by (17). Whenever (16) is fulfilled in , the total energy consumed by the device due to the feeding tasks at time , , is equal to and the total energy stored in the battery at time can be expressed as with expectationLet be replicas of the process for . Starting from , the system evolves as follows. For , we get that The harvested energy in is with density function given by (14). Right before , the energy stored in the battery is equal towhere is given by (19) and denotes the battery capacity with . The overflow in is produced when with probability given by where is given by (14). At time , if , the feeding task is performed properly and the device remains in ON mode. Otherwise, the feeding task is unfeasible and the device enters in OFF mode. Hence, the energy stored in the battery right after the feeding time is given bywhere is given by (25).

Let be the following event:that corresponds to a feasible feeding task at time . Given , next feeding time is scheduled at time and the energy harvested in is given by that follows a gamma distribution with density shown in (14). Hence, right before , the energy stored in the battery is given bywhere denotes the indicator function which equals 1 if the argument is true and 0 otherwise and is given by (29). At time , condition corresponds to a feasible feeding task. Otherwise, an unfeasible feeding task takes place consuming energy units. Hence, the energy stored in the battery right after is given by where is given by (29).

In a general setting, we assume that the feeding tasks performed at times have been feasible. Then, right before feeding time , the energy stored in the battery is given byand just after the feeding time, we get that the energy stored in the battery is equal to As before, in (33), the variable is gamma distributed with density given by (14).

#### 4. Time to the First Unfeasible Feeding Task

After the turn on of the device and before the first unfeasible feeding task, the device is in ON mode. After an unfeasible feeding task, the device enters in OFF mode and it remains in this mode until the stored energy reaches again the level . Let be the time from the turn on to the first unfeasible feeding task. Let be the index of the first unfeasible feeding task, that is,then .

Before computing the reliability function of , we define the mixture random variable given bywith and where the density of is shown in (14). The distribution of is given byNotice that corresponds to the energy stored in the battery between two feeding times and if the energy in the battery just after is equal to ().

Let be the reliability function of , that is, Next, the analytical formulation of is obtained.(i)For , by the condition given in (17).(ii)For , we get that where is given by (14). For , represents the probability of a feasible feeding time at time .(iii)For , we get that where and denote the density of probability of and given in (37).(iv)In a general setting, for with , we get that for where for , denotes the density of the variable given in (37).

#### 5. Renewal Process

As in previous sections, the device is switch on at time and the energy stored in the battery is equal to . The device remains in ON state until the first unfeasible feeding time (see Figure 4). After this unfeasible feeding time, the device enters in OFF mode and the system keeps harvesting energy until the energy stored in the battery exceeds the level , time in which the device again enters in ON mode. Let be the time in which the device is OFF. Next, we analyse the distribution of . By Assumption 10, at the time of the first unfeasible feeding time, that is, , where is given by (36), the energy stored in the battery is equal to After this unfeasible feeding task, the device remains in OFF mode until the energy stored in the battery reaches the level . Then, follows the same distribution as , where denotes the first hitting time to the level of an homogeneous gamma process. The distribution function of is given bywhere is given by (6) and denotes the incomplete gamma function given by (see [27, 28] for more details). Although the exact distribution of is given in (45), it is very difficult to compute in practice. However, some approximations to this distribution are given. A Birnbaum−Saunders (BS) distribution is used to approximate the distribution of the first hitting time in a gamma process [29]. So, we can approximate by a BirnbaumSaunders (BS) distribution with parameters Furthermore, if then the distribution of is approximated closely by the inverse Gaussian distribution with parameters Due to the jump discontinuities of the sample paths of a gamma process, the energy stored in the battery in the OFF period fulfills ; hence the overshoot is positive almost surely. To avoid complicated mathematical technicalities, it seems realistic from a practical point of view to disregard the overshoot and consider that the energy stored in the battery right after the end of the OFF period of the device is equal to , that is, As a matter of fact, in most studies, the overshoot of the gamma process is not mentioned at all [30].

A goal of this paper is the optimization of the functioning of this system. The renewal process and its generalization (the quasirenewal process) are useful tools to analyze the optimal transmission policy for a system.

To optimize the transmission policy, a renewal argument is used. We recall that a renewal process is an arrival process in which the interarrival intervals are independent and identically distributed random variables. In this paper, a* renewal cycle* takes place every time that the OFF mode of the device ends. So, let be the time between renewals, where are independent and identically distributed variables with and where is given by with expectation

Let be the number of feasible feeding tasks in a renewal cycle; then where is given by (36). Under the assumptions of the model since the initial consumption of the device is guaranteed. Following the previous results, we get that where is given in (14). In a general setting, for , we get that where is given by (43). It is easy to check thatHence, the expected time to a renewal cycle can be expressed aswhere this approximation has been obtained using the inverse Gaussian distribution if and where is given by (43). The expectation also can also be approximated byusing the Birnbaum−Saunders distribution.

Next, we describe an optimal transmission policy for this EHS. By optimal, we mean the policy that maximizes an objective function. Renewal techniques are used to approach the search of the optimal transmission policy. Let be the successive renewal cycles (as we explained before, a renewal cycle in this paper takes place every time the OFF mode of the device ends). Let be the total of feasible tasks performed in a renewal cycle. Denoting by the total of feasible tasks up to time and using the renewal-reward theorem, we get that Hence, the limit of the total of feasible tasks per unit time can be reduced to the behavior of this ratio in the first renewal cycle. To deal with the optimal transmission policy, the ratiowhere and are given by (59) and (58), respectively, is chosen as objective cost function. The optimal transmission policy corresponds to the value of (the time between feeding times) that maximizes the expected number of feasible feeding tasks per unit time in a renewal cycle, that is,The analytic formulation of given in (63) is not straightforward since it involves multidimensional integrals as we can check in (58) and (59). Traditional quadrature methods to solve these multidimensional integrals are outperformed and the Monte Carlo method is chosen to analyze the optimization problem formulated in (63). Monte Carlo method is an adequate method when integrals are multidimensional [31].

Next section shows some numerical examples of the search of the optimal transmission policy using Monte Carlo method. For that, replicas of the first renewal cycle are simulated and the expected values given in (62) are numerically estimated using these simulations. Optimal values for the transmission policy are obtained by inspection considering these estimations.

#### 6. Numerical Examples

We assume a battery with a finite capacity of milijoules () and an energy threshold of . This battery powers a device and we assume that its consumption is equal to in each feasible feeding task. Due to the transitory effects, at the device turn on, it consumes . The initial stored energy in the battery is equal to (see Figure 4 for a better understanding). We consider that the harvested energy comes from an RF channel with Nakagami -fading and a mean power harvested equals to .

The model used to describe the harvested energy is an homogeneous gamma process with parameters and . Parameter depends on fading rate. In this numerical example, three different rates have been considered , , and . Notice that is related to the fading length, so the examples show a range of possibilities: from moderate fast fading to slow fading depending on . A total of 50000 complete renewal cycles for each point have been performed to obtain the results.

In Figure 5, the expected number of feasible feedings per minute, is shown. For high values of , almost all the feeding times are feasible since there is enough energy in the battery to feed the device. Instead, high values of cause a lot of unfeasible feeding tasks and, therefore, the repetitive switch off of the device. The maximum value of feasible feeding tasks per unit time corresponds to the optimal policy for the system operation and it is reached at (7 feeding times per minute), , and according to the value of . It is relevant that, for high values of , the number of feasible feeding tasks decreases. So, for the maximum expected feasible feeding tasks per minute are equal to , for are , and for are . As one would hope, the expected number of feasible feeding tasks per unit time is closely related to the RF channel characteristics. In contrast, for low values of (that is, large number of feeding times per unit time), slow fading channel may be less disadvantageous, as in the present case.

Figure 6 shows the energy lost per unit time (power lost) versus . The lost power in the system is due to three factors: energy not stored in the battery due to its finite capacity (overflow), energy consumed in the device turn on and wasted energy due to the unfeasible feeding tasks. As we could expect, the minimum of the lost power is reached at the values described above () with values of , , and , respectively. Notice that the power received at the antenna is consequently, and the efficiency of the system is , , and for , , and , respectively.

Another interesting measure is the ratio of unfeasible feeding tasks, calculated as which is represented in Figure 7. In Figure 7, we can see that for small values of the ratio is reduced. As we expected, the lower (fast fading) the better ratio.

In general, fast fading channels (that correspond to small values of ) are more advantageous for harvesting systems since the storage battery fluctuations are much less susceptible to variations; i.e., the gamma process associated with the harvested energy presents a lower variance .

#### 7. Conclusions and Further Works

An EHS system that consists of an antenna, a RF rectifier and a battery used to store harvested energy and a device powered by the battery, is considered in this paper.

The harvested energy of ambient RF energy in a NLOS fading environment is modeled using a gamma process. As far as the authors are concerned, this stochastic process is new to the area of energy harvesting. Due to the properties of the gamma process, the mathematical model of the whole system is tractable and some probabilistic quantities are obtained. These quantities are useful for the construction and analysis of a transmission policy of this EHS. Since the functioning of the system can be described using renewal arguments, the optimal scheduling policy is numerically analyzed.

As we showed in the examples, the fading of the signal causes losses in the efficiency of the harvesting system but this loss can be reduced to reasonable limits in fast fading environments. Fast fading channels are more advantageous for harvesting systems since the storage battery are much less susceptible to variations rather than slow fading channels.

In this paper, we have considered a perfect and ideal energy storage. It means that the battery saves in a proper way the energy that it receives. However, energy storage inefficiency can be developed in further works. For example, an inefficiency which occurs at the time of the energy storage can be analyzed. Also, the energy demand of the device has been modeled as a ideal instantaneous process. Further works could consider a continuous demand of energy.

We have also assumed in this paper an uniform transmission policy that schedules the feeding times uniformly. However, a scheduling policy could be considered in future works where the future feeding times are planned taking into account the stored energy in previous feeding times.

#### Data Availability

No data were used to support this study.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.

#### Acknowledgments

This work is supported by the Gobierno de Extremadura, Spain (Project IB18073), and European Union (European Regional Development Funds, ERDF). For the first author, this research was also supported by Ministerio de Economía y Competitividad, Spain (Project MTM2015-63978-P). For the second and third author, this research was supported by Gobierno de Extremadura, ERDF and SferaOne (Project GlobalEnergy AA-16-0124-2), and by Ministerio de Economía y Competitividad, Spain (Project TEC2017-85376-C2-X-R). For the third author this research was supported also by Spanish Ministerio de Educación, Cultura y Deporte (FPU 00022/15).