4. ТЕОРИЯ МАССОВОГО ОБСЛУЖИВАНИЯ

4.1. Классификация систем массового обслуживания и их показатели эффективности

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

СМО могут быть классифицированы по признаку организации обслуживания следующим образом:

Системы с отказами не имеют очередей.

Системы с ожиданием имеют очереди.

Заявка, поступившая в момент, когда все каналы обслуживания заняты:

Покидает систему с отказами;

Становится в очередь на обслуживание в системах с ожиданием при неограниченной очереди или на свободное место при ограниченной очереди;

Покидает систему с ожиданием при ограниченной очереди, если в этой очереди нет свободного места.

В качестве меры эффективности экономической СМО рассматривают сумму потерь времени:

На ожидание в очереди;

На простои каналов обслуживания.

Для всех видов СМО используются следующие показатели эффективности :

- относительная пропускная способность - это средняя доля поступающих заявок, обслуживаемых системой;

- абсолютная пропускная способность - это среднее число заявок, обслуживаемых системой в единицу времени;

- вероятность отказа - это вероятность того, что заявка покинет систему без обслуживания;

- среднее число занятых каналов - для многоканальных СМО.

Показатели эффективности СМО рассчитываются по формулам из специальных справочников (таблиц). Исходными данными для таких расчетов являются результаты моделирования СМО.

4.2. Моделирование системы массового обслуживания:

основ­ные параметры, граф состояний

При всем многообразии СМО они имеют общие черты , которые позволяют унифицировать их моделирование для нахождения наиболее эффективных вариантов организации таких систем .

Для моделирования СМО необходимо иметь следующие исходные данные:

Основные параметры;

Граф состояний.

Результатами моделирования СМО являются вероятности ее состояний, через которые выражаются все показатели ее эффективности.

Основные параметры для моделирования СМО включают:

Характеристики входящего потока заявок на обслуживание;

Характеристики механизма обслуживания.

Рассмотрим характеристики потока заявок .

Поток заявок - последовательность заявок, поступающих на обслуживание.

Интенсивность потока заявок - среднее число заявок, поступающих в СМО в единицу времени.

Потоки заявок бывают простейшими и отличными от простейших.

Для простейших потоков заявок используются модели СМО.

Простейшим , или пуассоновским называется поток, являющийся стационарным , одинарным и в нем отсутствуют последействия .

Стационарность означает неизменность интенсивности поступления заявок с течением времени.

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

Отсутствие последействия заключается в том, что число заявок, поступивших в СМО за один интервал времени, не влияет на количество заявок, полученных за другой интервал времени.

Для отличных от простейших потоков заявок используются имитационные модели.

Рассмотрим характеристики механизма обслуживания .

Механизм обслуживания характеризуется:

- числом каналов обслуживания ;

Производительностью канала, или интенсивностью обслуживания - средним числом заявок, обслуживаемых одним каналом в единицу времени;

Дисциплиной очереди (например, объемом очереди , порядком отбора из очереди в механизм обслуживания и т. п.).

Граф состояний описывает функционирование системы обслуживания как переходы из одного состояния в другое под действием потока заявок и их обслуживания.

Для построения графа состояний СМО необходимо:

Составить перечень всех возможных состояний СМО;

Представить перечисленные состояния графически и отобразить возможные переходы между ними стрелками;

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

4.3. Вычисление вероятностей состояний

системы массового обслуживания

Граф состояний СМО со схемой "гибели и рождения" представляет собой линейную цепочку, где каждое из средних состояний имеет прямую и обратную связь с каждым из соседних состояний, а крайние состояния только с одним соседним:

Число состояний в графе на единицу больше, чем суммарное число каналов обслуживания и мест в очереди.

СМО может быть в любом из своих возможных состояний, поэтому ожидаемая интенсивность выхода из какого-либо состояния равна ожидаемой интенсивности входа системы в это состояние. Отсюда система уравнений для определения вероятностей состояний при простейших потоках будет иметь вид:

где - вероятность того, что система находится в состоянии

- интенсивность перехода, или среднее число переходов системы в единицу времени из состояния в состояние .

Используя эту систему уравнений, а также уравнение

вероятность любого -ого состояния можно вычислить по следующему общему правилу :

вероятность нулевого состояния рассчитывается как

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

Выводы по четвертому разделу

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

Моделирование систем массового обслуживания осуществляется для нахождения наиболее эффективных вариантов их организации и предполагает следующие исходные данные для этого: основные параметры, граф состояний. К таким данным относятся следующие: интенсивность потока заявок, количество каналов обслуживания, интенсивность обслуживания и объем очереди. Число состояний в графе на единицу больше, чем сумма числа каналов обслуживания и мест в очереди.

Вычисление вероятностей состояний системы массового обслуживания со схемой «гибели и рождения» осуществляется по общему правилу.

Вопросы для самопроверки

Какие системы называются системами массового обслуживания?

Как классифицируются системы массового обслуживания по признаку их организации?

Какие системы массового обслуживания называются системами с отказами, а какие – с ожиданием?

Что происходит с заявкой, поступившей в момент времени, когда все каналы обслуживания заняты?

Что рассматривают в качестве меры эффективности экономической системы массового обслуживания?

Какие используются показатели эффективности системы массового обслуживания?

Что служит исходными данными для расчетов показателей эффективности систем массового обслуживания?

Какие исходные данные необходимы для моделирования систем массового обслуживания?

Через какие результаты моделирования системы массового обслуживания выражают все показатели ее эффективности?

Что включают основные параметры для моделирования систем массового обслуживания?

Чем характеризуются потоки заявок на обслуживание?

Чем характеризуются механизмы обслуживания?

Что описывает граф состояний системы массового обслуживания

Что необходимо для построения графа состояний системы массового обслуживания?

Что представляет собой граф состояний системы массового обслуживания со схемой «гибели и рождения»?

Чему равно число состояний в графе состояний системы массового обслуживания?

Какой вид имеет система уравнений для определения вероятностей состояний системы массового обслуживания?

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

Примеры решения задач

1. Построить граф состояний системы массового обслуживания и привести основные зависимости ее показателей эффективности.

а) n-канальная СМО с отказами (задача Эрланга)

Основные параметры:

Каналов ,

Интенсивность потока ,

Интенсивность обслуживания .

Возможные состояния системы:

Все каналов заняты ( заявок в системе).

Граф состояний:

Относительная пропускная способность ,

Вероятность отказа ,

Среднее число занятых каналов .

б) n-канальная СМО с m-ограниченной очередью

Возможные состояния системы:

Все каналы свободны (ноль заявок в системе);

Один канал занят, остальные свободны (одна заявка в системе);

Два канала заняты, остальные свободны (две заявки в системе);

...................................................................................

Все каналы заняты, две заявки в очереди;

Все каналы заняты, заявок в очереди.

Граф состояний:

в) Одноканальная СМО с неограниченной очередью

Возможные состояния системы:

Все каналы свободны (ноль заявок в системе);

Канал занят, ноль заявок в очереди;

Канал занят, одна заявка в очереди;

...................................................................................

Канал занят, заявка в очереди;

....................................................................................

Граф состояний:

Показатели эффективности системы:

,

Среднее время пребывания заявки в системе ,

,

,

Абсолютная пропускная способность ,

Относительная пропускная способность .

г) n-канальная СМО с неограниченной очередью

Возможные состояния системы:

Все каналы свободны (ноль заявок в системе);

Один канал занят, остальные свободны (одна заявка в системе);

Два канала заняты, остальные свободны (две заявки в системе);

...................................................................................

Все каналов заняты ( заявок в системе), ноль заявок в очереди;

Все каналы заняты, одна заявка в очереди;

....................................................................................

Все каналы заняты, заявок в очереди;

....................................................................................

Граф состояний:

Показатели эффективности системы:

Среднее число занятых каналов ,

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

Среднее число заявок в очереди ,

Среднее время пребывания заявки в очереди .

2. Вычислительный центр имеет три ЭВМ. В центр поступает на решение в среднем четыре задачи в час. Среднее время решения одной задачи - полчаса. Вычислительный центр принимает и ставит в очередь на решение не более трех задач. Необходимо оценить эффективность центра.

РЕШЕНИЕ. Из условия ясно, что имеем многоканальную СМО с ограниченной очередью:

Число каналов ;

Интенсивность потока заявок (задача / час);

