Файл: Казаков В.А. Введение в теорию марковских процессов и некоторые радиотехнические задачи.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 28.06.2024

Просмотров: 205

Скачиваний: 3

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

Матрицы

такого

рода

в литературе называют с т о х а ­

с т и ч е с к и м и .

 

 

 

 

 

 

В том случае, когда вероятности перехода не зависят

от времени,

цепь Маркова называется о д н о р о д н о й и

матрица

перехода

записывается

в более

простом

виде

 

 

 

 

 

 

 

 

(1.2)

Рассмотрим

несколько примеров.

 

 

Пример

1.

Пусть

система

имеет

три .состояния

Хь х% ха\

процесс

в системе развивается

в соответствии

со схемой однородной цепи Маркова; возможные

пере­

ходы системы из состояния ів состояние задаются

матри­

цей перехода

 

 

 

 

 

 

- По виду матрицы можно заключить, что," если система

находится в состоянии х±, то через один шаг во

времени

она с равной вероятностью может остаться в

том же

состоянии или перейти в любое из

двух других. Вторая

строка

матрицы показывает, что система с малой вероят­

ностью

(Ѵб) останется в состоянии

х% и может перейти

в состояния xi и хз с вероятностями

Ѵг и Ѵз соответствен­

но. Если же система находится в состоянии х 3 ,

то наибо­

лее вероятно, что в следующий момент она

останется

вэтом состоянии—вероятность этого события равна 3 Д,

водном случае из четырех она перейдет в состояние хі;

переход же из состояния Хз в состояние х% невозможен. Пример 2. Обратимся к описанию движения блуждаю­ щей, под действием случайных сил частицы. Рассмотре­ ние этого движения позволяет проиллюстрировать основ­ ные разновидности марковских цепей, что определило широкое распространение этого примера в литературе

[1 - 4] .

Пусть движение частицы ограничено двумя экранами, установленными на границах. Зададим число возможных состояний N=5, совместив граничные точки с состоя­ ниями ХІ и л*. Предположим, что экраны являются п о-

9



г л о щ а ю щ ІІ м и . Это означает, что процесс, достигнув состояний ХІ, х5, не может перейти из них ни в какое другое состояние. Следовательно, рц—\ для /=1,5. Что касается внутренних состояний (1 = 2, 3, 4), то по усло­ вию в каждый момент времени частица должна переме­

щаться

.по направлению

к состоянию

Хі

или

«

состоя­

нию х5

с вероятностями

q и р

соответственно. Очевидно,

/? + <7=1, таік что /011 =

0,

если

і=2, 3, 4. Матрица

перехо­

да для

этого случая

 

имеет вид

 

 

 

 

 

 

PIT.

 

