Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

g5

.pdf
Скачиваний:
5
Добавлен:
25.03.2016
Размер:
387.94 Кб
Скачать

n

Pij 0, Pij =1, для любых i. j=1

Квадратные матрицы, для которых выполняются эти свойства, называют стохастическими. Зная матрицу переходных вероятностей, можно построить граф состояний с отмеченными на нем переходными вероятностями. Граф, на котором отмечены переходные вероятности, называется размеченным

графом.

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

шаге для дискретной марковской цепи.

Вероятности P(n)

P[ξ

n

= j /(ξ

0

= i)], i =1,2...r называют вероятностями

ij

 

 

 

перехода цепи Маркова за n шагов, соответствующие матрицы вида

 

 

P(n)

P(n)

P(n)

...P(n)

 

 

 

11

12

13

1r

P

(n)

P(n)

P(n)

P(n)

...P(n)

 

=

21

22

23

2r

 

 

 

...............................

 

 

P(n)

P(n)

P(n)

...P(n)

 

 

 

r1

r2

r3

rr

называются матрицей вероятностей перехода за n шагов. Имеет место следующая теорема.

Теорема. Матрица перехода за n шагов есть n-я степень матрицы перехода за один шаг, т.е.

P(n) = (P)n .

Доказательство. По определению P[ξn = j /(ξ0 = i)] вероятность перехода из состояния Si в состояние S j . По формуле полной вероятности

 

 

 

 

r

 

 

 

 

 

P[ξn = j /(ξ0 = i)]= P(ξ1 = k / ξ0 = i)P(ξn = j / ξ1 = k,ξ0 = i),

откуда с учетом равенства

 

k =1

 

 

 

 

 

 

 

 

 

 

 

 

P(ξ

n

= j /

ξ

= k,ξ

0

= i)= P(ξ

n

= j / ξ = k )= P(n1)

 

 

1

 

 

1

kj

получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

Pij(n)

r

 

 

 

 

 

 

 

= Pik Pkj(n1), n = 2,3..

(1)

 

 

 

 

 

 

k =1

 

 

 

В матричной записи доказанное соотношение имеет вид P(n) = P(P)(n1).

Откуда по индукции следует P(n) = (P)n , что и требовалось доказать.

Вектор a = {a1, a2,....ar }, где

ai = P(ξ0 = i), i =1,2..r

называется вектором

начальных вероятностей. Имеет место следующее утверждение.

Утверждение 1. Свойства однородных марковских цепей полностью

 

 

 

 

 

P

P

P

...P

 

 

 

 

 

 

 

11

12

13

 

1r

 

 

 

 

 

 

 

 

 

 

 

 

 

определяются вектором a = {a , a

 

,....a

 

}и матрицей

P

P

P

 

...P

 

2

r

P =

21

22

23

2r

.

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

...............................

 

 

 

 

 

 

P

P

P

...P

 

 

 

 

 

 

 

r1

r2

r3

 

rr

 

Доказательство. Пусть в начальный момент времени система находится в одном из состояний Si с вероятностью ai = P(ξ0 = i). Через n шагов система

будет

находиться в

одном

из

 

состояний Si

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

P(n) P[ξ

n

= j /(ξ

0

= i)]. Так как события

(ξ

0

= i) образуют полную группу, то

ij

 

 

 

 

 

 

 

 

 

 

 

 

 

по формуле полной вероятности

 

 

 

 

 

 

 

 

 

 

 

 

 

P(ξ

n

= j)= r

P(ξ

0

= i)P(ξ

n

=

j / ξ

0

= i)= r a P(n),

 

 

 

 

 

 

 

 

 

 

 

i ij

 

 

 

 

 

 

 

i=1

 

 

 

 

 

 

 

 

i=1

где вероятности Pij(n) вычисляются по рекуррентной формуле (1). Что и требовалось показать.

Утверждение 2. Для однородной цепи Маркова при любом m выполняется равенство

P[ξn+m = j /(ξm = i)]= P[ξn = j /(ξ0 = i)].

Доказательство.

Очевидно, что событие (ξn+m = j /(ξm = i)) равносильно

сумме событий следующих событий

(ξn+m = j /(ξm = i))=

r r

r

= ∑ ∑ .... [(ξm+1 = k1,ξm+2 = k2,ξm+3 = k3,...ξm+n = j)/(ξm = i)] k1 =1k2 =1 kn1