Время обслуживания одной заявки (час / задача), интенсивность обслуживания (задача / час);

Длина очереди .

Перечень возможных состояний:

Заявок нет, все каналы свободны;

Один канал занят, два свободны;

Два канала заняты, один свободен;

Три канала заняты;

Три канала заняты, одна заявка в очереди;

Три канала заняты, две заявки в очереди;

Три канала заняты, три заявки в очереди.

Граф состояний:

Рассчитаем вероятность состояния :

Показатели эффективности:

Вероятность отказа (все три ЭВМ заняты и три заявки стоят в очереди)

Относительная пропускная способность

Абсолютная пропускная способность

Среднее число занятых ЭВМ

3. (Задача с использованием СМО с отказами.) В ОТК цеха работают три контролера. Если деталь поступает в ОТК, когда все контролеры заняты обслуживанием ранее поступивших деталей, то она проходит непроверенной. Среднее число деталей, поступающих в ОТК в течение часа, равно 24, среднее время, которое затрачивает один контролер на обслуживание одной детали, равно 5 мин. Определить вероятность того, деталь пройдет ОТК необслуженной, насколько загружены контролеры и сколько их необходимо поставить, чтобы (* - заданное значение ).

РЕШЕНИЕ. По условию задачи , тогда .

1) Вероятность простоя каналов обслуживания:

,

3) Вероятность обслуживания:

4) Среднее число занятых обслуживанием каналов:

.

5) Доля каналов, занятых обслуживанием:

6) Абсолютная пропускная способность:

При . Произведя аналогичные расчеты для , получим

Так как , то произведя расчеты для , получим

ОТВЕТ. Вероятность того, что при деталь пройдет ОТК необслуженной, составляет 21%, и контролеры будут заняты обслуживанием на 53%.

Чтобы обеспечить вероятность обслуживания более 95%, необходимо не менее пяти контролеров.

4. (Задача с использованием СМО с неограниченным ожиданием.) Сберкасса имеет трех контролеров-кассиров () для обслуживания вкладчиков . Поток вкладчиков поступает в сберкассу с интенсивностью чел./ч. Средняя продолжительность обслуживания контролером-кассиром одного вкладчика мин.

Определить характеристики сберкассы как объекта СМО.

РЕШЕНИЕ. Интенсивность потока обслуживания , интенсивность нагрузки .

1) Вероятность простоя контролеров-кассиров в течение рабочего дня (см. предыдущую задачу №3):

.

2) Вероятность застать всех контролеров-кассиров занятыми:

.

3) Вероятность очереди:

.

4) Среднее число заявок в очереди:

.

5) Среднее время ожидания заявки в очереди:

мин.

6) Среднее время пребывания заявки в СМО:

7) Среднее число свободных каналов:

.

8) Коэффициент занятости каналов обслуживания:

.

9) Среднее число посетителей в сберкассе:

ОТВЕТ. Вероятность простоя контролеров-кассиров равна 21% рабочего времени , вероятность посетителю оказаться в очереди составляет 11,8%, среднее число посетителей в очереди 0,236 чел., среднее время ожидания посетителями обслуживания 0,472 мин.

5. (Задача с применением СМО с ожиданием и с ограниченной длиной очереди.) Магазин получает ранние овощи из пригородных теплиц. Автомобили с грузом прибывают в разное время с интенсивностью машин в день. Подсобные помещения и оборудование для подготовки овощей к продаже позволяют обрабатывать и хранить товар, привезенный двумя автомашинами (). В магазине работают три фасовщика (), каждый из которых в среднем может обрабатывать товар с одной машины в течение ч. Продолжительность рабочего дня при сменной работе составляет 12 ч.

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

РЕШЕНИЕ. Определим интенсивность загрузки фасовщиков:

Авт./дн.

1) Найдем вероятность простоя фасовщиков при отсутствии машин (заявок):

причем 0!=1,0.

2) Вероятность отказа в обслуживании:

.

3) Вероятность обслуживания:

Так как , произведем аналогичные вычисления для , получим), при этом вероятность полной обработки товара будет .

Задания для самостоятельной работы

Для каждой из следующих ситуаций определить:

a) к какому классу относится объект СМО;

b) число каналов ;

c) длину очереди ;

d)интенсивность потока заявок ;

e) интенсивность обслуживания одним каналом;

f) количество всех состояний объекта СМО.

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

a) ОО – одноканальная с отказами; МО – многоканальная с отказами; ОЖО – одноканальная с ожиданием с ограниченной очередью; ОЖН - одноканальная с ожиданием с неограниченной очередью; МЖО – многоканальная с ожиданием с ограниченной очередью; МЖН - многоканальная с ожиданием с неограниченной очередью;

b) =… (единиц);

c) =… (единиц);

d) =ххх/ххх (единиц /мин);

e) =ххх/ххх (единиц /мин);

f) (единиц).

1. Дежурный по администрации города имеет пять телефонов. Телефонные звонки поступают с интенсивностью 90 заявок в час, средняя продолжительность разговора составляет 2 мин.

2. На стоянке автомобилей возле магазина имеются 3 места, каждое из которых отводится под один автомобиль. Автомобили прибывают на стоянку с интенсивностью 20 автомобилей в час. Продолжительность пребывания автомобилей на стоянке составляет в среднем 15 мин. Стоянка на проезжей части не разрешается.

3. АТС предприятия обеспечивает не более 5 переговоров одновременно. Средняя продолжительность разговоров составляет 1 мин. На станцию поступает в среднем 10 вызовов в сек.

4. В грузовой речной порт поступает в среднем 6 сухогрузов в сутки. В порту имеются 3 крана, каждый из которых обслуживает 1 сухогруз в среднем за 8 ч. Краны работают круглосуточно. Ожидающие обслуживания сухогрузы стоят на рейде.

5. В службе «Скорой помощи» поселка круглосуточно дежурят 3 диспетчера, обслуживающие 3 телефонных аппарата. Если заявка на вызов врача к больному поступает, когда диспетчеры заняты, то абонент получает отказ. Поток заявок составляет 4 вызова в минуту. Оформление заявки длится в среднем 1,5 мин.

6. Салон-парикмахерская имеет 4 мастера. Входящий поток посетителей имеет интенсивность 5 человек в час. Среднее время обслуживания одного клиента составляет 40 мин. Длина очереди на обслуживание считается неограниченной.

7. На автозаправочной станции установлены 2 колонки для выдачи бензина. Около станции находится площадка на 2 автомашины для ожидания заправки. На станцию прибывает в среднем одна машина в 3 мин. Среднее время обслуживания одной машины составляет 2 мин.

8. На вокзале в мастерской бытового обслуживания работают три мастера. Если клиент заходит в мастерскую, когда все мастера заняты, то он уходит из мастерской, не ожидая обслуживания. Среднее число клиентов, обращающихся в мастерскую за 1 ч, равно 20. Среднее время, которое затрачивает мастер на обслуживание одного клиента, равно 6 мин.

9. АТС поселка обеспечивает не более 5 переговоров одновременно. Время переговоров в среднем составляет около 3 мин. Вызовы на станцию поступают в среднем через 2 мин.

10. На автозаправочной станции (АЗС) имеются 3 колонки. Площадка при станции, на которой машины ожидают заправку, может вместить не более одной машины, и если она занята, то очередная машина, прибывшая к станции, в очередь не становится, а проезжает на соседнюю станцию. В среднем машины прибывают на станцию каждые 2 мин. Процесс заправки одной машины продолжается в среднем 2,5 мин.

11. В небольшом магазине покупателей обслуживают два продавца. Среднее время обслуживания одного покупателя – 4 мин. Интенсивность потока покупателей – 3 человека в минуту. Вместимость магазина такова, что одновременно в нем в очереди могут находиться не более 5 человек. Покупатель, пришедший в переполненный магазин, когда в очереди уже стоят 5 человек, не ждет снаружи и уходит.

12. Железнодорожную станцию дачного поселка обслуживает касса с двумя окнами. В выходные дни, когда население активно пользуется железной дорогой, интенсивность потока пассажиров составляет 0,9 чел./мин. Кассир затрачивает на обслуживание пассажира в среднем 2 мин.

Для каждой из указанных в вариантах СМО интенсивность потока заявок равна и интенсивность обслуживания одним каналом . Требуется:

Составить перечень возможных состояний;

