This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Age-Oriented Opportunistic Relaying in Cooperative Status Update Systems with Stochastic Arrivals
thanks: This work is done when the first author Bohai Li is a visiting student at the Chinese University of Hong Kong.

Bohai Li1,2, He Chen2, Yong Zhou3, and Yonghui Li1 1School of Electrical and Information Engineering, The University of Sydney, Sydney, Australia 2Department of Information Engineering, The Chinese University of Hong Kong, Hong Kong SAR, China 3School of Information Science and Technology, ShanghaiTech University, Shanghai, China
1{bohai.li, yonghui.li}@sydney.edu.au, 2[email protected], 3[email protected]
Abstract

This paper considers a cooperative status update system with a source aiming to send randomly generated status updates to a designated destination as timely as possible with the help of a relay. We adopt a recently proposed concept, Age of Information (AoI), to characterize the timeliness of the status updates. We propose an age-oriented opportunistic relaying (AoR) protocol to reduce the AoI of the considered system. Specifically, the relay opportunistically replaces the source to retransmit the successfully received status updates that have not been correctly delivered to the destination, but the retransmission of the relay can be preempted by the arrival of a new status update at the source. By carefully analyzing the evolution of AoI, we derive a closed-form expression of the average AoI for the proposed AoR protocol. We further minimize the average AoI by optimizing the generation probability of the status updates at the source. Simulation results validate our theoretical analysis and demonstrate that the average AoI performance of the proposed AoR protocol is superior to that of the non-cooperative system.

I Introduction

With the rapid development of the Internet of Things (IoT) applications, timely status updates have become increasingly critical in many areas, such as telesurgery and self-driving cars [1]. However, the conventional performance metrics, e.g., throughput and delay, cannot properly characterize the timeliness of the status updates [2]. For example, throughput can be maximized by generating and transmitting the status updates as frequent as possible. However, excessive update rates may lead to network congestion, which makes the updates suffer from long transmission delay. Such long delay can be reduced by lowing the update rate. However, reducing the rate of updates overmuch can also make the monitor receive undesired outdated status updates. Motivated by this, age of information (AoI), defined as the time elapsed since the generation of the latest received status update, has been recently introduced to quantify the information freshness from the perspective of the receiver that monitors a remote process [3]. Unlike conventional metrics, AoI is related to both the transmission delay and the update generation rate [4]. As a result, AoI is a more comprehensive evaluation criterion for information freshness.

Since the AoI concept was first proposed to characterize the information freshness in a vehicular status update system [5], extensive work focusing on the analysis and optimization of AoI has appeared. Most of the existing work is concerned with the AoI of single-hop wireless networks [6, 7, 8, 9, 10, 11, 12, 13]. On the other hand, the AoI performance in multi-hop networks has also been studied in [14, 15, 16, 17]. However, all the above work overlooked the direct link between the source and the destination. Therefore, the updates from the source can only be transmitted to the destination via the relay. As far as we know, there is no existing work that designs and optimizes the average AoI of a cooperative status update system with the existence of a direct link. Such a design is indeed non-trivial. This is because delivering the status updates via the direct link will have shorter transmission time at the cost of a higher error probability, while the delivery of status updates through the two-hop relay link could be more reliable at the cost of longer transmission time.

Motivated by this gap, in this paper we consider a three-node cooperative status update system, where a source timely reports randomly generated status updates to its destination with the help of a relay. We concentrate on designing and analyzing an age-oriented opportunistic relaying (AoR) protocol to reduce the AoI of the considered system. By carefully observing the AoI evolution process, we first define some necessary time intervals to mathematically express the average AoI of the proposed protocol. By representing these time-interval definitions in terms of key system parameters, including the generation probability of the status updates and the transmission success probabilities of three links, we then attain a closed-form expression of the average AoI for the proposed AoR protocol. Given reasonable values of the transmission success probabilities, we further minimize the average AoI by optimizing the status generation probability at the source. Simulation results are then provided to validate the theoretical analysis. Furthermore, the simulation results demonstrate that the average AoI of the proposed AoR protocol is smaller than that of the non-cooperative system, especially when the source-destination link has bad channel conditions.

II System Model and Protocol Description

We consider a three-node cooperative status update system, in which the source (S) wants to send randomly generated status updates to a designated destination (D) as timely as possible with the help of a relay (R). There is a direct link between S and D. The transmission of the status updates can either go through the direct link or the S-R-D link. To quantify the timeliness of the status updates, we adopt a recently proposed metric, termed AoI, which is defined as the time elapsed since the generation of the last successfully received status update. It is assumed that the considered system is time slotted and the transmission of each status update by S/R takes exactly one time slot, whose length is normalized to one without loss of generality. All the channels suffer from block fading, i.e., the channels responses remain unchanged within one time slot but vary independently from one time slot to another.

II-A Protocol Description

We now propose an AoR protocol for the considered system. S transmits a newly generated status update to D. If the status update is successfully received by D, then an acknowledgement (ACK) packet is sent back to S to discard the status update at the end of the current time slot. Otherwise, either S or R retransmits the status update to reduce the average AoI. Specifically, if R successfully receives the status update, then R feeds back an ACK to inform S to discard the status update, and transmits the status update to D by itself in the following time slot(s). Note that the channel gain of the R-D link is generally better than that of the S-D link. Retransmitting the same status update packets by R instead of S can reduce the AoI. If R fails to receive the status update, then S keeps retransmitting until either D or R succeeds in receiving the said status update. To ensure that fresh status updates can be transmitted timely, we enforce S to preempt the retransmissions from R in case of a new status update arrival. Hence, if R has a status update to retransmit, it first senses the state of S at the beginning of one time slot. If S with a new status update transmits, then R discards its current status packet and attempts to decode the new one sent by S. Otherwise, R keeps retransmitting the update until the successful reception at D or the preemption by S. For simplicity, we assume that both the time for sensing and transmitting an ACK packet are negligible.

It is worth mentioning that in the conventional cooperation protocols that are designed mainly from a physical layer perspective, normally allocate dedicated channel resources for R to facilitate the cooperation [18]. However, if the transmission success probability of the S-D link is large, it is apparent that transmitting more status updates through the S-D link reduces the AoI. Therefore, conventional cooperation protocols become sub-optimal from the perspective of minimizing the AoI. We note that a novel cooperation protocol, where R utilizes the silence periods of S terminals to enable cooperation, was proposed in [19] to improve the cooperative system performance. However, the design in [19] focused on the throughout maximization rather than AoI minimization. As far as we know, the AoR protocol is the first effort towards the design and analysis for the considered system from the perspective of AoI.

Refer to caption

Figure 1: Sample staircase path of the AoI evolution. We use SDS\!\rightarrow\!D, SRS\!\rightarrow\!R, RDR\!\rightarrow\!D and Failure to denote the successful transmission through the S-D link, the successful transmission through the S-R link, the successful transmission through the R-D link and the failed transmission, respectively.

