Файл: Ганин, М. П. Прикладные методы теории марковских случайных процессов учебное пособие.pdf

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

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

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

Добавлен: 17.10.2024

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

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

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

ращению числа частиц за время Ат на отрезке координатной оси от у до у + Ду. Приращение плотности газа за время Ат при этом пропорционально разности fi(y, т) — fi(y, т + Дт), а потому прира­ щение числа частиц за время Дт на отрезке длиной Ду пропорцио­

нально \fi(y;

x) — fi(y; т -f- Ах)]Ау. Но тогда справедливо равенство

 

[f1{У, х) — f1(у; х + Дт) ]Ду =

 

 

=

[Si(y +

Ау,

х) — S1(y,

т)]Дт.

(30.27)

Разделив обе

части

этого

соотношения

на АуАх, после

перехода

к пределу при Дт -» 0 и Ау

0 получим равенство (30.26), кото­

рое является уравнением непрерывности потока, или, что то же са­

мое, уравнением сохранения числа частиц.

 

По аналогии с

(30.25)

введем условный поток вероятности, т. е.

функцию, определяемую формулой

 

у) —

S(t, х; т, у) =

а(х, y)f(t, х; т,

 

 

 

 

 

(30.28)

Тогда второе уравнение

Колмогорова

(30.23)

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

df{t, х ; т,

у)

dS(t,

х; х, у)

(30.29)

 

дх

 

'

ду

 

 

 

§ 31. НАЧАЛЬНЫЕ И ГРАНИЧНЫЕ УСЛОВИЯ.

ВЕРОЯТНОСТНЫЕ ХАРАКТЕРИСТИКИ

ПРЕБЫВАНИЯ ПРОЦЕССА В ЗАДАННОЙ ОБЛАСТИ

По принятой

классификации дифференциальных уравнений

в частных производных второго порядка уравнения (30.10), (30.11), (30.23) и (30.24) относятся к уравнениям параболического типа. Чтобы получить однозначное решение любого из указанных урав­ нений, необходимо задать начальное и граничные условия для со­ ответствующей искомой функции. Начальное условие представляет собой значение искомой функции в начальный момент времени. Для уравнения (30.24) начальный момент времени будем считать равным t, так что всегда т > /. Пусть в момент t ордината марков­ ского процесса является случайной величиной с известной плот­ ностью распределения f(x; t). Тогда начальное условие для иско­ мой одномерной плотности распределения f\{y, т), являющейся ре­

шением дифференциального уравнения

(30.24), записывается в виде

Ы у; "OIt- i

t).

(31.1)

Иногда начальная ордината рассматривается не как случайная ве­ личина с плотностью распределения /(х; /), а считается известной и равной заданному значению х. Так как это единственное значе-

231


нне начальная ордината принимает с вероятностью 1, то плотность распределения f(x, t) вырождается в дельта-функцию 6 ( у х). Начальное условие (31.1) для искомой функции fi(y; т) в этом случае принимает вид

/1 (у; т)|х=1 = о (у — х).

(31.2)

Для второго уравнения Колмогорова (30.23) начальный момент времени равен t. Еслп в искомой условной плотности распределе­ ния f(t, х\ х, у) принять х — t, то согласно (29.6) эта функция об­ ратится в дельта-функцию. Следовательно, начальное условие для функции /(/, х; х, у), являющейся решением дифференциального уравнения (29.23), записывается в виде

fit,

х: х, у)|т=1 = й (у

х).

(31.3)

При решении первого уравнения Колмогорова

(30.11) необхо­

димо знать значение

искомой условной

плотности

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

f(t, х; т, у) при t = т, т. е. в этом случае начальным является мо­ мент х. Начальное условие для условной плотности распределения f(t, х; х, у) при этом совпадает с (31.3).

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

ния F(t, х; т, у) в момент t x — 0.

Воспользовавшись равенством

(29.12), получаем

 

 

 

F(t, х- х, у ) = в (

у - х ) ,

(31.4)

где е х) — единичная

функция,

определяемая (29.13).

.

В дальнейшем будем

рассматривать только прямые дифферен­

