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

Age of Information with Hybrid-ARQ: A Unified Explicit Result

Aimin Li ID , Graduate Student Member, IEEE, Shaohua Wu ID , Member, IEEE, Jian Jiao ID , Member, IEEE, Ning Zhang ID , Senior Member, IEEE, and Qinyu Zhang ID Senior Member, IEEE. This work was supported in part by the National Natural Science Foundation of China under Grant nos. 61871147, 61831008, 91638204, and in part by the Shenzhen Municipal Science and Technology Plan under Grant nos. JCYJ20170811160142808, JCYJ20170811154309920. (Corresponding author: Shaohua Wu.) A. Li, S. Wu, J. Jiao and Q. Zhang are with the Department of Electronic Engineering, Harbin Institute of Technology (Shenzhen), Guangdong, China. N. Zhang is with the Department of Electrical and Computer Engineering, University of Windsor, Windsor, ON, N9B 3P4, Canada (e-mail: [email protected]; [email protected]; [email protected]; [email protected]; [email protected]).
Abstract

Delivering timely status updates in a timeliness-critical communication system is of paramount importance to assist accurate and efficient decision making. Therefore, the topic of analyzing Age of Information (AoI) has aroused new research interest. This paper contributes to new results in this area by systematically analyzing the AoI of two types of Hybrid Automatic Repeat reQuest (HARQ) techniques that have been newly standardized in the Release-16 5G New Radio (NR) specifications, namely reactive HARQ and proactive HARQ. Under a code-based status update system with non-trivial coding delay, transmission delay, propagation delay, decoding delay, and feedback delay, we derive unified closed-form average AoI and average Peak AoI expressions for reactive HARQ and proactive HARQ, respectively. Based on the obtained explicit expressions, we formulate an AoI minimization problem to investigate the age-optimal codeblock assignment strategy in the finite block-length (FBL) regime. Through case studies and analytical results, we provide comparative insights between reactive HARQ and proactive HARQ from a perspective of freshness of information. The numerical results and optimization solutions show that proactive HARQ draws its strength from both age performance and system robustness, thus enabling the potential to provide new system advancement of a freshness-critical status update system.

Index Terms:
5G NR, proactive HARQ, reactive HARQ, age of information, finite blocklength, real-time status update, low latency.

I Introduction

I-A Background

In recent ten years since Kaul et al. proposed a framework to quantify the timeliness of information in 2012 [1], one of the most popular ideas in timely update system design has been how to keep information as fresh as possible and ensure timely information delivery. For timely update systems such as vehicle networks where the vehicle’s velocity and location are disseminated to ensure safe transportation [2], environmental sensor networks where the updates of a time-varying phenomenon are collected for large-scale monitoring [3], and wireless communication networks where adaptive scheduling algorithms are adopted based on the time-varying channel state information [4], achieving timely delivery can freshen the monitor’s awareness of the sources and thus assist correct and efficient decision making.

This has aroused new interest in analyzing Age of Information (AoI) performance metrics. AoI has been broadly used to capture the freshness of a monitor’s knowledge of an entity or process. Different from conventional performance metrics such as delay and throughput, AoI comprehensively measures the effects of update rate, latency, and system utilization. Initial works on this issue were mainly based on queue analysis, which originated from single-source single-server queues [5, 6, 7, 1], and subsequently developed to multiple-sources single-server queues [8, 9, 10, 11] and wireless queuing networks [12, 13, 14, 15, 16, 17]. These works are based on an ideal assumption that the status update is transmitted through a perfect channel without packet errors and losses. In practice, however, packet errors and losses are inevitable due to ubiquitous noises, signal interference, and channel fading. As the incorrectly decoded message does not bring about fresh awareness, the packet errors and losses will result in staleness of information, leading to uncontrollable residual errors, system instability, and wrong decisions. Therefore, it is imperative to analyze the AoI over unreliable channels.

I-B Related Works

Some recent works have noticed the above limitation and have extended the AoI analyses to the physical (PHY) layer. One pioneering work concerning this issue was accomplished by Chen, et al. in 2016 [18], in which the update is delivered over an erasure channel and the Peak Age of Information (PAoI) is studied. This work has aroused extensive research interest in understanding the effect of system reliability on AoI. From then on, including but not limited to the follow-up works that also analyzed AoI over the erasure channel [19, 20, 21], various transmission protocols, ranging from conventional protocols like non-ARQ, classical ARQ and truncated ARQ protocols, to state-of-the-art protocols such as HARQ with Chase Combing (HARQ-CC) and HARQ with Incremental Redundancy (HARQ-IR) protocols, have been investigated under different types of noisy channels [22, 23, 24, 25, 26].

We notice that the above AoI analyses focus on the transmission delay, and neglect other types of system delay such as coding delay, propagation delay, decoding delay and feedback delay. An exception work is [27], which considers non-trivial propagation delay and studies the AoI of HARQ-IR with a fixed number of retransmitted packets m=2m=2 under Satellite-IoT Systems, but also assumes negligible coding delay and decoding delay. Nevertheless, in practical communication systems, especially the short-packet communication, the coding delay and decoding delay are also nontrivial compared to the transmission delay, resulting in the staleness of information by nature. Thus, we focus on a more realistic (or general) scenario where different types of delay elements naturally exist and the number of retransmitted packets is not fixed to m=2m=2. In this regard, we would like to provide a basic framework to comprehensively study the trade-off among coding complexity, decoding complexity, code length, number of retransmitted packets and error probability from the AoI perspective.

Up to this point, we have only introduced AoI research based on conventional reactive HARQ (also known as stop-and-wait HARQ), which allows for retransmissions only upon the reception of a Negative ACKnowledgment (NACK). As such, the retransmission process is not truly automatic. In the Release-16 5G NR specifications by the 3rd3^{rd} Generation Partnership Project (3GPP), a new HARQ protocol named proactive HARQ is designated for the Up-Link Grant-Free communication to enable the potential for meeting the stringent requirements for URLLC [28]. Some recent works have shown the superiority of proactive HARQ in terms of latency and throughput compared to reactive HARQ [29, 30, 31]. Inspiringly, these available studies also witness the potential for proactive HARQ to be applied in the freshness-critical status update system. To this end, we would like to theoretically analyze the AoI performance of reactive HARQ and proactive HARQ to investigate whether proactive HARQ will facilitate timeliness of information in the freshness-critical status update system.

TABLE I: Construsting The Novelty Of Our Work To The Literature
Contributions This Work [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17] [18, 19, 20, 22, 23] [25, 26] [21, 24] [27]
Age of Information (AoI) \checkmark \checkmark \checkmark \checkmark \checkmark \checkmark
Finite Block-Length Regime \checkmark \checkmark
Reactive Hybrid-ARQ \checkmark \checkmark \checkmark \checkmark \checkmark
Flexible Number of Retransmissions \checkmark \checkmark
Effect of Delay Elements
Other Than Transmission Delay
\checkmark \checkmark
Proactive Hybrid-ARQ \checkmark

I-C Contributions

The research on the HARQ-based timely status update system is still in the ascendant, and some open issues remain to be addressed. First, there have been a lot of works providing explicit average age results under different types of protocols and systems. Examples include the average AoI expressions under fixed-length non-ARQ protocols, truncated-ARQ, classical ARQ, and the explicit results of some advanced ARQ-based techniques like HARQ-CC and HARQ-IR. However, there has not been a unified expression that can unify the aforementioned expressions in a single closed-form formula. By providing such a unified result, the comparative insights and the intrinsic relationships among different protocols will be further investigated. Second, the existing literature only considers certain types of delay in the status update system and assumes others to be negligible. However, as the delay exists by nature and plays as a critical part in affecting the freshness of information, to comprehensively consider the coding delay (or processing delay), propagation delay, transmission delay, decoding delay and feedback delay in the status update system and provide a unified closed-form result will provide a systematic understanding in analyzing the age of a realistic freshness-critical status update system. Third, the age performance has been extensively studied over erasure channels. However, little research considers the short-packet AoI analysis over the AWGN channel. Finally, the majority of existing works mainly focus on AoI analysis of conventional reactive HARQ. Some recent works analyzing the performance of proactive HARQ are based on some conventional performance metrics, such as throughput and latency. Thus, to analyze the AoI of proactive HARQ will fill this research gap and may further facilitate new system advancement of a status update system.

Motivated by the above, this work achieves several key contributions and we summarize them as follows:

  • We derive unified closed-form average AoI and average Peak AoI expressions for reactive HARQ, wherein: ii) different kinds of delay elements (i.e., coding delay, transmission delay, propagation delay, decoding delay, and feedback delay) are comprehensively considered; iiii) the number of repeated packets is not fixed as m=2m=2, but is relaxed to a variable value; iiiiii) different types of protocols are unified to a single expression.

  • We investigate the AoI explicit expressions and comparative insights for proactive HARQ, which is the first work analyzing proactive HARQ from the AoI perspective. Theoretical and numerical comparisons are given to show the superiority of proactive HARQ in enabling timely information delivery.

  • We also try to further optimize the AoI for both reactive HARQ and proactive HARQ. By formulating an AoI minimization problem in the FBL regime, we solve out the age-optimal block assignment strategy for reactive HARQ and proactive HARQ, respectively. The results show that the optimal strategy for proactive HARQ turns out to be the finest grained symbol-by-symbol transmission, while that for reactive HARQ is highly dependent on the propagation delay and SNR.

I-D Organization

The rest of this paper is organized as follows. In Section II, we briefly introduce the considered system model. The generalized closed-form expressions of average AoI and average Peak AoI for reactive HARQ and proactive HARQ are provided in Section III, where the effect of different delay elements is added into the analysis. In Section IV, we design an optimization problem to reduce the average AoI of reactive HARQ and proactive HARQ, respectively. Numerical results and discussions are given in Section V, followed by conclusions in Section VI.

II System Model

We consider an end-to-end (E2E) code-based timely status update system. The update generator (source) is monitoring a time-varying phenomenon (t){0,,2k1}\mathcal{F}\left(t\right)\in\left\{0,\cdots,2^{k}-1\right\}, where the time tt is divided into some time slots in units of channel use such that tt\in\mathbb{N}111Here we consider the symbol-level AoI analysis. Some recent works focusing on PHY-layer AoI analysis also discretize the time into time slots to analyze symbol-level AoI [19, 21, 23, 22].. We assume that the monitored phenomenon is modeled as a sequence of independent and uniformly distributed symbols. In such a case, the size of the generated observation is kk information bits. The monitored data is transmitted through a noisy channel to a central location. We use the notation \mathbb{N} for non-negative integers and the notation +\mathbb{Z}^{+} for positive integers. Also, we define the notation [m]\left[m\right] as [m]{1,2,,m}\left[m\right]\triangleq\left\{1,2,\cdots,m\right\} for any positive integer m+m\in\mathbb{Z}^{+}.

Refer to caption
Figure 1: A general HARQ-based real-time status update system.

II-A Channel Model

We consider an E2E communication setup leveraging a power-limited AWGN model:

Y=PX+Z,Y=\sqrt{P}X+Z, (1)

where PP is the average transmit power, XX is the unit-variance coded symbol and Z𝒩(0,1)Z\sim\mathcal{N}\left(0,1\right) is the independent and identically distributed (i.i.d) AWGN.

Refer to caption
(a) Reactive HARQ
Refer to caption
(b) Proactive HARQ
Figure 2: Instantaneous age evolutions of reactive HARQ and proactive HARQ. Here the maximum retransmissions (the number of sub-codeblocks) is set as m=4m=4. The length of the yellow block represents the coding delay τc\tau_{\rm c}, the length of the green block denotes the decoding delay τd\tau_{\rm d}, the length of the blue block represents the transmission delay i\ell_{i}, the length of the oblique arrow projective on the timeline describes the propagation delay τp\tau_{\rm p} (or feedback delay τf\tau_{\rm f}).
Remark 1.

We notice that the AoI analysis over erasure channels has been extensively studied in the existing literature. However, the performance over AWGN channels is still not clear so far. As such, we consider the AWGN channel in this paper to reveal the AoI performance over the AWGN channel.

II-B Hybrid ARQ

