Chapter 8

Random processes

[49/98%]
Show
LaTeX
Problem 8.2.1

Let {Xn}\left\{ X_{n}\right\} be a Markov chain on the state space S={0,1}S = \left\{ 0,1\right\} with transition matrix

P=[1−ααβ1−β] \mathbf{P} =\left[\begin{smallmatrix} 1-\alpha & \alpha \\ \beta & 1-\beta \end{smallmatrix}\right]

where α+β>0\alpha +\beta > 0. Find:

?
(a)

the correlation ρ(Xm,Xm+n)\rho \left(X_{m}, X_{m+n}\right), and its limit as m→∞m \rightarrow \infty with nn remaining fixed,

(b)

lim⁡n→∞n−1∑r=1nP(Xr=1)\lim_{n \rightarrow \infty } n^{-1} \sum_{r = 1}^{n} \mathbb {P}\left(X_{r} = 1\right), and the condition under which the process is strongly stationary.

Problem 8.2.2

Let {N(t):t≥0}\left\{ N(t): t \geq 0\right\} be a Poisson process of intensity λ\lambda, and let T0T_{0} be an independent random variable such that P(T0=±1)=12\mathbb {P}\left(T_{0} = \pm 1\right) = \frac{1}{2}. Define T(t)=T0(−1)N(t)T(t) = T_{0}(-1)^{N(t)}. Show that {T(t):t≥0}\left\{ T(t): t \geq 0\right\} is stationary and find:

?
(a)

ρ(T(s),T(s+t))\rho (T(s), T(s+t)),

(b)

the mean and variance of X(t)=∫0tT(s)dsX(t) = \int_{0}^{t} T(s) d s. [The so-called Goldstein-Kac process X(t)X(t) denotes the position of a particle moving with unit speed, starting from the origin along the positive xx-axis, whose direction is reversed at the instants of a Poisson process.]

Problem 8.2.3

An integer-valued counting process {N(t):t≥0}\left\{ N(t): t \geq 0\right\} with N(0)=N(0) = 0 is called crudely stationary if pk(s,t)=P(N(s+t)−N(s)=k)p_{k}(s, t) = \mathbb {P}\left(N(s+t)-N(s) = k\right) depends only on the length t−st-s and not on the location ss. It is called simple if, almost surely, it has jump discontinuities of size 1 only. Show that, for a simple crudely stationary process N,lim⁡t↓0t−1P(N(t)>0)=E[N(1)]N, \lim_{t \downarrow 0} t^{-1} \mathbb {P}\left(N(t) > 0\right) = \mathbb {E}\left[N(1)\right].

?
Problem 8.2.4

Let X={Xn:n≥0}X = \left\{ X_{n}: n \geq 0\right\} be an ergodic Markov chain started in its stationary distribution π\pi. Let AA be a subset of states, with stationary probability π(A)\pi (A), and let TA=min⁡{n≥1:Xn∈A}T_{A} = \min \left\{ n \geq 1: X_{n} \in A\right\}. Show that E[TA∣X0∈A]=1/π(A)\mathbb {E}\left[T_{A} \mid X_{0} \in A\right] = 1 / \pi (A).

?
Problem 8.3.1

Let (fn:n≥1)\left(f_{n}: n \geq 1\right) be a probability distribution on the positive integers, and define a sequence (un:n≥0)(u_{n}: n \geq 0) by u0=1u_{0} = 1 and un=∑r=1nfrun−r,n≥1u_{n} = \sum_{r = 1}^{n} f_{r} u_{n-r}, n \geq 1. Explain why such a sequence is called a renewal sequence, and show that uu is a renewal sequence if and only if there exists a Markov chain UU and a state ss such that un=P(Un=s∣U0=s)u_{n} = \mathbb {P}\left(U_{n} = s \mid U_{0} = s\right).

?
Problem 8.3.2

Let {Xi:i≥1}\left\{ X_{i}: i \geq 1\right\} be the inter-event times of a discrete renewal process on the integers. Show that the excess lifetime BnB_{n} constitutes a Markov chain. Write down the transition probabilities of the sequence {Bn}\left\{ B_{n}\right\} when reversed in equilibrium. Compare these with the transition probabilities of the chain UU of your solution to Exercise (8.3.1).