циальные уравнения - (30.23) и (30.24), с помощью которых нахо­ дятся условная плотность распределения /(/, х; х, у) и одномерная плотность распределения fi(y, х) марковского случайного процесса при его развитии, т. е. в момент х > t. Граничные условия для ука­ занных искомых функций вытекают из физической сущности ре­ шаемой задачи и определяются ограничениями,накладываемыми на ординату у процесса яри любом т > t. Для обоих уравнений гранич­ ные условия записываются одинаково. Чтобы не повториться при

их

написании,

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

к

уравнению

(30.24). Для уравнения (30.23) граничные условия

получаются нз соответствующих условий для уравнения (30.24)

при замене одномерной плотности распределения fi (у; х)

условной

плотностью

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

f(/,

X; х, у)

и

потока вероятности

Si (г/; т) условным потоком вероятности S(t, х; х, у).

 

Еслп марковский случайный процесс X{t)

может принимать лю­

бые значения от

со до

со ,

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

уравнения

(30.23) и

(30.24)

справедливы при любых

значениях у.

Наиболее

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

232


процесса не может уходить в бесконечность и появляться из беско­ нечности. Вследствие этого поток вероятности *S] (г/; т) при боль­ ших по абсолютной величине значениях ординаты у должен рав­ няться нулю, а потому граничные условия записываются в виде

 

S t ( — с о ;

т) =

0,

Si ( с о ; х) =

0-

(31.5)

Наряду с

(31.5) при т >

t

для

одномерной плотности

распределе­

ния f1 (у;

т) выполняются более сильные условия:

 

 

/ х ( — с ° ; т ) =

0 ; / i ( c о ; т ) =

0 .

( 3 1 . 6 )

Равенства (31.6) выполняются всегда, так как они необходимы для того, чтобы неотрицательная функция f\(y; т) удовлетворяла усло­ вию нормировки

J fi(y, r )d y = 1 .

 

(31.7)

В некоторых случаях областью возможных

значений

марков­

ского случайного процесса X(t) является не вся

координатная

ось,

а ее часть от а до (5, причем может быть а = — ■

со или |3=

оо.

При

этом дифференциальные уравнения (30.23) и (30.24) справедливы только при а < у < |3. Так как все возможные значения случайного

процесса лежат в интервале от а до

(3, то кроме

(31.7)

справедливо

равенство

 

 

 

 

fi(y,

т )= 0 при у <

а н при у >

р.

(31.8)

Процесс, достигнув границы а пли

р, отражается внутрь области

A t? , расположенной

между границами у а и у =

р. Ордината

процесса не может выйти из области А з , а потому отсутствует по­ ток вероятности через ее границы. Вследствие этого граничные

условия для уравнения (30.24)

записываются в виде:

 

S ,(a; т) =

0; А (Р ; т) — 0.

(31.9)

При решении многих практических задач, связанных с опреде­ лением вероятности нахождения марковского случайного процесса

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

вобщем случае переменные, границы обозначим через Х(т) и р(т),

где Л(т) и р(т) — заданные функции, причем р(т) > Я (т), а об­ ласть, ограниченную этими линиями, — через Di^, Будем считать, что при достижении процессом любой границы сразу происходит его «поглощение». Если ордината процесса совпадает с какой-либо характеристикой физической системы, то каждое поглощение озна­ чает нарушение функционйрования системы. Результатом первого поглощения может быть прекращение дальнейшего функциониро­ вания системы или ее временное отключение с восстановлением

233


функционирования прп возвращении ординаты процесса в область Поэтому можно говорить о поглощении всего процесса или только его части. Примером системы первого типа является крыла­ тая ракета, совершающая полет на малой высоте над земной (вод­ ной) поверхностью. Если марковским случайным процессом при этом является высота полета ракеты, то полное поглощение этого процесса происходит при достижении нулевой высоты. Ко второму типу систем относится, например, цепь стрельбы, которая при не­ прерывно изменяющейся ошибке наводки отключается на время, в течение которого ошибка превосходит заданное значение. Подоб­ ные задачи возникают в теории надежности, когда для безаварий­ ной работы системы необходимо, чтобы определенная характери­

стика не выходила за допустимые пределы.

При наличии поглощающих границ Х(т) и д(т) одномерную плотность распределения fi (у; т) марковского случайного процесса

прп т >

t можно представить в виде

 

Л Су;

' ) = M

't) 3 [у — 4 х) 1 + р(1 (х) 5 [ у - 11(г)] + ® (у ;

х), (зп ю )

где рх(т) и

р (т) — вероятности поглощения процесса

границами

Х(т) п

|х(т)

соответственно, причем px(t) = p^t) — 0.