Построить граф состояний по схеме "гибели и размножения".

В ответе указать для каждой задачи:

Количество состояний системы;

Интенсивность перехода из последнего состояния в предпоследнее.

Вариант № 1

1. одноканальная СМО с очередью длиной в 1 заявку

2. 2-канальная СМО с отказами (задача Эрланга)

3. 31-канальная СМО с 1-ограниченной очередью

5. 31-канальная СМО с неограниченной очередью

Вариант № 2

1. одноканальная СМО с очередью длиной в 2 заявки

2. 3-канальная СМО с отказами (задача Эрланга)

3. 30-канальная СМО с 2-ограниченной очередью

4. Одноканальная СМО с неограниченной очередью

5. 30-канальная СМО с неограниченной очередью

Вариант № 3

1. одноканальная СМО с очередью длиной в 3 заявки

2. 4-канальная СМО с отказами (задача Эрланга)

3. 29-канальная СМО с 3-ограниченной очередью

4. Одноканальная СМО с неограниченной очередью

5. 29-канальная СМО с неограниченной очередью

Вариант № 4

1. одноканальная СМО с очередью длиной в 4 заявки

2. 5-канальная СМО с отказами (задача Эрланга)

3. 28-канальная СМО с 4-ограниченной очередью

4. Одноканальная СМО с неограниченной очередью

5. 28-канальная СМО с неограниченной очередью

Вариант № 5

1. одноканальная СМО с очередью длиной в 5 заявок

2. 6-канальная СМО с отказами (задача Эрланга)

3. 27-канальная СМО с 5-ограниченной очередью

4. Одноканальная СМО с неограниченной очередью

5. 27-канальная СМО с неограниченной очередью

Вариант № 6

1. одноканальная СМО с очередью длиной в 6 заявок

2. 7-канальная СМО с отказами (задача Эрланга)

3. 26-канальная СМО с 6-ограниченной очередью

4. Одноканальная СМО с неограниченной очередью

5. 26-канальная СМО с неограниченной очередью

Вариант № 7

1. одноканальная СМО с очередью длиной в 7 заявок

2. 8-канальная СМО с отказами (задача Эрланга)

3. 25-канальная СМО с 7-ограниченной очередью

4. Одноканальная СМО с неограниченной очередью

5. 25-канальная СМО с неограниченной очередью

Вариант № 8

1. одноканальная СМО с очередью длиной в 8 заявок

2. 9-канальная СМО с отказами (задача Эрланга)

3. 24-канальная СМО с 8-ограниченной очередью

4. Одноканальная СМО с неограниченной очередью

5. 24-канальная СМО с неограниченной очередью

Вариант № 9

1. одноканальная СМО с очередью длиной в 9 заявок

2. 10-канальная СМО с отказами (задача Эрланга)

3. 23-канальная СМО с 9-ограниченной очередью

4. Одноканальная СМО с неограниченной очередью

5. 23-канальная СМО с неограниченной очередью

Вариант № 10

1. одноканальная СМО с очередью длиной в 10 заявок

2. 11-канальная СМО с отказами (задача Эрланга)

3. 22-канальная СМО с 10-ограниченной очередью

4. Одноканальная СМО с неограниченной очередью

5. 22-канальная СМО с неограниченной очередью

Теория СМО посвящена разработке методов анализа, проектирования и рациональной организации систем, относящихся к различным областям деятельности, таким как связь, вычислительная техника, торговля, транспорт, военное дело. Несмотря на все свое разнообразие, приведенные системы обладают рядом типичных свойств, а именно.

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

Задача анализа СМО заключается в определении ряда показателей ее эффективности, которые можно разделить на следующие группы:

• показатели, характеризующие систему в целом: число n занятых каналов обслуживания, число обслуженных (λ b ), ожидающих обслуживание или получивших отказ заявок (λ c ) в единицу времени и т.д.;
• вероятностные характеристики : вероятность того, что заявка будет обслужена (P обс) или получит отказ в обслуживании (P отк), что все приборы свободны (p 0) или определенное число их занято (p k ), вероятность наличия очереди и т.д.;
• экономические показатели : стоимость потерь, связанных с уходом не обслуженной по тем или иным причинам заявки из системы, экономический эффект, полученный в результате обслуживания заявки, и т.д.

Часть технических показателей (первые две группы) характеризуют систему с точки зрения потребителей , другая часть – характеризует систему с точки зрения её эксплуатационных свойств . Часто выбор перечисленных показателей, может улучшать эксплуатационные свойства системы, но ухудшать систему с точки зрения потребителей и наоборот. Использование экономических показателей позволяет разрешить указанное противоречие и оптимизировать систему с учетом обеих точек зрения.
В ходе выполнения домашней контрольной работы изучаются простейшие СМО. Это системы разомкнутого типа, бесконечный источник заявок в систему не входит. Входной поток заявок, потоки обслуживания и ожидания этих систем являются простейшими. Приоритеты отсутствуют. Системы однофазные.

Многоканальная система с отказами

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

Смешанные системы

1. Система с ограничением на длину очереди .
   Состоит из накопителя (очереди) и узла обслуживания. Заявка покидает очередь и уходит из системы, если в накопителе к моменту ее появления уже находятся m заявок (m – максимально возможноечисло мест в очереди). Если заявка поступила в систему и застала, хотя бы один канал свободным, она мгновенно начинает обслуживаться. Если в момент поступления заявки в систему все каналы заняты, то заявка не покидает систему, а занимает место в очереди. Заявка покидает систему не обслуженной, если к моменту её поступления в систему заняты все каналы обслуживания и все места в очереди.
   Для каждой системы определяется дисциплина очереди. Это система правил, определяющих порядок поступления заявок из очереди в узел обслуживания. Если все заявки и каналы обслуживания равнозначны, то чаще всего действует правило «кто раньше пришел, тот раньше обслуживается».
2. Система с ограничением на длительность пребывания заявки в очереди .
   Состоит из накопителя (очереди) и узла обслуживания. От предыдущей системы она отличается тем, что заявка, поступившая в накопитель (очередь), может ожидать начала обслуживания лишь ограниченное время Т ож (чаще всего это случайная величина). Если её время Т ож истекло, то заявка покидает очередь и уходит из системы не обслуженной.

Математическое описание СМО

СМО рассматриваются как некоторые физические системы с дискретными состояниями х 0 , х 1 , …, х n , функционирующие при непрерывном времени t . Число состояний n может быть конечным или счетным (n → ∞). Система может переходить из одного состояния х i (i= 1, 2, … , n) в другое х j (j= 0, 1, … ,n) в произвольный момент времени t . Чтобы показать правила таких переходов, используют схему, называемую графом состояний . Для типов перечисленных выше систем графы состояний образуют цепь, в которой каждое состояние (кроме крайних) связано прямой и обратной связью с двумя соседними состояниями. Это схема гибели и размножения.
Переходы из состояния в состояние происходят в случайные моменты времени. Удобно считать, что эти переходы происходят в результате действия каких-то потоков (потоков входных заявок, отказов в обслуживании заявок, потока восстановления приборов и т.д.). Если все потоки простейшие, то протекающий в системе случайный процесс с дискретным состоянием и непрерывным временем будет марковским.
Поток событий - это последовательность однотипных событий, протекающих в случайные моменты времени. Его можно рассматривать как последовательность случайных моментов времени t 1 , t 2 , … появления событий.
Простейшим называют поток, обладающий следующими свойствами:

• Ординарность . События следуют по одиночке (противоположность потоку, где события следуют группами).
• Стационарность . Вероятность попадания заданного числа событий на интервал времени Т зависит только от длины интервала и не зависит от того, где на оси времени находиться этот интервал.
• Отсутствие последействия . Для двух непересекающихся интервалов времени τ 1 и τ 2 число событий, попадающих на один из них, не зависит от того, сколько событий попало на другой интервал.

