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

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

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

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

Добавлен: 17.10.2024

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

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

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

Обеспечить точное выполнение граничных условий (38.11) за счет выбора постоянных параметров Си С2, Х\ и х2 удается только в частных случаях. Как правило, поглощающие границы для функ­

ции 0(г/,

т) не совпадают с исходными поглохцающими границами

у — X(т)

и у = ц (т). Однако иногда эти границы получаются доста­

точно близкими, причем всегда их можно сделать пересекающимися по крайней мере в двух точках. В этих точках выполняются гра­ ничные условия (38.3) для функции w(y, т). Согласно методу коллокации искомую плотность распределения w(y, т) можно аппрок­ симировать функцией, которая удовлетворяет граничным условиям (38.3) не всюду, а только в отдельных точках. В качество такой аппроксимирующей функции иногда можно взять 0(г/, т), т. е. при­ нять

w(y,x)^Q(y, т).

(38.12)

Чтобы проверить возможность указанной приближенной за­ мены, нужно сравнить заданные поглощающие границы у = Х(т) и 1 /- ц (т ) с аналогичными границами для функции 0(г/, г). При фиксированных значениях параметров Си С2, х{ и х2 уравнение поглощающей границы для в (у, т) записывается в виде 0(1/, т) = 0, т. е.

f(;t, х\ т, y)=Cif(t,

х,

y) + c 2f(t, х2; т,

у).

(38.13)

Если, например,

д = с о ,

то С2 =

0.

Разрешая равенство

(38.13)

oj

носительно у,

получаем

уравнение поглощающей

границы

для

■0(г/, т) в виде г/ = ф(т; Сi, хл). Параметрами С; и Х\ можно распо­ рядиться так, чтобы кривая г/ = ф с границей поглощения У ~ Х(х) пересекалась в точках с заданными абсциссами x = t' и т — 1". Значения f и t" выбираются из условия наилучшего приближения кривых. В простейшем случае можно принять t' = t, t" = tu где / и t\— начальный и конечный моменты в рассматриваемом интер вале. Аналогично выбираются параметры С2 и х2, когда ограничена

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

а X — — оо. Если

ограничены обе функции X(т) и

д(т), то

из

(38.13)

следует

получить

два уравнения:

у —

= x])i(t; Ci,

С2,

хи х2)

и г/ = ф2(т; Сь С2,

хи х2). Параметры

Си

С2, Xj и х2 находятся из условий пересечения кривой У = ^h с ниж­ ней поглощающей границей у = Х(х) и кривой г/ = фг — с верхней поглощающей границей у = р(т) при x — t' и при x — t"\ значения t' и t" выбираются из условия наилучшего приближения кривых.

Если между заданными поглощающими границами и аналогии ными границами для функции 0 (у, т) получаются большие расхож­ дения, то интервал от t до т следует разбить на части. Условная плотность распределения w(t, х; th у) ординаты процесса в мо­ мент 11, конечный для первой части, находится изложенным выше способом. Если f\(x; t) — плотность распределения начальной орди­

296


наты, то плотность распределения W\(yi\ t{) ординаты непоглощен­ ного процесса в момент t\ находится по формуле

p-(t)

t) dx,

(38.14)

Щ(уй tx) = f w(t, х; tlt

X (t)

причем эта ордината может принимать любые значения от X{t\) до p(/i). Дальнейший расчет ведется, как для первого участка, когда начальный момент равен t\, а начальная ордината процесса У\.

Пример 38.1. Определить приближенные выражения для услов­ ной плотности распределения w{t, х\ т, у) непоглощенной части случайного процесса и для условной вероятности P(t, X; т) суще­ ствования процесса в момент т, если коэффициент сноса а (г, х) = = — ар2г/, а коэффициент диффузии Ь(т, у ) = (52, где а и р — за­ данные положительные постоянные. Граничные условия записы­

ваются в виде (38.3), причем

функция ц(т)

ограничена

при лю­

бом т > t. Рассмотреть два

случая, когда

Х = — со

и когда

Мт ) = - и М ,

Реш ен и е . При заданных коэффициентах сноса и диффузии

решением уравнения (38.1) с начальным условием (38.5) и гра­ ничными условиями (38.6) является функция

f(t,x ;

 

/ К

ехр

а (у xz)2

т,у) =

1 -

z2

 

 

1Лг(1 - Z 2)

где z —

со,

то согласно

(38.7)

 

 

Если X =

 

 

 

'9(1/,

т) = /(/, х;

г, y ) — C2f(t, х2; т,

у),

где С2 > 0 , а х2> р