?
Problem 8.3.3

Let (un:n≥1)(u_{n}: n \geq 1) satisfy u0=1u_{0} = 1 and un=∑r=1nfrun−ru_{n} = \sum_{r = 1}^{n} f_{r} u_{n-r} for n≥1n \geq 1, where (fr:r≥1)\left(f_{r}: r \geq 1\right) is a non-negative sequence. Show that:

?
(a)

vn=ρnunv_{n} = \rho^{n} u_{n} is a renewal sequence if ρ>0\rho > 0 and ∑n=1∞ρnfn=1\sum_{n = 1}^{\infty } \rho^{n} f_{n} = 1,

(b)

as n→∞,ρnunn \rightarrow \infty , \rho^{n} u_{n} converges to some constant cc.

Problem 8.3.4

Events occur at the times of a discrete-time renewal process NN (see Example (5.2.15)). Let unu_{n} be the probability of an event at time nn, with generating function U(s)U(s), and let F(s)F(s) be the probability generating function of a typical inter-event time. Show that, if ∣s∣<1\left|s\right| < 1 :

∑r=0∞E[N(r)]sr=F(s)U(s)1−s and ∑t=0∞E[(N(t)+kk)]st=U(s)k1−s for k≥0 \sum _{r = 0}^{\infty } \mathbb {E}\left[N(r)\right] s^{r} = \frac{F(s) U(s)}{1-s} \quad \text{ and } \quad \sum _{t = 0}^{\infty } \mathbb {E}\left[\binom {N(t)+k}{k}\right] s^{t} = \frac{U(s)^{k}}{1-s} \quad \text{ for } k \geq 0
?
Problem 8.3.5

Prove Theorem (8.3.5): Poisson processes are the only renewal processes that are Markov chains.

?
Problem 8.4.1

The two tellers in a bank each take an exponentially distributed time to deal with any customer; their parameters are λ\lambda and μ\mu respectively. You arrive to find exactly two customers present, each occupying a teller.

?
(a)

You take a fancy to a randomly chosen teller, and queue for that teller to be free; no later switching is permitted. Assuming any necessary independence, what is the probability pp that you are the last of the three customers to leave the bank?

(b)

If you choose to be served by the quicker teller, find pp.

(c)

Suppose you go to the teller who becomes free first. Find pp.

Problem 8.4.2

Customers arrive at a desk according to a Poisson process of intensity λ\lambda. There is one clerk, and the service times are independent and exponentially distributed with parameter μ\mu. At time 0 there is exactly one customer, currently in service. Show that the probability that the next customer arrives before time tt and finds the clerk busy is

λλ+μ(1−e−(λ+μ)t) \frac{\lambda }{\lambda +\mu }\left(1-e^{-(\lambda +\mu ) t}\right)
?
Problem 8.4.3

Vehicles pass a crossing at the instants of a Poisson process of intensity λ\lambda; you need a gap of length at least aa in order to cross. Let TT be the first time at which you could succeed in crossing to the other side. Show that E[T]=(eaλ−1)/λ\mathbb {E}\left[T\right] = \left(e^{a \lambda }-1\right) / \lambda, and find E[eθT]\mathbb {E}\left[e^{\theta T}\right].

Suppose there are two lanes to cross, carrying independent Poissonian traffic with respective rates λ\lambda and μ\mu. Find the expected time to cross in the two cases when:

?
(a)

there is an island or refuge between the two lanes,

(b)

you must cross both in one go. Which is the greater?

Problem 8.4.4

Customers arrive at the instants of a Poisson process of intensity λ\lambda, and the single server has exponential service times with parameter μ\mu. An arriving customer who sees nn customers present (including anyone in service) will join the queue with probability (n+1)/(n+2)(n+1) /(n+2), otherwise leaving for ever. Under what condition is there a stationary distribution? Find the mean of the time spent in the queue (not including service time) by a customer who joins it when the queue is in equilibrium. What is the probability that an arrival joins the queue when in equilibrium?

?
Problem 8.4.5