The overall system model is shown in Fig. 1. The considered system is in close-up fashion with perfect HARQ feedback222In this article, we assume that the feedback is error-free. The research with erroneous feedback can be extended following this work.. At the transmitter end, the update generator (source) generates a kk-bit short-packet update uu and encodes it to a parent codeword with length i=1mi\sum\nolimits_{i=1}^{m}{{\ell_{i}}} channel uses, which is then divided into mm sub-codeblocks with length i,i[m]\ell_{i},i\in\left[m\right] and stored in a buffer waiting to be transmitted. The coding process above will take up τc\tau_{\rm c} channel uses, and we call it as coding delay. Next, the stored sub-codeblocks are transmitted over a noisy channel sub-codeblock by sub-codeblock, with each transmission taking a transmission delay i,i[m]\ell_{i},i\in\left[m\right] channel uses. The transmitted sub-blocks will take τp\tau_{\rm p} channel uses to arrive at the receiver end. At the receiver end, we assume that the decoding process is conducted once receiving any complete sub-codeblock. The decoding delay in each transmission round is assumed to be the same and is denoted by τd\tau_{\rm d}. If the update is decoded correctly such that u^=u\hat{u}=u, an ACKnowledgment (ACK) will be fed back to the transmitter; otherwise, a NACK will be sent back. The feedback, similar to the forwarding information propagation, generally takes time and results in delay by nature, and we denote the delay as τf\tau_{\rm f} channel uses.

The generate-at-will model is adopted in the considered E2E status update system. That is, when the transmitter receives an ACK, the process of sensing and sampling will be performed, and a new update will be generated. In such a case, we mainly focus on two types of HARQ schemes: reactive HARQ and proactive HARQ. The detailed processes are shown in Fig. 2(a) and Fig. 2(b), respectively.

II-B1 Reactive HARQ

Reactive HARQ is also know as stop-and-wait HARQ. In Fig. 2(a), we demonstrate a detailed stop-and-wait retransmission process of reactive HARQ, wherein the maximum number of sub-codeblocks (or maximum retransmissions) is set as m=4m=4. The so-called reactive scheme implies that the transmitter allows for retransmissions only upon the reception of a NACK. As such, the transmitter should always wait for a feedback to decide whether to generate a new update or retransmit the old update’s sub-codeblocks. The waiting time, however, is referred ro as the HARQ round trip time (RTT) and will result in additional latency. Therefore, the reactive HARQ scheme allows for only a limited number of retransmissions in the URLLC application scenarios and thereby enables great potential to be further advanced [31].

II-B2 Proactive HARQ

The proactive HARQ scheme with maximum sub-codeblocks m=4m=4 is shown in Fig. 2(b). As its name indicates, the retransmission process is completely spontaneous and proactive, which is interrupted only when an ACK is received. The core idea of proactive HARQ is to eliminate the need for waiting for a feedback and implement consecutive retransmitting. By proactive retransmitting, the latency introduced by waiting for a feedback is reduced, and thus the issue of long HARQ round trip time (RTT) is resolved.

II-C Performance Metric

We focus on AoI analysis and optimization in this paper. Here we simply review the definition of instantaneous AoI as in Definition 1. For more intuitive and visualized results, Fig. 2 also gives the instantaneous age evolutions of reactive HARQ and proactive HARQ respectively.

Definition 1.

(AoI) Denote tiGt_{i}^{\rm G} as the generation time instant of the ithi^{\rm th} status update packet that can be correctly decoded, and denote tiSt_{i}^{\rm S} as the time instant at which this packet is correctly decoded. At a time instant τ\tau, the index of the most recently generated update can be given by N(τ)=min{i|tiS>τ}N(\tau)=\min\left\{i|t_{i}^{\rm S}>\tau\right\} and the time stamp is U(τ)=tN(τ)GU(\tau)=t_{N(\tau)}^{\rm G}. Then, the instantaneous AoI is defined as Δ(t)tU(t)\Delta\left(t\right)\triangleq t-U\left(t\right).

II-C1 Average AoI

Also known as time-average AoI, Average AoI is a statistical metric that measures the long-term average age of a status update system. In our considered discrete symbol-level system where the time is divided into some time slots in units of channel use, the average AoI is defined as follows.

Definition 2.

(Average AoI) The average AoI of a real-time status update system is defined as:

Δ¯limN1Nt=1NΔ(t).\bar{\Delta}\triangleq\mathop{\lim}\limits_{N\to\infty}\frac{1}{N}\sum\nolimits_{t=1}^{N}{\Delta\left(t\right)}. (2)

II-C2 Average Peak AoI

We also provide explicit expressions for the average Peak AoI in this paper. The peak age indicates the maximum value of age in each renewal process. In our considered system, the average Peak AoI is defined as follows.

Definition 3.

(Average Peak AoI) The average Peak AoI of a real-time status update system is

Δ¯PlimN1Nj=1NΔ(tjS1).\bar{\Delta}^{\rm P}\triangleq\lim\limits_{N\to\infty}\frac{1}{N}\sum_{j=1}^{N}\Delta\left(t^{S}_{j}-1\right). (3)

III Analytical Results

In this section, we study the symbol-level AoI of reactive HARQ and proactive HARQ. We first give the closed-form expressions for the AoI in Proposition 1 and Proposition 3, and then conduct a theoretical AoI comparison between the two considered transmission protocols in Corollary 1. The AoI expressions, given in (14) and (20), are functions of the block assignment vector 𝐧\bf n and its dependent error probability vector 𝐞\bf e, where the element nin_{i} in vector 𝐧\bf n denotes the number of cumulative transmitted symbols up to the ithi^{\rm th} transmission round with ni=j=1ijn_{i}=\sum\nolimits_{j=1}^{i}{{\ell_{j}}}, and the element ϵi\epsilon_{i} in vector 𝐞\bf e denotes the probability that the ithi^{\rm th} re-transmitted message remains incorrectly decoded.

By flexible choices of the vector 𝐧\bf n and the vector 𝐞\bf e, we also demonstrate that our derived expression for reactive transmission protocol also unifies the available AoI analyses in the existing literature. Moreover, by using the result of the achievable rate of finite-length codes, we can obtain the AoI closed-form expression under the finite block-length (FBL) regime.

III-A Reactive Scheme

III-A1 Average AoI

Denote the generation time of the jthj^{\rm th} collected message as tjt_{j}, and denote the jthj^{\rm th} collected message as Mj(tj)M_{j}\triangleq\mathcal{F}\left(t_{j}\right). The jthj^{\rm th} collected message MjM_{j} is encoded to a parent code 𝒞(Mj)\mathcal{C}\left(M_{j}\right) with size nm=a[m]an_{m}=\sum\nolimits_{a\in[m]}\ell_{a} being stored in a transmission buffer. Then, the transmission is evoked round by round until all the symbols stored in the buffer is transmitted or an ACK is received. In the ithi^{\rm th} round of transmission, the transmitter will transmit i\ell_{i} symbols, and the decoder will leverage the cumulatively received ni=a[i]an_{i}={\sum\nolimits_{a\in[i]}}\ell_{a} symbols to decode the message.

In such a case, we introduce ςj,i{0,1}\varsigma_{j,i}\in\left\{0,1\right\} to denote the feedback signal in the ithi^{\rm th} transmission round when transmitting message MjM_{j}: If ςj,i=1\varsigma_{j,i}=1 or i=mi=m, the transmitter will no longer transmit message MjM_{j}, instead, it will collect new update Mj+1M_{j+1} for transmission; If ςj,i=0\varsigma_{j,i}=0 and i<mi<m, the source will transmit additional i+1\ell_{i+1} encoded symbols from 𝒞(Mj)\mathcal{C}\left(M_{j}\right). As such, the probability that the ithi^{\rm th} transmission is not correctly decoded is given as ϵi=1𝔼ςj,i\epsilon_{i}=1-\mathbb{E}\varsigma_{j,i}.

Let Qjinf{i>Qj1:ςi,m=1},Q0=0Q_{j}\triangleq\inf\left\{i>Q_{j-1}:\varsigma_{i,m}=1\right\},Q_{0}=0 be the cumulative number of generated packets until the jthj^{\rm th} decoding success. Let RjQjQj11R_{j}\triangleq Q_{j}-Q_{j-1}-1 represent the decoding failures between two successful decoding and Vjinf{i[m]:ςQj,i=1}V_{j}\triangleq\inf\left\{i\in[m]:\varsigma_{Q_{j},i}=1\right\} denote the round in which the QjthQ_{j}^{\rm th} packet gets decoded, we have that

Lemma 1.

The random sequences RjR_{j} and VjV_{j} are independent, and they are i.i.d with distributions

(Rj=a)=(1ϵm)ϵma,a,\displaystyle\mathbb{P}\left(R_{j}=a\right)=\left(1-\epsilon_{m}\right)\epsilon^{a}_{m},\quad a\in\mathbb{N}, (4)
(Vj=i)=ϵi1ϵi1ϵm,i[m].\displaystyle\mathbb{P}\left(V_{j}=i\right)=\frac{\epsilon_{i-1}-\epsilon_{i}}{1-\epsilon_{m}},\quad i\in[m].
Proof.

Please refer to Appendix A. ∎

Lemma 2.

The first and second moments for RjR_{j} are

𝔼Rj\displaystyle\mathbb{E}R_{j} =ϵm1ϵm,\displaystyle=\frac{\epsilon_{m}}{1-\epsilon_{m}}, (5)
𝔼Rj2\displaystyle\mathbb{E}R_{j}^{2} =ϵm2+ϵm(1ϵm)2.\displaystyle=\frac{\epsilon_{m}^{2}+\epsilon_{m}}{\left(1-\epsilon_{m}\right)^{2}}.
Proof.

By utilizing the distributions in (4) and the definitions that 𝔼Rj=aa(Rj=a)\mathbb{E}R_{j}=\sum_{a\in\mathbb{N}}a\cdot\mathbb{P}\left(R_{j}=a\right) and 𝔼Rj2=aa2(Rj=a)\mathbb{E}R_{j}^{2}=\sum_{a\in\mathbb{N}}a^{2}\cdot\mathbb{P}\left(R_{j}=a\right), we obtain the results in (5). ∎

We then let τiReac\tau^{\rm Reac}_{i} represent the elapsed time form generating an update to receiving its ithi^{\rm th} feedback signal for reactive HARQ. For reactive HARQ, the transmitter should wait for an integral RTT to receive a feedback. As such, we have τiReac=ni+τc+i(τd+τf+τp)\tau^{\rm Reac}_{i}=n_{i}+\tau_{\rm c}+i\left(\tau_{\rm d}+\tau_{\rm f}+\tau_{\rm p}\right). Denote 𝒯=τd+τf+τp\mathcal{T}=\tau_{\rm d}+\tau_{\rm f}+\tau_{\rm p}, the first and second moments for τVjReac\tau^{\rm Reac}_{V_{j}} is given as follows.

Lemma 3.

The first and second moments for τVjReac\tau^{\rm Reac}_{V_{j}} are

𝔼τVjReac\displaystyle\mathbb{E}\tau^{\rm Reac}_{V_{j}} =n1+τc+𝒯+i[m1](ni+1ni+𝒯)ϵi1ϵm,\displaystyle=n_{1}+\tau_{\rm c}+\mathcal{T}+\sum_{i\in[m-1]}\frac{\left(n_{i+1}-n_{i}+\mathcal{T}\right)\epsilon_{i}}{1-\epsilon_{m}}, (6)
𝔼(τVjReac)2\displaystyle\mathbb{E}\left(\tau^{\rm Reac}_{V_{j}}\right)^{2} =(n1+τc+𝒯)2+i[m1](ni+1ni+𝒯)(ni+1+ni+2τc+(2i+1)𝒯)ϵi1ϵm.\displaystyle=\left(n_{1}+\tau_{\rm c}+\mathcal{T}\right)^{2}+\sum_{i\in[m-1]}\frac{\left(n_{i+1}-n_{i}+\mathcal{T}\right)\left(n_{i+1}+n_{i}+2\tau_{\rm c}+\left(2i+1\right)\mathcal{T}\right)\epsilon_{i}}{1-\epsilon_{m}}.
Proof.

By utilizing the distributions in (4) and the definitions that 𝔼τVjReac=i[m]τiReac(Vj=i)\mathbb{E}\tau^{\rm Reac}_{V_{j}}=\sum_{i\in[m]}\tau^{\rm Reac}_{i}\cdot\mathbb{P}\left({V_{j}}=i\right) and 𝔼(τVjReac)2=i[m](τiReac)2(Vj=i)\mathbb{E}\left(\tau^{\rm Reac}_{V_{j}}\right)^{2}=\sum_{i\in[m]}\left(\tau^{\rm Reac}_{i}\right)^{2}\cdot\mathbb{P}\left({V_{j}}=i\right), we obtain the results in (6). ∎

