Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

КМ-АХД-краткий конспект лекций

.pdf
Скачиваний:
50
Добавлен:
11.03.2016
Размер:
861.53 Кб
Скачать

Количественные метолы АХД

стр. 31 из 55

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

10.Литература

11.Основные понятия теории массового обслуживания

12.Показатели эффективности системы массового обслуживания (СМО)

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

14.Моделирование СМО с помощью графов состояний. Процессы «рождения-гибели»

15.СМО с отказами в обслуживании

16.СМО с очередью

17.СМО с ограничением на длину очереди

10. Литература

Использованная при подготовке темы:

3.Таха Хамди. Введение в исследование операций. М.: Мир, 1985

4.Фомин Г.П. Системы и модели массового обслуживания в коммерческой деятельности: Учеб. пособие. – М.: Финансы и статистика, 2000. -144 с.: ил.

Рекомендуется для самостоятельного обучения:

2.Лабскер Л.Г., Бабешко Л.О. Теория массового обслуживания в экономической сфере: учеб. пособие для вузов. М.: Банки и биржи, ЮНИТИ, 1998. – 319 с.

3.Исследование операций в экономике: Учебн. пособие для вузов/ Н.Ш.Кремер, Б.А.

Путко и др.; Под ред. проф. Н.Ш.Кремера. – М.: ЮНИТИ, 2000. – 407 с.

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

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

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

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

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

Количественные метолы АХД стр. 32 из 55

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

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

Основоположником теории массового обслуживания считается известный датский ученый А.К.Эрланг. Являясь сотрудником Копенгагенской телефонной компании, он опубликовал в 1909 г. работу «Теория вероятностей и телефонные переговоры», ы

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

Большой вклад в развитие теории массового обслуживания внесли и советские математики А.Я. Хинчин, Б.В. Гнезденко, А.Н.Колмогоков и др.

Каждая СМО включает в свою структуру некоторое число обслуживающих устройств,

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

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

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

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

Схема СМО

Количественные метолы АХД

 

 

стр. 33 из 55

 

 

 

 

 

 

 

 

 

 

 

Каналы

 

 

 

 

 

 

обслуживания

 

 

Входящий

 

 

 

Поток

поток

1

 

заявок

 

 

 

 

обслуженных

Очередь

 

 

 

 

 

 

 

 

 

заявок

 

 

 

 

 

 

 

2

Вход

Поток

необслуженных

заявок

Выход

 

n

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

4. Показатели эффективности системы массового обслуживания (СМО)

В качестве характеристик эффективности функционирования СМО используются 3

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

ееиспользования, качества работы СМО, доходности СМО:

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

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

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

3.Показатели доходности СМО: средний доход СМО, приносимый в единицу времени, прибыль СМО и пр.

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

в зависимости от характера потоков – на марковские и немарковские;

Количественные метолы АХД

стр. 34 из 55

по числу каналов – на одноканальные и многоканальные;

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

по дисциплине очереди – ПЕРПО, ПОСПО, СОЗ, КП;

по ограничению потока заявок – замкнутые и разомкнутые (открытые);

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

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

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

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

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

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

 

Λ01

S0

S1

Λ10

Система может находиться в одном из двух состояний: S0 – канал свободен и S1

канал занят. Переход системы из состояния S0 в состояние S1 происходит под воздействием пуассоновского (простейшего) потока заявок интенсивностью 01. Из состояния S1 в состояние S0 систему переводит уже поток обслуживания интенсивностью

10.

В связи с вероятностным характером состояний СМО можно сказать, что pi(t) – это вероятность того, что система будет находиться в состоянии Si в момент времени t. Такая вероятность определяется числом поступивших заявок на обслуживание (0 или 1 в

примере).

Случайный процесс, происходящий в системе, заключается в том, что в случайные моменты времени t0, t1,…,tn система оказывается в том или ином заранее известном дискретном состоянии последовательно.

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

Рассмотрим немногим более сложную схему СМО с состояниями S0, S1 и S2.

Предполагается, что происходящий процесс непрерывен во времени.

Количественные метолы АХД

 

 

стр. 35 из 55

 

 

Λ01

 

Λ12

 

 

S0

 

S1

 

S2

 

 

Λ10

 

Λ21

 

 

 

 

 

Для любого момента времени t справедливо записать нормировочное уравнение,

показывающее, что сумма вероятностей всех возможных состояний системы равна 1:

2

pi (t) p0 (t) p1 (t) p2 (t) 1

i 0