Customers enter a shop at the instants of a Poisson process of rate 2. At the door, two representatives separately demonstrate a new corkscrew. This typically occupies the time of a customer and the representative for a period which is exponentially distributed with parameter 11, independently of arrivals and other demonstrators. If both representatives are busy, customers pass directly into the shop. No customer passes a free representative without being stopped, and all customers leave by another door. If both representatives are free at time 00, show the probability that both are busy at time tt is 25−23e−2t+415e−5t\frac{2}{5}-\frac{2}{3} e^{-2 t}+\frac{4}{15} e^{-5 t}.

?
Problem 8.5.1

For a Wiener process WW with W(0)=0W(0) = 0, show that

P(W(s)>0,W(t)>0)=14+12πsin⁡−1st for s<t \mathbb {P}\left(W(s) > 0, W(t) > 0\right) = \frac{1}{4}+\frac{1}{2 \pi } \sin ^{-1} \sqrt{\frac{s}{t}} \text{ for } s < t

Calculate P(W(s)>0,W(t)>0,W(u)>0)\mathbb {P}\left(W(s) > 0, W(t) > 0, W(u) > 0\right) when s<t<us < t < u.

?
Problem 8.5.2

Let WW be a Wiener process. Show that, for s<t<us < t < u, the conditional distribution of W(t)W(t) given W(s)W(s) and W(u)W(u) is normal

N((u−t)W(s)+(t−s)W(u)u−s,(u−t)(t−s)u−s) N\left(\frac{(u-t) W(s)+(t-s) W(u)}{u-s}, \frac{(u-t)(t-s)}{u-s}\right)

Deduce that the conditional correlation between W(t)W(t) and W(u)W(u), given W(s)W(s) and W(v)W(v), where s<t<u<vs < t < u < v, is

(v−u)(t−s)(v−t)(u−s) \sqrt{\frac{(v-u)(t-s)}{(v-t)(u-s)}}
?
Problem 8.5.3

For what values of aa and bb is aW1+bW2a W_{1}+b W_{2} a standard Wiener process, where W1W_{1} and W2W_{2} are independent standard Wiener processes?

?
Problem 8.5.4

Show that a Wiener process WW with variance parameter σ2\sigma^{2} satisfies

∑j=0n−1{W((j+1)t/n)−W(jt/n)}2→ m.s. σ2t as n→∞ \sum _{j = 0}^{n-1}\left\{ W((j+1) t / n)-W(j t / n)\right\} ^{2} \xrightarrow {\text{ m.s. }} \sigma ^{2} t \quad \text{ as } n \rightarrow \infty

The process WW is said to have finite quadratic variation.

?
Problem 8.5.5

Let WW be a Wiener process. Which of the following define Wiener processes?

?
(a)

−W(t)-W(t),

(b)

tW(1)\sqrt{t} W(1),

(c)

W(2t)−W(t)W(2 t)-W(t).

Problem 8.5.6

Find the distribution of ∫0t[W(u)/u]du\int_{0}^{t}[W(u) / u] d u where WW is the Wiener process. [The integral is interpreted as the limit as ϵ↓0\epsilon \downarrow 0 of the integral from ϵ\epsilon to tt.]

?
Problem 8.5.7

Wiener process in nn dimensions. An nn-dimensional Wiener process is a process W(t)=(W1(t),W2(t),…,Wn(t))W(t) = \left(W_{1}(t), W_{2}(t), \ldots , W_{n}(t)\right) taking values in Rn\mathbb {R}^{n} with independent increments, started at the origin W(0)=0W(0) = 0, and such that W(t)−W(s)W(t)-W(s) is multivariate normal with means 0 and covariance matrix (t−s)I(t-s) \mathbf{I} where I\mathbf{I} is the n×nn \times n identity matrix. Let W(t)W(t) be an nn-dimensional Wiener process, and let AA be an orthonormal n×nn \times n matrix. Show that AW(t)A W(t) is an nn-dimensional Wiener process.

?
Problem 8.5.8

Let WW be a standard Wiener process. Show that, for p≥0p \geq 0,

E[∣W(t)−W(s)∣p]=cp∣t−s∣p/2,s,t∈R \mathbb {E}\left[|W(t)-W(s)|^{p}\right] = c_{p}|t-s|^{p / 2}, \quad s, t \in \mathbb {R}