With these notations, we can recursively write the time-instant of the jthj^{\rm th} successful decoding tjSt_{j}^{S} as follows.

tjS=tj1S+τmReacRj+τVjReac+τf.t_{j}^{S}=t_{j-1}^{S}+\tau^{\rm Reac}_{m}R_{j}+\tau^{\rm Reac}_{V_{j}}+\tau_{\rm f}. (7)

Therefore, the interval between the (j1)th{\left(j-1\right)}^{\rm th} and the jthj^{\rm th} successful decoding for reactive HARQ is given as follows.

TjReac\displaystyle T^{\rm Reac}_{j} =tjStj1S=τmReacRj+τVjReac+τf.\displaystyle=t_{j}^{S}-t_{j-1}^{S}=\tau^{\rm Reac}_{m}R_{j}+\tau^{\rm Reac}_{V_{j}}+\tau_{\rm f}. (8)

Since both the random sequences RjR_{j} and VjV_{j} are i.i.d. and have finite first and second moments, it turns out that the sequence TjReacT^{\rm Reac}_{j} is also i.i.d. with finite finite first and second moments. We give the first and second moments of TjReacT^{\rm Reac}_{j} as follows.

Lemma 4.

The first and second moments for TjReacT^{\rm Reac}_{j} are

𝔼TjReac=n1+τc+𝒯+i[m1](ni+1ni+𝒯)ϵi1ϵm,\displaystyle\mathbb{E}T^{\rm Reac}_{j}=\frac{n_{1}+\tau_{\rm c}+\mathcal{T}+\sum_{i\in[m-1]}{\left(n_{i+1}-n_{i}+\mathcal{T}\right)\epsilon_{i}}}{1-\epsilon_{m}}, (9)
𝔼(TjReac)2=(nm+τc+m𝒯)2(1+ϵm)(1ϵm)2\displaystyle\mathbb{E}\left(T^{\rm Reac}_{j}\right)^{2}=\frac{\left(n_{m}+\tau_{\rm c}+m\mathcal{T}\right)^{2}\left(1+\epsilon_{m}\right)}{\left(1-\epsilon_{m}\right)^{2}}
i[m1]1ϵi1ϵm(ni+1ni+𝒯)[ni+ni+1+2τc+(2i+1)𝒯+2(nm+τc+m𝒯)ϵm1ϵm].\displaystyle-\sum_{i\in[m-1]}\frac{1-\epsilon_{i}}{1-\epsilon_{m}}\left(n_{i+1}-n_{i}+\mathcal{T}\right)\left[n_{i}+n_{i+1}+2\tau_{\rm c}+\left(2i+1\right)\mathcal{T}+\frac{2\left(n_{m}+\tau_{\rm c}+m\mathcal{T}\right)\epsilon_{m}}{1-\epsilon_{m}}\right].
Proof.

From (8), we obtain

𝔼TjReac=τmReac𝔼Rj+𝔼τVjReac+τf,\displaystyle\mathbb{E}T^{\rm Reac}_{j}=\tau^{\rm Reac}_{m}\mathbb{E}R_{j}+\mathbb{E}\tau^{\rm Reac}_{V_{j}}+\tau_{\rm f}, (10)
𝔼(TjReac)2=𝔼(τmReacRj+τVjReac+τf)2.\displaystyle\mathbb{E}\left(T^{\rm Reac}_{j}\right)^{2}=\mathbb{E}\left(\tau^{\rm Reac}_{m}R_{j}+\tau^{\rm Reac}_{V_{j}}+\tau_{\rm f}\right)^{2}.

Note that the random variables RjR_{j} and Vj{V_{j}} are independent with each other and the first and second moments of RjR_{j} and τVjReac\tau^{\rm Reac}_{V_{j}} have been given in Lemma 2 and Lemma 3. Substitute them into (10), we can obtain the results in (9). ∎

As Definition 1 indicates, the generation time instant U(τ)U\left(\tau\right) is given by U(τ)=tN(τ)GU\left(\tau\right)=t^{G}_{N\left(\tau\right)}. Because N(tj1S)=min{i|tiS>tj1S}=jN\left(t^{S}_{j-1}\right)=\min\left\{i|t_{i}^{\rm S}>t^{S}_{j-1}\right\}=j, we obtain that

U(tj1S)=tjG=tj1S+τfτVj1Reac.U\left(t^{S}_{j-1}\right)=t^{G}_{j}=t^{S}_{j-1}+\tau_{\rm f}-\tau^{\rm Reac}_{V_{j-1}}.

Thus, for any time slots tt in the jthj^{\rm th} renewal interval j{tj1S,,tjS1}\mathcal{I}_{j}\triangleq\left\{t^{S}_{j-1},\cdots,t^{S}_{j}-1\right\}, we have

U(t)=U(tj1S)=tj1S+τfτVj1Reac, for tj.U\left(t\right)=U\left(t^{S}_{j-1}\right)=t^{S}_{j-1}+\tau_{\rm f}-\tau^{\rm Reac}_{V_{j-1}},\text{ for }t\in\mathcal{I}_{j}. (11)

As such, the instantaneous age is given as

Δ(t)\displaystyle\Delta\left(t\right) =tU(t)\displaystyle=t-U\left(t\right) (12)
=ttj1Sτf+τVj1Reac, for tj.\displaystyle=t-t^{S}_{j-1}-\tau_{\rm f}+\tau^{\rm Reac}_{V_{j-1}},\text{ for }t\in\mathcal{I}_{j}.
Lemma 5.

For the considered reactive HARQ model, the average AoI can be calculated by

Δ¯Reactive=𝔼tjΔ(t)𝔼TjReac=𝔼TjProac2𝔼(TjReac)2+𝔼τVjReacτf12.\bar{\Delta}_{\rm Reactive}=\frac{\mathbb{E}\sum_{t\in\mathcal{I}_{j}}\Delta{\left(t\right)}}{\mathbb{E}T^{\rm Reac}_{j}}=\frac{\mathbb{E}T^{\rm Proac}_{j}}{2\mathbb{E}\left(T^{\rm Reac}_{j}\right)^{2}}+\mathbb{E}\tau^{\rm Reac}_{V_{j}}-\tau_{\rm f}-\frac{1}{2}. (13)
Proof.

Please refer to Appendix B. ∎

Then, by adopting the available first and second moments of TjReacT^{\rm Reac}_{j} and τVjReac\tau^{\rm Reac}_{V_{j}} given in (6) and (10), we have the average AoI for reactive HARQ as follows.

Proposition 1.

(The Generalized Closed-form Average AoI Expression for Reactive Scheme) For reactive HARQ with maximum retransmissions mm, block assignment vector 𝐧=(n1,n2,,nm){\bf n}=\left(n_{1},n_{2},\cdots,n_{m}\right) and error probability vector 𝐞=(ϵ1,ϵ2,,ϵm){\bf e}=\left(\epsilon_{1},\epsilon_{2},\cdots,\epsilon_{m}\right), the average AoI can be calculated by

Δ¯Reactive=12τf+τc+n1+𝒯+i=1m1(ni+1ni+𝒯)ϵi1ϵm+\displaystyle{\bar{\Delta}_{\rm Reactive}}=-\frac{1}{2}-{\tau_{\rm f}}+\frac{{{\tau_{\rm c}}+{n_{1}}+\mathcal{T}+\sum\limits_{i=1}^{m-1}{\left({{n_{i+1}}-{n_{i}}+\mathcal{T}}\right){\epsilon_{i}}}}}{{1-{\epsilon_{m}}}}+ (14)
(τc+n1+𝒯)2+i=1m1(ni+1ni+𝒯)(2τc+ni+1+ni+(2i+1)𝒯)ϵi2(τc+n1+𝒯+i=1m1(ni+1ni+𝒯)ϵi),\displaystyle\frac{{{{\left({{\tau_{\rm c}}+{n_{1}}+\mathcal{T}}\right)}^{2}}+\sum\limits_{i=1}^{m-1}{\left({{n_{i+1}}-{n_{i}}+\mathcal{T}}\right)\left(2{{\tau_{\rm c}}+{n_{i+1}}+{n_{i}}+\left({2i+1}\right)\mathcal{T}}\right){\epsilon_{i}}}}}{{2\left({{\tau_{\rm c}}+{n_{1}}+\mathcal{T}+\sum\limits_{i=1}^{m-1}{\left({{n_{i+1}}-{n_{i}}+\mathcal{T}}\right){\epsilon_{i}}}}\right)}},

where 𝒯=τf+τp+τd\mathcal{T}=\tau_{\rm f}+\tau_{\rm p}+\tau_{\rm d} with τc\tau_{\rm c}, τp\tau_{\rm p}, τd\tau_{\rm d} and τf\tau_{\rm f} denoting the coding delay, propagation delay, decoding delay and feedback delay, respectively.

III-A2 Case Study: A Unified Result

With Proposition 1 in hand, we can conduct some case studies by flexibly considering the choices of the block assignment vector 𝐧\bf n and the error probability vector 𝐞\bf e. By this means, we theoretically show that the closed-form AoI expressions given in this paper is a unified result.333The average AoI expressions in Case 1, Case 2 and Case 3 are corresponding to Proposition 1, Proposition 3 and Proposition 2 of [22], respectively. The average AoI expression in Case 4 is a variant of the result in [21].. Though the given examples are not exhaustive in this paper, we can observe from these case studies that the unified expression given in (14) enables potential for exploring the intrinsic relationship and comparative insights among different types of transmission protocols.

Case 1.

(Average AoI for Fixed-rate Codes without ARQ) We show that the available average AoI expression for fixed-rate codes in [22] is a specific case of our unified result in (14). For fixed-rate codes without ARQ, the maximum retransmissions turns to m=1m=1. Substitute m=1m=1 into (14) and remove the effect of delay elements such that τc=τp=τd=τf=0\tau_{\rm c}=\tau_{\rm p}=\tau_{\rm d}=\tau_{\rm f}=0, we can obtain the average AoI as the Proposition 1 in [22]:

Δ¯NonARQ=12+n11ϵ1+n12,{\bar{\Delta}_{\rm Non-ARQ}}=-\frac{1}{2}+\frac{{{n_{1}}}}{{1-{\epsilon_{1}}}}+\frac{{{n_{1}}}}{2},

where n1n_{1} is the code length and ϵ1\epsilon_{1} is the error probability of the fixed-rate codes.

Case 2.

(Average AoI for Truncated ARQ (TARQ)) We demonstrate that the average AoI expression for TARQ is also a specific case of our unified result in (14). For truncated ARQ, the transmitter retransmits the same packet till the allowable maximum retransmissions mm is reached or this packet is successfully received. Since the retransmitted packet is the same as the first packet, the cumulative transmitted message length is ni=in1n_{i}=in_{1} and the corresponding error probability is ϵi=ϵ1i\epsilon_{i}={\epsilon_{1}}^{i}. Then, by substituting them back into (14) and similarly remove the effect of delay elements such that τc=τp=τd=τf=0\tau_{\rm c}=\tau_{\rm p}=\tau_{\rm d}=\tau_{\rm f}=0, we can obtain the average AoI as the Proposition 3 in [22]:

Δ¯TARQ=12+n1(21ϵ112mϵ1m1ϵ1m).{\bar{\Delta}_{\rm TARQ}}=-\frac{1}{2}+{n_{1}}\left({\frac{2}{{1-{\epsilon_{1}}}}-\frac{1}{2}-\frac{{m\epsilon_{1}^{m}}}{{1-\epsilon_{1}^{m}}}}\right).
Case 3.

(Average AoI for Classical ARQ) We also find that the average AoI expression for TARQ is a specific case of our unified result in (14). For classical ARQ, the transmitter re-transmits the same packet till the packet is successfully received, while the maximum retransmissions is not limited. The classical ARQ is a special case of TARQ where mm\to\infty. Then, by calculating the limit limmΔTARQ\mathop{\lim}\limits_{m\to\infty}{\Delta_{\rm TARQ}}, we can obtain the average AoI as the Proposition 2 in [22]:

Δ¯ClassicalARQ=12+n1(21ϵ112).{\bar{\Delta}_{\rm Classical-ARQ}}=-\frac{1}{2}+{n_{1}}\left({\frac{2}{{1-{\epsilon_{1}}}}-\frac{1}{2}}\right).
Case 4.

(Average AoI for HARQ-IR) By removing the effect of delay elements such that τc=τp=τd=τf=0\tau_{\rm c}=\tau_{\rm p}=\tau_{\rm d}=\tau_{\rm f}=0, we find that the result in (14) is transformed to a variant of that in [21], given as:

Δ¯HARQIR=12+n1+i=1m1(ni+1ni)ϵi1ϵm+n12+i=1m1(ni+12ni2)ϵi2(n1+i=1m1(ni+1ni)ϵi).\bar{\Delta}_{\rm HARQ-IR}=-\frac{1}{2}+\frac{{{n_{1}}+\sum\limits_{i=1}^{m-1}{\left({{n_{i+1}}-{n_{i}}}\right){\epsilon_{i}}}}}{{1-{\epsilon_{m}}}}+\frac{{n_{1}^{2}+\sum\limits_{i=1}^{m-1}{\left({n_{i+1}^{2}-n_{i}^{2}}\right){\epsilon_{i}}}}}{{2\left({{n_{1}}+\sum\limits_{i=1}^{m-1}{\left({{n_{i+1}}-{n_{i}}}\right){\epsilon_{i}}}}\right)}}.

III-A3 Average Peak AoI

Definition 3 has indicated that Δ¯PlimN1Nj=1NΔ(tjS1)\bar{\Delta}^{\rm P}\triangleq\lim\limits_{N\to\infty}\frac{1}{N}\sum_{j=1}^{N}\Delta\left(t^{S}_{j}-1\right). Note that for reactive scheme, the terms in the summation can be obtained from (12), given as Δ(tjS1)=TjReac1τf+τVj1Reac\Delta\left(t^{S}_{j}-1\right)=T^{\rm Reac}_{j}-1-\tau_{\rm f}+\tau^{\rm Reac}_{V_{j-1}}. Then, from the law of large numbers, we can obtain the following almost sure equality

Δ¯ReactiveP=1τf+𝔼TjReac+𝔼τVj1Reac.\bar{\Delta}^{\rm P}_{\rm Reactive}=-1-\tau_{\rm f}+\mathbb{E}T^{\rm Reac}_{j}+\mathbb{E}\tau^{\rm Reac}_{V_{j-1}}.

Then, by applying the available first moments of TjReacT^{\rm Reac}_{j} and τVjReac\tau^{\rm Reac}_{V_{j}} given in (6) and (10), we obtain the average Peak AoI for reactive HARQ in the following proposition.

Proposition 2.

(The Generalized Closed-form Average Peak AoI Expression for Reactive HARQ) For reactive HARQ with maximum retransmissions mm, block assignment vector 𝐧=(n1,n2,,nm){\bf n}=\left(n_{1},n_{2},\cdots,n_{m}\right) and error probability vector 𝐞=(ϵ1,ϵ2,,ϵm){\bf e}=\left(\epsilon_{1},\epsilon_{2},\cdots,\epsilon_{m}\right), the average Peak AoI can be calculated by

Δ¯ReactiveP=1τf(τc+nm+m𝒯)ϵm1ϵm+21ϵm(τc+n1+𝒯+i=1m1(ni+1ni+𝒯)ϵi).{\bar{\Delta}^{\rm P}_{\rm Reactive}}=-1-{\tau_{\rm f}}-\frac{{\left({{\tau_{\rm c}}+{n_{m}}+m\mathcal{T}}\right){\epsilon_{m}}}}{{1-{\epsilon_{m}}}}+\frac{2}{{1-{\epsilon_{m}}}}\left({{\tau_{\rm c}}+{n_{1}}+\mathcal{T}+\sum\limits_{i=1}^{m-1}{\left({{n_{i+1}}-{n_{i}}+\mathcal{T}}\right){\epsilon_{i}}}}\right).

III-B Proactive Scheme

III-B1 Average AoI

Let τiProac\tau^{\rm Proac}_{i} represent the elapsed time form generating an update to receiving its ithi^{\rm th} feedback signal for proactive HARQ, we can observe from Fig. 2(b) that τiProac=ni+τc+𝒯\tau^{\rm Proac}_{i}=n_{i}+\tau_{\rm c}+\mathcal{T}. As such, the first and second moments for τVjProac\tau^{\rm Proac}_{V_{j}} is correspondingly given as follows.

Lemma 6.

The first and second moments for τVjProac\tau^{\rm Proac}_{V_{j}} are

𝔼τVjProac\displaystyle\mathbb{E}\tau^{\rm Proac}_{V_{j}} =n1+τc+𝒯+i[m1](ni+1ni)ϵi1ϵm,\displaystyle=n_{1}+\tau_{\rm c}+\mathcal{T}+\sum_{i\in[m-1]}\frac{\left(n_{i+1}-n_{i}\right)\epsilon_{i}}{1-\epsilon_{m}}, (15)
𝔼(τVjProac)2\displaystyle\mathbb{E}\left(\tau^{\rm Proac}_{V_{j}}\right)^{2} =(n1+τc+𝒯)2+i[m1](ni+1ni)(ni+1+ni+2τc+2𝒯)ϵi1ϵm.\displaystyle=\left(n_{1}+\tau_{\rm c}+\mathcal{T}\right)^{2}+\sum_{i\in[m-1]}\frac{\left(n_{i+1}-n_{i}\right)\left(n_{i+1}+n_{i}+2\tau_{\rm c}+2\mathcal{T}\right)\epsilon_{i}}{1-\epsilon_{m}}.
Proof.

By utilizing the distributions in (4) and the definitions that 𝔼τVjProac=i[m]τiProac(Vj=i)\mathbb{E}\tau^{\rm Proac}_{V_{j}}=\sum_{i\in[m]}\tau^{\rm Proac}_{i}\cdot\mathbb{P}\left({V_{j}}=i\right) and 𝔼(τVjProac)2=i[m](τiProac)2(Vj=i)\mathbb{E}\left(\tau^{\rm Proac}_{V_{j}}\right)^{2}=\sum_{i\in[m]}\left(\tau^{\rm Proac}_{i}\right)^{2}\cdot\mathbb{P}\left({V_{j}}=i\right), we obtain the results in (15). ∎

Denote the interval between the (j1)th{\left(j-1\right)}^{\rm th} and the jthj^{\rm th} successful decoding for proactive HARQ as TjProacT^{\rm Proac}_{j}, we can similarly derive that

TjProac=τmProacRj+τVjProac+τf.T^{\rm Proac}_{j}=\tau^{\rm Proac}_{m}R_{j}+\tau^{\rm Proac}_{V_{j}}+\tau_{\rm f}. (16)

Then we have the first and second moments of TjProacT^{\rm Proac}_{j} as follows.

Lemma 7.

The first and second moments for TjProacT^{\rm Proac}_{j} are

𝔼TjProac=n1+τc+𝒯+i[m1](ni+1ni)ϵi1ϵm,\displaystyle\mathbb{E}T^{\rm Proac}_{j}=\frac{n_{1}+\tau_{\rm c}+\mathcal{T}+\sum_{i\in[m-1]}{\left(n_{i+1}-n_{i}\right)\epsilon_{i}}}{1-\epsilon_{m}}, (17)
𝔼(TjProac)2=(nm+τc+𝒯)2(1+ϵm)(1ϵm)2\displaystyle\mathbb{E}\left(T^{\rm Proac}_{j}\right)^{2}=\frac{\left(n_{m}+\tau_{\rm c}+\mathcal{T}\right)^{2}\left(1+\epsilon_{m}\right)}{\left(1-\epsilon_{m}\right)^{2}}
i[m1]1ϵi1ϵm(ni+1ni)[ni+ni+1+2τc+2𝒯+2(nm+τc+𝒯)ϵm1ϵm].\displaystyle-\sum_{i\in[m-1]}\frac{1-\epsilon_{i}}{1-\epsilon_{m}}\left(n_{i+1}-n_{i}\right)\left[n_{i}+n_{i+1}+2\tau_{\rm c}+2\mathcal{T}+\frac{2\left(n_{m}+\tau_{\rm c}+\mathcal{T}\right)\epsilon_{m}}{1-\epsilon_{m}}\right].
Proof.

From (16), we obtain that

𝔼TjProac=τmProac𝔼Rj+𝔼τVjProac+τf,\displaystyle\mathbb{E}T^{\rm Proac}_{j}=\tau^{\rm Proac}_{m}\mathbb{E}R_{j}+\mathbb{E}\tau^{\rm Proac}_{V_{j}}+\tau_{\rm f}, (18)
𝔼(TjProac)2=𝔼(τmProacRj+τVjProac+τf)2.\displaystyle\mathbb{E}\left(T^{\rm Proac}_{j}\right)^{2}=\mathbb{E}\left(\tau^{\rm Proac}_{m}R_{j}+\tau^{\rm Proac}_{V_{j}}+\tau_{\rm f}\right)^{2}.

Substitute (4) and (5) into (18), we obtain the results in (17). ∎

For the considered proactive HARQ model, the average AoI can be similarly expressed as

Δ¯Proactive=𝔼tjΔ(t)𝔼TjProac=𝔼TjProac2𝔼(TjProac)2+𝔼τVjProacτf12.\bar{\Delta}_{\rm Proactive}=\frac{\mathbb{E}\sum_{t\in\mathcal{I}_{j}}\Delta{\left(t\right)}}{\mathbb{E}T^{\rm Proac}_{j}}=\frac{\mathbb{E}T^{\rm Proac}_{j}}{2\mathbb{E}\left(T^{\rm Proac}_{j}\right)^{2}}+\mathbb{E}\tau^{\rm Proac}_{V_{j}}-\tau_{\rm f}-\frac{1}{2}. (19)

Hence, applying the first and second moments given in Lemma 6 and Lemma 7 leads to the explicit expression in the following proposition.

Proposition 3.

(The Generalized Closed-form Average AoI Expression for Proactive HARQ) For proactive HARQ with maximum retransmissions mm, block assignment vector 𝐧=(n1,n2,,nm){\bf n}=\left(n_{1},n_{2},\cdots,n_{m}\right) and error probability vector 𝐞=(ϵ1,ϵ2,,ϵm){\bf e}=\left(\epsilon_{1},\epsilon_{2},\cdots,\epsilon_{m}\right), the average AoI can be calculated by

Δ¯Proactive=12τf+τc+n1+𝒯+i=1m1(ni+1ni)ϵi1ϵm+\displaystyle{\bar{\Delta}_{\rm Proactive}}=-\frac{1}{2}-{\tau_{\rm f}}+\frac{{{\tau_{\rm c}}+{n_{1}}+\mathcal{T}+\sum\limits_{i=1}^{m-1}{\left({{n_{i+1}}-{n_{i}}}\right){\epsilon_{i}}}}}{{1-{\epsilon_{m}}}}+ (20)
(τc+n1+𝒯)2+i=1m1(ni+1ni)(2τc+2𝒯+ni+1+ni)ϵi2(τc+n1+𝒯+i=1m1(ni+1ni)ϵi).\displaystyle\frac{{{{\left({{\tau_{\rm c}}+{n_{1}}+\mathcal{T}}\right)}^{2}}+\sum\limits_{i=1}^{m-1}{\left({{n_{i+1}}-{n_{i}}}\right)\left({2{\tau_{\rm c}}+2\mathcal{T}+{n_{i+1}}+{n_{i}}}\right){\epsilon_{i}}}}}{{2\left({{\tau_{\rm c}}+{n_{1}}+\mathcal{T}+\sum\limits_{i=1}^{m-1}{\left({{n_{i+1}}-{n_{i}}}\right){\epsilon_{i}}}}\right)}}.

III-B2 Average Peak AoI

As the definition indicates, we have Δ¯PlimN1Nj=1NΔ(tjS1)\bar{\Delta}^{\rm P}\triangleq\lim\limits_{N\to\infty}\frac{1}{N}\sum_{j=1}^{N}\Delta\left(t^{S}_{j}-1\right). For proactive HARQ, Δ(tjS1)\Delta\left(t^{S}_{j}-1\right) is given by Δ(tjS1)=TjProac1τf+τVj1Proac\Delta\left(t^{S}_{j}-1\right)=T^{\rm Proac}_{j}-1-\tau_{\rm f}+\tau^{\rm Proac}_{V_{j-1}}. Then, from the law of large numbers, we get the following equality as

