Renewal-reward processes

[20/100%]
Show
LaTeX
Problem 5.36 from Gallager R.G. "Stochastic Processes: Theory for Applications" [Gal13]More problems in Chapter 5 "Renewal processes"

A large system is controlled by nn identical computers. Each computer independently alternates between an operational state and a repair state. The duration of the operational state, from completion of one repair until the next need for repair, is a random variable XX with finite expected duration E[X]\mathbb {E}\left[X\right]. The time required to repair a computer is an exponentially distributed random variable with density λe−λt\lambda e^{-\lambda t}. All operating durations and repair durations are independent. Assume that all computers are in the repair state at time 0.

?
(a)

For a single computer, say the iith, do the epochs at which the computer enters the repair state form a renewal process? If so, find the expected inter-renewal interval.

(b)

Do the epochs at which it enters the operational state form a renewal process?

(c)

Find the fraction of time over which the iith computer is operational and explain what you mean by fraction of time.

(d)

Let Qi(t)Q_{i}(t) be the probability that the iith computer is operational at time tt and find lim⁡t→∞Qi(t)\lim_{t \rightarrow \infty } Q_{i}(t).

(e)

The system is in failure mode at a given time if all computers are in the repair state at that time. Do the epochs at which system failure modes begin form a renewal process?

(f)

Let P(t)\mathbb {P}\left(t\right) be the probability that the the system is in failure mode at time tt. Find lim⁡t→∞P(t)\lim_{t \rightarrow \infty } \mathbb {P}\left(t\right).

(g)

For δ\delta small, find the probability that the system enters failure mode in the interval (t,t+δ](t, t+\delta ] in the limit as t→∞t \rightarrow \infty.

(h)

Find the expected time between successive entries into failure mode.

(i)

Next assume that the repair time of each computer has an arbitrary density rather than exponential, but has a mean repair time of 1/λ1/\lambda. Do the epochs at which system failure modes begin form a renewal process?

(j)

Repeat (f) for the assumption in (i).

This problem provides another way of treating ensemble averages for renewal-reward problems. Assume for notational simplicity that XX is a continuous random variable.

?
(a)

Show that P(one or more arrivals in (τ,τ+δ))=m(τ+δ)−m(τ)−o(δ)\mathbb {P}\left(\text{one or more arrivals in } (\tau , \tau +\delta )\right) = m(\tau +\delta )-m(\tau )-\mathrm{o}(\delta ) where o(δ)≥0\mathrm{o}(\delta ) \geq 0 and lim⁡δ→0o(δ)/δ=0\lim_{\delta \rightarrow 0} \mathrm{o}(\delta )/\delta = 0.

(b)

Show that P(Z(t)∈[z,z+δ),X~(t)∈(x,x+δ))\mathbb {P}\left(Z(t) \in [z, z+\delta ), \widetilde{X}(t) \in (x, x+\delta )\right) is equal to [m(t−z)−m(t−z−δ)−o(δ)][FX(x+δ)−FX(x)][m(t-z)-m(t-z-\delta )-\mathrm{o}(\delta )][\mathsf{F}_{X}(x+\delta )-\mathsf{F}_{X}(x)] for x≥z+δx \geq z+\delta.

(c)

Assuming that m′(τ)=dm(τ)/dτm'(\tau ) = dm(\tau )/d\tau exists for all τ\tau, show that the joint density of Z(t),X~(t)Z(t), \widetilde{X}(t) is fZ(t),X~(t)(z,x)=m′(t−z)fX(x)\mathsf{f}_{Z(t), \widetilde{X}(t)}(z, x) = m'(t-z) \mathsf{f}_{X}(x) for x>zx>z.

(d)

Show that E[R(t)]=∫z=0t∫x=z∞R(z,x)fX(x) dx m′(t−z) dz\mathbb {E}\left[R(t)\right] = \int_{z=0}^{t} \int_{x=z}^{\infty } \mathcal{R}(z, x) \mathsf{f}_{X}(x) \, dx \, m'(t-z) \, dz.