В простейшем потоке интервалы времени Т 1 , Т 2 ,… между моментами t 1 , t 2 , … появления событий случайны, независимы между собой и имеют показательное распределение вероятностей f(t)=λe -λt , t≥0, λ=const, где λ - параметр показательного распределения, являющийся одновременно интенсивностью потока и представляющий собой среднее число событий, происходящих в единицу времени. Таким образом, t =M[T]=1/λ.
Марковские случайные события описываются обыкновенными дифференциальными уравнениями . Переменными в них служат вероятности состояний р 0 (t), p 1 (t),…,p n (t) .
Для очень больших моментов времени функционирования систем (теоретически при t → ∞) в простейших системах (системы, все потоки в которых – простейшие, а граф – схема гибели и размножения) наблюдается установившийся, или стационарный режим работы. В этом режиме система будет изменять свое состояние, но вероятности этих состояний (финальные вероятности ) р к , к= 1, 2 ,…, n, не зависят от времени и могут рассматриваться как среднее относительное время пребывания системы в соответствующем состоянии. Показатели эффективности СМО

• абсолютная и относительная пропускная способность системы;
• коэффициенты загрузки и простоя;
• среднее время полной загрузки системы;
• среднее время пребывания заявки в системе.

Показатели, характеризующие систему с точки зрения потребителей :

• P обс – вероятность обслуживания заявки,
• t сист – время пребывания заявки в системе.

Показатели, характеризующие систему с точки зрения её эксплуатационных свойств :

• λ b – абсолютная пропускная способность системы (среднее число обслуженных заявок в единицу времени),
• P обс – относительная пропускная способность системы,
• k з – коэффициент загрузки системы.

см. также Параметры экономической эффективности СМО

Задача . В вычислительный центр коллективного пользования с тремя ЭВМ поступают заказы от предприятий на вычислительные работы. Если работают все три ЭВМ, то вновь поступающий заказ не принимается, и предприятие вынуждено обратиться в другой вычислительный центр. Среднее время работы с одним заказом составляет 3 ч. Интенсивность потока заявок 0,25 (1/ч). Найти предельные вероятности состояний и показатели эффективности работы вычислительного центра.
Решение. По условию n=3, λ=0,25(1/ч), t об. =3 (ч). Интенсивность потока обслуживаний μ=1/t об. =1/3=0,33. Интенсивность нагрузки ЭВМ по формуле (24) ρ=0,25/0,33=0,75. Найдем предельные вероятности состояний:
по формуле (25) p 0 =(1+0,75+0,75 2 /2!+0,75 3 /3!) -1 =0,476;
по формуле (26) p 1 =0,75∙0,476=0,357; p 2 =(0,75 2 /2!)∙0,476=0,134; p 3 =(0,75 3 /3!)∙0,476=0,033 т.е. в стационарном режиме работы вычислительного центра в среднем 47,6% времени нет ни одной заявки, 35,7% - имеется одна заявка (занята одна ЭВМ), 13,4% - две заявки (две ЭВМ), 3,3% времени - три заявки (заняты три ЭВМ).
Вероятность отказа (когда заняты все три ЭВМ), таким образом, P отк. =p 3 =0,033.
По формуле (28) относительная пропускная способность центра Q = 1-0,033 = 0,967, т.е. в среднем из каждых 100 заявок вычислительный центр обслуживает 96,7 заявок.
По формуле (29) абсолютная пропускная способность центра A= 0,25∙0,967 = 0,242, т.е. в один час в среднем обслуживается 0,242 заявки.
По формуле (30) среднее число занятых ЭВМ k =0,242/0,33 = 0,725, т.е. каждая из трех ЭВМ будет занята обслуживанием заявок в среднем лишь на 72,5/3 =24,2%.
При оценке эффективности работы вычислительного центра необходимо сопоставить доходы от выполнения заявок с потерями от простоя дорогостоящих ЭВМ (с одной стороны, у нас высокая пропускная способность СМО, а с другой стороны - значительный простой каналов обслуживания) и выбрать компромиссное решение.

