CC BY
13
ВЕСТНИК ТОМСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА
2016 Управление, вычислительная техника и информатика № 2 (35)
ОБРАБОТКА ИНФОРМАЦИИ
УДК 519.21
DOI: 10.17223/19988605/35/2
М.А. Бахолдина, А.М. Горцев
ОЦЕНКА МАКСИМАЛЬНОГО ПРАВДОПОДОБИЯ ДЛИТЕЛЬНОСТИ НЕПРОДЛЕВАЮЩЕГОСЯ МЕРТВОГО ВРЕМЕНИ В МОДУЛИРОВАННОМ ОБОБЩЕННОМ ПОЛУСИНХРОННОМ ПОТОКЕ СОБЫТИЙ
Рассматривается модулированный обобщенный полусинхронный поток событий, являющийся одной из математических моделей информационных потоков заявок, функционирующих в телекоммуникационных и информационно-вычислительных сетях связи, и относящийся к классу дважды стохастических потоков событий (DSPPs). Функционирование потока рассматривается в условиях непродлевающегося мертвого времени. В статье приводятся аналитические результаты по нахождению оценки максимального правдоподобия длительности мертвого времени по наблюдениям за моментами наступления событий в потоке.
Ключевые слова: модулированный обобщенный полусинхронный поток событий; дважды стохастический поток событий (DSPP); MAP (Markovian Arrival Process)-поток событий; непродлевающееся мертвое время; функция правдоподобия; оценка максимального правдоподобия; длительность мертвого времени.
Условия функционирования реальных систем массового обслуживания таковы, что если в отношении параметров обслуживающих устройств, как правило, можно утверждать, что они известны и с течением времени не меняются, то в отношении интенсивностей входящих потоков этого сказать во многих случаях нельзя. Более того, интенсивности входящих потоков заявок обычно меняются со временем; часто эти изменения носят случайный характер, что приводит к рассмотрению математических моделей дважды стохастических потоков событий (DSPPs) [1-9]. Интерес к рассмотрению дважды стохастических потоков событий проявляется неслучайно. Все это находит широкое применение в различных отраслях науки и техники, таких как теория сетей, P2P^ra и адаптивное вещание видео, системы оптической связи, статистическое моделирование, финансовая математика и др. [10-16]. Как было отмечено выше, в реальных ситуациях параметры, задающие входящий поток событий, известны либо частично, либо вообще неизвестны, либо, что еще более ухудшает ситуацию, изменяются со временем случайным образом. Поэтому при реализации адаптивного управления системой массового обслуживания возникают, в частности, следующие задачи: 1) задача фильтрации интенсивности потока (или задача оценивания состояний потока по наблюдениям за моментами наступления событий) [17-26]; 2) задача оценивания параметров потока по наблюдениям за моментами наступления событий [27-33].
Одним из искажающих факторов при оценке состояний и параметров потока выступает мертвое время регистрирующих приборов. Необходимость рассмотрения случая мертвого времени вызвана тем, что на практике любое регистрирующее устройство затрачивает на измерение и регистрацию события некоторое конечное время, в течение которого оно не способно правильно обработать следующее событие, т.е. событие, поступившее на обслуживающий прибор, порождает период так называемого мертвого времени [34-42], в течение которого другие наступившие события потока недоступны наблюдению (теряются). Можно считать, что этот период продолжается некоторое фиксированное время (непродле-вающееся мертвое время). В частности, подобные ситуации встречаются в компьютерных сетях, например, при использовании протокола случайного множественного доступа с обнаружением конфликта (протокол CSMA/CD). В момент регистрации (обнаружения) конфликта на входе некоторого узла сети по сети рассылается сигнал «заглушки»; в течение времени рассылки сигнала «заглушки» заявки, поступившие в данный узел сети, получают отказ в обслуживании и направляются в источник повторных
вызовов. Здесь время, в течение которого узел сети закрыт для обслуживания заявок, поступающих в него после обнаружения конфликта, можно трактовать как мертвое время прибора, регистрирующего конфликт в узле сети.
В настоящей работе рассматривается модулированный обобщенный полусинхронный поток событий, являющийся обобщением полусинхронного потока [43-46] и обобщенного полусинхронного потока событий [47-49] и относящийся к классу дважды стохастических потоков событий с кусочно-постоянной интенсивностью. Достаточно обширная литература по исследованию подобных потоков событий (асинхронных, синхронных и полусинхронных) приведена в [21, 22, 50, 51], при этом в [50] показано, что данные потоки могут быть представлены в виде моделей МАР-потоков событий. Настоящая статья является непосредственным развитием работ [51-55], где решается задача нахождения совместной плотности вероятности длительности интервалов модулированного обобщенного полусинхронного потока событий при непродлевающемся мертвом времени. В данном исследовании приводятся аналитические результаты по нахождению оценки максимального правдоподобия длительности мертвого времени в модулированном обобщенном полусинхронном потоке событий по наблюдениям за моментами наступления событий в потоке.
1. Постановка задачи
Рассматривается модулированный обобщенный полусинхронный поток событий (далее - поток или поток событий), интенсивность которого является кусочно-постоянным стационарным случайным процессом ) с двумя состояниями Х1 и Х2 > Х2 > 0). Будем говорить, что имеет место первое состояние процесса (потока), если Х(?) = Х1, и второе состояние процесса (потока), если Х(1) = Х 2. В течение временного интервала случайной длительности, когда процесс Х^) находится в состоянии Xi (Х(?) = Хг-), имеет место пуассоновский поток событий с интенсивностью X^, I = 1,2. Переход из первого состояния процесса Х^) во второе возможен в момент наступления события пуассоновского потока интенсивности Х1, при этом переход осуществляется с вероятностью р (0 < р < 1); с вероятностью 1 - р процесс Х^) остается в первом состоянии. Переход из первого состояния процесса Х(() во второе также возможен в произвольный момент времени, не совпадающий с моментом наступления события, при этом длительность пребывания процесса Х^) в первом состоянии распределена по экспоненциальному закону с параметром р: ^(т) = 1 - е-рт, т > 0 . Тогда длительность пребывания процесса Х^) в первом состоянии есть случайная величина с экспоненциальной функцией распределения ^](т) = 1 - е-(рХ+р)т, т > 0 . Переход из второго состояния процесса Х(^) в первое в момент наступления события пуассоновского потока интенсивности Х 2 невозможен и может осуществляться только в произвольный момент времени. При этом длительность пребывания процесса Х(^) во втором состоянии распределена по экспоненциальному закону с параметром а : ^2(т) = 1 - е-ат, т > 0 . В момент окончания второго состояния процесса Х(^) при его переходе из второго состояния в первое инициируется с вероятностью 5 (0 < 5 < 1) дополнительное событие. Отметим, что события пуассоновских потоков и дополнительные события неразличимы для наблюдателя. В сделанных предпосылках Х(^) - скрытый марковский процесс. При этом матрицы инфинитезимальных характеристик принимают вид
^0 =
-(*■! +Р) Р
а(1 - 5) - (Х2 + а)
А =
(1 - р)Х1 РХ1 а5 Х 2
Элементами матрицы О^ являются интенсивности переходов процесса Х(^) из состояния в состояние с наступлением события. Недиагональные элементы матрицы О0 - это интенсивности переходов из состояния в состояние без наступления события. Диагональные элементы матрицы О0 - это интенсивности выхода процесса Х(^) из своих состояний, взятые с противоположным знаком. Отметим, что если р = 0 , то имеет место обобщенный полусинхронный поток событий [47-49].
После каждого зарегистрированного в момент времени ¿к события наступает период мертвого времени фиксированной длительности Т , в течение которого другие события потока недоступны наблюдению. По окончании периода мертвого времени первое наступившее событие снова создает период мертвого времени длительности Т и т.д. (непродлевающееся мертвое время). Вариант возникающей ситуации приведен на рис.1, где 1, 2 - состояния процесса ); дополнительные события, которые могут наступать при переходе процесса ) из второго состояния в первое, помечены буквами 5; периоды мертвого времени длительности Т помечены штриховкой; ненаблюдаемые события отображены черными кружками, наблюдаемые 11, ¿2,... - белыми.
Рассматривается установившийся (стационарный) режим функционирования потока событий, поэтому переходными процессами на полуинтервале наблюдения (¿0,1], где 10 - начало наблюдений, I -окончание наблюдений (момент вынесения решения), пренебрегаем. Тогда без потери общности можно положить t0 = 0 . Поскольку процесс ) является принципиально ненаблюдаемым, то говорить о состоянии потока можно только в вероятностном смысле. Вся доступная информация о потоке - это моменты наступления событий ¿1,¿2,...,^ с начала наблюдения t0 до момента t.
Рис. 1. Формирование наблюдаемого потока событий
Параметры потока Л,1 > А, 2 > 0, 0 < р < 1, р> 0, а> 0, 0 < 5 < 1 полагаются известными, длительность мертвого времени Т неизвестна. Необходимо в момент окончания наблюдений (в момент времени t) на основании выборки наблюденных моментов наступления событий потока осуществить методом максимального правдоподобия оценку Т длительности мертвого времени.
2. Построение функции правдоподобия
Обозначим через тк = ¿к+1 - ¿к , к = 1,2,..., значение длительности к-го интервала между соседними событиями наблюдаемого потока (тк > 0 ). Так как рассматривается стационарный режим, то плотность вероятности значений длительности к-го интервала есть Рт (тк) = Рт (т), т > 0 , для любого к (индекс Т подчеркивает, что плотность вероятности зависит от длительности мертвого времени). В силу этого момент времени ¿к без потери общности можно положить равным нулю, т.е. момент наступления события есть т = 0 . Тогда одномерная плотность вероятности рт(т), т > 0 , примет вид [51, 55]:
Pt (т) =<
0, 0 <т<T,
Pi +P:
-f (T)
g-X_T ) . ^
Pi +P:
■f (T )
g- ), T> T, (i)
-(Pi + P2T
f (T) = X1a + (pX1 + P)(X2 +a5) + a(X1 -X2 -a5){(pX1 + Р)[(1 -p + p5)-X2 + 5p]-p^a}
F(T) = zz -qe-(P1+P2)T, P1 = pX1 +P, P2 =a , q = X1[X2 -p(X2 + a5)],
F (T)
z12 = 2(x1 +X2 +a+P + yj(X1 -X2 -a + P)2 + 4aP(1 -5)), 0 <z1 <z2.
Пусть т = t2 - , X2 = % - ¿2, . * *, т^ = tk+1 - , Т > 0, X2 > 0, ..., тк > 0, - последовательность измеренных в результате наблюдения за потоком в течение интервала наблюдения (0, t] значений длительностей интервалов между соседними событиями потока. Упорядочим величины т1, ..., тк по возрастанию: тт1п = т(1) < т(2) <... < т(к). В силу предпосылок последовательность моментов наступления событий t1, t2,..., tк,... образует вложенную цепь Маркова {X(tk)}, т.е. поток обладает марковским свойством, если его эволюцию рассматривать с момента наступления события tk, к = 1,2,.... Тогда функция правдоподобия с учетом (1) запишется в виде
L(T|t(1),..., т(к)) = Прт (т( 1)), 0 < T <Tmin; L( T | т(1),..., т(к)) = 0, T >ТтШ.
j=1
Поскольку задача заключается в построении оценки T длительности мертвого времени в предположении, что все параметры потока X1 >X 2 > 0, 0 < p < 1, P> 0, a> 0, 0 <5< 1 известны, то согласно методу максимального правдоподобия ее реализация есть решение оптимизационной задачи
L(T | т(1),...,т(к)) = ППрт(т( 1)) ^ max, 0 < T < Tmm, ТтШ > 0, (2)
1=1
где pT(т( 1)) определена в (1) для т = т1. Значение T , при котором функция правдоподобия (2) достигает своего глобального максимума, есть оценка T длительности мертвого времени.
3. Решение оптимизационной задачи (2)
Произведем переобозначение: тт = тт1п. В силу того что функция правдоподобия (2) отличается от нуля при 0 < T < тт , положим pT (т(j)) = 0, j = 2,к , при T >тт (тт > 0). В дальнейшем изложенная ситуация, когда принимается тт = 0, означает доопределение изучаемых функций в граничной точке. Перейдем к исследованию pT (тт) как функции переменной T (0 < T < тт). Отметим, что pT (тт) > 0, так как pT (тт) есть одномерная плотность вероятности. Исследуем производную p'T (тт) по T функции pT (тт). Имеем
PT(Тт) = 7-' P ,{z,e-Z1T- [Z1Z2 (P1 +P2)-zf(T) - f '(T)]e-Z1T -Z2e-Z2Tm [Z1Z2( +P2)-
(Z2 - Z1 ))P1 +P2 Г (3)
2f (T) - f '(T)] e -Z2T },
f '(T) = -a(X1 - X2 - a5)(p1 + P2 fe {(^^ + P)[ (1 - p + p5) - X2 + 5p]- p^a}} , 0 < T < т„
"СРх + Р2)Г
^Чг)
где /(Т), F(Т) определены в (1); /'(Т) - производная функции /(Т) по Т .
Лемма 1. Производная р'Т (хт) является положительной функцией переменной тт при Т = 0,
о <тт С рТ=оСтт) > 0).
1
1
z
Z2 Z1
Z2 Z1
Доказательство. Так как тт - любое неотрицательное число (тт > 0), то рТ=0(^т) можно рассматривать как функцию переменной тт . Подставляя Т = 0 в (3) и проделывая необходимые преобразования, получаем
С
) =
А\z2 - z1)
[(C - z A)e-ZA - Zi (C - Z2 A)e-z'x" ], xm > ü,
(4)
C = Я1а[А,1(1 - p + p5) + 5ß]+(X2 + a5)[X2(pX1 +ß) + pX^] = -ZjZ2(ß +ß2)+(Z1 + z2)А > ü ,
А = А^а + (+ Р)(Х2 + а5) = 2Х22 - д > 0 . Рассмотрим на предмет существования корней уравнение р'т=0(тт) = 0, которое с учетом (4) пре-
образуется к виду
Из (4) находим
— Z2 - ZlK
B = e
B =
Z1( Z2 A - C) Z2(Z1A - C)
P;=ü(^m = ü) = (C/ A)2 > ü, lim p'T==ü(xm) = ±0 .
(5)
(6)
Подставляя в (5) выражения для Z1, Z2, определенные в (1), получаем
B = -
■2A + ( +ß2)(Z1 + Z2) + ( +ß2-X2 -a + ß)2 + 4aß(1 -5)
-2A + (ß1 +ß2)(z1 + Z2)-(ß1 +ß2)(X1 -X2 - a+ ß)2 + 4aß(1 -5)
Bl
B2
Тогда имеем
B = Z12
-2A + (ß1 +ß2)( + X2 + a + ß)+(ß1 +ß2)(X1 - X2 - a + ß)2 + 4aß(1 - 5) > Z12[- 2A +
+ (ß1 +ß2)X1 +X2 +a + ß)+(ß1 +ß2)(X1 -X2 -a + ß)2 > 2z12[pX1a(1 -5) + aß(2-5)]> ü вне зависимости от знака выражения X1 -X 2 -a+ß (X1 -X 2 -a + ß> ü, X1 -X 2 -a+ß = ü, X1 - X2 -a + ß < ü ). Для выражения (5) имеем два варианта:
1) B1 = Z1(z2A-C) > ü; B2 = Z2(Z1A -C) > ü . Тогда B > ü и разность B1 -B2 = C(Z2 - Z1) > 0 . Следовательно, B1 > B2 и B > 1. Тогда уравнение (5) решения не имеет, следовательно, pT=ü(^m) > ü, xm > ü, так как в силу (6) pT=ü(^m = ü) > ü, при этом lim pT=ü(^m) = +ü ;
Xm ^да
2) B1 = Z1(Z2A - C) > ü; B2 = Z2(Z1A - C) < ü . Тогда B < ü и уравнение (5) решения не имеет, следовательно, pT=ü(^m) > ü, xm > ü, аналогично варианту 1.
Осталось рассмотреть особый случай Z1A - C = ü . Подставляя Z1A - C = ü в (4), получаем
pT=ü(^m ) =■
C
(Z2A - C)Z1e-Z1X" > ü, xm > ü,
-А2(^ -
так как В1= г1 (22А - С) > 0 всегда и, следовательно, 22А - С > 0 всегда.
Таким образом, если 22 (11А - С) > 0 либо 22 (г1А - С) < 0 , то рТ=0 (тт) > 0, хт > 0. Лемма 1 доказана. Замечание 1. д > д1 всегда.
Доказательство. В лемме 1 показано, что г2А - С > 0 всегда. Можно показать, что 12А - С = 11(д - д1), следовательно, д > д1 всегда.
Лемма 2. Производная р'Т (тт) строго больше нуля при Т = тт, 0 < хт <<х> (р'т=г (хт) > 0). Доказательство. Подставляя Т = тт в (3), получаем
p'O^m) = pt=X, (Tm) =
1
ß1 +ß2
!C + )"
-((+ß2
F2(Tm)
Y(Tm ) = C + (Z1 + Z2 )?[
- e
"(ß1 + ß2
(7)
т... -^да
2
Z
2
Z
2
ф(?) = (01 +Р2 ) - Л2 = -?2 +[2^2 —(01 +Р2 + ^2 )] - ^1^2 {21 -Р1 -Р2 )2 -Р1 -Р2 ) ,
где д, F(хт) определены в (1), А, С - в (4). Можно показать, что ф(д) =-(? - )(? - ?2), где
= 22 (21 - 01 - 02), д2 = 21 (22 - 01 - 02) - корни уравнения ф(?) = 0 . Отметим, что
С С
р'СТт = 0) = (С / А)2 > 0, НтрК) =——> 0, р'СТт = 0) - Шп рК) = . ( 0 ) фС?). (8)
Хт ^ 01 +02 Хт ^ А (01 +02 )
В силу того что С > 0, у(хт) > 0 для хт > 0, то знак производной р'(хт) определяется знаком ф(?). Тогда если ф(д) > 0 , то производная (7) строго больше нуля ( р'(Хт ) > 0) для Хт > 0.
Пусть ф(?) < 0 . Введем в рассмотрение вторую производную р''(хт) по переменной Т в точке Т = хт . Используя (3), находим
р''(Хт ) = -2122 ^ф?- *-(01+02 [[122 (21 + 22 -01 -02 )-(( + 22 +01 +02 )-(0+02 К ], Хт > 0,
Р (Хт)
где д , Р(хт) определены в (1), ф(д) определена в (7). По условию имеем
ф(?) *-(0! +02
- 2122 Р 3(Хт ) *
>0,
следовательно, знак производной р '' (хт) определяется знаком функции
у(Хт ) = 2122 (21 + 22 - 01 -02 ) -^ + 22 + 01 + 02 ^(0^ К , Хт > 0 .
Исследуем функцию у(хт) как функцию переменной хт (хт > 0). Имеем
У(0) = С - ?(01 +02 ), 11т У (Хт ) = = 2122 (21 + 22 - 01 -02 )= 2122 ^(1 - р) + Х 2 ]> 0 , где д определена в (1), С - в (4). Здесь возможны следующие варианты:
1) д < 0 и, следовательно, у(0) = С - ?(01 +02 )> 0. Тогда функция у(хт) - убывающая функция переменной хт (хт > 0); убывает от у(0) = С - ?(01 +02 )> 0 до у(да) = 2122 (21 + 22 - 01 - 02 )> 0. Следовательно, у(хт) > 0 , хт > 0 . Отсюда следует, что р''(хт) > 0, хт > 0, и функция р'(хт) является возрастающей функцией переменной хт (хт > 0). Таким образом, в силу (8) р' (хт) > 0 , хт > 0;
2) д = 0. Тогда у(хт) = 2122 (21 + 22 - 01 - 02)> 0, хт > 0. Результат идентичен результату предыдущего пункта;
3) д > 0. Тогда у(хт) - возрастающая функция переменной хт (хт > 0); возрастает от у(0) = С - +02) до у(да) = 2122 (21 + 22 - 01 -02 )> 0. В данном случае имеют место три варианта: а) У(0) > 0 ; б) у(0) = 0 ; в) у(0) < 0. Рассмотрим более подробно каждый из вариантов.
Вариант а.
Пусть у(0) = С - д( +02) > 0. Тогда у(хт) > 0, хт > 0. Отсюда следует, что р' (хт) > 0, хт > 0, и функция р' (хт) является возрастающей функцией переменной хт (хт > 0). Таким образом, в силу (8) р' (хт) > 0 , хт > 0 .
Вариант б. Пусть у(0) = С - +02 ) = 0. Тогда у(хт) > 0, хт > 0, причем равенство нулю (у(хт) = 0 ) достигается в единственной точке хт = 0. Результат идентичен результату варианта а.
Вариант в. Пусть у(0) = С - д( +02 )< 0. Тогда: 1) у(Хт) < 0, 0 < Хт <хт ; 2) у(Хт) = 0, Хт =х*п ; 3) У(хп) > 0 , хП < хп < да . Следовательно, в точке хп = хП производная р' (хп) достигает минимального значения р' (хП ). Точка минимума хП находится из уравнения р' (хт) = 0 :
Х* =- 1 1п 2122 (21 + 22 - 01 -02 ) 2122 (21 + 22 - 01 -02 ) > 0 т 01 +02 (21 + 22 +01 +02' (21 + 22 + 01 +02 Вычислим значение р' (хт) в точке хП :
р'(хП) = \с + ф(д) (21 +(2-01 -02)2 1 .
РКт) 01 +02 I 4(01 +02)д I
Тогда имеем
p'(<) > p-1ß2 |c + 9(q) (Z1 + ^ "ß 2 )2 } = 4C (ß1 + ß 2) 2 + ( +ß 2 )[c (2(1- p)2 +X 22 + 2pX^)-
-q2( +ß2))> 4C(ß1 + ß ){3C2 + ( +ß2)[(X12a(1 -p) + X22ß + X1X2p(X2 -p)2 +X22)-
- (X1X2 (1 - p))2 ( + ß2)]} = X14a(1 - p)3 + X1X2pa X1X2 (1 - p) + X12(1 - p)2 + X22 ]+ X24p1 > ü.
Отсюда следует, что p'(xm) > ü , xm > ü . Таким образом, если 9(q) < ü , то производная p' (xm) строго больше нуля (p' (xm) > ü) для xm > ü . Лемма 2 доказана.
Перейдем к исследованию производной p'T (xm) как функции переменной T , ü < T < xm . Рассмотрим (на предмет существования корней) уравнение p'T (xm) = ü , которое с учетом (3) приводится к виду
e-(Z2-Z1 )(x"-T '= f1(T)/ f2(T), ü < T <xm , (9)
f1(T) = Z1[Z1Z2 (ß1 +ß2 )-Zf(T) - f '(T)], ] = Z2 [Z1Z2 (ß1 +ß2 )-Z2f(T) - f'(T)] .
-(z -z )(x -t)
При этом функция ^(T) = e-Z2-Z1Ax"- ' есть возрастающая функция переменной T , ü < T < xm ; возрастает от x(ü) = e-(Z2-Z1 )x" до x(T = xm) = 1, ее максимальное значение есть 1.
Так как xm, в принципе, может равняться бесконечности, то изучим функции f1(T), f2(T) как функции переменной T , T > ü .
Утверждение 1. Для функции f1(T) = z1[z1z2(ß1 +ß2)-Z1 f (T) - f '(T)], T > ü, справедливо:
1) f1(ü) = ( +ß2))£(A-C)=( +ß2^(q-q)> ü, q1 = Z2(Z1 -ß1 -ß2), q > q1;
2) f (да) = lim f (T) = Z1 (Z2 A - C) = z2 (q - q) > ü;
t ^да
f'(T) = 9(q)e-(^+ß2)T[Z1Z2(Z1 -ß1 -ß2)-( +ß1 +ß2)qe-(ß)], T> ü ;
3)
' ' р 3(Т)
4) точка Т1 (д) - точка экстремума (либо точка перегиба) функции / (Т):
T1(q)=пг+г in Z( +ЛТ2}, ^^^ > ü; (ю)
ß1 +ß2 (Z1 +ß1 +ß2 )q q
5) f1(T1(q))=)2
4( +ß2 )q
12
" Z1 (Z1 -ß1 -ß2 )2 (Z2 -ß1 -ß2 )] Z1 -ß1 -ß2
(1 +ß1 +ß2 )2 J' q
r-Z2 (z1 -ß1 ß2 )
> ü,
где Я11 = Я1 = ^( -Р1 -Р2), Я12 = ^ (( Р| Р-) Р Р2 ) - корни уравнения /ДТДд)) = 0 . Здесь д,
(( +Р1 +Р2)
^(Т) определены в (1); А, С - в (4); ф(д) - в (7).
Лемма 3. Функция /1(Т) является положительной функцией (/1(Т) > 0) переменной Т , Т > 0 . Доказательство. Пусть д = д1, тогда /1(Т, д = д1) = г12[г2(Р1 +Р2)-г1г2 + г2( -Р1 -Р2)] = 0 ; пусть
я = я2, тогда /1 (Т, я = Я2) = (2 - г1 )((31 + Р2) > 0 . Точки д = Я1, Я = Я2 из дальнейшего рассмотрения исключаем, так как в этих точках ф(д) = 0 . Рассмотрим далее возможные случаи.
Случай 1. ф(я) > 0, 0 <я1 <Я <Я2. Знак производной /1(Т) функции /1 (Т) определяется знаком функции у1(Т) = 2Х22 ( -Р1 -Р2 )-(( +Р1 +Р2)-(Р1+Р2)Т, Т > 0 . В данном случае функция у1(Т) является возрастающей функцией переменной Т, Т > 0; возрастает от у1(Т = 0, д) < 0 до у1(да) = 2121 ( -Р1 -Р2 )> 0, т.е. один раз пересекает ось абсцисс в точке Т1(д), определенной выражени-
ем (10). Таким образом, в точке T = T1(q) реализуется единственный минимум функции f (T). Можно показать, что f1(T1(q)) > 0 , тогда f1(T) > 0 , T > 0 . Итак, если q<q) > 0, 0 < q1 < q < q2, то f1(T) > 0 , T > 0 .
Случай 2. 9(q) > 0, q1 = 0 < q < q2. Имеем y1 (T) = -(z1 + P1 + P2 )qe-(P+P2 )T < 0 , T > 0 . Тогда f/(T) < 0, T > 0, т.е. f1(T) - убывающая функция переменной T , T > 0; убывает от /1(0) > 0 до У1(да) > 0 . Тогда f1(T) > 0 , T > 0 . Итак, если q<q) > 0, q1 = 0 < q < q2, то f1(T) > 0 , T > 0 .
Случай 3. ф(q) > 0, q1 <0 <q <q2. Имеем y1(T) = zz(z -P1 -P2)-(z1 +P1 +P2)qe-(P +P2)T <0, T>0. Тогда f1(T) < 0 , T > 0 ; далее ход доказательства аналогичен случаю 2. Итак, если ф^) > 0, q1 < 0 < q < q2, то f1(T) > 0 , T > 0 .
Случай 4. ф(д) > 0, q1 < q = 0 < q2. Имеем y1(T) = z1 z2 (z1 -P1 -P2 ) = const < 0, T > 0. Тогда f1(T) < 0 , T > 0 ; далее ход доказательства аналогичен случаю 2. Итак, если ф(д) > 0, q1 < q = 0 < q2, то f1(T) > 0, T > 0.
Случай 5. ф(q) > 0, q1 < q < 0 < q2. Имеем y1 (T) = z1z2 (z1 - P1 - P2) - (z1 + P1 + P2 )qe-(P1+P2 )T, T > 0 . В данном случае функция y1 (T) является убывающей функцией переменной T , T > 0 ; убывает от y1(T = 0,q) до уДда) = ^2((-P1 P2)<0. Рассмотрим y1(T = 0,q) = z1z2(z1 -P1 -P2)-(z +P1 +P2)q,
q1 < q < 0 . Тогда: а) y^T = 0,q = q1) = -z2(P1 +P2) -P1 -P2)> 0; б) ^(T = 0,q = 0) = zz(z1 -P1 -P2)< 0;
* * / \ / \ * *
в) y1(T = 0,q = q1 ) = 0, q1 = z1z2(z1 -P1 -P2)/(z1 +P1 +P2), q1 < q1 < 0. Тогда для q1 < q < q1 < 0 имеем: 1) y1 (T) > 0, 0 < T < T1(q); 2) y1 (T) = 0, T = T1(q); 3) y1 (T) < 0, T > T1(q), где T1 (q) определяется выражением (10). Следовательно, для q1 < q < q1 < 0 имеем: 1) f/(T) > 0, 0 < T < T1(q); 2) f/(T) = 0, T = T1(q); 3) f1(T) < 0 , T > T1 (q) . Таким образом, в точке T = T1 (q) реализуется единственный максимум функции f1(T), в силу этого f1(T) > 0 , T > 0 . Итак, если q>(q) > 0, q1 < q < q1 < 0, то f1(T) > 0 , T > 0 .
Для q = q1* имеем y1(T) = z1 z2(z1 -P1 -P2)[1 -e-(P1+P2)T]< 0, T > 0. Тогда f'(T) < 0, T > 0, т.е. f1 (T) - убывающая функция переменной T , T > 0 ; убывает от f1 (0) > 0 до f1 (да) > 0 . Тогда f1 (T) > 0 ,
T > 0 . Итак, если q>(q) > 0, q1 < q = q1 < 0 , то f1(T) > 0 , T > 0 .
*
Для q1 < q < 0 экстремальной точки T = T1 (q) не существует (T > 0), при этом y1(T = 0, q) < 0, следовательно, y1(T) < 0, T > 0 . Далее ход доказательства аналогичен случаю 2. Итак, если ф^) > 0, q1 < q < 0, то f1(T) > 0, T > 0.
Случай 6. ф(q) > 0, q1 < q < q2 = 0. Имеем y1 (T) = - z1z2 (z2 - z1) - (z1 + z2 )qe-hT , T > 0 . В данном случае функция y1 (T) является убывающей функцией переменной T , T > 0 ; убывает от y1(T = 0, q) до y1 (да) = -z1z2 (z2 - z1) < 0 . Рассмотрим y1 (T = 0, q) = -[z1z2 (z2 - z1)+(z1 + z2 )q], q1 < q < q2 = 0 . Тогда: а) y1(T = 0, q = q1) =-z22 (z1 -P1 -P2 )> 0; б) y1(T = 0, q = q2 = 0) = z2 (z2 - z1 )< 0; в) y1(T = 0, q = q^) = 0, q1* =-z1z2 (z2 - z1 )/(z1 + z2), q1 < q1 < 0. Далее ход доказательства аналогичен случаю 5. Итак, если ф^) > 0, q1 < q < q2 = 0 , то f1(T) > 0 , T > 0 .
Случай 7. ф^) > 0, q1 < q < q2 < 0. Имеем y1 (T) = z1z2 (z1 - P1 - P2) - (z1 + P1 + P2 )qe-(P1+P2 )T, T > 0 . В данном случае функция y1 (T) является убывающей функцией переменной T , T > 0 ; убывает от y1(T = 0, q) до У1 (да) = z1z2 (z1 -P1 -P2 )< 0. Рассмотрим y1(T = 0, q) = zjz2 (z1 -Pj-P2 )-(z1 +Pj +P2) q, q1 < q < q2 < 0. Тогда: 1) y,(T = 0,q = q^ =-z2(P1 +P2)x (zt-P1 -P2)> 0; 2) ^(T = 0,q = q2) = = -z1 (P1 + P2)[z2 - z1 + (z2 - P1 - P2)] , здесь z2 - P1 - P2 < 0. Далее возможны варианты: а) y1(T = 0, q = q2) > 0; б) y1(T = 0, q = q2) = 0; в) y1(T = 0, q = q2) < 0.
Вариант а: у1(Т = 0,q = q2) > 0. В данном варианте функция у1(Т) один раз пересекает ось абсцисс в точке T1(q), определенной выражением (10). Тогда поведение производной //(Т) будет выглядеть следующим образом: 1) //(Т) > 0 , 0 < Т < T1(q); 2) //(Т) = 0 , Т = T1(q); 3) //(Т) < 0 , Т > T1(q) . Таким образом, в точке Т = T1 (q) реализуется единственный максимум функции /1 (Т). Тогда /1 (Т) > 0 , Т > 0 , так как /1 (0) > 0, /1 (да) > 0 , /1 (0) > /1 (да). Итак, если 9(q) > 0, q1 < q < q2 < 0, У1(Т = 0, q = q2) > 0, то /1(Т) > 0, T > 0.
Вариант б: у1(Т = 0, q = q2) = 0. Результат идентичен результату варианта а.
Вариант в: у1(Т = 0, q = q2) < 0 . Результат идентичен результату в случае 5.
Итак, если 9(q) > 0, q1 < q < q2 < 0, то /1(Т) > 0 , Т > 0 .
Суммируя результаты случаев 1-7, получаем: если 9(q) > 0, q1 < q < q2, где q1 = z2 ( -P1 -P2), q2 = z(-P1 -P2), то /1(Т) > 0, T > 0.
Перейдем к рассмотрению случаев, когда 9(q) < 0 . Так как q > q1, область q < q1 отсекается.
Случай 8. 9(q) < 0, 0 < q1 < q2 < q. Имеем y1(T) = z^z2(z1 -P1 -P2)-(z1 +P1 +P2)qe-(P1+P2)T, T > 0. В данном случае функция y1 (Т) является возрастающей функцией переменной Т , Т > 0 ; возрастает от У1(Т = 0, q) до У1 (да) = Z1Z2 ( -P1 -P2 )> 0. Рассмотрим у1(Т = 0, q) = Z1Z2(z-P1 -P2)- (z +Pj +P2) q, 0 < q < q2 < q. Тогда y1(T = 0, q = q2) < 0 и: 1) уДТ) < 0, 0 < T < T^q); 2) y1(T) = 0, T = T1(q); 3) y1(T) > 0, T > T1(q), где T1(q) определяется выражением (10). Следовательно: 1) f(T) > 0; 0 < T < T1(q); 2) /(T) = 0 , T = T1(q); 3) /(T) < 0 , T > T1(q) . Таким образом, в точке Т = T1(q) реализуется максимум функции / (Т), при этом / (0) > 0, / (да) > 0 , /(0) < / (да). Отсюда следует, что /1(Т) > 0 , Т > 0 . Итак, если q<q) < 0 , 0 < q1 < q2 < q, то /1(Т) > 0 , Т > 0 .
Случай 9. 9(q) < 0, q1 = 0 < q2 < q . Имеем y1 (T) = -(z1 + P1 + P2 )qe-(P +P2 )T < 0 , T > 0 . Тогда /1(Т) > 0 , T > 0 , т.е. /1(Т) - возрастающая функция переменной Т , Т > 0 ; возрастает от /1(0) > 0 до /1(да) > 0 . Тогда /1(Т) > 0 , Т > 0 . Итак, если q<q) < 0 , q1 = 0 < q2 < q , то /1(Т) > 0 , Т > 0 .
Случай 10. 9(q) < 0, q1 < 0 < q2 < q. Имеем y1(T) = z1z2 (z1 -P1 -P2 )-(z1 +P1 +P2 )qe-(P1+P2 < 0, T > 0 . Тогда /1(Т) > 0, T > 0 . Далее ход доказательства аналогичен случаю 9. Итак, если 9(q) < 0 , q < 0 < q2 < q, то /1(Т) > 0, Т > 0.
Случай 11. 9(q) < 0, q1 < q2 < 0 < q. Имеем y1(T) = z1z2 (z1 -P1 -P2 )-(z1 +P1 +P2 )qe-(P1+P2 < 0, T > 0 . Тогда /1(Т) > 0, T > 0 . Далее ход доказательства аналогичен случаю 9. Итак, если 9(q) < 0 , q1 < q2 < 0 < q, то /1(Т) > 0 , Т > 0 .
Случай 12. 9(q) < 0, q1 < q2 < q = 0. Имеем yj(T) = zjz2 (zj -Pj -P2 ) = const < 0, T > 0. Тогда /1(Т) > 0 , T > 0 . Далее ход доказательства аналогичен случаю 9. Итак, если 9(q) < 0 , q1 < q2 < q = 0 , то /1(Т) > 0, Т > 0.
Случай 13. 9(q) < 0, q1 < q2 < q < 0. Имеем y1(T) = z1z2(z1 -P1 -P2)-(z1 +P1 +P2)qe-(P1+P2)T, T > 0. В данном случае функция y1 (Т) является убывающей функцией переменной Т , Т > 0 ; убывает от у1(Т = 0,q) до уДда) = z1z2(z1 -P1 -P2)<0. Рассмотрим у1(Т = 0,q) = z1z2(z1 -P1 -P2)- (z1 +P1 +P2)q, q2 < q < 0. Тогда: 1) y1(T = 0,q = q2) = - z1 (P1 +P2\z2-z1 +(z2-P1 -P2)]; 2) y1(T = 0,q = 0) = = z1 z2 (z1 -P1 P2 )< 0.
*
Подслучай 13.1. y1(T = 0, q = q2) > 0. Тогда возможны варианты: а) у1(Т = 0, q) > 0, q2 < q < q1 < 0 ; б) y1(T = 0, q) = 0, q = q/ < 0; в) Л(Т = 0, q) < 0, q/ < q < 0, q/ = z1z2 (z-P1 -P2)/ (z1 +P1 +P2) .
Вариант а: yj(T = 0,q) > 0, q2 < q < qj < 0 . В данном варианте функция yj (T) один раз пересекает ось абсцисс в точке Tj(q), определенной выражением (10). Тогда поведение производной fj(T) будет выглядеть следующим образом: 1) fJ(T) < 0, 0 < T < Tj(q); 2) fJ'(T) = 0, T = Tj(q); 3) fJ'(T) > 0,
T > Tj (q) . Таким образом, в точке T = Tj (q) реализуется единственный минимум функции f (T), T > 0 .
*
Можно показать, что f (Tj (q)) > 0, q2 < q < qj < 0, тогда f (T) > 0 , T > 0 . Итак, если p(q) < 0 ,
*
q2 < q < qj < 0 , то fJ(T) > 0 , T > 0 .
Вариант б: yJ(T = 0,q) = 0, q = qj < 0 . В данном варианте yj(T) является убывающей функцией переменной T , T > 0 ; убывает от yJ(T = 0, q = qj ) до yj (да) = zjz2(zj - ßj -ß2) < 0. Тогда: J) yj (T) < 0, T > 0 ; 2) yJ(T) = 0, T = 0 . Таким образом: J) fJ'(T) > 0 , T > 0 ; 2) fJ'(T) = 0 , T = 0 . Отсюда следует, что
f (T) - возрастающая функция переменной T , T > 0, и неубывающая функция в точке T = 0 ; f (T)
*
возрастает от f (0) > 0 до f (да) > 0 . Тогда f (T) > 0 , T > 0 . Итак, если p(q) < 0 , q = qJ < 0 , то fj(T) > 0, T > 0.
Вариант в: yJ(T = 0,q) < 0, qJ < q < 0 . В данном варианте yJ(T) является убывающей функцией переменной T , T > 0 ; убывает от yJ (T = 0, q) < 0 до yJ (да) = zJz2 (zJ - ßJ - ß2) < 0. Тогда yJ (T) < 0, T > 0 , следовательно, fJ'(T) > 0 , T > 0 . Отсюда следует, что f (T) - возрастающая функция переменной T ,
T > 0 ; f (T) возрастает от f (0) > 0 до f (да) > 0 . Тогда f (T) > 0 , T > 0 . Итак, если p(q) < 0 ,
*
qJ < q < 0, то fJ(T) > 0, T > 0.
Таким образом, если p(q) < 0 , qJ < q2 < q < 0, yJ(T = 0, q = q2) > 0, то f (T) > 0 , T > 0 . Подслучай 13.2. yJ(T = 0,q = q2) = 0. В данном подслучае функция yJ(T) является убывающей функцией переменной T , T > 0; убывает от yJ(T = 0,q) < 0 до yJ(<^) = zJz2 (zj-ßj-ß2 )< 0. Тогда yJ(T) < 0, T > 0, и дальнейший ход доказательства аналогичен случаю 9. Итак, если p(q) < 0, qj < q2 < q < 0, yj(T = 0, q = q2) = 0, то fj(T) > 0, T > 0.
Подслучай 13.3. yJ(T = 0,q = q2) < 0 . Дальнейший ход доказательства аналогичен подслучаю J3.2. Итак, если p(q) < 0, qJ < q2 < q < 0, yJ(T = 0, q = q2) < 0, то fJ(T) > 0, T > 0.
Суммируя подслучаи J3.J-J3.3, получаем: если p(q) < 0, qJ < q2 < q < 0, то fj(T) > 0, T > 0.
Суммируя случаи 8-J3, получаем: если 9(q) < 0 , qj < q2 < q, где qj = z2(zj -ßj -ß2), q2 = zj(z2 - ßj -ß2 ), то fj(T) > 0, T > 0.
Таким образом, j) если p(q) > 0, qj < q < q2, то f (T) > 0 , T > 0 ; 2) если p(q) = 0, q = q2, то
fj(T) = const > 0, T > 0; 3) если p(q) < 0, qj < q2 < q, то fj(T) > 0, T > 0, где qj = z2(zj-ßj-ß2),
q2 = zj (z2 - ßj - ß2) - любые вещественные числа. Лемма 3 доказана.
Утверждение 2. Для функции f2(T) = z2[zjz2(ßj +ß2)- z2f (T) - f '(T)], T > 0 , справедливо:
z C z^C
j) f2(0) = (ßj +ß2))cr(zjA - C) = (ßj +ß2))r(q - q2), q2 = ^2-ßj-ß2); при этом /2(0) может быть:
а) f2(0) > 0; б) f2(0) = 0; в) f2(0) < 0; q > qj, qj < q2;
2) У2(да) = lim f2(T) = z2 (zj A - C) = z^ (q - q2); при этом У2(да) может быть: а) f2(<») > 0; б) f2(<») = 0; в) У2(да) < 0;
3) f'(T) = zjzi(ß;7!)ß2}cp(qK(ft+ß2[zjz2(z2 -ßj - ß2)-(z2 +ßj +ß2)qe-(( +ß2], T > 0 ;
4) точка T2(q) - точка экстремума (либо точка перегиба) функции f2(T):
T2(q) _-
1 ln ZiZ2 (Z2 -ßl -ß2 ) Z2 -ßl -ß2
ßi +ß2 (( +ßi +ß2 )q
> 0:
5) f2(T2(q)) _
_ Z2 ( +ßi +ß2 )2
4(ßi +ß2 )q
r-zi(z2 -ßi -ß2 )
Z2 (Zi -ßi -ß2 )(Z2 -ßi -ß2 )2 (Z2 +ßi +ß2 )2
Z2 -ßi -ß2
> 0,
где q22 _ q2 _ Zi(Z2 - ßi -ß2), q2i _ Z2(zi ßi ß)z2 ßi ß2) - корни уравнения f2(T2(q)) _ 0 . Здесь q,
(Z2 +ßi +ß2)
F(T) определены в (i); A, C - в (4); 9(q) - в (7).
Лемма 4. Функция f2(T) является отрицательной функцией (f2(T) < 0) переменной T , T > 0, если 9(q) > 0, qi < q < q2, и положительной функцией (f2 (T) > 0) переменной T , T > 0 , если 9(q) < 0 , qi < q2 < q, где qi _ z2 (z1 -ßi -ß2), q2 _ z1(z2 -ßi -ß2) - любые вещественные числа. Доказательство аналогично доказательству леммы 3.
Рассмотрим функцию 0(T) _ fi(T) - f2(T) как функцию переменной T , T > 0 . Утверждение 3. Для функции 0(T) _ fi(T) -f2(T), T > 0 , справедливо:
1) O(T) _(Z2 - Zi )[-ZiZ2 (ßi +ß2 )+(zi + Z2 )f(T) + f'(T)], T > 0;
2) Ф(0) _ (z2 - zi )(j3i + ß2)C / A)2 > 0, следовательно, f (0) > f2 (0);
3) Ф(<х>) _ lim Ф(T) _ (z2 - zi )С > 0, следовательно, fi(x>) > f2(<»);
t ^вд
4) Ф(0) -Ф(вд) _(z2 - Zi )cp(q)C / A2;
5) Ф'(T) _ -(Z2 - zi)q) ZiZ2(ßl FßT)(ft+ß2[ziZ2(zi + Z2 -ßi -ß2)-(( + Z2 +ßi +ß2)+ß2)) ], T > 0 ;
6) точка T* (q) - точка экстремума (либо точка перегиба) функции Ф^):
T*(q) _-
i_ln ZiZ2 (zi + Z2 -ßi ß2 ) Zi + Z2 -ßi ß2
ßi +ß2 (Zi + Z2 +ßi +ß2 )q
>0;
7) Ф(T*(q)) _(Z2 - zj
C + 9(q)
(zi + Z2-ßi-ß2 )2
4(ßi +ß2 )q
Zi + Z2 -ßi ß2 q
>0.
Здесь д , F(Т) определены в (1); А, С - в (4); ф(д) - в (7).
Лемма 5. Функция Ф(Т) является положительной функцией (Ф(Т) > 0) переменной Т , Т > 0 , т.е. Л(Т) > /2(Т), Т > 0.
Доказательство аналогично доказательству леммы 3. Леммы 1-5 позволяют сформулировать следующую теорему.
Теорема 1. Производная р'Т (тт) есть положительная функция (р'Т (тт) > 0) переменной Т , 0 < Т <тт , 0 <тт <да, при любых значениях параметров Л,1 >Х 2 > 0, 0 < р < 1, р> 0, а> 0, 0 <5< 1. Доказательство. Доказательство осуществляется последовательным применением лемм 1-5:
1. Если ф(д) > 0, ^ < д < д2, д ф д2, то /1(Т) > 0 , /2(Т) < 0 , Т > 0 . Тогда имеем /1(Т)//2(Т) < 0, Т > 0, в том числе и для 0 < Т <тт . Таким образом, уравнение (9) решения не имеет, т.е. функция рТ(Тт) не достигает нуля для 0 < Т <Тт . Так как р'Т=0(Тт) > 0, р'Т=%т (Хт) > 0, то р'Т(Тт) > 0 , 0 < Т < Тт .
2. Если ф(д) = 0, д = д1, д Ф д2, то /1 (Т) = 0 , /2(Т) < 0 , Т > 0 . Тогда имеем /1 (Т) //2(Т) = 0, Т > 0 , в том числе и для 0 < Т < тт . Таким образом, уравнение (9) решения не имеет. Тогда р'Т (тт) > 0 ,
0 < Т <Тт .
3. Если ф(д) = 0, д Ф д1, д = д2, то /1 (Т) > 0 , /2 (Т) = 0, Т > 0 . Тогда в уравнении (9) имеет место деление на ноль. Преобразуем (9) к виду
)(Т" -Т]= /2 (Т) / /1(Т), 0 < Т <Тт . (11)
q
Тогда имеем f2(T)/ f1(T) = 0, T > 0 , в том числе и для 0 < T < xm . Следовательно, уравнение (11) решения не имеет, уравнение (9) также не имеет решения. Тогда p'T (xm) > 0 , 0 < T < xm .
4. Если q>(q) < 0, q > q2, то f1(T) > 0, f2(T) > 0, f1(T) > f2(T), T > 0. Тогда имеем f1 (T) / f2 (T) > 1, T > 0 , в том числе и для 0 < T < xm . Таким образом, уравнение (9) решения не имеет. Тогда p'T (xm) > 0 , 0 < T < xm . Теорема 1 доказана.
Теорема 2. Плотность вероятности pT (xm) есть возрастающая функция переменной T, 0 < T <xm , 0 <xm <да, при любых значениях параметров Л,1 > А, 2 > 0, 0 < p < 1, P> 0, а> 0, 0 <5< 1 и достигает своего максимального значения в точке T = xm, 0 < xm < да . Доказательство вытекает из результата теоремы 1.
Следствие 1. Из теоремы 1 вытекает, что функции pT(x(j)), j = 2,k , являются возрастающими функциями переменной T , 0 < T < x(j), 0 < x(j) <да , j = 2, k , при любых значениях параметров X1 >А, 2 > 0, 0 < p < 1, P> 0, а > 0, 0 < 5 < 1.
Следствие 2. Из теоремы 2 вытекает, что функция правдоподобия L(T | x(1),...,x(k)) достигает своего глобального максимума в точке T = xm, т.е. решением оптимизационной задачи (2) является оценка длительности мертвого времени T = xm .
Заключение
Полученный результат делает возможным решение задачи оценивания длительности мертвого времени без привлечения численных методов: в процессе наблюдения за потоком событий (в течение временного интервала (0, t]) вычисляются величины xk, k = 1, n , после чего находится xmin = min xk (k = 1, n ), и полагается T = xmin . Подчеркнем, что по определению оценка максимального правдоподобия длительности мертвого времени при конечных t будет всегда смещенная (xmin > T ); ее несмещен-ность реализуется только в асимптотическом случае при t ^ да .
ЛИТЕРАТУРА
1. Cox D.R. Some Statistical Methods Connected with Series of Events // J. Royal Statistical Society B. 1955. V. 17. P. 129-164.
2. Kingman Y.F.C. On doubly stochastic Poisson process // Proceedings of Cambridge Phylosophical Society. 1964. V. 60, No 4.
P. 923-930.
3. Basharin G.P., Kokotushkin V.A., Naumov V.A. Method of equivalent substitutions for calculating fragments of communication
networks for digital computer // Engineering cybernetics. 1979. V. 17(6). P. 66-73.
4. Башарин Г.П., Кокотушкин В.А., Наумов В.А. О методе эквивалентных замен расчета фрагментов сетей связи // Известия
АН СССР. Техн. кибернетика. 1980. № 1. С. 55-61.
5. Neuts M.F. A versatile Markov point process // Journal of Applied Probability. 1979. V. 16. P. 764-779.
6. Cox D. R., Isham V. Point Processes. London : Chapman & Hall, 1980.
7. Bremaud P. Point Processes and Queues: Martingale Dynamics. N.Y. : Springer-Verlag, 1981.
8. Last G., Brandt A. Marked Point Process on the Real Line: The Dynamic Approach. N.Y. : Springer-Verlag, 1995.
9. Gortsev A.M., Nezhelskaya L.A. An asynchronous double stochastic flow with initiation of superfluous events // Discrete Mathemat-
ics and Applications. 2011. V. 21, No. 3. P. 283-290.
10. Башарин Г.П., Гайдамака Ю.В., Самуйлов К.Е. Математическая теория телетрафика и ее приложения к анализу мультисер-висных сетей связи следующих поколений // Автоматика и вычислительная техника. 2013. № 2. С. 11-21.
11. Adamu A., Gaidamaka Y., Samuylov A. Discrete Markov Chain Model for Analyzing Probability Measures of P2P Streaming Network // Lecture Notes in Computer Science: Proc. of the 11-th International Conference on Next Generation Wired/Wireless Networking NEW2AN-2011 (August 23-25, 2011, St. Petersburg, Russia). 2011. P. 428-439.
12. Bouzas P.R., Valderrama M.J., Aguilera A.M., Ruiz-Fuentes N. Modelling the mean of a doubly stochastic Poisson process by functional data analysis // Computational Statistics and Data Analysis. 2006. V. 50(10). P. 2655-2667.
13. Centanni S., Minozzo M. A Monte Carlo approach to filtering for a class of marked doubly stochastic Poisson processes // Journal of the American Statistical Association. 2006. V. 101. P. 1582-1597.
14. Dubois J.-P. Traffic estimation in wireless networks using filtered doubly stochastic point processes (Conference Paper) // Proceedings - 2004 International Conference on Electrical, Electronic and Computer Engineering, ICEEC'04 2004. 2004. P. 116-119.
15. Hossain M.M., Lawson A.B. Approximate methods in Bayesian point process spatial models // Computational Statistics and Data Analysis. 2009. V. 53(8). P. 2831-2842.
16. Snyder D.L., Miller M.I. Random Point Processes in Time and Space. Springer-Verlag, Heidelberg, 1991.
17. Горцев А.М., Нежельская Л. А., Шевченко Т.И. Оценивание состояний MC-потока событий при наличии ошибок измерений // Известия вузов. Физика. 1993. № 12. С. 67-85.
18. Gortsev A.M., Shmyrin I.S. Optimal estimation of states of a double stochastic flow of events in the presence of measurement errors of time instants // Automation and Remote Control. 1999. V. 60, No. 1. P. 41-51.
19. Горцев А.М., Шмырин И.С. Оптимальная оценка состояний дважды стохастического потока событий при наличии ошибок в измерениях моментов времени // Автоматика и телемеханика. 1999. № 1. С. 52-66.
20. Горцев А.М., Ниссенбаум О.В. Оптимальная оценка состояний асинхронного альтернирующего потока с инициированием лишних событий // Вестник ТюмГУ. 2008. № 6. С. 107-119.
21. Горцев А.М., Зуевич В. Л. Оптимальная оценка состояний асинхронного дважды стохастического потока событий с произвольным числом состояний // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2010. № 2(11). С. 44-65.
22. Горцев А.М., Леонова М.А. Оптимальная оценка состояний обобщенного асинхронного дважды стохастического потока // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2010. № 1(10). С. 33-47.
23. Gortsev A.M., Nezhelskaya L.A., Solovev A.A. Optimal State Estimation in MAP Event Flows with Unextendable Dead Time // Automation and Remote Control. 2012. V. 73, No. 8. P. 1316-1326.
24. Бахолдина М.А. Оптимальная оценка состояний модулированного обобщенного полусинхронного потока событий // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2013. № 2(23). С. 10-21.
25. Бахолдина М.А., Горцев А.М. Оптимальная оценка состояний модулированного обобщенного полусинхронного потока событий при непродлевающемся мертвом времени // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2014. № 1(26). С. 13-24.
26. Bakholdina M.A., Gortsev A.M. Optimal estimation of the states of modulated semi-synchronous integrated flow of events in condition of its incomplete observability // Applied Mathematical Sciences. 2015. V. 9, No. 29. P. 1433-1451.
27. Горцев А. М., Завгородняя М. Е. Оценка параметров альтернирующего потока событий при условии его частичной наблюдаемости // Оптика атмосферы и океана. 1997. Т. 10. № 3. С. 273-280.
28. Васильева Л.А., Горцев А.М. Оценивание параметров дважды стохастического потока событий в условиях его неполной наблюдаемости // Автоматика и телемеханика. 2002. № 3. С 179-184.
29. Горцев А.М., Нежельская Л.А. Оценивание длительности «мертвого времени» и интенсивностей синхронного дважды стохастического потока событий // Радиотехника. 2004. № 10. С. 8-16.
30. Горцев А. М., Ниссенбаум О. В. Оценивание длительности мертвого времени и параметров асинхронного альтернирующего потока событий с инициированием лишнего события // Вестник Томского государственного университета. 2004. № 284. С. 137-145.
31. Горцев А.М., Зуевич В. Л. Оптимальная оценка параметров асинхронного дважды стохастического потока событий с произвольным числом состояний // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2011. № 4(17). С. 25-40.
32. Леонова М. А., Нежельская Л. А. Оценка максимального правдоподобия длительности мертвого времени в обобщенном асинхронном потоке событий // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2013. № 2(23). С. 54-63.
33. Горцев А. М., Соловьев А. А. Оценка максимального правдоподобия длительности непродлевающегося мертвого времени в потоке физических событий // Известия вузов. Физика. 2015. Т. 58, № 11. С. 141-149.
34. Феллер В. Введение в теорию вероятностей и ее приложения. М. : Мир, 1967. Т. 1.
35. Баруча-Рид А.Т. Элементы теории марковских процессов и их приложения. М. : Наука, 1969.
36. Апанасович В.В., Коляда А.А., Чернявский А.Ф. Статистический анализ случайных потоков в физическом эксперименте. Минск : Университетское, 1988.
37. Normey-Rico J.E. Control of dead-time processes. (Advanced textbooks in control and signal processing). London : SpringerVerlag, 2007.
38. Горцев А. М., Климов И. С. Оценка интенсивности пуассоновского потока событий в условиях частичной его ненаблюдаемости // Радиотехника. 1991. № 12. С. 3-7.
39. Горцев А. М., Климов И. С. Оценивание периода ненаблюдаемости и интенсивности пуассоновского потока событий // Радиотехника. 1996. № 2. С. 8-11.
40. Горцев А.М., Паршина М.Е. Оценивание параметров альтернирующего потока событий в условиях «мертвого времени» // Известия вузов. Физика. 1999. № 4. С. 8-13.
41. Горцев А. М., Ниссенбаум О. В. Оценивание длительности мертвого времени и параметров асинхронного альтернирующего потока событий при непродлевающемся мертвом времени // Известия вузов. Физика. 2005. № 10. С. 35-49.
42. Gortsev A.M., Nissenbaum O.V. Estimation of the dead time period and parameters of an asynchronous alternative flow of events with unextendable dead time period. Russian Physics Journal. 2005. V. 48(10). Р. 1039-1054.
43. Нежельская Л.А. Оптимальное оценивание состояний полусинхронного потока событий в условиях его частичной наблюдаемости // Вестник Томского государственного университета. 2000. № 269. С. 95-98.
44. Горцев А.М., Нежельская Л.А. Оценивание параметров полусинхронного дважды стохастического потока событий методом моментов // Вестник Томского государственного университета. 2002. № 1. С. 18-23.
45. Gortsev A.M., Nezhelskaya L.A. Estimation of the dead-time period and parameters of a semi-synchronous double-stochastic stream of events // Measurement Techniques. 2003. V. 46, No. 6. P. 536-545.
46. Горцев А.М., Нежельская Л.А. Полусинхронный дважды стохастический поток событий при продлевающемся мертвом времени // Вычислительные технологии. 2008. Т. 13, № 1. С. 31-41.
47. Горцев А.М., Калягин А. А. Оптимальная оценка состояний обобщенного полусинхронного потока событий в условиях непродлевающегося мертвого времени // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2010. № 4(13). С. 50-60.
48. Горцев А.М., Калягин А.А., Нежельская Л.А. Совместная плотность вероятностей длительности интервалов обобщенного полусинхронного потока событий при непродлевающемся мертвом времени // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2014. № 2(27). С. 19-29.
49. Горцев А.М., Калягин А.А., Нежельская Л.А. Оценка максимального правдоподобия длительности мертвого времени в обобщенном полусинхронном потоке // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2015. № 1(30). С. 27-37.
50. Горцев А.М., Нежельская Л.А. О связи MC-потоков и MAP-потоков событий // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2011. № 1(14). С. 13-21.
51. Бахолдина М. А. Совместная плотность вероятностей длительности интервалов модулированного обобщенного полусинхронного потока событий при непродлевающемся мертвом времени и условия его рекуррентности // Вестник Томского государственного университета. Управление, вычислительная техника и информатика. 2015. № 2(31). С. 4-17.
52. Bakholdina M., Gortsev A. Joint probability density of the intervals length of the modulated semi-synchronous integrated flow of events and its recurrence conditions // Communications in Computer and Information Science. 2014. V. 487. P. 18-25.
53. Бахолдина М. А., Горцев А. М. Совместная плотность вероятностей длительности интервалов модулированного обобщенного полусинхронного потока событий и условия его рекуррентности // Информационные технологии и математическое моделирование (ИТММ-2014): материалы XIII Междунар. науч.-практ. конф. им. А.Ф. Терпугова (20-22 ноября 2014 г.). Томск : Изд-во Том. ун-та, 2014. Ч. 2. С. 137-143.
54. Бахолдина М. А., Горцев А. М. Плотность вероятностей длительности интервала между соседними событиями модулированного обобщенного полусинхронного потока событий при непродлевающемся мертвом времени // Теория вероятностей, случайные процессы, математическая статистика и приложения: материалы Междунар. науч. конф., посвящ. 80-летию проф., д-ра физ.-мат. наук Г.А. Медведева, Минск, 23-26 фев. 2015 г. / редкол.: Н.Н. Труш [и др.]. Минск : РИВШ, 2015. С. 17-22.
55. Bakholdina M., Gortsev A. Joint probability density of the intervals length of modulated semi-synchronous integrated flow of events in conditions of a constant dead time and the flow recurrence conditions // Communications in Computer and Information Science. 2015. V. 564. P. 13-27.
Бахолдина Мария Алексеевна. E-mail: [email protected]
Горцев Александр Михайлович, д-р техн. наук, профессор. E-mail: [email protected]
Томский государственный университет
Поступила в редакцию 19 февраля 2016 г.
Bakholdina Maria A., Gortsev Alexander M. (Tomsk State University, Russian Federation).
Maximum likelihood estimation of unextendable dead time period duration in the modulated semi-synchronous generalized flow of events.
Keywords: modulated semi-synchronous generalized flow of events; doubly stochastic Poisson process (DSPP); Markovian arrival process (MAP); maximum likelihood estimation; likelihood function; dead time period duration.
DOI: 10.17223/19988605/35/2
In this paper we consider the modulated semi-synchronous generalized flow of events, which is one of the mathematical models for incoming streams of events in computer communication networks and which is related to the class of doubly stochastic Poisson processes (DSPPs). The flow intensity process is a piecewise constant stationary random process X(t) with two states 1, 2 (first, second correspondingly). In the first state X(t) = Xj and in the second state X(t) = X2 (X1 >X2 > 0 ). During the time interval of a random duration when the process X(t) is in state (X(t) = ), a Poisson flow of events with intensity , i = 1,2, arrives. The transition of the process X(t) from the first state to the second state is possible at any moment of a Poisson event occurrence in state 1 of the process X(t) , herewith the process X(t) can change its state to the second one with probability p (0 < p < 1) or continue to stay in state 1 with complementary probability 1 — p . The transition of the process X(t) from state 1 to state 2 is also possible at any moment that does not coincide with the moment of a Poisson event occurrence, herewith the duration of the process X(t) staying in the first state is distributed according to the exponential law with parameter P : F(x) = 1 — e—Px, x > 0 . Then the duration of the process X(t) staying
in the first state is distributed according to the exponential law with distribution function F1(x) = 1 - e-(pX1 +P'x, x> 0 . The transition of the process X(t) from the second state to the first state at the moment of a Poisson event occurrence in state 2 is impossible and can be done only at a random time moment. In this case the duration of the process X(t) staying in state 2 is distributed according to the exponential law with parameter а : F2 (x) = 1 — e—ax, x> 0 . At the moment when the state changes from the second to the first one, an additional event is assumed to be initiated with probability 5 (0 < 5 < 1) .
The registration of the flow events is considered in conditions of a constant (unextendable) dead time. The dead time period of a constant duration T begins after every registered at the moment tk , k > 1, event. During this period no other events are observed. When the dead time period is over, the first coming event causes the next interval of a dead time of duration T and so on.
This paper contains analytical results that are devoted to finding the maximum likelihood estimate T of the dead time period duration on monitoring the time moments of the events occurrence. We assume that the flow parameters Xj > > 0, 0 < p < 1, P > 0, a> 0, 0 <5< 1 are known and the duration of the dead time period T is not known. According to the maximum-likelihood technique the likelihood function L(T | x(1),...,x(k)) is maximized and the following task of optimization is solved:
L(T | x(1),..., x(k)) = f\pT (j ^ max, 0 < T <xmm, xmm > 0, j=1
where pT (x(j-1) is the one-dimensional probability density function of the interval length between two consecutive flow events. Finally, we obtain that the likelihood function L(T |x(1),..., x(k > ) reaches its global maximum at the point T = x min, where xmin = min xk (k = 1 , n ); %k — tk+1 — tk , x > 0 , k — 1, n - the sequence of the values of the intervals lengths between consecutive flow events measured during the interval of observation (0, t] , i.e. the solution of optimization problem is the estimate of the dead time period duration: T = xmm .
REFERENCES
1. Cox, D.R. (1955) Some Statistical Methods Connected with Series of Events. Journal of Royal Statistical Society B. 17. pp. 129-164.
DOI: 10.2307/2983950
2. Kingman, Y.F.C. (1964) On doubly stochastic Poisson process. Proceedings of Cambridge Phylosophical Society. 60(4). pp. 923-930.
3. Basharin, G.P., Kokotushkin, V.A. & Naumov, V.A. (1979) Method of equivalent substitutions for calculating fragments of commu-
nication networks for digital computer. Engineering cybernetics. 17(6). pp. 66-73. DOI: 10.1016/0166-5316(84)90009-9
4. Basharin, G.P., Kokotushkin, V.A. & Naumov, V.A. (1980) O metode ekvivalentnykh zamen rascheta fragmentov setey svyazi. Izv.
ANSSSR. Tekhn. kibernetika. 1. pp. 55-61.
5. Neuts, M.F. (1979) A versatile Markov point process. Journal of Applied Probability. 16. pp. 764-779. DOI: 10.2307/3213143
6. Cox, D.R. & Isham, V. (1980) Point Processes. London: Chapman & Hall.
7. Bremaud, P. (1981) Point Processes and Queues: Martingale Dynamics. New York: Springer-Verlag.
8. Last, G. & Brandt, A. (1995) Marked Point Process on the Real Line: The Dynamic Approach. New York: Springer-Verlag.
9. 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.1515/dma.2011.017
10. Basharin, G.P., Gaidamaka, Y.V. & Samouylov, K.E. (2013) Mathematical Theory of Teletraffic and Its Application to the Analysis of Multiservice Communication of Next Generation Networks. Automatic Control and Computer Sciences. 47(2). pp. 62-69. DOI: 10.3103/S0146411613020028
11. Adamu, A., Gaidamaka, Y. & Samuylov, A. (2011) Discrete Markov Chain Model for Analyzing Probability Measures of P2P Streaming Network. Lecture Notes in Computer Science. Proc. of the 11-th International Conference on Next Generation Wired/Wireless Networking NEW2AN-2011. August 23-25, 2011. St. Petersburg, Russia. pp. 428-439.
12. Bouzas, P.R., Valderrama, M.J., Aguilera, A.M. & Ruiz-Fuentes, N. (2006) Modelling the mean of a doubly stochastic Poisson process by functional data analysis. Computational Statistics and Data Analysis. 50(10). pp. 2655-2667. DOI: 10.1016/j.csda.2005.04.015
13. Centanni, S. & Minozzo, M. (2006) A Monte Carlo approach to filtering for a class of marked doubly stochastic Poisson processes. Journal of the American Statistical Association. 101. pp. 1582-1597. DOI: 10.1198/016214506000000276
14. Dubois, J.-P. (2004) Traffic estimation in wireless networks using filtered doubly stochastic point processes (Conference Paper). Proceedings - 2004 International Conference on Electrical, Electronic and Computer Engineering, ICEEC'04 2004. pp. 116-119.
15. Hossain, M.M. & Lawson, A.B. (2009) Approximate methods in Bayesian point process spatial models. Computational Statistics and Data Analysis. 53(8). pp. 2831-2842. DOI: 10.1186/1471-2288-11-133
16. Snyder, D.L. & Miller, M.I. (1991) Random Point Processes in Time and Space. Heidelberg: Springer-Verlag.
17. Gortsev, A.M., Nezhelskaya, L.A. & Shevchenko, T.I. (1993) Optimal states estimation of asynchronous doubly stochastic flow of events with arbitrary number of states. Russian Physics Journal. 12. pp. 67-85. (In Russian).
18. Gortsev, A.M. & Shmyrin, I.S. (1999) Optimal estimation of states of a double stochastic flow of events in the presence of measurement errors of time instants. Automation and Remote Control. 60(1). pp. 41-51.
19. Gortsev, A.M. & Shmyrin, I.S. (1999) Optimal estimation of states of a double stochastic flow of events in the presence of measurement errors of time instants. Automation and Remote Control. 1. pp. 52-66. (In Russian).
20. Gortsev, A.M. & Nissenbaum, O.V. (2008) Optimal'naya otsenka sostoyaniy asinkhronnogo al'terniruyushchego potoka s initsiiro-vaniem lishnikh sobytiy [Optimal estimation of the states of an asynchronous alternative flow of events with additional events initiation]. Vestnik Tyumenskogo gosudarstvennogo universiteta - The Bulletin of Tyumen State University. 6. pp. 107-119.
21. Gortsev, A.M. & Zuevich, V.L. (2010) Optimal states estimation of asynchronous doubly stochastic flow of events with arbitrary 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).
22. 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).
23. Gortsev, A.M., Nezhelskaya, L.A. & Solovev, 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
24. Bakholdina, M.A. (2013) Optimal states estimation of the modulated semi-synchronous integrated flow of events in conditions of constant dead time. Vestnik Tomskogo gosudarstvennogo universiteta. Upravlenie, vychislitel'naya tekhnika i informatika - Tomsk State University Journal of Control and Computer Science. 2(23). pp. 10-21. (In Russian).
25. Bakholdina, M.A. & Gortsev, A.M. (2014) Optimal states estimation of the modulated semi-synchronous integrated flow of events in the condition of constant dead time. Vestnik Tomskogo gosudarstvennogo universiteta. Upravlenie, vychislitel'naya tekhnika i informatika - Tomsk State University Journal of Control and Computer Science. 1(26). pp. 13-24. (In Russian).
26. Bakholdina, M.A. & Gortsev, A.M. (2015) Optimal estimation of the states of modulated semi-synchronous integrated flow of events in condition of its incomplete observability. Applied Mathematical Sciences. 9(29). pp. 1433-1451. DOI: 10.12988/ams.2015.5135
27. Gortsev, A.M. & Zavgorodnyaya, 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.
28. Vasileva, L.A. & Gortsev, A.M. (2002) Estimation of parameters of a double-stochastic flow of events under conditions of its incomplete observability. Automation and Remote Control 63(3). pp. 511-515. DOI: 10.1023/A:1014718921138
29. Gortsev, A.M. & Nezhelskaya, L.A. (2004) Estimation of the dead time period and intensities of the synchronous double stochastic event flow. Radiotechnics. 10. pp. 8-16.
30. Gortsev, A.M. & Nissenbaum, O.V. (2004) Otsenivanie dlitel'nosti mertvogo vremeni i parametrov asinkhronnogo al'terniruyush-chego potoka sobytiy s initsiirovaniem lishnego sobytiya [Estimation of the dead time period and parameters of an asynchronous alternative flow of events with additional events initiation]. Vestnik Tomskogo gosudarstvennogo universiteta - Tomsk State University Journal. 284. pp. 137-145.
31. Gortsev, A.M. & Zuevich, V.L. (2011) Optimal estimation of parameters of an asynchronous doubly stochastic flow of events with arbitrary number of the states. Vestnik Tomskogo gosudarstvennogo universiteta. Upravlenie, vychislitel'naya tekhnika i informatika - Tomsk State University Journal of Control and Computer Science. 4(17). pp. 25-40. (In Russian).
32. 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).
33. Gortsev, A.M. & Solovev, A.A. (2015) Otsenka maksimal'nogo pravdopodobiya dlitel'nosti neprodlevayushchegosya mertvogo vremeni v potoke fizicheskikh sobytiy [Maximum likelihood estimation of fixed dead time in physical flow of events]. Izvestiya vuzov. Fizika. 58(11). pp. 141-149.
34. Feller, W. (1950) Vvedenie v teoriyu veroyatnostey i ee prilozheniya [An introduction to probability theory and its applications]. Translated from English. Vol. I. New York-London-Sydney: John Wiley and Sons, Inc.
35. Bharucha-Reid, A.T. (1960) Elementy teorii markovskikh protsessov i ikh prilozheniya [Elements of the theory of Markov processes and their applications]. Translated from English by V. Kalashnikov, Yu. Krutovskiy, M. Epelman. Moscow: Nauka.
36. Apanasovich, V.V., Kolyada, A.A. & Chernyavskiy, A.F. (1988) Statisticheskiy analiz sluchaynykh potokov v fizicheskom eksperi-mente [The statistical analysis of series of random events in physical experiment]. Minsk: University Press.
37. Normey-Rico, J.E. (2007) Control of dead-time processes. (Advanced textbooks in control and signal processing). London: Springer-Verlag.
38. Gortsev, A.M. & Klimov, I.S. (1991) Otsenka intensivnosti puassonovskogo potoka sobytiy v usloviyakh chastichnoy ego nena-blyudaemosti [Intensity estimation of the Poisson flow of events in condition of its incomplete observability]. Radiotekhnika - Ra-diotechnics. 12. pp. 3-7.
39. Gortsev, A.M. & Klimov, I.S. (1996) Otsenivanie perioda nenablyudaemosti i intensivnosti puassonovskogo potoka sobytiy [The estimation of intensity process and period of unobservability of the Poisson flow of events]. Radiotekhnika - Radiotechnics. 2. pp. 8-11.
40. Gortsev, A.M. & Parshina, M.E. (1999) Estimation of parameters of an alternate stream of events in "dead" time conditions. Russian Physics Journal. 4. pp. 8-13. DOI: 10.1007/BF02509672
41. Gortsev, A.M., Nissenbaum, O.V. (2005) Otsenivanie dlitel'nosti mertvogo vremeni i parametrov asinkhronnogo al'terniruyushchego potoka sobytiy pri neprodlevayushchemsya mertvom vremeni [Estimation of the dead time period and parameters of an asynchronous alternative flow of events with unextendable dead time period]. Izvestiya vuzov. Fizika. 10. pp. 35-49.
42. Gortsev, A.M. & Nissenbaum, O.V. (2005) Estimation of the dead time period and parameters of an asynchronous alternative flow of events with unextendable dead time period. Russian Physics Journal. 48(10). pp. 1039-1054. DOI: 10.1007/s11182-006-0023-y
43. 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.
44. Gortsev, A.M. & Nezhelskaya, L.A. (2002) Otsenivanie parametrov polusinkhronnogo dvazhdy stokhasticheskogo potoka sobytiy metodom momentov [Parameters estimation of a semi-synchronous doubly stochastic flow of events using method of moments]. Vestnik Tomskogo gosudarstvennogo universiteta - Tomsk State University Journal. 1. pp. 18-23.
45. 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
46. 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).
47. Gortsev, A.M. & Kalyagin, A.A. (2010) Optimal states estimation of generalized semi-synchronous flow of events in conditions of constant dead time. Vestnik Tomskogo gosudarstvennogo universiteta. Upravlenie, vychislitel'naya tekhnika i informatika - Tomsk State University Journal of Control and Computer Science. 4(13). pp. 50-60. (In Russian).
48. 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).
49. Gortsev, A.M., Kalyagin, A.A. & Nezhelskaya, L.A. (2015) Maximum likelihood estimation of dead time value at a generalized semi-synchronous flow of events. Vestnik Tomskogo gosudarstvennogo universiteta. Upravlenie, vychislitel'naya tekhnika i in-formatika - Tomsk State University Journal of Control and Computer Science. 1(30). pp. 27-37. (In Russian).
50. 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).
51. Bakholdina, M.A. (2015) Joint probability density of the intervals length of the modulated semi-synchronous integrated flow of events in conditions of constant dead time and flow recurrence conditions. Vestnik Tomskogo gosudarstvennogo universiteta. Up-ravlenie, vychislitel'naya tekhnika i informatika - Tomsk State University Journal of Control and Computer Science. 2(31). pp. 4-17. (In Russian).
52. Bakholdina, M. & Gortsev, A. (2014) Joint probability density of the intervals length of the modulated semi-synchronous integrated flow of events and its recurrence conditions. Communications in Computer and Information Science. 487. pp. 18-25. DOI: 10.1007/978-3-319-13671-4_3
53. Bakholdina, M.A. & Gortsev, A.M. (2014) [Joint probability density of the intervals length of the modulated semi-synchronous integrated flow of events and its recurrence conditions]. Informatsionnye tekhnologii i matematicheskoe modelirovanie (ITMM-2014) [Information Technologies and Mathematical Modeling (ITMM-2014)]. Proc. of the 13th International Scientific Research and Practice Conference named after A.F. Terpugov. November 20-22, 2014. Tomsk: Tomsk State University. pp. 137-143. (In Russian).
54. Bakholdina, M.A. & Gortsev, A.M. (2015) [Probability density of the interval length between consecutive events of the modulated semi-synchronous integrated flow of events in condition of a constant dead time]. Teoriya veroyatnostey, sluchaynye protsessy, ma-tematicheskaya statistika i prilozheniya [Probability Theory, Random Processes, Mathematical Statistics and Applications]. Proc. of International Scientific Conference dedicated to 80th anniversary of Professor, Doctor of Physical and Mathematical Sciences G.A. Medvedev. Minsk. Feburary 23-26, 2015. Minsk: RIVSh. pp. 17-22. (In Russian)
55. Bakholdina, M., Gortsev, A. (2015) Joint probability density of the intervals length of modulated semi-synchronous integrated flow of events in conditions of a constant dead time and the flow recurrence conditions. Communications in Computer and Information Science. 564. pp. 13-27. DOI: 10.1007/978-3-319-13671-4_3