This problem is designed to give you an alternate way of looking at ensemble averages for renewal-reward problems. First we find an exact expression for P(SN(t)>s)\mathbb {P}\left(S_{N(t)}>s\right). We find this for arbitrary ss and tt, 0<s<t0<s<t.

?
(a)

By breaking the event {SN(t)>s}\left\{ S_{N(t)}>s\right\} into subevents {SN(t)>s,N(t)=n}\left\{ S_{N(t)}>s, N(t) = n\right\}, explain each of the following steps:

P(SN(t)>s)=∑n=1∞P(t≥Sn>s,Sn+1>t)=∑n=1∞∫y=stP(Sn+1>t∣Sn=y) dFSn(y)=∫y=stFXc(t−y) d∑n=1∞FSn(y)=∫y=stFXc(t−y) dm(y)where m(y)=E[N(y)]. \begin{aligned} \mathbb {P}\left(S_{N(t)}>s\right) & = \sum _{n=1}^{\infty } \mathbb {P}\left(t \geq S_{n}>s, S_{n+1}>t\right) \\ & = \sum _{n=1}^{\infty } \int _{y=s}^{t} \mathbb {P}\left(S_{n+1}>t \mid S_{n} = y\right) \, d\mathsf{F}_{S_{n}}(y) \\ & = \int _{y=s}^{t} \mathsf{F}_{X}^{c}(t-y) \, d\sum _{n=1}^{\infty } \mathsf{F}_{S_{n}}(y) \\ & = \int _{y=s}^{t} \mathsf{F}_{X}^{c}(t-y) \, dm(y) \qquad \text{where } m(y) = \mathbb {E}\left[N(y)\right]. \end{aligned}
(b)

Show that for 0<s<t<u0<s<t<u,

P(SN(t)>s,SN(t)+1>u)=∫y=stFXc(u−y) dm(y).(A.116) \mathbb {P}\left(S_{N(t)}>s, S_{N(t)+1}>u\right) = \int _{y=s}^{t} \mathsf{F}_{X}^{c}(u-y) \, dm(y). \tag {A.116}
(c)

Draw a two dimensional sketch, with age and duration as the axes, and show the region of (age, duration) values corresponding to the event {SN(t)>s,SN(t)+1>u}\left\{ S_{N(t)}>s, S_{N(t)+1}>u\right\}.

(d)

Assume that for large tt, dm(y)dm(y) can be approximated (according to Blackwell) as (1/X‾) dy(1/\overline{X}) \, dy, where X‾=E[X]\overline{X} = \mathbb {E}\left[X\right]. Assuming that XX also has a density, use the result in parts (b) and (c) to find the joint density of age and duration.

Consider a ferry that carries cars across a river. The ferry holds an integer number kk of cars and departs the dock when full. At that time, a new ferry immediately appears and begins loading newly arriving cars ad infinitum. The ferry business has been good, but customers complain about the long wait for the ferry to fill up.

?
(a)

Assume that cars arrive according to a renewal process. The IID interarrival times have mean X‾\overline{X}, variance σ2\sigma^{2} and moment generating function gX(r)\mathsf{g}_{X}(r). Does the sequence of departure times of the ferries form a renewal process? Explain carefully.

(b)

Find the expected time that a customer waits, starting from its arrival at the ferry terminal and ending at the departure of its ferry. Note 1: Part of the problem here is to give a reasonable definition of the expected customer waiting time. Note 2: It might be useful to consider k=1k = 1 and k=2k = 2 first.

(c)