Задача . В порту имеется один причал для разгрузки судов. Интенсивность потока судов равна 0,4 (судов в сутки). Среднее время разгрузки одного судна составляет 2 суток. Предполагается, что очередь может быть неограниченной длины. Найти показатели эффективности работы причала, а также вероятность того, что ожидают разгрузки не более чем 2 судна.
Решение. Имеем ρ = λ/μ = μt об. =0,4∙2=0,8. Так как ρ = 0,8 < 1, то очередь на разгрузку не может бесконечно возрастать и предельные вероятности существуют. Найдем их.
Вероятность того, что причал свободен, по (33) p 0 = 1 - 0,8 = 0,2, а вероятность того, что он занят, P зан. = 1-0,2 = 0,8. По формуле (34) вероятности того, что у причала находятся 1, 2, 3 судна (т.е. ожидают разгрузки 0, 1, 2 судна), равны p 1 = 0,8(1-0,8) = 0,16; p 2 = 0,8 2 ∙(1-0,8) = 0,128; p 3 = 0,8 3 ∙(1-0,8) = 0,1024.
Вероятность того, что ожидают разгрузку не более чем 2 судна, равна
P=p 1 +p 2 +p 3 = 0,16 + 0,128 + 0,1024 = 0,3904
По формуле (40) среднее число судов, ожидающих разгрузки
L jч =0,8 2 /(1-0,8) = 3,2
а среднее время ожидания разгрузки по формуле (15.42)
T оч =3,2/0,8 = 4 сутки.
По формуле (36) среднее число судов, находящихся у причала, L сист. = 0,8/(1-0,8) = 4 (сутки) (или проще по (37) L сист. = 3,2+0,8 = 4 (сутки), а среднее время пребывания судна у причала по формуле (41) T сист = 4/0,8 = 5 (сутки).
Очевидно, что эффективность разгрузки судов невысокая. Для ее повышения необходимо уменьшение среднего времени разгрузки судна t об либо увеличение числа причалов n .

Задача . В универсаме к узлу расчета поступает поток покупателей с интенсивностью λ = 81 чел. в час. Средняя продолжительность обслуживания контролером-кассиром одного покупателя t об = 2 мин. Определить:
а. Минимальное количество контролеров-кассиров п min , при котором очередь не будет расти до бесконечности, и соответствующие характеристики обслуживания при n=n min .
б. Оптимальное количество n опт. контролеров-кассиров, при котором относительная величина затрат С отн., связанная с издержками на содержание каналов обслуживания и с пребыванием в очереди покупателей, задаваемая, например, как , будет минимальна, и сравнить характеристики обслуживания при n=n min и n=n опт.
в. Вероятность того, что в очереди будет не более трех покупателей.
Решение.
а. По условию l = 81(1/ч) = 81/60 = 1,35 (1/мин.). По формуле (24) r = l/ m = lt об = 1,35×2 = 2,7. Очередь не будет возрастать до бесконечности при условии r/n < 1, т.е. при n > r = 2,7. Таким образом, минимальное количество контролеров-кассиров n min = 3.
Найдем характеристики обслуживания СМО при п = 3.
Вероятность того, что в узле расчета отсутствуют покупатели, по формуле (45) p 0 = (1+2,7+2,7 2 /2!+2,7 3 /3!+2,7 4 /3!(3-2,7)) -1 = 0,025, т.е. в среднем 2,5% времени контролеры-кассиры будут простаивать.
Вероятность того, что в узле расчета будет очередь, по (48) P оч. = (2,7 4 /3!(3-2,7))0,025 = 0,735
Среднее число покупателей, находящихся в очереди, по (50) L оч. = (2,7 4 /3∙3!(1-2,7/3) 2)0,025 = 7,35.
Среднее время ожидания в очереди по (42) T оч. = 7,35/1,35 = 5,44 (мин).
Среднее число покупателей в узле расчета по (51) L сист. = 7,35+2,7 = 10,05.
Среднее время нахождения покупателей в узле расчета по (41) T сист. = 10,05/1,35 = 7,44 (мин).
Таблица 1

Характеристика обслуживания	Число контролеров-кассиров
	3	4	5	6	7
Вероятность простоя контролеров-кассиров p 0	0,025	0,057	0,065	0,067	0,067
Среднее число покупателей в очереди T оч.	5,44	0,60	0,15	0,03	0,01
Относительная величина затрат С отн.	18,54	4,77	4,14	4,53	5,22

Среднее число контролеров-кассиров, занятых обслуживанием покупателей, по (49) k = 2,7.
Коэффициент (доля) занятых обслуживанием контролеров-кассиров
= ρ/n = 2,7/3 = 0,9.
Абсолютная пропускная способность узла расчета А = 1,35 (1/мин), или 81 (1/ч), т.е. 81 покупатель в час.
Анализ характеристик обслуживания свидетельствует о значительной перегрузке узла расчета при наличии трех контролеров-кассиров.
б. Относительная величина затрат при n = 3
C отн. = = 3/1,35+3∙5,44 = 18,54.
Рассчитаем относительную величину затрат при других значениях п (табл. 1).
Как видно из табл. 2, минимальные затраты получены при n = n опт. = 5 контролерах-кассирах.
Определим характеристики обслуживания узла расчета при n = n опт. =5. Получим P оч. = 0,091; L оч. = 0,198; Т оч. = 0,146 (мин); L сист. = 2,90; T снст. = 2,15 (мин); k = 2,7; k 3 = 0,54.
Как видим, при n = 5 по сравнению с n = 3 существенно уменьшились вероятность возникновения очереди P оч. , длина очереди L оч. и среднее время пребывания в очереди T оч. и соответственно среднее число покупателей L сист. и среднее время нахождения в узле расчета T сист., а также доля занятых обслуживанием контролеров k 3. Но среднее число занятых обслуживанием контролеров-кассиров k и абсолютная пропускная способность узла расчета А естественно не изменились.
в. Вероятность того, что в очереди будет не более 3 покупателей, определится как
= 1- P оч. + p 5+1 + p 5+2 + p 5+3 , где каждое слагаемое найдем по формулам (45) – (48). Получим при n=5:

Заметим, что в случае n=3 контролеров-кассиров та же вероятность существенно меньше: P(r ≤ 3) =0,464.

Рассмотренный в предыдущей лекции марковский случайный процесс с дискретными состояниями и непрерывным временем имеет место в системах массового обслуживания (СМО).

Системы массового обслуживания – это такие системы, в которые в случайные моменты времени поступают заявки на обслуживание, при этом поступившие заявки обслуживаются с помощью имеющихся в распоряжении системы каналов обслуживания.

Примерами систем массового обслуживания могут служить:

• расчетно-кассовые узлы в банках, на предприятиях;
• персональные компьютеры, обслуживающие поступающие заявки или требования на решение тех или иных задач;
• станции технического обслуживания автомобилей; АЗС;
• аудиторские фирмы;
• отделы налоговых инспекций, занимающиеся приёмкой и проверкой текущей отчетности предприятий;
• телефонные станции и т. д.

	Узлы	Требования
Больница	Санитары	Пациенты
Производство
Аэропорт	Выходы на взлетно-посадочные полосы Пункты регистрации	Пассажиры

Рассмотрим схему работы СМО (рис. 1). Система состоит из генератора заявок, диспетчера и узла обслуживания, узла учета отказов (терминатора, уничтожителя заявок). Узел обслуживания в общем случае может иметь несколько каналов обслуживания.

Рис. 1

1. Генератор заявок – объект, порождающий заявки: улица, цех с установленными агрегатами. На вход поступает поток заявок (поток покупателей в магазин, поток сломавшихся агрегатов (машин, станков) на ремонт, поток посетителей в гардероб, поток машин на АЗС и т. д.).
2. Диспетчер – человек или устройство, которое знает, что делать с заявкой. Узел, регулирующий и направляющий заявки к каналам обслуживания. Диспетчер:

• принимает заявки;
• формирует очередь, если все каналы заняты;
• направляет их к каналам обслуживания, если есть свободные;
• дает заявкам отказ (по различным причинам);
• принимает информацию от узла обслуживания о свободных каналах;
• следит за временем работы системы.

1. Очередь – накопитель заявок. Очередь может отсутствовать.
2. Узел обслуживания состоит из конечного числа каналов обслуживания. Каждый канал имеет 3 состояния: свободен, занят, не работает. Если все каналы заняты, то можно придумать стратегию, кому передавать заявку.
3. Отказ от обслуживания наступает, если все каналы заняты (некоторые в том числе могут не работать).

Кроме этих основных элементов в СМО в некоторых источниках выделяются также следующие составляющие:

терминатор – уничтожитель трансактов;

склад – накопитель ресурсов и готовой продукции;

счет бухгалтерского учета – для выполнения операций типа «проводка»;

менеджер – распорядитель ресурсов;

Классификация СМО

Первое деление (по наличию очередей):

• СМО с отказами;
• СМО с очередью.

В СМО с отказами заявка, поступившая в момент, когда все каналы заняты, получает отказ, покидает СМО и в дальнейшем не обслуживается.

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

СМО с очередями подразделяются на разные виды в зависимости от того, как организована очередь, – ограничена или не ограничена . Ограничения могут касаться как длины очереди, так и времени ожидания, «дисциплины обслуживания».

Итак, например, рассматриваются следующие СМО:

• СМО с нетерпеливыми заявками (длина очереди и время обслуживания ограничено);
• СМО с обслуживанием с приоритетом, т. е. некоторые заявки обслуживаются вне очереди и т. д.

Типы ограничения очереди могут быть комбинированными.

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

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

Кроме этого СМО делятся на открытые СМО и замкнутые СМО.

В открытой СМО характеристики потока заявок не зависят от того, в каком состоянии сама СМО (сколько каналов занято). В замкнутой СМО – зависят. Например, если один рабочий обслуживает группу станков, время от времени требующих наладки, то интенсивность потока «требований» со стороны станков зависит от того, сколько их уже исправно и ждет наладки.

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

По количеству каналов СМО делятся на:

• одноканальные;
• многоканальные.

Характеристики системы массового обслуживания

Основными характеристиками системы массового обслуживания любого вида являются:

• входной поток поступающих требований или заявок на обслуживание;
• дисциплина очереди;
• механизм обслуживания.

Входной поток требований

Для описания входного потока требуется задать вероятностный закон, определяющий последовательность моментов поступления требований на обслуживание, и указать количество таких требований в каждом очередном поступлении. При этом, как правило, оперируют понятием «вероятностное распределение моментов поступления требований». Здесь могут поступать как единичные, так и групповые требования (количество таких требований в каждом очередном поступлении ). В последнем случае обычно речь идет о системе обслуживания с параллельно-групповым обслуживанием.

А i – время поступления между требованиями – независимые одинаково распределенные случайные величины;

E(A) – среднее (МО) время поступления;

λ=1/E(A) – интенсивность поступления требований;

Характеристики входного потока:

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

Дисциплина очереди

Очередь – совокупность требований, ожидающих обслуживания.

Очередь имеет имя.

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

• первым пришел – первый обслуживаешься;

first in first out (FIFO)

самый распространенный тип очереди.

Какая структура данных подойдет для описания такой очереди? Массив плох (ограничен). Можно использовать структуру типа СПИСОК.

Список имеет начало и конец. Список состоит из записей. Запись – это ячейка списка. Заявка поступает в конец списка, а выбирается на обслуживание из начала списка. Запись состоит из характеристики заявки и ссылки (указатель, за кем стоит). Кроме этого, если очередь с ограничением на время ожидания, то еще должно быть указано предельное время ожидания.

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

Действия со списком:

• вставить в хвост;
• взять из начала;
• удалить из списка по истечении времени ожидания.
• пришел последним - обслуживаешься первым LIFO (обойма для патронов, тупик на железнодорожной станции, зашел в набитый вагон).

Структура, известная как СТЕК. Может быть описан структурой массив или список;

• случайный отбор заявок;
• отбор заявок по критерию приоритетности.

Каждая заявка характеризуется помимо прочего уровнем приоритета и при поступлении помещается не в хвост очереди, а в конец своей приоритетной группы. Диспетчер осуществляет сортировку по приоритету.

Характеристики очереди

• ограничение времени ожидания момента наступления обслуживания (имеет место очередь с ограниченным временем ожидания обслуживания, что ассоциируется с понятием «допустимая длина очереди»);
• длина очереди.

Механизм обслуживания

Механизм обслуживания определяется характеристиками самой процедуры обслуживания и структурой обслуживающей системы. К характеристикам процедуры обслуживания относятся:

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

Для аналитического описания характеристик процедуры обслуживания оперируют понятием «вероятностное распределение времени обслуживания требований».

S i – время обслуживания i -го требования;

E(S) – среднее время обслуживания;

μ=1/E(S) – скорость обслуживания требований.

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

Коэффициент использования СМО

N ·μ – скорость обслуживания в системе, когда заняты все устройства обслуживания.

ρ=λ/(N μ) – называется коэффициентом использования СМО , показывает, насколько задействованы ресурсы системы.

Структура обслуживающей системы

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

Пример. Кассы в магазине.

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

Пример. Медицинская комиссия.

Комбинированное обслуживание – обслуживание вкладов в сберкассе: сначала контролер, потом кассир. Как правило, 2 контролера на одного кассира.

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

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

Основные критерии эффективности функционирования СМО

В качестве основных критериев эффективности функционирования систем массового обслуживания в зависимости от характера решаемой задачи могут выступать:

• вероятность немедленного обслуживания поступившей заявки (Р обсл =К обс /К пост);
• вероятность отказа в обслуживании поступившей заявки (P отк =К отк /К пост);

Очевидно, что Р обсл + P отк =1.

Потоки, задержки, обслуживание. Формула Поллачека–Хинчина

Задержка – один из критериев обслуживания СМО, время проведенное заявкой в ожидании обслуживания.

D i – задержка в очереди требования i ;

W i =D i +S i – время нахождения в системе требования i .

(с вероятностью 1) – установившаяся средняя задержка требования в очереди;

(с вероятностью 1) – установившееся среднее время нахождения требования в СМО (waiting).

Q(t) – число требований в очереди в момент времени t;

L(t) – число требований в системе в момент времени t (Q(t) плюс число требований, которые находятся на обслуживании в момент времени t.

Тогда показатели (если существуют)

(с вероятностью 1) – установившееся среднее по времени число требований в очереди;

(с вероятностью 1) – установившееся среднее по времени число требований в системе.

Заметим, что ρ<1 – обязательное условие существования d, w, Q и L в системе массового обслуживания.

Если вспомнить, что ρ= λ/(N μ), то видно, что если интенсивность поступления заявок больше, чем N μ, то ρ>1 и естественно, что система не сможет справиться с таким потоком заявок, а следовательно, нельзя говорить о величинах d, w, Q и L.

К наиболее общим и нужным результатам для систем массового обслуживания относятся уравнения сохранения

Следует обратить внимание, что упомянутые выше критерии оценки работы системы могут быть аналитически вычислены для систем массового обслуживания M/M/N (N >1), т. е. систем с Марковскими потоками заявок и обслуживания. Для М/G/ l при любом распределении G и для некоторых других систем. Вообще распределение времени между поступлениями, распределение времени обслуживания или обеих этих величин должно быть экспоненциальным (или разновидностью экспоненциального распределения Эрланга k-го порядка), чтобы аналитическое решение стало возможным.

Кроме этого можно также говорить о таких характеристиках, как:

• абсолютная пропускная способность системы – А=Р обсл *λ;
• относительная пропускная способность системы –

Еще один интересный (и наглядный) пример аналитического решения – вычисление установившейся средней задержки в очереди для системы массового обслуживания M/G/ 1 по формуле:

.

В России эта формула известна как формула Поллачека– Хинчина, за рубежом эта формула связывается с именем Росса (Ross).

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

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

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

Курсовая работа

«Имитационное моделирование системы массового обслуживания»

по курсу «Исследование операций»

Введение

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

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

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

– Абсолютная пропускная способность системы (А

Q

– вероятность отказа обслуживания заявки ();

k );

– среднее число заявок в очереди ();

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

Одним из методов расчета показателей эффективности СМО является метод имитационного моделирования. Практическое использование компьютерного имитационного моделирования предполагает построение соответствующей математической модели, учитывающей факторы неопределенности, динамические характеристики и весь комплекс взаимосвязей между элементами изучаемой системы. Имитационное моделирование работы системы начинается с некоторого конкретного начального состояния. Вследствие реализации различных событий случайного характера, модель системы переходит в последующие моменты времени в другие свои возможные состояния. Этот эволюционный процесс продолжается до конечного момента планового периода, т.е. до конечного момента моделирования.

1. Основные характеристики CМОи показатели их эффективности

1.1 Понятие марковского случайного процесса

Пусть имеется некоторая система, которая с течением времени изменяет свое состояние случайным образом. В этом случае говорят, что в системе протекает случайный процесс.

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

Процесс работы СМО – это случайный процесс с дискретными состояниями и непрерывным временем.

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

При анализе процессов работы СМО удобно пользоваться геометрической схемой – графом состояний . Обычно состояния системы изображаются прямоугольниками, а возможные переходы из состояния в состояние – стрелками. Пример графа состояний приведен на рис. 1.

Поток событий – последовательность однородных событий, следующих одно за другим в случайные моменты времени.

Поток характеризуется интенсивностью λ – частотой появления событий или средним числом событий, поступающих в СМО в единицу времени.

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

Поток событий называется стационарным, если его вероятностные характеристики не зависят от времени. В частности, интенсивность стационарного потока есть величина постоянная: .

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

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

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

1.2 Уравнения Колмогорова

Все переходы в системе из состояния в состояние происходят под некоторым потоком событий. Пусть система находится в некотором состоянии , из которого возможен переход в состояние , тогда можно считать, что на систему воздействует простейший поток с интенсивностью , переводящий ее из состояния в . Как только появляется первое событие потока, происходит ее переход . Для наглядности на графе состояний у каждой стрелки, соответствующей переходу, указывается интенсивность . Такой размеченный граф состояний позволяет построить математическую модель процесса, т.е. найти вероятности всех состояний как функции времени. Для них составляются дифференциальные уравнения, называемые уравнениями Колмогорова.

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

Например, для графа состояний, приведенного на рис. 1, уравнения Колмогорова имеют вид:

Т.к. в правой части системы каждое слагаемое входит 1 раз со знаком и 1 раз со знаком , то, складывая все уравнений, получим, что

,

,

Следовательно, одно из уравнений системы можно отбросить и заменить уравнением (1.2.1).

Чтобы получить конкретное решение надо знать начальные условия, т.е. значения вероятностей в начальный момент времени.

1.3 Финальные вероятности и граф состояний СМО

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

Смысл финальных вероятностей состоит в том, что они равны среднему относительному времени нахождения системы в данном состоянии.

Т.к. в стационарном состоянии производные по времени равны нулю, то уравнения для финальных вероятностей получаются из уравнений Колмогорова путем приравнивания нулю их правых частей.

Графы состояний, используемые в моделях систем массового обслуживания, называются схемой гибели и размножения. Такое название обусловлено тем, что эта схема используется в биологических задачах, связанных с изучением численности популяции. Его особенность состоит в том, что все состояния системы можно представить в виде цепочки, в которой каждое из состояний связано с предыдущим и последующим (рис 2).

Рис. 2. Граф состояний в моделях СМО

Предположим, что все потоки, переводящие систему из одного состояния в другое, простейшие. По графу, представленному на рис. 2, составим уравнения для финальных вероятностей системы. Они имеют вид:

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

,

.

Подставляя эти выражения в последнее уравнение системы, находим , затем находим остальные вероятности состояний СМО.

1.4 Показатели эффективности СМО

Цель моделирования СМО состоит в том, чтобы рассчитать показатели эффективности системы через ее характеристики. В качестве показателей эффективности СМО используются:

– абсолютная пропускная способность системы (А ), т.е. среднее число заявок, обслуживаемых в единицу времени;

– относительная пропускная способность (Q ), т.е. средняя доля поступивших заявок, обслуживаемых системой;

– вероятность отказа (), т.е. вероятность того, что заявка покинет СМО не обслуженной;

– среднее число занятых каналов (k );

– среднее число заявок в СМО ();

– среднее время пребывания заявки в системе ();

– среднее число заявок в очереди () – длина очереди;

– среднее число заявок в системе ();

– среднее время пребывания заявки в очереди ();

– среднее время пребывания заявки в системе ()

– степень загрузки канала (), т.е. вероятность того, что канал занят;

– среднее число заявок, обслуживаемых в единицу времени;

– среднее время ожидания обслуживания;

– вероятность того, что число заявок в очереди превысит определенное значение и т.п.

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

(1.4.1)

Формулы (1.4.1) и (1.4.2) называются формулами Литтла. Они вытекают из того, что в предельном стационарном режиме среднее число заявок, прибывающих в систему, равно среднему числу заявок, покидающих ее, т.е. оба потока заявок имеют одну и ту же интенсивность .

Формулы для вычисления показателей эффективности приведены в таб. 1.

Таблица 1.

Показатели	Одноканальная СМО с ограниченной очередью	Многоканальная СМО с ограниченной очередью
Финальные вероятности
Вероятность
Абсолютная пропускная способность
Относительная пропускная способность
Среднее число заявок в
Среднее число заявок под обслуживанием
Среднее число заявок в системе

1.5 Основные понятия имитационного моделирования

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

Компьютерное имитационное моделирование следует рассматривать как статический эксперимент.

Из теории функций случайных величин известно, что для моделирования случайной величины с любой непрерывной и монотонно возрастающей функцией распределения достаточно уметь моделировать случайную величину , равномерно распределенную на отрезке . Получив реализацию случайной величины , можно найти соответствующую ей реализацию случайной величины , так как они связаны равенством

Предположим, что в некоторой системе массового обслуживания время обслуживания одной заявки распределено по экспоненциальному закону с параметром , где – интенсивность потока обслуживания. Тогда функция распределения времени обслуживания имеет вид

Пусть - реализация случайной величины , равномерно распределенной на отрезке , а – соответствующая ей реализация случайного времени обслуживания одной заявки. Тогда, согласно (1.5.1)

1.6 Построение имитационных моделей

Первый этап создания любой имитационной модели – этап описания реально существующей системы в терминах характеристик основных событий. Эти события, как правило, связаны с переходами изучаемой системы из одного возможного состояния в другое и обозначаются как точки на временной оси. Для достижения основной цели моделирования достаточно наблюдать систему в моменты реализации основных событий.

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

Характеристики самого процесса массового обслуживания могут изменять свои значения либо в момент поступления новой заявки на обслуживание, либо при завершении обслуживания очередной заявки. К обслуживанию очередной заявки СМО может приступить немедленно (канал обслуживания свободен), но не исключена необходимость ожидания, когда заявке придется занять место в очереди (СМО с очередью, канал обслуживания занят). После завершения обслуживания очередной заявки СМО может сразу приступить к обслуживанию следующей заявки, если она есть, но может и простаивать, если таковая отсутствует. Необходимую информацию можно получить, наблюдая различные ситуации, возникающие при реализациях основных событий. Так, при поступлении заявки в СМО с очередью при занятом канале обслуживания длина очереди увеличивается на 1. Аналогично длина очереди уменьшается на 1, если завершено обслуживание очередной заявки и множество заявок в очереди не пусто.

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

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

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

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

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

1) увеличить длительность времени имитационного моделирования процесса функционирования изучаемой системы. В этом случае не только увеличивается вероятность достижения системой стационарного режима функционирования, но и возрастает число используемых псевдослучайных чисел, что также положительно влияет на качество получаемых результатов.

2) при фиксированной длительности времени Т имитационного моделирования провести N вычислительных экспериментов, называемых еще прогонами модели, с различными наборами псевдослучайных чисел, каждый из которых дает одно наблюдение. Все прогоны начинаются при одном и том же начальном состоянии моделируемой системы, но с использованием различных наборов псевдослучайных чисел. Преимуществом этого метода является независимость получаемых наблюдений , показателей эффективности системы. Если число N модели достаточно велико, то границы симметричного доверительного интервала для параметра определяются следующим образом:

, , т.е. , где

Исправленная дисперсия, ,

N – число прогонов программы, – надежность, .

2. Аналитическое моделирование СМО

2.1 Граф состояний системы и уравнения Колмогорова

Рассмотрим двухканальную систему массового обслуживания (n = 2) с ограниченной очередью равной шести (m = 4). В СМО поступает простейший поток заявок со средней интенсивностью λ = 4,8 и показательным законом распределения времени между поступлением заявок. Поток обслуживаемых в системе заявок является простейшим со средней интенсивностью μ = 2 и показательным законом распределения временем обслуживания.

Данная система имеет 7 состояний, обозначим их:

S 0 – система свободная, нет заявок;

S 1 – 1 заявка на обслуживании, очередь пуста;

S 2 – 2 заявки на обслуживании, очередь пуста;

S 3 – 2 заявки на обслуживании, 1 заявка в очереди;

S 4 – 2 заявки на обслуживании, 2 заявки в очереди;

S 5 – 2 заявки на обслуживании, 3 заявки в очереди;

S 6 – 2 заявки на обслуживании, 4 заявки в очереди;

Вероятности прихода системы в состояния S 0 , S 1 , S 2 , …, S 6 соответственно равны Р 0 , Р 1 , Р 2 , …, Р 6 .

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

