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

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

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

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

Добавлен: 17.10.2024

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

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

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

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

 

 

[у-х(т)Р

h(y, *)

1

2'2хМ

3x( t ) /

(34.28)

 

Математическое ожидание лг(т) и среднее квадратическое отклоне­ ние зх(х) для Х(ъ) могут быть определены с помощью равенства (34.20), если в нем заменить характеристическую функцию Е(и, т)

выражением вида (34.21) с параметрами х(%) и ох(т), а характе­ ристическую функцию начальной ординаты записать в виде

iu х (t) -

2

- iArx(t)

E(u, t) — e

 

(34.29)

После указанной подстановки из сравнения коэффициентов при

одинаковых степенях параметра и находим:

 

 

 

X

 

*

X

 

 

 

J’ аДЩИтг!

J аДтрсЬ)

 

х (т )=

x[t) е1

-J- I aQ{i)e%

d\\

(34.30)

 

X

 

t

т

 

 

 

 

*

a,(T|)dT

 

 

2 Jа,(т|) di]

2 j

 

З Д =

3 (* )е *

 

+ \ ь 0(1) е 1

di.

(34.31)

 

 

 

t

 

 

 

Для условного среднего квадратического отклонения нормаль­ ной случайной величины справедливо выражение

 

Зу/Х =

3 /

1 —

где Гху=

kXy

 

 

. Поэтому

 

 

 

ху

 

 

 

kху=

°х

~

ху

(34.32)

:2

(34.33)

у/х

При X — X(t) и У= X (т) по этой формуле с учетом (34.23) и (34.30) получаем следующее выражение для корреляционной функ­

ции марковского случайного процесса в моменты времени t

и т:

 

 

 

j а,(т,) dr,

 

 

KAt,

*) = *хУ)е'

(34.34)

Пусть исходное уравнение (34.1) имеет вид

 