Is there a `slow truck' phenomenon (a dependence on E[X2]\mathbb {E}\left[X^{2}\right]) here? Give an intuitive explanation. Hint: Look at k=1k = 1 and k=2k = 2 again.

(d)

In an effort to decrease waiting, the ferry managers institute a policy where no customer ever has to wait more than one hour. Thus, the first customer to arrive after a ferry departure waits for either one hour or the time at which the ferry is full, whichever comes first, and then the ferry leaves and a new ferry starts to accumulate new customers. Does the sequence of ferry departures form a renewal process under this new system? Does the sequence of times at which each successive empty ferry is entered by its first customer form a renewal process? You can assume here that t=0t = 0 is the time of the first arrival to the first ferry. Explain carefully.

(e)

Give an expression for the expected waiting time of the first new customer to enter an empty ferry under this new strategy.

Planes land at Heathrow airport at the times of a renewal process with interarrival time distribution function FF. Each plane contains a random number of people with a given common distribution and finite mean. Assuming as much independence as usual, find an expression for the rate of arrival of passengers over a long time period.

?

If X(t)X(t) is an irreducible positive recurrent Markov chain, and u(⋅)u(\cdot ) is a bounded function on the integers, show that

1t∫0tu(X(s))ds→ a.s. ∑i∈Sπiu(i) \frac{1}{t} \int _{0}^{t} u(X(s)) d s \xrightarrow {\text{ a.s. }} \sum _{i \in S} \pi _{i} u(i)

where π\mathbf{\pi } is the stationary distribution of X(t)X(t).

?

Let M(t)M(t) be an alternating renewal process, with interarrival pairs {Xr,Yr:r≥1}\left\{ X_{r}, Y_{r}: r \geq 1\right\}. Show that

1t∫0tI{M(s) is even }ds→ a.s. E[X1]E[X1]+E[Y1] as t→∞ \frac{1}{t} \int _{0}^{t} I_{\left\{ M(s) \text{ is even }\right\} } d s \xrightarrow {\text{ a.s. }} \frac{\mathbb {E}\left[X_{1}\right]}{\mathbb {E}\left[X_{1}\right]+\mathbb {E}\left[Y_{1}\right]} \quad \text{ as } t \rightarrow \infty

Is this limit valid for an arbitrary joint distribution of independent pairs (Xi,Yi)\left(X_{i}, Y_{i}\right)?

?

Let C(s)C(s) be the current lifetime (or age) of a renewal process N(t)N(t) with a typical interarrival time XX. Show that

1t∫0tC(s)ds→ a.s. E[X2]2E[X] as t→∞ \frac{1}{t} \int _{0}^{t} C(s) d s \xrightarrow {\text{ a.s. }} \frac{\mathbb {E}\left[X^{2}\right]}{2 \mathbb {E}\left[X\right]} \quad \text{ as } t \rightarrow \infty

Find the corresponding limit for the excess lifetime.

?

Let jj and kk be distinct states of an irreducible discrete-time Markov chain XX with stationary distribution π\pi. Show that

P(Tj<Tk∣X0=k)=1/πkE[Tj∣X0=k]+E[Tk∣X0=j] \mathbb {P}\left(T_{j} < T_{k} \mid X_{0} = k\right) = \frac{1 / \pi _{k}}{\mathbb {E}\left[T_{j} \mid X_{0} = k\right]+\mathbb {E}\left[T_{k} \mid X_{0} = j\right]}

where Ti=min⁡{n≥1:Xn=i}T_{i} = \min \left\{ n \geq 1: X_{n} = i\right\} is the first passage time to the state ii. [Hint: Consider the times of return to jj having made an intermediate visit to kk.]

?

Use the result of Exercise (10.5.2) to show that the limiting distribution function as t→∞t \rightarrow \infty of the total life D(t)D(t) of a renewal process is

FD(y)=∫0y1μxf(x)dx F_{D}(y) = \int _{0}^{y} \frac{1}{\mu } x f(x) d x

where ff is the density function of a typical interarrival time XX, and μ=E[X]\mu = \mathbb {E}\left[X\right]. The integrand is called the size-biased (or length-biased) density for XX.

?
?
(a)

Let XX be exponentially distributed with parameter λ\lambda. Show that

E[min⁡{X,d}]=1λ(1−e−λd),d≥0 \mathbb {E}\left[\min \left\{ X, d\right\} \right] = \frac{1}{\lambda }\left(1-e^{-\lambda d}\right), \quad d \geq 0
(b)

John's garage offers him the choice of two tyre replacement plans.

  1. The garage undertakes to replace all the tyres of his car at the normal price whenever one of the tyres requires replacing.

  2. The garage undertakes to replace all the tyres of his car at 5%5 \% of the normal price two years after they have last been replaced. However, if any tyre needs replacing earlier, the garage will then replace all the tyres at a price that is 5%5 \% higher than the normal price. Assuming that a new tyre has an exponentially distributed lifetime with mean 8 years, determine the long run average cost per year under the two options. Which option should John choose if his car is new?

A machine MM is repaired at time t=0t = 0, and its uptime after any repair is exponentially distributed with parameter λ\lambda, at which point it breaks down (assume the usual independence). Following any repair at time TT, say, it is inspected at times T,T+m,T+2m,…T, T+m, T+2 m, \ldots, and instantly repaired if found to be broken (the inspection schedule is then restarted). Show that the long run proportion of time that MM is working (the 'uptime ratio') is m−1∫0me−λxdxm^{-1} \int_{0}^{m} e^{-\lambda x} d x.

?

The interarrival times of this process are Z0,Y1,Z1,Y2,…Z_{0}, Y_{1}, Z_{1}, Y_{2}, \ldots, where the YiY_{i} and ZjZ_{j} are independent with respective common moment generating functions MYM_{Y} and MZM_{Z}. Let p(t)p(t) be the probability that the epoch tt of time lies in an interval of type ZZ. Show that the Laplace-Stieltjes transform p∗p^{*} of pp satisfies

p∗(θ)=1−MZ(−θ)1−MY(−θ)MZ(−θ) p^{*}(\theta ) = \frac{1-M_{Z}(-\theta )}{1-M_{Y}(-\theta ) M_{Z}(-\theta )}
?

Bears arrive in a village at the instants of a renewal process; they are captured and confined at a cost of $c\$ c per unit time per bear. When a given number BB bears have been captured, an expedition (costing $d\$ d) is organized to remove and release them a long way away. What is the long-run average cost of this policy?

?

Let X={Xn:n≥0}X = \left\{ X_{n}: n \geq 0\right\} be an ergodic Markov chain on a finite state space SS, with stationary distribution π\pi. Let k∈Sk \in S, and let TT be a strictly positive stopping time with finite mean such that XT=kX_{T} = k. For j∈Sj \in S, let Vj(k)V_{j}(k) be the number of visits to jj by the chain started in kk and stopped at time TT (set Vk(k)=1V_{k}(k) = 1). Use the renewal-reward theorem to show that E[Vj(k)]=πjEk[T]\mathbb {E}\left[V_{j}(k)\right] = \pi_{j} \mathbb {E}_{k}\left[T\right].

?

For a given type of car, costing cc new, the number of years between its manufacture and being broken up for scrap is X(≥1)X( \geq 1), where XX is a random variable with distribution function FF. In its kk th year after manufacture (with k≥1k \geq 1), depreciation has reduced its resale value to cλkc \lambda^{k}, and in that year it costs rμk−1r \mu^{k-1} in repairs and maintenance, where μ≠1\mu \neq 1. (Usually, λ<1\lambda < 1 and μ>1\mu > 1.)

You buy a new car immediately your old car either is scrapped or reaches the age of mm years, whichever is the sooner. Show that, if you continue the policy indefinitely, in the long run you minimize your expected average cost by choosing mm such that

1E[Y]{c+r1−μ−cG(λ)−r1−μG(μ)} \frac{1}{\mathbb {E}\left[Y\right]}\left\{ c+\frac{r}{1-\mu }-c G(\lambda )-\frac{r}{1-\mu } G(\mu )\right\}

is as small as possible, where G(s)=E[sY]G(s) = \mathbb {E}\left[s^{Y}\right] and Y=min⁡{m,X}Y = \min \left\{ m, X\right\}.

?

Consider the M(λ)/G/1M(\lambda ) / G / 1 queue with no waiting room. Customers who arrive while the server is busy are lost. Show that the long run proportion of arrivals lost is 1/(1+ρ−1)1 /\left(1+\rho^{-1}\right) where ρ=λE[S]\rho = \lambda \mathbb {E}\left[S\right] is the traffic intensity.

?

Consider an M(λ)/G/∞M(\lambda ) / G / \infty queue. Use the renewal-reward theorem to show that the expected duration of a busy period is (eρ−1)/λ\left(e^{\rho }-1\right) / \lambda where ρ=λE[S]\rho = \lambda \mathbb {E}\left[S\right].

?

After the level of an industrial process has been set at its desired value, it wanders in a random fashion. To counteract this the process is periodically reset to this desired value, at times 0,T,2T,…0, T, 2 T, \ldots. If WtW_{t} is the deviation from the desired level, tt units of time after a reset, then {Wt:0≤t<T}\left\{ W_{t}: 0 \leq t < T\right\} can be modelled by a standard Wiener process. The behaviour of the process after a reset is independent of its behaviour before the reset. While WtW_{t} is outside the range (−a,a)(-a, a) the output from the process is unsatisfactory and a cost is incurred at rate CC per unit time. The cost of each reset is RR. Show that the period TT which minimises the long-run average cost per unit time is T∗T^{*}, where

R=C∫0T∗a(2πt)exp⁡(−a22t)dt R = C \int _{0}^{T^{*}} \frac{a}{\sqrt{(2 \pi t)}} \exp \left(-\frac{a^{2}}{2 t}\right) d t
?

Let X1,X2,…X_1, X_2, \ldots be the i.i.d. lifetimes of some component in some large machine. The simplest replacement policy is to change a component as soon as it fails. In this case it may be necessary to call a repairman at night, which might be costly. Another policy, called replacement based on age, is to replace at failure or at some given age, aa, say, whichever comes first, in which case the interreplacement times are

Wk=min⁡{Xk,a},k≥1. W_k = \min \left\{ X_k, a\right\} , \quad k \geq 1.

Suppose that c1c_1 is the cost for replacements due to failure and that c2c_2 is the cost for replacements due to age. In addition, let YkY_k be the cost attached to replacement kk, k≥1k \geq 1, and let N(t)N(t) be the number of replacements made in the time interval (0,t](0,t], where {N(t),t≥0}\left\{ N(t), t \geq 0\right\} is a Poisson process, which is independent of X1,X2,…X_1, X_2, \ldots. This means that

Z(t)=∑k=1N(t)Yk Z(t) = \sum _{k=1}^{N(t)} Y_k

is the total cost caused by the replacements in the time interval (0,t](0,t] (with Z(t)=0Z(t)=0 whenever N(t)=0N(t)=0).

?
(a)

Compute E[Y1]\mathbb {E}\left[Y_1\right] and Var⁡[Y1]\operatorname {Var}\left[Y_1\right].

(b)

Compute E[Z(t)]\mathbb {E}\left[Z(t)\right] and Var⁡[Z(t)]\operatorname {Var}\left[Z(t)\right].

References

[Gal13] Gallager R.G. Stochastic Processes: Theory for Applications. Cambridge University Press, 2013.
[GrSt20] Grimmett G.R., Stirzaker D.R. Probability and Random Processes. 4th ed., Oxford University Press, 2020.
[Gut09] Gut A. An Intermediate Course in Probability. 2nd ed., Springer, 2009.