Опубликован: 20.04.2011 | Уровень: для всех | Доступ: платный
Лекция 14:

Прикладная теория организации очередей

Распределение времени ожидания: M/D/n, FCFS

Распределение времени ожидания определяется распределением Кроммелина:

p\{W \le t \}=1-\sum_{i=0}^{n-1} \sum_{k=0}^i p(k)*\sum_{j=1}^{\infty} \frac{\{A(j- \tau)\}^{(T+j+1)n-1-i}}{\{(T+j+1)n-1-i\}!}, ( 13.56)

где A - предложенная нагрузка и:

t=T*h+ \tau, 0 \le \tau < h. ( 13.57)

В компактной форме по аналогии с (13.50):

p\{W \le t\}=\sum_{i=0}^{n-1} \sum_{k=0}^i p(k)\sum_{j=0}^T \frac{\{A(j-t)\}^{j*n+n-1-i}}{\{j*n+n-1-i)!}*e^{-A(j-t)}. ( 13.58)

Для составных значений времени ожидания t мы имеем:

p\{W \le t \}=\sum_{j=0}^{n(t+1)-1}p(j) ( 13.59)

Для несоставных времен ожидания при t = T+ \tau, T - целое число, 0 \le \tau < 1 можно выразить распределение времени ожидания в элементах составных времен ожидания, поскольку для M/D/1:

p\{W \le t\}=p\{W \le T+ \tau\}=e^{ \lambda \tau}\sum_{j=0}^k \left \{ \frac{(- \lambda \tau)^j}{j!}*\sum_{i=0}^{k-j}p(i) \right \}, ( 13.60)

где n(T+1)-1 и p(i) - вероятность состояния (13.55).

Точное среднее время ожидания всех клиентов получить W трудно. Приближенное значение получено Молина (Molina)

W \approx \frac{n}{n+1}*E_{2,n}(A)*\frac{n}{n-A}*\frac{1-(\frac An)^{n+1}}{1-(\frac An)^n}. ( 13.61)

Для любой системы с бесконечной очередью мы имеем (3.20):

w=\frac WD

где:

D=1-\sum_{j=0}^{n-1}p(j)

k -Эрланговский процесс поступления вызовов: Ek /D/r

Рассмотрим систему организации очереди к системе, в которой: n = r * k обслуживающих приборов ( r, k - целые числа), процесс поступления вызовов - общий, время обслуживания - постоянное и дисциплина организации очереди - FCFS. Клиенты, прибывающие в течение, когда обслуживающие приборы свободны, выбирают их в циклическом порядке:

1,2, \dots ,п- 1,п,1,2, \dots

Тогда некоторый обслуживающий прибор обслужит только каждого n -ого клиента, так как клиенты из-за постоянного времени обслуживания освобождают обслуживающие приборы в том же порядке, в каком они занимают обслуживающие приборы. Ни один клиент не может догнать другого клиента в очереди.

Пусть обслуживающие приборы разбиты на группы, тогда

x,x+k,x+2*2, \dots, x+(r-1)*k, \quad 0 < x \le k. ( 13.62)

прибор в группе обслужит только каждого k -ого клиента. Если мы рассматриваем обслуживающие приборы в (13.62) как одну группу, то они эквивалентны системе организации очереди GI^{k*}/D/r, где процесс поступления вызовов GI^{k*} есть свертка распределения времени прибытия отдельных k времен.

То же самое относится к другим k-1 системам. Нагрузка в этих k системах является взаимно коррелированной, но если рассматривать только одну систему одновременно, то она - GI^{k*}/D/n, FCFS -система организации очереди.

Предположение о циклическом поиске обслуживающих приборов в пределах отдельных систем обязательно (13.62). Вероятности состояния, средние времена ожидания и т.д. не зависят от дисциплины организации очереди, которая играет важную роль только для распределения времени ожидания.