II-B AoI Definition

We follow [6], [8] and adopt a Bernoulli process to model the stochastic arrivals at S. Particularly, a new status update is generated with probability pp at the beginning of each time slot. We denote tit_{i} as the generation time of the iith update at S. Taking into account channel fading, we define by P1P_{1}, P2P_{2}, and P3P_{3} to denote the transmission success probabilities through the S-D link, S-R link and R-D link, respectively. Due to the transmission error and the aforementioned preemption, the status updates generated at S may not always be successfully received by D. Hence, we further denote tkt^{\prime}_{k} as the arrival time of the kkth status update successfully received by D. We denote U(t)U(t) as the generation time of the most recently received status update at D until time slot tt. The AoI at time slot tt can then be defined as

Δ(t)=tU(t).\Delta(t)=t-U(t).\vspace{-0.1cm} (1)

For the ease of understanding the AoI evolution, we illustrate an example staircase path for 10 consecutive time slots with an initial value of one in Fig. 1. Note that we also define some time intervals in Fig. 1 to facilitate the calculation of the AoI, which will be explained one by one in the following. We denote by NkN_{k} the index of the most recently generated update at S before tkt^{\prime}_{k}, i.e.,

Nk=max{iti<tk},N_{k}=\max\{i\!\mid\!t_{i}<t^{\prime}_{k}\}, (2)

and denote by NkN^{\prime}_{k} the index of the first generated update at S since the last successfully received update at D, i.e.,

Nk=min{iti>tk1}.N^{\prime}_{k}=\min\{i\!\mid\!t_{i}>t^{\prime}_{k-1}\}. (3)

We further define SkS_{k} as the service time of the kkth successfully received update at D, given by

Sk=tktNk.S_{k}=t^{\prime}_{k}-t_{N_{k}}. (4)

In addition, WkW_{k} is defined as the waiting time starting from the arrival of the (k1)(k-1)th successfully received update at D until the generation of the first status update at S after tk1t^{\prime}_{k-1}, which is given by

Wk=tNktk1.W_{k}=t_{N^{\prime}_{k}}-t^{\prime}_{k-1}. (5)

We define TkT_{k} as the time duration between the generation time of the first status update after tk1t^{\prime}_{k-1} and the arrival time of the kkth successfully received update at D, which is given by

Tk=tktNk,T_{k}=t^{\prime}_{k}-t_{N^{\prime}_{k}}, (6)

and ZkZ_{k} is defined as interdeparture time between two consecutive successfully received status updates at D, given by

Zk=tktk1.Z_{k}=t^{\prime}_{k}-t^{\prime}_{k-1}. (7)

In the following section, we shall analyze the average AoI of the proposed AoR protocol for the considered cooperative status update system.

III Analysis and Optimization of Average AoI

In this section, we first express the average AoI by the defined time intervals given in Section II. By representing these time-interval definitions in terms of key system parameters, (i.e., pp, P1P_{1}, P2P_{2} and P3P_{3}), we then manage to attain a closed-form expression of the average AoI for the proposed AoR protocol. Given reasonable value sets of P1P_{1}, P2P_{2} and P3P_{3}, we further minimize the average AoI by optimizing the status generation probability pp.

III-A Analysis of Average AoI

We define XtX_{t} as the number of successfully received updates at D till time tt and it can be expressed as

Xt=max{ktk<t}.X_{t}=\max\{k\!\mid\!t^{\prime}_{k}<t\}.\vspace{-0.075cm} (8)

As depicted in Fig. 1, the average AoI can be expressed by the area under the AoI curve, i.e., QkQ_{k}, and mathematically we have

Δ¯=limtXtt1Xtk=1XtQk=𝔼[Qk]𝔼[Zk].\bar{\Delta}=\lim\limits_{t\to\infty}{\frac{X_{t}}{t}\frac{1}{X_{t}}\sum_{k=1}^{X_{t}}Q_{k}}=\frac{\mathbb{E}[Q_{k}]}{\mathbb{E}[Z_{k}]}. (9)

For further simplification, we represent QkQ_{k} in terms of the time intervals defined in Section II and it can be expressed as

Qk\displaystyle Q_{k} =Sk1+(Sk1+1)++(Sk1+Zk1)\displaystyle=S_{k-1}+(S_{k-1}+1)+\cdots+(S_{k-1}+Z_{k}-1) (10)
=Sk1Zk+Zk2Zk2.\displaystyle=S_{k-1}Z_{k}+\frac{Z_{k}^{2}-Z_{k}}{2}.

By substituting (10) into (9), we can express the average AoI as

Δ¯=𝔼[Sk1Zk]𝔼[Zk]+𝔼[Zk2]2𝔼[Zk]12.\bar{\Delta}=\frac{\mathbb{E}[S_{k-1}Z_{k}]}{\mathbb{E}[Z_{k}]}+\frac{\mathbb{E}\Big{[}Z_{k}^{2}\Big{]}}{2\mathbb{E}[Z_{k}]}-\frac{1}{2}.\vspace{-0.05cm} (11)

In order to further simplify (11), we have the following lemma.

Lemma 1: The service time of the (k1)(k-1)th successfully received update at D is independent of the interdeparture time between the (k1)(k-1)th and the kkth successfully received updates at D such that:

𝔼[Sk1Zk]=𝔼[Sk1]𝔼[Zk].\mathbb{E}[S_{k-1}Z_{k}]=\mathbb{E}[S_{k-1}]\mathbb{E}[Z_{k}]. (12)

The lemma can be proved by following a similar process as that for [6, Lemma 1]

By applying Lemma 1, the expression of the average AoI can be simplified as

Δ¯=𝔼[Sk1]+𝔼[Zk2]2𝔼[Zk]12.\bar{\Delta}=\mathbb{E}[S_{k-1}]+\frac{\mathbb{E}\Big{[}Z_{k}^{2}\Big{]}}{2\mathbb{E}[Z_{k}]}-\frac{1}{2}. (13)

To obtain the average AoI, we now derive the terms 𝔼[Sk1]\mathbb{E}[S_{k-1}], 𝔼[Zk]\mathbb{E}[Z_{k}] and 𝔼[Zk2]\mathbb{E}[Z_{k}^{2}] one by one in the following. We start with the expectation of the service time, i.e., 𝔼[Sk1]\mathbb{E}[S_{k-1}]. In the considered system, there are two types of status updates that can be successfully received by D without being preempted. One are these packets that are successfully received by D through the direct link, and the other are these packets that are successfully delivered through the two-hop S-R-D link. Therefore, 𝔼[Sk1]\mathbb{E}[S_{k-1}] can be expressed as a weighted sum of the average service time of each type. For the first type of successful updates, the probability that the update is successfully received after ll times of transmission by S can be expressed as

