CC BY
48
УДК 519.248
МОДЕЛИРОВАНИЕ КОНЕЧНЫМИ АВТОМАТАМИ СИСТЕМ МАССОВОГО ОБСЛУЖИВАНИЯ С РАЗЛИЧИМЫМИ КАНАЛАМИ
© 2014 А.П. Котенко, М.Б. Букаренко
Самарский государственный аэрокосмический университет имени академика С.П.Королёва (национальный исследовательский университет)
Поступила в редакцию 17.12.2014
Изложена постановка задачи исследования систем массового обслуживания с различимыми каналами, обладающими различной пропускной способностью и/или отдельными очередями. В этом случае нельзя применить результаты классической теории массового обслуживания с линейным графом состояний процесса "гибели и размножения". Предложено описание системы с различимыми каналами конечными детерминированными и недетерминированными автоматами. Это позволило провести имитационное моделирование поведения системы при различных вариантах диспетчеризации входных заявок. Ключевые слова: система массового обслуживания, неоднородные каналы обслуживания, раздельные очереди каналов обслуживания, конечные автоматы.
В нотификации [1, 2, 3] рассмотрим систему массового обслуживания (СМО) с различимыми каналами (имеющими, к примеру, разную пропускную способность или раздельные очереди) сигнатуры Т=Т(м1,м2,...,м1!;т1,т2,...,т1!), где мг-пропускная способность, тг - число мест в очереди ¿-го канала, 1 е 1,к ; к>0. Оптимизация работы системы по среднему времени обслуживания заявок и минимизации вероятности отказа достигается с помощью следующего протокола диспетчеризации входных заявок [4, 5, 6].
Пусть очередная входная заявка обнаруживает систему в состоянии (х1,х2,_,хк;у1,у2,_,ук), не являющемся состоянием отказа (1,1,_,1;т1,т2,_,тк), где х.= 1, если ¿-й канал занят, х=0, если свободен, у соответствует наполненности очереди этого канала. В случае, если существует один и только один канал (с номером ¿), способный принять заявку, то заявка направляется к нему. В противном случае оптимальным считаем выбор ¿-го канала обслуживания, способного обработать заявку с минимальным средним суммарным временем Т обслуживания попавших в него заявок.
СМО с детерминированной диспетчеризацией и недетерминированной выработкой сигналов на освобождение приборов. Далее рассмотрим 2-канальную СМО с различимыми каналами пропускной способности т1>т2 без очереди Т=Т(м1,м2;0,0). Представим её поведение в дискретном времени недетерминированным конечным автоматом (НКА) К(Б,А) с алфавитом
Котенко Андрей Петрович, кандидат физико-математических наук, доцент кафедры "Математические методы в экономике". Б-таИ: ako1959@maU.ru Букаренко Максим Борисович, ассистент кафедры "Математические методы в экономике". Б-таИ: max¿m.bukarenko@gma¿l.com
внутренних состояний 5={(00),(01),(10),(11)| без
выделенных начального и конечного состояний, входным алфавитом А={0,1} и пустым выходным алфавитом. Буква 1 входного алфавита А соответствует приходу заявки в систему, а 0 - выработке сигнала освобождения заявки одним из каналов.
Рассмотрим все сочетания стохастической природы дуг орграфа, представленные (не)де-терминированными конечными автоматами.
При детерминированной диспетчеризации входных заявок получим для а=а(п), 51=51(п), ^2=^2(п) нелинейные нестационарные рекурсивные стохастические булевы функции в правой части уравнений состояний НКА К(Б,А):
(п +1) =
а © 5152 ©а5152 а
Я1 Я 2
( +1) =
а © а^ © а^ 5 2 а © 5 2 © а^1
Я1 Я 2
Тогда рекурсивную (то есть разрешимую последовательно) систему линейных рекурсивных нестохастических уравнений состояний автомата К(Б,А 1) получим из соответствующей таблицы истинности при входном сигнале
а1 (п)а2(п) = а]а2 е А1:
5 (п + 1) = а1 © а1а2 © а25 © а1а251 = = ап ©РАА 52 (п +1) = а1я1 © аха^\ © а252 © s2 © а^^ =
= Гп ©ЗА(и),_ _ _ где а := аа2, Д :=аа, г„ := аал, 3 :=а©ат.
5
Отсюда получим явный вид уравнений состояний
г-1
* (и + 1 )= ХхПД-,© * (0)ПД
X ^ ПДи-г© * 1
г=0 т=о
« г-1
* 2 (« + 1) = X Г - г П 8п - © *2 (0 )П ^ - .
г=0 т=о т=о
Здесь для универсальности обозначений по--1
лагается равенство
ПД«- := 1.
т=0
Поскольку,
П=Д = 0 , если 3к е 0,«: Дк = 0(>2 (к) = 0; П=А-г = 0, если3ке 0,и: 4 =ак©а(к)а2(к)1(к)=0,
то явные уравнения состояний примут окончательный вид:
* (и +1) =
*2 (и + 1) =
X а - г © * (0)^
— (к е 0, и ^ Дк = а1а2 = 1),
т -1
XX - г —
— (т < и)л(Д, = 0)л _
л (ук е 0, т -1 ^ Дк = а1а2 = 1);
XXX Ги - г © *2 (0)^
— (ук е 0, и ^ 8к = а2 © а1а2* = 1),
т -1
ХГи - г —
г=0
— (т < и)л
л (ук е 0, т -1 ^ 8к = 1)л л( = а2 © а1а2* =
*1(и + 1) =
а ©*1*2 ©а**2 а *1*2 ©а* ©а* а*©а* ©а**2
ад ЬР ЧР2 Ь2Р2
2 (и + 1) =
а © а* © а**2 а © *1*2 © а*1 а а © **2 © а**2
Ь Р1 ЧР Ь Р2
Построим изоморфный ДКА K(S,A2) со входным алфавитом A1={00,01,10,11}, обозначив буквой 00 сигнал освобождения канала 1 (в том числе "холостой ход" при простое (01) этого канала), буквой 01 - сигнал освобождения канала 2 (в том числе "холостой ход" при простое (10) этого канала), буквой 10 - сигнал прихода входной заявки, которую может обработать только канал 1, буквой 11 - сигнал прихода входной заявки, которую может обработать только канал 2. Отношение вероятностей появления букв 00 и 01 во входном потоке сигналов ДКА K(S,A2) равно Ч\/ Ч2 , а отношение вероятностей появления букв 10 и 11 во входном потоке сигналов равно р^р2.
Аналогично первому случаю, составим таблицу истинности детерминированных булевых функций 51(и+1) и 52(и+1) ДКА K(S,A2) и выведем из их СДНФ линейные независимые рекурсивные детерминированные уравнения состояний
* (и + 1) = а © а1а2 © а2* = а(и) © Д(и)* (и), *2 (и + 1) = **2 © а2**2 © а1а2 = у(и) © $(и)*2 (и), где а(и):= а1 а2 , Д(и):= а2 , у(и):= а1а2 , 8(и):=
а2 .
Отсюда получим явный вид уравнений состояний:
Ха(и - г)© * () — г=0 / _ х
— (Ук е 0, и ^ Д(к) = а2 (к) = 1),
т -1
Ха(и - г) —
г=0
— (т < и)л
л (Ук е 0, т -1 ^ Д(к) = а2 (к) = 1)л л (Д(т) = а2 (т) = 0);
* (и +1) =
СМО со стохастической диспетчеризацией входных заявок и выработкой сигналов на освобождение приборов. В общем случае стохастического поведения пропускной способности дуг, выходящих из вершин (00) и (11), уравнения перехода зависят от совместного распределения вероятностей переходов по этим четырём дугам. В этом случае получим рекурсивные стохастические булевы функции в правых частях уравнений состояний НКА
*2 (и +1) =
Хг(и - г )© *2 (0) —
г=0 , __,
— (Ук е 0, и ^ 8(к) = аЦ) = 1),
т -1
Хг(и - г)—
г=0
— (т < и)л
л (Ук е 0,т -1 ^ 8(к) = X) = 1)л л ((т) = а2 (т) = 0).
=0
СМО со стохастической диспетчеризацией и детерминированной выработкой сигналов на освобождение приборов. При недетерминированной диспетчеризации входных заявок с вероятностью р1<1 перехода по оптимальной дуге (00)>(10) и вероятностью р2=1-р1>0 перехода по неоптимальной дуге (00)>(01) введём детерминированную диспетчеризацию выходных заявок: ц1=1, ?2=°.
Из таблицы истинности стохастических булевых функций 51(п+1) и 52(п+1) НКА К(Б,А) с двумя стохастическими дугами (00) ^ (10) и (00) ^ (01) получим нелинейные нестационарные рекурсивные стохастические булевы функции в правой части уравнений состояний НКА К(Б,А) с двумя оставшимися недетерминированными переходами (00) ^ (10) и (00) ^ (01):
S (n + l) =
a Ф s1s2 Ф as1s2 as1 Ф as2 Ф s1s2
Р1 Р2
,(n + l) =
as1 Ф as2 Ф as^2 a
Р1 Р2
нии состоянии
sj (n +1) =
t=0 —
m—1
n
— t )Ф s1 () —
(vk e 0,n ^
p(k)=ak [ (k) ф a()¡2 (k)]=i —t ) —
t=0
— (m < п)л
A(vke 0,m — 1 ^¡3(k) = 1 л ((m) = 0);
л
s2 (n + 1) =
¿ y{n—t )Ф S2 (0) —
— (vk e 0,и ^ S(k) = a (k )a2 (k) = 1),
m—1
£r(n—t) —
t=0
— (m < п)л
л (vk e 0, m — 1 ^ S(k) = 1) A((m) = 0).
л
Тогда уравнения состоянии автомата K(S,A3) получим из таблицы истинности детерминированных булевых функции s^n+i) и s2(n+i) при входном сигнале а1 (n)a2 (n) = а1а2 e A3:
s1(n+1) = = а1а2 Ф a1s1 Ф a1a2s1 Ф s1s2 Ф Ф a1s1s2 Ф а2 s1s2 Ф a1a2s1s2 = а(п)ф /3(n)s1(n) , s2 (n +1) = a1 Ф a1a2 Ф a1a2s2 = y(n) Ф S(n)s2 (n),
где a(n) := a1 (n )a2 (n),
P(n) := a2 (n )[a1 (n) Ф a!(ñ]s 2 (n)],
y(n):= a^n))^, S(n):= a1(n>?2(n).
Таким образом, получим явныИ вид уравне-
Заключение. Выкладки аналогично переносятся на случай большего числа каналов обслуживания, наличия очередей, а также совместного стохастического поведения параметров каналов обслуживания. Несмотря на значительное усложнение неклассического орграфа СМО явный вид уравнений переходов состояний позволяет проводить имитационное моделирование как в случае простейших, так и произвольных случайных процессов поступления и обработки заявок.
СПИСОК ЛИТЕРАТУРЫ
1. Котенко А.П., Букаренко М.Б. Аналитическое описание систем массового обслуживания с использованием колец вычетов // Труды VII Всероссийской научной конференции "Математическое моделирование и краевые задачи". Самара: Изд-во СамГТУ, 2010. С. 136-140.
2. Букаренко М.Б. Система массового обслуживания с раздельными очередями к каналам // Тезисы 42-й Всероссийской конференции "Современные проблемы математики". Екатеринбург: Институт математики и механики УрО РАН, 2011. С. 11-13.
3. Комарова Е.Д., Котенко А.П. Исследование транспортной задачи линейного программирования как системы массового обслуживания // Материалы VI Международной научно-практической конференции "Научная дискуссия: вопросы технических наук". М.: Изд-во "Международный центр науки и образования", 2013. С. 17-22.
4. Котенко А.П., Букаренко М.Б. Комплекс программ имитационного моделирования работы системы массового обслуживания с неоднородными приборами и раздельными накопителями // Вестник СамГТУ. Серия "Физ.-мат. науки". 2013. №2 (31). С. 50-57.
5. Котенко А.П., Букаренко М.Б. Система массового обслуживания с различимыми каналами как конечный автомат // Вестник СамГТУ Серия "Физ.-мат. науки". 2012. №3 (28). С. 114-124.
6. Котенко А.П., Букаренко М.Б. Математическое моделирование систем массового обслуживания с ис-
S
пользованием колец вычетов в управлении органи- ние организационно-экономическими системами.
зационно-экономическими системами // Управле- Вып. 7. Самара: Изд-во СГАУ, 2010. С. 31-34.
MODELING SYSTEMS FINITE AUTOMATA QUEUING VISIBLE CHANNELS
© 2014 A.P. Kotenko, M.B. Bukarenko
Samara State Aerospace University named after academician S.P. Korolyov (National Research University)
Distinct channels queuing systems with servers of different processing capacities and separate queues are considered. In this case the results of classical queuing theory, which suppose that a state graph is a linear "birth and death" process, graph. System with distinct channels is suggested to be described by deterministic and stochastic finite automatons. This approach gives an opportunity to carry out simulation of the system behavior under different variants of entries scheduling.
Key words: queuing system, heterogeneous servers, separate queues servers, finite automaton.
Andrey Kotenko, Candidate of Physics and Mathematical Sciences, docent of Department "Mathematical methods in the economy". E-mail: ako1959@mail.ru
Maxim Bukarenko, assistant of Department "Mathematical
methods in the economy".
E-mail: maxim.bukarenko@gmail.com