Если мы предположим, что GI -процесс поступления вызовов - Пуассоновский процесс, то GI^{k*} становится k-Эрланговским процессом поступления вызовов. Таким образом, мы получили, что следующие системы эквивалентны относительно распределения времени ожидания:

M/D/r-k, FCFS \equiv  E_k/D/r, FCFS

Поэтому система Ek/D/r может работать с таблицами для M/D/n.

Вообще мы знаем, что нагрузка на один обслуживающий прибор и среднее времени ожидания уменьшается, когда число обслуживающих приборов увеличивается. По той же причине среднее время ожидания уменьшается, когда процесс поступления вызовов становится более регулярным. Это замечено непосредственно из вышеупомянутого разложения, где процесс поступления вызовов для E_k /D/r становится более регулярным при увеличении k ( r константа). Для = 0,9 Эрл. на один обслуживающий прибор ( L = средняя длина очереди) находим:

E_4/E_1/2:\qquad L = 4.5174 ,\\
Е_4/Е_2/2: \qquad L = 2.6607 ,\\
Е_4/Е_3/2: \qquad L = 2.0493 ,\\
E_4/D/2:\qquad L = 0.8100

Система с конечной очередью: M/D/1/k

В реальных системах мы всегда имеем конечную очередь. В компьютерных системах размер памяти конечен, и в системах ATM существуют конечные буфера. То же самое происходит для мест ожидания в гибких производственных системах ( FMS - Flexible Manufacturing Systems ).

Как показано в секции 13.3.4, вероятности состояния p_k(i) для конечной буферной системы получены из вероятностей состояний p(i) бесконечной буферной системы, используя (13.10) и (13.11).

Составные времена ожидания получены из вероятностей состояний и несоставных времен ожидания - из составных времен ожидания, как показано выше (секция 13.5.4).

Для бесконечной буферной системы вероятности состояния существуют только, когда предложенная нагрузка меньше, чем емкость (A < n) . Но для конечной буферной системы вероятности состояния существуют и для A > n, однако мы не можем получить их вышеупомянутым методом.

В M/D/1/k конечные буферные вероятности состояния p_k (i) могут быть получены для любой предложенной нагрузки следующим способом.

В системе с одним сервером и (k- 1) местами ожидания мы имеем (k + 1) состояний (0, 1 \dots , k) . Уравнения равновесия для вероятностей состояния p_k(i), i = 0, 1 \dots , k-2, составляя k-1 линейных уравнений между состоянием \{p_k(0), p_k(1), \dots , p_k(k-1)\} , могут быть установлены с помощью уравнений состояний Фрея. Но невозможно записать простые уравнения времени для состояния k- 1 и k. Однако первые (k- 1) уравнения (13.41) вместе с требованием нормализации:

\sum_{j=0}^kp_k(j)=1 ( 13.63)

и фактом, что предложенная нагрузка равняется обслуженной нагрузке плюс отклоненная нагрузка ( свойство PASTA ) дают

A=1-p_k(0)+A*p_k(k) ( 13.64)

в результате (k+1) независимых линейных уравнений, которые легко решить.

Два подхода дают, конечно, один и тот же результат. Первый метод справедлив только для A<1, тогда как второй - для любой предложенной нагрузки.

Пример 13.5.2: "Дырявое ведро"

"Дырявое ведро" - механизм для управления процессами поступления ячеек (пакетов) от пользователя (источника) в ATM -системе. Механизм соответствует системе организации очереди с постоянным временем обслуживания (размер ячейки) и конечным буфером. Если процесс поступления вызовов - Пуассоновский процесс, то мы имеем систему M/D/1/k. Размер "вытекания" из дыры соответствует долгосрочной средней интенсивности поступления, тогда как размер "ведра" описывает допустимый избыток (число пакетов). Механизм работает как виртуальная система организации очереди, где ячейки или принимаются немедленно, или отклоняются согласно значению счетчика, который является составным значением функции нагрузки (рис.13.2). В контракте между пользователем и сетью соглашение заключается при заданном размере утечки и размере "ведра". На этом основании сеть способна гарантировать заданный уровень обслуживания.