Pl=(1p)l1(1P2)l1(1P1)l1P1.P_{l}=(1-p)^{l-1}(1-P_{2})^{l-1}(1-P_{1})^{l-1}P_{1}.\vspace{-0.15cm} (14)

The corresponding expectation of the service time can be calculated by

𝔼l=l=1Plll=1Pl.\mathbb{E}_{l}=\frac{\sum_{l=1}^{\infty}P_{l}\cdot l}{\sum_{l=1}^{\infty}P_{l}}. (15)

Similarly, for the second type of successful status update, the expectation of its service time can be expressed by

𝔼mn=n=1m=1Pmn(m+n)n=1m=1Pmn,\mathbb{E}_{mn}=\frac{\sum_{n=1}^{\infty}\sum_{m=1}^{\infty}P_{mn}\cdot(m+n)}{\sum_{n=1}^{\infty}\sum_{m=1}^{\infty}P_{mn}}, (16)

where mm and nn denote the number of transmission times via the S-R link and the R-D link, respectively. PmnP_{mn} is the probability that the update is successfully received by D after transmitting (m+n)(m+n) times, which can be expressed as

Pmn=(1p)m1(1P2)m1(1P1)mP2(1p)n(1P3)n1P3.\leavevmode\resizebox{460.35788pt}{}{$P_{mn}\!=\!(1-p)^{m-1}(1-P_{2})^{m-1}(1-P_{1})^{m}P_{2}(1-p)^{n}(1-P_{3})^{n-1}P_{3}$}. (17)

Since the above two types of updates make up all the updates that can be successfully received by D, the expectation of the service time of the considered system can be evaluated as

𝔼[Sk1]=𝔼ll=1Pll=1Pl+n=1m=1Pmn+𝔼mnn=1m=1Pmnl=1Pl+n=1m=1Pmn,\begin{aligned} \mathbb{E}[S_{k-1}]=&\ \mathbb{E}_{l}\cdot\frac{\sum_{l=1}^{\infty}P_{l}}{\sum_{l=1}^{\infty}P_{l}+\sum_{n=1}^{\infty}\!\sum_{m=1}^{\infty}\!P_{mn}}\\ &\,+\mathbb{E}_{mn}\cdot\frac{\sum_{n=1}^{\infty}\!\sum_{m=1}^{\infty}\!P_{mn}}{\sum_{l=1}^{\infty}P_{l}+\sum_{n=1}^{\infty}\!\sum_{m=1}^{\infty}\!P_{mn}},\end{aligned}

(18)

which follows according to the law of total probability.

By applying the finite sum equations given in [20, Eqs. (0.112) and (0.113)], we can further simplify (18) to

𝔼[Sk1]=11β+11αγP1(1α)+γ,\mathbb{E}[S_{k-1}]=\frac{1}{1-\beta}+\frac{1}{1-\alpha}\cdot\frac{\gamma}{P_{1}(1-\alpha)+\gamma}, (19)

where α=(1p)(1P3)\alpha=(1-p)(1-P_{3}), β=(1p)(1P1)(1P2)\beta=(1-p)(1-P_{1})(1-P_{2}) and γ=P2P3(1p)(1P1)\gamma=P_{2}P_{3}(1-p)(1-P_{1}).

Based on (5), (6) and (7), we have Zk=Wk+TkZ_{k}=W_{k}+T_{k}. Therefore, the expectation of the interdeparture time can be written as 𝔼[Zk]=𝔼[Wk]+𝔼[Tk]\mathbb{E}[Z_{k}]=\mathbb{E}[W_{k}]+\mathbb{E}[T_{k}]. Recall that S can preempt the transmission of R when there is a new arrival. This indicates that the waiting time WkW_{k} only depends on the generation probability pp. Since the status updates are generated according to a Bernoulli process, WkW_{k} follows a geometric distribution with parameter pp and its expectation can be readily given by 𝔼[Wk]=(1p)/p\mathbb{E}[W_{k}]=(1-p)/p. We then move to the calculation of term 𝔼[Tk]\mathbb{E}[T_{k}].

Note that TkT_{k} behaves differently for the following four possible cases: 1) The update is successfully received by D through the S-D link without being preempted; 2) The update is successfully received by D through the S-R-D link without being preempted; 3) The update is preempted by a new update before being successfully received by either R or D, and the new update may be preempted by multiple new updates; 4) The update is preempted by a new update generated at S after being successfully received by R, and the new update may be preempted by multiple new updates at S. We note that the number of preemptions can approach infinity in the third and fourth cases, which generally makes the expectation of TkT_{k} difficult to derive mathematically. To tackle this issue, we resort to the recursive method applied in [6],[21] and evaluate the expectation of TkT_{k} as (20) at the top of the next page.

𝔼[Tk]=l=1Pll+n=1m=1Pmn(m+n)+l=1(1P1)l(1P2)l(1p)l1p(l+𝔼[Tk])+m=1n=0[(1P1)m(1P2)m1P2(1p)m1×(1p)n(1P3)np(m+n+𝔼[Tk])]\begin{aligned} \mathbb{E}[T_{k}]=&\ \sum_{l=1}^{\infty}P_{l}\cdot l+\sum_{n=1}^{\infty}\sum_{m=1}^{\infty}P_{mn}\cdot(m+n)+\sum_{l=1}^{\infty}(1-P_{1})^{l}(1-P_{2})^{l}(1-p)^{l-1}p\cdot\big{(}l+\mathbb{E}\big{[}T^{\prime}_{k}\big{]}\big{)}\\ &\,+\sum_{m=1}^{\infty}\sum_{n=0}^{\infty}\Big{[}(1-P_{1})^{m}(1-P_{2})^{m-1}P_{2}(1-p)^{m-1}\times(1-p)^{n}(1-P_{3})^{n}p\cdot\big{(}m+n+\mathbb{E}\big{[}T^{\prime}_{k}\big{]}\big{)}\!\Big{]}\\ \end{aligned}

(20)

 

Note that the four terms on the right hand side of (20) correspond to the above four cases, respectively, and TkT^{\prime}_{k} is the time duration between the preemption time (i.e., the generation time of the second status update after tk1t^{\prime}_{k-1}) and the arrival time of the kkth successfully received update at D. As TkT_{k} is defined as the time duration between the generation time of the first status update after tk1t^{\prime}_{k-1} and the arrival time of the kkth successfully received update at D, following the idea of recursion, we have 𝔼[Tk]=𝔼[Tk]\mathbb{E}[T_{k}]=\mathbb{E}\big{[}T^{\prime}_{k}\big{]}. After some manipulations by applying [20, Eqs. (0.112) and (0.113)], we have

𝔼[Tk]=(1α)+βP2P1(1α)+γ,\mathbb{E}[T_{k}]=\frac{(1-\alpha)+\beta\cdot P^{\prime}_{2}}{P_{1}(1-\alpha)+\gamma},\vspace{-0.1cm} (21)

