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

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

Пример 13.4.3: Система с программным управлением (см. пример 13.4.1)

Предположим, что компьютерная система в Примере 13.4.1 работает с дисциплиной "приоритетное возвращение к работе". Находим: Тип 1 - самый высокий приоритет:

W_1=\frac{\frac 12(0.1)^2}{1-0.1}+0=0.0056c.,\\
W_2=\frac{1.2850}{(1-0.1)(1-0.9}}+\frac{0.1}{1-0.1}*1.6=14.46c.

Тип приоритета 2:

W_2=\frac{\frac 12*0.5*(1.6)^2}{1-0.8}+0=6.40c,\\
W_1=\frac{1.2850}{(1-0.8)(1-0.9}}+\frac{0.8}{1-0.8}*0.1=64.65c.

Это показывает что, изменяя состав клиентов типа 1, мы можем обеспечить этим клиентам очень короткое время ожидания, не нарушая характеристик очереди для клиентов типа 2, но обратное преобразование не обеспечивает такого свойства.

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

M/M/n с приоритетным возвращением к работе

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

Системы организации очереди с постоянными временами занятия

В этой секции мы сосредотачиваем свое внимание на системе организации очереди M/D/n, FCFS. Системы с постоянными временами обслуживания имеют определяющее свойство: клиенты покидают обслуживающие приборы в том же самом порядке, в котором они приняты для обслуживания.

Исторические замечания по M/D/n

Системы организации очереди с Пуассоновским потоком вызовов и постоянными временами обслуживания были проанализированы первыми. Интуитивно можно было бы думать, что анализ более систем с постоянными временами обслуживания проще, чем с экспоненциально распределенными временами обслуживания, но это явно не так. С экспоненциальным распределением иметь дело проще из-за отсутствия у него памяти: остающееся время "жизни" имеет то же самое распределение, как полное время "жизни" (секция 4.1), и поэтому мы можем забыть о периоде или моменте времени, когда начинается время обслуживания. Постоянные времена занятия требуют, чтобы мы помнили точное время начала. Эрланг первым проанализировал систему с постоянным временем обслуживания M/D/n, FCFS (Brockmeyer, 1948 [11] ):

Эрланг: 1909 n = 1 ошибка для n > 1,

Эрланг: 1917 n = 1, 2, 3 без доказательства,

Эрланг: 1920 n - произвольное, решения для n = 1, 2, 3.

Эрланг получил распределение времени ожидания, но не рассматривал вероятности состояния. Фрей (Fry (1928 [30] )) также исследовал M/D/1 и получил вероятности состояния (уравнения состояния Фрея), используя принцип статистического равновесия Эрланга. Сам Эрланг применял другие теоретические методы.

Кроммелин (1932 [20] , 1934 [21] ), британский телефонный инженер, представил общее решение для M/D/n. Он обобщил уравнения состояний Фрея для произвольного n, и получил распределение времени ожидания, теперь называемое распределением Кроммелина.

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

Хинчин (1932 [63] ) исследовал M/D/n и получил распределение времени ожидания.

Вероятности состояния M/D/1

Вероятности состояний для M/D/1 могут быть получены простым способом, исходя из предположения о статистическом равновесии.

Пусть интенсивность поступления вызовов равна \lambda, а постоянное время занятия - h. Рассматриваем чистую систему с ожиданием и с единственным обслуживающим прибором.

\mbox{Предложенная нагрузка = Обслуженная нагрузка} =\lambda * h < 1,      (13.39) ( 13.39)

то есть

А= Y=\lambda *h = 1 -р(0),

в каждом состоянии, кроме нулевого, обслуженная нагрузка равна 1 Эрл.

Мы рассматриваем два периода (момента времени) t и t+h на расстоянии h. Каждый клиент, обслуживаемый в период t (самое большее один период), может покинуть обслуживающий прибор в период t+h. Клиенты, прибывающие в течение интервала (t, t+ h) , находятся все еще в очереди в период t+h (ожидание или обслуживание).