Δ¯ProactiveP=1τf+𝔼TjProac+𝔼τVj1Proac.\bar{\Delta}^{\rm P}_{\rm Proactive}=-1-\tau_{\rm f}+\mathbb{E}T^{\rm Proac}_{j}+\mathbb{E}\tau^{\rm Proac}_{V_{j-1}}.

Finally, by applying the available first of TjProacT^{\rm Proac}_{j} and τVjProac\tau^{\rm Proac}_{V_{j}} given in (15) and (18), we obtain the average Peak AoI for reactive HARQ in the following proposition.

Proposition 4.

(The Generalized Closed-form Average Peak AoI Expression for Proactive HARQ) For proactive HARQ with maximum retransmissions mm, block assignment vector 𝐧=(n1,n2,,nm){\bf n}=\left(n_{1},n_{2},\cdots,n_{m}\right) and error probability vector 𝐞=(ϵ1,ϵ2,,ϵm){\bf e}=\left(\epsilon_{1},\epsilon_{2},\cdots,\epsilon_{m}\right), the average Peak AoI can be calculated by

Δ¯ProactiveP=1τf(τc+nm+𝒯)ϵm1ϵm+21ϵm(τc+n1+𝒯+i=1m1(ni+1ni)ϵi).{\bar{\Delta}^{\rm P}_{\rm Proactive}}=-1-{\tau_{\rm f}}-\frac{{\left({{\tau_{\rm c}}+{n_{m}}+\mathcal{T}}\right){\epsilon_{m}}}}{{1-{\epsilon_{m}}}}+\frac{2}{{1-{\epsilon_{m}}}}\left({{\tau_{\rm c}}+{n_{1}}+\mathcal{T}+\sum\limits_{i=1}^{m-1}{\left({{n_{i+1}}-{n_{i}}}\right){\epsilon_{i}}}}\right).

III-B3 Case Study: Rateless Codes

For rateless codes, the encoder can generate as many symbols as possible to achieve error-free transmission. As such, rateless codes can be regarded as a type of proactive HARQ with infinite code-length setup. By leveraging the obtained results regarding to proactive HARQ, the average AoI and average Peak AoI of rateless code are give in the following Propositions.

Proposition 5.

(The Generalized Closed-form Average AoI Expression for Rateless Codes) For rateless codes transmitted over a noisy channel with non-trivial coding delay τc\tau_{\rm c}, propagation delay τp\tau_{\rm p}, decoding delay τd\tau_{\rm d} and feedback delay τf\tau_{\rm f}, the average AoI can be calculated by

Δ¯Rateless=12+𝒯+n1++i=1(ni+1ni)ϵi\displaystyle{\bar{\Delta}_{\rm Rateless}}=-\frac{1}{2}+\mathcal{T}+{n_{1}}++\sum\limits_{i=1}^{\infty}{\left({{n_{i+1}}-{n_{i}}}\right){\epsilon_{i}}} (21)
+(τc+n1+𝒯)2+i=1(ni+1ni)(2τc+2𝒯+ni+1+ni)ϵi2(τc+n1+𝒯+i=1(ni+1ni)ϵi).\displaystyle+\frac{{{{\left({{\tau_{\rm c}}+{n_{1}}+\mathcal{T}}\right)}^{2}}+\sum\limits_{i=1}^{\infty}{\left({{n_{i+1}}-{n_{i}}}\right)\left({2{\tau_{\rm c}}+2\mathcal{T}+{n_{i+1}}+{n_{i}}}\right){\epsilon_{i}}}}}{{2\left({{\tau_{\rm c}}+{n_{1}}+\mathcal{T}+\sum\limits_{i=1}^{\infty}{\left({{n_{i+1}}-{n_{i}}}\right){\epsilon_{i}}}}\right)}}.
Proof.

Rateless codes is a special type of proactive HARQ where mm\to\infty. Since limmϵm=0\lim\limits_{m\to\infty}\epsilon_{m}=0, we can obtain the average AoI of rateless codes as in (21) by calculating the limit limmΔ¯Proactive\mathop{\lim}\limits_{m\to\infty}{\bar{\Delta}_{\rm Proactive}}. ∎

Proposition 6.

(The Generalized Closed-form Average Peak AoI Expression for Rateless Codes) For rateless codes transmitted over the channel with non-trivial coding delay τc\tau_{\rm c}, propagation delay τp\tau_{\rm p}, decoding delay τd\tau_{\rm d} and feedback delay τf\tau_{\rm f}, the average AoI can be calculated by

Δ¯RatelessP=1τf+2(τc+n1+𝒯+i=1(ni+1ni)ϵi).{\bar{\Delta}^{\rm P}_{\rm Rateless}}=-1-{\tau_{\rm f}}+{2}\left({{\tau_{\rm c}}+{n_{1}}+\mathcal{T}+\sum\limits_{i=1}^{\infty}{\left({{n_{i+1}}-{n_{i}}}\right){\epsilon_{i}}}}\right). (22)
Proof.

Rateless codes is a special case of proactive HARQ where mm\to\infty. By calculating the limit limmΔ¯Proactive\mathop{\lim}\limits_{m\to\infty}{\bar{\Delta}_{\rm Proactive}}, we can obtain the average AoI of rateless codes as in (22). ∎

Remark 2.

Note that there are infinite series in (21) and (22), which are

i=1(ni+1ni)ϵi,\displaystyle\sum\limits_{i=1}^{\infty}{\left({{n_{i+1}}-{n_{i}}}\right){\epsilon_{i}}}, (23)
i=1(ni+1ni)(2τc+2𝒯+ni+1+ni)ϵi.\displaystyle\sum\limits_{i=1}^{\infty}{\left({{n_{i+1}}-{n_{i}}}\right)\left({2{\tau_{\rm c}}+2\mathcal{T}+{n_{i+1}}+{n_{i}}}\right){\epsilon_{i}}}.

In this regard, the sufficient condition that Proposition 5 and Proposition 6 exist is that the infinite series in (23) converge to some finite values. In the following Lemma 8, we would like to discuss this issue.

Lemma 8.

The infinite series in (23) are always bounded (less than or equal to some finite number).

Proof.

Please refer to Appendix C. ∎

III-C Reactive HARQ vs. Proactive HARQ

Corollary 1.

(Reactive HARQ vs. Proactive HARQ) The average age performance of reactive HARQ would not exceed that of proactive HARQ under the same block assignment vector 𝐧\bf n and the same error probability vector 𝐞\bf e. The necessary and sufficient condition for their equivalence is m=1m=1 or 𝒯=0\mathcal{T}=0.

Proof.

Please refer to Appendix D. ∎

Corollary 1 demonstrates that Δ¯ReactiveΔ¯Proactive\bar{\Delta}_{\rm Reactive}\geq\bar{\Delta}_{\rm Proactive}, where the equivalence happens only if ii) m=1m=1, in such a case, both the reactive HARQ and the proactive HARQ turns to a open-loop fashion non-ARQ system, and the system does not send any incremental redundancy; iiii) 𝒯=0\mathcal{T}=0, this condition infers to an ideal assumption where the propagation delay, decoding delay and feedback delay are negligible. In this regard, the RTT issue of reactive HARQ does not exist any more, and thus the considered reactive scheme is the same as the proactive one.

III-D Average Age in the FBL Regime

With the above closed-form results, we observe that the average AoI of a HARQ-based system can be directly evaluated by determining the error probability vector e and the block assignment vector n. The error probability ϵi\epsilon_{i} is affected by three factors, which are, ii) the channel condition; iiii) the coding and decoding technique; iiiiii) the message length kk and the code length nin_{i}. As such, the framework given in this paper is general-purpose, enabling potential AoI research under different coding schemes and channel conditions.

For instance, the given generalized expressions allow us to adopt the FBL results in [32] to evaluate the AoI of the considered HARQ protocols. Over the power-limited AWGN channel with SNR γ\gamma, the error probability ϵi\epsilon_{i} can be approximated by the Theorem 54 of [32] as:444Here we focus on FBL analysis as a case study. Note that ϵi\epsilon_{i} can also be characterized by other specific error-correcting techniques.

ϵiQ(C(γ)k/ni12log2ni/niV(γ)/ni).{\epsilon_{i}}\approx Q\left({\frac{C\left({\gamma}\right)-k/{n_{i}}-\frac{1}{2}{\log_{2}{n_{i}}}/{n_{i}}}{{\sqrt{V\left(\gamma\right)/{n_{i}}}}}}\right). (24)

where CC is the channel capacity with C=log2(1+γ)C=\log_{2}\left(1+\gamma\right), VV is the channel dispersion with V=(11(1+γ)2)log22eV=(1-\frac{1}{\left(1+\gamma\right)^{2}})\log_{2}^{2}e, and Q()Q\left(\cdot\right) denotes the QQ function with

Q(x)=x12πexp(t22)𝑑t.Q\left(x\right)=\int_{x}^{\infty}\frac{1}{\sqrt{2\pi}}\exp\left(-\frac{t^{2}}{2}\right)dt.

Substitute (24) into (14) and (20), we obtain the average AoI closed-form expressions in the FBL regime.

IV Age-optimal Block Assignment

In addition to the error probability 𝐞\bf e, which has been determined by the available finite-length results in (24), the other factor that can significantly affect the age performance is the block assignment vector 𝐧\bf n. The block assignment vector 𝐧\bf n is an important system parameter regarding the design of the transmission strategy. To this end, this section provides design guidelines for the system parameters selection for an status update system to improve the average information timeliness. By this means, we would like to answer how many retransmissions should be and what lengths they are in an age-optimal system.

IV-A Problem Formulation

Input: The signal-to-noise ratio γ\gamma; The message length kk; The lower bound of the range of block length nminn_{\min}; The upper bound of the range of block length nmaxn_{\max}; The system delay τc\tau_{\rm c}, τp\tau_{\rm p}, τd\tau_{\rm d} and τf\tau_{\rm f};
Output: The optimal block assignment vector 𝐧optimal{\bf n}_{\rm optimal}; The minimum average age Δ¯min\bar{\Delta}_{\min};
1 Initialization: Δ¯min=\bar{\Delta}_{\min}=\infty; 𝒮={0,1}nmaxnmin+1\mathcal{S}=\left\{0,1\right\}^{n_{\max}-n_{\min}+1};
2 for 𝐩{\bf p} in 𝒮\mathcal{S}  do
3       Map vector 𝐩\bf p to the block assignment vector 𝐧\bf n;
4       According to the obtained 𝐧\bf n, calculate the average age Δ¯\bar{\Delta} by using (14) or (20) ;
5       if Δ¯<Δ¯min\bar{\Delta}<\bar{\Delta}_{\min} then
6             Update Δ¯min=Δ¯\bar{\Delta}_{\min}=\bar{\Delta};
7             Update 𝐧optimal=𝐧{\bf n}_{\rm optimal}={\bf n};
8            
9      
return 𝐧optimal{\bf n}_{\rm optimal} and Δ¯min\bar{\Delta}_{\min}
Algorithm 1 The algorithm for solving Problem 1.

We establish an average AoI minimization problem here to further explore the age-optimal transmission mechanism in the FBL regime with non-trivial delay:

1) Objective function: To minimize the average age Δ¯\bar{\Delta}.

2) Decision variable: The block assignment vector 𝐧=(n1,n2,,nm){\bf n}=\left(n_{1},n_{2},\cdots,n_{m}\right).

Problem 1.

Age-optimal block assignment for reactive HARQ (or proactive HARQ)

min𝐧\displaystyle\min_{\bf n} Δ¯ReactiveorΔ¯Proactive\displaystyle{\bar{\Delta}_{\rm Reactive}}\quad{\rm or}\quad{\bar{\Delta}_{\rm Proactive}}
s.t. c1:nminn1<n2<<nmnmax,\displaystyle c_{1}:{{n_{\min}}\leq{n_{1}}<{n_{2}}<\ldots<{n_{m}}\leq{n_{\max}}},
c2:1mnmaxnmin+1,\displaystyle c_{2}:{1\leq m\leq n_{\max}-n_{\min}+1},
c3:ϵi=Q(C(γ)k/ni12log2ni/niV(γ)/ni),\displaystyle c_{3}:{\epsilon_{i}}=Q\left({\frac{C\left({\gamma}\right)-k/{n_{i}}-\frac{1}{2}{\log_{2}n_{i}}/{n_{i}}}{{\sqrt{V\left(\gamma\right)/{n_{i}}}}}}\right),
c4:m,ni+,i=1,,m.\displaystyle c_{4}:{m,{n_{i}}\in{\mathbb{Z}^{+}},i=1,\ldots,m}.

