CC BY
15
ВЕСТНИК ТОМСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА
2016 Управление, вычислительная техника и информатика № 3 (36)
УДК 519.21
DOI: 10.17223/19988605/36/7
М.Н. Сиротина
СРАВНЕНИЕ МП- И ММ-ОЦЕНОК ДЛИТЕЛЬНОСТИ МЕРТВОГО ВРЕМЕНИ В МОДУЛИРОВАННОМ СИНХРОННОМ ДВАЖДЫ СТОХАСТИЧЕСКОМ
ПОТОКЕ СОБЫТИЙ
Работа выполнена в рамках государственного заказа Минобрнауки РФ на проведение научных исследований в Национальном исследовательском Томском государственном университете на 2014-2016 годы.
Рассматривается модулированный синхронный дважды стохастический поток событий. Поток функционирует в условиях непродлевающегося мертвого времени, т.е. после каждого зарегистрированного события наступает время фиксированной длительности, в течение которого другие события исходного модулированного синхронного потока недоступны наблюдению. Полагается, что длительность мертвого времени - неизвестная величина. Проводится сравнение качества получаемых (по наблюдениям за моментами наступления наблюдаемых событий потока) оценок длительности мертвого времени методом максимального правдоподобия (МП-оценки) и методом моментов (ММ-оценки).
Ключевые слова: модулированный синхронный поток событий; непродлевающееся мертвое время; МП-оценки; ММ-оценки; длительность мертвого времени.
Настоящая работа является непосредственным продолжением исследований модулированного синхронного потока событий, начатых в статьях [1-4].
Математические модели систем массового обслуживания широко применяются при описании реальных физических, технических и других процессов и систем. В связи с бурным развитием компьютерной техники и информационных технологий появилась важная сфера приложений теории массового обслуживания - проектирование и создание информационно-вычислительных сетей, компьютерных сетей связи, спутниковых сетей, телекоммуникационных сетей, объединенных термином «цифровые сети интегрального обслуживания» (ЦСИО) [5, 6].
На практике параметры, определяющие входящий поток событий, изменяются со временем, при этом изменения часто носят случайный характер, последнее приводит к рассмотрению дважды стохастических потоков событий. По-видимому, одной из первых работ в этом направлении явилась статья [7], где дважды стохастический поток определяется как поток, интенсивность которого есть случайный процесс. Дважды стохастические потоки можно разделить на два класса: к первому классу относятся потоки, интенсивность которых есть непрерывный случайный процесс; ко второму классу относятся потоки, интенсивность которых есть кусочно-постоянный случайный процесс с конечным числом состояний. Подчеркнем, что потоки второго класса впервые введены в рассмотрение практически одновременно в 1979 г. в [8, 9]. В [8] введенные потоки названы MC (Markov сЬат)-потоками, в [9] - MVP (Markov versatile ргосе88е8)-потоками. Отечественные и зарубежные авторы в своих работах начала 90-х гг. [1015] называют данные потоки событий либо дважды стохастическими потоками событий, либо MAP-потоками, либо MC-потоками. В свою очередь, в зависимости от того, каким образом происходит переход из состояния в состояние, MC-потоки можно разделить на три типа: 1) синхронные потоки событий [16-21]; 2) асинхронные и обобщенные асинхронные потоки событий [22-27]; 3) полусинхронные и обобщенные полусинхронные потоки событий [28-33]. В [34] введены в рассмотрение MAP-потоки событий первого порядка (собственно MAP-потоки, введенные в [9]) и MAP-потоки событий второго порядка (суперпозиция двух синхронизированных MAP-потоков первого порядка, отличающихся друг от друга исходными параметрами). В [34] показывается, что синхронный MC-поток является частным случаем MAP-потока первого порядка, асинхронный и полусинхронный MC-потоки являются частными случаями MAP-потока второго порядка.
В реальных ситуациях параметры, задающие входящий поток событий, известны либо частично, либо вообще неизвестны, либо (что еще более ухудшает ситуацию) изменяются со временем. Вследствие этого возникают задачи: 1) оценки состояний потока (задача фильтрации интенсивности потока) по наблюдениям за моментами наступления событий [24, 27, 29, 33, 35-38]; 2) оценка параметров потока по наблюдениям за моментами наступления событий [16-22, 26, 28, 30-32].
Одним из искажающих факторов при оценке состояний и параметров потока событий выступает мертвое время регистрирующих приборов [39-46], которое порождается зарегистрированным событием. Другие же события, наступившие в течение периода мертвого времени, недоступны наблюдению (теряются). Можно считать, что этот период продолжается некоторое фиксированное время (непродле-вающееся мертвое время). В качестве примера приведем протокол С8МЛ/СБ - протокол случайного множественного доступа с обнаружением конфликта, широко используемого в компьютерных сетях. В момент регистрации (обнаружения) конфликта на входе некоторого узла сети по ней рассылается сигнал «заглушки» («пробки»); в течение времени рассылки сигнала «заглушки» заявки, поступившие в данный узел сети, получают отказ в обслуживании и направляются в источник повторных вызовов. Здесь время, в течение которого узел сети закрыт для обслуживания заявок, поступающих в него после обнаружения конфликта, можно трактовать как мертвое время прибора, регистрирующего конфликт в узле сети. Для того чтобы оценить потери заявок в узле сети, необходимо оценить длительность мертвого времени, которым выступает в данном случае длительность сигнала «заглушки».
В работах [1-4] введен в рассмотрение модулированный синхронный поток событий, являющийся обобщением синхронного потока и относящийся к классу МЛР-потоков второго порядка. Синхронный поток событий систематически исследовался в работах [16-21, 47-52]. В статье [4] приведено решение задачи оценивания длительности мертвого времени методом максимального правдоподобия [53-59]. В настоящей статье предложено решение задачи оценивания длительности мертвого времени методом моментов, а также производится сравнение оценок, полученных методом максимального правдоподобия (МП-оценки) и методом моментов (ММ-оценки).
1. Постановка задачи
Рассматривается модулированный синхронный поток событий (далее - поток), интенсивность которого есть кусочно-постоянный стационарный случайный процесс X(?) с двумя состояниями: X2 > X2 > 0). Длительность пребывания процесса X(?) (потока) в 7-м состоянии распределена по экспоненциальному закону с параметром а7, 7 = 1,2 . Если процесс X(?) в момент времени ? находится в 7-м состоянии, то на полуинтервале [?, t + Д?), где Д? - достаточно малая величина, с вероятностью а7Д? + о(Д?) пребывание процесса X(?) в 7-м состоянии закончится и процесс X(?) с вероятностью, равной единице, перейдет из 7-го состояния в /-е (7, / = 1,2; 7 ф/). В течение временного интервала случайной длительности, когда X(?) = X 7, имеет место пуассоновский поток событий с интенсивностью X 7, 7 = 1,2. Кроме того, переход из первого состояния процесса X(?) во второе возможен в момент наступления события пуассоновского потока интенсивности X!; переход осуществляется с вероятностью р (0 < р < 1); с вероятностью 1-р процесс X(?) остается в первом состоянии. Переход из второго состояния процесса X (?) в первое возможен также в момент наступления события пуассоновского потока интенсивности X2; переход осуществляется с вероятностью q (0<д<1); с вероятностью 1-д процесс X (?) остается во втором состоянии. В сделанных предпосылках X(?) - марковский процесс.
Матрицы инфинитезимальных коэффициентов рассматриваемого потока примут вид
(1 - Р^ pX1 II
- (X, + а,) а,
А = 1
а 2 - (X 2 + а 2)
А =
qX2 (1 - q)X2
Элементами матрицы А1 являются интенсивности переходов процесса X (?) из состояния в состояние с наступлением события. Недиагональные элементы матрицы А0 - интенсивности переходов из
состояния в состояние без наступления события. Диагональные элементы матрицы О0 — интенсивности выхода процесса X (г) из своих состояний, взятые с противоположным знаком. Отметим, что если а г = 0, г = 1,2, то имеет место обычный синхронный поток событий [16—21].
После каждого зарегистрированного события наступает время фиксированной длительности Т (далее - мертвое время), в течение которого другие события потока недоступны наблюдению. Рассматривается непродлевающееся мертвое время, т.е. события, наступившие в течение интервала мертвого времени не вызывают его продления. По окончании длительности периода мертвого времени первое наступившее событие вновь порождает период мертвого времени длительности Т и т.д. Вариант возникающей ситуации приведен на рис.1, где Х1,X2 - состояния процесса X(г), г2,... - моменты наступления наблюдаемых событий потока, штриховка - периоды мертвого времени длительности Т, ось под номером 1 отображает исходный модулированный поток событий, под номером 2 - схему создания мертвого времени, под номером 3 - наблюдаемые события модулированного синхронного потока, функционирующего в условиях мертвого времени.
40
1
2
1-Р
1-Р
1-ч
-ЫШИщЩ—ЫППШЩ
' ЛГ 1 ЛГ
1-Я 1-я
-уиннщ
1-Р
1 Р
-ЦШПМЦ
Рис. 1. Формирование наблюдаемого потока событий
я
г
г
г
Процесс X(г) и типы событий (события пуассоновских потоков с интенсивностями Х1 либо X 2) являются принципиально ненаблюдаемыми (в сделанных предпосылках X (г) - скрытый марковский процесс), а наблюдаемыми являются только моменты наступления наблюдаемых событий на временной оси г1,г2,.... Последовательность моментов г1,г2...,гк,... наступления событий наблюдаемого потока порождает вложенную цепь Маркова {X(tk)}, т.е. поток обладает марковским свойством, если его эволюцию рассматривать с момента гк (момент наступления события потока), к = 1,2,... . Рассматривается стационарный режим функционирования наблюдаемого потока событий, поэтому переходными процессами на интервале наблюдения (г0, г), где г0 - момент начала наблюдений, г - момент окончания наблюдений, можно пренебречь, при этом полагается г0 = 0 .
Главная задача работы - в момент окончания наблюдений (в момент времени г) построить оценки Т длительности мертвого времени методом максимального правдоподобия и методом моментов (в предположении, что остальные параметры потока X!, X2, ах, а2, р, я известны) и произвести сравнение получаемых оценок.
2. МП-оценка длительности мертвого времени
Обозначим %к = гк+1 - гк, к = 1,2,..., - значение длительности к-го интервала между моментами наступления соседних событий наблюдаемого потока (рис. 2).
-т-*—Г
Рис. 2. Интервал между соседними наблюдаемыми событиями потока
Поскольку поток функционирует в стационарном режиме, то плотность вероятности значений длительности интервала между моментами наступления соседних событий наблюдаемого потока Рт (т к) = РТ (т),т — 0, для любого к (индекс Т подчеркивает, что плотность вероятности зависит от длительности мертвого времени). Таким образом, без потери общности момент наступления события ?к можно положить равным нулю, т.е. т = 0 .
Тогда плотность вероятностей рт (т) значений длительности интервала между наступлениями соседних событий модулированного синхронного потока, функционирующего в условиях мертвого времени, запишется в виде [4]:
[0,0 < т < Т,
РТ (т) = 1 i ТЧ i тч (1)
[у(Т)г1е-21(х-Т) + (1-у(Т))г2е-22(т-Т), т — Т,
У(т) =-1-(г2 ) ^2п2 (Т)),
29 - 2
2Х = ((ах + А1 +а 2 + А2) -
а1 + А1 -а 2 -А 2 ) + 4аха 2)/2,
2 2 = ((а1 + А1 +а 2 + А 2) + >/( а1 + А1 -а 2 -А 2) + 4а1а 2)/2, 21 < 22,
п,(Т) = П, -(П, -П,(0| Т))е-(а1+а2+pXl2)Т, п2(Т) = п2 -(п2 -п2(0| Т))е-(а1+а2+pXl+qX2)Т, qX2а1 + (1 - p)Xl (X2 + а2) - ^2(1 - р - q)(п2 + п1е-(а1+а2+pXl+qX2)Т)
п1(0|Т) = п 2(0|Т) =
2122 - X1X2 (1 - р - q)e-(а1+а2+pXl+qX2)Т (1- q)X2(X1 + а1) + pX1а2 -X1X2(1- p-q)(п1 + п2е-(а1+"2+pXl+^2)Т) 2122 - X1X2 (1 - р - q)e-(а1+а2+^+qX2)T :
а 2 + qX 2 ах + pX1 п, =-, п2 =
1 _-' "'2
ах + pX1 + а 2 + qX 2 ах + pX1 + а 2 + qX 2
Пусть т, = ?2 - ?1, т2 = ?3 - ?2,..., тк = ?к+1 - ?к, к = 1,2,..., - последовательность измеренных в результате наблюдения за потоком на интервале наблюдения (?0, ?) значений длительностей интервалов между соседними событиями наблюдаемого потока. Упорядочим величины т,, т2,...,тк по возрастанию: тт1п = т(1) < т(2) < ... < т(к). Тогда функция правдоподобия, с учетом (1), примет вид
ДА7,а7,р,q,T | т(1),...,т(к}) = 0,0 < Тт1П <Т,
ДА7,а7,р,q,T|т(1),.,т(к}) = П рт(т(;}),0 < Т <ТтШ.
/=1
Так как часть задачи исследования заключается в построении оценки длительности мертвого времени Т, то согласно методу максимального правдоподобия ее решение есть нахождение глобального максимума функции правдоподобия по неизвестному параметру Т:
Щ | т(1), т(к}) = П рт (т(;}) = П {у(Т)2!е"21(т(')-Т} + (1 - у(Т))22е-2^')}Ц тах,0 < Т < V», (2)
/=1 /=11 ) т
где 2,, 22, у(Т) определены в (1).
Значение T, при котором (2) достигает своего глобального максимума, есть МП-оценка T длительности мертвого времени.
В [5] аналитически строго решена оптимизационная задача (2): при любых значениях параметров рассматриваемого потока событий Xj > 0, X2 > 0 (Xj > X2), ах > 0, а2 > 0, 0 < p < 1, 0 < q < 1 МП-оценка
длительности мертвого времени ТМП = т min . Таким образом, в процессе наблюдения потока событий (в
течение временного интервала (t0, t)) вычисляются величины т к, k = 1, n , после чего находится
Т min = min Т k , k = 1 n и полагается ТМП = Т min .
3. ММ-оценка длительности мертвого времени
В [3] показано, что модулированный синхронный поток событий, функционирующий в условиях мертвого времени, в общем случае является коррелированным потоком. Только в частных случаях поток становится рекуррентным.
Пусть т k = tk+j - tk, т k+i = tk+2 - tk+j, k = 1,2,..., - значения длительностей смежных k-го и k+1-го интервалов между соседними событиями наблюдаемого потока (k = 1,2,...). В силу стационарности потока можно положить k = 1 и рассматривать длительности интервалов тх = 12 - tl, т 2 = t3 -12, тх > 0, т 2 > 0. Тогда тх = 0 соответствует моменту tl наступления события наблюдаемого потока, а т 2 = 0 соответствует моменту 12 наступления следующего события наблюдаемого потока. При этом соответствующая совместная плотность вероятностей имеет вид pT (т k, т k+1) = pT (т j, т 2), т 1 > 0, т 2 > 0 [2]:
pT (тх, т2) = 0,0 < ^ <T, 0 < т2 <T,
pt (т1, т2) = PT 0i)PT (т2) +
+CT (z1e-z1(X1-T} - z2e-z2(X1-T})(z1e-z1(x2-T} - z2e-z2(x2-T>),T1 > T, x2 > T, c = e-(ß1+ß,)T P(X -X2)(ß1^1(0)-ß2^2(0))(P + P2) T [(z2 - z1)(zz2 - Pe-(ß1+ß2)T )(ß1 + ß2)"
x{z1z2 - e-(ß1+ß2)T (2z1 z2 - (ß +ß2)(z1 + z2)) + e-2(ß1+ß2)T (z1z2 -(ß1 +ß2)(^1(1 - p) + M1 -q)))}, (3) ß1 =a1 + pXb ß2 =a2 + q^2, P = ^(1 - p - q),
P1 = (1 - p) A1a2 + qX2( X1 +a1), P2 = (1 - q) X 2a1 + pX1( X 2 +a2), Я1 (0) = P / (P1 + P2 ), (0) = P2 / (P1 + P2 ), где z1, z2, pT (тk) определены в (1) для т = тk, k = 1,2 .
Теоретическая ковариация значений т1, т 2 запишется в виде
2
2x
J тpT (т)dт
T
(4)
СОУ(х1; X2) = | |Т^Т2Pт (т1, Т2)<яМТ2 -
TT
Подставляя плотность (1) и совместную плотность (3) в интегральную формулу для ковариации (4), находим ее явный вид
еоу(т1, т 2) =
z 2 z1 V z1 z2 J
CT , (5)
где z1, z2 определены в (1), CT определена в (3).
Пусть за время наблюдения (в течение временного интервала (^0, 1)) реализовалось п интервалов (tk, tk+1), k = 1,п, длительности тк, k = 1, п . Введем статистику:
Л 1 п—1 ( 1 П \2
С°У( Т1, Т 2 ) =--2 ТkТk+1 — 1-Х Тk I , (6)
П — 1 k=1 I nk=1 )
2
которая представляет собой оценку теоретической ковариации (5). Тогда согласно методу моментов уравнение моментов, учитывающее коррелированность потока событий, может быть записано в виде
/ - ^ 2 Л 2 1 Ст = еоу(т1, т2). (7)
V ZiZ2 )
Подставляя в (7) значение Ст , определенное в (3), вводя новую переменную х = е Ф^Р2^ и проделывая необходимые выкладки, выражение (7) преобразуется в кубическое уравнение относительно переменной х вида
ах3 + Ьх 2 + сх + ё = 0, а = И[22 -(Р1 + в2 xx1с - р) + X2(1- q))]
Ь = -{¿[2212^(Р1 + в2 )(21 + 22)] + Р2оОУ(т,, т2)}, (8)
Л
И + 2Р соу(ть т2)
с — zi z 2
л
, d — -(ziz 2)2cov(ii, X 2),
h — -X 2)(ßiPi -ß 2^2)/(Zl Z 2 (ßl + ß 2 )) 2 • При решении уравнения (8) определяются три корня xi, i — i,2,3 , которые, в свою очередь, задают три ММ-оценки длительности мертвого времени:
Тмм —--i-ln xt, i —1,2,3 .
ßi + ß 2 '
Используя полученные выше ММ-оценки Т^, i — i, 2,3, можно определить единственную ММ-оценку ТММ по алгоритму, описанному ниже:
1. Для заданного набора параметров Xi, X 2, ai, а 2, p, q, Т осуществляется в течение длительности
Tm — t ед. времени интервала наблюдения (интервала моделирования) имитационное моделирование наблюдаемого потока событий.
л
2. В ходе имитационного моделирования находится оценка ковариации cov(xi, т2), определенная в
(6), где n > 2 (n - количество наблюдаемых событий потока).
3. Используя метод Виета, решается кубическое уравнение (8), т.е. находятся три корня xi, i — i,2,3 , ко-
торые могут быть как вещественными, так и комплексными.
4. Если все корни x{, i — i, 2,3, - комплексные, то ТММ — т min .
5. Если среди трех корней xi, i — i, 2,3, выделяется хотя бы один вещественный, то возможны три варианта:
5.i. Вещественный корень один - xi, тогда:
5.П. Если xi < 0, то Тмм — т min. 5.i.2. Если xi > 0, то:
5.1.2.1. ТММ — т min , если ТММ > Tmin .
5.1.2.2. ТММ — ТМЫ , если 0 < Тми < Tmin . 5Л.2.3. ТММ — Тmin , если ТММ < 0.
5.2. Вещественных корня два - xi, x2 (xi < x2), тогда:
5.2.L Если xi < x2 < 0 , то ТММ — Tmin .
5.2.2. Если xi < 0 < x2, то:
5.2.2.1. ТММ — тmin , если Тмм > Tmin . 5 222. ТММ — Тми , если 0 < TMM < Tmin .
5.2.2.3. ТММ = Tmin , если TMM < 0.
5.2.3. Если 0 < Xj < x2, то:
5.2.3.1. ТММ =Т min, если Tmin < TMM < TMM ■ 5232. Тмм = TMM , если 0 < TMM < Tmin < TMM .
5.2.3.3. Тмм = 1 min , если TMM < 0 < Tmin < TMM ■
5.2.3.4. TMM = (TMM + TMM>)/ 2 , если 0 < TMM < TMM < Tmin . 5-2-3.5. TMM = TMM , если TMM < 0 < TMM < Tmin .
5.2.3.6. Tmm = Tmin , если TMM < TMM < 0 .
5.3. Вещественных корня три - Xj, x2, x3 (Xj < x2 <x3), тогда:
5.3.1. Если Xj < X 2 < X3 < 0, то Tmm = T min.
5.3.2. Если Xj < X2 < 0 < x3 , то:
5.3.2.j. TMM = Tmin , если TMM > Tmin .
5 322. TMM = TMM , если 0 < TMM < Tmin .
5.3.2.3. TMM = Tmin , если TMm < 0.
5.3.3. Если Xj < 0 < X2 < x3 , то:
23
5.3.3.j. TMM =T min, если Tmin < TMM < TMM .
5.3.3.2. TMM = TMM , если 0 < TMM < Tmin < TMM .
5.3.3.3. TMM = Tmin , если TMM < 0 < Tmin < TMM ■
5.3.3.4. TMM = (TMM + TMM)/ 2 , если 0 < TMM < TMM < Tmin . 5-3-3.5. TMM = TMM , если TMM < 0 < TMM < Tmin .
5.3.3.6. TMM = Tmin , если TMm < TMM < 0 .
5.3.4. Если 0 < Xj < X2 < x3 , то:
T min , если Tmin < TMM < TMM < TMM .
Ä(3) если 0 < f(3) <T < f(2) < f(1)
MM MM min MM MM
г если T(3) < 0 <T < T(2) < T(1) min MM min MM MM
7"(2) + 7Ч3) )/2 если 0 < Т(3) < Т(2) <T < Т(1)
MM MM MM MM min MM
"(2) если 743) < 0 < f(2) <T < 741)
MM MM MM min MM
г если Т(3) < Т(2) < 0 <T < Т(1) min MM MM min MM
741) + 742) + 743) )/3 если 0 < 743) < 742) < f(1) <T
MM MM MM MM MM MM min
7"(1) + 742) ) / 2 если Т(3) < 0 < Т(2) < Т(1) <T
MM MM MM MM MM min
"(1) если 7^(3) < Г® < 0 < 7^(1) <T
MM MM MM MM min
T если т(3) < т(2) < г(1) < 0 <T
min MM MM MM min
В результате работы алгоритма выбирается один из вышеописанных вариантов, который определяет единственную ММ-оценку TMM длительности мертвого времени.
4. Численное сравнение МП- и ММ-оценок
Для получения статистических результатов разработан алгоритм вычисления МП- и ММ-оценок, который реализован на языке программирования C# в среде программирования Microsoft Visual Studio C#.
5.3.4.1. TMM =
5.3.4.2. =
5.3.4.3. TMM =
5.3.4.4. =
5.3.4.5. 7JMM =
5.3.4.6. TMM =
5.3.4.7. 7JMM =
5.3.4.8. 7JMM =
5.3.4.9. =
5.3.4.10 . =
Программа расчета состоит из трех этапов. На первом этапе осуществляется имитационное моделирование рассматриваемого модулированного синхронного дважды стохастического потока событий, функционирующего в условиях мертвого времени (при заданных параметрах моделирования Xi, X 2, ai, а 2, p, q, Т). Сам алгоритм имитационного моделирования в данной работе не приводится; хотя он достаточно трудоемок, но никаких принципиальных сложностей не представляет. Результатом работы имитационной модели является последовательность значений длительности интервалов между моментами наступления соседних событий наблюдаемого потока ti,т2,...тn , n — 2,3,...(тn — tn+i -tn, где tn+i - последнее наблюдаемое событие на интервале моделирования Tm).
Второй этап расчета заключается в нахождении численных значений МП- и ММ-оценок рассматриваемого потока. Опишем сокращенный алгоритм программной реализации второго этапа расчета:
1. На основании последовательности ti,т2,...тп ,n — 2,3,., находится МП-оценка
ТМП — тmin ( тmin — min тk , k — ^ n ).
2. Вычисляется оценка ковариации (6).
3. Методом Виета находятся корни xi, i — i, 2,3, кубического уравнения (8), на основании ко-
торых вычисляются три ММ-оценки ТММ, i — i, 2,3 .
4. Реализуется алгоритм нахождения единственной ММ-оценки ТММ .
Третий этап расчета заключается в сравнении качества МП- и ММ-оценок путем проведения статистического эксперимента:
1. Для заданного набора параметров Xi, X 2, ai, а 2, p, q, Т осуществляется моделирование
наблюдаемого потока событий на интервале моделирования Tm для каждого j-го эксперимента, j — i, N .
2. В каждом j-м эксперименте вычисляются значения МП- и ММ-оценок ТМП, TMM , j — i, N .
3. В каждом j-м эксперименте вычисляются значения величин АТМП — (ТМП - Т)2,
АТММ — (ТММ -Т)2, j — i, N, где Т - истинное значение мертвого времени, известное из имитационной модели.
4. Осуществляется повтор шагов i-3 N раз.
Результатом третьего этапа расчета являются две выборки: ^АТ\Мтг,АТм77,."АТмП), ^А ТММ, А ТММ"А ТММ ), используя которые вычисляются выборочные вариации полученных оценок:
N „ N
Умп — (i/N) S шМП> Умм — (i/N) ЕАТММ j—i j—i
Сравнивая значения выборочных вариаций, решается, какая из оценок при заданных параметрах Xi, X2, ai, а2, p, q, Т лучше, а какая хуже: если !/МП < VMM, то МП-оценка лучше ММ-оценки, в противном случае ММ-оценка лучше МП-оценки. Стоит отметить, что по определению при заданном времени моделирования Тт МП-оценка будет всегда смещенная (т min > Т). Вариант, когда МП-оценка несмещенная, реализуется только в асимптотическом случае, когда Тт ^ да .
Результаты статистического эксперимента приведены в табл. i—6. В первой строке таблиц указана длительность имитационного моделирования Тт, причем Тт изменяется от i0 до 50 с шагом i0 ед. времени в табл. i—3, в табл. 4—6 Тт изменяется от 600 до i000 с шагом i00 ед. времени. Во второй и третьей сроках таблиц для каждой длительности имитационного моделирования Тт приведены численные значения для выборочных вариаций МП- и ММ-оценок УМП и VMM соответственно. В четвертой строке
каждой из таблиц приведены численные значения разности !/МП - VMM . Численные результаты во всех таблицах получены для N = i00 (N - количество экспериментов).
Т а б л и ц а 1
Результаты эксперимента при X1 = 2, X 2 = 1, р = 0,5, q = 0,4, а 1 = 0,03, а 2 = 0,03, Т = 0,5
т т 10 20 30 40 50
V у МП 0,08936555 0,02139580 0,00190583 0,00114900 0,00128956
V у ММ 0,08919526 0,02131419 0,00190530 0,00114900 0,00445481
V -V у МП у ММ 0,00170293 0,00008161 0,00000052 0 -0,00316525
Т а б л и ц а 2
Результаты эксперимента при X1 = 1, X 2 = 0,5, р = 0,4, q = 0,2, а 1 = 0,05, а 2 = 0,05, Т = 0,7
т т 10 20 30 40 50
V у МП 0,90196699 0,12758867 0,040312912 0,02378204 0,00815925
V ММ 0,90034098 0,12616512 0,038268349 0,02355205 0,00812320
V -\г у МП ММ 0,00162600 0,00142355 0,002044563 0,00229987 0,00003605
Т а б л и ц а 3
Результаты эксперимента при X1 = 1,5, X 2 = 0,9, р = 0,4, q = 0,3, а 1 = 0,02, а 2 = 0,02, Т = 0,2
т т 10 20 30 40 50
V у МП 0,04937792 0,00691690 0,00299561 0,00132811 0,00104303
V ММ 0,04930285 0,00690571 0,00296747 0,00126077 0,00103060
V -\г у МП ММ 0,00007507 0,00001118 0,00002814 0,00006734 0,00001243
Т а б л и ц а 4
Результаты эксперимента при X1 = 2, X 2 = 1, р=0,5, q=0,4, а 1 = 0,03, а 2 = 0,03, Т = 0,5
т т 600 700 800 900 1000
V у МП 0,00000410 0,00000405 0,00000304 0,00000225 0,00000151
V ММ 0,00314312 0,00460004 0,00594022 0,00305194 0,00614238
V -\г у МП ММ -0,00313901 -0,00459599 -0,00593717 -0,00304969 -0,00614086
Т а б л и ц а 5
Результаты эксперимента при X1 = 1, X 2 = 0,5 , р = 0,4, д = 0,2, а 1 = 0,05 , а 2 = 0,05 , Т = 0,7
т т 600 700 800 900 1000
V у МП 0,00009986 0,00006699 0,00003027 0,00003759 0,00002732
V ММ 0,00875014 0,00404040 0,00710656 0,00359553 0,00014134
V -\г у МП ММ -0,00865027 -0,00397328 -0,00707628 -0,00355793 -0,00011401
Т а б л и ц а 6
Результаты эксперимента при X1 = 1,5, X 2 = 0,9 , р = 0,4, д = 0,3, а 1 = 0,02, а 2 = 0,02, Т = 0,2
т т 600 700 800 900 1000
V у МП 0,00000638 0,00000349 0,00000335 0,00000166 0,00000137
V ММ 0,00000638 0,00000349 0,00103220 0,00002244 0,00002016
V -\г у МП ММ 0 0 -0,00102884 -0,00002078 -0,00001879
Анализ результатов, полученных в ходе статистического эксперимента, показывает, что при достаточно малом времени моделирования Tm (Tm = 10, 20, 30, 40,50 ед. времени) ММ-оценки длительности мертвого времени лучше МП-оценок длительности мертвого времени (см. табл. 2, 3) либо, по крайней мере, не хуже их (см. табл. 1). Это вполне объяснимо, так как при малом времени наблюдения МП-
оценка ТМП может быть достаточно сильно смещенной относительно Т. При достаточно большом времени моделирования Tm (Tm = 600, 700, 800, 900, 1000 ед. времени) МП-оценки лучше ММ-оценок (табл. 4-6), что вполне естественно, так как при увеличении времени наблюдения за потоком смещение МП -оценки Тмп относительно Т уменьшается.
Заключение
Анализ проведенного исследования оценок длительности мертвого времени Т на базе метода максимального правдоподобия и метода моментов позволяет сделать вывод о том, что при малом времени
наблюдения за потоком предпочтительнее применять метод моментов и использовать ММ-оценку ТММ , при большом времени наблюдения за потоком - МП-оценку ТМП. Границу применения той или иной оценки при заданных параметрах потока Х1, X 2, а1, а 2, p, q, T следует определять для каждого случая отдельно путем имитационного моделирования.
ЛИТЕРАТУРА
1. Сиротина М.Н. Оптимальная оценка состояний модулированного синхронного дважды стохастического потока событий в
условиях непродлевающегося мертвого времени // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2014. № 1(26). C. 63-74.
2. Сиротина М.Н., Горцев А.М. Плотность вероятностей длительности интервала между соседними событиями в модулирован-
ном синхронном потоке при непродлевающемся мертвом времени // Материалы Международной научной конференции, посвященной 80-летию профессора, доктора физико-математических наук Геннадия Алексеевича Медведева. Белоруссия. Минск. 23-26 февраля 2015. Минск : Изд-во РИВШ, 2015. С. 299-304.
3. Сиротина М.Н. Совместная плотность вероятностей значений длительности интервалов между моментами наступления со-
седних событий в модулированном синхронном дважды стохастическом потоке при непродлевающемся мертвом времени и условия рекуррентности потока // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2015. № 2(31). C. 53-67.
4. Сиротина М.Н., Горцев А.М. Оценка максимального правдоподобия длительности мертвого времени в модулированном
синхронном дважды стохастическом потоке событий // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2016. № 1(34). C. 50-64.
5. Дудин А.Н., Клименок В.Н. Системы массового обслуживания с коррелированными потоками. Минск : Изд-во БГУ, 2000.
175 с.
6. Башарин Г.П., Гайдамака Ю.В., Самуйлов К.Е. Математическая теория телетрафика и ее приложения к анализу мультисер-
висных сетей связи следующих поколений // Автоматика и вычислительная техника. 2013. № 2. С. 11-21.
7. Kingman J.F.C. On doubly stochastic Poisson process // Proceedings Cambridge Phylosophical Society. 1964. V. 60, No. 4. P. 923-
930.
8. Башарин Г.П., Кокотушкин В.А., Наумов В.А. О методе эквивалентных замен расчета фрагментов сетей связи // Изв. АН
СССР. Техн. кибернетика. 1979. № 6. С. 92-99.
9. Neuts M.F. A versatile Markov point process // Journal of Applied Probability. 1979. V. 16. P. 764-779.
10. Lucantoni D.M. New results on the single server queue with a batch markovian arrival process // Communication in Statistics Stochastic Models. 1991. V. 7. P. 1-46.
11. Lucantoni D.M., Neuts, M.F. Some steady-state distributions for the MAP/SM/1 queue // Communication in Statistics Stochastic Models. 1994. V. 10. P. 575-598.
12. Лившиц К.И., Бублик Я.С. Вероятность разорения страховой компании при дважды стохастическом потоке страховых выплат // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2010. № 1(10). С. 66-77.
13. Yang Y.W. Woo Shin. BMAP/G/1 queue with correlated arrivals of customers and disasters // Operation Research Letters. 2004. V. 32, Issue 4. P. 364-373.
14. Artalejo Jesus R., Chakravarthy Srinivas R. Computational analysis of the maximal queue length in the MAP/M/c retrival queue // Applied Mathematics and Computation. 2006. V. 183, Issue 2. P. 1399-1409.
15. Best J. Doubly Stochastic Processes: an Approach for Understanding Central Nervous System Activity // Selected Topics on Applied Mathematics, Circuits, Systems and Signals. WSEAS Press, 2009. P. 155-158.
16. Bushalnov I. V., Gortsev A.M., Nezhel'skaya L.A. Estimating parameters of the synchronous twofold-stochastic flow of events // Automation and Remote Control. 2008. V. 69, No. 9. P. 1517-1533.
17. Горцев А.М., Нежельская Л.А. Оценивание длительности «мертвого времени» и интенсивностей синхронного дважды стохастического потока событий // Радиотехника. 2004. №. 10. С. 8-16.
18. Василевская Т.П., Горцев А.М., Нежельская Л.А. Оценивание длительности мертвого времени и параметров синхронного альтернирующего потока с проявлением либо непроявлением событий // Вестник Томского государственного университета. 2004. № 9. С. 129-138.
19. Горцев А.М., Нежельская Л.А. Оценивание длительности мертвого времени и параметров синхронного альтернирующего потока событий // Вестник Томского государственного университета. 2003. № 6. С. 232-239.
20. Горцев А.М., Нежельская Л.А. Оцениваение параметров синхронного дважды стохастического потока событий методом моментов // Вестник Томского государственного университета. 2002. № 1. С. 24-29.
21. Горцев А.М., Нежельская Л.А. Оценка параметров синхронного альтернирующего пуассоновского потока событий методом моментов // Радиотехника. 1995. №. 7-8. С. 6-10.
22. Леонова М.А., Нежельская Л.А. Оценка максимального правдоподобия длительности мертвого времени в обобщенном асинхронном потоке событий // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2013. № 2 (23). С. 54-63.
23. Горцев А.М., Леонова М.А., Нежельская Л.А. Совместная плотность вероятностей длительности интервалов обобщенного асинхронного потока событий при непродлевающемся мертвом времени // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2012. № 4 (21). С. 14-25.
24. Леонова М.А., Нежельская Л.А. Вероятность ошибки при оценивании состояний обобщенного синхронного потока событий // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2012. № 2(19). С. 88-101.
25. Gortsev A.M., Nezhel'skaya L.A. An asynchronous double stochastic flow with initiation of superfluous events // Discrete Mathematics and Applications. 2011. V. 21, No. 3. P. 283-290.
26. Горцев А.М., Нежельская Л.А. Оценивание параметров асинхронного потока с инициированием лишних событий методом моментов // Вестник Томского государственного университета. 2006. № 18. С. 267-273.
27. Горцев А.М., Нежельская Л.А., Шевченко Т.И. Оценивание состояний MC-потока событий при наличии ошибок измерений // Известия высших учебных заведений. Физика. 1993. № 12. С. 67-85.
28. Калягин А.А., Нежельская Л.А. Оценка длительности мертвого времени в обобщенном полушнхронном потоке событий // Новые информационные технологии в исследовании сложных структур : материалы Десятой российской конференции с международным участием (9-13 июня 2014 г.). Томск : Издательский Дом Томского государственного университета, 2014. С. 96-97.
29. Горцев А.М., Калягин А.А., Нежельская Л.А. Оптимальная оценка состояний обобщенного полусинхронного потока событий // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2010. № 2(11). С. 66-81.
30. Горцев А. М., Нежельская Л. А. Полусинхронный дважды стохастический поток событий при продлевающемся мертвом времени // Вычислительные технологии. 2008. Т. 13, № 1. С. 31-41.
31. Gortsev A.M., Nezhel'skaya L.A. Estimation of the dead-time period and parameters of a semi-synchronous double-stochasic stream of events // Measurement Techniques. 2003. V. 46, No. 6. P. 536-545.
32. Горцев А. М., Нежельская Л. А. Оценивание параметров полусинхронного дважды стохастического потока событий методом моментов // Вестник Томского государственного университета. 2002. № 1. С. 18-23.
33. Нежельская Л. А. Оптимальное оценивание состояний полусинхронного потока событий в условиях его частичной наблюдаемости // Вестник Томского государственного университета. 2000. № 269. С. 95-98.
34. Горцев А.М., Нежельская Л.А. О связи MC-потоков и MAP-потоков событий // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2011. № 1(14). С. 13-21.
35. Gortsev A.M., Nezhel'skaya L.A., Solov'ev A.A. Optimal State Estimation in MAP Event Flows with Unextendable Died Time // Automation and Remote Control. 2012. V. 73, No. 8. P. 1316-1326.
36. Горцев А.М., Шмырин И.С. Оптимальная оценка состояний дважды стохастического потока событий при наличии ошибок в измерениях моментов времени // Автоматика и телемеханика. 1999. № 1. С. 52-66.
37. Горцев А.М., Леонова М.А. Оптимальная оценка состояний обобщенного асинхронного дважды стохастического потока // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2010. № 1. С. 33-47.
38. Горцев А.М., Зуевич В.Л. Оптимальная оценка состояний асинхронного дважды стохастического потока событий с произвольным числом состояний // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2010. № 2. С. 44-65.
39. Апанасович В.В., Коляда А.А., Чернявский А.Ф. Статистический анализ случайных потоков в физическом эксперименте. Минск : Университетское, 1988. 254 с.
40. Горцев А.М., Васильева Л.А. Оценивание параметров дважды стохастического потока событий в условиях его неполной наблюдаемости // Автоматика и телемеханика. 2002. № 3. С. 179-184.
41. Горцев А.М., Завгородняя М.Е. Оценка параметров альтернирующего потока событий при условии его частичной наблюдаемости // Оптика атмосферы и океана. 1997. Т. 10, № 3. С. 273-280.
42. Горцев А.М., Климов И.С. Оценка интенсивности пуассоновского потока событий в условиях частичной его ненаблюдаемости // Радиотехника. 1991. № 12. С. 3-7.
43. Горцев А.М., Паршина М.Е. Оценивание параметров альтернирующего потока событий в условиях «мертвого времени» // Известия высших учебных заведений. Физика. 1999. № 4. С. 8-13.
44. Горцев А.М., Ниссенбаум О.В. Оценивание длительности мертвого времени и параметров асинхронного альтернирующего потока событий при непродлевающемся мертвом времени // Известия высших учебных заведений. Физика. 2005. № 10. С. 35-40.
45. Горцев А.М., Ниссенбаум О.В. Оценивание длительности мертвого времени и параметров асинхронного альтернирующего потока событий с инициированием лишних событий // Вестник Томского государственного университета. 2004. № 284. С. 137-145.
46. Горцев А. М., Климов И. С. Оценивание периода ненаблюдаемости и интенсивности пуассоновского потока событий // Радиотехника. 1996. № 2. С. 8-11.
47. Gortsev A.M., Nezhel'skaya L.A. Estimate of parameters of synchronously alternating Poisson stream of events by the moment method // Telecommunications and Radio Engineering (English translation of Elektrosvyaz and Radiotekhnika). 1996. V. 50, No. 1. P. 56-63.
48. Нежельская Л.А. Оптимальная оценка состояний синхронного MC-потока событий // Избранные доклады международной конференции «Всесибирские чтения по математике и механике». Томск : Изд-во ТГУ, 1997. Т. 1. С. 97-102.
49. Gortsev A.M., Shmyrin I.S. Optimal estimation of the states of a double stochastic flow of events in the presence oof measurement errors of time instants // Automation and Remote Control. 1999. V. 60, No. 1. P. 41-51.
50. Горцев А. М., Шмырин И. С. Оптимальная оценка параметров дважды стохастического пуассоновского потока событий при наличии ошибок в измерениях моментов наступления событий // Известия высших учебных заведений. Физика. 1999. № 4. С. 19.
51. Горцев А. М., Нежельская Л. А. Синхронный дважды стохастический поток событий при продлевающемся мертвом времени // Теория вероятностей, случайные процессы, математическая статистика и приложения: материалы международной конференции. Минск : Изд-во БГУ, 2005. С. 60-69.
52. Бушланов И.В., Горцев А.М. Оптимальная оценка состояний синхронного дважды стохастического потока событий // Автоматика и телемеханика. 2004. № 9. С. 40-51.
53. Горцев А. М., Калягин А. А., Нежельская Л. А. Оценка максимального правдоподобия длительности мертвого времени в обобщенном полусинхронном потоке // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2015. № 1. С. 27-37.
54. Горцев А.М., Соловьев А.А. Совместная плотность вероятностей длительности интервалов потока физических событий при непродлевающемся мертвом времени // Известия высших учебных заведений. Физика. 2014. Т. 57, № 7. С. 103-111.
55. Горцев А. М., Леонова М. А., Нежельская Л. А. Сравнение МП- и ММ-оценок длительности мертвого времени в обобщенном асинхронном потоке событий // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2013. № 4 (25). С. 32-42.
56. Калягин А. А., Нежельская Л. А. Сравнение МП- и ММ-оценок длительности мертвого времени в обобщенном полусинхронном потоке событий // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2015. № 3(32). С. 23-32.
57. Горцев А. М., Калягин А. А., Нежельская Л. А. Совместная плотность вероятностей длительности интервалов обобщенного полусинхронного потока событий при непродлевающемся мертвом времени // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2014. № 2(27). С. 19-29.
58. Горцев А. М., Соловьев А. А. Оценка максимального правдоподобия длительности непродлевающегося мертвого времени в потоке физических событий // Известия высших учебных заведений. Физика. 2015. Т. 58, № 11. С. 141-149.
59. Горцев А. М., Соловьев А. А. Сравнение МП- и ММ-оценок длительности непродлевающегося мертвого времени в MAP-потоке событий // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2015. № 4. С. 13-22.
Сиротина Мария Николаевна. E-mail: mashuliagol@mail.ru Томский государственный университет
Поступила в редакцию 14 мая 2016 г.
Sirotina Maria N. (Tomsk State University, Russian Federation).
Comparison of ML- and MM-estimators of period duration of dead time in modulated synchronous double stochastic flow of events.
Keywords: modulated synchronous flow; nonextendable dead time; ML-estimation; MM-estimation; dead time duration. DOI: 10.17223/19988605/36/7
One considers the modulated synchronous twice stochastic flow of events, which rate is piecewise constant random process X(t) with two states: X1, X2 (X1 > X2 > 0) . The time when process X(t) is staying in i-state has exponential probability distribution function with the parameter ai, i = 1,2 . During the time random interval when X(t) = Xi, there is the Poisson flow with the rate Xi, i = 1,2 .
A state transition of process X(t) occurs in the moment of the Poisson flow event arrival, moreover, the passing from the first to the
second state is realized with probability p, the passing from the second to the first state is realized with probability q.
The flow functions in conditions of dead time, which is one of the falsifying factors of state and parameter estimation of the flow. After each registered event, there is the time of fixed duration T (dead time), during which another flow events is inaccessible for observation. Now consider the nonextendable dead time, that is all events, which happen during dead time interval without its prolongation. At the end of the period duration of the dead time, the first happened event creates dead time period of duration T again and etc. Note that for the flow, which functions in dead time conditions, events are observable if they did not get into dead time interval (events which cause its approach). For the flow which functions in dead time absence all events are observable.
One of confounding factor during the flow state and parameter estimation is dead time of recording device, which is generated by observable flow event occurrence. All other events occurred during the dead time interval are not accessible for observation.
The main purpose of the research work is to obtain dead time period estimator T using the maximum likelihood method and method of moments (ML-estimators and MM-estimators respectively) and to compare these estimators.
During the investigation it was obtained that ML-estimator of dead time period duration takes a value T = %m = min %k, k = 1, n . It
means that the solution of optimization problem is the dead time estimate T = xm . Also, it was designed an algorithm of the single MM-estimator obtaining. On the base of these methods it was realized a computing program and it was made a statistical experiment during which for different values of input flow parameters X1,X2, a1, a2, p, q, T there was computed the values of ML-estimates
TMJj and MM-estimates Tj, j = 1, N (N is a number of experiment steps). For the values obtained there are calculated the values of
^ N ^ ^ N ^ _
sample variances Vun = (1/N)^(Tj - T)2, TV]UM = (1/N)^ OmM - T)2 , j = 1, N . Then it was chosen which of estimates is
j=i j=1
better: if Vun < Vuu then ML-estimate is better than MM-estimate, otherwise MM-estimate is better than ML-estimate.
REFERENCES
1. Sirotina, M.N. (2014) Optimal state estimation of modulated synchronous twice-stochastic flow of events in conditions of fixed dead
time. Vestnik Tomskogo gosudarstvennogo universiteta. Upravlenie, vychislitel'naya tekhnika i informatika - Tomsk State University Journal of Control and Computer Science. 1(26). pp. 63-74. (In Russian).
2. Sirotina, M.N. & Gortsev, A.M. (2015) [Probability density function of interval duration between two adjacent events of modulated synchronous flow in conditions of fixed dead time]. Proceedings International Scientific Conference devoted to the 80th anniversary of Prof. G.A. Medvedev. Belorussia, 23-26 February, 2015. Minsk: RIVSH. pp. 299-304.
3. Sirotina, M.N. & Gortsev, A.M. (2015) Joint probability density function of interval duration between two adjacent events of modulated synchronous flow in conditions of fixed dead time and recurrence flow conditions. Vestnik Tomskogo gosudarstvennogo uni-versiteta. Upravlenie, vychislitel'naya tekhnika i informatika - Tomsk State University Journal of Control and Computer Science. 2(31). pp. 53-67. (In Russian)
4. Sirotina, M.N. & Gortsev, A.M. (2016) Maxumum likelihood function estimation of dead time duration in modulated synchronous
doubly stochastic flow of events. Vestnik Tomskogo gosudarstvennogo universiteta. Upravlenie, vychislitel'naya tekhnika i informat-ika - Tomsk State University Journal of Control and Computer Science. 1(34). pp. 50-64. (In Russian)
5. Dudin, A.N. & Klimenok, V.N. (2000) Sistemy massovogo obsluzhivaniya s korrelirovannymi potokami [Queueing systems with correlated flows]. Minsk: BSU.
6. Kingman, J.F.C. (1964) On doubly stochastic Poisson process. Proceedings Cambridge Phylosophical Society. 60(4). pp. 923-930.
7. Basharin, G.P., Kokotushkin, V.A. & Naumov, V.A. (1979) O metode ekvivalentnykh zamen rascheta fragmentov setey svyazi
[About the method of renewal of subnetwork computation]. Izv. ANSSSR. Tekhn. kibernetika. 6. pp. 92-99.
8. Basharin, G.P., Kokotushkin, V.A. & Naumov, V.A. (1980) Basharin, G.P., Kokotushkin, V.A. & Naumov, V.A. (1980) O metode ekvivalentnykh zamen rascheta fragmentov setey svyazi [About the method of renewal of subnetwork computation]. Izv. AN SSSR. Tekhn. kibernetika. 1. pp. 55-61.
9. Neuts, M.F. (1979) A versatile Markov point process. Journal of Applied Probability. 16. pp. 764-779. DOI: 10.2307/3213143
10. Lucantoni, D.M. (1991) New results on the single server queue with a batch markovian arrival process. Communication in Statistics Stochastic Models. 7. pp. 1-46. DOI: 10.1080/15326349108807174
11. Lucantoni, D.M. & Neuts, M.F. (1994) Some steady-state distributions for the MAP/SM/1 queue. Communication in Statistics Stochastic Models. 10. pp. 575-598. DOI: 10.1080/15326349408807311
12. Livshits, K.N. & Bublik, Ya.S. (2010) The probability of insurance company collapse with twice stochastic flow of insurance payouts. Vestnik Tomskogo gosudarstvennogo universiteta. Upravlenie, vychislitel'naya tekhnika i informatika - Tomsk State University Journal of Control and Computer Science. 1(10). pp. 66-77. (In Russian).
13. Yang, Y.W. & Woo Shin. (2004) BMAP/G/1 queue with correlated arrivals of customers and disasters. Operation Research Letters. 32(4). pp. 364-373. DOI: 10.1016/j.orl.2003.09.005
14. Artalejo, J.R. & Chakravarthy, S.R. (2006) Computational analysis of the maximal queue length in the MAP/M/c retrival queue. Applied Mathematics and Computation. 183(2). pp. 1399-1409.
15. Best, J. (2009) Doubly Stochastic Processes: an Approach for Understanding Central Nervous System Activity. In: Pardalos, P. et al. (eds) Selected Topics on Applied Mathematics, Circuits, Systems and Signals. WSEAS Press. pp. 155-158.
16. Bushalnov, I.V., Gortsev, A.M. & Nezhelskaya, L.A. (2008) Estimating parameters of the synchronous twofold-stochastic flow of events. Automation and Remote Control. 69(9). pp. 1517-1533. DOI: 10.1134/S0005117908090075
17. Gortsev, A.M. & Nezhelskaya, L.A. (2004) Estimation of the dead time period and intensities of the synchronous double stochastic event flow. Radioengineering. 10. pp. 8-16.
18. Vasilevskaya, T.P., Gortsev, A.M. & Nezhelskaya, L.A. (2004) Otsenivanie dlitel'nosti mertvogo vremeni i parametrov sinkhron-nogo al'terniruyushchego potoka s proyavleniem libo neproyavleniem sobytiy [Dead time and parameters estimation of synchronous alternating flow with or without event manifestation]. Vestnik Tomskogo gosudarstvennogo universiteta - Tomsk State University Journal. 9. pp. 129-138.
19. Gortsev, A.M. & Nezhelskaya, L.A. (2003) Otsenivanie dlitel'nosti mertvogo vremeni i parametrov sinkhronnogo al'terniruyushchego potoka sobytiy [Dead time period and parameter estimation of synchronous alternating flow of events]. Vestnik Tomskogo gosudarstvennogo universiteta - Tomsk State University Journal. 6. pp. 232-239.
20. Gortsev, A.M. & Nezhelskaya, L.A. (2002) Otsenivaenie parametrov sinkhronnogo dvazhdy stokhasticheskogo potoka sobytiy metodom momentov [Parameter estimation of synchronous twice-stochastic flow of events using the method of moments]. Vestnik Tomskogo gosudarstvennogo universiteta - Tomsk State University Journal. 1. pp. 24-29.
21. Gortsev, A.M. & Nezhelskaya, L.A. (1995) Estimation of the parameters of a synchro-alternating Poisson event flow by the method of moments. Radiotechnics. 40(7-8). pp. 6-10.
22. Leonova, M.A. & Nezhelskaya, L.A. (2013) Maximum likelihood estimation of dead time value at a generalized asynchronous flow of events. Vestnik Tomskogo gosudarstvennogo universiteta. Upravlenie, vychislitel'naya tekhnika i informatika - Tomsk State University Journal of Control and Computer Science. 2(23). pp. 54-63. (In Russian).
23. Gortsev, A.M., Leonova, M.A. & Nezhelskaya, L.A. (2012) Joint probability density function of interval duration of generic asynchronous event flow in conditions of fixed dead time. Vestnik Tomskogo gosudarstvennogo universiteta. Upravlenie, vychislitel'naya tekhnika i informatika - Tomsk State University Journal of Control and Computer Science. 4(21). pp. 14-25. (In Russian).
24. Leonova, M.A. & Nezhelskaya, L.A. (2012) Error probability of state estimation of generic synchronous flow. Vestnik Tomskogo gosudarstvennogo universiteta. Upravlenie, vychislitel'naya tekhnika i informatika - Tomsk State University Journal of Control and Computer Science. 2(19). pp. 88-101. (In Russian).
25. Gortsev, A.M. & Nezhelskaya, L.A. (2011) An asynchronous double stochastic flow with initiation of superfluous events. Discrete Mathematics and Applications. 21(3). pp. 283-290. DOI: 10.4213/dm1141
26. Gorstev, A.M. & Nezhelskaya, L.A. (2006) Otsenivanie parametrov asinkhronnogo potoka s initsiirovaniem lishnikh sobytiy metodom momentov [Parameter estimation of asynchronous flow with additional event initiation using the method of moments]. Vestnik Tomskogo gosudarstvennogo universiteta. 18. pp. 267-273.
27. Gortsev, A.M., Nezhelskaya, L.A. & Shevchenko, T.I. (1993) Estimation of the states of an MC-stream of events in the presence of measurement errors. Russian Physics Journal. 36(12). pp. 1153-1167.
28. Kalyagin, A.A. & Nezhelskaya, L.A. (2014) [Dead time period estimation of generic semi-synchronous flow of events]. Novye in-formatsionnye tekhnologii v issledovanii slozhnykh struktur [New information technologies in complex structure investigation]. Proceedings of 10th Russian conference with international participation (9-13 of June, 2014). Tomsk: Tomsk State University. pp. 9697. (In Russian).
29. Gortsev, A.M., Kalyagin, A.A. & Nezhelskaya, L.A. (2010) Optimal state estimation of generic semi-synchronous flow of events. Vestnik Tomskogo gosudarstvennogo universiteta. Upravlenie, vychislitel'naya tekhnika i informatika - Tomsk State University Journal of Control and Computer Science. 2(11). pp. 66-81. (In Russian).
30. Gortsev, A.M. & Nezhelskaya, L.A. (2008) Semi-synchronous twice-stochastic event flow in conditions of prolonged dead time. Vychislitel'nye tekhnologii - Computational Technologies. 13(1). pp. 31-41. (In Russian).
31. Gortsev, A.M. & Nezhelskaya, L.A. (2003) Estimation of the dead-time period and parameters of a semi-synchronous double-stochasic stream of events. Measurement Techniques. 46(6). pp. 536-545. DOI: 10.1023/A:1025499509015
32. Gortsev A.M. & Nezhelskaya, L.A. (2002) Otsenivanie parametrov polusinkhronnogo dvazhdy stokhasticheskogo potoka sobytiy metodom momentov [Parameter estimation of semi-synchronous twice-stochastic event flow using the method of moments]. Vestnik Tomskogo gosudarstvennogo universiteta - Tomsk State University Journal. 1. pp. 18-23.
33. Nezhelskaya, L.A. (2000) Optimal'noe otsenivanie sostoyaniy polusinkhronnogo potoka sobytiy v usloviyakh ego chastichnoy na-blyudaemosti [Optimal state estimation of semi-synchronus flow in conditions of its incomplete observability]. Vestnik Tomskogo gosudarstvennogo universiteta - Tomsk State University Journal. 269. pp. 95-98.
34. Gortsev, A.M. & Nezhelskaya, L.A. (2011) About connectivity of MC-flows and MAP-flows of events. Vestnik Tomskogo gosudarstvennogo universiteta. Upravlenie, vychislitel'naya tekhnika i informatika - Tomsk State University Journal of Control and Computer Science. 1(14). pp. 13-21. (In Russian).
35. Gortsev, A.M., Nezhelskaya, L.A. & Solovyev, A.A. (2012) Optimal State Estimation in MAP Event Flows with Unextendable Dead Time. Automation and Remote Control. 73(8). pp. 1316-1326. DOI: 10.1134/S000511791208005X
36. Gortsev, A.M. & Shmyrin, I.S. (1999) Optimal state estimation of doubly stochastic flow of events under conditions of time moments measurement errors. Automation and Remote Control. 1. pp. 52-66.
37. Gortsev, A.M. & Leonova, M.A. (2010) Optimal state estimation of generic asynchronous doubly stochastic flow. Vestnik Tomskogo gosudarstvennogo universiteta. Upravlenie, vychislitel'naya tekhnika i informatika - Tomsk State University Journal of Control and Computer Science. 1(10). pp. 33-47. (In Russian).
38. Gortsev, A.M. & Zuevich, V.L. (2010) Optimal state estimation of generic asynchronous doubly stochastic flow with random number of states. Vestnik Tomskogo gosudarstvennogo universiteta. Upravlenie, vychislitel'naya tekhnika i informatika - Tomsk State University Journal of Control and Computer Science. 2(11). pp. 44-65. (In Russian).
39. Apanasovich, V.V., Kolyada, A.A. & Chernyavskiy, A.F. (1988) Statisticheskiy analiz sluchaynykh potokov v fizicheskom eksperi-mente [Statistical analysis of stochastic flows in physical experiment]. Minsk: Universitetskoe.
40. Gortsev, A.M. & Vasilieva, L.A. (2002) Parameter estimation of a doubly stochastic flow of events under incomplete observability. Avtomatika i Telemekhanika. 3. pp. 179-184.
41. Gortsev, A.M. & Zavgorodnaya, M.E. (1997) Otsenka parametrov al'terniruyushchego potoka sobytiy pri uslovii ego chastichnoy nablyudaemosti [Parameter estimation of alternating flow of events under conditions of particulate observabiity]. Optika atmosfery i okeana - Atmospheric and Oceanic Optics. 10(3). pp. 273-280.
42. Gortsev, A.M. & Klimov, I.S. (1991) An estimate for intensity of Poisson flow of events under the condition of its partial missing. Radiotekhnika. 12. pp. 3-7.
43. Gortsev, A.M. & Parshina, M.E. (1999) Otsenivanie parametrov al'terniruyushchego potoka sobytiy v usloviyakh "mertvogo vremeni" [Parameter estimation of alternating flow of events under conditions of dead time]. Russian Physics Journal. 4. pp. 8-13.
44. Gortsev, A.M. & Nissenbaum, O.V. (2005) Dead time and parameter estimation of asynchronous alternating flow in conditions of a fixed dead time. Higher Educational Establishment Journal. Physics. 10. pp. 35-40.
45. Gortsev, A.M. & Nissenbaum, O.V. (2004) Dead time and parameter estimation of asynchronous alternating flow with additional event initiation. Vestnik Tomskogo gosudarstvennogo universiteta - Tomsk State University Journal. 284. pp. 137-145. (In Russian).
46. Gortsev, A.M. & Klimov, I.S. (1996) Estimation of the non-observability period and intensity of Poisson event flow. Radiotekhnika. 2. pp. 8-11.
47. Gortsev, A.M. & Nezhelskaya, L.A. (1996) Estimate of parameters of synchronously alternating Poisson stream of events by the moment method. Telecommunications and Radio Engineering. 50(1). pp. 56-63.
48. Nezhelskaya, L.A. (1997) [Optimal state estimation of synchronous MC-flow of events]. Vsesibirskie chteniyapo matematike i mek-hanike [Siberian readings on mathematics and mechanics]. Proc. of the International Conference. Vol. 1. Tomsk: Tomsk State University. pp. 97-102. (In Russian).
49. Gortsev, A.M. & Shmyrin, I.S. (1999) Optimal estimation of states of a double stochastic flow of events under conditions of time moments measurement errors. Automation and Remote Control. 60(1). pp. 41-51.
50. Gortsev, A.M. & Shmyrin, I.S. (1999) Optimal parameter estimation of a double stochastic Poisson flow of events in presence of measurement errors of time instants. Higher Educational Establishment Journal. Physics. 4. p. 19. (In Russian).
51. Gortsev, A.M. & Nezhelskaya, L.A. (2005) [Synchronous twice-stochastic flow of events in conditions of fixed dead time]. Teoriya veroyatnostey, sluchaynye protsessy, matematicheskaya statistika i prilozheniya [Probability theory, stochastic processes, mathemat-ic statistics and applications]. Proc. of the International Conference. Minsk: BSU. pp. 60-69. (In Russian).
52. Bushlanov, I.V. & Gortsev, A.M. (2004) Optimal estimation of the states of a synchronous double stochastic flow of events. Automation and Remote Control. 65(9). pp. 1389-1399. DOI: 10.1023/B:AURC.0000041418.09187.63
53. Gortsev, A.M., Kalyagin, A.A. & Nezhelskaya, L.A. (2015) Maximum likelihood estimation of dead time of generic semichronous flow. Vestnik Tomskogo gosudarstvennogo universiteta. Upravlenie, vychislitel'naya tekhnika i informatika - Tomsk State University Journal of Control and Computer Science. 1(30). pp. 27-37. (In Russian).
54. Gortsev, A.M. & Solovyev, A.A. (2014) Joint Probability Density of Interarrival Interval of a Flow of Physical Events with Unex-tendable Dead Time Period. Russian Physics Journal. 57(7). pp. 973-983. DOI: 10.1007/s11182-014-0333-4
55. Gortsev, A.M., Leonova, M.A. & Nezhelskaya, L.A. (2013) Comparison of MP- and MM-estimations of dead time in generic asynchronous flow. Vestnik Tomskogo gosudarstvennogo universiteta. Upravlenie, vychislitel'naya tekhnika i informatika - Tomsk State University Journal of Control and Computer Science. 4(25). pp. 32-42. (In Russian).
56. Kalyagin, A.A. & Nezhelskaya, L.A. (2015) Comparison of MP- and MM-estimations of dead time in generic semisynchronous flow of events. Vestnik Tomskogo gosudarstvennogo universiteta. Upravlenie, vychislitel'naya tekhnika i informatika - Tomsk State University Journal of Control and Computer Science. 3(32). pp. 23-32. (In Russian).
57. Gortsev, A.M., Kalyagin, A.A. & Nezhelskaya, L.A. (2014) The joint probability density of duration of the intervals in a generalized semi-synchronous flow of events with unprolonging dead time. Vestnik Tomskogo gosudarstvennogo universiteta. Upravlenie, vychislitel'naya tekhnika i informatika - Tomsk State University Journal of Control and Computer Science. 2(27). pp. 19-29. (In Russian).
58. Gortsev, A.M. & Solovyev, A.A. (2015) Maximum likelihood estimation of fixed dead time in physical flow of events. Russian Physics Journal. 58(11). pp. 141-149.
59. Gortsev, A.M. & Solovyev, A.A. (2015) Comparison of MP- and MM-estimations of fixed dead time in MAP-flow of events. Vestnik Tomskogo gosudarstvennogo universiteta. Upravlenie, vychislitel'naya tekhnika i informatika - Tomsk State University Journal of Control and Computer Science. 4(33). pp. 13-22. (In Russian).