where P2=P2/(1P2)P^{\prime}_{2}=P_{2}/(1-P_{2}). By combining 𝔼[Wk]=(1p)/p\mathbb{E}[W_{k}]=(1-p)/p and (21), we now attain a closed-form expression for the term 𝔼[Zk]\mathbb{E}[Z_{k}], given by

𝔼[Zk]=𝔼[Wk]+𝔼[Tk]=(1α)(1β)p[P1(1α)+γ].\displaystyle\mathbb{E}[Z_{k}]=\mathbb{E}[W_{k}]+\mathbb{E}[T_{k}]=\frac{(1-\alpha)(1-\beta)}{p\big{[}P_{1}(1-\alpha)+\gamma\big{]}}. (22)

Finally, we work on the expectation of Zk2Z_{k}^{2}, which can be expressed as

𝔼[Zk2]=𝔼[(Wk+Tk)2]=𝔼[Wk2]+2𝔼[WkTk]+𝔼[Tk2].\displaystyle\mathbb{E}\Big{[}Z_{k}^{2}\Big{]}=\mathbb{E}\Big{[}(W_{k}+T_{k})^{2}\Big{]}=\mathbb{E}\Big{[}W_{k}^{2}\Big{]}+2\mathbb{E}[W_{k}T_{k}]+\mathbb{E}\Big{[}T_{k}^{2}\Big{]}. (23)

It is readily to find that WkW_{k} and TkT_{k} are independent. Thus, we have 𝔼[WkTk]=𝔼[Wk]𝔼[Tk]\mathbb{E}[W_{k}T_{k}]=\mathbb{E}[W_{k}]\mathbb{E}[T_{k}]. As WkW_{k} follows a geometric distribution with parameter pp, we have 𝔼[Wk2]=(p23p+2)/p2\mathbb{E}[W_{k}^{2}]=(p^{2}-3p+2)/p^{2}. In order to evaluate 𝔼[Zk2]\mathbb{E}[Z_{k}^{2}], the only remaining task is to calculate 𝔼[Tk2]\mathbb{E}[T_{k}^{2}].

Similar to (20), the expectation of Tk2T_{k}^{2} can be evaluated as (24) at the top of the next page.

𝔼[Tk2]=l=1Pll2+n=1m=1Pmn(m2+2mn+n2)+l=1[(1P1)l(1P2)l(1p)l1p(l2+2l𝔼[Tk]+𝔼[Tk2])]+m=1n=0[(1P1)m(1P2)m1P2(1p)m1(1p)n(1P3)np(m2+n2+𝔼[Tk2]+2mn+2m𝔼[Tk]+2n𝔼[Tk])]\begin{aligned} \mathbb{E}\Big{[}T_{k}^{2}\Big{]}\!=\!\!&\ \sum_{l=1}^{\infty}P_{l}\!\cdot\!l^{2}\!+\!\sum_{n=1}^{\infty}\!\sum_{m=1}^{\infty}\!P_{mn}\!\cdot\!(m^{2}\!+\!2mn\!+\!n^{2})\!+\!\sum_{l=1}^{\infty}\!\bigg{[}\!(1\!-\!P_{1})^{l}(1\!-\!P_{2})^{l}(1\!-\!p)^{l-1}p\!\cdot\!\Big{(}l^{2}\!+\!2l\mathbb{E}\big{[}T_{k}\big{]}\!+\!\mathbb{E}\Big{[}T_{k}^{2}\Big{]}\Big{)}\!\bigg{]}\\ &\,\!+\!\sum_{m=1}^{\infty}\!\sum_{n=0}^{\infty}\!\bigg{[}\!(1\!-\!P_{1})^{m}(1\!-\!P_{2})^{m-1}P_{2}(1\!-\!p)^{m-1}(1\!-\!p)^{n}(1\!-\!P_{3})^{n}p\!\cdot\!\Big{(}\!m^{2}\!+\!n^{2}\!+\!\mathbb{E}\Big{[}T_{k}^{2}\Big{]}\!+\!2mn\!+\!2m\mathbb{E}[T_{k}]\!+\!2n\mathbb{E}[T_{k}]\!\Big{)}\!\bigg{]}\end{aligned}

(24)

 

Similar to the process of obtaining (21) from (20), we can rewrite 𝔼[Tk2]\mathbb{E}[T_{k}^{2}] as

𝔼[Tk2]=1(1α)(1β)(1α)βpβpP2×1(1α)(1β)×{(1α)2(1+β)+(3αβαβ)βP2+2𝔼[TK][(1α)2βp+(1αβ)βpP2]},\begin{aligned} \mathbb{E}\Big{[}T_{k}^{2}\Big{]}\!=&\,\frac{1}{(1-\alpha)(1-\beta)-(1-\alpha)\beta\cdot p^{\prime}-\beta\cdot p^{\prime}\cdot P^{\prime}_{2}}\\ &\times\!\frac{1}{(1-\alpha)(1-\beta)}\\ &\times\!\bigg{\{}(1-\alpha)^{2}(1+\beta)+(3-\alpha-\beta-\alpha\beta)\beta\cdot P^{\prime}_{2}\\ &\qquad\!\!+2\mathbb{E}[T_{K}]\!\cdot\!\Big{[}(1\!-\!\alpha)^{2}\beta\!\cdot\!p^{\prime}\!+\!(1\!-\!\alpha\beta)\beta\!\cdot\!p^{\prime}\!\cdot\!P^{\prime}_{2}\Big{]}\bigg{\}},\end{aligned}

(25)

where p=p/(1p)p^{\prime}=p/(1-p).

By substituting the terms derived in (19), (22), (23) and (25) into (13), we can obtain the exact closed-form expression of the average AoI for the proposed AoR protocol, given by

Δ¯=[1(1p)(1P3)][1(1p)(1P1)(1P2)]p[pP1+(1p)P3(1p)(1P1)(1P2)P3].\bar{\Delta}\!=\!\frac{\big{[}1\!-\!(1\!-\!p)(1\!-\!P_{3})\big{]}\!\cdot\!\big{[}1\!-\!(1\!-\!p)(1\!-\!P_{1})(1\!-\!P_{2})\big{]}}{p\!\cdot\!\big{[}pP_{1}\!+\!(1\!-\!p)P_{3}\!-\!(1\!-\!p)(1\!-\!P_{1})(1\!-\!P_{2})P_{3}\big{]}}. (26)

III-B Optimization of Average AoI

In this subsection, given reasonable value sets of P1P_{1}, P2P_{2} and P3P_{3}, we will minimize the average AoI by optimizing the generation probability pp at S.

By assuming that 0<P1<P2<10<P_{1}<P_{2}<1 and 0<P1<P3<10<P_{1}<P_{3}<1, the optimal generation probability pp that minimizes the average AoI is presented in the following theorem.