Процесс поступления вызовов - Пуассоновский процесс. Следовательно, мы имеем Пуассоновское распределение прибытия во временном интервале (t, t+h) :

p(i,h)=p\{j \quad calls \quad in \quad h\}=\frac{(\lambda h)^j}{j!}*e^{- \lambda h}, j=0,1,2, \dots. ( 13.40)

Вероятность того, что данное состояние в период t+h может быть получено из состояния в период t, определяется, принимая во внимание общее количество вызовов и выходы из состояния в течение (t, t+h) . Рассматривая эти периоды, получаем марковскую цепь, внедренную в первоначальный процесс обслуживания нагрузки (рис.13.4).

Мы получаем уравнения состояний Фрея для n = 1 (Fry, 1928 [30] ):

p(j,h)=p\{j \quad calls \quad h\}=\frac{(\lambda h)^j}{j!}*e^{- \lambda h}, j=0,1,2,\dots. ( 13.40)

Выше было установлено:

р(0) = 1-А

и согласно предположению о статистическом равновесии p_t(i) = p_{t +h}(i), мы последовательно находим:

p(1)=(1-A)*\{e^A-1\},\\
p(2)=(1-A)*\{-e^A*(1+A)+e^{2A}\},

и далее

p(i)=(1-A)*\sum_{j=1}^i(-1)^{i-j}*e^{jA}* \left \{\frac{(jA)^{i-j}}{(i-j)!}+\frac{(jA)^{i-j-1}}{(i-j-1)!} \right \}, i=2,3,\dots. ( 13.42)

Последнее слагаемое, соответствующее j = i, всегда равняется e^{iA}, поскольку (-1)! \equiv \infty

В принципе p(0) может быть получено из условия, что все вероятности состояния должны в сумме давать единицу.

 Иллюстрация уравнений состояния Фрея для системы организации очереди M/D/1

Рис. 13.4. Иллюстрация уравнений состояния Фрея для системы организации очереди M/D/1

Средние времена ожидания и период занятости M/D/1

Для Пуассоновского процесса поступления вызовов вероятность задержки D равна вероятности того, что он не находится в нулевом состоянии (свойство PASTA ):

D=A=1-p(0) ( 13.43)

W обозначает среднее время ожидания для всех клиентов, и w обозначает среднее время ожидания для клиентов, которые стоят в очереди (положительное значение времени ожидания). Мы имеем для любой системы организации очереди (3.20):

w=\frac WD ( 13.44)

W и w могут быть легко получены, используя формулу Поллячека-Хинчина (13.2):

W=\frac{A*h}{2(1-A)}, ( 13.45)
w=\frac{h}{1(1-A)}. ( 13.46)

Средняя величина периода занятости была получена для M/G/1 в (13.7) и иллюстрирована для постоянных времен обслуживания на рис.13.1:

m_{T_1}=\frac{h}{1-A}. ( 13.47)

Среднее время ожидания для задержанных клиентов, таким образом, - половина периода занятости. Это показывает, что вызовы от клиентов прибывают случайно в течение периода занятости, но мы знаем, что ни один вызов не поступает в течение последнего времени обслуживания периода занятости.

Можно показать, что распределение числа вызовов, прибывающих в течение периода занятости, можно выразить распределением Бореля (Borel ):

B(i)=\frac{(iA)^{i-1}}{i!} e^{-iA}, \quad i=1,2, \dots ( 13.48)

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

Можно показать, что:

p\{W \le t\}=1-(1- \lambda)* \sum_{j=1}^{\infty} \frac{\{\lambda (j- \tau)\}^{T+j}}{(T+j)!}*e^{-\lambda(j - \tau)}, ( 13.49)

где h=1 выбрано как единица времени, t=T+ \tau , T - целое число, и 0 \le \tau < 1.