where

cp=12π∫−∞∞∣y∣pe−12y2dy c_{p} = \frac{1}{\sqrt{2 \pi }} \int _{-\infty }^{\infty }|y|^{p} e^{-\frac{1}{2} y^{2}} d y
?
Problem 8.6.1

Prove that the characteristic function ϕ(t,θ)=E[eiθX(t)]\phi (t, \theta ) = \mathbb {E}\left[e^{i \theta X(t)}\right] of a Lévy process XX is a continuous function of tt.

?
Problem 8.6.2

Show that a compound Poisson process is a Lévy process, and find its Lévy symbol.

?
Problem 8.6.3

Verify the formula of Example (8.6.6) for the Laplace exponent of the Moran gamma process.

?
Problem 8.6.4

Prove Theorem (8.6.7): the property of being a Lévy process is preserved under time change by an independent subordinator.

?
Problem 8.6.5

Let NN be a Poisson process of rate 11, and TT an independent random variable with the Γ(1,t)\Gamma (1, t) distribution. Show that Y=N(T)Y = N(T) has probability generating function E[sY]=(2−s)−t\mathbb {E}\left[s^{Y}\right] = (2-s)^{-t} for s<2s < 2.

?
Problem 8.6.6

Let XX be a Lévy process with finite variances. Show that the following are martingales:

?
(a)

X(t)−E[X(t)]X(t)-\mathbb {E}\left[X(t)\right]

(b)

Z(t)2−E[Z(t)2]Z(t)^{2}-\mathbb {E}\left[Z(t)^{2}\right], where Z(t)=X(t)−E[X(t)]Z(t) = X(t)-\mathbb {E}\left[X(t)\right]

(c)

eiθX(t)/ϕ(t,θ)e^{i \theta X(t)} / \phi (t, \theta ), where ϕ(t,θ)=E[eiθX(t)]\phi (t, \theta ) = \mathbb {E}\left[e^{i \theta X(t)}\right].

Problem 8.6.7

Let XX be a continuous-time martingale and TT an independent subordinator. Show that, subject to the moment condition E[∣Y(t)∣]<∞,Y(t)=X(T(t))\mathbb {E}\left[\left|Y(t)\right|\right] < \infty , Y(t) = X(T(t)) defines a martingale with respect to a suitable filtration. Show that the moment condition is valid whenever XX is a positive martingale.

?
Problem 8.7.1

Let XX be stable, and let X1,X2,…,XnX_{1}, X_{2}, \ldots , X_{n} be a random sample distributed as XX. Show that, for n≥1n \geq 1, there exist An,BnA_{n}, B_{n} such that X1+X2+⋯+Xn=DBnX+AnX_{1}+X_{2}+\cdots +X_{n} \stackrel{\mathrm{D}}{ = } B_{n} X+A_{n}.

?
Problem 8.7.2

Prove that, for a non-degenerate, self-similar process XX with stationary increments, finite variances, and Hurst exponent H>0H > 0,

Var⁡[X(t)]=t2HVar⁡[X(1)],t≥0E[X(s)X(t)]=12(s2H+t2H−∣t−s∣2H)E[X(1)2],s,t≥0 \begin{aligned} \operatorname {Var}\left[X(t)\right] & = t^{2 H} \operatorname {Var}\left[X(1)\right], \quad t \geq 0 \\ \mathbb {E}\left[X(s) X(t)\right] & = \frac{1}{2}\left(s^{2 H}+t^{2 H}-|t-s|^{2 H}\right) \mathbb {E}\left[X(1)^{2}\right], \quad s, t \geq 0 \end{aligned}
?
Problem 8.7.3

Show that a stable distribution is infinitely divisible.

?
Problem 8.7.4

Show that a self-similar Lévy process either is a Wiener process or has infinite variances.

?
Problem 8.7.5

Prove that a distribution with characteristic function ϕ(θ)=exp⁡(−∣θ∣α)\phi (\theta ) = \exp \left(-|\theta |^{\alpha }\right), where α∈(0,2]\alpha \in (0,2], is strictly α\alpha-stable and symmetric.