Theorem 1: The optimal generation probability that minimizes the average AoI is given by:

p={ψ+ψ24χω2χ,if 0<P1<P2+P2P3(P2P2P3)2+4P2P32(P21),1,ifP2+P2P3(P2P2P3)2+4P2P32(P21)<P1<1,p=\left\{\begin{aligned} &\frac{-\!\psi\!+\!\sqrt{\psi^{2}\!-\!4\chi\omega}}{2\chi},\\ &\quad\mathrm{if}\ 0<P_{1}<\frac{P_{2}\!+\!P_{2}P_{3}\!-\!\sqrt{(P_{2}\!-\!P_{2}P_{3})^{2}\!+\!4P_{2}P_{3}}}{2(P_{2}\!-\!1)},\\ &1,\\ &\quad\mathrm{if}\ \frac{P_{2}\!+\!P_{2}P_{3}\!-\!\sqrt{(P_{2}\!-\!P_{2}P_{3})^{2}\!+\!4P_{2}P_{3}}}{2(P_{2}\!-\!1)}<P_{1}<1,\end{aligned}\right. (27)

where χ=P2P3(1P2P3)P12(1P2)(1P3)2P1P2P32(1P2)(2P1)P1P2(1P3)\chi\!=\!P_{2}P_{3}(1\!-\!P_{2}P_{3})\!-\!P_{1}^{2}(1\!-\!P_{2})(1\!-\!P_{3})^{2}\!-\!P_{1}P_{2}P_{3}^{2}(1\!-\!P_{2})(2\!-\!P_{1})\!-\!P_{1}P_{2}(1\!-\!P_{3}), ψ=2P3(P1+P2P1P2)(P1P1P3P2P3+P1P2P3)\psi\!=\!-2P_{3}(P_{1}\!+\!P_{2}\!-\!P_{1}P_{2})(P_{1}\!-\!P_{1}P_{3}\!-\!P_{2}P_{3}\!+\!P_{1}P_{2}P_{3}) and ψ24χω=4P2P32(P3P1)(1P1)(P1+P2P1P2)2\psi^{2}\!-\!4\chi\omega\!=\!4P_{2}P_{3}^{2}(P_{3}\!-\!P_{1})(1\!-\!P_{1})(P_{1}\!+\!P_{2}\!-\!P_{1}P_{2})^{2}.

Proof: See Appendix A.

Remark 1: For the second case in Theorem 1, the average AoI is minimized by setting p=1p=1 and it is given by

Δ¯min=1P1.\bar{\Delta}_{min}=\frac{1}{P_{1}}. (28)

Note that in the considered system, the minimum possible time required for a successful status update through the S-D link is one time slot, while that through the S-R-D link is at least two time slots. Therefore, when the transmission success probability of the S-D link (i.e., P1P_{1}) is relatively large, having more status updates successfully delivered through the S-D link is beneficial to reduce the average AoI. If the generation probability of the updates is 1, S always has a new status update to transmit (i.e., the generate-at-will model [22]), which means that the updates successfully received by R in the previous time slot is always preempted. As a result, all the successfully received status updates by D can only come from the S-D link, which will minimize the average AoI.

On the other hand, for the first case in Theorem 1, the transmission success probability of the S-D link is relatively small compared to the second case. Hence, the updates successfully delivered through the S-R-D link also have a significant contribution on minimizing the average AoI. As mentioned above, if the generation probability is 1, the updates can only be successfully delivered to D through the S-D link. Obviously, this is not optimal for minimizing the average AoI in this case, which means that p=1p=1 should not be the optimal status generation probability to minimize the average AoI.

IV Simulation Results and Discussions

In this section, we present simulation results to validate our analytical results, and demonstrate the superiority of the considered cooperative status update system over its non-cooperative counterpart in terms of achieving lower average AoI. For the non-cooperative status update system, we plot its average AoI curves according to the result given in [6, Eqs. (35)]. Each curve presented in the section is averaged over 10710^{7} time slots.

We first study the average AoI of both cooperative and non-cooperative systems against the the transmission success probability of the S-D link (i.e., P1P_{1}), as shown in Fig. 2, for two different generation probabilities of the status updates. In both cases, we assume that the transmission success probabilities of the S-R link and the R-D link are P2=P3=0.8P_{2}\!\!=\!\!P_{3}\!\!=\!\!0.8. It can be observed that in both cases, the performance of the considered cooperative status update system is obviously superior to that of the non-cooperative system, especially when P1P_{1} is small. Specifically, when p=0.5p\!=\!0.5 (or p=0.9p\!=\!0.9), the considered cooperative system has a significant performance improvements over the non-cooperative system for the case P1<0.5P_{1}\!<\!0.5 (or P1<0.3P_{1}\!<\!0.3), respectively. This is because when the S-D link has bad channel conditions, the generated updates at S can be transmitted to D via R in our system while those in the non-cooperative system can only be retransmitted via the direct link. However, we can also observe that the AoI performance gain decreases as P1P_{1} increases, and approaches 1 when P1P_{1} is close to 0.8 (i.e., equal to P2P_{2} and P3P_{3}). The rationality behind this is that when P1P_{1} is large, most of the successfully received updates at D are transmitted via the S-D link. Therefore, R is not much helpful in reducing the average AoI. We also note that when P1P_{1} is less than 0.27, a lower generation probability, i.e., p=0.5p\!=\!0.5, results in a better AoI performance. On the contrary, the AoI performance of a higher generation probability, i.e., p=0.9p\!=\!0.9, is always better when P1P_{1} is greater than 0.27. Motivated by this interesting observation, we explore the relationship between the generation probability and the average AoI in the following.

We now investigate the average AoI against the generation probability pp by fixing the transmission success probabilities of the S-R link and the R-D link as P2=P3=0.8P_{2}\!=\!P_{3}\!=\!0.8. As depicted in Fig. 3, we also consider two cases with the transmission success probability of the S-D link being 0.25 and 0.5, respectively. In the first case of P1=0.25P_{1}=0.25, it is apparent that the average AoI first decreases and then increases as the generation probability increases. This is because if the S-D link has bad channel conditions, the updates successfully delivered to D by R have a significant contribution on minimizing the average AoI. However, as S can preempt the transmission of R whenever there is a new arrival, if the updates are generated too frequently, the updates at R will be always preempted by S and cannot be transmitted to D to effectively reduce the AoI. Therefore, too high generation probabilities in turn increase the average AoI. We can also find that the average AoI when p=0.5p=0.5 is smaller than that when p=0.9p=0.9, which confirms our mentioned observation when P1P_{1} is less than 0.27 in Fig. 2. In the second case of P1=0.5P_{1}=0.5, the average AoI always decreases as the generation probability increases. The monotone decreasing property demonstrates our mentioned phenomenon when P1P_{1} is greater than 0.27 in Fig. 2. It is also worth mentioning that the optimal generation probabilities to minimize the average AoI in the case of P1=0.25P_{1}\!=\!0.25 and P1=0.5P_{1}\!=\!0.5 are p=0.6770p=0.6770 and p=1p=1, respectively, which is consistent with the results calculated by using the formula in Theorem 1.

Refer to caption
Figure 2: Average AoI of the considered cooperative system and the non-cooperative system versus the transmission success probability P1P_{1} for different generation probabilities.
Refer to caption
Figure 3: Average AoI of the considered cooperative system and the non-coopertive system versus the generation probability pp for different transmission success probabilities, i.e., (P1,P2,P3)(P_{1},P_{2},P_{3}).

V Conclusions

In this paper, we proposed an age-oriented opportunistic relaying (AoR) protocol to reduce the age of information (AoI) of a cooperative status update system, where the source timely reports randomly generated status updates to the destination with the help of a relay. By analyzing the evolution of AoI, we derived a closed-form expression of the average AoI of the proposed AoR protocol as a function of the generation probability of the status updates and the transmission success probability of each link. Given reasonable assumptions about the transmission success probabilities, we further figured out the optimal generation probability of the status updates that minimizes the average AoI. Simulation results validated our theoretical analysis and demonstrated the superiority of the proposed AoR protocol over the non-cooperative system in terms of minimizing the average AoI, especially when the source-destination link is not as good as the other two links.

Appendix

A. Proof of Theorem 1

Recall that Δ¯=[1(1p)(1P3)][1(1p)(1P1)(1P2)]p[pP1+(1p)P3(1p)(1P1)(1P2)P3]\bar{\Delta}=\frac{[1-(1-p)(1-P_{3})]\cdot[1-(1-p)(1-P_{1})(1-P_{2})]}{p[pP_{1}+(1-p)P_{3}-(1-p)(1-P_{1})(1-P_{2})P_{3}]}. To proceed, we derive the first-order derivative of Δ¯\bar{\Delta} with respect to (w.r.t) pp. After some algebra manipulations, we have

Δ¯p=χp2+ψp+ωp2[pP1+(1p)P3(1p)(1P1)(1P2P3)]2,\frac{\partial\bar{\Delta}}{\partial p}=\frac{\chi\cdot p^{2}+\psi\cdot p+\omega}{p^{2}\!\cdot\!\big{[}pP_{1}\!+\!(1\!-\!p)P_{3}\!-\!(1\!-\!p)(1\!-\!P_{1})(1\!-\!P_{2}P_{3})\big{]}^{2}}, (29)

where

χ=(P12P22P322P12P2P32+2P12P2P3P12P2+P12P322P12P3+P122P1P22P32+2P1P2P32P1P2P3+P1P2+P22P32P2P3),\begin{aligned} \chi\!=\!-&(P_{1}^{2}P_{2}^{2}P_{3}^{2}\!-\!2P_{1}^{2}P_{2}P_{3}^{2}\!+\!2P_{1}^{2}P_{2}P_{3}\!-\!P_{1}^{2}P_{2}\!+\!P_{1}^{2}P_{3}^{2}\!-\!2P_{1}^{2}P_{3}\!\\ &\,+\!P_{1}^{2}\!-\!2P_{1}P_{2}^{2}P_{3}^{2}\!+\!2P_{1}P_{2}P_{3}^{2}\!-\!P_{1}P_{2}P_{3}\!+\!P_{1}P_{2}\!+\!P_{2}^{2}P_{3}^{2}\!-\!P_{2}P_{3}),\end{aligned}

(30a)

ψ=2P12P22P324P12P2P32+2P12P2P3+2P12P322P12P34P1P22P32+4P1P2P322P1P2P3+2P22P32,\begin{aligned} \psi\!=\,&2P_{1}^{2}P_{2}^{2}P_{3}^{2}\!-\!4P_{1}^{2}P_{2}P_{3}^{2}\!+\!2P_{1}^{2}P_{2}P_{3}\!+\!2P_{1}^{2}P_{3}^{2}\!-\!2P_{1}^{2}P_{3}\!\\ &\qquad\qquad\qquad-\!4P_{1}P_{2}^{2}P_{3}^{2}\!+\!4P_{1}P_{2}P_{3}^{2}\!-\!2P_{1}P_{2}P_{3}\!+\!2P_{2}^{2}P_{3}^{2},\end{aligned}

(30b)

ω=(P12P22P322P12P2P32+P12P322P1P22P32+2P1P2P32+P22P32),\begin{aligned} \omega\!=\!-(P_{1}^{2}P_{2}^{2}P_{3}^{2}\!-\!2P_{1}^{2}P_{2}P_{3}^{2}\!+\!P_{1}^{2}P_{3}^{2}\!-\!2P_{1}P_{2}^{2}P_{3}^{2}\!+\!2P_{1}P_{2}P_{3}^{2}\!+\!P_{2}^{2}P_{3}^{2}),\end{aligned}

(30c)

are defined for the notation simplicity.

After a careful observation on the right hand side (RHS) of (29), we can deduce that the sign of Δ¯p\frac{\partial\bar{\Delta}}{\partial p} is only determined by the numerator

κ(p)=χp2+ψp+ω,\kappa(p)=\chi\!\cdot\!p^{2}+\psi\!\cdot\!p+\omega, (31)

since the denominator is always large than zero. To determine the monotonicity of the function Δ¯\bar{\Delta}, we need to investigate the properties of the quadratic function κ(p)\kappa(p) on the feasible set of pp (i.e., [0, 1]).

Firstly, after some algebra manipulations, we note that

ω=P32[P12(P21)2+P2(P2+2P1(1P2))],\omega=-P_{3}^{2}\Big{[}P_{1}^{2}(P_{2}-1)^{2}+P_{2}\big{(}P_{2}+2P_{1}(1-P_{2})\big{)}\Big{]}, (32a)
ψ24χω=4P32P2(P3P1)(1P1)(P1+P2P1P2)2.\psi^{2}-4\chi\omega=4P_{3}^{2}P_{2}(P_{3}-P_{1})(1-P_{1})(P_{1}+P_{2}-P_{1}P_{2})^{2}. (32b)

As 0<P1<P2<10<P_{1}<P_{2}<1 and 0<P1<P3<10<P_{1}<P_{3}<1, we have that ω<0\omega<0 and ψ24χω>0\psi^{2}-4\chi\omega>0. This means that the curve of κ(p)\kappa(p) always intersects the Y-axis at the negative half of the Y-axis and it always has two roots on the X-axis, which can be given by

x1=ψψ24χω2χ,x_{1}=\frac{-\psi-\sqrt{\psi^{2}-4\chi\omega}}{2\chi}, (33a)
x2=ψ+ψ24χω2χ,x_{2}=\frac{-\psi+\sqrt{\psi^{2}-4\chi\omega}}{2\chi}, (33b)

where x2>x1x_{2}>x_{1}.

To further characterize the shapes of the function κ(p)\kappa(p), we now investigate χ\chi and ψ\psi for the following four possible cases:

1) When χ>0\chi>0 and ψ>0\psi>0, we always have x1<0x_{1}<0 and x2>0x_{2}>0. Since pp is in the range of [0,1][0,1], we have two sub-cases: (a) x2<1x_{2}<1 and (b) x2>1x_{2}>1. We draw the possible shapes of κ(p)\kappa(p) versus pp for the sub-case (a) and sub-case (b) in Fig. 4 (a) and Fig. 4 (b), respectively.

From Fig. 4 (a), we can see that in the sub-case (a), κ(p)<0\kappa(p)<0 holds for p[0,x2]p\in[0,x_{2}], and κ(p)>0\kappa(p)>0 holds for p(x2,1]p\in(x_{2},1]. This means that the function Δ¯\bar{\Delta} is decreasing for p[0,x2]p\in[0,x_{2}] and is increasing for p(x2,1]p\in(x_{2},1]. As a result, the minimum value of Δ¯\bar{\Delta} is achieved at p=x2p=x_{2}.

From Fig. 4 (b), we can see that in the sub-case (b), κ(p)<0\kappa(p)<0 holds for p[0,1]p\in[0,1]. Therefore, the function Δ¯\bar{\Delta} is always decreasing for p[0,1]p\in[0,1], which indicates that the value of Δ¯\bar{\Delta} is minimized at p=1p=1.

2) When χ>0\chi>0 and ψ<0\psi<0, we always have x1<0x_{1}<0 and x2>0x_{2}>0. Since pp is in the range of [0,1][0,1], we also have two sub-cases: (a) x2<1x_{2}<1 and (b) x2>1x_{2}>1. The possible shapes of κ(p)\kappa(p) versus pp for the sub-case (a) and sub-case (b) are depicted as Fig. 5 (a) and Fig. 5 (b), respectively.