Откуда по теореме сложения вероятностей несовместных событий получим

P(ξn+m = j /(ξm = i))=

= r r .... r P[(ξm+1 = k1,ξm+2 = k2,ξm+3 = k3,...ξm+n = j)/(ξm = i)](1) k1 =1k2 =1 kn1

где по теореме умножения

P[(ξm+1 = k1,ξm+2 = k2,ξm+3 = k3,...ξm+n = j)/(ξm = i)]= P(ξm+1 = k1 /(ξm

= i))

(2)

P(ξm+2 = k2 /(ξm+1 = k1))....P(ξm+n = j /(ξm+n1 = kn1))

 

 

 

 

 

По свойству однородности

 

 

 

 

 

 

 

 

 

 

 

P(ξm+1 = k1 /(ξm = i))= P(ξ1 = k1 /(ξ0 = i)),

 

 

P(ξm+2 = k2 /(ξm+1 = k1))= P(ξ2 = k2 /(ξ1 = k1)),

(3)

………………………………………………………..,

 

 

P(ξm+n = j /(ξm+n1 = kn1))

= P(ξn = j /(ξn1 = kn1))

 

Так как события S1, S2,.....Sr образуют полную группу,

то по формуле полной

вероятности

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P[ξ1 = k1 / ξ0 = i]P[ξ2 = k2 / ξ1 = k1]= P[ξ2 = k2 / ξ0 = i],

(4)

k1 =1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P[ξ2 = k2 / ξ0 = i]P[ξ3 = k3 / ξ2 = k2 ]= P[ξ3 = k3 / ξ0 = i],

(5)

k2 =1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

……………………………………………………………………….,

 

 

r

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P[ξn1 = kn1 / ξ0 = i]P[ξn = j / ξn1 = kn1]= P[ξn = j / ξ0 = i],

(6)

kn1 =1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставляя (2) и (3) в (1) и учитывая (4-6), получим

 

 

 

P[ξ

n+m

= j /(ξ

m

= i)]= P[ξ

n

= j /(ξ

0

= i)]= P(n) ,

 

 

 

 

 

 

 

 

 

ij

 

 

что и требовалось показать.

 

 

 

 

 

 

 

 

P(n0 )матрицы P(n0 )

Теорема. Если при

некотором n

все

элементы

 

 

 

 

 

0

 

 

 

 

 

 

ij

 

 

положительны, то существуют пределы

lim Pij(n) = b j , i, j =1,2...r , n→∞

где предельные вероятности b j не зависят от начального состояния и являются единственным решением системы уравнений

r

b j =1,=

j 1

r

bk Pkj = b j.

=k 1

Физический смысл этой теоремы заключается в том, что нахождения

системы в состоянии S j

практически не зависит от того,

в каком состоянии

она находилась в далеком прошлом.

 

 

 

 

 

Исследование случайных процессов зависит от вида рассматриваемых

состояний. Состояние

называется

несущественным,

если

существует

такое и такое , что

 

, но

 

для всех

, где

 

– число

тактов перехода. Таким образом0,

несущественное состояние характеризуется

 

0

 

 

,

 

тем, что из него можно попасть в некоторое другое состояние, но вновь вернуться в первоначальное состояние уже нельзя. Все состояния, отличные от несущественных, называются существенными. Говорят, что система обладает эргодическим свойством, если ее состояния принадлежат к одному существенному классу состояний, т.е. это свойство заключается в том, что объект, находящийся в момент t, в состоянии i, через достаточно большой промежуток времени возвращается в это состояние, т.е. если существуют вероятности b j . Эргодическое свойство играет большую роль при

исследовании стационарных случайных процессов. Эффективным инструментом исследования таких процессов являются методики, опирающиеся на теорему Биркхофа-Хинчина.

Теорема Биркхофа-Хинчина. Если непрерывный стационарный процесс имеет конечное математическое ожидание, то с единичной

вероятностью существует предел

lim

1

.

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

Уравнения Колмогорова для вероятностей состояний

Исследование непрерывной цепи Маркова основывается на уравнениях Колмогорова. Пусть имеется непрерывная цепь Маркова, т.е. система может

находиться в

дискретных состояниях S1, S2,....Sr

переход в

которые

осуществляется в любой случайный момент времени. Пусть

 

вероятность того, что в момент t система находится в состоянии

,

и пусть

требуется найти

алгоритм, описывающий изменение

