Главная
База знаний "Allbest"
Математика
Стационарное распределение в системах массового обслуживания с разнотипными заявками и ограниченным числом требований
Стационарное распределение в системах массового обслуживания с разнотипными заявками и ограниченным числом требований
Стационарное распределение вероятностей. Построение математических моделей, графов переходов. Получение уравнения равновесия систем массового обслуживания с различным числом приборов, требованиями различных типов и ограниченными очередями на приборах.
| Рубрика | Математика |
| Вид | дипломная работа |
| Язык | русский |
| Дата добавления | 23.12.2012 |
| Размер файла | 2,4 M |
Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже
Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.
Размещено на http://www.allbest.ru/
Министерство образования Республики Беларусь
Учреждение образования
"Гомельский государственный университет
имени Франциска Скорины"
математический факультет
Кафедра экономической кибернетики и теории вероятности
Дипломная работа
Стационарное распределение в СМО с разнотипными заявками и ограниченным числом требований
Исполнитель:
студент группы М-51 Колядова Ю.В.
Научный руководитель:
к. физ. - мат. н., доцент Бураковский В.В.
Рецензент
к. физ. - мат. н., старший преподаватель Вересович П.П.
Гомель 2008
Реферат
Дипломная работа _____ страницы, 5 источников, 4 рисунка.
Ключевые слова:
Система массового обслуживания (СМО), математическая модель СМО, уравнения равновесия, заявка, прибор, очередь.
Объект исследования: стационарное распределение вероятностей в СМО с различным числом приборов, требованиями различных типов и ограниченными очередями на приборах.
Цель работы: исследование стационарного распределения вероятностей в СМО с различным числом приборов, требованиями различных типов и ограниченными очередями на приборах.
Метод исследования: теория массового обслуживания.
Полученные результаты: в данной работе было исследовано стационарное распределение вероятностей, построены математические модели, графы переходов, получены уравнения равновесия СМО с различным числом приборов, требованиями различных типов и ограниченными очередями на приборах.
Область применения: системы массового обслуживания.
Содержание
- Введение
- 1 Стационарное распределение вероятностей состояний СМО с двумя приборами, разнотипными требованиями и ограниченной очередью
- 1.1 Построение математической модели СМО с двумя приборами, разнотипными требованиями и ограниченной очередью
- 1.2 Построение уравнений равновесия СМО с двумя приборами, заявками двух типов и ограниченной очередью
- 1.3 Построение уравнений равновесия СМО с двумя приборами, заявками n типов и ограниченной очередью
- 2. Исследование стационарного распределения вероятностей состояний СМО с несколькими приборами и ограниченными очередями на приборах
- 2.1 Построение математической модели СМО с N приборами, разнотипными требованиями n типов, у которой на i-том приборе могут находиться не более mi требований (i=)
- 2.2 Построение уравнений равновесия для частного случая СМО с N=2 приборами, заявками n=3 типов, у которой на i-том приборе могут находиться не более mi=1 требований (i=1,2,3)
- 2.3 Построение математической модели СМО с N приборами, разнотипными требованиями n типов, которых может находиться в СМО не более m
- 2.4 Построение уравнений равновесия для частного случая СМО с N=2 приборами, заявками n=3 типов, которых может находиться в СМО не более m=2
- 3. Стационарное распределение СМО с N приборами, заявками n типов, у которой на каждом приборе имеется буфер ёмкости M, причём заявки i-того типа обслуживаются на i-том приборе
- 3.1 Построение математической модели СМО с N=2 приборами, заявками N=2 типов, у которой на каждом приборе имеется буфер ёмкости m=2, причём заявки i-того типа обслуживаются на i-том приборе
- 3.2 Построение уравнений равновесия для частного случая СМО с N=2 приборами, заявками N=2 типов, у которой на каждом приборе имеется буфер ёмкости m=2, причём заявки i-того типа обслуживаются на i-том приборе
- 3.3 Построение математической модели СМО с N=3 приборами, заявками N=3 типов, у которой на каждом приборе имеется буфер ёмкости m=2, причём заявки i-того типа обслуживаются на i-том приборе
- 3.4 Построение уравнений равновесия для частного случая СМО с N=3 приборами, заявками N=3 типов, у которой на каждом приборе имеется буфер ёмкости m=2, причём заявки i-того типа обслуживаются на i-том приборе
- 3.5 Построение математической модели СМО с N приборами, заявками N типов, у которой на каждом приборе имеется буфер ёмкости m, причём заявки i-того типа обслуживаются на i-том приборе
- 3.6 Построение уравнений равновесия для СМО с N приборами, заявками N типов, у которой на каждом приборе имеется буфер ёмкости m, причём заявки i-того типа обслуживаются на i-том приборе
- 4. Построение векторно-матричных уравнений СМО с двумя приборами, требованиями трех типов и конечными буферами
- Заключение
- Список использованных источников
Введение
Практические требования телефонного дела, физики и рациональной организации массового обслуживания (билетные кассы, магазины и т.д.) выдвинули в начале ХХ столетия ряд интересных математических задач нового типа. Первоначально эти задачи касались преимущественно вопросов обслуживания абонентов телефонной станции, расчета запасов магазинов для бесперебойного снабжения покупателей, а также установления наиболее рационального числа продавцов и касс в торговых предприятиях. На первичное развитие этой теории особое влияние оказали работы известного датского ученого А.К. Эрланга (1878 - 1929) - многолетнего сотрудника Копенгагенской телефонной компании. С того времени интерес к проблемам, выдвинутым Эрлангом, необычайно возрос. Оказалось, что задачи типа телефонных возникают в самых разнообразных направлениях исследований: в естествознании, в технике, экономике, транспорте, военном деле, организации производства.
Рассмотрим несколько областей применения теории массового обслуживания.
Предположим, что на телефонную станцию в случайном порядке поступают вызовы. Если в момент поступления вызова на станции имеются свободные линии, то происходит подключение абонента к свободной линии и начинается разговор в течение того времени, которое необходимо для его завершения. Если же на станции все линии заняты, то возможны различные системы обслуживания абонентов. В настоящее время особенно хорошо разработаны две системы обслуживания: система с ожиданием и система с потерями.
Ситуация, которая создается около театральной кассы, когда в нее обращаются за билетами, весьма напоминает описание системы обслуживания абонентов на телефонной станции. Если только в первоначальной постановке задачи шла речь о телефонных линиях, то теперь вопрос касается занятости кассира. Стремление рационально обслуживать потребителей приводит к необходимости изучения закономерностей образования очередей. Знание этих закономерностей должно, в частности, помочь решению вопроса о числе касс, которые рационально установить для продажи билетов на железнодорожной станции или в магазине. Содержание каждой кассы вызывает некоторые расходы, но и потеря требований также наносит определенный ущерб. Возникает вопрос о разыскании некоторого оптимума.
Для многих реальных задач научного, производственного и экономического характера естественны не только задачи, в которых рассматриваются обслуживание с потерями и ожидание без ограничения времени. Зачастую мы отказываемся от обслуживания только из-за возможной длительной задержки с началом обслуживания. Так, если в очереди к продавцу имеется более пяти покупателей, то порой мы уходим из магазина и откладываем предполагаемую покупку. Несколько иная ситуация может создаваться, когда ограничено не время ожидания, а время пребывания в системе обслуживания. С такой постановкой задачи приходится встречаться при продаже скоропортящихся продуктов: от момента изготовления до употребления должно пройти не более чем единиц времени, т.к. иначе эти продукты теряют свои ценные качества и могут представлять угрозу для здоровья потребителя.
Таким образом, естественна постановка следующей группы близких по характеру задач. Требования, поступающие для обслуживания, остаются в очереди, если число ранее прибывших и ожидающих обслуживания требований не превосходит заданного числа k, в противном случае требование теряется. Каждое требование остается в системе не более чем время , даже если началось его обслуживание. Возможна также и третья постановка вопроса: ограничено время ожидания величиной , но если до истечения этого срока обслуживание началось, то оно доводиться до конца. Во всех трех случаях особый интерес представляет вычисление среднего числа потерь за определенный промежуток времени, среднего времени ожидания начала обслуживания или потерянного времени на ожидание. Во второй постановке задачи естественно различать среднее число потерь не совсем обслуженных требований, обслуживание которых все-таки было начато.
В поставленных задачах мы исходили из предположения, что обслуживающие приборы обладают абсолютной надежностью, и сами никогда не выходят из рабочего состояния. Такое положение несколько идеализирует реальные системы. В результате возникает естественная и важная задача учета влияния на эффективность системы обслуживания порчи обслуживающих приборов. Изучение этой задачи началось сравнительно недавно. Возможное разнообразие практически интересных вопросов здесь совершенно не ограничено.
Можно указать множество других постановок задач реального содержания, которые в своей математической части сводится к вопросам теории массового обслуживания. Задача математической теории состоит в первую очередь в выработке общих методов, применимых не только к решению тех частных задач, на базе которых была начата ее разработка, но и множества других, быть может, даже очень далеких по своей формулировке от первоначальных.
В подавляющем большинстве работ по теории массового обслуживания, как тех из них, которые послужили базой построения теории, так и современных, рассматривается простейший случай потоков, когда вероятность поступления в промежуток времени t ровно k требований задается формулой
где >0 - среднее число заявок, поступающих в единицу времени. Поступающий поток при этом считается таким, что для любой конечной группы непересекающихся отрезков времени числа появившихся на их протяжении требований представляют собой взаимно независимые случайные величины.
Стационарность потока означает, что для любых попарно непересекающихся промежутков времени 1,2,.,n для любого n=1,2,. вероятность поступления в этих промежутках k1,k2,.,kn заявок соответственно зависит только от этих чисел и длин промежутков 1, 2,., n, и не зависит от их взаимного расположения. В частности, вероятность появления k заявок в промежутке времени [Т, Т+t) не зависит от Т и является функцией только переменных k и t.
Отсутствие последействия состоит в том, что вероятность поступления k требований в течение промежутка времени (Т, Т+t) не зависит от того, сколько требований и как поступали до этого промежутка.
Ординарность потока требований выражает собой условие практической невозможности появления двух или нескольких требований в один и тот же момент времени. Это условие точнее сформулируем следующим образом: обозначим через P>1 (h) вероятность появления в промежутке длины h двух или более требований. Условие ординарности потока состоит в том, что при h0
Простейшим или пуассоновским потоком называется стационарный, ординарный поток без последействия.
Экспериментальная проверка показала, что простейший поток наблюдается не так часто, как это предполагалось первоначально. Действительно, предположение стационарности в реальной обстановке является довольно сильной абстракцией. Поток вызовов, поступающий на телефонную станцию, не может считаться вполне стационарным, т.к. в течение суток режим работы станции существенно меняется. Однако если рассматривать явления в сравнительно ограниченные промежутки времени, то предположение стационарности может служить достаточно удовлетворительным первым приближением.
Гипотеза отсутствия последействия во многих случаях также должна считаться недостаточно обоснованной. Имеются многочисленные явления, в которых наступление одного события влечет за собой появление других. Один телефонный звонок может повлечь за собой большое число звонков к другим абонентам.
Предположение ординарности потока во многих случаях оказывается выполненным далеко не с полной строгостью. Известно, например, что в магазины и в билетные кассы обращаются сразу группами.
система массовое обслуживание вероятность
1 Стационарное распределение вероятностей состояний СМО с двумя приборами, разнотипными требованиями и ограниченной очередью
1.1 Построение математической модели СМО с двумя приборами, разнотипными требованиями и ограниченной очередью
Рассмотрим СМО с двумя приборами, у которой входящий поток является простейшим с интенсивностью л. В СМО поступают заявки n типов. С вероятностью P (k) k=, каждая заявка независимо от предыдущих оказывается заявкой типа k. Если поступающая в СМО заявка застаёт приборы занятыми, то она становится в очередь. Всего в СМО могут находиться не более m заявок. Обслуживание производится в порядке поступления. Заявки разных типов отличаются тем, что требуют различного обслуживания. Время обслуживания в СМО распределено по показательному закону с интенсивностью мk k=, где k-тип заявки. Будем предполагать для определённости, что P (k) >0 k=; . Состоянием СМО будем называть тройку чисел (l,j,k), первое из которых показывает, сколько заявок находится в СМО в данный момент времени (l ? m), второе - номер типа заявки, обслуживаемой на первом приборе, третье - номер типа заявки, обслуживаемой на втором приборе. Если в СМО заявок нет, то для унификации обозначений будем говорить, что она находится в состоянии (0,0,0).
1.2 Построение уравнений равновесия СМО с двумя приборами, заявками двух типов и ограниченной очередью
Предположим, что в СМО существует стационарный режим. Через P (l,j,k) обозначим стационарную вероятность того, что СМО находится в состоянии (l,j,k). Процесс изменения состояния в СМО будет марковским. При поступлении на пустую станцию, сообщение любого типа будет обслуживаться первым прибором. Получим следующие уравнения для вероятностей состояния (учитываем, что P (1) +P (2) =1):
Уравнения для состояний нулевого уровня имеют вид:
2лP (0,0,0) =м1P (1,1,0) +м2P (1,2,0).
Уравнения для состояний первого уровня имеют вид:
(м1+л) P (1,1,0) =лP (1) P (0,0,0) + м1P (2,1,1) + м2P (2,1,2),
(м2+л) P (1,2,0) =лP (2) P (0,0,0) + м1P (2,2,1) + м2P (2,2,2).
Уравнения для состояний второго уровня имеют вид:
(м1+л) P (2,1,1) =лP (1) P (1,1,0) + 2м1P (1) P (3,1,1) + м2P (1) P (3,1,2) + м2P (1) P (3,2,1),
(м2+л) P (2,1,2) =лP (2) P (1,1,0) + (м1P (1) + м2P (2)) P (3,2,1) + м2P (2) P (3,1,1) + +м2P (1) P (3,2,2).
Уравнения для состояний m-того уровня имеют вид:
P (m,1,1) (2м1P (1) + м1P (2) +м2P (2)) = лP (m-1,1,1),
P (m,1,2) (м1+м2) = лP (m-1,1,2),
P (m,2,1) (м1+м2) = лP (m-1,2,1,P (m,2,2) 2м2= лP (m-1,2,2).
1.3 Построение уравнений равновесия СМО с двумя приборами, заявками n типов и ограниченной очередью
Для общего случая n типов заявок уравнения для стационарных вероятностей имеют вид:
Уравнения для состояний нулевого уровня имеют вид:
2лP (0,0,0) = . P (1,k,0).
Уравнения для состояний первого уровня имеют вид:
(мk+л) P (1,k,0) =лP (k) P (0,0,0) + P (2,k,j), k=.
Уравнения для состояний второго уровня имеют вид:
(мj+л) P (2,k,j) =лP (j) P (1,k,0) + P (j) P (3,k,r) мr+P (k) P (3,m,j) мm k,j=.
Уравнения для состояний l-го уровня имеют вид:
(л+мk+мj) P (l,k,j) =лP (l-1,k,j) + P (j) P (l+1,k,r) мr+P (k) P (l+1,m,j) мm, k,j=, (мk+мj) P (m,k,j) =лP (m-1,k,j), k,j=.
Рисунок 1-Граф переходов для двух приборов и n типов заявок.
Таким образом, мы получили стационарное распределение вероятностей состояний СМО с двумя приборами, разнотипными требованиями и ограниченной очередью.
В результате исследования были построены математические модели, графы переходов, уравнения равновесия для данных СМО.
2. Исследование стационарного распределения вероятностей состояний СМО с несколькими приборами и ограниченными очередями на приборах
2.1 Построение математической модели СМО с N приборами, разнотипными требованиями n типов, у которой на i-том приборе могут находиться не более mi требований (i=)
Рассмотрим СМО с N приборами. Поступающий поток требований является простейшим с интенсивностью л. В СМО поступают требования (заявки) n типов, причём интенсивность поступления заявок i-го типа лi, лi=л. Интенсивности обслуживания заявок i-го типа на j-том приборе мji, j=, i=. C вероятностью P (k) >0, k=1,…,n, каждая заявка независимо от предыдущих оказывается заявкой типа k. Если поступающая в СМО заявка застаёт приборы занятыми, то она становится в очередь.
Всего в СМО могут находиться на i-том приборе не более mi требований (i=). Вероятность того, что поступающая в СМО заявка становится на i-ый прибор равна Qi (i=). Обслуживание производится в порядке поступления заявок. Qi=1.
Состоянием СМО будем называть набор чисел (l,j,1…,jN), где l-число заявок в СМО в данный момент времени (0?l?m), ji, i= - номер типа заявки, обслуживаемой на i-том приборе, ji{1,2,…n}. Если в СМО заявок нет, то состояние СМО полагаем (0,0,…,0). Если заявок нет, то ji=0.
2.2 Построение уравнений равновесия для частного случая СМО с N=2 приборами, заявками n=3 типов, у которой на i-том приборе могут находиться не более mi=1 требований (i=1,2,3)
Предположим, что в СМО существует стационарный режим. Через P (l,j,k) обозначим стационарную вероятность того, что СМО находится в состоянии (l,j,k). Процесс изменения состояний в СМО является марковским, получим следующие уравнения для вероятностей состояния:
Уравнения для состояний нулевого уровня имеют вид:
2лP (0,0,0) =м11P (1,1,0) +м12P (1,2,0) +м13P (1,3,0) +м21P (1,0,1) +м22P (1,0,2) + +м23P (1,0,3).
Уравнения для состояний первого уровня имеют вид:
(лQ2+м11) P (1,1,0) = м21P (2,1,1) +м22P (2,1,2) +м23P (2,1,3) +л1Q1P (0,0,0),
(лQ2+м12) P (1,2,0) = м21P (2,2,1) +м22P (2,2,2) +м23P (2,2,3) +л2Q1P (0,0,0),
(лQ2+м13) P (1,3,0) = м21P (2,3,1) +м22P (2,3,2) +м23P (2,3,3) +л3Q1P (0,0,0),
(лQ1+м21) P (1,0,1) = м11P (2,1,1) +м12P (2,2,1) +м13P (2,3,1) +л1Q2P (0,0,0),
(лQ1+м22) P (1,0,2) = м11P (2,1,2) +м12P (2,2,2) +м13P (2,3,2) +л2Q2P (0,0,0),
(лQ1+м23) P (1,0,3) = м11P (2,1,3) +м12P (2,2,3) +м13P (2,3,3) +л3Q2P (0,0,0).
Уравнения для состояний второго уровня имеют вид:
(м11+ м21) P (2,1,1) =л1Q1P (1,0,1) +л1Q2P (1,1,0),
(м11+ м22) P (2,1,2) =л1Q1P (1,0,2) +л2Q2P (1,1,0),
(м11+ м23) P (2,1,3) =л1Q1P (1,0,3) +л3Q2P (1,1,0),
(м21+ м12) P (2,2,1) =л2Q1P (1,0,1) +л1Q2P (1,2,0),
(м22+ м12) P (2,2,2) =л2Q1P (1,0,2) +л2Q2P (1,2,0),
(м23+ м12) P (2,2,3) =л2Q1P (1,0,3) +л3Q2P (1,2,0),
(м21+ м13) P (2,3,1) =л3Q1P (1,0,1) +л1Q2P (1,3,0),
(м22+ м13) P (2,3,2) =л3Q1P (1,0,2) +л2Q2P (1,3,0),
(м23+ м13) P (2,3,3) =л3Q1P (1,0,3) +л3Q2P (1,3,0).
Уравнения для состояния нулевого уровня в общем виде:
2лP (0,0,0) =м1jP (1,j,0) +м2kP (1,0,k).
Уравнения для состояния первого уровня в общем виде:
(лQ2+м1j) P (1,j,0) =м2kP (2,j,k) +лjQ1P (0,0,0), j={1,2,3},
(лQ1+м2j) P (1,0,j,) =м1kP (2,k,j) +лjQ2P (0,0,0), j={1,2,3}.
Уравнения для состояния второго уровня в общем виде:
(м2j+ м1k) P (2,j,k) =лkQ1P (1,0,k) + лjQ2P (1,j,0), k,j={1,2,3}.
Рисунок 2 - Граф переходов для СМО с N=2 приборами, заявками n=3 типов
2.3 Построение математической модели СМО с N приборами, разнотипными требованиями n типов, которых может находиться в СМО не более m
Рассмотрим СМО с N приборами. Поступающий поток требований является простейшим с интенсивностью л. В СМО поступают требования (заявки) n типов, причём интенсивность поступления заявок i-го типа лi, лi=л. Интенсивности обслуживания заявок i-го типа на j-том приборе мji, j=, i=. C вероятностью P (k) >0, k=, каждая заявка независимо от предыдущих оказывается заявкой типа k. Если поступающая в СМО заявка застаёт приборы занятыми, то она становится в очередь.
Всего в СМО могут находиться не более m требований. Вероятность того, что поступающая в СМО заявка становится на i-ый прибор равна Qi (i=). Обслуживание производится в порядке поступления заявок. Qi=1.
Состоянием СМО будем называть набор чисел (l,j,1…,jN), где l-число заявок в СМО в данный момент времени (0?l?m), ji, i= - номер типа заявки, обслуживаемой на i-том приборе, ji{1,2,…n}. Если в СМО заявок нет, то состояние СМО полагаем (0,0,…,0). Если заявок нет, то ji=0.
2.4 Построение уравнений равновесия для частного случая СМО с N=2 приборами, заявками n=3 типов, которых может находиться в СМО не более m=2
Предположим, что в СМО существует стационарный режим. Через P (l,j,k) обозначим стационарную вероятность того, что СМО находится в состоянии (l,j,k). Процесс изменения состояний в СМО является марковским, получим следующие уравнения для вероятностей состояния:
Уравнения для состояний нулевого уровня имеют вид:
2лP (0,0,0) =м11P (1,1,0) +м12P (1,2,0) +м13P (1,3,0) +м21P (1,0,1) +м22P (1,0,2) + +м23P (1,0,3).
Уравнения для состояний первого уровня имеют вид:
(л+м11) P (1,1,0) =м11P (1) P (2,1,0) +м12P (1) P (2,2,0) +м13P (1) P (2,3,0) +м21P (2,1,1) + +м22P (2,1,2) +м23P (2,1,3) +л1Q1P (0,0,0),
(л+м12) P (1,2,0) =м11P (2) P (2,1,0) +м12P (2) P (2,2,0) +м13P (2) P (2,3,0) +м21P (2,2,1) + +м22P (2,2,2) +м23P (2,2,3) +л2Q1P (0,0,0),
(л+м13) P (1,3,0) =м11P (3) P (2,1,0) +м12P (3) P (2,2,0) +м13P (3) P (2,3,0) +м21P (2,3,1) + +м22P (2,3,2) + м23P (2,3,3) +л3Q1P (0,0,0),
(л+м21) P (1,0,1) =м21P (1) P (2,0,1) +м22P (1) P (2,0,2) +м23P (1) P (2,0,3) +м11P (2,1,1) м12P (2,2,1) +м13P (2,3,1) +л1Q2 P (0,0,0),
(л+м22) P (1,0,2) =м21P (2) P (2,0,1) +м22P (2) P (2,0,2) +м23P (2) P (2,0,3) +м11P (2,1,2) + +м12P (2,2,2) +м13P (2,3,2) +л2Q2 P (0,0,0),
(л+м23) P (1,0,3) =м21P (3) P (2,0,1) +м22P (3) P (2,0,2) +м23P (3) P (2,0,3) +м11P (2,1,3) + +м12P (2,2,3) + м13P (2,3,3) +л3Q2 P (0,0,0).
Уравнения для состояний второго уровня имеют вид:
м11P (2,1,0) =лQ1P (1,1,0),
м12P (2,2,0) =лQ1P (1,2,0),
м13P (2,3,0) =лQ1P (1,3,0),
(м11+ м21) P (2,1,1) =л1Q1P (1,0,1) +л1Q2P (1,1,0),
(м11+ м22) P (2,1,2) =л1Q1P (1,0,2) +л2Q2P (1,1,0),
(м11+ м23) P (2,1,3) =л1Q1P (1,0,3) +л3Q2P (1,1,0),
(м21+ м12) P (2,2,1) =л2Q1P (1,0,1) +л1Q2P (1,2,0),
(м22+ м12) P (2,2,2) =л2Q1P (1,0,2) +л2Q2P (1,2,0),
(м23+ м12) P (2,2,3) =л2Q1P (1,0,3) +л3Q2P (1,2,0),
(м21+ м13) P (2,3,1) =л3Q1P (1,0,1) +л1Q2P (1,3,0),
(м22+ м13) P (2,3,2) =л3Q1P (1,0,2) +л2Q2P (1,3,0),
(м23+ м13) P (2,3,3) =л3Q1P (1,0,3) +л3Q2P (1,3,0),
м21P (2,0,1) =лQ2P (1,0,1),
м22P (2,0,2) =лQ2P (1,0,2),
м23P (2,0,3) =лQ2P (1,0,3).
Уравнения для состояния нулевого уровня в общем виде:
2лP (0,0,0) =м1jP (1,j,0) +м2kP (1,0,k).
Уравнения для состояния первого уровня в общем виде:
(л+м1j) P (1,j,0) =P (j) м1iP (2, i,0) +м2kP (2,j,k) +лjQ1P (0,0,0), j={1,2,3},
(л+м2j) P (1,0,j) =P (j) м2iP (2,0, i) +м1kP (2,k,j) +лjQ2P (0,0,0), j={1,2,3},
Уравнения для состояния второго уровня в общем виде:
м1jP (2,j,0) =лQ1P (1,j,0), j={1,2,3},
(м2j+ м1k) P (2,j,k) =лkQ1P (1,0,k) + лjQ2P (1,j,0), k,j={1,2,3},
м2jP (2,0,j) =лQ2P (1,0,j), j={1,2,3}.
Таким образом, нами было исследовано стационарное распределение вероятностей состояний СМО с несколькими приборами и ограниченными очередями на приборах.
В частности, были рассмотрены следующие СМО:
СМО с N приборами, разнотипными требованиями n типов, у которой на i-том приборе может находиться не более mi требований (i=).
СМО с N=2 приборами, заявками n=3 типов, которых может находиться в СМО не более m=2.
В результате исследования были построены математические модели, графы переходов, уравнения равновесия для данных СМО.
3. Стационарное распределение СМО с N приборами, заявками n типов, у которой на каждом приборе имеется буфер ёмкости M, причём заявки i-того типа обслуживаются на i-том приборе
3.1 Построение математической модели СМО с N=2 приборами, заявками N=2 типов, у которой на каждом приборе имеется буфер ёмкости m=2, причём заявки i-того типа обслуживаются на i-том приборе
Рассмотрим СМО с N=2 приборами. Поступающий поток требований является простейшим с интенсивностью л. В СМО поступают заявки N=2 типов, причём интенсивность поступления заявок i-го типа лi, i=, . Заявки i-того типа обслуживаются на i-том приборе. Интенсивность обслуживания заявок i-того типа - мi, i=. На каждом приборе имеется буфер ёмкости m=2, в который становятся заявки соответствующих типов. С вероятностью p (k) >0, k= каждая заявка независимо от предыдущих оказывается заявкой типа k. Если поступающая в СМО заявка застаёт приборы занятыми, то она становится в очередь. Обслуживание производится в порядке поступления заявок. Состоянием СМО будем называть тройку чисел (l,j,k), где l - число заявок в СМО в данный момент времени (0?l ? 4), j - число заявок 1-го типа, обслуживаемых на первом приборе, k - число заявок 2-го типа, обслуживаемых на втором приборе. Если в СМО заявок нет, то для унификации обозначений будем говорить, что она находится в состоянии (0,0,0).
Рисунок 3 - Граф переходов для СМО с двумя приборами, заявками двух типов, в которой заявки i-го типа обслуживаются на i-ом приборе
3.2 Построение уравнений равновесия для частного случая СМО с N=2 приборами, заявками N=2 типов, у которой на каждом приборе имеется буфер ёмкости m=2, причём заявки i-того типа обслуживаются на i-том приборе
Предположим, что в СМО существует стационарный режим.
Через P (l,j,k) обозначим стационарную вероятность того, что СМО находится в состоянии (l,j,k).
Процесс изменения состояния в СМО будет марковским. Получим следующие уравнения для вероятностей состояния:
Уравнения для состояний нулевого уровня имеют вид:
2лP (0,0,0) =м1P (1,1,0) +м2P (1,0,1).
Уравнения для состояний первого уровня имеют вид:
(л +м1) P (1,1,0) =л1P (0,0,0) + м1P (2,2,0) + м2P (2,1,1),
(л +м2) P (1,0,1) =л2P (0,0,0) + м1P (2,1,1) + м2P (2,0,2).
Уравнения для состояний второго уровня имеют вид:
(л+ м1) P (2,2,0) =л1P (1,1,0) + м1P (3,3,0) + м2P (3,2,1),
(л +м1+ м2) P (2,1,1) =л1P (1,0,1) + л2P (1,1,0) + м1P (3,2,1) +м2P (3,1,2),
(л+ м2) P (2,0,2) =л2P (1,0,1) + м1P (3,1,2) + м2P (3,0,3).
Уравнения для состояний третьего уровня имеют вид:
(л+ м1) P (3,3,0) =л1P (2,2,0) + м1P (4,4,0) + м2P (4,3,1),
(л +м1+ м2) P (3,1,2) =л1P (2,0,2) + л2P (2,1,1) + м1P (4,2,2) +м2P (4,1,3),
(л +м1+ м2) P (3,2,1) =л1P (2,1,1) + л2P (2,2,0) + м1P (4,3,1) +м2P (4,2,2),
(л+ м2) P (3,0,3) =л2P (2,0,2) + м1P (4,1,3) + м2P (4,0,4).
Уравнения для состояний четвёртого уровня имеют вид:
м1P (4,4,0) =л1P (3,3,0),
(м1+ м2) P (4,1,3) =л1P (3,0,3) + л2P (3,1,2),
(м1+ м2) P (4,2,2) =л1P (3,1,2) + л2P (3,2,1),
(м1+ м2) P (4,3,1) =л1P (3,2,1) + л2P (3,3,0),
м2P (4,0,4) =л1P (3,3,0).
3.3 Построение математической модели СМО с N=3 приборами, заявками N=3 типов, у которой на каждом приборе имеется буфер ёмкости m=2, причём заявки i-того типа обслуживаются на i-том приборе
Рассмотрим СМО с N=3 приборами. Поступающий поток требований является простейшим с интенсивностью л. В СМО поступают заявки N=3 типов, причём интенсивность поступления заявок i-го типа лi, i=, . Заявки i-того типа обслуживаются на i-том приборе. Интенсивность обслуживания заявок i-того типа - мi, i=, . На каждом приборе имеется буфер ёмкости m=2, в который становятся заявки соответствующих типов. С вероятностью p (k) >0, k= каждая заявка независимо от предыдущих оказывается заявкой типа k. Если поступающая в СМО заявка застаёт приборы занятыми, то она становится в очередь. Обслуживание производится в порядке поступления заявок. Состоянием СМО будем называть четвёрку чисел (l,j,k,s), где l - число заявок в СМО в данный момент времени (0?l ? 6), j - число заявок 1-го типа, обслуживаемых на первом приборе, k - число заявок 2-го типа, обслуживаемых на втором приборе, k - число заявок 3-го типа, обслуживаемых на третьем приборе. Если в СМО заявок нет, то для унификации обозначений будем говорить, что она находится в состоянии (0,0,0,0).
3.4 Построение уравнений равновесия для частного случая СМО с N=3 приборами, заявками N=3 типов, у которой на каждом приборе имеется буфер ёмкости m=2, причём заявки i-того типа обслуживаются на i-том приборе
Предположим, что в СМО существует стационарный режим. Через P (l,j,k,s) обозначим стационарную вероятность того, что СМО находится в состоянии (l,j,k,s). Процесс изменения состояния в СМО будет марковским. Получим следующие уравнения для вероятностей состояния:
Уравнения для состояний нулевого уровня имеют вид:
лP (0,0,0,0) =м1P (1,1,0,0) +м2P (1,0,1,0) + м3P (1,0,0,1).
Уравнения для состояний первого уровня имеют вид:
(л +м1) P (1,1,0,0) =л1P (0,0,0,0) + м1P (2,2,0,0) + м2P (2,1,1,0) + м3P (2,1,0,1),
(л +м2) P (1,0,1,0) =л2P (0,0,0,0) + м1P (2,1,1,0) + м2P (2,0,2,0) + м3P (2,0,1,1),
(л +м3) P (1,0,0,1) =л3P (0,0,0,0) + м1P (2,1,0,1) + м2P (2,0,1,1) + м3P (2,0,0,2),
Уравнения для состояний второго уровня имеют вид:
(л+ м1) P (2,2,0,0) =л1P (1,1,0,0) + м1P (3,3,0,0) + м2P (3,2,1,0) + м3P (3,2,0,1),
(л+ м2) P (2,0,2,0) =л2P (1,0,1,0) + м1P (3,1,2,0) + м2P (3,0,3,0) + м3P (3,1,1,1),
(л+ м3) P (2,0,0,2) =л3P (1,0,0,1) + м1P (3,0,1,2) + м2P (3,0,2,1) + м3P (3,0,0,3),
(л +м1+м2) P (2,1,1,0) =л1P (1,0,1,0) + л2P (1,1,0,0) + м1P (3,2,1,0) +м2P (3,1,2,0) + +м2P (3,1,1,1),
(л +м2+м3) P (2,0,1,1) =л2P (1,0,0,1) + л3P (1,0,1,0) + м1P (3,1,1,1) +м2P (3,0,2,1) + +м3P (3,0,1,2),
(л +м1+ м3) P (2,1,0,1) =л1P (1,0,0,1) + л3P (1,1,0,0) + м1P (3,2,0,1) +м3P (3,1,0,2) + +м2P (3,1,1,1).
Уравнения для состояний третьего уровня имеют вид:
(л+ м1) P (3,3,0,0) =л1P (2,2,0,0) + м1P (4,4,0,0) + м2P (4,3,1,0) + м3P (4,3,0,1),
(л+ м2) P (3,0,3,0) =л2P (2,0,2,0) + м1P (4,1,3,0) + м2P (4,0,4,0) + м3P (4,0,3,1),
(л+ м3) P (3,0,0,3) =л3P (2,0,0,2) + м1P (4,1,0,3) + м2P (4,0,1,3) + м3P (4,0,0,4),
(л +м1+ м2) P (3,1,2,0) =л1P (2,0,2,0) + л2P (2,1,1,0) + м1P (4,2,2,0) +м2P (4,1,3,0) + +м3P (4,1,2,1),
(л +м1+ м2) P (3,1,0,2) =л1P (2,0,0,2) + л3P (2,1,0,1) + м1P (4,2,0,2) +м2P (4,1,1,2) + +м3P (4,1,0,3),
(л +м1+ м3) P (3,2,1,0) =л1P (2,1,1,0) + л2P (2,2,0,0) + м1P (4,3,1,0) +м2P (4,2,2,0) + +м3P (4,2,1,1),
(л +м1+ м3) P (3,2,0,1) =л1P (2,1,0,1) + л3P (2,0,0,2) + м1P (4,3,0,1) +м2P (4,2,1,1) + +м3P (4,2,0,2),
(л +м2+ м3) P (3,0,1,2) =л2P (2,0,0,2) + л3P (2,0,1,1) + м1P (4,1,1,2) +м2P (4,0,2,2) + +м3P (4,0,1,3),
(л +м2+ м3) P (3,0,2,1) =л2P (2,0,1,1) + л3P (2,0,2,0) + м1P (4,1,2,1) +м2P (4,0,3,1) + +м3P (4,0,2,2),
(л +м) P (3,1,1,1) =л1P (2,0,1,1) + л2P (2,1,0,1) + л3P (2,1,1,0) + +м1P (4,2,1,1) +м2P (4,1,2,1) + м3P (4,1,1,2).
Уравнения для состояний шестого уровня имеют вид:
м1P (6,6,0,0) =л1P (5,5,0,0),
м2P (6,0,6,0) =л1P (5,0,5,0),
м3P (6,0,0,6) =л1P (5,0,0,5),
(м1+ м2) P (6,5,1,0) =л1P (5,4,1,0) + л2P (5,5,0,0),
(м1+ м3) P (6,5,1,0) =л1P (5,4,0,1) + л3P (5,5,0,0),
(м2+ м3) P (6,0,5,1) =л2P (5,0,4,1) + л3P (5,0,5,0),
(м1+ м2) P (6,1,5,0) =л1P (5,0,5,0) + л2P (5,1,4,0),
(м1+ м3) P (6,1,0,5) =л1P (5,0,0,5) + л3P (5,1,0,4),
(м2+ м3) P (6,0,1,5) =л2P (5,0,0,5) + л3P (5,0,1,4),
(м1+ м2) P (6,3,3,0) =л1P (5,2,3,0) + л2P (5,3,2,0),
(м1+ м3) P (6,3,0,3) =л1P (5,2,0,3) + л3P (5,3,0,2),
(м2+ м3) P (6,0,3,3) =л2P (5,0,2,3) + л3P (5,0,3,2),
(м1+ м2) P (6,2,4,0) =л1P (5,1,4,0) + л2P (5,2,3,0),
(м1+ м3) P (6,2,0,4) =л1P (5,1,0,4) + л3P (5,2,0,3),
(м2+ м3) P (6,0,2,4) =л2P (5,0,1,4) + л3P (5,0,2,3),
(м1+ м2) P (6,4,2,0) =л1P (5,3,2,0) + л2P (5,4,1,0),
(м1+ м3) P (6,4,0,2) =л1P (5,3,0,2) + л3P (5,4,0,1),
(м2+ м3) P (6,0,5,1) =л2P (5,0,3,2) + л3P (5,0,4,1),
мP (6,2,2,2) =л1P (5,1,2,2) + л2P (5,2,1,2) + л3P (5,2,2,1),
мP (6,1,2,3) =л1P (5,0,2,3) + л2P (5,1,1,3) + л3P (5,1,2,2),
мP (6,3,1,2) =л1P (5,2,1,2) + л2P (5,3,0,2) + л3P (5,3,1,1),
мP (6,1,3,2) =л1P (5,0,3,2) + л2P (5,1,2,2) + л3P (5,1,3,1),
мP (6,2,1,3) =л1P (5,1,1,3) + л2P (5,2,0,3) + л3P (5,2,1,2),
мP (6,3,2,1) =л1P (5,2,2,1) + л2P (5,3,1,1) + л3P (5,3,2,0),
мP (6,2,3,1) =л1P (5,1,3,1) + л2P (5,2,2,1) + л3P (5,2,3,1),
мP (6,4,1,1) =л1P (5,3,1,1) + л2P (5,4,0,1) + л3P (5,4,1,0),
мP (6,1,4,1) =л1P (5,0,4,1) + л2P (5,1,3,1) + л3P (5,1,4,0),
мP (6,1,1,4) =л1P (5,0,1,4) + л2P (5,1,0,4) + л3P (5,1,1,3).
Составим уравнения равновесия в общем виде (будем учитывать, что :
l) Уравнения для состояний нулевого уровня в общем виде:
лP (0,0,0,0) =мiP (1,) I{}.
Уравнения для состояний первого уровня имеют вид:
(л+мi) I{} P (1,) =лiP (0,0,0,0) +мj P (2,) I{}, i=1,2,3.
Уравнения для состояний второго уровня имеют вид:
(л+мi) I{} P (2,) = лiP (1, ) I{} +
+мj P (3,) I{}, i=1,2,3, (л+мi+мj) I{,} P (2,) =
= лiP (1, ) I{} + лjP (1, ) I{} +
+мk P (3,) I{}I{}, i=1,2,3, j=1,2,3.
Уравнения для состояний третьего уровня имеют вид:
(л+мi) I{,} P (3,) = лiP (2, ) I{} +
+мj P (4,) I{}, i=1,2,3, (л+мi+мj) I{} P (3,) =
=лiP (2, ) I{} + лjP (2, ) I{}+ +мk P (4,) I{}I{}, i=1,…,3, j=1,2,3, (л+м) I{} P (3,) =
=лnP (2, ) I{} + +мmP (4,) I{}, i=1,2,3, j=1,2,3, k=1,2,3.
Уравнения для состояний шестого уровня имеют вид:
мiP (6,) I{} = лiP (5, ) I{}, i=1,2,3, (мi+ мj) P (6,) I{}=
=лiP (5, ) I{} +
+лjP (5, ) I{}, i=1,…,3, k=1,…,3, j=1,2,3,мP (6,) = лiP (5,) I{}, i,j=.
3.5 Построение математической модели СМО с N приборами, заявками N типов, у которой на каждом приборе имеется буфер ёмкости m, причём заявки i-того типа обслуживаются на i-том приборе
Рассмотрим СМО с N приборами. Поступающий поток требований является простейшим с интенсивностью л. В СМО поступают заявки N типов, причём интенсивность поступления заявок i-го типа лi, i=, . Заявки i-того типа обслуживаются на i-том приборе. Интенсивность обслуживания заявок i-того типа - мi, i=. На каждом приборе имеется буфер ёмкости m, в который становятся заявки соответствующих типов. С вероятностью p (k) >0, k= каждая заявка независимо от предыдущих оказывается заявкой типа k. Если поступающая в СМО заявка застаёт приборы занятыми, то она становится в очередь. Обслуживание производится в порядке поступления заявок. Состоянием СМО будем называть набор чисел (l,j1,…,jN), где l - число заявок в СМО в данный момент времени (0?l ? mN), ji - число заявок i-го типа (i=), обслуживаемых на i-том приборе. Если в СМО заявок нет, то для унификации обозначений будем говорить, что она находится в состоянии (0,0,…,0). Если заявок нет, то ji =0.
3.6 Построение уравнений равновесия для СМО с N приборами, заявками N типов, у которой на каждом приборе имеется буфер ёмкости m, причём заявки i-того типа обслуживаются на i-том приборе
Предположим, что в СМО существует стационарный режим. Через P (l, ) обозначим стационарную вероятность того, что СМО находится в состоянии (l, ). Процесс изменения состояния в СМО будет марковским. Получим следующие уравнения для вероятностей состояния (будем учитывать, что l):
Уравнения для состояний нулевого уровня:
лP (0,0,0,0) =мiP (1,) I{}.
Уравнения для состояний первого уровня имеют вид:
(л+мi) I{} P (1, ) =лiP (0,0,0,0) +
+мj P (2, ) I{}, i=.
Уравнения для состояний второго уровня имеют вид:
(л+мi) I{} P (2, ) = лiP (1, ) I{} +
+мj P (3, ) I{}, i=,
(л+мi+мj) I{,} P (2, ) = лiP (1, ) I{} + +лjP (1, ) I{} +
+мk P (3, ) I{}I{}, i=, j=.
Уравнения для состояний третьего уровня имеют вид:
(л+мi) I{,} P (3,) = лiP (2, ) I{} +
+мj P (4, ) I{}, i=,
(л+мi+мj) I{} P (3, ) =
=лiP (2, ) I{} +лjP (2, ) I{}+ +мk P (4, ) I{}I{},
i=, j=,
(л+м) I{} P (3, ) =
=лnP (2, ) I{} +
+мm P (4, ) I{}, i=, j=, k=.
Уравнения для состояний N-го уровня имеют вид:
мiI{} P (Nm,) = лiP (Nm-1, ) I{}, i=,
(мi+мj) I{} P (Nm,) =
=лiP (Nm-1, ) I{}+ лjP (Nm-1, ) I{}, i=, j=.
…
I{} P (Nm,) =
=P (Nm-1, ) I{}, i=, j=,
мP (Nm,) =P (Nm-1, ) I{}, i=,
j=.
Таким образом, нами было исследовано стационарное распределение вероятностей состояний СМО с N приборами, заявками N типов, у которой на каждом приборе имеется буфер ёмкости m, причём заявки i-того типа обслуживаются на i-том приборе.
В частности, были рассмотрены следующие СМО:
СМО с N=2 приборами, заявками N=2 типов, у которой на каждом приборе имеется буфер ёмкости m=2, причём заявки i-того типа обслуживаются на i-том приборе.
СМО с N=3 приборами, заявками N=3 типов, у которой на каждом приборе имеется буфер ёмкости m=2, причём заявки i-того типа обслуживаются на i-том приборе.
В результате исследования были построены математические модели, графы переходов, уравнения равновесия для данных СМО.
4. Построение векторно-матричных уравнений СМО с двумя приборами, требованиями трех типов и конечными буферами
Состояние СМО, вообще говоря, задается тройкой чисел (l,j,k), где l - число заявок в СМО в данный момент времени (l0), j, - номер типа заявки, обслуживаемой на первом приборе, k - номер типа заявки, обслуживаемой на втором приборе. Однако, если перенумеровать все пары (j,k) через m (в случае, когда j{1,2,3}, k{1,2,3}), то состояние СМО можно задавать парой (l,m), где l>1, m{1,.,9}. В случае, когда l=0 имеем пару (0,0), т.е. m=0 (j=0, k=0).
Будем называть множество { (l,m) mL} l-тым уровнем. Обозначим через P (l) вектор стационарных вероятностей l-ого уровня:
P (l) = (P (l,1), P (l,2),., P (l,9)), l>1;
P (1) = (P (1,1), P (1,2),., P (1,6));
P (0) = P (0,0).
Очевидно, что возможные переходы для l-ого уровня схематически можно изобразить следующим образом:
Рисунок 4
Так как на нулевом уровне имеем лишь одно состояние (0,0), то при l=0 L={0}. На первом уровне получаем множество состояний (1,m), где m {1,.,6}, т.е. L={1,2,3,4,5,6}, поскольку в парах (j,k) при l=1 либо j=0 и k{1,2,3}, либо j{1,2,3} и k=0.
На l-том уровне (l>1) имеем множество состояний { (l,k) k=1,.,9}, следовательно, L={1,2,.,9}.
Кроме того, для стационарных вероятностей выполняются следующие соотношения:
Q=0; 1=1,где= (Р (0), Р (1), Р (2), Р (3),.),
Q - инфинитезимальная матрица, имеющая следующий блочный вид:
Для l 2 введем следующие обозначения:
Аlj=A1 при l > j; All=A2; Alj=A3 при l < j.
Таким образом, из соотношения Q=0 имеем с учетом обозначений:
Р (0) A00+ Р (1) A10=0,Р (0) A01+ Р (1) A11+ Р (2) A21=0,Р (1) A12+ Р (2) A2+ Р (3) A1=0, (4.1)
Р (2) A3+ Р (3) A2+ Р (4) A1=0,.
Р (l-1) A3+ Р (l) A2+ Р (l+1) A1=0,.
P (m-2) A3+ P (m-1) A2+ P (m) A1=0
P (m-1) A3+P (m) A2=0
Состояние (l,j,k) при l>2 переходит в состояние (l+1,j,k) с интенсивностью , в состояние (l-1,j',k) - с 1jp (j'), в состояние (l-1,j,k') - с 2kp (k'), где 1j - интенсивность обслуживания заявки j-ого типа на первом приборе; 2k - интенсивность обслуживания заявки k-ого типа на втором приборе; p (j') - вероятность появления заявки типа j'; p (k') - вероятность появления заявки типа k'; j, k, j', k'{1,2,3}.
Введем следующие обозначения:
1
2
= 3
вектор-столбец интенсивностей обслуживания заявок различных типов;
где 1= 11+21+31,2= 12+22+32,3= 13+23+33,
рТ= (р (1), р (2), р (3)),
I - единичная (nn) матрица.
С учетом введенных обозначений получаем:
A00 = - 2.
Для блоков A01 имеем, что состояние (0,0,0) переходит в сотояние (1,j,0) с интенсивностью p (j) и в состояние (1,0,k) с интенсивностью p (k). Следовательно, блок A01 - строка, состоящая из шести элементов и имеющая вид:
A01= (p (1) p (2) p (3) p (1) p (2) p (3)) или A01= (pТ pТ).
Для блоков A10 получаем, что состояние (1,j,0) переходит в состояние (0,0,0) с интенсивностью 1j, а состояние (1,0,k) переходит в состояние (0,0,0) с интенсивностью 2k. Таким образом, A10 - столбец, состоящий из шести элементов, который имеет вид:
Далее, для блоков A12 имеем, что состояние (1,j,0) переходит в состояние (2,j,k) с интенсивностью р (k), а состояние (1,0,k) переходит в состояние (2,j,k) с интенсивностью p (j). Следовательно, блоки A12 размерности (69) и вида:
или с учетом обозначений
Учитывая вид построенных блоков A10 и A12, блоки A11 размерности (66) и также имеют диагональный вид:
Для блоков A21 имеем, что состояние (2,j,k) переходит в состояние (1,j,0) с интенсивностью 2k, а в состояние (1,0,k) с интенсивностью 1j. Следовательно, блоки A21 размерности (96) и вида:
или с учетом обозначений:
С учетом введенных ранее обозначений блоки A1 имеют вид:
Где
Блоки A3 также имеют размерность (99) и диагональный вид:
т.е.
Блоки A2 также диагонального вида и размерности (99):
Для стационарных вероятностей можно записать:
P (l) = P (l-1) R, l>2.
Таким образом, последнее соотношение из (4.1) будет иметь вид:
P (l-1) (A3+R A2+R2 A1) =0 (4.2)
Введем следующие обозначения:
k=k/,
D= A2/.
Учитывая вид блоков A1, получим:
A1=рТI+IрТ, A1/=рТI+ IрТ,
где - символ кронекеровского произведения.
Кроме того, последнее соотношение из (4.1) примет вид:
P (l-1) +P (l) D+P (l+1) (рТI+ IрТ) = 0
С учетом соотношения P (l) =P (l-1) R, l>2 (4.3) получим:
P (l-1) +P (l) D+P (l) R (рТI+ IрТ) = 0 или
P (l-1) +P (l) (D+R (рТI+ IрТ)) = 0.
Следовательно, принимая во внимание соотношение (4.3), имеем:
R = - (D+R (рТI+ IрТ)) - 1.
Полагая R (0) =0, R (N+1) = - (D+R (N) (рТI+ IрТ)) - 1, получим, что R (N+1) R при N.
Таким образом, получены векторно-матричные уравнения для определения стационарных вероятностей состояний СМО с двумя приборами, тремя типами заявок и конечными очередями на каждом приборе.
Заключение
В ходе написания дипломной работы было исследовано стационарное распределение вероятностей в СМО с различным числом приборов, требованиями различных типов и ограниченными очередями на приборах.
В частности, были рассмотрены следующие СМО:
1) СМО с двумя приборами, разнотипными требованиями и ограниченной очередью.
2) СМО с несколькими приборами и ограниченными очередями на приборах:
СМО с N приборами, разнотипными требованиями n типов, у которой на i-том приборе может находиться не более mi требований (i=).
СМО с N=2 приборами, заявками n=3 типов, которых может находиться в СМО не более m=2.
3) СМО с N приборами, заявками N типов, у которой на каждом приборе имеется буфер ёмкости m, причём заявки i-того типа обслуживаются на i-том приборе:
СМО с N=2 приборами, заявками N=2 типов, у которой на каждом приборе имеется буфер ёмкости m=2, причём заявки i-того типа обслуживаются на i-том приборе.
СМО с N=3 приборами, заявками N=3 типов, у которой на каждом приборе имеется буфер ёмкости m=2, причём заявки i-того типа обслуживаются на i-том приборе.
В результате исследования были построены математические модели, графы переходов, уравнения равновесия для данных СМО.
Список использованных источников
1. Климов, Г.П. Стохастические системы обслуживания [Текст]: главная редакция физико-математической литературы / Г.П. Климов. - М.: Наука, 1966. - 265 с.
2. Медведев, Г.А. Приближенный анализ процессов передачи разнотипных сообщений по сети связи [Текст]: Управляемые системы массового обслуживания/ Г.А. Медведев. - Томск, 1986. - 187с.
3. Назаров, А.А. Асимптотический анализ марковизируемых систем [Текст] / А. А. Назаров. - Томск: Изд-во Том. Ун-та, 1991. - 158с.
4. Гнеденко, Б.В. Введение в теорию массового обслуживания [Текст]: учеб. пособие для вузов / Б.В. Гнеденко, И.Н. Коваленко. - М.: Наука. Главная редакция физико-математической литературы, 1987. - 574 с.
5. Мишина, А.П. Высшая алгебра [Текст]: учеб. пособие для вузов /
А.П. Мишина, И.В. Проскуряков. - М.: Наука. Главная редакция
физико-математической Литературы, 1965. - 283 с.
6. Малинковский, Ю.В. Методические указания по курсу "Теория массового обслуживания" для студентов 4 курса специальности "Математика" Часть 1 [Текст]: учеб. пособие для вузов / Ю.В. Малинковский, Е.А. Ковалев. - Гомель, ГГУ, 1985. - 25 с.
Размещено на Allbest.ru
Подобные документы
Основные понятия теории массового обслуживания: марковский процесс, простой поток, сеть Джексона. Исследование стационарного распределения сети с ромбовидным контуром: для марковских и немарковских процессов, а также для сети с отрицательными заявками.
дипломная работа [957,4 K], добавлен 17.12.2012Оптимизация управления потоком заявок в сетях массового обслуживания. Методы установления зависимостей между характером требований, числом каналов обслуживания, их производительностью и эффективностью. Теория графов; уравнение Колмогoрова, потоки событий.
контрольная работа [35,0 K], добавлен 01.07.2015Исследование стационарного распределения сетей массового обслуживания и доказательство инвариантности. Уравнения глобального равновесия и понятие эргодичности. Доказательство инвариантности стационарного распределения, а также определение его вида.
дипломная работа [439,7 K], добавлен 12.12.2009Понятие системы массового обслуживания, ее сущность и особенности. Теория массового обслуживания как один из разделов теории вероятностей, рассматриваемые вопросы. Понятие и характеристика случайного процесса, его виды и модели. Обслуживание с ожиданием.
курсовая работа [1,4 M], добавлен 15.02.2009Общая структура системы массового обслуживания. Каналы и линии связи, вычислительные машины, объединенные общей структурой, число каналов обслуживания. Регулярный поток с ограниченным последействием. Применение различных величин и функций в системе.
курсовая работа [199,4 K], добавлен 13.11.2011Математическая теория массового обслуживания как раздел теории случайных процессов. Системы массового обслуживания заявок, поступающих через промежутки времени. Открытая марковская сеть, ее немарковский случай, нахождение стационарных вероятностей.
курсовая работа [374,3 K], добавлен 07.09.2009Теория массового обслуживания – область прикладной математики, анализирующая процессы в системах производства, в которых однородные события повторяются многократно. Определение параметров системы массового обслуживания при неизменных характеристиках.
курсовая работа [439,6 K], добавлен 08.01.2009Примеры процессов размножения и гибели в случае простейших систем массового обслуживания. Математическое ожидание для системы массового обслуживания. Дополнительный поток и бесконечное число приборов. Система с ограничением на время пребывания заявки.
курсовая работа [1003,1 K], добавлен 26.01.2014Составление имитационной модели и расчет показателей эффективности системы массового обслуживания по заданны параметрам. Сравнение показателей эффективности с полученными путем численного решения уравнений Колмогорова для вероятностей состояний системы.
курсовая работа [745,4 K], добавлен 17.12.2009Определение случайного процесса и его характеристики. Основные понятия теории массового обслуживания. Понятие марковского случайного процесса. Потоки событий. Уравнения Колмогорова. Предельные вероятности состояний. Процессы гибели и размножения.
реферат [402,0 K], добавлен 08.01.2013