Проведем анализ системы в момент времени t, задав малое приращение времениt, и найдем вероятность p1(t+ t) того, что система в момент времени t+ t будет находится в состоянии S1, которое достигается несколькими возможными путями:

1. система в момент времени t с вероятностью p1(t) находилась в состоянии S1 и за малое приращение времени t не успела перейти в другое состояние (ни в S0, ни в S2).

Вывести систему

из

состояния S1

можно суммарным

простейшим

потоком

с

интенсивностью

( 10

+ 12). Тогда

вероятность выхода

из состояния

S1 за малый

промежуток времени t приближенно равна ( 10 + 12)* t, а вероятность невыхода из этого состояния равна (1-( 10 + 12)* t). В соответствии с этим вероятность того, что система останется в состоянии S1, на основании теоремы умножения вероятностей равна:

p1(t) * (1-( 10 + 12)* t);

 

2. система находилась в соседнем состоянии S0

и за малое время t перешла в

состояние S1. Переход произошел под воздействием

потока 01 с вероятностью,

приближенно равной 01* t. Вероятность того, что система будет находиться в состоянии

S1, в этом варианте равна p0(t)* 01* t.

 

3. система находилась в соседнем состоянии S2

и за малое время t перешла в

состояние S1. Переход произошел под воздействием

потока 21 с вероятностью,

приближенно равной 21* t. Вероятность того, что система будет находиться в состоянии

S1, в этом варианте равна p2(t)* 21* t.

Применим теперь теорему сложения вероятностей для этих вариантов: p1(t+ t) = p1(t) * (1-( 10 + 12)* t) + p0(t)* 01* t + p2(t)* 21* t

которое можно записать иначе:

p1 (t t) p1 (t) p0 (t) 01 p2 (t) 21 p1 (t)( 10 12 )

t

Переходя к пределу при t 0, приближенные равенства перейдут в точные,

получаем производные первого порядка

Количественные метолы АХД

стр. 36 из 55

dpdt1 p0 01 p2 21 p1 ( 10 12 )

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

dpo p1 10 po 01

dt

dpdt1 p0 01 p2 21 p1 ( 10 12 )

dpdt2 p1 12 p2 21

Уравнения Колмогорова позволяют вычислить все вероятности состояний СМО Si

в функции времени pi(t). Из теории случайных процессов известно, что если число состояний системы конечно, а из каждого состояния можно перейти в любое другое состояние, то существуют предельные (финальные) вероятности состояний, которые дают среднюю относительную величину времени пребывания системы в этом состоянии. Если предельная вероятность состояния S0 равна p0=0,2, то в среднем 20% времени система находится в состоянии S0.

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

умноженная на суммарную интенсивность всех потоков, выводящих (выходящие стрелки)

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

n

pi 1

i 1

Составим систему предельных вероятностей для СМО с тремя состояниями:

Состояния

Входящие

 

Выходящие

 

 

 

 

S0

p0* 01

=

p1* 10

 

 

 

 

S1

p1*( 10+ 12)

=

p0* 01+p2* 21

 

 

 

 

S2

p2* 21

=

p1* 12

 

 

 

 

Количественные метолы АХД

 

стр. 37 из 55

 

 

 

 

 

p0+p1+p2

=

1

 

 

 

 

Среди однородных марковских процессов существует класс случайных процессов,

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

 

Λ0

 

Λ1

Λn-1

 

S0

S1

Sn

 

 

 

 

μ0

 

μ1

μn-1

 

 

 

 

Величина k отображает интенсивность рождения нового представителя некоторой популяции, причем текущий объем популяции равен k; величина является интенсивностью гибели одного представителя этой популяции, если текущий объем популяции равен k.

Найдем финальное распределение для такого процесса и составим уравнения для каждого состояния (S1, S2,…, Sk,…, Sn):

Для S0: 0p0= 0p1

S1: ( 1+ 0)p1= 0p0+ 1p2 или после преобразований 1p1= 1p2

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

0p0= 0p11p1= 1p2

……………….

n-1pn-1= n-1pn

p0+p1+…+pn=1

Решая эту систему уравнений, можно получить выражения, определяющие финальные состояния системы массового обслуживания:

p

 

(1

0

 

0

1

...

0 1... n 1

) 1

0

 

 

 

 

 

 

0

 

0

1

 

0 1... n 1

 

 

 

 

 

p

0

p

;

p

 

 

0 1

p

; ...; p

 

 

0 1... n 1

 

2

 

n

 

1

0

0

 

 

 

