CC BY
32
МАТЕМАТИКА
УДК 519.872
ИССЛЕДОВАНИЕ RQ-СИСТЕМ, ФУНКЦИОНИРУЮЩИХ В ПОЛУМАРКОВСКОЙ СРЕДЕ
В. А. Вавилов
RESEARCH OF THE RETRIAL QUEUEING SYSTEMS OPERATING IN SEMI-MARKOV ENVIRONMENT
V. A. Vavilov
Работа выполнена при поддержке Российского Фонда Фундаментальных Исследований. Проект № 11-
01-90720-моб_ст.
В работе предлагаются математические модели систем с повторными вызовами (RQ-систем, Retrial Queueing Systems), функционирующих в полумарковской среде. Исследуются асимптотические средние характеристики, величины отклонения количества заявок в источнике повторных вызовов от их асимптотического среднего значения. Проводится глобальная аппроксимация процесса изменения числа заявок в системе и исследуется плотность распределения вероятностей значений этого процесса.
The paper is devoted to the mathematical models of retrial queueing systems concerning semi-Markov environment influence on the repeated calls source and services. The asymptotical average characteristics of the considered systems and deviation numbers of applications in the repeated calls source depending on their asymptotical average are investigated. The author provides the global approximationof the process of changing applications number in the repeated calls source, and investigates the density of distribution of the probabilities of the process values.
Ключевые слова: системы с повторными вызовами, RQ-системы, случайная среда, источник повторных вызовов, вероятностные характеристики систем массового обслуживания, метод асимптотического анализа.
Keywords: retrial queueing systems, random environment, source of repeated calls, probability characteristics of mass service system, method of asymptotic analysis.
Введение
В настоящее время исследованию систем с повторными вызовами [7] (RQ-систем, Retrial Queueing Systems) посвящается достаточно много трудов [1 - 2; 6 - 9]. Причинами являются необходимость оценки производительности и разработки новых, более эффективных локальных компьютерных сетей, радиосетей, сетей сотовой связи. Общим моментом является наличие единого передающего ресурса со случайным доступом к нему со стороны абонентских станций и, как следствие, необходимость повторной передачи в случае возникновения коллизий.
Проектирование и построение наиболее эффективно функционирующих сетей связи невозможно без учёта случайной среды, то есть неконтролируемых внешних воздействий в процессе передачи данных по каналу [1].
В данной работе проводится построение и исследование математических моделей систем с повторными вызовами, в которых изменение параметров функционирования происходит под влиянием внешнего фактора - случайной среды. Инструментом моделирования такого рода систем является аппарат теории массового обслуживания [5], позволяющий изучить вероятностно-временные характеристики функционирования систем.
Постановка задачи и построение математической модели
Рассмотрим однолинейную RQ-систему, на вход которой поступает простейший с параметром X поток заявок. Прибор этой системы может находиться в одном из двух состояний: к = 0, если он свободен;
к = 1, если он занят обслуживанием заявки. Заявка, заставшая в момент поступления прибор свободным, начинает немедленно обслуживаться. Если в течение обслуживания этой заявки другие требования на прибор не поступают, то исходная заявка по завершении обслуживания покидает систему. Если во время обслуживания одной заявки поступает новая заявка, то она переходит в источник повторных вызовов. Число заявок в источнике повторных вызовов обозначим i.
RQ-система функционирует в случайной среде. В качестве математической модели случайной среды рассмотрим полумарковский процесс s (t) с непрерывным временем t, то есть такой дискретный случайный процесс, который принимает значения из конечного множества состояний s = 1,2,...,S и для которого вложенная по моментам времени tn изменения состояний цепь s(tn ) является марковской. Времена пребывания этого процесса в различных состояниях являются условно независимыми случайными величинами, распределение вероятностей значений которых зависит лишь от номера состояния полумарковского процесса.
Для определения полумарковского процесса [4] s (t) зададим стохастическую матрицу одношаговых
вероятностей ^ переходов вложенной цепи Маркова:
Ps1s2 = P(s(tn+1) = S2 I s(tn ) = s1^
при этом будем полагать, что pss = 0.
Заметим, что
S
ZPss2 = Ь s = 1,2,...>S . (1)
•2=1
В. А. Вавилов, 2014
99
МАТЕМАТИКА
100
Также зададим набор функций распределения Gs (x) значений времени пребывания полумарковского процесса в s -м состоянии.
Будем полагать, что влияние случайной среды на функционирование сети определяется зависимостью интенсивности у обслуживания заявок в источнике повторных вызовов от состояний s случайной среды, то есть у = у(s). Вероятность обращения заявок на прибор из источника повторных вызовов за бесконечно малый промежуток времени At равна
y(s)At + o(t), при условии, что среда находится в состоянии s . Также влияние среды сказывается на интенсивности ц обслуживания заявок на приборе, то есть ц = jll(s), где s - текущее состояние случайной среды. Вероятность окончания обслуживания заявки на приборе за бесконечно малый промежуток времени At равна ^s)At + o(At).
В силу свойств приведенной математической модели трехмерный случайный вектор {к (t),i(t\s(t)} изменения во времени состояний к(t),i(t)} математической модели ссистемы и состояний {(t)} математической модели случайной среды является полумарковским процессом [4].
Для исследования описанной математической модели марковизируем процесс {к { УОХ s(t)} методом дополнительной переменной [5]. Введём переменную C(t), имеющую смысл длины интервала времени от момента t до момента смены текущего состояния случайной среды, тогда процесс изменения значений вектора {k(t),i(tX s(tX C(t)} является марковским процессом.
Обозначим
P(k(t) = k,i(t) = i, s(t) = s, C(t) < C) =
= P(k, i, s, C, t).
В любой момент времени должно выполняться
1 да S
условие нормировки: IIIP(k, i, s, даt)=1.
k=0 i=0 s=1
Представим интенсивность y(s) обращения заявок на прибор из источника повторных вызовов в виде y(s) = yc(s). Для распределения вероятностей P(k, i, s, C, t) можно составить следующую систему дифференциальных уравнений Колмогорова:
CP(0, i, s, C, t)
■( + iy^(s))P(0, i, s, C, t) =
dt
_ cp(0, i, s,c, t) cp(0,i, s,0, t)
+
CC cC
+ цф P(1, i, s, C, t) + Gs (C) iCP(0,!;cs1,0, t} Pss
s1=1
CC
CP(1, i, s, C, t) Ct
( + ,x(s))P(1, i, s, C, t) =
CP(1, i, s, C, t) CP(1, i, s,0, t)
+ Xp(0, i, s, C, t) +
+ (i + 1)yc(s) P(0, i +1, s, C, t) + + Щ1, i -1, s, C, t) +
S CP(1, i, s1,0, t)
+
Gs (C) I-
-Ps1s
(2)
s1=1 CC
Решение P(k,i,s,C,t) системы (1) будем искать методом асимптотического анализа [3] в условиях большой задержки у — 0 .
Обозначим
2
s2t = X
(3)
у = S
и рассмотрим предельный процесс
x(x) = lim (s 2i (х / s 2) I, характеризующий асим-
s—0V '
птотическое среднее нормированного числа заявок в источнике повторных вызовов.
Рассмотрим также процесс
y(x) = lim((s 2i(x / s 2) - x(x
который характеризует изменение величин отклонения нормированного числа заявок в источнике повторных вызовов от их асимптотического среднего и покажем, что он является диффузионным процессом авторегрессии. Процесс изменения состояний канала
k (х / S2) при S —— 0 является дискретным марковским процессом, независимым от процесса у(х).
Используя предельные процессы х(х) и у(х) для достаточно малых значений параметра S , рассмотрим процесс z(x) = х(х) + sy, который аппроксимирует
процесс изменения числа заявок в системе S2i(x/ s2) , и покажем, что он является однородным диффузионным процессом.
В системе (2) выполним замены 21
S i = x + sy,-P(k,i,s,C,t) = H(k, y, s,C,X,s) ,
будем иметь:
S2 CH(0,y,s,C,X,s)
(4)
-SX'(X)CH (0y.sCX-S) +
Cx Cy
+ (A,+ c(s)( x + sy)) H (0, y, s, C, X, s) = CH (0, y, s, C, X, s) CH (0, y, s,0, X, s)
" cc cc
+ M'(s)H1( y, s, C, X, s) +
S CH (0, y, s1,0, X, s) p
+ Gs (C)I
s1=1
CC
+
CC CC
Вестник Кемеровского государственного университета 2014 № 3 (59) Т. 3
2 CH(1, y, s, C, X, s) , CH(1, y, s, C, X, s) +
Cx Cy
+ (^ + a(s)( x + sy) + ц( s)) H (1, y, s, C, X, s) = CH (1, y, s, C, X, s) CH (1, y, s,0, X, s) +
cc cc
+ Ш (1, y, s, C, X, s) +
+ c(s)( x + s( y + s)) H (0, y + S, s, C, X, s) +
+ Ш (1, y-S, s, C, X, s) +
S
+
s1 s
МАТЕМАТИКА
+ g (Q £ М(1,т,.) pss. (J)
Si =1
8Q
+ Gs «) i p
S1 =1
Дальнейшие исследования будем проводить, основываясь на этой системе.
Исследование асимптотических средних характеристик RQ-систем
Под асимптотическими средними характеристиками будем понимать распределение вероятностей Rk (x) состояний канала сети связи и функцию х(т).
Теорема 1. Асимптотическое при у ^ 0 среднее значение нормированного числа заявок в источнике повторных вызовов х(т) есть детерминированная функция, определяемая обыкновенным дифференциальным уравнением вида
х' (т) = -xy0 R0 (x) + XR1 (x),
здесь Rk (x) определяются равенством
l4s)Q1( x, s, Q) =
8Q
8Qj( x, s, Q)
8Q
8Qj( x, s,0)
8Q
+ (A, + a(s) x)Q0( x, s, Q) +
+G,(Q)± 821 xs,0> Pv,
(11)
(12)
Rk (x) = lim i Qk (x, s,Q), (7)
s=1
в котором функции Qk (x, s, Q) определяются системой (11) и условием нормировки(12).
Доказательство. В системе (5) перейдем к пределу приs ^ 0 и, полагая, что существуют конечные пределы:
limH(k,у,s,Q,т,s) = H(k,у,s,Q,т), (8)
s^-0
получим систему:
8H (0, у, s, Q, т)
и условием нормировки:
iiQk(x, s, да) = 1.
k=0 s=1
Будем полагать, что Qk (x, s, Q) известны, если (6) удается решить систему (11).
Далее покажем, что x = x^) является детерминированной функцией.
В системе (5) функции Hk (у + s, s, т,.) разложим в ряд по приращениям аргумента у с точностью
до o(s), получим:
(A, + a(s) x) H (0, у, s, Q, т) = 8H (0, у, s,0, т)
8Q
8Q
+ Gs (Q) i
s1=1
+ l(s)H (1, у, s, Q, т) + 8H(0,у,st,0,т) p
l(s)H (1, у, s, Q, т) =
8H (1, у, s,0, т)
8Q ^ ’
8H (1, у, s, Q, т)
8Q
8Q
+ (A, + c(s) x) H (0, у, s, Q, т)
+ Gs (Q) i 8H (1, •8QV1,0, т) Ps1s
(9)
(X + a(s) x)Q0( x, s, Q) =-
8Q
8Q0( x, s,0)
8Q
+ l(s)Q1( x, s, Q)+
-sx'(t)8H (0- уs’Qтs) +
8у
+ (Я + a(s)( x + sy))H (0, у, s, Q, т,s) = 8H (0, у, s, Q, т,s) 8H (0, у, s,0, т,s)
" 8Q 8Q
+ l(s)H (1, у, s, Q, т,s) +
S 8H (0, у, sb0, т,s)
+
+
Gs (Q)i
s1 =1
8Q
ps1s
8H (1, у, s, Q, т,s)
- .x (т)----------------+ |(s)H(1, у, s, Q, т,s) =
+
8у
8H (1, у, s, Q, т,s) 8H (1у, s,0, т,s)
8Q 8Q
+ (Я + a(s)( x + .у)) H (0, у, s, Q, т,s) +
8
+ {xa(s) H (0, у, s, Q, т,s) - Ш (1, у, s, Q, т,s)}+
8у
s1=1 8Q
Решение H (k, у, s, Q, т) системы (9) будем искать в виде:
H (k, y, s, С т) = Qk (x, s Q)H (y, т). (10)
С учётом (10) функция Qk (x, s, Q), имеющая смысл условного совместного распределения вероятностей состояний k канала и s среды при условии x(^ = x , как следует из (9), определяется системой вида:
8Qq( x, s, Q)
+
Gs (Q i
8H (1, у, s1,0, т,s)
Ps1s + o(s). (13)
s1=1 8Q
Все уравнения системы (13) просуммируем по k и по s , учтем (1) и при Q ^ да получим:
-ex'W 8^1 iiH (k, у, s, да, т, s)1 =
8у Ik=0 s =1 J
8 Г s
= s^4 xic(s)H (0, у, s, да, т, s) -
8у l s=1
- ^iH(1, у, s, да, т, s) 1 + o(s). s=1 J
Вестник Кемеровского государственного университета 2014 № 3 (59) Т. 3
101
МАТЕМАТИКА
Поделим на s обе части полученного уравнения, выполним предельный переход (8), учтем (10), получим:
- x'
(т) YYQ(x, s, да)
к=0 s=1
дН (у, т)
ду
= \ x£a(s)Qo(x,s,да) -
s=1
-xY Qi( x, s,да)
дН (у, т)
s=1 J дУ
Учтем условие нормировки (12), обозначим:
Rk (x) = lim Y Qk (x, s, С), к = 0,1,
С^да
s=1
yRo(x) = lim Yc(s)Qo(x,s,С), (14)
^да s=1
заметим, что Rk(x) имеет смысл распределения вероятностей состояний канала, получим:
|x' (т) + yxR0 (x) - XR1 (x)}^^ = 0.
ду
Поскольку производная плотности распределения H (у, т) не может тождественно равняться нулю, следовательно, функция x = x(t) является решением обыкновенного дифференциального уравнения
x'(t) = -^xR0( x) + XR1( x), (15)
здесь Rk (x) есть распределение вероятностей состояний канала, определяемое равенством (14), в котором Qk (x, s, С) есть двумерное распределение вероятностей состояний k канала и состояний s случайной среды, которое определяется системой (11) и условием нормировки (12). Таким образом, (15) совпадает с (6). Теорема доказана.
Исследование величин отклонения числа заявок в источнике повторных вызовов от их асимптотического среднего
Докажем следующую теорему.
Теорема 2. Асимптотически при у ^ 0 случайный процесс у(т) определяется стохастическим дифференциальным уравнением вида:
dy(т) = Ax (x) у(т)л(т + B( x)dw(x), (16)
где w(t) есть стандартный винеровский процесс, A'(x) определяется производной по x от правой части дифференциального уравнения (6), а функция B( x) определяется равенством:
B2( x) = yxR0 (x) + XR1 (x) + 2{r\xh<01 (x) -
- ^h1(1)( x) + x'(T)(h01)( x) + h1(1)(x))), (17)
если выражение в правой части больше нуля, здесь параметры a и X заданы, Rk (x) определяются равенствами (7), x определяется дифференциальным уравнением (6), h^( x) определяются равенством:
lim Yh(()(x, s, С) = hk1}(x),
С^» ,
^ s=1
в котором hk\ x, s, С) есть решение системы (21).
Доказательство. Будем искать решение
Н(k, у, s, С,т, s) системы (13) в виде следующего разложения:
Н (k, ^ s, £ т, s) = Qk (x, s, С)Н (У, т) +
+ shk (у, s, С, т) + o(s). (18)
Подставим в систему (13) разложение (18), учтем (11) и запишем полученную систему, сократив на s все уравнения в следующем виде:
- (X + a(s)x)h) (у, s, С, т) + ju.(s)h1 (у, s, С, т) + dh0( у, s, С, т) dh0( у, s,0, т)
+ -
дС
+Gs (С) Y'"10
s1=1
дС
dh0( у, s1,0, т)
■ +
дС
Ps1s =
= a(s)Q0( x, s, С) уН (y, т) -дН (у, т)
- x'Cr)Q0( x, s, С)-
ду
- M-(s)h1 (у, s, С, т)) + (X + a(s) x)h (у, s, С, т) + + dht( у, s, С, т) dht( у, s,0, т) +
дС дС
S dh1 (у, sb0, т)
+Gs (С) Y
s1=1
дС
Ps,s =
= -^(s)Q0(x, s С) уН (y, т) -- (X'(т)Q1(x S, С) + °(s)xQ0 (x, s С) -
-XQ1(x, s, С))^Нду^ . (19)
ду
Будем искать решение системы (19) в следующем виде:
hk(y, s, С, т) = hkk) (x s, С) дЯдy, т) +
ду
+ hf}( x, s, С) уН (у, т). (20)
Подставим (20) в (19) и представим систему в виде двух систем:
- (X + a(s) x)^^ (x, s, С) + M.(s)/?1:1) (x, s, С) + dh01)( x, s, С) dh^1 (x, s,0)
+ -
дС
дС
+
G (r)'S dh0 ^(X, s1,0) '( )Q ( )
+ Gs(С) Y-------77-----Ps1s =- x (т)Q0( x, s).
s1 =1 дС 1
- M.(s)h1(1) (x, s, С) + (X + a(s) x)h01) (x, s, Q
) +
+ dhYCxsC) dh1(1)( x, s,0) +
дС
+Gs (С) Y
s1=1
дС
dh!(1)( x, s1;0)
дС
Ps1s =
102 | Вестник Кемеровского государственного университета 2014 № 3 (59) Т. 3
МАТЕМАТИКА
= -x'(x)Qi (х, , Q) - ^(s)xQ0 (х, , 0(21)
- (k + a(s) x)h02)( x, s, Q) + ц(.)^(2)( x, s, Q) + dh02)( x, s, Q) dh02)( x, s,0)
+
Gs (С) X
8Q
5h02)( x, si,0)
ас
+
Ps. =V(s)Qo( X, S, Q)■.
1=1 ас
- M.(s)hi(2) (x, s, Q) + (k + a(s) x)h02) (x, s, Q) + + 8h i(2)( x, s, Q) 5hi(2)( x, s,0) +
ас
ас
ax
8x
+ a fag0(x, s, С) I 8 [8Q0(x, s,0) | + ad 8x } 8Q1 8x J
+ G. (С) } Pv =ФШ x, s, С),
- ц(.) 8a<xsQ) + (k + „(.) x) 8Q0( x sC) +
8x
8x
+ 8 \ 8Qi(x,s,Q)J 8 j 8Qi(x,s,0) J +
8Q 1 8x
8Q 1 8x
?0(
+ Ц( s) H (i, y, s, Q, x, s) +
+ ^ (Q) £ 8H (0 y.C.A x-s) Psis
si=l
s2 8H(i, y, s, С, x, s)
8Q
-sx'(x)8H <■■ys-Qxs) +
+
+ G (С) X 8hP)(x, si,0) p =
+ G. (Q) nr psis
si=i 8Q i
., Q). (22)
Продифференцируем систему (ii) по x , получим: - (k + a(s)x) 8Q0(x's'С> + ц(.)ШЫЫ. +
8x 8y
+ (k + |x(s))H (i, y, s, Q, x, s) =
8H (i, y, s, Q, x, s) 8H (i, y, s,0, x, s)
" 8Q 8Q
+ (k + a(s)( x + sy))H (0, y, s, Q, x, s) +
+ kH (i, y, s, Q, x, s) +
8
+ s—{c(s)( x + sy)H (0, y, s, Q, x, s) -
8y
-kH (i, y, ., Q, x, s)}+
s2 82
+ — —y fc(s) xH (0, y, s, Q, x, s) +
2 8y
+ kH (i, y, ., Q, x, s)}+
Gs(Q)£ 8H(ixs) p + o(s2). (26)
t=i ^ '
Сложим все уравнения системы (26) по k, полу-
+
чим:
+g, (Q) X | Л. =
.“i 8Q1 8x rsis
= -a(.)Q0(x,., С). (23)
Из (22) и (23) следует, что решение h;(2)(x,s,Q) системы (22) имеет вид:
h(2)(x,.,q) = 8Qk8.Q). (24)
ox
С учетом (24) и (i0) разложение (i8) примет вид: H (к, y,., Q, x, s) = Qk (x,., Q)H (y, x) +
+ shki)( x,, С) ^8^ +
8y
+ syH(y, x) 8Qk 8’1 C) + o(s). (25)
ox
s2 ^1x H (k, y,., с, x s)|-
-sx'(x) 8- jx H (k, y, s, С, x, s)J =
=iQ{sH (k,y.,Qxs)}--^QjX H (k,y s■0■xs)}-
d
- s—{- a(s)( x + sy) H (0, y, s, Q, x, s) +
8y
+
22 s 8
kH (i, y, .,Q, x, s)}
+
+ — —у fc(s) xH (0, y, s, Q, x, s) +
+
2 8y2
+ kH (i, y, s, Q, x, s) +
i S 8H(k,y,.i,0,x,s)
G. (С) XX
k=0 .i=i
8Q
Psis + o(s ).
Теперь найдем вид функции H (y, x). Для этого Подставим в полученную систему разложение
функции в правой части системы (5) разложим в ряд функций H (k,y,., х,s) в виде (25), получим:
по приращениям аргумента y с точностью до o(s ), получим:
2 8H(0,y,s,Q, x,s) '( ) 8H(0, y,s,Q,x,s) +
s--------------------sx (x)---------------+
8x 8y
+ (k + a(s)( x + sy)) H (0, y, s, С, x, s) =
8H (0, y, s, Q, x, s) 8H (0, y, s,0, x, s)
\
( i
X Qk (x, s, С)
V k=0 ( i
8H (y, x)
- sx'(x)
X Qk (x, s, С)
V k=0
8x
8H (y, x)
8Q
+
-s2x'(x)8x jX Qk(x,.,C)
8y
8{yH (y, x)}
8Q ''8x {t=T...........J 8у
Вестник Кемеровского государственного университета 2014 № 3 (59) Т. 3
и
2
s
i03
МАТЕМАТИКА
-в2 х'(т)
f 1
Z hk1( X,^ Q)
V k=0
^8 2h (у, т)
8у2
=8т{х Hk(у, s ?т, в)[-
8? [к=0 J
-J?{|q H (уsA т,в)}-
-в(-а(5)хбо( x, s, Q)
+
+ Q x, s, gj^M-в 2 (-a(s)Qo( x, s, ?) -
ау
-„( s) x rnxso+
dx
dQi( x, s, ?) ']8{yH (у, т)}
+ X-
8x
8j
+ y [o(s) xQo( x, s, ?) + ^Qi( x, s, ?) +
+ 2((s) xh01}( x, s, ?) - ^hi(1)( x, s, 0)]8Н(2У,Т)
8у 2
1 S
+
G (?) ZZ
8H(k, у, s,,0, т, в) x 2 ч
---v ^ u’ ’ > Ps1s + о(в2). (27)
к=0 s1=1 8?
Просуммируем уравнения системы (27) по s , выполним предельный переход при ? ^ да, воспользуемся условием нормировки (12), обозначением (14), также обозначим:
lim Z hf)(x, s, ?) = hf)(x),
^ s=1
lim Z hci1)(x, s, ?) = nh01}(x)
s=1
учтем (1), получим
в 2 8H (у, т) ,(,8H (у, т)
8т
-- bx'^)-
8у
-в 2 x'(т)
Z hf( x)
V k=0
^ 8 2 H (у, т)
8у2
= -B(-yxR0(x) + Rx)) 8H(у’т)
8у
- в 2 f- VR0( x) - W 8RXX1 +
+ ?8R,( x) )_8(уЯ(у,т)} +
8x
8у
(28)
2
[yxR0 (x) + XR1 (x) -
+ 2(nxh01} (x) - ^h1(1) (x) j]8^^ + o(b2) . (29)
8у
В силу дифференциального уравнения (6) уничтожим в (29) слагаемые порядка о(в), поделим обе
2
части полученного уравнения на в , выполним несложные преобразования, будем иметь:
2
Б
+
8H (у, т) 8т
yR0 (x) - yx
8Rq( x) 8x
+
+ ?8k,(x) ^ S^H(у, т)} +
8x J 8у
1
+ —
2
(yxR0 (x) + XR1 (x) + 2 \\xhQ') (x) -
- ^h1(1) (x) + x^^h11 (x) +
+ h»( x))])8^. (30)
8у
Получили уравнение Фоккера-Планка для плотности распределения вероятностей H( у, т) значений диффузионного процесса авторегрессии у(т). Обозначим коэффициент переноса уравнения (30) как A (x). Заметим, что A (x) является производной по x от правой части дифференциального уравнения (6), то есть:
8
A'x(x) = — {-¥xRq(x) + Rx)}. (31)
8x
Коэффициент диффузии обозначим следующим образом:
Б2 (x) = yxR0 (x) + XR1 (x) + 2(qxhQ1:) (x) -
- ^h1(1) (x) + x'(т)(h01) (x) + h1(1) (x))j, (32)
если выражение в правой части больше нуля.
Получили, что (32) совпадает с (17). Из (30) следует, что H( у, т) является плотностью распределения вероятностей некоторого диффузионного процесса у(т) , который удовлетворяет стохастическому дифференциальному уравнению
dу(т) = Ax (x) у(т)Ат + B(x)dw(x), (33)
где w^) является стандартным винеровским процессом, Ax (x) определяется равенством (31), Б( x) - равенством (32), совпадающим с (17), следовательно, уравнение (33) совпадает с уравнением (16), а процесс у(т) является процессом авторегрессии.
Теорема доказана.
Следствие 2.1. Решение у(т) стохастического дифференциального уравнения (33) имеет вид:
т и
JЛ'х(x(s))ds т -JЛ'х(x(s))ds
у(т) = e0 J Б(x(u))e 0 dw(u), (34)
о
где A'(x) определяется равенством (31), Б(x) -равенством (17), x^) - дифференциальным уравнением (6), а w^) есть стандартный винеровский процесс.
Доказательство. Представим процесс у(т) в следующем виде:
т
J A'x (x(s)) ds
у (т) = e0 f (т ), (35)
тогда
104 | Вестник Кемеровского государственного университета 2014 № 3 (59) Т. 3
МАТЕМАТИКА
х
-J A X (x (^ )) ds
f (х) = e 0 y (х) , (36)
здесь A (x) определяется равенством (31), а функция х(х) - дифференциальным уравнением (6).
Продифференцируем (36), используя формулы Ито, получим:
-J A'x (x (s)) ds
df (х) = -A'x(x(x)) f (x)dx + e 0 dy .
Учитывая (35), получим:
х
-J A'x (x( s)) ds
df (х) = -A'x(x(x))e 0 У(x)dx +
х
-J A'x (x( s ))ds
+ e 0 [A'x(x(x))y(x)dx + B(x(x))dw(x)].
Выполним преобразования, будем иметь:
х
-J A'x (x(s))ds
df (х) = B(x(x))e 0 dw(x). (37)
Проинтегрируем (37), положив f (0) = 0, тогда исходное представление у(х) в виде (35) примет вид:
х и
J A'x ( x(s))ds х -J A'x ( x( s))ds
у(х) = e0 JB(x(u))e 0 dw(u), (38)
0
где A (x) определяется равенством (31), B( x) -равенством (17), x^) - дифференциальным уравнением (6), а w^) есть стандартный винеровский процесс, то есть (38) совпадает с (34). Следствие доказано.
Глобальная аппроксимация процесса изменения состояний RQ-систем в полумарковской среде
Покажем, что для достаточно малых значений параметра s случайный процесс z(х) = x^) + sy, аппроксимирующий процесс изменения числа заявок в ИПВ s2i(r / s2), является однородным диффузионным процессом. Докажем следующую теорему.
Теорема 3. С точностью до o(s) случайный процесс z^) является решением стохастического дифференциального уравнения:
dz(х) = A(z)dх + sB(z)dw(х), (39)
где w^) есть стандартный винеровский процесс, функция A(z) определяется правой частью дифференциального уравнения (6), а функция B( z) - равенством (17), то есть z^) является однородным диффузионным процессом с коэффициентом переноса A( z)
и коэффициентом диффузии s2B2(z) .
Доказательство. Поскольку z^) = x^) + sy , то дифференцируя z^) по х получаем:
dz(%) = x'(%)d% + sdy. (40)
В силу (6) и (16) имеем:
dz (х) = [-yxR0 (x) + XR (x)]dx +
д
+ sy—{- {R0 (x) +XR (x)}dх + sB(x)dw(х).
dx
Так как правая часть содержит разложение в ряд по приращениям sy аргумента x, то можно записать:
dz(x) = [-y(x + sy)R (x+sy) + XR (x + sy)]dr + + sB(z - sy)dw^).
Заметим, что z^) = x(х) + sy, тогда с точностью до o(s) имеем:
dz(%) = [-yzR (z) + XR (z)]dT +
+ sB(z)dw^) + o(s).
С учетом (6) уравнение (40) окончательно примет вид: dz(х) = A(z)dх + sB(z)dw(х) + o(s).
Таким образом, z^) является однородным диффузионным процессом с коэффициентом переноса A(z) и коэффициентом диффузии s2B2(z) и определяется с точностью до o(s) стохастическим дифференциальным уравнением вида (39). Теорема доказана.
Следствие 3.1. Плотность распределения вероятностей значений процесса z^) имеет вид:
1
F (z) =-
B 2( z)
4 )AUUL du
s20 B2(u )
B2(z)
4 )AuuX du
3s20 B2(u )
(41)
dz
где A(z) определяется правой частью дифференциального уравнения (6), B( z) - равенством (17).
Доказательство. Обозначим F(z, х) плотность распределения вероятностей значений процесса z^), тогда можно записать уравнение Фоккера-Планка для плотности этого процесса: dF (z, х) д
дх
= -—{(z) F (z, х)}+
dz
+ ■
:2 д2
2 dz2
{B 2( z) F (z, х)},
где A(z) определяется правой частью дифференциального уравнения (6), B(z) - равенством (17). Рассмотрим функционирование процесса z(х) в стационарном режиме, то есть F(z, х) = F(z), тогда стабильное распределение можно найти из уравнения Фоккера-Планка:
д
0 = -—{(z) F (z)}-
dz
22 s d
{ 2( z) F (z)}.
2 ^ ^ v-/,- (42)
2 dz2
Уравнение [42] является однородным дифференциальным уравнением второго порядка. Обозначим
B2( z)F (z) = G(z). (43)
Понизим порядок уравнения (43) и, положив константу, возникшую в результате интегрирования, равной нулю, запишем:
dG(z) 2 A(z)
dz
s2 B2 (z)
G (z).
Проинтегрируем последнее уравнение
да
1
Вестник Кемеровского государственного университета 2014 № 3 (59) Т. 3 | 105
МАТЕМАТИКА
f ЗШ du=4 f f G(u) s2 f
выполним преобразования:
2 Z A(u)
du + C1,
82J„ г 2(u) 1
lnlG (z)=h f du+4 C
s 0 B (u)
4 du
G( z) = C • eS0 B(u) .
Учтем замену (43) и перепишем последнее урав нение в виде:
2 f A(u )
F (z) = C es2 0 Bl(u)
-du
В 2( z)
(44)
Константу C найдем из условия нормировки
f F (z )dz = 1
=1,тогда
c = I f
1 sH
2 z A(u)
s20 B2(u)
du
B2( z)
(45)
Подставим (45) в (44), получим плотность распределения вероятностей для процесса z(т) в виде (41). Следствие доказано.
Таким образом, в данной работе найдено распределение вероятностей Rk (x) состояний k канала в виде (7). Получено дифференциальное уравнение (6), определяющее асимптотическое среднее значение х(т) нормированного числа заявок в источнике повторных вызовов. Исследованы величины отклонения от этого среднего, показано, что процесс их изменения у(т) определяется стохастическим дифференциальным уравнением вида (16). Найдено решение данного уравнения в виде (34). Доказано, что для достаточно малых значений параметра s случайный процесс z(t) = х(т) + sy, аппроксимирующий процесс
изменения числа заявок в системе s2i(T / s2), является однородным диффузионным процессом. Найдена важнейшая из вероятностно-временных характеристик этого процесса - плотность распределения вероятностей F (z) в виде (41).
Полученные результаты могут быть использованы при проведении анализа существующих систем с повторными вызовами, а также при проектировании новых сетей связи, реализующих более производительные протоколы передачи данных.
Литература
1. Вавилов В. А. Исследование RQ-систем с конечным числом обслуживающих приборов, функционирующих в случайной среде // Информационные технологии и математическое моделирование (ИТММ-2011): матер. X Всерос. науч.-практ. конф. с междунар. участием. Томск: Изд-во Том. ун-та, 2011. Ч. 1. С. 102 - 107.
2. Гарайшина И. Р., Моисеева С. П., Назаров А. А. Методы исследования коррелированных потоков и специальных систем массового обслуживания. Томск: Изд-во НТЛ, 2010. 204 с.
3. Назаров А. А., Моисеева С. П. Метод асимптотического анализа в теории массового обслуживания. Томск: Изд-во НТЛ, 2006. 112 с.
4. Назаров А. А., Терпугов А. Ф. Теория вероятностей и случайных процессов. Томск: Томск: Изд-во НТЛ,
2006. 204 с.
5. Назаров А. А., Терпугов А. Ф. Теория массового обслуживания. Томск: Изд-во НТЛ, 2010. 228 с.
6. Судыко Е. А., Назаров А. А. Исследование марковской RQ-системы с конфликтами заявок и простейшим входящим потоком // Вестник Томского государственного университета. Управление, вычислительнаятехни-каиинформатика. 2010. № 4(12). С. 79 - 90.
7. Artalejo J. R., Gomez-Corral A. Retrial Queueing Systems: a Computational Approach // Springer. 2008. 309 p.
8. Dudin A., Klimenok V. Queueing System BMAP|G| 1 with Repeated Calls // Mathematical and Computer Modelling. 1999. № 30. P. 115 - 128.
9. Falin G. I. A diffusion approximation for retrial queueing systems // Theory of Probability and Its Application.
1991. Vol. 36. № 1. P. 149 - 152.
Информация об авторе:
Вавилов Вячеслав Анатольевич - кандидат физико-математических наук, доцент, доцент кафедры информатики и математики филиала КемГУ в г. Анжеро-Судженске, vavilovv@yandex. ru.
Vyacheslav A. Vavilov - Candidate of Physics and Mathematics, Associate Professor, Assistant Professor at the Department of Computer Science and Mathematics, Anzhero-Sudjensk Branch of Kemerovo State University.
Статья поступила в редколлегию 28 июля 2014 г.
106 | Вестник Кемеровского государственного университета 2014 № 3 (59) Т. 3