Similar to case 1), we can prove that the minimum value of Δ¯\bar{\Delta} is achieved by setting p=x2p=x_{2} and p=1p=1 for the case x2<1x_{2}<1 and x2>1x_{2}>1, respectively.

3) When χ<0\chi<0 and ψ<0\psi<0, we always have x1<x2<0x_{1}<x_{2}<0. Therefore, there only exists one possible shape of κ(p)\kappa(p) versus pp in this case, which is depicted as Fig. 6.

From Fig. 6, it is straightforward to see that κ(0)<0\kappa(0)<0 holds for p[0,1]p\in[0,1]. This means that the function Δ¯\bar{\Delta} is always decreasing for p[0,1]p\in[0,1] in this case. As a result, the minimum value of Δ¯\bar{\Delta} is achieved at p=1p=1.

4) When χ<0\chi<0 and ψ>0\psi>0, we always have x2>x1>0x_{2}>x_{1}>0. Therefore, the possible shapes of κ(p)\kappa(p) versus pp can be depicted as Fig. 7. According to (30b), after some algebra manipulations, ψ\psi can be rewritten as

ψ=2P3(P1+P2P1P2)(P1P1P3P2P3+P1P2P3).\psi\!=-2P_{3}(P_{1}\!+\!P_{2}\!-\!P_{1}P_{2})(P_{1}\!-\!P_{1}P_{3}\!-\!P_{2}P_{3}\!+\!P_{1}P_{2}P_{3}). (34)