д[_

д

 

Ё

 

д-

— а ду

( ) '/ ) -

(34.35)

2

264


где

а и

р — но равные нулю постоянные. При

этом

коэффициент

сноса я(т, у) — чу, коэффициент диффузии

Ь(т,

у) — |32,

т. е.

а0 =

0;

щ — —а; & о=Р2. На основании вышеизложенного

полу­

чаем, что при заданном значении х начальной ординаты случайная величина Y= Х (т) имеет нормальное распределение; условный за­

кон распределения случайной величины Y=

X (") при любом фик­

сированном значении х начальной ординаты

X (/)

также

нормаль­

ный, причем согласно (34.24)

 

 

 

 

 

у = xe-« x~v;

а2 =

_ | l

[ 1— е-2*(т-ч]

_

(34.36)

В этом случае

 

 

 

 

 

 

Л (у ;

*=) = /(< ,

■*;

у ) =

 

 

 

________ / а ___

р|Л:[1 — e - 2*6 -‘>j

а [у—хе "(т 6 )2

Э» [1 — е -2а(т—t)]

(34.37)

Если начальная ордината X(t) является нормальной случайной величиной с математическим ожиданием x(t) и средним квадрати­

ческим

отклонением зх (/), то

плотность

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

нормаль­

ной случайной величины Y= X (-) записывается

в

виде

(34.28),

где согласно (34.30)

и (34.31)

 

 

 

 

 

 

 

 

л(т) = х Й

Г И |

;

 

 

(34.38)

 

<%(*) =

ъ 2У)е~ъъ +

[1 — е—5ta(T-t)] .

 

 

 

 

 

 

 

 

 

 

Корреляционная функция

(34.34) при этом определяется формулой

 

 

Kx{t,

T) =

a2 (^)e-“6- ‘) .

 

 

(34.39)

Рассмотрим стохастическое

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

уравнение

 

 

;f(f) + a * (/)= p g (0

 

 

 

(34.40)

с постоянными коэффициентами а и (5. Согласно

(32.16) или

(32.23)

и (32.24) коэффициенты сноса и диффузии для этого урав

нения

следующие:

a(t, х ) = ах;

b(t, х )— р2,

а потому второе

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

При a > 0 после затухания переходного процесса решение урав­ нения (34.40) можно считать стационарной случайной функцией.

265


Ее спектральная плотность определяется с помощью равенства

=

(34.41)

где Ф(гсо) — передаточная функция уравнения (34.40); 5£(со) — спектральная плотность нормального белого шума

Так как Ф (to) = ,-m^

, а

5£(ш) =

~ - ,

то

 

 

 

го» -j- а

 

 

 

 

 

 

 

 

5х(со)

2к (о)2 -f- а2)

 

(34.42)

 

 

 

 

 

 

 

Соответствующая корреляционная функция

 

 

 

 

 

оо

 

 

 

 

 

 

 

/Сх(т)'=

Гe,a>rSx (ш) dw =

 

. е-ф|

(34.43)

 

 

 

 

 

 

j / 2о

 

Последнее

равенство

получается также

из

(34.39),

если заменить

т — t на

т

и дисперсию

з2 (/)

начальной ординаты на дисперсию

Кх (0) =

- у = - стационарной случайной функции.

 

 

у

 

 

 

 

 

 

Таким образом, корреляционная функция стационарного нор­ мального случайного процесса имеет вид (34.43). Справедливо и обратное утверждение о том, что нормальный процесс является марковским, если его корреляционная функция записывается в виде

Кх(т) = и2е~а1т 1.

Так как эта функция не

имеет производной

в точке т = 0, то

нормальный

марковский процесс недифференци­

руемый.

 

 

 

 

 

Если в (34.35)

коэффициент а =

0, то это уравнение принимает

вид

 

 

 

 

 

 

дз

_1_

 

 

(34.44)

 

2

 

 

 

 

 

 

Так как при этом ах= 0, то согласно

(34.25)

 

 

у — х ; ау =

(3 тf x — t .

(34.45)

Плотность распределения f\(y;

т)

случайной

величины Y— X(t)

при заданном значении х начальной ординаты X{t) и условная

плотность распределения f(t,

х\ т, у) при любом фиксированном х

совпадают и определяются формулой

 

 

(у—хр

М у ;

Х\ t, у ) -

_______ 1______ е 2p,(T_t) . (34.46)

0 / ^ 7 т — t)

266


Если начальная ордината X(t) является^ нормальной случайной ве­ личиной с математическим ожиданием x(t) и средним квадратиче­ ским отклонением ах (^), то одномерная плотность распределения fi(y; т) определяется формулой (34.28), в которой

x(x) = x ( t ) ;

З Д = У<*(*) + Р2( * - * ) ■

(34-47)

Пример 34.1. Отклонение X(t) руля высоты, производимое авто­ пилотом для ликвидации пульсаций ветра, характеризуемых слу­ чайной функцией ~{%(t) , приближенно описывается стохастическим дифференциальным уравнением

TX(t) + X(t) =

n( t),

 

где Т и Т — известные постоянные,

а

— нормальный белый

шум, для которого 1 (t) — 0; Кх(т) = б (£).

Определить условную плотность распределения f(t, х; т, у) ор­

динаты Y = X(т), если в начальный момент t отклонение X(t)

руля

высоты равно х.

Исходное

стохастическое

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

урав-

Р еш ен и е .

нение совпадает

с (34.40)

при а = -1у

и (3 = т

Воспользовав­

шись формулой (34.37), для искомой условной плотности распре­ деления получаем следующее выражение:

/ Г

X

 

 

t - t 1

'г'

у — хе

Т

Т

 

Xехр

§35. РЕШЕНИЕ УРАВНЕНИЙ КОЛМОГОРОВА МЕТОДОМ ФУРЬЕ

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

а(т, у) = а(у)с(х) ; Ь(т, у ) = Ь ( у ) с ( т).

(35.1)

267


В частности, если c ( T ) = const, то коэффициенты сноса и диффузии не зависят от времени т. Уравнение (34.1) при этом имеет вид

7 W ~ 7 T + ~ h lal:,Uly'

Т,1 “ °- <35'2>

Согласно методу разделения переменных Фурье предполагаем, что искомая плотность распределения f(y, т) является произведе­ нием двух функций, одна из которых зависит только от у, а дру­ гая — только от т. Тогда

 

 

!(У, 1) =

Х(У)'Ф),

 

 

(35.3)

где х(У) и % (т)— неизвестные функции. Подставив (35.3)

в (35.2)

и разделив результат

такой

подстановки

на %(г/)%(т),

приходим

к равенству

 

 

 

 

 

 

 

 

__

1

d x __

 

 

 

 

 

с (т) х (т) di

 

 

 

=

- ^ у Г { т ^ т 1 1’ ( У У г А у ) 1 - ^ г

1‘‘ Ш

т } -

(35.4)

Левая

часть данного

равенства

зависит только

от т, а

правая —

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

которую обозначим через

Тогда получим:

 

/Уу

 

(35.5)

- ^ - + А ф ) х ( т ) = 0 .

- § f

(у) /. (у)] -

2-|г

м у- (у)] + 2^ у) =

°- ■ (35-6^

Из (35.5)

находим

 

 

 

 

 

 

-* L .=

- l c ( x)d*,

 

 

 

У,

 

 

 

поэтому функцию х(т)

из (35.3)

можно принять равной

 

,

 

Т

 

 

 

 

. — X j' с(Е) d;

 

 

 

 

*W = e

1

,

(35.7)

где t — начальный момент времени т.

 

задаются гра­

При определении плотности распределения f(y, т)

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

268