Рис. 3. Граф состояний двухканальной СМО

Для построенного графа запишем уравнения Колмогорова:

Чтобы решить данную систему зададим начальные условия:

Систему уравнений Колмогорова (систему дифференциальных уравнений) решим численным методом Эйлера с помощью программного пакета Maple 11 (см. Приложение 1).

Метод Эйлера

где- в нашем случае, это правые части уравнений Колмогорова, n=6.

Выберем шаг по времени . Предположим , где Т – это время, за которое система выходит на стационарный режим. Отсюда получаем число шагов . Последовательно N раз вычисляя по формуле (1) получим зависимости вероятностей состояний системы от времени, приведенной на рис. 4.

Значения вероятностей СМО при равны:

Рис. 4. Зависимости вероятностей состояний системы от времени

P 0

P 5

P 4

P 3

P 2

P 1

2.2 Финальные вероятности системы

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

Т.к. в стационарном состоянии производные по времени равны 0, то уравнения для финальных вероятностей получаются из уравнений Колмогорова путем приравнивания правых частей 0. Запишем уравнения для финальных вероятностей для нашей СМО.

Решим данную систему линейных уравнений с помощью программного пакета Maple 11 (см. Приложение 1).

Получим финальные вероятности системы:

Сравнение вероятностей, полученных из системы уравнений Колмогорова при , с финальными вероятностями показывает, что ошибки равны:

Т.е. достаточно малы. Это подтверждает правильность полученных результатов.

2.3 Расчет показатели эффективности системы по финальным вероятностям

Найдем показатели эффективности системы массового обслуживания.

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

1) Вероятность отказав обслуживании заявки, т.е. вероятность того, что заявка покидает систему не обслуженной.В нашем случае заявке отказывается в обслуживании, если все 2 канала заняты, и очередь максимально заполнена (т.е. 4 человек в очереди), это соответствует состоянию системы S 6 . Т.к. вероятность прихода системы в состояние S 6 равна Р 6 , то

4) Средняя длина очереди, т.е. среднее число заявок в очереди, равна сумме произведений числа заявок в очереди на вероятность соответствующего состояния.

5) Среднее время пребывания заявки в очередиопределяется формулой Литтла:

3. Имитационное моделирование СМО

3.1 Алгоритм метода имитационного моделирования СМО (пошаговый подход)

Рассмотрим двухканальную систему массового обслуживания (n = 2) с максимальной длиной очереди равной шести (m = 4). В СМО поступает простейший поток заявок со средней интенсивностью λ = 4,8 и показательным законом распределения времени между поступлением заявок. Поток обслуживаемых в системе заявок является простейшим со средней интенсивностью μ = 2 и показательным законом распределения временем обслуживания.

Для имитации СМО воспользуемся одним из методов статистического моделирования – имитационным моделированием. Будем использовать пошаговый подход. Суть этого подхода в том, что состояния системы рассматриваются в последующие моменты времени, шаг между которыми является достаточно малым, чтобы за его время произошло не более одного события.

Выберем шаг по времени (). Он должен быть много меньше среднего времени поступления заявки () и среднего времени ее обслуживания (), т.е.

Где (3.1.1)

Исходя из условия (3.1.1) определим шаг по времени .

Время поступления заявки в СМО и время ее обслуживания являются случайными величинами. Поэтому, при имитационном моделировании СМО их вычисление производится с помощью случайных чисел.

Рассмотрим поступление заявки в СМО. Вероятность того, что на интервале в СМО поступит заявка, равна: . Сгенерируем случайное число , и, если , то будем считать, что заявка на данном шаге в систему поступила, если , то не поступила.

В программе это осуществляет isRequested () . Интервал времени примем постоянным и равным 0,0001, тогда отношение будет равно 10000. Если заявка поступила, то она принимает значение «истина», в противном случае значение «ложь».

bool isRequested()

double r = R. NextDouble();

if (r < (timeStep * lambda))

Рассмотрим теперь обслуживание заявки в СМО. Время обслуживания заявки в системе определяется выражением , где – случайное число. В программе время обслуживания определяется с помощью функции GetServiceTime () .

double GetServiceTime()