Граф распределения времени ожидания имеет нарушения каждый раз, когда время ожидания превышает кратное число постоянного времени занятия. Пример для распределения вероятности дополнения времени ожидания (p(W > t)) показан на рис.13.5.

 Распределение дополнения времени ожидания для всех клиентов в системе организации очереди M/M/1 и M/D/1 для дисциплины очереди (FCFS). Единица времени = среднее время обслуживания. Заметим, что среднее время ожидания для M/D/1 - только половина времени для M/M/1.

Рис. 13.5. Распределение дополнения времени ожидания для всех клиентов в системе организации очереди M/M/1 и M/D/1 для дисциплины очереди (FCFS). Единица времени = среднее время обслуживания. Заметим, что среднее время ожидания для M/D/1 - только половина времени для M/M/1.

Формула (13.49) неудобна для числовой оценки. Можно показать (Iversen, 1982 [39] ), что формула времени ожидания может быть записана в более компактной форме, например, введенной Эрлангом в 1909 г.:

p\{W \le t \}=(1- \lambda )*\sum_{j=0}^T \frac{\{\lambda (j-t)\}^j}{j!}*e^{- \lambda (j-t)} ( 13.50)

и может применяться для числовой оценки при малых временах ожидания.

Для больших времен ожидания мы обычно интересуемся составными значениями t. Можно показать (Iversen, 1982 [39] ), что для составного значения t

p\{W \le t\}=p(0)+p(1)+ \dots p(t). ( 13.51)

Вероятности состояния p(i) вычисляются наиболее точно, с помощью рекурсивной формулы, основанной на уравнениях состояний (13.42) Фрея:

p(i+1)=\frac{1}{p(0,h)} \left \{ p(i)-\{p(0)+p(1)\}*p(i,h)-\sum_{j=2}^ip(j)*p(i-j+1,h) \right \} ( 13.52)

Для несоставных времен ожидания можно выразить распределение времени ожидания с помощью комбинации составных времен ожидания.

Если мы примем h = 1, то (13.50) может быть биноминальным разложением, записанным с помощью степени \tau, где t=T+ \tau; T - целое число, 0 \le \tau < 1.

p\{W \le T+ \tau\}=e^{\lambda \tau} \sum_{j=0}^T \frac{(- \lambda \tau)^j}{j!}*p\{W \le T-j\}, ( 13.53)

где p\{W \le T-j\} определяется в (13.51).

Очень точная числовая оценка может быть получена при использовании (13.51), (13.52) и (13.53).

Вероятности состояния: M/D/n

При выводе уравнений состояния Фрея (13.41) мы получаем больше комбинаций:

p_{t+h}(i)=\left \{ \sum_{j=0}^n p_t(j) \right \} p(i,h)+ \sum_{j=n+1}^{n+j}p_t(j)*p(n+i-j,h) ( 13.54)

При условии статистического равновесия (A < n) можно вычислить абсолютные моменты времени выхода из состояния:

p(i)=\left \{ \sum_{j=0}^n p(j) \right \} p(i,h)+ \sum_{j=n+1}^{n+1}p(j)*p(in+i-j,h), i=0,1, \dots ( 13.55)

Если мы знаем первые n вероятности состояний \{p(0), p(1)\dots, p (n-1)\} , система уравнений (13.55) может быть решена непосредственно подстановкой. Практически мы можем получить числовые значения, подставляя приблизительный набор значений для \{p(0), p(1) \dots , p(n - 1)\} . Затем, заменяя эти значения, согласно рекурсивной формуле (13.55), получаем новые значения. После нескольких приближений мы получим точные результаты.

Явное математическое решение может быть получено с помощью производящих функций (Erlang [11] стр. 75 (83)).

Нияз Сабиров
Нияз Сабиров

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

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

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

алексей оглы
алексей оглы
Россия
рафич Салахиев
рафич Салахиев
Россия