10

o

 

 

0

1... n 1

 

 

 

 

 

 

 

 

 

 

Пример:

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

Размеченный граф их работы выглядит следующим образом:

Количественные метолы АХД

 

 

стр. 38 из 55

 

 

Λ0

 

Λ1

 

 

S0

 

S1

 

S2

 

 

μ0

 

μ1

 

 

 

 

 

S0 – обе кассы свободны;

S1 – одна касса занята (любая из двух);

S2 – обе кассы заняты.

Найдем предельные вероятности p0, p1, p2 при следующих исходных данных:

0 =1 пок/мин; 1 =0,8 пок/мин; 1= 2=0,6 пок/мин.

0p0= 0p1

p0=0,6p1

p0= 0,204

1p1= 1p2

0,8p1=0,6p2

p1=0,341

p0+p1+p2=1

p0+p1+p2=1

p2=0,455

Следовательно, доля времени простоя касс составляет 20,4% от всего рабочего времени, и т.д.

6.СМО с отказами в обслуживании Рассмотрим одноканальную СМО с отказами в обслуживании.

На вход системы поступает пуассоновский поток заявок с интенсивностью , а

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

 

 

Λ

 

 

S0

 

 

 

S1

 

 

 

 

 

μ

 

 

 

Примерами одноканальных СМО с отказами в обслуживании являются:

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

Запишем систему дифференциальных уравнений Колмогорова:

dpo (t) p1 (t) po (t) dt

dp1 (t) p1 (t) p0 (t) dt

p0(t) + p1(t) = 1

Отсюда дифференциальное уравнение для определения вероятности p0(t) состояния

S0:

Количественные метолы АХД

стр. 39 из 55

dp0 (t) ( ) p0 (t) dt

Пусть в момент времени t=0 система находилась в состоянии S0, тогда p0(0) = 1, p1(0)=0. В этом случае решение дифференциального уравнения позволит определить вероятность того, что канал свободен и не занят обслуживанием:

 

(t)

 

 

 

( )t

 

 

 

p0

 

 

e

 

 

 

 

 

 

p (t) 1 p

 

(t)

 

(1

e ( )t ) - вероятность занятости канала.

0

 

 

1

 

 

 

 

 

 

 

 

 

С течением времени вероятность того, что канал свободен уменьшается и в пределе при t= стремится к величине:

p0

 

 

, тог да как p1

 

 

 

 

 

 

 

 

 

 

Функции p0(t) и p1(t) определяют переходный процесс в одноканальной СМО и описывают процесс экспоненциального приближения СМО к своему предельному

состоянию с постоянной времени

1

, характерной для рассматриваемой системы. С

 

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

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

Абсолютная пропускная способность, определяющая число заявок, обслуживаемых в единицу времени в пределе при t= , равна

A *

Соответственно доля заявок, получивших отказ при этих предельных условиях:

pотк p1

 

, а общее число необслуженных заявок равно

2

.

 

 

 

 

 

Пример. Статистическими исследованиями в результате наблюдения установлено,

что интенсивность потока телефонных звонков коммерческому директору фирмы =1,2

вызова в минуту, средняя продолжительность разговора (обслуживания заявки) tобсл =2,5

мин и все потоки событий имеют характер простейших пуассоновских потоков.

Количественные метолы АХД стр. 40 из 55

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

течение 1 часа работы СМО.

Размеченный граф состояний этой СМО приведен выше.

Начнем с определения интенсивности потока обслуживания (заявок в минуту) –

 

1

0,4заявок / мин

 

 

 

2,5

 

 

 

 

 

 

 

 

Постоянная времени

1

0,63

мин, тогда переходный процесс в СМО

 

 

 

 

 

 

 

завершается через 3 =1,9 мин;

относительная пропускная способность или вероятность того, что канал свободен -

p0

Q

 

0,25;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Абсолютная пропускная способность – А= Q=0,3 заявки в минуту;

 

Вероятность занятости канала -

p1 1 p0

 

0,75.

 

 

 

 

 

 

 

 

 

 

 

Число

заявок, обслуженных

в течение

часа составляет: 60А=18 заявок, а

получивших отказ – 60* *p1=54 заявки.

Сравним фактическую пропускную способность СМО с номинальной, т.е.

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

перерыва.

 

 

 

Номинальная пропускная способность равна

Q

60

24заявки в час , т.е.

 

 

н

2,5

 

 

 

 

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

происходящего в СМО, составляет только 75%=18:24*100% от номинальной.