Note that the decision variable 𝐧\bf n is a variable-length vector with infinite solution space, we introduce nminn_{\min} and nmaxn_{\max} as constraints of the solution space, which denotes the lower bound and the upper bound of the range of block length, respectively.

IV-B Solutions and Discussions

Problem 1 is a nonlinear integer problem. To solve the optimal solution of Problem 1, an auxiliary vector 𝐩𝒮{0,1}nmaxnmin+1{\bf p}\in\mathcal{S}\triangleq\left\{0,1\right\}^{n_{\max}-n_{\min}+1} can be introduced here.

Lemma 9.

There exists an one-to-one mapping between vectors 𝐧\bf n and 𝐩\bf p.

Proof.

Please refer to appendix E, where we construct a specific one-to-one mapping function between. ∎

Lemma 9 illustrates that the introduced auxiliary vector 𝐩\bf p can be regarded as an index of the solution space of Problem 1, which can help us traverse the entire solution space efficiently and find the optimal solution. The detailed algorithm process is provided in Algorithm 1.

Refer to caption
Figure 3: Age-optimal block assignment of reactive HARQ and proactive HARQ. Here the length of message is k=100k=100, the minimum code length nmin=100n_{\min}=100, the maximum code length nmax=120n_{\max}=120, the coding delay is τc=20\tau_{\rm c}=20, the decoding delay is τd=30\tau_{\rm d}=30.

Fig. 3 gives some detailed examples of the solved optimal block assignment vector 𝐧optimal{\bf n}_{\rm optimal} under different protocols, SNRs, and propagation delays. For example, under SNR=0.7=0.7 dB and τp=0\tau_{\rm p}=0, the optimal block assignment vector for reactive HARQ is 𝐧optimal=(105,120){\bf n}_{\rm optimal}=\left(105,120\right); under SNR=1.9=1.9 dB and τp=20\tau_{\rm p}=20, the optimal block assignment vector for reactive HARQ is 𝐧optimal=(100,112,120){\bf n}_{\rm optimal}=\left(100,112,120\right).

Fig. 3. leads to the following conclusions:

  • For proactive HARQ, the finest grained symbol-by-symbol strategy always minimizes the average AoI.

  • For reactive HARQ, the age-optimal block assignment varies among different SNRs and propagation delays. As the propagation delay increases, the number of retransmissions will monotonically decrease and finally converge to m=1m=1. In such a case, the transmission scheme turns an open-loop fashion without any retransmission. This indicates that there exists a threshold of the propagation delay, only within which retransmission is beneficial to AoI.

  • From a perspective of channel coding, we can see that the trade-off between reliability and effectiveness can be well evaluated by the new metric, AoI. It is well known that a longer code length can improve the reliability while sacrificing the effectiveness; however, what is not fully explored is that an appropriate choice of code length can minimize the AoI.

IV-C A Heuristic Algorithm for Reactive HARQ

For proactive HARQ, it has been empirically shown that the finest grained strategy minimizes the average AoI; however, for reactive HARQ, the age-optimal strategies vary along with channel conditions and propagation delay. Thus, to repeatedly determine the age-optimal scheme requires amounts of calculations. It is also pertinent to note that the implementation of Algorithm 1 aims to exhaustively search the whole solution space to find an age-optimal block assignment strategy. As such, the complexity of such an Algorithm is exponentially increasing with the range of code length nmaxnmin+1n_{\max}-n_{\min}+1. For a broader range of code length, here we heuristically provide a sub-optimal algorithm to circumvent the high-complexity issue. Specifically, the heuristic sub-optimal algorithm is given in Algorithm 2.

Input: The signal-to-noise ratio (SNR); The message length kk; The lower bound of the range of block length nminn_{\min}; The upper bound of the range of block length nmaxn_{\max}; The system delay τc\tau_{\rm c}, τp\tau_{\rm p}, τd\tau_{\rm d} and τf\tau_{\rm f};
Output: The optimal block assignment vector 𝐧optimal{\bf n}_{\rm optimal}; The minimum average age Δ¯min\bar{\Delta}_{\min};
1 Initialization: Δ¯min=\bar{\Delta}_{\min}=\infty;
2 for mm\leftarrow 11 to nmaxnmin+1n_{\max}-n_{\min}+1  do
3       Δ¯min,m=\bar{\Delta}_{\min,m}=\infty;
4       Construct the sub set 𝒮m\mathcal{S}_{m};
5       for 𝐩\bf p in 𝒮m\mathcal{S}_{m} do
6             Map vector 𝐩\bf p to the block assignment vector 𝐧\bf n;
7             According to the obtained 𝐧\bf n, calculate the average age Δ¯\bar{\Delta} by using (14) or (20) ;
8             if Δ¯<Δ¯min,m\bar{\Delta}<\bar{\Delta}_{\min,m} then
9                   Update Δ¯min,m=Δ¯\bar{\Delta}_{\min,m}=\bar{\Delta};
10                   Update 𝐧optimal,m=𝐧{\bf n}_{{\rm optimal},m}={\bf n};
11                  
12            
13      if Δ¯min,m>Δ¯min,m1\bar{\Delta}_{\min,m}>\bar{\Delta}_{\min,m-1} then
14             Update Δ¯min=Δ¯min,m1\bar{\Delta}_{\min}=\bar{\Delta}_{\min,m-1};
15             Update 𝐧optimal=𝐧optimal,m1{\bf n}_{{\rm optimal}}={\bf n}_{{\rm optimal},m-1};
16             break;
17            
18      else
19            Update Δ¯min=Δ¯min,m\bar{\Delta}_{\min}=\bar{\Delta}_{\min,m};
20             Update 𝐧optimal=𝐧optimal,m{\bf n}_{{\rm optimal}}={\bf n}_{{\rm optimal},m};
21      
return 𝐧optimal{\bf n}_{\rm optimal} and Δ¯min\bar{\Delta}_{\min}
Algorithm 2 The sub-optimal algorithm for solving Problem 1

The design of such an algorithm is based on the empirical observation that the age-optimal mm is always monotonically decreasing with propagation delay (see Fig. 3). Meanwhile, the age-optimal mm for reactive HARQ tends to remain small since large mm will lead to multiple RTT and thus result in staleness of information. Attributed to the above factors, heuristically, we denote the minimal age under a fixed mm as Δ¯min,m\bar{\Delta}_{\min,m} and recursively search the solution space 𝒥m{(n1,,ni):i=m,nmin<n1<n2<ni<nmax}\mathcal{J}_{m}\triangleq\left\{\left(n_{1},\cdots,n_{i}\right):i=m,n_{\min}<n_{1}<n_{2}\cdots<n_{i}<n_{\max}\right\} with increasing value of mm. The on-the-fly searching process will terminate only if Δ¯min,m>Δ¯min,m1\bar{\Delta}_{\min,m}>\bar{\Delta}_{\min,m-1}. By this means, the algorithm outputs Δ¯min,m1\bar{\Delta}_{\min,m-1} as the optimal age. In such a case, this algorithm eliminates the need for searching the sub-space i=m+1nmaxnmin+1𝒮i𝒮\cup_{i=m+1}^{n_{\max}-n_{\min}+1}\mathcal{S}_{i}\subseteq\mathcal{S}, thereby bypassing the calculations required for searching for the whole solution space.

Note that an auxiliary set 𝒮m{𝐩{0,1}nmaxnmin+1:𝐩1=m}\mathcal{S}_{m}\triangleq\left\{{\bf p}\in\left\{0,1\right\}^{n_{\max}-n_{\min}+1}:\left\|{\bf p}\right\|_{1}=m\right\} is also introduced in Algorithm 2 to assist high-efficiency searching, where 1\left\|\cdot\right\|_{1} represents the implementation of 1\ell_{1}-norm.

Lemma 10.

There exists an one-to-one mapping between 𝐧𝒥m{\bf n}\in\mathcal{J}_{m} and 𝐩𝒮m{\bf p}\in\mathcal{S}_{m}.

Proof.

The one-to-one mapping function between 𝐧𝒥m{\bf n}\in\mathcal{J}_{m} and 𝐩𝒮m{\bf p}\in\mathcal{S}_{m} is the same as that of Lemma 9, which has been discussed in Appendix E in detail. ∎

V Numerical Results

V-A The Closed-Form Results

Refer to caption
Figure 4: Instantaneous age evolution and the statistic characterizations among reactive HARQ, proactive HARQ and rateless codes. Here the encoding delay is τc=2\tau_{\rm c}=2, the decoding delay is τd=3\tau_{\rm d}=3, the propagation delay is τp=5\tau_{\rm p}=5 and the feedback delay is τf=6\tau_{\rm f}=6.

In addition to the case studies given in Section III. B, we also carry out Monte Carlo simulations to verify our closed-form expressions. For the simulation setup, we leverage an i.i.d uniformly distributed random sequence 𝒳j𝒰(0,1)\mathcal{X}_{j}\sim\mathcal{U}\left(0,1\right) to generate the feedback signal sequence ςj,i,i[m]\varsigma_{j,i},i\in[m] when transmitting message MjM_{j}. Specifically, the feedback signal sequence is generated by

ςj,i=sign(𝒳jϵi), for i[m],\varsigma_{j,i}=sign\left(\mathcal{X}_{j}-\epsilon_{i}\right)\text{, for }i\in[m],

where ϵi\epsilon_{i} is obtained by (24) and sign()sign\left(\cdot\right) is defined as

sign(x)={1 , for x00 , for x<0.sign\left(x\right)=\left\{{\begin{array}[]{*{20}{c}}1\text{ , for }x\geq 0\\ 0\text{ , for }x<0\end{array}}\right..
Refer to caption
Figure 5: Reactive HARQ vs proactive HARQ. Here the message length is k=100k=100, the encoding delay is τc=20\tau_{\rm c}=20, the decoding delay is τc=30\tau_{\rm c}=30, the propagation delay range is τp=0:200\tau_{\rm p}=0:200 and the feedback delay is calculated by τf=τp+1\tau_{\rm f}=\tau_{\rm p}+1. For fairness, here we consider the finest grained block assignment vectors for both reactive HARQ and proactive HARQ, with n1=100,ni=ni1+1n_{1}=100,n_{i}=n_{i-1}+1.

Then, with the feedback signal sequence ςj,i\varsigma_{j,i} in hand, the transmission-decoding model is almost sure, and we can recursively obtain the instantaneous age evolution as shown in Fig. 4. For reactive HARQ and proactive HARQ, we set k=100k=100, m=11m=11 and 𝐧=100:110\mathbf{n}=100:110. For rateless codes, we find that a sufficiently large value mm will directly lead to an almost convergent AoI. Thus, we set m=10000m=10000 for the simulation setup of rateless codes.

Fig. 4(a) demonstrates the instantaneous age evolution for reactive HARQ, proactive HARQ, and rateless codes, respectively. Intuitively, we can observe that the age of reactive HARQ tends to exhibit a number of large sawtooth waveforms, while that of proactive HARQ and rateless codes cut off the large sawtooth waveforms and keep at a relatively low level.

Refer to caption
Figure 6: Average age comparisons among the finest grained reactive HARQ, the optimal reactive HARQ and the finest grained proactive HARQ. Here the message length is k=100k=100, the encoding delay is τc=20\tau_{\rm c}=20, the decoding delay is τc=200\tau_{\rm c}=200, the propagation delay is τp=50\tau_{\rm p}=50 and the feedback delay is τf=51\tau_{\rm f}=51.

Fig. 4(b) and Fig. 4(c) depict the average AoI and average peak AoI comparisons between the simulation results and the analytical closed-form results, wherein the discrete orange points are obtained through Monte Carlo simulations, while the blue curves are plotted by utilizing the available closed-form results given in Section III. It can be seen that the simulation results fit well with the analytical results, verifying that our provided closed-form expressions enable exact and efficient AoI evaluations.

V-B Reactive HARQ vs. Proactive HARQ

Fig. 5 demonstrates the average AoI comparison between reactive HARQ and proactive HARQ from a multi-dimensional perspective. The comparisons are conducted among different settings of mm, τp\tau_{\rm p} and SNR. It is shown that the proactive HARQ surface remains below the reactive HARQ surface. Also, they intersects with each other at m=1m=1. These numerical results are consistent with Corollary 1. In addition, Fig. 5 also illustrates the impact that τp\tau_{\rm p} and mm exert on average AoI. On the one hand, the average AoI is monotonically increasing with respect to the propagation delay τp\tau_{\rm p}. On the other hand, the impact of mm on average AoI could be complex: ii) for proactive HARQ, retransmitting redundancy remains beneficial for AoI performance metric; iiii) for reactive HARQ, however, retransmitting redundancy naturally brings about RTT and thus results in staleness of information when the SNR is high enough to achieve reliable communication; in contrast, if the channel condition is poor, retransmitting redundancy is essential for reliable delivery, and thus may even compensate the AoI losses due to RTT.

