image

11. 

(b) Follows from the hint about using the lack of memory property and the fact that εiimage, the minimum of j-(i-1)image independent exponentials with rate λimage, is exponential with rate (j-i-1)λimage.

(c) From parts (a) and (b)

P{T1++Tjt}=Pmax1ijXit=(1-e-λt)j

image

(d) With all probabilities conditional on X(0)=1image,

P1j(t)=P{X(t)=j}=P{X(t)j}-P{X(t)j+1}=P{T1++Tjt}-P{T1++Tj+1t}

image

(e) The sum of iimage independent geometrics, each having parameter p=e-λtimage, is a negative binomial with parameters i,pimage. The result follows since starting with an initial population of iimage is equivalent to having iimage independent Yule processes, each starting with a single individual.

16. Let the state be

2: an acceptable molecule is attached

0: no molecule attached

1: an unacceptable molecule is attached.

Then, this is a birth and death process with balance equations

μ1P1=λ(1-α)P0μ2P2=λαP0

image

Since 02Pi=1image, we get

P2=1+μ2λα+1-ααμ2μ1-1=λαμ1λαμ1+μ1μ2+λ(1-α)μ2

image

where P2image is the percentage of time the site is occupied by an acceptable molecule. The percentage of time the site is occupied by an unacceptable molecule is

P1=1-ααμ2μ1P1=λ(1-α)μ2λαμ1+μ1μ2+λ(1-α)μ2

image

19. There are four states. Let state 0 mean that no machines are down, state 1 that machine 1 is down and 2 is up, state 2 that machine 1 is up and 2 is down, and state 3 that both machines are down. The balance equations are as follows:

(λ1+λ2)P0=μ1P1+μ2P2(μ1+λ2)P1=λ1P0(λ1+μ2)P2=λ2P0+μ1P3μ1P3=λ2P1+λ1P2P0+P1+P2+P3=1

image

The equations are easily solved and the proportion of time machine 2 is down is P2+P3image.

24. We will let the state be the number of taxis waiting. Then, we get a birth and death process with λn=1,μn=2image. This is an M/M/1image. Therefore:

(a) Average number of taxis waiting =1μ-λ=12-1=1image.

(b) The proportion of arriving customers that gets taxis is the proportion of arriving customers that find at least one taxi waiting. The rate of arrival of such customers is 2(1-P0)image. The proportion of such arrivals is therefore

2(1-P0)2=1-P0=1-1-λμ=λμ=12

image

28. Let Pijx,viximage denote the parameters of the X(t)image and Pijy,viyimage of the Y(t)image process; and let the limiting probabilities be Pix,Piyimage, respectively. By independence we have that for the Markov chain {X(t),Y(t)}image its parameters are

v(i,l)=vix+vly,P(i,l)(j,l)=vixvix+vlyPijx,P(i,l)(i,k)=vlyvix+vlyPlky,

image

and

limtP{(X(t),Y(t))=(i,j)}=PixPjy

image

Hence, we need to show that

PixPlyvixPijx=PjxPlyvjxPjix

image

(That is, the rate from (i,l)image to (j,l)image equals the rate from (j,l)image to (i,l)image.) But this follows from the fact that the rate from iimage to jimage in X(t)image equals the rate from jimage to iimage; that is,

PixvixPijx=PjxvjxPjix

image

The analysis is similar in looking at pairs (i,l)image and (i,k)image.