?
Problem 8.7.6

Let XX and YY be independent stable random variables, where XX is symmetric with exponent α∈(0,2]\alpha \in (0,2], and YY is non-negative with Laplace transform M(θ)=exp⁡(−kθβ)M(\theta ) = \exp \left(-k \theta^{\beta }\right) for θ∈[0,∞)\theta \in [0, \infty ), where 1≠β∈(0,2]1 \neq \beta \in (0,2]. Show that Z=XY1/αZ = X Y^{1 / \alpha } is symmetric and αβ\alpha \beta-stable. Deduce that, for independent N(0,1)N(0,1) random variables U,VU, V, the ratio Z=U/VZ = U / V has the Cauchy distribution.

?
Problem 8.8.1

Let ZZ be a continuous-time Markov chain, and let TT be an independent subordinator. Show that X(t)=Z(T(t))X(t) = Z(T(t)) defines a Markov chain.

?
Problem 8.8.2

Let Z={Zn:n≥0}Z = \left\{ Z_{n}: n \geq 0\right\} be a discrete-time Markov chain with nn-step transition probabilities zij(n)z_{i j}(n). Let NN be a Poisson process of rate λ\lambda that is independent of ZZ, and set X(t)=ZN(t)X(t) = Z_{N(t)}. Show that XX is a continuous-time Markov chain with transition probabilities

pij(t)=∑n=0∞e−λt(λt)nn!zij(n) p_{i j}(t) = \sum _{n = 0}^{\infty } e^{-\lambda t} \frac{(\lambda t)^{n}}{n!} z_{i j}(n)
?
Problem 8.8.3

Let WW be the standard Wiener process, and let TT be an independent stable subordinator with E[e−θT(t)]=exp⁡(−tθa/2)\mathbb {E}\left[e^{-\theta T(t)}\right] = \exp \left(-t \theta^{a / 2}\right) where a∈(0,2)a \in (0,2) and θ≥0\theta \geq 0. Show that the time-changed process Y(t)=W(T(t))Y(t) = W(T(t)) is a symmetric Lévy process.

?
Problem 8.8.4

Markov, but not strong Markov. Let VV be exponentially distributed with parameter 1, and define X={X(t):t≥0}X = \left\{ X(t): t \geq 0\right\} by: X(t)=0X(t) = 0 for t≤Vt \leq V, and X(t)=t−VX(t) = t-V for t≥Vt \geq V.

?
(a)

Show that XX is a Markov process.

(b)

Show that VV is a stopping time for XX.

(c)

Show that the two processes X1(t)=X(t)X_{1}(t) = X(t) and X2(t)=X(t+V)X_{2}(t) = X(t+V) satisfy X1(0)=X2(0)=0X_{1}(0) = X_{2}(0) = 0, but have different fdds.

Problem 8.10.1

Let {Zn}\left\{ Z_{n}\right\} be a sequence of uncorrelated real-valued variables with zero means and unit variances, and define the 'moving average'

Yn=∑i=0rαiZn−i Y_{n} = \sum _{i = 0}^{r} \alpha _{i} Z_{n-i}

for constants α0,α1,…,αr\alpha_{0}, \alpha_{1}, \ldots , \alpha_{r}. Show that YY is stationary and find its autocovariance function.

?
Problem 8.10.2

Let {Zn}\left\{ Z_{n}\right\} be a sequence of uncorrelated real-valued variables with zero means and unit variances. Suppose that {Yn}\left\{ Y_{n}\right\} is an 'autoregressive' stationary sequence in that it satisfies Yn=αYn−1+ZnY_{n} = \alpha Y_{n-1}+Z_{n}, −∞<n<∞-\infty < n < \infty, for some real α\alpha satisfying ∣α∣<1\left|\alpha \right| < 1. Show that YY has autocovariance function c(m)=α∣m∣/(1−α2)c(m) = \alpha^{\left|m\right|} /\left(1-\alpha^{2}\right).

?
Problem 8.10.3

Let {Xn}\left\{ X_{n}\right\} be a sequence of independent identically distributed Bernoulli variables, each taking values 0 and 1 with probabilities 1−p1-p and pp respectively. Find the mass function of the renewal process N(t)N(t) with interarrival times {Xn}\left\{ X_{n}\right\}.

