EEL6507sp09L38
From BoykinWiki
Contents |
EEL6507 Spring 2009, Lecture 38, Monday 2009/04/20 (Notes created by Xiaoyuan Li)
- Conditional Distribution on Queue size, G/M/m
Conditional Distribution on Queue size
From last time,
![\begin{align}
P_{i,j}& =Prob[q_{n+1}'=j|q_n'=i] \\
& = Prob[i+1-j\ departs|q_n'=i] \\
& = \int_{0}^{\infty}\binom{i+1}{j}(1-e^{-\mu t})^{i+1-j}e^{-\mu tj}a(t)\,dt, for j\le i+1 \le m ,
\end{align}](/wiki/images/math/9/9/5/995f77d9f0322688b4ee27d361ef20dd.png)
Refer to the Regions of pij
- All busy:
![Prob[k\ customer\ served|time=t, m\ busy] = \frac{(m\mu t)^k}{k!}e^{-m\mu t}](/wiki/images/math/3/3/1/33119f70cf06dcae57ae2da9664fc463.png)
When



- Transition Period: All Busy
Some Busy
![P_{i,j}= \int_{0}^{\infty}\binom{m}{j}e^{-j\mu t} \times \left[\int_{0}^{\infty}\frac{(m\mu t)^{i-m}}{(i-m)!}(e^{-\mu y}-e^{-\mu t})^{m-j}m\mu\,dy\right]a(t)dt](/wiki/images/math/f/0/7/f0716ca977e1fac6f5303da0f7d95009.png)
We want to get
![Prob[k\ customer\ served|time=t, m\ busy] = \frac{(m\mu t)^k}{k!}e^{-m\mu t}](/wiki/images/math/3/3/1/33119f70cf06dcae57ae2da9664fc463.png)
When
![\gamma_k\ =\ Prob[\ arrival\ finds\ k\ in\ system]](/wiki/images/math/4/b/8/4b8016761fd6c2328c5becccb6f8dcdf.png)


Given arrival finds m, there is queue
Conditional Queue Length (Arrival Queues)

Refer to the [Chain Diagram ]
number of arrivals in (0,t) that finds k in the system
![\begin{align}
\sigma_k &= E[number\ of\ visits\ to\ E_{k+1}\ between\ visits\ to\ E_{k}]\\
&= \beta_0\sum n\gamma^{n-1}(1-\gamma)\\
&=\frac{\beta_0}{1-\gamma}\\
&=\sigma_k
\end{align}](/wiki/images/math/e/9/6/e96dc5f0bbd7d71d51899dd7b733bf72.png)
Observation
![\begin{align}
\gamma&=Prob[leave\ E_{k+1}\ and\ return\ without\ passing\ E_j\ where\ j\le k]\\
&=Prob[leave\ E_{k+1}\ and\ return\ without\ passing\ E_k]
\end{align}](/wiki/images/math/a/8/d/a8d9396c51afbccf9c844afb060c118a.png)

![Prob[n-1 \ visits\ to\ E_{k+1}\ without\ E_{k}, 1\ with\ E_k ]=\gamma^{n-1}(1-\gamma)\beta_0](/wiki/images/math/5/f/a/5fa16e061484f0e4f38e20eeb9cbdc3a.png)



![\vec \gamma=[\gamma_0,\gamma_1,...\gamma_{m},\sigma \gamma_m,\sigma^2 \gamma_m,... ]](/wiki/images/math/9/1/4/914b6fdf17ed6a1247fa077a11525ed8.png)