As ψ>0\psi>0, 0<P1<P2<10<P_{1}<P_{2}<1 and 0<P1<P3<10<P_{1}<P_{3}<1, we thus have P1P1P3P2P3+P1P2P3<0P_{1}-P_{1}P_{3}-P_{2}P_{3}+P_{1}P_{2}P_{3}<0, which can be rewritten as

P1<P2P31P3+P2P3.P_{1}<\frac{P_{2}P_{3}}{1-P_{3}+P_{2}P_{3}}. (35)

To further validate this case, we investigate the value of χ-\chi under the constraint of (35). We rewrite χ-\chi as a quadratic function of P1P_{1}, which can be expressed as

χ(P1)=λP12+μP1+ξ,-\chi(P_{1})=\lambda P_{1}^{2}+\mu P_{1}+\xi, (36)

where

λ=P22P322P2P32+2P2P3P2+P322P3+1,\lambda=P_{2}^{2}P_{3}^{2}-2P_{2}P_{3}^{2}+2P_{2}P_{3}-P_{2}+P_{3}^{2}-2P_{3}+1, (37a)
μ=2P2P32(1P2)+P2(1P3),\mu=2P_{2}P_{3}^{2}(1-P_{2})+P_{2}(1-P_{3}), (37b)
ξ=P2P3(P2P31),\xi=P_{2}P_{3}(P_{2}P_{3}-1), (37c)

are defined for the simplicity of notations.

As 0<P1<P2<10<P_{1}<P_{2}<1 and 0<P1<P3<10<P_{1}<P_{3}<1, according to (37b) and (37c), we have μ>0\mu>0 and ξ<0\xi<0. In addition, if we set P1=1P_{1}=1, after some manipulations, we can attain that χ(P1)=(P31)2>0-\chi({P_{1}})=(P_{3}-1)^{2}>0. Similarly, if we substitute P1=P2P31P3+P2P3P_{1}=\frac{P_{2}P_{3}}{1-P_{3}+P_{2}P_{3}} into the expression of χ(P1)-\chi(P_{1}), we can have χ(P1)=P2P3(P21)(P31)2(P2P3P3+1)2<0-\chi(P_{1})=\frac{P_{2}P_{3}(P_{2}-1)(P_{3}-1)^{2}}{(P_{2}P_{3}-P_{3}+1)^{2}}<0. Based on the above analysis, we can draw the possible shapes of χ(P1)-\chi(P_{1}) versus P1P_{1} for λ>0\lambda>0 and λ<0\lambda<0 in Fig. 8 (a) and Fig. 8 (b), respectively. Note that in these figures, we also show the possible positions of the point with X-coordinate being P2P31P3+P2P3\frac{P_{2}P_{3}}{1-P_{3}+P_{2}P_{3}}.

Refer to caption
(a) When x2<1x_{2}<1
Refer to caption
(b) When x2>1x_{2}>1
Figure 4: Possible shapes for the function κ(p)\kappa(p) versus pp when χ>0\chi>0 and ψ>0\psi>0.
Refer to caption
(a) When x2<1x_{2}<1
Refer to caption
(b) When x2>1x_{2}>1
Figure 5: Possible shapes for the function κ(p)\kappa(p) versus pp when χ>0\chi>0 and ψ<0\psi<0.

Refer to caption