((). Из условия в(у,

т) = 0 следует, что

ехр

а (у xz)2

= С2ехр

а (у X2Z)2

_ _ _ _ _ _

T ^ z 2

 

 

 

 

или, что то же самое,

 

 

 

— (У -

x z f + (у -

x2z)2=

1 - Z 2 In С2.

 

 

 

 

а

 

С помощью этого равенства получаем уравнение поглощающей гра­ ницы для функции '0(у, т) в виде

z -f- х 2)

( \ - z 2)\nC2

У = 2

z (х — лг2)

297


Выберем параметры х2 и

С2 так, чтобы

эта граница пересекалась

с поглощающей

границей

у — р (т)

при

х = ¥

и при т = t". Если

z' = e~aP й'~*>

а

z" = е ~а^ <r-t>,. то должшно

быть:

 

 

 

z'

+

 

+

[ Ю

Т ) 2] 1пС2

 

 

0 (^') — ~2

 

 

2o.z' (х — х 2)

 

 

 

 

 

 

 

 

 

 

 

(t") = у

( х + Х 2)

[1 — (z")2\1пС2

 

 

z" (х х 2)

 

 

 

 

 

 

 

 

Исключая 1пС2, приходим к равенству

 

 

 

 

 

 

2 p ( t ' ) - z ' ( x

+

x 2)

_

z” [1 — (г')2]

поэтому

2 v { t " ) - z " { x

+

x 2)

~

z'[\

- О Т ] ’

 

 

 

 

 

 

 

 

 

 

 

 

*2 =

■- .*+

 

 

 

 

{2' 11-

СТ1 И(О -

 

 

 

- 2 " [ 1 - О Т 1 И Г ) } ;

 

 

1пС2 =

I — ( z ' j 2 ^ 2 ~

12'(-*2 +

- * ) - 2 ц (0 1 =

 

=

l

- {z"f

(Хз ~~ х) ^

(А'2 +

~

(О] •

В

частном случае, когда

функция

р(т) равна

нулю при любом

х >

t, получается х2 = — х;

С2= 1. При этом уравнение поглощаю­

щей границы для функции 0 (у,

х) записывается в виде у — 0. Дан­

ная граница совпадает с исходной поглощающей границей, а по­ тому искомая плотность распределения точно совпадает с функ­ цией В (у, т), т. е.

w{t, х; х, y) = f(t, х; т, y) — f(t, — X; х, у).

В общем случае, полагая t ' = t и t " = t u где ^ — конечное зна­

чение аргумента т, получаем:

 

 

 

 

Х2=

х + 2\x(t) ;

 

ln С* =

 

Iх - t1

( * i ) ~

(01,

где Zi — е~а<Р

При этом уравнение поглощающей границы

для функции В (У, т)

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

 

 

У =

(*) + -^

-1 * 2)- (*t) ~

(01 =

= Sh [ а р Н .- ^ ТГ ^ (^ } Sh 1а?2 (Т “ 01

+ |l (t) Sh l*P! - Т)1) •

298


Если данная граница близка к у — р (т), то искомая плотность рас­ пределения непоглощенной части процесса

 

 

w (*, х;

х,

у) « йf (t, х;

т, у) -

C2f(t,

х2; т, у ) .

 

 

 

Условная вероятность существования процесса

 

 

 

 

 

P(t,x;

т ) =

Г w (t,x) т, у)

 

 

С2)~\-

 

I

 

/ / 2 а[р(т)-х2 :]\

_ Л( /2а [р (т)—хаг] \)

 

+ 2-|Ф^ т / 1 - z ’

) ~ С‘ (

 

 

I ) '

 

Когда ограничены обе функции >„(т) и р(т), согласно

(38.7)

б( У, x) = f(t,

х;

х, у)

 

хГ, х, y ) - C 2f(t, Х2; х, у).

 

Найдем постоянные С\

и С2 из условия обращения в нуль функции

6 (у, т)

на

границах

поглощения у =

%(т) =

— р (т)

и

У =

р (т)

при х =

t\.

В

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

с

равенствами

0[— p(^i),

^i] =

0 и

0[p('/i), t{\ — 0 получаются следующие уравнения:

 

 

 

 

Ci exp | aZ'^x_

 

[2p (tx) +

Zi (x + jct)]j +

 

 

 

 

+ ' C2exp j

az,(x__ ^

 

[2p (tx) +

2 , (x -f x2)] J-

 

1;

 

 

C\ exp

 

azx(x — x x\

[2p(A)

 

— (A'-f-X1)H-f-

 

 

 

 

 

1 — z\

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

C2exp

 

■azx(x x2)

|2p {tx) zx( x - f x 2)]

=

1.

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

Решение этих уравнений записывается в виде:

sh

" 2az, (X2 — X)li(t1)

j

 

" az\ (x2 x 2)

С,:

1 — z\

exp

~2az1(х2~ х ,) ц (^ )

1 — z\

sh

1 — z2

 

 

 

 

 

 

 

sh

"2a2, (x — Xx) p(^ )

 

 

' az\ ( x \ X 2)

 

l — z2

 

exp

sh

"2aZj (x2 — x x)i).(t,)~

1 — z\

1 — z\

 

 

 

 

 

 

 

Так как Л(т) = — р(т), то уравнения поглощающих границ для функции Q(y, т) можно записать в виде p = pi(t) и у = — pi(x).

299