V-C The Age-Optimal Block Assignment

Fig. 6 shows a comparison among the finest grained reactive HARQ, the optimal reactive HARQ, and the finest grained proactive HARQ. It is shown that the optimized reactive HARQ approaches proactive HARQ in average AoI performance. Notice that this gain lies in an adaptive block assignment strategy which requires accurate channel status information. In this regard, we find that adopting proactive HARQ for freshness-critical status update systems would be a robust and timeliness-efficient approach.

VI Conclusion and Future Work

In this paper, we have comprehensively considered different types of nontrivial system delay and derived unified closed-form average AoI and average Peak AoI expressions for both reactive HARQ and proactive HARQ. The unifying characteristic of our result has been shown by several case studies, wherein some existing PHY-layer AoI expressions in the literature are shown to be only some specific cases of our result. With these closed-form results in hand, we have theoretically proven that under the same communication conditions, the proactive scheme always outperforms the reactive scheme in terms of the average AoI. Also, a block assignment design framework at the PHY layer has been provided to further achieve timely delivery in a status update system. The simulation and analytical results demonstrate that the age-optimal block assignment strategy of reactive schemes is sensitive to both channel conditions and propagation delays, while that of proactive scheme exhibits both strategy robustness and age superiority. In this regard, we witness the potential for the proactive HARQ to be applied in the freshness-critical system.

The research in this paper also leaves some open challenges and issues for future research. First, it will be an interesting work to carry out AoI analyses and comparisons for some specific state-of-the-art channel coding techniques, such as polar codes, LDPC codes, Turbo codes, and rateless Raptor codes, etc. As such, from the AoI perspective, the trade-off among coding complexity, decoding complexity, codelength, the number of retransmitted packets, and the error probability can be explored. Second, since this work is based on an ideal assumption of perfect feedback, the analysis considering lossy feedback can be further conducted. Third, notwithstanding the AoI superiority of proactive HARQ compared to reactive HARQ, proactive HARQ may consume more energy due to the consecutive retransmissions. To this end, to further investigate the trade-off of proactive HARQ between timeliness and energy efficiency would be an interesting topic.

Appendix A Proof of Lemma 1

The event {Rj=a}\left\{R_{j}=a\right\} is equivalent to

{Rj=a}={ςQj,m=1}i[a]{ςi+Qj1,m=0}.\left\{R_{j}=a\right\}=\left\{\varsigma_{Q_{j},m}=1\right\}\bigcap_{i\in[a]}\left\{\varsigma_{i+Q_{j-1},m}=0\right\}.

Note that the AWGN is i.i.d, we have that

(Rj=a)\displaystyle\mathbb{P}\left(R_{j}=a\right) =(ςQj,m=1)i[a](ςi+Qj1,m=0)\displaystyle=\mathbb{P}\left(\varsigma_{Q_{j},m}=1\right)\cdot\prod_{i\in[a]}\mathbb{P}\left(\varsigma_{i+Q_{j-1},m}=0\right)
=(1ϵm)ϵma.\displaystyle=\left(1-\epsilon_{m}\right)\epsilon^{a}_{m}.

Also, the probability of event {Vj=i}\left\{V_{j}=i\right\} can be expressed as

(Vj=i)=({ςQj,m=1,ςQj,i=1}r[i1]{ςQj,r=0})(ςQj,m=1).\mathbb{P}\left(V_{j}=i\right)=\frac{\mathbb{P}\left(\left\{\varsigma_{Q_{j},m}=1,\varsigma_{Q_{j},i}=1\right\}\bigcap_{r\in[i-1]}\left\{\varsigma_{Q_{j},r}=0\right\}\right)}{\mathbb{P}\left(\varsigma_{Q_{j},m}=1\right)}. (25)

As the variable ςQj,i\varsigma_{Q_{j},i} follows the monotonic property such that

{ςQj,i=1}{ςQj,m=1},for im,\displaystyle\left\{\varsigma_{Q_{j},i}=1\right\}\subseteq\left\{\varsigma_{Q_{j},m}=1\right\},\text{\rm for }i\leq m,
{ςQj,i=0}{ςQj,r=0},for ri.\displaystyle\left\{\varsigma_{Q_{j},i}=0\right\}\subseteq\left\{\varsigma_{Q_{j},r}=0\right\},\text{\rm for }r\leq i.

The event in (25) can be simplified as

{ςQj,m=1,ςQj,i=1}r[i1]{ςQj,r=0}\displaystyle\left\{\varsigma_{Q_{j},m}=1,\varsigma_{Q_{j},i}=1\right\}\bigcap_{r\in[i-1]}\left\{\varsigma_{Q_{j},r}=0\right\} (26)
={ςQj,i=1,ςQj,i1=0}\displaystyle=\left\{\varsigma_{Q_{j},i}=1,\varsigma_{Q_{j},i-1}=0\right\}
={ςQj,i=1}/{ςQj,i1=1}.\displaystyle=\left\{\varsigma_{Q_{j},i}=1\right\}/\left\{\varsigma_{Q_{j},i-1}=1\right\}.