( 1

0

0

0

0

 

 

0

р

0

0

 

Ргг

 

ч

 

 

 

 

 

 

0

ч

0

р

0

 

 

 

 

 

 

0

0

ч

0

р

 

Ріі

 

V о

0

0

0

1

 

 

 

 

 

 

 

 

 

 

(1.3)

Пример 3. Запишем матрицу перехода для блуждаю­ щей частицы между о т р а ж а ю щ и м и экранами.. Условия отражения могут быть весьма разнообразными. Сначала мы предположим, что частица при достижении граничной точки на следующем шаге отражается в то состояние, из которого она пришла. Для данной задачи такими состояниями могут быть только х-ь и х і и поэтому матрица перехода выглядит следующим образом:

 

0

1

0

0

0

 

 

 

q

0

р

0

0

 

 

Р =

0

яq

0

р

0

(1.4)

0

0

 

0

0

ч

0

р

 

 

 

0

0

0

1

0

 

 

Условия отражения могут быть такими, что частица из граничной точки переходит в любое другое заданное

состояние, например в х3.

В этом

счучае

 

 

 

9

0

1

0

0

 

 

 

ч

0

р

0

0

 

 

р =

0

чg

0

р

0

(1.5)

0

0

 

0

0

ч

0

р

 

 

 

0

0

1

0

0

 

 

Матрицы (1.4), (1.5) описывают эффект полного от­ ражения от границы. Возможны и другие варианты. Рас­ смотрим пример полупоглощающих (или полуотражаю­ щих) границ, для которых характерно, что частица на

10


следующем шаге с вероятностью и отражается в какоето состояние (пусь в л'з) и с вероятностью 1—и остается на границе. Тогда

— и

0

а

0

0

 

Q

0

Р

0

0

(1.6)

0

q

' 0

Р

0

0

0

я

0

р

 

0

0

и

0

1 —

 

Аналогично можно составить матрицы перехода для дру­ гих граничных условий. Матрицы перехода удобно пред­ ставлять в виде сигнальных графов, интерпретируя пере­ ходы системы из состояния ХІ В состояние Xj кяк переда-

р г г і

Р32=Ч

Д«=£

PSM=f

 

Рис.

!.2.

 

чу сигнала от одного узла к другому с коэффициентом передачи ветви рц. Сигнальные графы для матриц (1.3) — (1.6) изображены на рис. 1.1—1.4 соответственно.

Вероятности перехода — важнейшие

характеристики

любой марковской цепи, однако они по

определению

являются у с л о в н ы м и , и поэтому знание матрицы пе­

рехода не полностью определяет цепь Маркова. Если отнести матрицу перехода к первому шагу, определяю­ щему н а ч а л о работы системы, то для исключения условности необходимо задать еще вероятности началь­ ных состояний — начальные условия.

Вероятности начальных состояний р[°\ р%\ ..., р^" являются безусловными вероятностями и образуют ма­ трицу-строку Ро = (/?{0), р{°\ ... , р(®), сумма элементов

11


которой по условию нормировки должна быть равна еди­ нице. Характер распределения вероятностей начальных состояний целиком определяется условиями задачи. В ча­ стности, в начальный момент система может находиться

Рис. 1.4.

в каждом из N С О С Т О Я Н И Й С равной вероятностью. При этом

pi°) = llN. для всех і.

В том случае, когда начальное состояние системы извест­ но точно, то

/?<0 ) =0 для всех іф\ и р(°) = 1.

1.3. Уравнение Колмогорова—Чепмена. Классификация состояния и цепей

Матрица перехода дает исчерпывающее представле­ ние о вероятностях возможных переходов за о д и н шаг. Естественно возникает вопрос: как рассчитать вероятно­ сти того, что система, находящаяся в данный момент в состоянии Xï, перейдет в состояние х,- за 2, 3, ... , п ша­ гов. Иными словами, требуется найти матрицы перехода для заранее заданного числа шагов, если известна ма­ трица перехода за один шаг *).

*' Эта и все последующие зада™ и данной главе непосредствен­ но связаны с матричными операциями. Изложение теории матриц можно .найти в курсах высйей алгебры (см., например, [5, 6]).

12

Для однородной цепи Маркова вероятность перехода системы из г'-го состояния в j - e за два шага pij{2) опре­ деляется следующим очевидным соотношением, которое учитывает все возможные пути перехода

 

N

 

 

 

 

 

 

 

 

 

 

PU ( 2 ) = S

PuPij,

і,

У =

1, .2,..., ff,

 

(1.7)

 

f=i

 

 

 

 

 

 

 

 

 

 

где pu, pij — элементы

заданной

 

матрицы

перехода за

один шаг. Совокупность вероятностей перехода

за

два

шага pij(2) составляет

'матрицу

перехода

за два

шага

Р(2):

 

 

 

 

 

 

 

 

 

 

 

 

 

/Рп

(2)

Р» (2)

 

. . .

р

ш

(2)

 

 

р (2) _

Р п

( 2 )

Р а ( 2 )

• •

р

(

2

)

 

 

 

W i ( 2 )

РтЖ

 

• •

 

PNNW

 

 

Соотношение

(1.7)

позволяет

заключить,

 

что

матрица

Р(2) является произведением двух одинаковых матриц

перехода (1.2),

т. е. Р(2) = Р - Р =

Р2 .

 

 

Аналогично

вероятность

перехода системы из і-го

в /-е состояние за три шага

p,j(3)

можно

вычислить

по

формуле

 

 

 

 

 

/>«(3) = 2 PÜ(2)PI5

=

T,PÜPIJ

(2).

(1.8)

 

i-i

г=і

 

 

Это означает, что матрица перехода за три шага Р(3) равна 'произведению матриц перехода за один и два шага: Р(3) = Р(2) Р = Р • Р(2) = Р3 .

Ясно, что матрица перехода за п шагов Р(/г) вычи­ сляется как /г-я степень матрицы перехода за один шаг

Р ( л ) = Р " .

(1.9)

Вычисление вероятностей перехода

за 2, 3,*..., п ша­

гов производится путем суммирования произведений ве­ роятностей перехода, относящихся к начальному, про­ межуточному и конечному моментам времени. Если для определенности предположить, что процесс в системе

наблюдается с нулевого момента, то

при

вычислении

вероятностей перехода за

два

шага

промежуточным

является дискретный

момент времени

h. Действительно,

в

соотношении

(1.7)

фигурируют

вероятности перехода

за

один шаг: из

нулевого

момента в первый

и затем из

13