?
Problem 8.10.4

Customers arrive in a shop in the manner of a Poisson process with parameter λ\lambda. There are infinitely many servers, and each service time is exponentially distributed with parameter μ\mu. Show that the number Q(t)Q(t) of waiting customers at time tt constitutes a birth-death process. Find its stationary distribution.

?
Problem 8.10.5

Let X(t)=Ycos⁡(θt)+Zsin⁡(θt)X(t) = Y \cos (\theta t)+Z \sin (\theta t) where YY and ZZ are independent N(0,1)N(0,1) random variables, and let X~(t)=Rcos⁡(θt+Ψ)\tilde{X}(t) = R \cos (\theta t+\Psi ) where RR and Ψ\Psi are independent. Find distributions for RR and Ψ\Psi such that the processes XX and X~\widetilde{X} have the same fdds.

?
Problem 8.10.6

Customers arrive at the entrance to a queueing system at the instants of an non-homogeneous Poisson process with rate function λ(t)\lambda (t). Their subsequent service histories are independent of each other, and a customer arriving at time ss is in state AA at time s+ts+t with probability p(s,t)p(s, t). Show that the number of customers in state AA at time tt is Poisson with parameter ∫−∞tλ(u)p(u,t−u)du\int_{-\infty }^{t} \lambda (u) p(u, t-u) d u.

?
Problem 8.10.7

An insurance company receives premiums (net of costs) at unit rate per unit time, and the claims of size X1,X2,…X_{1}, X_{2}, \ldots are independent random variables with common distribution function FF, arriving at the instants of a Poisson process of rate 1. Show that the probability r(y)r(y) of ruin (that is, assets becoming negative) for an initial capital yy satisfies

dr(y)dy=λr(y)−λP(X1>y)−λ∫0yr(y−x)dF(x) \frac{d r(y)}{d y} = \lambda r(y)-\lambda \mathbb {P}\left(X_{1} > y\right)-\lambda \int _{0}^{y} r(y-x) d F(x)
?
Problem 8.10.8

Let XX be fBMH\mathrm{fBM}_{H} with Var⁡[X(1)]=1\operatorname {Var}\left[X(1)\right] = 1, where H∈(0,1)H \in (0,1). Show that:

?
(a)

X(t)−X(s)X(t)-X(s) has the N(0,∣t−s∣2H)N\left(0,|t-s|^{2 H}\right) distribution,

(b)

for r>0r > 0,

E[∣X(t)−X(s)∣r]=C∣t−s∣rH \mathbb {E}\left[|X(t)-X(s)|^{r}\right] = C|t-s|^{r H}

where CC depends on rr.

(c)

XX is a Wiener process if and only if H=12H = \frac{1}{2}.

Problem 8.10.9

Holtsmark distribution for stellar gravity, Problem (6.15.57) revisited. Suppose stars are distributed in R3\mathbb {R}^{3} in the manner of a Poisson process with intensity λ\lambda, and let GλG_{\lambda } be the xx-coordinate of their aggregated gravitational force exerted on a unit mass at the origin.

?
(a)

Show that Gλ+μ=DGλ+GμG_{\lambda +\mu } \stackrel{\mathrm{D}}{ = } G_{\lambda }+G_{\mu }, the sum of two independent variables.

(b)

Show by the inverse square law that Gλ=Dλ2/3G1G_{\lambda } \stackrel{\mathrm{D}}{ = } \lambda^{2 / 3} G_{1}.

(c)

Comment on the relevance of the above in studying gravitational attraction in three dimensions.

Problem 8.10.10

In a Prague teashop (U Myšáka) before the Velvet Revolution of 1989, customers queue at the entrance for a blank bill. In the shop there are separate counters for coffee, sweetcakes, pretzels, milk, drinks, and ice cream, and queues form at each of these. At each service point the customers' bills are marked appropriately. There is a restricted number NN of seats, and departing customers have to queue in order to pay their bills. If interarrival times and service times are exponentially distributed and the process is in equilibrium, find how much longer a greedy customer must wait if he insists on sitting down. Answers on a postcard to the authors, please.

?