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

Системы с ожиданием

< Лекция 12 || Лекция 13: 123 || Лекция 14 >
Состояния терминалов и характеристики нагрузки

Характеристики качества работы легко могут быть получены на основе аналогии с классической системой с потерями Эрланга (12.37). Заменяя p(0) на E_{1,S} (\varrho) , получаем, что компьютер работает с вероятностью \{1 - E_{1,S} (\varrho)\} . Тогда среднее число терминалов, обслуживаемых компьютером, равно:

n_s=1-E_{1,S}(\varrho). ( 12.38)

Среднее число ожидающих терминалов равно:

n_w=S-n_s-n_t=S-\{1-E_{1,S}(\varrho)- \varrho *\{1-E_{1,S}(\varrho)\}\\

\quad = S-\{1-E_{1,S}(\varrho)\}\{1+ \varrho\}. ( 12.40)

Если мы рассматриваем случайный терминал в случайный момент времени, то получаем:

p(\mbox{обслуживающий терминал})p_s=\frac{n_s}{S}=\frac{1-E_{1,S}(\varrho}{S}, ( 12.41)
p(\mbox{терминал в паузе})p_t=\frac{n_t}{S}=\frac{\varrho(1-E_{1,S}(\varrho)}{S}, ( 12.42)
p(\mbox{ожидающий терминал})=p_w=\frac{n_w}{S}=1-\frac{\{1-E_{1,S}(\varrho)\}\{1+ \varrho\}}{S} ( 12.43)

Мы также интересуемся временем реакции R, которая имеет среднюю величину m_r = m_w + m_s. Применяя формулу Литла L = \lambda W к терминалам, установленным на ожидание, и к компьютеру, мы, соответственно, получаем (обозначая скорость обращения заявок \lambda ):

\frac{1}{\lambda}=\frac{m_t}{n_t}=\frac{m_w}{n_w}=\frac{m_s}{n_s}=\frac{m_r}{n_w+n_s}, ( 12.44)

или

m_r=\frac{n_w+n_s}{n_s}*m_s=\frac{S_n_t}{n_s}*m_s.

Используя (12.38) и (12.44) \frac{n_t}{n_s}=\frac{m_t}{m_s}, мы получим:

m_r=\frac{S}{n_s}*m_s-m_t\\
m_r=\frac{S}{1-E_{1,S}(\varrho)}*m_s-m_t. ( 12.45)

Таким образом, среднее время ответа не зависит от типа распределения времени, поскольку оно базируется на формуле Литла (12.38) и (12.44). Однако, E_{1, S} (varrho) будет зависеть от типов распределений, так же, как это было в B-формуле Эрланга. Если время обслуживания компьютера является экспоненциально распределенным (средняя величина ms = 1/ \mu ), то E_{1,S}(\varrho) будет определяться формулой (12.37). Pис.12.8 в этом случае показывает время реакции как функцию числа терминалов.

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

K=\frac{m_t+m_s}{m_s}\\

= \varrho+1 ( 12.46)

K - подходящий параметр для описания точки насыщения системы. Среднее время ожидания для произвольного терминала может быть получено из (12.45):

m_w = m_r- m_s
 Фактическое среднее время реакции, определяется как функция числа терминалов

Рис. 12.8. Фактическое среднее время реакции, определяется как функция числа терминалов

Сервисный коэффициент - \varrho = 30. Среднее время реакции переходит в прямую линию, пересекающую ось x при S = 30 терминалов. Среднее виртуальное время реакции для системы с S терминалами равно фактическому среднему времени реакции для системы с S +1 терминалами (теорема моментов поступления, теорема 8.1).

Пример 12.5.2: Компьютер в режиме разделения времени

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

Pис.12.9 показывает нагрузку времени ожидания в Эрл для компьютера и для одного терминала. Соответствующая надбавка затрат и сумма времен ожидания для компьютера и для всех терминалов дает стоимость ожидания.


Рис. 12.9.

Нагрузка времени ожидания (соотношение времени, потраченного на ожидание), измеренная в Эрл для компьютера и соответственно для терминалов при интерактивной системе организации очереди (коэффициент обслуживания - \varrho = 30).

Например, на рис.12.9 мы получаем минимальные полные затраты ожидания приблизительно для 45 терминалов. Стоимость ожидания компьютера - в сотню раз больше стоимости времени одного терминала. При 31 терминале и одном компьютере каждый терминал тратит 11,4 % времени для ожидания. Если отношение стоимости - 31, то 31 - оптимальное число терминалов. Однако должны быть учтены несколько других факторов.

Пример 12.5.3: Потери по нагрузке

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

Удельная нагрузка будет (8.8):

a=\frac{\beta}{1+\beta}=\frac{m_s}{m_t+m_s}

На один источник обслуженная нагрузка:

y=\frac{m_s}{m_t+m_w+m_s}

Потери по нагрузке равны:

C=\frac{a-y}{a}\\
=1-\frac{m_t-m_s}{m_t+m_w+m_s}=\frac{m_w}{m_t+m_w+m_s},\\
C=p_w

Потери по нагрузке становятся равными соотношению времени, потраченного на ожидание. Для системы Эрланга с ожиданием потери по нагрузке являются нулевыми, потому что вся предложенная нагрузка обслуживается.

Модель восстановления машин с п серверами

Вышеупомянутая модель легко может быть обобщена на п компьютеров.

Диаграмма переходов показана на рис.12.10 Вероятности устойчивых состояний равны:

p(i)={S\choose i}(\frac{\gamma}{\mu})^ip(0), 0 \le i \le n,\\
p(i)=\frac{(S-n)!}{(S-i)!}(\frac{\gamma}{n \mu})^{i-n}*p(n), n \le i \le S. ( 12.47)

где мы имеем нормировочное ограничение:

\sum_{i=0}^Sp(i)=1 ( 12.48)
 Диаграмма переходов состояний для модели восстановления машин с S терминалами и n компьютерами

Рис. 12.10. Диаграмма переходов состояний для модели восстановления машин с S терминалами и n компьютерами

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

Произвольный терминал - в случайный момент времени может находиться в одном из трех возможных состояний:

p_s = p\{\mbox{ терминал, обслуживаемый компьютером}\};p_s = p\{ \mbox{терминал, ожидающий обслуживания}\}; p_w = p\{ \mbox{терминал в паузе}\}.

Мы имеем:

p_s=\frac 1S\left \{\sum_{i=0}^ni*p(i)+\sum_{i=n+1}^Sn*p(i) \right \}, ( 12.49)
p_t=p_s*\frac{\mu}{\gamma}, ( 12.50)
p_w=1-p_s-p_t. ( 12.51)

Среднее использование компьютеров равно:

\alpha=\frac{p_s}{n}*S=\frac{n_s}{n}. ( 12.52)

Среднее время ожидания для терминала равно:

W=\frac{p_w}{p_s}*\frac{1}{\mu}. ( 12.53)

Иногда p_w называют коэффициентом потерь терминалов, и аналогично (1 - \alpha) называют коэффициентом потерь компьютеров (рис.12.9).

Пример 12.5.4: Числовой пример

Следующие числовые примеры иллюстрируют то, что мы получаем самое высокое использование для больших значений nS ). Рассмотрим систему с S/n = 30 и \mu / \gamma = 30 для большого числа компьютеров (в этом случае p_t = \alpha ).

n 1 2 4 8 16
p_s 0.0289 0.0300 0.0307 0.0313 0.0316
p_w 0.1036 0.0712 0.0477 0.0311 0.0195
p_t 0.8675 0.8989 0.9215 0.9377 0.9489
a 0.8675 0.8989 0.9215 0.9377 0.9489
W[\mu^1] 3.5805 2.3754 1.5542 0.9945 0.6155

Оптимизация модели восстановления машин

В этой секции мы оптимизируем модель восстановления машин тем же способом, как сделал это Пальм в 1947 г. Заметим, что модель для единственного ремонтника (одного обслуживающего прибора) похожа на систему с потерями Эрланга, которую мы оптимизировали в Лекции 7. Мы можем видеть, что одна и та же модель может быть оптимизирована несколькими способами.

Рассмотрим терминальную систему с одним компьютером и S терминалами и найдем оптимальное значение S. Примем следующую структуру затрат:

c_t - стоимость на терминал в единицу времени при паузе (размышление или работа только с терминалом);

c_w - стоимость на один терминал в единицу времени при ожидании;

c_s - стоимость на один терминал в единицу времени при обслуживании;

c_a - стоимость на компьютер в единицу времени.

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

Общая стоимость c_0 в единицу времени, когда терминал в паузе:

p_t*c_0=p_t*c_t+p_s*c_s+p_w*c_w+\frac 1S*c_a ( 12.54)

Рис. 12.11.

Модель восстановления машины. Общие стоимости, данные в (12.57), показывают общую стоимость как функцию числа терминалов для сервисного отношения \varrho = 25 и отношения стоимости r = 1/25 (см. рис. 7.6)

Мы хотим минимизировать с_0. Сервисное отношение \varrho = m_t /m_s равно p_t /p_s. Вводя отношение стоимости r = с_w /с_a, получаем:

c_0=c_t+\frac{p_s}{p_t}*c_s+\frac{p_w*c_w+\frac 1S*c_a}{p_t}\\
=c_t+\frac{1}{\varrho}*c_s+c_a*\frac{r*p_w+(1/S)}{p_t}, ( 12.55)

Это выражение должно быть свернуто как функция S. Учитывая, что только последний член зависит от числа терминалов, получаем:

min_S\{c_0\}=min_S \left \{\frac{r*p_w+(1/S)}{p_t} \right \}\\
=min_S \left \{\frac{r*(n_w/S)+(1/S)}{n_t/S} \right \}\\
=min_S \left \{ \frac{r*n_w+1}{n_t} \right \} ( 12.56)
=min_S \left \{ \frac{r[S-\{1-E_{1,S}(\varrho)\}\{1+ \varrho\}]+1}{\{1-E_{1,S}(\varrho)\}*\varrho}\right \}\\
=min_s \left \{ \frac{r*S+1}{\{1-E_{1,S}(\varrho)\}* \varrho}+1+\frac{1}{\varrho} \right \}, ( 12.57)

Мы замечаем, что минимум не зависит от c_t и c_s и что только отношение r = c_w /c_a содержит характеристики стоимости. Числитель соответствует (7.31), тогда как знаменатель соответствует обслуженной нагрузке в соответствующей системе с потерями. Таким образом мы минимизируем стоимость обслуживания в Эрл. для соответствующей системе с потерями. На рис.12.10 показан пример. Мы замечаем, что результат отличается от результата, полученного с использованием принципа Мо для системы с потерями Эрланга (рис.7.6), где мы оптимизируем прибыль.

Краткие итоги

  • В этой лекции рассматривается система с ожиданием, имеющая бесконечное число мест ожидания. Когда все n обслуживающих приборов заняты, поступивший вызовов ставится в очередь и ждет, пока не освободится хотя бы один обслуживающий прибор.
  • Когда хотя бы один обслуживающий прибор свободен, клиенты не могут оставаться в очереди (полная доступность). o Рассмотрены два случая потоков нагрузки, Пуассоновский поток вызовов ( бесконечное число источников ) и экспоненциально распределенное время обслуживания ( PCT-I ) - Эрланговская система с ожиданием ( M/M/n ).
  • Также рассмотрены ограниченное число источников и экспоненциально распределенное время обслуживания ( PCT-II ). Это - модель Пальма, называемая моделью восстановления машин.
  • Система с ожиданием типа M/M/n обслуживает Пуассоновский поток вызовов ( M ) и имеет экспоненциальные время обслуживания ( M ) и n обслуживающих приборов при бесконечном числе мест ожидания.
  • Состояние системы определяется как общее количество пользователей в системе (или в обслуживании, или ожидающих в очереди).
  • Вероятность, что система c ожиданием находится в состоянии i, равна

    p(i)=\begin{cases}
p(o)* \frac{A^i}{i!}, & 0 \le i \le n\\
p(n)*(\frac An)^{i-n}=p(0)*\frac{A^i}{n!*n^{i-n}} & i \ge n
\end{cases}.

    где

    p(0)=\frac{1}{\sum_{i=0}^{n-1} \frac{A^i}{i!}+\frac{A^n}{n!} \frac{n}{n-A}}, A < n.
  • Вероятность ожидания p\{W > 0\} = E_{2,n} ( A ) в системе c ожиданием равна:

    E_{2,n}(A)=\frac{\frac{A^n}{n!}\frac{n}{n-A}}{1+\frac A1+\frac{A^2}{2!}+ \dots +\frac{A^{n-1}}{(n-1)!}+\frac{A^n}{n!}\frac{n}{n-A}}, A < n.

    ( формула Эрланга для систем с ожиданием ).

  • Поскольку мы имеем очень точную рекурсивную формулу для числовой оценки B-формулы (7.29) Эрланга, можно использовать следующие отношения для того, чтобы получить числовые значения для C-формулы:

    E_{2,n}(A)=\frac{E_{1,n}(A)}{1-A\{1-E_{1,n}(A)\}/n}
  • Мы должны отличать длину очереди в произвольный момент времени и длину очереди, когда есть клиенты, стоящие в очереди.
  • Длина очереди L в произвольной момент времени называется виртуальной длиной очереди. Эта длина очереди определяется для произвольного клиента:

    L_n=E_{2,n}(A)*\frac{A}{n-A} \qquad A < n
  • Средняя длина очереди со временем ожидания больше нуля (длина очереди, когда есть клиенты, стоящие в очереди):

    L_{nq}=\frac{L_n}{p\{L > 0\}

    где L - случайная переменная, обозначающая длину очереди.

  • Теорема Литла говорит, что средняя длина очереди в произвольный момент времени равна интенсивности прибытия, умноженной на среднее время ожидания:

    L_n=\lambda W_n
  • Представляют интерес две характеристики: среднее время ожидания W для всех клиентов и среднее времени ожидания w для клиентов, для которых время ожидания имеет положительное значение.
  • Среднее время обслуживания для системы равно:

    W_n=E_{2,n}(A)*\frac{s}{n-A}.
  • Среднее время ожидания w для задержанных вызовов:

    W_n=\frac{s}{n-A}.
  • Mo предложил свой принцип для систем организации очереди. Согласно этому принципу оптимальное решение равно:

    F_{W,n_i-1}(A) > \frac{c_i}{\vartheta} \ge F_{W, n_i}(A), \quad i=1,2, \dots, k

    где c_i - переменная стоимость канала;

    F_{W,n}(A) - функция увеличения времени ожидания;

    \vartheta (тета) - множитель Лагранжа.

  • Системы организации очереди, где дисциплина обслуживания зависит от времени поступления вызовов, все имеют одни и те же средние времена ожидания. В этом случае стратегия оказывает влияние только на распределение времен ожидания для каждого отдельного клиента.
  • Вероятность p (W \le t) = F (t) , что положительное время ожидания t превосходит заданную величину, равна вероятности, что в Пуассоновском потоке вызовов с интенсивностью n \mu, по крайней мере, (k+1) вызовов поступят в течение интервала t.
  • Исследование распределения времени ожидания упрощается в случае дисциплины (First Come First Served - "Первый прибыл - Первый обслужен"). Та же дисциплина обозначается FIFO (First In First Out). Эта дисциплина также называется дисциплина обслуживания в порядке поступления.
  • При поступлении вызова в систему с ожиданием и дисциплиной обслуживания FCFS со всеми ранее занятыми обслуживающими приборами можно:

    1. рассчитать число k ждущих клиентов. Полное время ожидания тогда будет (k+1) - Эрланговское распределение;
    2. если не учитывать этих ждущих клиентов, то время ожидания становится экспоненциально распределенным.
  • Среднее время пребывания заявки в системе в случае одного обслуживающего прибора:

    m=\frac{1}{\mu - \lambda}
  • Режим разделения времени - это лучшее решение для оптимального обслуживания большой группы источников нагрузки, использующих, например, терминалы, подключенные к универсальному компьютеру.
  • Вероятности состояния модели ремонта машин с одним компьютером и S терминалами справедливы в течение произвольных времен пауз (размышления и работа с компьютером), когда времена обслуживания компьютером являются экспоненциально распределенными.
  • Модель ремонта машин легко может быть обобщена на n компьютеров и может быть оптимизирована для терминальной системы с одним компьютером и S терминалами.
< Лекция 12 || Лекция 13: 123 || Лекция 14 >
Нияз Сабиров
Нияз Сабиров
Стоимость "обучения"
Елена Сапегова
Елена Сапегова
диплом
Александр Орешин
Александр Орешин
Россия, Ростов-на-Дону
Андрей Лучицкий
Андрей Лучицкий
Россия