Substituting (26) into (25 results in the probability as

(Vj=i)=ϵi1ϵi1ϵm.\mathbb{P}\left(V_{j}=i\right)=\frac{\epsilon_{i-1}-\epsilon_{i}}{1-\epsilon_{m}}.

Appendix B Proof of Lemma 5

Define GjReac=tjΔ¯(t)G^{\rm Reac}_{j}=\sum_{t\in\mathcal{I}_{j}}\bar{\Delta}\left(t\right), we have that

GjReac=t=tj1StjS1(ttj1Sτf+τVj1Reac)\displaystyle G^{\rm Reac}_{j}=\sum_{t=t^{S}_{j-1}}^{t^{S}_{j}-1}\left(t-t^{S}_{j-1}-\tau_{\rm f}+\tau^{\rm Reac}_{V_{j-1}}\right)
=x=ttj1Sx=0TjReac1(xτf+τVj1Reac)\displaystyle\overset{x=t-t^{S}_{j-1}}{=}\sum_{x=0}^{T^{\rm Reac}_{j}-1}\left(x-\tau_{\rm f}+\tau^{\rm Reac}_{V_{j-1}}\right)
=TjReac(TjReac1)2τfTjReac+τVj1ReacTjReac.\displaystyle=\frac{T^{\rm Reac}_{j}\left(T^{\rm Reac}_{j}-1\right)}{2}-\tau_{\rm f}T^{\rm Reac}_{j}+\tau^{\rm Reac}_{V_{j-1}}T^{\rm Reac}_{j}.

Thus, we have the first moment of SjS_{j} as

𝔼GjReac=𝔼(TjReac)22+𝔼TjReac(𝔼τVj1Reacτf12).\mathbb{E}G^{\rm Reac}_{j}=\frac{\mathbb{E}\left(T^{\rm Reac}_{j}\right)^{2}}{2}+\mathbb{E}T^{\rm Reac}_{j}\left(\mathbb{E}\tau^{\rm Reac}_{V_{j-1}}-\tau_{\rm f}-\frac{1}{2}\right). (27)

With (27), we can obtain the average age as

Δ¯Reactive=𝔼GjReac𝔼TjReac=𝔼TjProac2𝔼(TjReac)2+𝔼τVjReacτf12.\bar{\Delta}_{\rm Reactive}=\frac{\mathbb{E}G^{\rm Reac}_{j}}{\mathbb{E}T^{\rm Reac}_{j}}=\frac{\mathbb{E}T^{\rm Proac}_{j}}{2\mathbb{E}\left(T^{\rm Reac}_{j}\right)^{2}}+\mathbb{E}\tau^{\rm Reac}_{V_{j}}-\tau_{\rm f}-\frac{1}{2}.

Appendix C Proof of Lemma 8

Recall that ϵi\epsilon_{i} a monotonically decreasing infinite sequence with ϵ1>ϵ2>ϵ3>\epsilon_{1}>\epsilon_{2}>\epsilon_{3}>\cdots, we can prove Lemma 8 by adopting the Dirichlet’s test. With Dirichlet’s test, this proof is equivalent to proving that the partial sums

i=1N(ni+1ni),\displaystyle\sum\limits_{i=1}^{N}{\left({{n_{i+1}}-{n_{i}}}\right)}, (28)
i=1N(ni+1ni)(2τc+2𝒯+ni+1+ni).\displaystyle\sum\limits_{i=1}^{N}{\left({{n_{i+1}}-{n_{i}}}\right)\left({2{\tau_{\rm c}}+2\mathcal{T}+{n_{i+1}}+{n_{i}}}\right)}.

are bounded. Evidently, we have the solutions to (28) as

i=1N(ni+1ni)=nN+1n1,\displaystyle\sum\limits_{i=1}^{N}{\left({{n_{i+1}}-{n_{i}}}\right)}=n_{N+1}-n_{1}, (29)
i=1N(ni+1ni)(2τc+2𝒯+ni+1+ni)=(nN+1+τc+𝒯)2(n1+τc+𝒯)2.\displaystyle\sum\limits_{i=1}^{N}{\left({{n_{i+1}}-{n_{i}}}\right)\left({2{\tau_{\rm c}}+2\mathcal{T}+{n_{i+1}}+{n_{i}}}\right)}=\left(n_{N+1}+\tau_{\rm c}+\mathcal{T}\right)^{2}-\left(n_{1}+\tau_{\rm c}+\mathcal{T}\right)^{2}.

Thus, we have that the infinite series in (23) are bounded.

Appendix D Proof of Corollary 1

Subtract Δ¯Proactive{\bar{\Delta}_{\rm Proactive}} from Δ¯Reactive{\bar{\Delta}_{\rm Reactive}}, we have

Δ¯ReactiveΔ¯Proactive=𝒯i=1m1ϵi1ϵm+\displaystyle{\bar{\Delta}_{\rm Reactive}}-{\bar{\Delta}_{\rm Proactive}}=\frac{{\mathcal{T}\sum\limits_{i=1}^{m-1}{{\epsilon_{i}}}}}{{1-{\epsilon_{m}}}}+ (30)
(τc+n1+𝒯)2+i=1m1(ni+1ni+𝒯)(2τc+ni+1+ni+(2i+1)𝒯)ϵi2(τc+n1+𝒯+i=1m1(ni+1ni+𝒯)ϵi)\displaystyle\frac{{{{\left({{\tau_{\rm c}}+{n_{1}}+\mathcal{T}}\right)}^{2}}+\sum\limits_{i=1}^{m-1}{\left({{n_{i+1}}-{n_{i}}+\mathcal{T}}\right)\left({2{\tau_{\rm c}}+{n_{i+1}}+{n_{i}}+\left({2i+1}\right)\mathcal{T}}\right){\epsilon_{i}}}}}{{2\left({{\tau_{\rm c}}+{n_{1}}+\mathcal{T}+\sum\limits_{i=1}^{m-1}{\left({{n_{i+1}}-{n_{i}}+\mathcal{T}}\right){\epsilon_{i}}}}\right)}}
(τc+n1+𝒯)2+i=1m1(ni+1ni)(2τc+2𝒯+ni+1+ni)ϵi2(τc+n1+𝒯+i=1m1(ni+1ni)ϵi)\displaystyle-\frac{{{{\left({{\tau_{\rm c}}+{n_{1}}+\mathcal{T}}\right)}^{2}}+\sum\limits_{i=1}^{m-1}{\left({{n_{i+1}}-{n_{i}}}\right)\left({2{\tau_{\rm c}}+2\mathcal{T}+{n_{i+1}}+{n_{i}}}\right){\epsilon_{i}}}}}{{2\left({{\tau_{\rm c}}+{n_{1}}+\mathcal{T}+\sum\limits_{i=1}^{m-1}{\left({{n_{i+1}}-{n_{i}}}\right){\epsilon_{i}}}}\right)}}
(a)𝒯i=1m1ϵi1ϵm+𝒯i=1m1(ni+1ni)(2i1)ϵi2(τc+n1+𝒯+i=1m1(ni+1ni)ϵi)+\displaystyle\overset{(a)}{\geq}\frac{{\mathcal{T}\sum\limits_{i=1}^{m-1}{{\epsilon_{i}}}}}{{1-{\epsilon_{m}}}}+\frac{{\mathcal{T}\sum\limits_{i=1}^{m-1}{\left({{n_{i+1}}-{n_{i}}}\right)\left({2i-1}\right){\epsilon_{i}}}}}{{2\left({{\tau_{\rm c}}+{n_{1}}+\mathcal{T}+\sum\limits_{i=1}^{m-1}{\left({{n_{i+1}}-{n_{i}}}\right){\epsilon_{i}}}}\right)}}+
𝒯i=1m1(2τc+ni+1+ni+(2i+1)𝒯)ϵi2(τc+n1+𝒯+i=1m1(ni+1ni)ϵi)(b)0.\displaystyle\frac{{\mathcal{T}\sum\limits_{i=1}^{m-1}{\left({2{\tau_{\rm c}}+{n_{i+1}}+{n_{i}}+\left({2i+1}\right)\mathcal{T}}\right){\epsilon_{i}}}}}{{2\left({{\tau_{\rm c}}+{n_{1}}+\mathcal{T}+\sum\limits_{i=1}^{m-1}{\left({{n_{i+1}}-{n_{i}}}\right){\epsilon_{i}}}}\right)}}\overset{(b)}{\geq}0.

D-1 Proof of Sufficiency

If 𝒯=0\mathcal{T}=0, the equal signs at both (a)(a) and (b)(b) in (30) are established ; if m=1m=1, all those sums i=1m1xi0\sum\nolimits_{i=1}^{m-1}{x_{i}}\equiv 0, and thus we have Δ¯ReactiveΔ¯Proactive=0{\bar{\Delta}_{\rm Reactive}}-{\bar{\Delta}_{\rm Proactive}}=0. Therefore, we prove the sufficiency that m=1m=1 or 𝒯=0Δ¯Reactive=Δ¯Proactive{\mathcal{T}=0\Rightarrow\bar{\Delta}_{\rm Reactive}}={\bar{\Delta}_{\rm Proactive}}.

D-2 Proof of Necessity

Assume m>1m>1 and 𝒯0\mathcal{T}\neq 0, we can easily obtain from (30) that Δ¯ReactiveΔ¯Proactive>0{\bar{\Delta}_{\rm Reactive}}-{\bar{\Delta}_{\rm Proactive}}>0. This is equivalent to the necessity that Δ¯Reactive=Δ¯Proactivem=1{\bar{\Delta}_{\rm Reactive}}={\bar{\Delta}_{\rm Proactive}}\Rightarrow m=1 or 𝒯=0\mathcal{T}=0.

Appendix E Constructing the Mapping in Lemma 9

A constructing mapping from 𝐩\bf p to 𝐧\bf n is shown below:

Step 1: Find the indexes of all the zero-value positions of vector 𝐩=(p1,p2,,pnmaxnmin+1){\bf p}=\left(p_{1},p_{2},\cdots,p_{n_{\max}-n_{\min}+1}\right), and store them in an empty set 𝒜\mathcal{A}. For example, if pi=0p_{i}=0, then ii is stored into 𝒜\mathcal{A}.

Step 2: Sort the elements in set 𝒜\mathcal{A} in the ascending order and denote the ordered elements as a vector (a1,a2,,a|𝒜|)\left(a_{1},a_{2},\dots,a_{\left|{\mathcal{A}}\right|}\right), with a1<a2<<a|𝒜|a_{1}<a_{2}<\cdots<a_{\left|{\mathcal{A}}\right|}.

Step 3: The vector 𝐧\bf n can be obtained by

ni=ai+nmin1,i=1,2,,|𝒜|.n_{i}=a_{i}+n_{\min}-1,i=1,2,\cdots,\left|\mathcal{A}\right|. (31)

The above process is reversible, and the mapping from 𝐧\bf n to 𝐩\bf p is elaborated below:

Step 1: Construct the vector 𝐚=(a1,a2,,am){\bf a}=\left(a_{1},a_{2},\cdots,a_{m}\right) by

ai=ninmin+1,i=1,2,,m.a_{i}=n_{i}-n_{\min}+1,i=1,2,\cdots,m. (32)

Step 2: Initialize 𝐩=𝟏nmaxnmin+1{\bf p}={\bf 1}^{n_{\max}-n_{\min}+1} and let pai=0,i=1,2,,mp_{a_{i}}=0,i=1,2,\dots,m, the vector 𝐩\bf p is obtained.

References

  • [1] S. Kaul, R. Yates, and M. Gruteser, “Real-time status: How often should one update?” in 2012 Proceedings IEEE INFOCOM, 2012, pp. 2731–2735.
  • [2] P. Papadimitratos, A. D. La Fortelle, K. Evenssen, R. Brignolo, and S. Cosenza, “Vehicular communication systems: Enabling technologies, applications, and future outlook on intelligent transportation,” IEEE Communications Magazine, vol. 47, no. 11, pp. 84–95, 2009.
  • [3] P. Corke, T. Wark, R. Jurdak, W. Hu, P. Valencia, and D. Moore, “Environmental wireless sensor networks,” Proceedings of the IEEE, vol. 98, no. 11, pp. 1903–1917, 2010.
  • [4] A. A. Reddy, S. Banerjee, A. Gopalan, S. Shakkottai, and L. Ying, “On distributed scheduling with heterogeneously delayed network-state information,” Queueing Systems, vol. 72, no. 3, pp. 193–218, 2012.
  • [5] M. Costa, M. Codreanu, and A. Ephremides, “On the Age of Information in Status Update Systems With Packet Management,” IEEE Transactions on Information Theory, vol. 62, no. 4, pp. 1897–1910, 2016.
  • [6] S. K. Kaul, R. D. Yates, and M. Gruteser, “Status updates through queues,” in 2012 46th Annual Conference on Information Sciences and Systems (CISS), 2012, pp. 1–6.
  • [7] 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 Transactions on Information Theory, vol. 63, no. 11, pp. 7492–7508, 2017.
  • [8] R. D. Yates and S. K. Kaul, “The Age of Information: Real-Time Status Updating by Multiple Sources,” IEEE Transactions on Information Theory, vol. 65, no. 3, pp. 1807–1827, 2019.
  • [9] L. Huang and E. Modiano, “Optimizing age-of-information in a multi-class queueing system,” in 2015 IEEE International Symposium on Information Theory (ISIT), 2015, pp. 1681–1685.
  • [10] E. Najm and E. Telatar, “Status updates in a multi-stream M/G/1/1 preemptive queue,” in IEEE INFOCOM 2018 - IEEE Conference on Computer Communications Workshops (INFOCOM WKSHPS), 2018, pp. 124–129.
  • [11] M. Moltafet, M. Leinonen, and M. Codreanu, “On the Age of Information in Multi-Source Queueing Models,” IEEE Transactions on Communications, vol. 68, no. 8, pp. 5003–5017, 2020.
  • [12] C. Kam, S. Kompella, G. D. Nguyen, and A. Ephremides, “Effect of Message Transmission Path Diversity on Status Age,” IEEE Transactions on Information Theory, vol. 62, no. 3, pp. 1360–1374, 2016.
  • [13] A. M. Bedewy, Y. Sun, and N. B. Shroff, “Minimizing the Age of Information Through Queues,” IEEE Transactions on Information Theory, vol. 65, no. 8, pp. 5215–5232, 2019.
  • [14] ——, “The Age of Information in Multihop Networks,” IEEE/ACM Transactions on Networking, vol. 27, no. 3, pp. 1248–1257, 2019.
  • [15] E. Najm, R. Nasser, and E. Telatar, “Content Based Status Updates,” IEEE Transactions on Information Theory, vol. 66, no. 6, pp. 3846–3863, 2020.
  • [16] S. K. Kaul and R. D. Yates, “Age of Information: Updates with Priority,” in 2018 IEEE International Symposium on Information Theory (ISIT), 2018, pp. 2644–2648.
  • [17] R. D. Yates, “Status Updates through Networks of Parallel Servers,” in 2018 IEEE International Symposium on Information Theory (ISIT), 2018, pp. 2281–2285.
  • [18] K. Chen and L. Huang, “Age-of-information in the presence of error,” in 2016 IEEE International Symposium on Information Theory (ISIT), 2016, pp. 2579–2583.
  • [19] R. D. Yates, E. Najm, E. Soljanin, and J. Zhong, “Timely updates over an erasure channel,” in 2017 IEEE International Symposium on Information Theory (ISIT), 2017, pp. 316–320.
  • [20] E. Najm, E. Telatar, and R. Nasser, “Optimal Age over Erasure Channels,” in 2019 IEEE International Symposium on Information Theory (ISIT), 2019, pp. 335–339.
  • [21] S. C. Bobbili, P. Parag, and J.-F. Chamberland, “Real-Time Status Updates With Perfect Feedback Over Erasure Channels,” IEEE Transactions on Communications, vol. 68, no. 9, pp. 5363–5374, 2020.
  • [22] M. Xie, Q. Wang, J. Gong, and X. Ma, “Age and Energy Analysis for LDPC Coded Status Update With and Without ARQ,” IEEE Internet of Things Journal, vol. 7, no. 10, pp. 10 388–10 400, 2020.
  • [23] J. You, S. Wu, Y. Deng, J. Jiao, and Q. Zhang, “An Age Optimized Hybrid ARQ Scheme for Polar Codes via Gaussian Approximation,” IEEE Wireless Communications Letters, vol. 10, no. 10, pp. 2235–2239, 2021.
  • [24] A. Arafa, K. Banawan, K. G. Seddik, and H. V. Poor, “On Timely Channel Coding with Hybrid ARQ,” in 2019 IEEE Global Communications Conference (GLOBECOM), 2019, pp. 1–6.
  • [25] B. Yu, Y. Cai, D. Wu, and Z. Xiang, “Average Age of Information in Short Packet Based Machine Type Communication,” IEEE Transactions on Vehicular Technology, vol. 69, no. 9, pp. 10 306–10 319, 2020.
  • [26] R. Devassy, G. Durisi, G. C. Ferrante, O. Simeone, and E. Uysal-Biyikoglu, “Delay and Peak-Age Violation Probability in Short-Packet Transmissions,” in 2018 IEEE International Symposium on Information Theory (ISIT), 2018, pp. 2471–2475.
  • [27] D. Li, S. Wu, Y. Wang, J. Jiao, and Q. Zhang, “Age-Optimal HARQ Design for Freshness-Critical Satellite-IoT Systems,” IEEE Internet of Things Journal, vol. 7, no. 3, pp. 2066–2076, 2020.
  • [28] Discussion on Explicit HARQ-ACK Feedback for Configured Grant Transmission, document R1-1903079, 3GPP TSG RAN WG1 96, 2019.
  • [29] R. Abreu, T. Jacobsen, G. Berardinelli, K. Pedersen, I. Z. Kovács, and P. Mogensen, “Power control optimization for uplink grant-free URLLC,” in 2018 IEEE Wireless Communications and Networking Conference (WCNC), 2018, pp. 1–6.
  • [30] N. H. Mahmood, R. Abreu, R. Böhnke, M. Schubert, G. Berardinelli, and T. H. Jacobsen, “Uplink Grant-Free Access Solutions for URLLC services in 5G New Radio,” in 2019 16th International Symposium on Wireless Communication Systems (ISWCS), 2019, pp. 607–612.
  • [31] Y. Liu, Y. Deng, M. Elkashlan, A. Nallanathan, and G. K. Karagiannidis, “Analyzing Grant-Free Access for URLLC Service,” IEEE Journal on Selected Areas in Communications, vol. 39, no. 3, pp. 741–755, 2021.
  • [32] Y. Polyanskiy, H. V. Poor, and S. Verdu, “Channel Coding Rate in the Finite Blocklength Regime,” IEEE Transactions on Information Theory, vol. 56, no. 5, pp. 2307–2359, 2010.