MathIsimple

Stochastic Processes – Problem 8: Prove that

Question

Let X1,XniidXX_{1},\cdots X_{n}\stackrel{iid}{\sim}X, and denote Sn=i=1nXiS_{n}=\sum_{i=1}^{n}X_{i}, N(t)=inf{n:Sn>t}N(t)=\inf\{n:S_{n}>t\}, so that N(t)N(t) is a renewal process.

(1) Prove that P(XN(t)>x)P(X>x),x0\mathbb{P}(X_{N(t)}>x)\geq\mathbb{P}(X>x),\forall x\geq0.

(2) If the distribution function of XX is F(x)=1exF(x)=1-e^{-x}, find the exact analytical expression for P(XN(t)>x)\mathbb{P}(X_{N(t)}>x).

(3) Using the Key Renewal Theorem, find limtEXN(t)\lim_{t\to\infty}\mathbb{E}X_{N(t)}.

Step-by-step solution

(1) Step 1. Denote H(t,x)=P(XN(t)>x)H(t, x) = \mathbb{P}(X_{N(t)} > x). By conditioning on the first renewal epoch X1X_1, we can establish a renewal equation for H(t,x)H(t, x). When X1>tX_1 > t, N(t)=1N(t) = 1, so XN(t)=X1X_{N(t)} = X_1, giving P(XN(t)>xX1>t)=P(X1>xX1>t)\mathbb{P}(X_{N(t)} > x | X_1 > t) = \mathbb{P}(X_1 > x | X_1 > t). When X1tX_1 \leq t, the process restarts from tX1t - X_1, so XN(t)X_{N(t)} has the same distribution as XN(tX1)X_{N(t-X_1)}. Thus: H(t,x)=P(X1>max(t,x))+0tH(ty,x)dF(y).H(t, x) = \mathbb{P}(X_1 > \max(t, x)) + \int_0^t H(t-y, x) dF(y).

Step 2. Denote Fˉ(x)=1F(x)=P(X>x)\bar{F}(x) = 1 - F(x) = \mathbb{P}(X > x). Let Δ(t)=H(t,x)Fˉ(x)\Delta(t) = H(t, x) - \bar{F}(x). Substituting H(t,x)=Δ(t)+Fˉ(x)H(t, x) = \Delta(t) + \bar{F}(x) into the renewal equation and simplifying: Δ(t)=z(t)+0tΔ(ty)dF(y),\Delta(t) = z(t) + \int_0^t \Delta(t-y) dF(y), where the source term is z(t)=Fˉ(max(t,x))Fˉ(x)+Fˉ(x)F(t)z(t) = \bar{F}(\max(t, x)) - \bar{F}(x) + \bar{F}(x)F(t).

Step 3. Analyze the sign of z(t)z(t). When t<xt < x: z(t)=Fˉ(x)F(t)0z(t) = \bar{F}(x)F(t) \geq 0. When txt \geq x: z(t)=Fˉ(t)F(x)0z(t) = \bar{F}(t)F(x) \geq 0. Thus z(t)0z(t) \geq 0 for all t0t \geq 0.

Step 4. By the renewal equation solution formula, Δ(t)=0tz(ty)dU(y)\Delta(t) = \int_0^t z(t-y) dU(y), where U(t)=n=0Fn(t)U(t) = \sum_{n=0}^\infty F^{*n}(t) is the renewal function. Since z(t)0z(t) \geq 0 and dU(y)dU(y) is a nonnegative measure, Δ(t)0\Delta(t) \geq 0. Thus H(t,x)Fˉ(x)H(t, x) \geq \bar{F}(x), i.e., P(XN(t)>x)P(X>x)\mathbb{P}(X_{N(t)} > x) \geq \mathbb{P}(X > x). QED.

(2) Since F(x)=1exF(x) = 1 - e^{-x}, XX follows an exponential distribution with parameter 11. The renewal process corresponding to the exponential distribution is a Poisson process, with renewal function U(t)=1+tU(t) = 1 + t (including the unit step at n=0n=0). From the derivation in (1), P(XN(t)>x)\mathbb{P}(X_{N(t)} > x) satisfies a renewal equation whose solution is: H(t,x)=Fˉ(max(t,x))+0tFˉ(max(ty,x))dy.H(t, x) = \bar{F}(\max(t, x)) + \int_0^t \bar{F}(\max(t-y, x)) dy. When x>tx > t: H(t,x)=ex(1+t).H(t, x) = e^{-x}(1 + t). When xtx \leq t: H(t,x)=ex(1+x).H(t, x) = e^{-x}(1 + x). In summary, P(XN(t)>x)=ex(1+min(x,t))\mathbb{P}(X_{N(t)} > x) = e^{-x}(1 + \min(x, t)).