Функция

w(y; т)

из (31.10) является одномерной плотностью распределения

«пепоглощснной» части процесса при х > t, т. е. это плотность рас­

пределения вероятности нахождения процесса в области

D\^.

Для

значений т и у,

соответствующих внутренним точкам области

Div.,

функция w (у, т)

удовлетворяет уравнению вида

(30.24),

т. е.

 

dwf£'

+ ^

И Х, У)™(у; хЯ -

j

-

у ) ™( г >г)1 = ° -

 

 

 

 

 

 

(31.11)

Если для момента t известна плотность распределения /«(*; t) ор­ динаты процесса, то начальное условие для функции w(y; т) запи­ сывается в виде

иду; T)|T_t = Ы * ; t),

(31.12)

причем, так как ордината х лежит в области D ^,

 

f М * ; t ) d x = \ .

(31.13)

MO

 

При известной начальной ордпнате х процесса, такой, что

k(t) <

< х < р ( / ) , начальное условие для w(y, т) принимает вид

 

w(y, x)|x_ t = 3 (у — х).

(31.14)

Граничные условия для функции w(y\ г) отличаются от выра­ жений вида (31.9), так как на границах Л(т) и ц(т) имеется от­ личный от нуля поток вероятности. Известно, что при марковском

234


случайном процессе мгновенная скорость изображающей частицы бесконечная. Ограниченность смещения частицы за конечное время объясняется тем, что скорость беспрестанно с неограниченной ча­ стотой меняет знак, а потому движение происходит вперед и назад. Если случайная функция X(t) в момент т имеет значение у, то непосредственно до этого момента данная функция принимала зна­ чения, расположенные слева и справа вблизи значения у. Вслед­ ствие этого практически каждая траектория, проходящая вблизи границы А(т) пли р(т), приводит к поглощению процесса. Так как плотность распределения w(y, т) вероятностей нахождения про­ цесса в области А ,, при т > t отлична от нуля только для непогло­ щенных траекторий, то эта функция удовлетворяет следующим гра­ ничным условиям:

ш[Х(т); т] = 0;

йу[ц (т); т] = 0.

 

(31.15)

При заданном начальном условии (31.12) пли

(31.14)

и гранич­

ных условиях (31.15) плотность

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

w(y; т)

вероятно­

стей нахождения марковского случайного процесса в области Diа находится однозначно как решение дифференциального уравнения (31.11). Зная эту функцию, можно определить различные вероят­ ности и числовые характеристики случайного процесса.

Обозначим через Р(т) вероятность недостижения нп разу про­ цессом границ Я(т) и р(т) области Оха во* временном интервале от t до т. Данная вероятность получается с помощью w(y; т) в резуль­ тате интегрирования этой функции по всем возможным значе­ ниям у из области Dy.. Следовательно,

г-М

(31.16)

Р (т )= Г w(y; i)dy.

Х(т)

Проинтегрируем обе части равенства (31.10) по у от —оо до оо. Так как функция w(y: т) равна нулю вне области Di,,., то с учетом

(31.7) и (31.16) получаем

1 = Рх (*) + /V М 4- р М -

Тогда

 

 

 

 

 

 

.

(31.17)

Так

как рк(t) = р

(t) 0,

то P ( t ) = l . Однако по мере

поглоще­

ния

процесса с

ростом т

вероятность Р(т) нахождения

процесса

в области Охр. уменьшается, а потому условие нормировки вида (31.7) для функции ш(г/; т) не выполняется. В пределе получается

Р (оо) = 0.

Обозначим через fz(z) плотность распределения времени Z пре­ бывания случайного марковского процесса в области Ар. после мо­ мента t (до первого поглощения). Чтобы получить аналитическое

235

\