всех

в

любой

r

момент времени. В любой момент времени Pi (t)=1 . Введем понятие

i=1

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

Плотностью вероятности перехода из состояния в состояние называется величина

lim ,

где – вероятность того, что система, находящаяся в момент времени t

в состоянии

 

, за

время

 

перейдет

в

состояние

. С

точностью

до

бесконечно малых высшего порядка

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

Будем

 

рассматривать

однородные

непрерывные

марковские

цепи

(

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

),

характеризующие

процессы

гибели

и

 

 

 

 

размножения.

Получим

дифференциальные

 

уравнения

 

для

вероятностей

. Прежде всего, найдем дифференциальное уравнение для

вероятности

 

начального состояния .

Дадим t малое приращение

и

найдем вероятность того, что в момент

 

 

система будет находиться в

состоянии

 

Это событие может произойти двумя способами:

 

 

 

– в момент .

система была уже в состоянии

и за время

не перешла в

состояние

;

 

 

 

 

 

и за время

перешла в состояние

– в момент

система была в состоянии

 

.

Вероятность

первого

события

 

по

теореме

умножения

равна

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

1

на условную вероятность не

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

, т.е.

 

 

 

 

 

 

 

 

 

 

1

 

 

.

 

 

 

 

 

Аналогично, вероятность второго события равна

.

Тогда по теореме сложения

1 ,

так что искомое дифференциальное уравнение имеет вид

lim

 

 

 

 

 

.

(1)

 

 

 

 

 

Найдем дифференциальное уравнение для вероятности

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

состояния .

Дадим t малое приращение и найдем вероятность того, что в

момент

система будет находиться в состоянии

 

.

Это событие может

произойти тремя способами:

 

 

и за время

не перешла ни в

– в момент

система была в состоянии

 

состояние

, ни в состояние

;

и

за

время

перешла

в

– в момент

система

была

в

состоянии

состояние

 

 

 

 

 

 

 

 

 

 

 

– в момент.

система

была

в

состоянии

и

за

время

перешла

в

состояние

 

 

 

 

 

 

 

 

 

 

 

Как и.

в предыдущем случае, вероятность первого события

определяется как

 

1

,

,

,

 

 

 

 

 

 

 

 

 

 

 

 

 

а вероятности второго и третьего событий равны

 

,

 

 

 

 

 

 

 

 

, ,

,

 

 

 

 

 

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

возможность получения функций

, имеет вид

 

,

 

 

 

 

,

 

,

,

,

(2)

где i =1,2...r .

 

 

 

 

 

 

 

 

Пример. Пусть система состоит из основного элемента

и двух резервных:

При отказе элемента

в работу включается

, при отказе

До, включения.

каждый из

резервных

элементов находится

в холодном.

резерве и отказать не может. Интенсивность потока отказов основного

элемнта

Все потоки отказов пуассоновские. Требуется

определить надежность; .

системы.

Решение. Представим процесс, протекающий в системе, как марковский случайный процесс с непрерывным временем и с дискретными состояниями:

– работает элемент

,

– работает элемент

– работает элемент

,

– не работает ни один

,элемент.

Система уравнений Колмогорова для таких состояний имеет вид

,

,

,

 

1.

 

 

 

,

 

 

0

1

из первого

 

 

причем,

При начальном условии

уравнения находим

.

 

 

 

 

 

0

 

0 дает

Интегрирование второго уравнения с начальным условием

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Аналогично,

.

Из условия полноты находим

1.

Тогда надежность системы равна сумме соответствующих вероятностей:

.

§3. Пуассоновский случайный процесс

Рассмотрим события, которые могут происходить в каждый момент непрерывно меняющегося времени. Например, к таким событиям можно отнести регистрацию частиц счетчиком Гейгера.т.е. рассмотрим марковский процесс с дискретным множеством состояний и непрерывным временем, называемый процессом Пуассона.

Пусть ξ(t) число появлений событий на промежутке [0,t]. Будем считать, что ξ(t) определено для всех t и является неотрицательной

целочисленной случайной величиной. Тогда для t1 < t2 число появлений событий

на (t1,t2] равно ξ(t2)ξ(t1).

Пусть А – состояние системы, в котором за промежуток времени (t + t) наступит m событий. Пусть H1, H2, H3,...., Hm состояния системы, в которых

за промежуток времени t наступает 0,1,2….m событий соответственно. Введем некоторые предположения относительно совокупности рассматриваемых событий.