(3) Denote g(t)=EXN(t)g(t) = \mathbb{E} X_{N(t)}. Conditioning on X1X_1, we establish a renewal equation for g(t)g(t): g(t)=tydF(y)+0tg(ty)dF(y).g(t) = \int_t^\infty y dF(y) + \int_0^t g(t-y) dF(y). Let h(t)=tydF(y)h(t) = \int_t^\infty y dF(y). This is a standard renewal equation g(t)=h(t)+0tg(ty)dF(y)g(t) = h(t) + \int_0^t g(t-y) dF(y). Verify integrability: 0h(t)dt=E[X2].\int_0^\infty h(t) dt = \mathbb{E}[X^2]. Assuming E[X2]<\mathbb{E}[X^2] < \infty, h(t)h(t) is directly Riemann integrable. By the Key Renewal Theorem, limtg(t)=1μ0h(t)dt=E[X2]E[X],\lim_{t \to \infty} g(t) = \frac{1}{\mu} \int_0^\infty h(t) dt = \frac{\mathbb{E}[X^2]}{\mathbb{E}[X]}, where μ=E[X]\mu = \mathbb{E}[X].

Final answer

(1) QED.; (2) P(XN(t)>x)=ex(1+min(x,t))\mathbb{P}(X_{N(t)} > x) = e^{-x}(1 + \min(x, t)); (3) E[X2]E[X]\dfrac{\mathbb{E}[X^2]}{\mathbb{E}[X]}

Marking scheme

The following is the rubric designed for this problem and its official solution.

1. Checkpoints (Total 7 pts)

Part (1) (3 pts)

  • Establish the renewal equation for H(t,x)H(t,x) [1 pt]: Correctly write the renewal equation satisfied by P(XN(t)>x)\mathbb{P}(X_{N(t)}>x), namely H(t,x)=P(X1>max(t,x))+0tH(ty,x)dF(y)H(t, x) = \mathbb{P}(X_1 > \max(t, x)) + \int_0^t H(t-y, x) dF(y). Note: The source term P(X1>max(t,x))\mathbb{P}(X_1 > \max(t, x)) must correctly reflect the relationship between the first renewal epoch X1X_1 and both tt and xx.
  • Analyze the sign of the difference function or source term [1 pt]: Construct the difference Δ(t)=H(t,x)P(X>x)\Delta(t) = H(t,x) - \mathbb{P}(X>x) and derive its renewal equation, OR directly analyze the iteration properties of the original equation. Must show the key algebraic steps proving the effective source term (such as z(t)z(t) in the solution) is nonnegative in both cases t<xt<x and txt \ge x.
  • Derive the inequality conclusion [1 pt]: Based on the nonnegativity of the source term and the sign-preserving property of the renewal function (or renewal equation solution), rigorously conclude Δ(t)0\Delta(t) \ge 0 or H(t,x)P(X>x)H(t,x) \ge \mathbb{P}(X>x).

Part (2) (2 pts)

  • Exponential distribution setup and integral formulation [1 pt]: Identify the renewal process properties corresponding to the exponential distribution (such as U(t)=1+tU(t)=1+t or dU(y)=δ(y)dy+dydU(y) = \delta(y)dy + dy), and correctly establish the specific integral expression. If the specific form of U(t)U(t) is not substituted and only the general formula is given, no credit.
  • Compute the piecewise analytical expression [1 pt]: Correctly compute the integrals for both cases x>tx > t and xtx \le t, and write the final piecewise result (or unified form ex(1+min(x,t))e^{-x}(1+\min(x,t))). Both parts must be correct for credit; only one case correct receives no credit.

Part (3) (2 pts)

  • Mean renewal equation and source term integral [1 pt]: Write the renewal equation for g(t)=EXN(t)g(t) = \mathbb{E}X_{N(t)}, and verify that the integral of its source term h(t)=tydF(y)h(t)=\int_t^\infty y dF(y) over (0,)(0,\infty) equals E[X2]\mathbb{E}[X^2]. Key point: must show the computation of 0h(t)dt\int_0^\infty h(t)dt or explain why it equals the second moment.
  • Apply the Key Renewal Theorem (KRT) [1 pt]: Cite the Key Renewal Theorem, derive the limit as h(t)dtE[X]\frac{\int h(t)dt}{\mathbb{E}[X]}, and give the final result E[X2]E[X]\frac{\mathbb{E}[X^2]}{\mathbb{E}[X]}.

Total (max 7)

2. Zero-credit items

  • Part (1): Only intuitively describing that "the interval covering tt tends to be longer than a typical interval" (length-biased sampling intuition) without providing a rigorous mathematical proof (such as renewal equation or coupling argument), receives 0 pts.
  • Part (1): Only verifying the cases t=0t=0 or tt \to \infty without proving for all tt, receives 0 pts.
  • Part (3): Without deriving via the renewal theorem, directly writing the expectation formula for the "length-biased distribution" without derivation, receives 0 pts (the problem explicitly requires "using the Key Renewal Theorem").

3. Deductions

  • Concept confusion: Mistakenly computing XN(t)X_{N(t)} (the total lifetime covering tt) as the residual life (SN(t)tS_{N(t)} - t) or the current age, leading to incorrect results in Part (2) or (3) (e.g., obtaining the residual life limit E[X2]2E[X]\frac{\mathbb{E}[X^2]}{2\mathbb{E}[X]}): that part receives 0 pts.
  • Lack of existence verification: In Part (3), if E[X2]<\mathbb{E}[X^2] < \infty or Riemann integrability conditions are not mentioned, no deduction.
  • Sign error: If the logic is reversed when proving the inequality (e.g., assuming the conclusion first then deriving), deduct 1 pt.
  • Integration computation error: Pure arithmetic/integration errors in Part (2) or (3) with completely correct approach, deduct 1 pt.
Ask AI ✨