33. Suppose first that the waiting room is of infinite size. Let Xi(t)image denote the number of customers at server i,i=1,2image. Then since each of the M/M/1image processes {X1(t)}image is time reversible, it follows from Exercise 28 that the vector process {(X1(t),(X(t)),t0}image is a time reversible Markov chain. Now the process of interest is just the truncation of this vector process to the set of states Aimage where

A={(0,m):m4}{(n,0):n4}{(n,m):nm>0,n+m5}

image

Hence, the probability that there are nimage with server 1 and mimage with server 2 is

Pn,m=kλ1μ1n1-λ1μ1λ2μ2m1-λ2μ2=Cλ1μ1nλ2μ2m,(n,m)A

image

The constant Cimage is determined from

Pn,m=1

image

where the sum is over all (n,m)image in Aimage.

37. 

(a) The state is (n1,,nk)image if there are niimage type iimage patients in the hospital, for all i=1,,kimage.

(b) It is a M/M/image birth and death process, and thus time reversible.

(c) Because Ni(t),t0image are independent processes for i=1,,kimage, the vector process is a time reversible continuous time Markov chain.

(d) 

P(n1,,nk)=i=1ke-λi/μi(λi/μi)ni/ni!

image

(e) As a truncation of a time reversible continuous time Markov chain, it has stationary probabilites

P(n1,,nk)=Ki=1ke-λi/μi(λi/μi)ni/ni!,(n1,,nk)A

image

where A={(n1,,nk):i=1kniwiC}image, and Kimage is such that

K(n1,,nk)Ai=1ke-λi/μi(λi/μi)ni/ni!=1

image

(f) With riimage equal to the rate at which type iimage patients are admitted,

ri=(n1,,ni+1,,nk)AλiP(n1,,nk)

image

(g) i=1kri/;i=1kλiimage

38. 

(a) The state is the set of idle servers.

(b) For iS,jSimage, the infinitesimal rates of the chain are

qS,S-i=λ/S,qS,S+j=μj

image

where Simage is the number of elements in Simage. The time reversibility equations are

P(S)λ/S=P(S-i)μi

image

which has a solution

P(S)=P0S!kS(μk/λ)

image

where P0image, the probability there are no idle servers, is found from

P0[1+SS!kS(μk/λ)]=1

image

where the preceding sum is over all nonempty subsets of {1,,n}image.

39. 

(a) The state is (i1,,ik)image if {i1,,ik}image is the set of idle servers, with i1image having been idle the longest, i2image the second longest, and so on.

(b) and

(c) For j{i1,,ik}image, the infinitesimal rates of the chain are

q(i1,,ik),(i1,,ik-1)=λ,q(i1,,ik),(i1,,ik,j)=μj

image

The time reversibility equations are

P(i1,,ik)λ=P(i1,,ik-1)μik

image

giving the solution

P(i1,,ik)=μi1μikλkP(0)

image

where P(0)image is the probability there are no idle servers.

40. The time reversible equations are

P(i)vin-1=P(j)vjn-1

image

yielding the solution

P(j)=1/vji=1n1/vi

image

Hence, the chain is time reversible with long run proportions given by the preceding.

41. Show in Example 6.22 that the limiting probabilities satisfy Equations (6.33), (6.34), and (6.35).

42. Because the stationary departure process from an M/M/1image queue is a Poisson process it follows that the number of customers with server 2image is the stationary probability of an M/M/1image system.

43. We make the conjecture that the reverse chain is a system of same type, except that the Poisson arrivals at rate λimage arrive at server 3image, then go to server 2image, then to server 1image, and then depart the system. Let ekimage be the 3-vector with 1image in position kimage and 0image elsewhere. With the state being i=(i1,i2,i3)image when that there are ijimage customers at server jimage for j=1,2,3image, the instantaneous transition rates of the chain are

q(i,j,k),(i+1,j,k)=λq(i,j,k),(i-1,j+1,k)=μ1,i>0q(i,j,k),(i,j-1,k+1)=μ2,j>0q(i,j,k),(i,j,k-1)=μ3,k>0

image

whereas the conjectured instantaneous rates for the reversed chain are

q(i,j,k),(i,j,k+1)=λq(i,j,k),(i,j+1,k-1)=μ3,k>0q(i,j,k),(i+1,j-1,k)=μ2,j>0q(i,j,k),(i-1,j,k)=μ1,i>0

image

The conjecture is correct if we can find probabilities P(i,j,k)image that satisfy the reverse time equations when the preceding are the instantaneous rates for the reversed chain, and it is easy to check that

P(i,j,k)=Kλμ1iλμ2jλμ3k

image

satisfy.

49. 

(a) The matrix Pimage can be written as

P=I+R/v

image

and so Pijnimage can be obtained by taking the i,jimage element of (I+R/v)nimage, which gives the result when v=n/timage.

(b) Uniformization shows that Pij(t)=E[PijN]image, where Nimage is independent of the Markov chain with transition probabilities Pijimage and is Poisson distributed with mean vtimage. Since a Poisson random variable with mean vtimage has standard deviation (vt)1/2image, it follows that for large values of vtimage it should be near vtimage. (For instance, a Poisson random variable with mean 106image has standard deviation 103image and thus will, with high probability, be within 3000 of 106image.) Hence, since for fixed iimage and j,Pijmimage should not vary much for values of mimage about vtimage where vtimage is large, it follows that, for large vtimage,

E[PijN]Pijnwheren=vt

image

Chapter 7

3. 

(a) 

P(N(t)=nSn=y)=1-F(t-y),ifyt0,ify>t

image

(b) 

P(N(t)=n)=0P(N(t)=nSn=y)fSn(y)dy=0te-λ(t-y)λe-λy(λy)n-1/(n-1)!dy=e-λtλn(n-1)!0tyn-1dy=e-λt(λt)nn!

image

6. 

(a) Consider a Poisson process having rate λimage and say that an event of the renewal process occurs whenever one of the events numbered r,2r,3r,image of the Poisson process occurs. Then

P{N(t)n}=P{nror more Poisson events byt}=i=nre-λt(λt)i/i!

image

(b) 

E[N(t)]=n=1P{N(t)n}=n=1i=nre-λt(λt)i/i!=i=rn=1[i/r]e-λt(λt)i/i!=i=r[i/r]e-λt(λt)i/i!

image

8. 

(a) The number of replaced machines by time timage constitutes a renewal process. The time between replacements equals Timage, if the lifetime of new machine is T;ximage, if the lifetime of new machine is x,x<Timage. Hence,

E[time between replacements]=0Txf(x)dx+T[1-F(T)]

image

and the result follows by Proposition 7.1.

(b) The number of machines that have failed in use by time timage constitutes a renewal process. The mean time between in-use failures, E[F]image, can be calculated by conditioning on the lifetime of the initial machine as E[F]=E[E[Fimage lifetime of initial machine]]. Now

E[Flifetime of machine isx]=x,ifxTT+E[F],ifx>T

image

Hence,

E[F]=0Txf(x)dx+(T+E[F])[1-F(T)]

image

or

E[F]=0Txf(x)dx+T[1-F(T)]F(T)

image

and the result follows from Proposition 7.1.

13. With Wiimage equal to your winnings in game i,i1image, and Nimage the number of games played, Wald’s equation yields

E[X]=E[N]E[W]=0

image

With pi=P(X=i),p1=1/2+1/8=5/8,p-1=1/4,p-3=1/8image, verifying that E[X]=0image.

18. We can imagine that a renewal corresponds to a machine failure, and each time a new machine is put in use its life distribution will be exponential with rate μ1image with probability pimage, and exponential with rate μ2image otherwise. Hence, if our state is the index of the exponential life distribution of the machine presently in use, then this is a two-state continuous-time Markov chain with intensity rates

q1,2=μ1(1-p),q2,1=μ2p

image

Hence,

P11(t)=μ1(1-p)μ1(1-p)+μ2pexp{-[μ1(1-p)+μ2p]t}+μ2pμ1(1-p)+μ2p

image

with similar expressions for the other transition probabilities (P12(t)=1-P11(t)image, and P22(t)image is the same with μ2pimage and μ1(1-p)image switching places). Conditioning on the initial machine now gives

E[Y(t)]=pE[Y(t)X(0)=1]+(1-p)E[Y(t)X(0)=2]=pP11(t)μ1+P12(t)μ2+(1-p)P21(t)μ1+P22(t)μ2

image

Finally, we can obtain m(t)image from

μ[m(t)+1]=t+E[Y(t)]

image

where

μ=p/μ1+(1-p)/μ2

image

is the mean interarrival time.

22 

(a) Let Ximage denote the length of time that J keeps a car. Let Iimage equal 1image if there is a breakdown by time Timage and equal 0image otherwise. Then

E[X]=E[XI=1](1-e-λT)+E[XI=0]e-λT=T+1μ(1-e-λT)+T+1λe-λT=T+1-e-λTμ+e-λTλ

image

1/E[X]image is the rate that J buys a new car.

(b) Let Wimage equal to the total cost involved with purchasing a car. Then, with Yimage equal to the time of the first breakdown

E[W]=0E[WY=y]λe-λydy=C+0Tr(1+μ(T-y)+1)λe-λydy+Trλe-λydy=C+r(2-e-λT)+r0Tμ(T-y)λe-λydy

image

J’s long run average cost is E[W]/E[X]image.

23. 

(a) Say that a new cycle begins each time A wins a point. With Nimage equal to the number of points in a cycle

E[N]=1+qa/pb

image

where the preceding used that, starting with Bimage serving, the number of points played until A wins a point is geometric with parameter pbimage. Hence, by renewal reward, the proportion of points won by A is 1/E[N]=pbpb+qaimage.

(b) (pa+pb)/2image

(c) pbpb+qa>(pa+pb)/2image is equivalent to paqa>pbqbimage.

30. 

A(t)t=t-SN(t)t=1-SN(t)t=1-SN(t)N(t)N(t)t

image

The result follows since SN(t)/N(t)μimage (by the strong law of large numbers) and N(t)/t1/μimage.

35. 

(a) We can view this as an M/G/image system where a satellite launching corresponds to an arrival and Fimage is the service distribution. Hence,

P{X(t)=k}=e-λ(t)[λ(t)]k/k!

image

where λ(t)=λ0t(1-F(s))dsimage.

(b) By viewing the system as an alternating renewal process that is on when there is at least one satellite orbiting, we obtain

limP{X(t)=0}=1/λ1/λ+E[T]

image

where Timage, the on time in a cycle, is the quantity of interest. From part (a)

limP{X(t)=0}=e-λμ

image

where μ=0(1-F(s))dsimage is the mean time that a satellite orbits. Hence,

e-λμ=1/λ1/λ+E[T]

image

so

E[T]=1-e-λμλe-λμ