1)Предположим, что H1, H2, H3,...., Hm образуют полную группу событий

исобытие А может наступить при условии появления только одного из них. Тогда по формуле полной вероятности

m

 

p(A)= p(Hi )p(A / Hi ).

(1)

i=0

Здесь предполагается, что условные вероятности p(A / Hk ) не зависят от

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

2) Предположим, что рассматриваемый марковский процесс обладает свойством однородности, т.е. вероятность наступления k событий за любой промежуток времени t зависит только от числа k и от длительности этого промежутка времени t и не зависит от начала его отсчета

p(A)= pm(t + t),

(2)

p(Hk )= pk (t),

 

k= 0,1,2...m.

3)Предположим, что условные вероятности p(A / Hk ) удовлетворяют условию

p(A / Hk )= pmk ( t).

 

(3)

Подставляя (2) и (3) в выражение (1), получим

 

 

m

 

 

pm(t + t)= pmi (t)pi ( t).

 

(4)

i=0

 

 

4) Предположим, что за малый промежуток времени

t

вероятность

наступления 1 события приближенно пропорциональна

t

, а наступление 2

или более событий можно пренебречь, т.е.

 

 

p1( t)= λ t,

 

(5)

p0( t)=1 λ t,

 

 

lim pm( t)= 0, m 2, t0

где λ = const > 0.

Учитывая (5), из (4) получим

pm(t + t)= pm(t)(1λ t)+ pm1(t)λ t ,

и, значит,

 

pm(t +

t)pm(t)

= −λp

(t)+λp

m1

(t) .

(6)

 

 

 

 

 

t

 

m

 

 

 

 

 

 

 

 

 

 

 

Переходя в (6) к пределу

 

при t 0 ,

получим последовательность

дифференциальных уравнений первого порядка

 

 

 

p0(t)+λp0(t)= 0,

p1(t)+λp1(t)= λp0(t), p2(t)+λp2(t)= λp1(t),

…………………….

pm(t)+λpm(t)= λpm1(t).

Из (5) вытекает

(7)

(8)

(9)

(10)

lim p0( t)= p0(0)=1.

(11)

t0

 

Решая уравнение (7) при начальном условии (11), находим вероятность отсутствия событий за промежуток времени t

p

(t)= eλt .

 

0

 

 

Согласно (5)

 

 

lim p1(

t)= p1(0)= 0 .

(12)

t0

Решим уравнение (8) с начальным условием (12) p1(t)+λp1(t)= λp0(t)= λeλt .

Решение этого линейного уравнения равно p1(t)= u(t)v(t) , где

v(t)= eλdt = eλt ,

u(t)= λdt +C1 = λt +C1.

Из начального условия (12) находим константу интегрирования p1(0)= 0 = C1. Окончательно имеем

t =1c

p (t)= λteλt .

(13)

1

 

Для уравнения (9) в соответствие с (5) начальное условие имеет тот же вид p2(0)= 0 . Решим уравнение (9) с этим условием с учетом (13).

 

 

 

 

2

 

λt

,

 

 

p2

(t)+λp2(t)= λp1(t)= λ te

 

 

 

 

v(t)= eλdt = eλt ,

 

 

 

 

 

u(t)=

λ2 tdt +C2 =

(λt)2

+C2,

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

p

(t)

=

(λt)2

eλt .

 

 

 

(14)

 

 

 

 

 

 

2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Поступая аналогично, находим

 

 

 

 

 

 

 

 

 

 

p

(t)= (λt)m eλt .

 

 

 

(15)

 

m

 

m!

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Формула (15) описывает закон

 

распределения

(СВ)

числа событий,

наступающих за промежуток времени t, который называют случайным процессом Пуассона.

Пример. В наблюдениях Резерфорда и Гейгера радиоактивное вещество за промежуток времени t = 7.5c испускало в среднем n t = 3.87 α частицы.

Найти вероятность того, что за это вещество испустит хотя бы одну частицу.

Решение. Найдем вероятность того, что за произвольное время t вещество не испустит ни одной частицы

p0(t)= (λ0!t)0 eλt = eλt ,

где по условию задачи

λ = n tt = 37.87.5 = 0.516 .

Тогда вероятность того, что это вещество за t =1 c испустит хотя бы одну частицу равна

pm1(t =1)=1p0(t =1)=1e0.516 0.403 .

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]