CC BY
34
ВЕСТНИК ТОМСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА
2010 Управление, вычислительная техника и информатика № 3(12)
УДК 519.21
А.М. Горцев, М.А. Леонова
ОПТИМАЛЬНАЯ ОЦЕНКА СОСТОЯНИЙ ОБОБЩЕННОГО АСИНХРОННОГО ПОТОКА В УСЛОВИЯХ НЕПРОДЛЕВАЮЩЕГОСЯ МЕРТВОГО ВРЕМЕНИ
Рассматривается задача оптимальной оценки состояний асинхронного дважды стохастического потока с инициированием дополнительных событий (обобщенный асинхронный поток событий) с двумя состояниями, являющегося математической моделью информационных потоков событий, функционирующих в цифровых сетях интегрального обслуживания - ЦСИО (Integrated Service Digital Networks - ISDN). Условия наблюдения за потоком таковы, что каждое событие порождает период мертвого времени, в течение которого другие события потока недоступны наблюдению и не вызывают продления его периода (непродлевающееся мертвое время). Находится явный вид апостериорных вероятностей состояний потока. Решение о состоянии потока выносится по критерию максимума апостериорной вероятности. Приведены численные результаты, полученные с
использованием расчетных формул и имитационного моделирования.
Ключевые слова: обобщенный асинхронный поток событий, состояние потока, апостериорная вероятность состояния, оценка состояния, мертвое время.
В последние два десятилетия в связи с бурным развитием компьютерной техники и информационных технологий появилась важная сфера приложений теории массового обслуживания - проектирование и создание информационно-вычислительных сетей, компьютерных сетей связи, спутниковых сетей связи, телекоммуникационных сетей и т.п., которые можно объединить единым термином - ЦСИО. Вполне естественно, что это дало толчок к необходимости построения новых математических моделей входящих потоков событий, достаточно адекватно описывающих реальные информационные потоки, функционирующие в ЦСИО. Отметим, что одними из первых работ в этом направлении были статьи [1 - 3]. Подчеркнем, однако, что в литературе по теории массового обслуживания и ее приложениям в целом, как и в литературе, посвященной исследованию ЦСИО, в частности, довольно незначительное количество работ посвящено адаптивным системам обслуживания, т.е. системам, функционирующим в условиях полной или частичной неопределенности. Более того, подавляющее число авторов рассматривают ситуации, когда все параметры, характеризующие входящий поток событий, априорно известны, хотя в реальных ситуациях дело обстоит, как правило, иначе. На практике параметры, определяющие входящий поток событий, известны либо частично, либо вообще неизвестны, либо (что ещё больше ухудшает ситуацию) они изменяются со временем, при этом изменения часто носят случайный характер, что приводит к рассмотрению дважды стохастических потоков событий. С другой стороны, очевидно, что функционирование системы обслуживания непосредственно зависит от параметров входящего потока событий. Потоки событий с интенсивностью, зависящей от времени и являющейся случайным процессом
(дважды стохастические потоки событий), можно разделить на два класса. К первому классу относятся потоки с интенсивностью, являющейся непрерывным случайным процессом. Ко второму классу относятся потоки, у которых интенсивность есть кусочно-постоянный случайный процесс с конечным числом состояний. Последние (потоки с переключениями или МС-потоки событий [1, 2]) являются наиболее характерными для реальных телекоммуникационных сетей. В свою очередь, в зависимости от того, каким образом происходит переход из состояния в состояние, МС-потоки можно разделить на три типа: 1) синхронные дважды стохастические потоки событий - потоки с интенсивностью, для которой переход из состояния в состояние происходит в случайные моменты времени, являющиеся моментами наступления событий [4 - 6]; 2) асинхронные дважды стохастические потоки событий - потоки с интенсивностью, для которой переход из состояния в состояние происходит в случайные моменты времени и не зависит от моментов наступления событий [7, 8]; 3) полусинхронные дважды стохастические потоки событий - потоки, у которых для одного множества состояний справедливо 1, а для остальных состояний справедливо 2 [9 - 11]. Наконец, подчеркнем, что отмеченные синхронные, асинхронные и полусинхронные потоки событий возможно представить в виде моделей МАР (Markovian Arrival Ргосе$5)-потоков событий [12, 13] с определенными ограничениями на параметры последних [14]. При исследовании дважды стохастических потоков событий могут быть выделены два класса задач: 1) задача фильтрации интенсивности потока (или задача оценивания состояний потока событий) по наблюдениям за потоком (по наблюдениям за моментами наступления событий) [7, 15]; 2) задача оценивания параметров потока по наблюдениями за моментами наступления событий [4 - 6, 8 - 11].
Одними из первых работ по оценке состояний дважды стохастических потоков, по-видимому, являются [7,16,18], в которых рассматриваются асинхронные дважды стохастические потоки с двумя и произвольным числом состояний, по оценке параметров - работа [17].
В большинстве публикаций авторы рассматривают математические модели потоков событий, когда все события доступны наблюдению. Однако на практике возникают ситуации, когда наступившее событие может повлечь за собой нена-блюдаемость последующих событий. Одним из искажающих факторов при оценке состояний и параметров потока событий выступает мертвое время регистрирующих приборов [19], которое порождается зарегистрированным событием. Другие же события, наступившие в течение периода мертвого времени, недоступны наблюдению (теряются). Можно считать, что этот период продолжается некоторое фиксированное время Т. Все устройства регистрации делятся на две группы. Первую группу составляют устройства с непродлевающимся мертвым временем, вторую - устройства с продлевающимся мертвым временем. Отметим, что одними из первых по оценке параметров случайных потоков событий, функционирующих в условиях мертвого времени, являются работы [20 - 22].
В работах [23, 24] решена задача об оптимальной оценке состояний асинхронного дважды стохастического потока с инициированием дополнительных событий (обобщенный асинхронный поток) в условиях отсутствия мертвого времени. В настоящей статье, являющейся непосредственным развитием работы [24], решается задача об оптимальной оценке состояний обобщенного асинхронного потока событий в условиях его неполной наблюдаемости (при непродлевающемся мертвом времени). Предлагается алгоритм оптимальной оценки состояний, когда ре-
шение о состоянии обобщенного асинхронного потока выносится по критерию максимума апостериорной вероятности, представляющей наиболее полную характеристику состояния потока, которую можно получить, располагая только выборкой наблюдений, и обеспечивающей минимум полной (безусловной) вероятности ошибки вынесения решения [25].
1. Постановка задачи
Рассматривается асинхронный дважды стохастический поток с инициированием дополнительных событий (далее обобщенный асинхронный поток), интенсивность которого есть кусочно-постоянный случайный процесс Х(г) с двумя состояниями Х1 и Х2 (^1 > Х2). В течение временного интервала, когда Х(г) = X, , имеет место пуассоновский поток событий с интенсивностью X, , ,=1,2. Переход из первого состояния процесса Х(г) во второе (из второго в первое) может осуществляться в произвольный момент времени. При этом длительность пребывания процесса Х(г) в ,-м состоянии распределена по экспоненцильному закону с параметром а, , , = 1,2. При переходе процесса Х(*) из первого состояния во второе инициируется с вероятностью р (0<р <1) дополнительное событие во втором состоянии (т.е. сначала осуществляется переход, а затем инициируется дополнительное событие). Наоборот, при переходе процесса Х(г) из второго состояния в первое инициируется с вероятностью я (0< я <1) дополнительное событие в первом состоянии. Очевидно, что в сделанных предпосылках Х(г) - марковский процесс. После каждого зарегистрированного в момент времени г, события наступает время фиксированной длительности Т (мертвое время), в течение которого другие события исходного потока недоступны наблюдению. События, наступившие в течение мертвого времени, не вызывают продления его периода (непродлевающееся мертвое время). По окончании мертвого времени первое наступившее событие снова создает период мертвого времени длительности Т и т. д. Вариант возникающей ситуации показан на рис.1, где 1, 2 - состояния случайного процесса Х(г); дополнительные события, которые могут наступать в момент перехода процесса Х(г) из состояния в состояние, помечены буквами р либо q; штриховка - периоды мертвого времени длительности Т; г1, г2, ... - моменты наступления событий в наблюдаемом потоке.
1 ; 1а 1 а2, 2 : Т ,,а1 а2| 1а1 а2,, „а1
[
: \р : ; Процесс Х(г) : я \ \р г
; ; ; Обобщенный асинхронный поток |
; Т Т ■ Т Т \ Схема создания непродлевающегося мёртвого времени —ё ё е ё ► *1 *2 *3 *4 *
Наблюдаемый поток событий
Рис. 1. Формирование наблюдаемого потока событий
Подчеркнем, что если р = я = 0, то имеет место обычный асинхронный поток [7, 8]. Так как процесс Х(г) и типы событий (события пуассоновских потоков и дополнительные события) являются ненаблюдаемыми, а наблюдаемыми являются только временные моменты наступления событий наблюдаемого потока г1, г2, ..., то необходимо по этим наблюдениям оценить состояние процесса (потока) Х(г) в момент окончания наблюдений.
Рассматривается установившийся (стационарный) режим функционирования наблюдаемого потока событий, поэтому переходными процессами на интервале наблюдения (г0, г), где г0 - начало наблюдений, г - окончание наблюдений (момент вынесения решения), пренебрегаем. Тогда без потери общности можно положить г0 = 0. Для вынесения решения о состоянии ненаблюдаемого процесса Х(г) в момент времени г необходимо определить апостериорные вероятности ю(Х,| г) = = ю(Х,| гь..., гт , г), , = 1, 2, того, что в момент времени г значение процесса Х(г) = Х, (т - количество наблюденных событий за время г), при этом ю(Х1| г) + ю(Х2| г) = 1. Решение о состоянии процесса Х(1) выносится путем сравнения апостериорных вероятностей: если ю(Х,| г) > ю(Х,| г) , ,,] = 1, 2, , Ф], то оценка состояния процесса есть Х(г) = X.
2. Апостериорная вероятность состояния наблюдаемого потока событий
Рассмотрим интервал времени (г,, г,+1) между двумя соседними событиями в наблюдаемом потоке, значение длительности которого есть т, = г,+1 - г,,, = 1, 2, ... (длительность интервала - случайная величина). С другой стороны, значение длительности этого интервала есть т, = Т + п, , где п, - значение длительности интервала между моментом окончания периода мертвого времени и моментом г, +1, т.е. интервал (г,, г, +1) разбивается на два смежных интервала: первый - (г,, г, + Т), второй - (г, + Т, г, +1). Подчеркнем одно важное обстоятельство: так как первое наступившее после окончания периода мертвого времени событие обобщенного асинхронного потока снова порождает период мертвого времени фиксированной длительности Т, в течение которого последующие события обобщенного асинхронного потока недоступны наблюдению (поток отсутствует), то условия нахождения апостериорной вероятности ю(Х1| г) на интервале (г,, г, + Т) длительности Т и интервале (г, + Т, г, +1), значение длительности которого есть п,, принципиально разные. Кроме того, для нахождения апостериорной вероятности ю(Х1| г) необходимо точно знать длительность Т мертвого времени либо, по крайней мере, предварительно осуществить ее оценку Т . В противном случае отсутствие информации о значении длительности Т мертвого времени делает попытку строгого нахождения апостериорной вероятности ю(Х1| г) невозможной. Здесь предполагается, что значение Т известно точно.
В [23] сформулирован алгоритм расчета апостериорной вероятности ю(Х1| г) для случая отсутствия мертвого времени (Т=0). При этом поведение вероятности ю(Х1| г) на полуинтервале [г,, г,+1) между соседними наблюдавшимися событиями обобщенного асинхронного потока определяется выражением
| г) Ю! [ю2 -Ю(Х! | г+0)]-Ю2 [<Й1 -Ю(Х! | г, + 0)] е~в(г ) (Х)
1 ю2 -ю(Х1 | гi + 0) -[ -ю(Х1 | гi + 0)]в(г-г)
где гi < г < гi+1,, = 0,1...;р = (ю2 -ю1)(Х1 -X2 + ра1 - яа2)Ф 0 ,
[ - Х2 + а1 + (1 - 2д) а2 ] - - Х2 + а1 - а2 )2 + 4а1а2 (1 - р(1 - д)
2 ( -X2 + ра1 - да2 )
2 (Х1 -X2 + ра1 - да2 )
В момент времени ґ (в момент наступления события обобщенного асинхронного потока) апостериорная вероятность (1) претерпевает разрыв 1-го рода (і = 1,2,...), поэтому в момент времени ґі имеет место формула пересчета
где ю(Х1| г, -0) вычисляется по формуле (1) в момент времени г = гi , когда г изменяется в полуинтервале [г,-1, г,), соседнем с полуинтервалом [г,, г,+1). В качестве начального значения ю(Х1| г0+0) = ю(Х1| г0=0) в (1) выбирается априорная финальная вероятность первого состояния процесса Х(г): п1=а2 /( а1+ а2), которая находится из уравнений а1п1-а2п2 = 0, п1+п2 = 1.
Таким образом, вычисление апостериорной вероятности ю(Х1| г) по формуле (1) в условиях, когда длительность мертвого времени Т Ф 0, справедливо на полуинтервале [г,+Т, гт), значение длительности которого есть п,. При этом начальное условие для ю(Х1| г) привязывается к моменту времени г,+Т, т.е. в формуле (1), во-первых, нужно ю(Х1| г, +0) заменить на ю(Х1| г, +Т), во-вторых, г,+Т <г < г,+1, , = 1,2,. . Формула пересчета (2) остается при этом без изменения, так как она предназначена для вычисления апостериорной вероятности в момент г, наступления наблюдаемого события, которое порождает мертвое время.
Рассмотрим отрезок [г,, г,+Т ], ,= 1,2,. . Так как на этом отрезке длительности Т событие наблюдаемого потока имеет место в граничной точке г, , а на самом отрезке события отсутствуют, то необходимо определить поведение апостериорной вероятности ю(Х1| г) на отрезке [г,, г,+Т ]. Определим апостериорную вероятность ю(Х11 г +Лг) того, что в момент времени г +Аг (г, < г+Лг < г,+Т), где дг - достаточно малая величина, процесс Х(г) находится в первом состоянии. Пусть в момент времени г процесс Х(г) находится в первом состоянии и на интервале (г, г +Аг) процесс Х(г) не перешел во второе состояние. Вероятность этого события есть ю(Х1| г)(1-а1Лг) + о(Лг). Пусть в момент времени г процесс Х(г) находится во втором состоянии и на интервале (г, г +Аг) процесс Х(г) перешел в первое состояние. Вероятность этого события есть ю(Х2| г) а2Лг +о(Лг). Другие возможности имеют вероятность о(Лг). Тогда
Производя в (3), (4) необходимые преобразования, после чего переходя к пределу при Л1—*0, получаем систему дифференциальных уравнений для апостериорных вероятностей ю(Х1| г) и ю(Х2| г):
с граничными условиями: ю(Х1| ґ=ґ)= ю(Х1| ґі +0), ю(Х2| ґ=ґ)= ю(Х2| ґі +0), і = 1, 2,. .
, і = 1,2,..., (2)
Х2 + да 2 + (Х1 -X 2 + ра1 - да 2) ю(Х1 | ґі - 0)
ю(Х1 | ґ + Дґ) = (1 -а1Дґ)ю(Х1 | ґ) + а2Дґю(Х2 | ґ) + о(Дґ). Аналогично находим ю(Х2| ґ +Аґ):
ю(Х2 | ґ + Дґ) = (1 -а2Дґ)ю(Х2 | ґ) + а1Дґю(Х1 | ґ) + о(Дґ).
(3)
(4)
сіґ 1 1 2 2 Ж
Последнее вытекает из того, что на полуинтервале [г,-1+Т, г,), , = 2,3,., смежном отрезку [г,, г,+Т ], апостериорная вероятность рассчитывается по формуле (1), где вместо ю(Х1| г,+0) стоит ю(Х1| г,+Т); в точке г=г, происходит пересчет апостериорной вероятности по формуле (2), так что её значение в этой точке есть ю(Х11 г,+0). Для граничного полуинтервала [г0, г1) расчет апостериорной вероятности ю(Х1| г) производится по формуле (1) с ее последующим пересчетом по формуле (2) в точке г = г1. Решая систему (5), находим
ю(Х 1| г) = п1 +[ю(Х 1| гг + 0)-п1 ](а1+“2)(г-г,-), (6)
где г, < г < г,+Т, ,=1,2,.; вероятность п1 определена в (1). Тогда из (6) следует, что ю(Х 1| гг + Т) = п1 +[со(Х 1| г, + 0)-п1 ]еч<Х1+“2)Т , , = 1,2,.... (7)
Полученные формулы позволяют сформулировать алгоритм расчета апостериорной вероятности ю(Х1| г) и алгоритм принятия решения о состоянии процесса Х(г) в любой момент времени г:
1) в момент времени г0=0 задается ю(Х1|г0+0)=ю(Х1|г0=0)=л1;
2) по формуле (1) для ,=0 рассчитывается вероятность ю(Х1|г) в любой момент времени г (0 < г < г1), где г1 - момент наблюдения первого события наблюдаемого потока;
3) по формуле (1) для ,=0 рассчитывается вероятность ю(Х1|г) в момент времени г1: ю(Х1|г1)=ю(Х1|г1-0);
4) увеличивается на единицу и по формуле (2) для =1 производится пересчет апостериорной вероятности ю(Х1|г) в момент времени г= г1, при этом ю(Х1|г1+0) является начальным значением для ю(Х1|г) в формуле (6);
5) по формуле (6) для ,=1 рассчитывается вероятность ю(Х1|г) в любой момент времени г (г1 < г < г1+Т);
6) по формуле (7) для , =1 рассчитывается вероятность ю(Х1|г) в момент времени г= г1+Т, т.е. ю(Х1|г1+Т); при этом ю(Х1|г1+7) является начальным условием для ю(Х1|г) на следующем шаге алгоритма;
7) для , =1 по формуле
| г) Ю1 [аз2 -ю(Х1 | г, +Т)]-ю2 [Ю1 -ю(Х1 | г, +Т)]е~т--Т) (8)
1 ю2 -ю(Х11 г, + т) -[ -ю(Х11 гг + т)]е-Р(г-г'-Т)
(г1+Т < г < г+1, ,= 1,2,. ; ю1, ю2, в определены в (1)) рассчитывается вероятность ю(Х1|г) в любой момент времени г (г1+Т< г < г2), где г2 - момент наблюдения второго события наблюдаемого потока;
8) по формуле (8) для ,=1 рассчитывается вероятность ю(Х1|г) в момент времени г = г2: ю(Х1|г2) = ю(Х1|г2-0);
9) алгоритм переходит на шаг 4, после чего шаги 4 - 8 повторяются для , = 2 и т.д.
Параллельно по ходу вычисления апостериорной вероятности ю(Х1|г) в любой момент времени г выносится решение о состоянии процесса Х(г): если
ю(Х1|г) >ю(Х2|г) ( ю(Х1|г)>1/2 ), то оценка к(г) = Х1, в противном случае Х(г) = Х2.
Для частных и особых случаев, приведенных в [24] и возникающих при определенном соотношении параметров Х1 , Х2 , а1 , а2 , р, д, алгоритмы оптимальной оценки состояний обобщенного асинхронного потока идентичны приведенному алгоритму для общего случая (изменяются только формулы (1), (2), (6), (7)). Вследствие этого алгоритмы для частных и особых случаев здесь отдельно не рассматриваются.
3. Результаты численных расчетов
Для получения численных результатов разработан алгоритм вычисления апостериорной вероятности ra(Xi|0 по формулам (2), (6), (7), (8). Программа расчета реализована на языке программирования C++ в среде Builder 6. Первый этап расчета предполагает имитационное моделирование обобщенного асинхронного потока событий, схемы создания непродлевающегося мертвого времени и наблюдаемого потока событий. Описание алгоритма имитационного моделирования здесь не приводится, так как никаких принципиальных трудностей алгоритм не содержит. Второй этап расчета - непосредственное вычисление вероятностей ю(Х1| t), t0 < t < t1; ra(X1|t,+0); ra(X1|t), ti < t < t+T; ra(X1|t,+T); ra(X1|t), t+T < t < ti+1, i = 1, 2, ..., по формулам (1), (2), (6) - (8) соответственно и определение оценки X(t). Расчеты произведены для следующих значений параметров: Х1 = 4, Х2 = 1, а1 = 0,05, а2 = 0,08, p = 0,4, q = 0,7, Т = 2 ед. времени и времени моделирования Tm = 20 ед. времени. В качестве иллюстрации на рис. 2 приведена траектория (нижняя часть рис. 2) случайного процесса X(t) (истинная траектория), полученная путем имитационного моделирования, где 1, 2 - состояния процесса X(t), и траектория (верхняя часть рис. 2) оценки X(t), где 1, 2 - состояния оценки X(t). Вынесение решения о том или ином состоянии процесса X(t) производилось с шагом At = 0,05. На рис. 2 штриховкой на оси времени обозначены временные промежутки, на которых оценка состояния не совпадает с истинным значением процесса X(t) (области ошибочных решений). На рис. 3 приведена траектория поведения апостериорной вероятности ra(X1|t), соответствующая полученной при имитационном моделировании последовательности моментов наступления событий t1, t2,. и последовательности моментов времени окончания периодов мертвого времени t1+T, t2+T,....
4t)
1
2
i(0
1 2
0
ю(Х1|г)
1
0,5
4 8 12 16
Рис. 2. Траектории процесса X(t) и оценки k(t)
t
t
0 4 8 12 16 t
Рис. 3. Траектория апостериорной вероятности ro(X1|t)
Для установления частоты ошибочных решений о состоянии случайного процесса Х(0 по наблюдениям за обобщенным асинхронным потоком, функционирующим в условиях мертвого времени (наблюдаемый поток), проведен статистический эксперимент, состоящий из следующих этапов: 1) для определенного набора параметров Хь Х2, а1, а2, р, д, Т ед.времени осуществляется моделирование наблюдаемого потока событий на заданном отрезке времени [0, Тт] (отдельный 1-й эксперимент, I = 1, 2,. ); 2) осуществляется расчет апостериорных вероятностей ю(Х1|/) первого состояния процесса Х(0 на заданном отрезке времени [0, Тт] по формулам (1), (2), (6) - (8); 3) осуществляется оценивание траектории процесса Х(/) (оценивание на отрезке [0, Тт] интервалов, когда оценка Х(/) принимает то
или иное значение); 4) осуществляется определение (для отдельного г-го эксперимента) - суммарной протяженности интервалов, на которых значение процесса
Х(/) не совпадает с его оценкой Х(/); 5) осуществляется вычисление доли ошибочных решений рг = й^Т; 6) осуществляется повторение N раз (г = 1, N) шагов 1 - 5 для расчета оценки безусловной вероятности ошибки оценивания состояний процесса Х(0 на отрезке [0, Тт].
Результатом выполнения описанного алгоритма является выборка Р1, р2,...,PN долей ошибочных решений в N экспериментах. По этому набору вычисляются выборочное среднее безусловной вероятности ошибочного решения
1 N 1 N
Р0 = ^ XРг и выборочная дисперсия I = ^^ ХР - Рг )2.
Результаты статистического эксперимента приведены в табл. 1 - 5. В первой строке таблиц указана длительность мертвого времени Т (Т = 1, 2, ..., 9 ед. времени). Во второй и третьей строках таблиц для каждой длительности мертвого времени Т приведены численные значения для Р0 и I соответственно.
Результаты получены при следующих значениях параметров, общих для всех таблиц: Х2 = 1, а! = 0,1, а2 = 0,15, р = 0,4, д = 0,7, Тт = 100, N = 100. При этом результаты в табл. 1 получены для ^ = 4, в табл. 2 - для ^ = 5, в табл. 3 - для ^ = 6, в табл. 4 - для = 7, в табл. 5 - для ^ = 8.
Т аблица 1
Результаты статистического эксперимента (11 = 4)
Т 1 2 3 4 5 6 7 8 9
р 0,2482 0,2797 0,2945 0,2988 0,3062 0,3315 0,3436 0,3441 0,3451
5 0,0048 0,0074 0,0079 0,0097 0,0067 0,0118 0,0123 0,0136 0,0150
Т аблица 2
Результаты статистического эксперимента (>^ = 5)
Т 1 2 3 4 5 6 7 8 9
Р 0,2228 0,2631 0,2731 0,3009 0,3190 0,3304 0,3316 0,3339 0,3403
5 0,0039 0,0054 0,0129 0,0097 0,0109 0,0115 0,0166 0,0148 0,0193
Т аблица 3
Результаты статистического эксперимента (11 = 6)
т 1 2 3 4 5 6 7 8 9
Р0 0,2134 0,2487 0,2625 0,2807 0,2955 0,3071 0,3063 0,3187 0,3201
Г 0,0044 0,0051 0,0067 0,0120 0,0091 0,0101 0,0124 0,0080 0,0094
Т аблица 4
Результаты статистического эксперимента (>^ = 7)
т 1 2 3 4 5 6 7 8 9
Р0 0,2003 0,2251 0,2551 0,2572 0,2841 0,3011 0,3050 0,3113 0,3189
Г 0,0024 0,0080 0,0080 0,0063 0,0079 0,0081 0,0101 0,0111 0,0109
Т аблица 1
Результаты статистического эксперимента (11 = 8)
т 1 2 3 4 5 6 7 8 9
Р0 0,1920 0,2220 0,2496 0,2531 0,2749 0,2933 0,3006 0,3117 0,3167
Г 0,0027 0,0043 0,0066 0,0061 0,0090 0,0119 0,0061 0,0063 0,0106
Во-первых, отметим, что анализ проведенных многочисленных вариантов численных расчетов по оценке безусловной вероятности ошибочного решения Р0 показывает, что оценка Р0 является достаточно стабильной для Тт > 100 ед.времени. Вследствие этого время моделирования Тт для всех экспериментов , результаты которых представлены в табл. 1 - 5, было выбрано равным 100 ед.времени. Во-вторых, анализ численных результатов, приведенных в табл. 1 - 5, показывает:
1) значение оценки безусловной вероятности ошибочного решения Р0 увеличивается с увеличением длительности мертвого времени Т (Т = 1, 2, ..., 9 ед. времени); последнее является вполне естественным, так как при увеличении длительности мертвого времени происходит увеличение потерь полезной информации о потоке событий, что в конечном итоге отрицательно сказывается на качестве оценивания;
2) при фиксированной длительности мертвого времени Т значения оценки Р0 уменьшаются в зависимости от ^ (^ = 4, 5, 6, 7, 8), что является естественным, так как при увеличении разности ^ - Х2 условия различимости состояний потока улучшаются; 3) оценка дисперсии Г для всех вариантов расчета достаточно мала.
4. Заключение
Полученные результаты показывают возможность оценивания состояний обобщенного асинхронного потока событий в условиях непродлевающегося мертвого времени по результатам текущих наблюдений (в течение некоторого временного интервала) за потоком. Это, в свою очередь, позволяет изменять режимы функционирования системы массового обслуживания в зависимости от того или иного состояния обобщенного асинхронного потока событий (адаптироваться к изменяющейся интенсивности входящего потока событий). Выражения
апостериорных вероятностей для оценки состояний потока получены в явном виде, что позволяет производить вычисления без привлечения численных методов. Сам же алгоритм оценки состояний потока обеспечивает минимум полной вероятности ошибки вынесения решения.
Наконец, отметим, что рассмотренный обобщенный асинхронный поток событий охватывает ранее изученные модели потоков, вытекающие из него как частные случаи [7, 8, 16, 17, 26 - 30].
ЛИТЕРАТУРА
1. Башарин Г.П., Кокотушкин В.А., Наумов В.А. О методе эквивалентных замен расчета фрагментов сетей связи. Ч.1 // Изв. АН СССР. Техн. кибернетика. 1979. № 6. С. 92 - 99.
2. Башарин Г.П., Кокотушкин В.А., Наумов В.А. О методе эквивалентных замен расчета фрагментов сетей связи. Ч. 2 // Изв. АН СССР. Техн. кибернетика. 1980. № 1. С. 55 - 61.
3. NeutsM.F. A versatile Markov point process // J. Appl. Probab. 1979. V. 16. P. 764 - 779.
4. Горцев А.М., Нежельская Л.А. Оценка параметров синхронного альтернирующего пуас-соновского потока событий методом моментов // Радиотехника. 1995. № 7 - 8. С. 6 - 10.
5. Горцев А.М., Нежельская Л.А. Оценивание длительности «мертвого времени» и интенсивностей синхронного дважды стохастического потока событий // Радиотехника. 2004. № 10. С. 8 - 16.
6. Бушланов И.В., Горцев А.М., Нежельская Л.А. // Оценка параметров синхронного дважды стохастического потока событий // АиТ. 2008. № 9. С. 76 - 93.
7. Горцев А.М., Нежельская Л.А. Оптимальная нелинейная фильтрация марковского потока событий с переключениями // Техника средств связи. Сер. Системы связи. 1989. Вып. 7. С. 46 - 54.
8. Васильева Л.А, Горцев А.М. Оценивание длительности мертвого времени асинхронного дважды стохастического потока событий в условиях его неполной наблюдаемости // АиТ. 2003. № 12. С. 69 - 79.
9. Горцев А.М., Нежельская Л.А. Оценивание параметров полусинхронного дважды стохастического потока событий методом моментов // Вестник ТГУ. Приложение. 2002. № 1(1). С. 18 - 23.
10. Горцев А.М., Нежельская Л.А. Оценивание периода мертвого времени и параметров полусинхронного дважды стохастического потока событий // Измерительная техника. 2003. № 6. С. 7 - 13.
11. Горцев А.М., Нежельская Л.А. Полусинхронный дважды стохастический поток событий при продлевающемся мертвом времени // Вычислительные технологии. 2008. Т. 13. № 1. С. 31 - 34.
12. Lucantoni D.M., Neuts M.F. Some steady-state distributions for the MAP/SM/1 queue // Communications in Statistics Stochastic Models. 1994. V. 10. P. 575 - 598.
13. Machihara F.A. A MAP/SM/1 queue with service times depending on the arrival process // Symposium on Perfomance Models for Information Communication Networks: Proc. Conf., Tokyo. 1997. P. 180 - 191.
14. Василевская Т.П., Завгородняя М.Е., Шмырин И.С. О соотношении моделей МАР-потока событий и асинхронного, полусинхронного и синхронного дважды стохастических потоков событий // Вестник ТГУ. Приложение. 2004. № 9(II). С. 138 - 144.
15. Бушланов И.В., Горцев А.М. Оптимальная оценка состояний синхронного дважды стохастического потока событий // АиТ. 2004. № 9. С. 40 - 51.
16. Горцев А.М., Нежельская Л.А. Оптимизация параметров адаптера при наблюдениях за МС-потоком // Стохастические и детерминированные модели сложных систем: сб. статей. Новосибирск: Изд-во ВЦ СО АН СССР, 1988. С. 20 - 32.
17. Горцев А.М., Коротаева Н.И. Оценка параметров МС-потока событий методом моментов // Распределенные микропроцессорные управляющие системы и локальные вычислительные сети: материалы Всесоюзной науч.-технич. конф., Томск, июнь 1991 г. Томск: Изд-во Том. ун-та, 1991. С. 69 - 72.
18. Горцев А.М., Полетавкин Ю.М. Оценка состояния дважды стохастического потока с произвольным числом состояний // Математические методы исследования сетей связи и сетей ЭВМ: тез. докл. Шестой Белорусской зимней школы-семинара по теории массового обслуживания, Витебск, январь-февраль 1990. Минск: Изд-во Белорусского унта, 1990. С. 36 - 37.
19. Апанасович В.В., Коляда А.А., Чернявский А.Ф. Статистический анализ случайных потоков в физическом эксперименте. Минск: Изд-во «Университетское», 1988. 254 с.
20. Горцев А.М., Климов И.С. Оценка интенсивности пуассоновского потока событий в условиях частичной его ненаблюдаемости // Радиотехника. 1991. № 12. С. 3 - 7.
21. Горцев А.М., Климов И.С. Оценивание параметров знакопеременного пуассоновского потока событий // Радиотехника. 1994. № 8. С. 3 - 9.
22. Горцев А.М., Климов И.С. Оценивание периода ненаблюдаемости и интенсивности пуассоновского потока событий // Радиотехника. 1996. № 2. С. 8 - 11.
23. Горцев А.М., ЛеоноваМ.А., Нежельская Л.А. Оптимальная оценка состояний асинхронного дважды стохастического потока с инициированием лишних событий // Массовое обслуживание: потоки, системы, сети: материалы Междунар. науч. конф. «Современные математические методы анализа и оптимизации информационно-телекоммуникационных сетей», Минск, 26 - 29 янв. 2009 г. Вып. 20. Минск: РИВШ, 2009. С. 90 - 96.
24. Горцев А.М., Леонова М.А. Оптимальная оценка состояний обобщенного асинхронного дважды стохастического потока // Вестник Томского госуниверситета. Управление, вычислительная техника и информатика. 2010. № 1(10). С. 33 - 47.
25. Хазен Э.М. Методы оптимальных статистических решений и задачи оптимального управления. М.: Сов. радио, 1968. 256 с.
26. Васильева Л.А., Горцев А.М. Оценивание параметров дважды стохастического потока событий в условиях его частичной наблюдаемости // АиТ. 2002. № 3. С. 179 - 184.
27. Горцев А.М., Шмырин И.С. Оптимальная оценка состояний дважды стохастического потока событий при наличии ошибок в измерениях моментов времени // АиТ. 1999. № 1. С. 52 - 66.
28. Горцев А.М., Завгородняя М.Е. Оценка параметров альтернирующего потока событий при условии его частичной наблюдаемости // Оптика атмосферы и океана. 1997. Т. 10. № 3. С. 273 - 280.
29. Горцев А.М., Паршина М.Е. Оценивание параметров альтернирующего потока событий в условиях «мертвого времени» // Изв. вузов. Физика. 1999. № 4. С. 8 - 13.
30. Горцев А.М., Ниссенбаум О.В. Оценивание длительности мертвого времени и параметров асинхронного альтернирующего потока событий с инициированием лишнего события // Вестник ТГУ. 2004. № 284. С. 137 - 145.
Горцев Александр Михайлович
Леонова Мария Алексеевна
Томский государственный университет
E-mail: [email protected] Поступила в редакцию 3 марта 2010 г.