Система организации очереди с одним обслуживающим прибором: GI/G/1

В секции 13.3 мы показали, что среднее время ожидания для всех клиентов в системе организации очереди M/G/1 определяется формулой Полячека-Хинчина:

W=\frac{A*s}{2(1-A)}* \varepsilon ( 13.65)

где \varepsilon - коэффициент формы распределения времени пребывания в системе.

Ранее мы анализировали следующие случаи.

M/M/1 (секция 12.2.4): \varepsilon = 2:

W=\frac{A*s}{(1-A)}, \quad Erlang \quad 1917. ( 13.66)

M/M/1 (секция 13.5.3): \varepsilon    1:

W=\frac{A*s}{2(1-A)}, \quad Erlang \quad 1909. ( 13.67)

Показано, что чем регулярнее распределение времени пребывания в системе, тем меньше становится нагрузка времени ожидания (для систем с потерями и с ограниченной доступностью этот эффект противоположен: чем больше коэффициент формы, тем меньше потери).

В системах с не пуассоновским поступлением моменты более высокого порядка будут также влиять на среднее время ожидания.

Общие результаты

Мы до настоящего времени принимали, что процесс поступления вызовов - Пуассоновский процесс. Для других процессов поступления вызовов редко можно найти точное выражение среднего времени ожидания, кроме случаев, где времена пребывания в системе дают основания предлагать, что или процесс поступления вызовов, или процесс обслуживания - марковские. До настоящего времени нет никаких общих точных формул, например, для M/G/n.

Для GI/G/1 можно получить теоретические верхние пределы для среднего времени ожидания. Обозначая дисперсию интервалов поступления va и дисперсию распределения времени пребывания в системе - v_d, с помощью неравенства Кингмана (1961) мы получаем верхний предел для среднего времени ожидания:

GI/G/1^ W \le \frac{A*s}{2(1-A)}* \left \{ \frac{v_a+v_d}{s^2} \right \}. ( 13.68)

Эта формула показывает, что мы имеем дело со стохастическими вариациями времен ожидания.

Формула (13.68) дает теоретическую верхнюю границу. Реалистическая оценка фактического среднего времени ожидания получена приближением Марчала (Marchal, 1976)

W \approx \frac{A*s}{2(1-A)}* \left \{ \frac{v_a+v_d}{s^2} \right \}* \left \{ \frac{s^2+v_d}{a^2+v_d} \right \}. ( 13.69)

где a - средний интервал поступления ( A = s/a ). Это приближение получается масштабированием приближенного неравенства Кингмана, что согласуется с Формулой Полячека-Хинчина для случая M/G/1.

Как пример непуассоновского потока вызовов, мы проанализируем систему организации очереди GI/M/1, где распределение интервалов поступления - общее распределение, данное плотностью распределения f (t) . Времена обслуживания являются экспоненциально распределенными со скоростью \mu

Если систему рассматривают в произвольном момент времени, то вероятности состояния не будут отображаться марковским процессом, потому что вероятность прибытия будет зависеть от временного интервала, начиная с последнего прибытия, и свойство PASTA не будет выполняться.

Однако если систему рассматривать непосредственно до (или после) момента поступления, то она будет независима, так как интервалы поступления стохастическая независимы от времени пребывания в системе и являются экспоненциально распределенными. Моменты поступления - это точки равновесия (пункты регенерации, секция 5.2.2), и мы рассматриваем так называемую вложенную цепь Маркова.

Вероятность, что перед моментом поступления система находится в состоянии j, будем обозначать \pi(j) . Можно показать, что при статистическом равновесии будет получен следующий результат D.G. Kendall, 1953 [62] :

\pi(i)=(1-\alpha) \alpha^i, \quad i=0,1,2 \dots ( 13.70)

где \alpha - положительный реальный корень, удовлетворяющий уравнению:

\alpha=\int_0^{\infty}e^{-\mu (1- \alpha)t}f(t)dt. ( 13.71)

Устойчивые вероятности состояния могут быть получены при рассмотрении двух последовательных моментов поступления t_1 и t_2 (подобно уравнениям состояния Фрея, секция 13.5.5).

Так как процесс отклонения - Пуассоновский процесс с постоянной интенсивностью \mu, когда есть клиенты в системе, то вероятность p (j) , что для j клиентов обслуживание будет закончено между двумя моментами поступления, может быть выражена числом событий в Пуассоновском процессе в течение стохастического интервала (интервал поступления). Существуют следующие уравнения состояний:

\pi_{t_2}(0)=\sum_{j=0}^{\infty} \pi_{t_1}(j)* \left \{ 1- \sum_{i=0}^j p(i) \right \},\\
\pi_{t_2}(1)=\sum_{j=0}^{\infty} \pi_{t_1}(j)*p(j),\\
\vdots \qquad \vdots\\
\pi_{t_2}(i)=\sum_{j=0}^{\infty} \pi_{t_1}(j)*p(j-i+1) ( 13.72)

Условие нормализации определяется обычно так:

\sum_{i=0}^{\infty}\pi_{t_1}(i)=\sum_{j=0}^{\infty} \pi_{t_2}(j)=1 ( 13.73)

Можно показать, что для вышеупомянутого геометрического распределения это единственное решение данной системы уравнений (Kendall, 1953 [62] ).

В принципе, система организации очереди GI/M/n может быть рассмотрена с помощью того же способа, что и предыдущая. Тогда вероятность состояния p(j) становится более сложной, так как интенсивность освобождения зависит от числа занятых каналов.

Заметим, что \pi (i) - не вероятность нахождения системы в состоянии i в произвольный момент времени (математическое ожидание по времени), а вероятность нахождения системы в состоянии i непосредственно перед прибытием (вызвать математическое ожидание вызова).

Характеристики GI/M/1

Вероятность непосредственного обслуживания становится:

p\{непосредственно\}=\pi(0)=1- \alpha ( 13.74)

Соответствующая вероятность того, чтобы выполнение заявки было отсрочено, становится:

D=p\{ задержка \}=\alpha ( 13.75)

Среднее число занятых обслуживающих приборов в случайный момент времени (математическое ожидание по времени) равно обслуженной нагрузке (= предложенная нагрузка A < 1 ).

Среднее число ждущих клиентов, непосредственно перед прибытием клиента, может быть получено с помощью вероятности состояний:

L_1=\sum_{j=1}^{\infty}(1- \alpha) \alpha^i(i-1),\\
L_1=\frac{\alpha^2}{1-\alpha}. ( 13.76)

Среднее число клиентов в системе перед моментом поступления вызова:

L_2=\sum_{i=0}^{\infty}(1-\alpha) \alpha^i*i\\
=\frac{\alpha}{1-\alpha} ( 13.77)

Среднее время ожидания для всех клиентов тогда равно:

W=\frac{1}{\mu}*\frac{\alpha}{1-\alpha}. ( 13.78)

Средняя длина очереди, по всей оси (виртуальная длина очереди) поэтому равна (формула Литтла):

L=A*\frac{\alpha}{1-\alpha}. ( 13.79)

Среднее время ожидания для клиентов, которые стоят в очереди, равно:

w=\frac WD,\\
w=\frac{1}{\mu}*\frac{1}{1-\alpha}. ( 13.80)
Нияз Сабиров
Нияз Сабиров

Здравствуйте. А уточните, пожалуйста, по какой причине стоимость изменилась? Была стоимость в 1 рубль, стала в 9900 рублей.

Елена Сапегова
Елена Сапегова

для получения диплома нужно ли кроме теоретической части еще и практическую делать? написание самого диплома требуется?

Иван Бузмаков
Иван Бузмаков
Россия, Сарапул
Никита Сомов
Никита Сомов
Россия, Удмуртская республика