4
$\begingroup$

I am almost finished reading The Goal, which emphasizes the importance of managing bottlenecks in production. I'm more into service networks than production, but there is a lot to carry over from one to the other. One the lessons is that meeting capacity with demand on average is not enough due to a combination of uncontrolled/predicted variation in demand/service and dependency between events. This brings us to queuing, as an example, which frequently occurs in the services I am interested in modelling.

I have simulated a M/M/1 queue in which $\mu = \lambda$. It is well-known that this is an unstable system; the wait list will tend to get longer and so will the wait times. But at what rate? I don't know an analytic result around this.

Here is an example output where $\lambda = \mu = 10$:

enter image description here

Here is an example output where $\lambda = \mu = 1$:

enter image description here

Prima facie it looks something like the growth is logarithmic (first guess) or like a square root (second guess). I am not sure if this hold true at every quantile, or just specific statistics like the maximum or average. It also appears (based on the two cases that were run) that it dependends on the values of $\lambda$ (or equivalent $\mu$ since they are set to be equal).

$\endgroup$
2
  • 1
    $\begingroup$ Note: on the Wikipedia page, it claims that for $\lambda / \mu \approx 1$, the process is approximately $|(\lambda - \mu)t + \sqrt{\lambda + \mu}W_t|$ for Brownian $W_t$; this of course is just $\sqrt{(\lambda + \mu) t}$ in expectation when $\lambda = \mu$ in your case. This seems to fit your simulations very well. $\endgroup$
    – daisies
    Commented Feb 12 at 20:21
  • $\begingroup$ @daisies That looks more like an answer than a comment to me! $\endgroup$
    – Galen
    Commented Feb 12 at 20:31

1 Answer 1

2
$\begingroup$

On the Wikipedia page, it claims that for $\lambda / \mu \approx 1$, the process is approximately $|(\lambda - \mu)t + \sqrt{\lambda + \mu}W_t|$ for Brownian $W_t$; this of course is $\propto \sqrt{(\lambda + \mu) t}$ in expectation when $\lambda = \mu$ in your case. This seems to fit your simulations very well.

More precisely, you can check this (different) Wikipedia link for details (and recall that $|W_t| = \sup_{t} W_t$ in distribution).

$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.