11.
(b) Follows from the hint about using the lack of memory property and the fact that εi, the minimum of j-(i-1) independent exponentials with rate λ, is exponential with rate (j-i-1)λ.
(c) From parts (a) and (b)
P{T1+⋯+Tj⩽t}=Pmax1⩽i⩽jXi⩽t=(1-e-λt)j
(d) With all probabilities conditional on X(0)=1,
P1j(t)=P{X(t)=j}=P{X(t)⩾j}-P{X(t)⩾j+1}=P{T1+⋯+Tj⩽t}-P{T1+⋯+Tj+1⩽t}
(e) The sum of i independent geometrics, each having parameter p=e-λt, is a negative binomial with parameters i,p. The result follows since starting with an initial population of i is equivalent to having i 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
Since ∑02Pi=1, we get
P2=1+μ2λα+1-ααμ2μ1-1=λαμ1λαμ1+μ1μ2+λ(1-α)μ2
where P2 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
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
The equations are easily solved and the proportion of time machine 2 is down is P2+P3.
24. We will let the state be the number of taxis waiting. Then, we get a birth and death process with λn=1,μn=2. This is an M/M/1. Therefore:
(a) Average number of taxis waiting =1μ-λ=12-1=1.
(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). The proportion of such arrivals is therefore
2(1-P0)2=1-P0=1-1-λμ=λμ=12
28. Let Pijx,vix denote the parameters of the X(t) and Pijy,viy of the Y(t) process; and let the limiting probabilities be Pix,Piy, respectively. By independence we have that for the Markov chain {X(t),Y(t)} its parameters are
v(i,l)=vix+vly,P(i,l)(j,l)=vixvix+vlyPijx,P(i,l)(i,k)=vlyvix+vlyPlky,
and
limt→∞P{(X(t),Y(t))=(i,j)}=PixPjy
Hence, we need to show that
PixPlyvixPijx=PjxPlyvjxPjix
(That is, the rate from (i,l) to (j,l) equals the rate from (j,l) to (i,l).) But this follows from the fact that the rate from i to j in X(t) equals the rate from j to i; that is,
PixvixPijx=PjxvjxPjix
The analysis is similar in looking at pairs (i,l) and (i,k).
33. Suppose first that the waiting room is of infinite size. Let Xi(t) denote the number of customers at server i,i=1,2. Then since each of the M/M/1 processes {X1(t)} is time reversible, it follows from Exercise 28 that the vector process {(X1(t),(X(t)),t⩾0} is a time reversible Markov chain. Now the process of interest is just the truncation of this vector process to the set of states A where
A={(0,m):m⩽4}∪{(n,0):n⩽4}∪{(n,m):nm>0,n+m⩽5}
Hence, the probability that there are n with server 1 and m with server 2 is
Pn,m=kλ1μ1n1-λ1μ1λ2μ2m1-λ2μ2=Cλ1μ1nλ2μ2m,(n,m)∈A
The constant C is determined from
∑Pn,m=1
where the sum is over all (n,m) in A.
37.
(a) The state is (n1,…,nk) if there are ni type i patients in the hospital, for all i=1,…,k.
(b) It is a M/M/∞ birth and death process, and thus time reversible.
(c) Because Ni(t),t⩾0 are independent processes for i=1,…,k, the vector process is a time reversible continuous time Markov chain.
(d)
P(n1,…,nk)=∏i=1ke-λi/μi(λi/μi)ni/ni!
(e) As a truncation of a time reversible continuous time Markov chain, it has stationary probabilites
P(n1,…,nk)=K∏i=1ke-λi/μi(λi/μi)ni/ni!,(n1,…,nk)∈A
where A={(n1,…,nk):∑i=1kniwi⩽C}, and K is such that
K∑(n1,…,nk)∈A∏i=1ke-λi/μi(λi/μi)ni/ni!=1
(f) With ri equal to the rate at which type i patients are admitted,
ri=∑(n1,…,ni+1,…,nk)∈AλiP(n1,…,nk)
(g) ∑i=1kri/;∑i=1kλi
38.
(a) The state is the set of idle servers.
(b) For i∈S,j∉S, the infinitesimal rates of the chain are
qS,S-i=λ/∣S∣,qS,S+j=μj
where ∣S∣ is the number of elements in S. The time reversibility equations are
P(S)λ/∣S∣=P(S-i)μi
which has a solution
P(S)=P0∣S∣!∏k∈S(μk/λ)
where P0, the probability there are no idle servers, is found from
P0[1+∑S∣S∣!∏k∈S(μk/λ)]=1
where the preceding sum is over all nonempty subsets of {1,…,n}.
39.
(a) The state is (i1,…,ik) if {i1,…,ik} is the set of idle servers, with i1 having been idle the longest, i2 the second longest, and so on.
(b) and
(c) For j∉{i1,…,ik}, the infinitesimal rates of the chain are
q(i1,…,ik),(i1,…,ik-1)=λ,q(i1,…,ik),(i1,…,ik,j)=μj
The time reversibility equations are
P(i1,…,ik)λ=P(i1,…,ik-1)μik
giving the solution
P(i1,…,ik)=μi1⋯μikλkP(0)
where P(0) is the probability there are no idle servers.
40. The time reversible equations are
P(i)vin-1=P(j)vjn-1
yielding the solution
P(j)=1/vj∑i=1n1/vi
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/1 queue is a Poisson process it follows that the number of customers with server 2 is the stationary probability of an M/M/1 system.
43. We make the conjecture that the reverse chain is a system of same type, except that the Poisson arrivals at rate λ arrive at server 3, then go to server 2, then to server 1, and then depart the system. Let ek be the 3-vector with 1 in position k and 0 elsewhere. With the state being i=(i1,i2,i3) when that there are ij customers at server j for j=1,2,3, 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
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
The conjecture is correct if we can find probabilities P(i,j,k) 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
satisfy.
49.
(a) The matrix P∗ can be written as
P∗=I+R/v
and so Pij∗n can be obtained by taking the i,j element of (I+R/v)n, which gives the result when v=n/t.
(b) Uniformization shows that Pij(t)=E[Pij∗N], where N is independent of the Markov chain with transition probabilities Pij∗ and is Poisson distributed with mean vt. Since a Poisson random variable with mean vt has standard deviation (vt)1/2, it follows that for large values of vt it should be near vt. (For instance, a Poisson random variable with mean 106 has standard deviation 103 and thus will, with high probability, be within 3000 of 106.) Hence, since for fixed i and j,Pij∗m should not vary much for values of m about vt where vt is large, it follows that, for large vt,
E[Pij∗N]≈Pij∗nwheren=vt