Figure 6: Possible shapes for the function κ(p)\kappa(p) versus pp when χ<0\chi<0 and ψ<0\psi<0.
Refer to caption
(a) When x2<1x_{2}<1
Refer to caption
(b) When x2>1x_{2}>1
Figure 7: Possible shapes for the function κ(p)\kappa(p) versus pp when χ<0\chi<0 and ψ>0\psi>0.
Refer to caption
(a) When λ>0\lambda>0
Refer to caption
(b) When λ<0\lambda<0
Figure 8: Possible shapes for the function χ(P1)-\chi(P_{1}) versus P1P_{1}.

From Fig. 8 (a) and Fig. 8 (b), we can see that if P1<P2P31P3+P2P3P_{1}<\frac{P_{2}P_{3}}{1-P_{3}+P_{2}P_{3}} and 0<P1<10<P_{1}<1, χ<0-\chi<0 always holds. This means that we always have χ>0\chi>0 if ψ>0\psi>0. This is contrary to the initial assumptions of this case (i.e., χ<0\chi<0 and ψ>0\psi>0). Therefore, the case 4) is not valid.

In conclusion, we can find that in all of the possible case 1), case 2) and case 3), the minimum value of Δ¯\bar{\Delta} is achieved by setting p=x2p=x_{2} and p=1p=1 when x2<1x_{2}<1 and x2>1x_{2}>1, respectively. Therefore, based on (30a), (30b), (30c) and (33b), after some algebra manipulations, the optimal pp that minimizes Δ¯\bar{\Delta} can be given by (27).

This completes the proof of Theorem 1.

References

  • [1] M. A. Abd-Elmagid, N. Pappas and H. S. Dhillon, “On the Role of Age of Information in the Internet of Things,” IEEE Commun. Mag., vol. 57, no. 12, pp. 72-77, December 2019.
  • [2] Y. Sun, I. Kadota, R. Talak, E. Modiano, and R. Srikant, “Age of Information: A New Metric for Information Freshness”, in Synthesis Lectures on Communication Networks, Morgan &\& Claypool Publishers, 2019.
  • [3] S. Kaul, R. Yates and M. Gruteser, “Real-time status: How often should one update?,” in Proc. IEEE INFOCOM, 2012, pp. 2731-2735.
  • [4] A. Kosta, N. Pappas, V. Angelakis, “Age of Information: A New Concept, Metric, and Tool,” Found. Trends Netw., vol. 12, no. 3, pp. 162-259, 2017.
  • [5] S. Kaul, M. Gruteser, V. Rai and J. Kenney, “Minimizing age of information in vehicular networks,” in Proc. SECON, 2011, pp. 350-358.
  • [6] Y. Gu, H. Chen, Y. Zhou, Y. Li and B. Vucetic, “Timely Status Update in Internet of Things Monitoring Systems: An Age-Energy Tradeoff,” IEEE Internet Things J., vol. 6, no. 3, pp. 5324-5335, June 2019.
  • [7] M. Costa, M. Codreanu and A. Ephremides, “On the Age of Information in Status Update Systems With Packet Management,” IEEE Trans. Inf. Theory, vol. 62, no. 4, pp. 1897-1910, April 2016.
  • [8] Y. Gu, H. Chen, C. Zhai, Y. Li and B. Vucetic, “Minimizing Age of Information in Cognitive Radio-Based IoT Systems: Underlay or Overlay?,” IEEE Internet Things J., vol. 6, no. 6, pp. 10273-10288, Dec. 2019.
  • [9] Q. Wang, H. Chen, Y. Li, Z. Pang and B. Vucetic, “Minimizing Age of Information for Real-Time Monitoring in Resource-Constrained Industrial IoT Networks,” arXiv preprint arXiv:1912.07186, 2019.
  • [10] Y. Sun, E. Uysal-Biyikoglu, R. D. Yates, C. E. Koksal and N. B. Shroff, “Update or Wait: How to Keep Your Data Fresh,” IEEE Trans. Inf. Theory , vol. 63, no. 11, pp. 7492-7508, Nov. 2017.
  • [11] R. Wang, Y. Gu, H. Chen, Y. Li, and B. Vucetic, “On the Age of Information of Short-Packet Communications with Packet Management,” arXiv preprint arXiv:1908.05447, 2019.
  • [12] I. Kadota, A. Sinha and E. Modiano, “Optimizing Age of Information in Wireless Networks with Throughput Constraints,” in IEEE INFOCOM, 2018, pp. 1844-1852.
  • [13] Q. Wang, H. Chen, Y. Gu, Y. Li and B. Vucetic, “Minimizing the Age of Information of Cognitive Radio-Based IoT Systems Under A Collision Constraint,” arXiv preprint arXiv: 2001.02482, 2020.
  • [14] A. A. Simiscuka and G. Muntean, “Age of Information as a QoS Metric in a Relay-Based IoT Mobility Solution,” in Proc. 14th Int. Wireless Commun. Mobile Comput. Conf. (IWCMC), 2018, pp. 868-873.
  • [15] A. Arafa and S. Ulukus, “Timely Updates in Energy Harvesting Two-Hop Networks: Offline and Online Policies,” IEEE Trans. Wireless Commun., vol. 18, no. 8, pp. 4017-4030, Aug. 2019.
  • [16] R. Talak, S. Karaman and E. Modiano, “Minimizing age-of-information in multi-hop wireless networks,” in Proc. 55th Annual Allerton Conf. Commun., Control, Comput. (Allerton), 2017, pp. 486-493.
  • [17] A. Maatouk, M. Assaad and A. Ephremides, “The Age of Updates in a Simple Relay Network,” in Proc. IEEE Inf. Theory Workshop (ITW), 2018, pp. 1-5.
  • [18] J. N. Laneman, D. N. C. Tse and G. W. Wornell, “Cooperative diversity in wireless networks: Efficient protocols and outage behavior,” IEEE Trans. Inf. Theory, vol. 50, no. 12, pp. 3062-3080, Dec. 2004.
  • [19] A. K. Sadek, K. J. R. Liu and A. Ephremides, “Cognitive multiple access via cooperation: Protocol design and performance analysis,” IEEE Trans. Inf. Theory, vol. 53, no. 10, pp. 3677-3696, Oct. 2007.
  • [20] A. Jeffrey and D. Zwillinger, Table of Integrals, Series, and Products. Amsterdam, The Netherlands: Elsevier, 2007.
  • [21] K. Chen and L. Huang,“Age-of-information in the presence of error,” in Proc. IEEE Int. Symp. Inf. Theory (ISIT), 2016, pp. 2579-2583.
  • [22] E. T. Ceran, D. Gündüz and A. György, “Average Age of Information With Hybrid ARQ Under a Resource Constraint,” IEEE Trans. Wireless Commun., vol. 18, no. 3, pp. 1900-1913, March 2019.