double r = R. NextDouble();

return (-1/mu*Math. Log (1-r, Math.E));

Алгоритм метода имитационного моделирования можно сформулировать следующим образом. Время работы СМО (Т ) разбивается на шаги по времени dt , на каждом из них выполняется ряд действий. Вначале определяются состояния системы (занятость каналов, длина очереди), затем, с помощью функции isRequested () , определяется, поступила ли на данном шаге заявка или нет.

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

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

По истечении времени Т , т.е., после моделирования работы СМО, вычисляются показатели эффективности работы системы и результаты выводятся на экран.

3.2 Блок-схема программы

Блок-схема программы, реализующей описанный алгоритм, приведена на рис. 5.

Рис. 5. Блок-схема программы

Распишем некоторые блоки более подробно.

Блок 1. Задание начальных значений параметров.

Random R; // Генератор случайных чисел

public uint maxQueueLength; // Максимальная длина очереди

public uint channelCount; // Число каналов в системе

public double lambda; // Интенсивность потока поступления заявок

public double mu; // Интенсивность потока обслуживания заявок

public double timeStep; // Шагповремени

public double timeOfFinishProcessingReq; // Время окончания обслуживания заявки во всех каналах

public double timeInQueue; // Время пребывания СМО в состояниях с очередью

public double processingTime; // Времяработысистемы

public double totalProcessingTime; // Суммарноевремяобслуживаниязаявок

public uint requestEntryCount; // Числопоступившихзаявок

public uint declinedRequestCount; // Числоотказанныхзаявок

public uint acceptedRequestCount; // Числообслуженныхзаявок

uint queueLength; // Длина очереди //

Тип, описывающий состояния СМО

enum SysCondition {S0, S1, S2, S3, S4, S5, S6};

SysCondition currentSystemCondition; // Текущее состояние системы

Задание состояний системы. Выделим у данной 2-х канальной системы 7 различных состояний: S 0 , S 1 . S 6 . СМО находится в состоянии S 0 , когда система свободна; S 1 – хотя бы один канал свободен; в состоянии S 2 , когда все каналы заняты, и есть место в очереди; в состоянии S 6 – все каналы заняты, и очередь достигла максимальной длины (queueLength = 4).

Определяем текущее состояние системы с помощью функции GetCondition()

SysCondition GetCondition()

SysCondition p_currentCondit = SysCondition.S0;

int busyChannelCount = 0;

for (int i = 0; i < channelCount; i++)

if (timeOfFinishProcessingReq[i] > 0)

busyChannelCount++;

p_currentCondit += k * (i + 1);

if (busyChannelCount > 1)

{p_currentCondit ++;}

return p_currentCondit + (int) QueueLength;

Изменение времени пребывания СМО в состояниях с длиной очереди 1, 2,3,4. Это реализуется следующим программным кодом:

if (queueLength > 0)

timeInQueue += timeStep;

if (queueLength > 1)

{timeInQueue += timeStep;}

Присутствует такая операция, как помещение заявки на обслуживание в свободный канал. Просматриваются, начиная с первого, все каналы, когда выполняется условие timeOfFinishProcessingReq [ i ] <= 0 (канал свободен), в него подается заявка, т.е. генерируется время окончания обслуживания заявки.

for (int i = 0; i < channelCount; i++)

if (timeOfFinishProcessingReq [i] <= 0)

timeOfFinishProcessingReq [i] = GetServiceTime();

totalProcessingTime+= timeOfFinishProcessingReq [i];

Обслуживаниезаявоквканалахмоделируетсякодом:

for (int i = 0; i < channelCount; i++)

if (timeOfFinishProcessingReq [i] > 0)

timeOfFinishProcessingReq [i] -= timeStep;

Алгоритм метода имитационного моделирования реализован на языке программирования C#.

3.3 Расчет показателей эффективности СМО на основе результатов ее имитационного моделирования

Наиболее важными являются такие показатели, как:

1) Вероятность отказа в обслуживании заявки, т.е. вероятность того, что заявка покидает систему не обслуженной.В нашем случае заявке отказывается в обслуживании, если все 2 канала заняты, и очередь максимально заполнена (т.е. 4 человек в очереди). Для нахождения вероятности отказа разделим время пребывания СМО в состоянии с очередью 4 на общее время работы системы.

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

3) Абсолютная пропускная способность– это среднее число заявок, обслуживаемых в единицу времени.

4) Длина очереди, т.е. среднее число заявок в очереди. Длина очереди равна сумме произведений числа человек в очереди на вероятность соответствующего состояния. Вероятности состояний найдем как отношение времени нахождения СМО в этом состоянии к общему времени работы системы.

5) Среднее время пребывания заявки в очереди определяется формулой Литтла

6) Среднее число занятых каналовопределяется следующим образом:

7) Процент заявок, которым было отказано в обслуживании, находится по формуле

8) Процент обслуженных заявок находится по формуле

3.4 Статистическая обработка результатов и их сравнение с результатами аналитического моделирования

Т.к. показатели эффективности получаются в результате моделирования СМО в течение конечного времени, они содержат случайную компоненту. Поэтому, для получения более надежных результатов нужно провести их статистическую обработку. С этой целью оценим доверительный интервал для них по результатам 20 прогонов программы.

Величина попадает в доверительный интервал, если выполняется неравенство

, где

математическое ожидание (среднее значение), находится по формуле

Исправленная дисперсия,

,

N =20 – число прогонов,

– надежность. При и N =20 .

Результат работы программы представлен на рис. 6.

Рис. 6. Вид программы

Для удобства сравнения результатов, полученных различными методами моделирования, представим их в виде таблицы.

Таблица 2.

Показатели эффективности СМО	Результаты аналитического моделирования	Результаты имитационного моделирования (послед. шаг)	Результаты имитационного моделирования
			Нижняя граница доверительного интервала	Верхняя граница доверительного интервала
Вероятность отказа		0,174698253017626	0,158495148639101	0,246483801571923
Относительная пропускная способность		0,825301746982374	0,753516198428077	0,841504851360899
Абсолютная пропускная способность		3,96144838551539	3,61687775245477	4,03922328653232
Средняя длина очереди		1,68655313447018	1,62655862750852	2,10148609204869
Среднее время пребывания заявки в очереди	0,4242558575	0,351365236347954	0,338866380730942	0,437809602510145
Среднее число занятых каналов		1,9807241927577	1,80843887622738	2,01961164326616

Из табл. 2 видно, что результаты, полученные при аналитическом моделировании СМО, попадают в доверительный интервал, полученный по результатам имитационного моделирования. Т.е., результаты, полученные разными методами, согласуются.

Заключение

В данной работе рассмотрены основные методы моделирования СМО и расчета показателей их эффективности.

Проведено моделирование двухканальной СМО с максимальной длиной очереди равной 4 с помощью уравнений Колмогорова, а также, найдены финальные вероятности состояний системы. Рассчитаны показатели ее эффективности.

Проведено имитационное моделирование работы такой СМО. На языке программирования C# составлена программа, имитирующая ее работу. Проведена серия расчетов, по результатам которых найдены значения показателей эффективности системы и выполнена их статистическая обработка.

Полученные при имитационном моделировании результаты согласуются с результатами аналитического моделирования.

Литература

1. Вентцель Е.С. Исследование операций. – М.: Дрофа, 2004. – 208 с.

2. Волков И.К., Загоруйко Е.А. Исследование операций. – М.: Изд.-во МГТУ им. Н.Э. Баумана, 2002. – 435 с.

3. Волков И.К., Зуев С.М., Цветкова Г.М. Случайные процессы. – М.: Изд.-во МГТУ им. Н.Э. Баумана, 2000. – 447 с.

4. Гмурман В.Е. Руководство к решению задач по теории вероятностей и математической статистике. – М.: Высшая школа, 1979. – 400 с.

5. Ивницкий В.Л. Теория сетей массового обслуживания. – М.: Физматлит, 2004. – 772 с.

6. Исследование операций в экономике/ под ред. Н.Ш. Кремера. – М.: Юнити, 2004. – 407 с.

7. Таха Х.А. Введение в исследование операций. – М.: ИД «Вильямс», 2005. – 902 с.

8. Харин Ю.С., Малюгин В.И., Кирлица В.П. и др. Основы имитационного и статистического моделирования. – Минск: Дизайн ПРО, 1997. – 288 с.