также нормальная, а потому ее плотность распределения записы вается в виде
|
|
[у-х(т)Р |
h(y, *) |
1 |
2'2хМ |
3x( t ) / |
(34.28) |
|
2« |
Математическое ожидание лг(т) и среднее квадратическое отклоне ние зх(х) для Х(ъ) могут быть определены с помощью равенства (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 |
где |
а и |
р — но равные нулю постоянные. При |
этом |
коэффициент |
сноса я(т, у) — —чу, коэффициент диффузии |
Ь(т, |
у) — |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)] . |
|
|
|
|
|
|
2а |
|
|
|
|
Корреляционная функция |
(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) можно считать стационарной случайной функцией.
Ее спектральная плотность определяется с помощью равенства
где Ф(гсо) — передаточная функция уравнения (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) = |
- у = - стационарной случайной функции. |
|
|
у |
2а |
|
|
|
|
|
|
Таким образом, корреляционная функция стационарного нор мального случайного процесса имеет вид (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)
Если начальная ордината 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
Xехр
§35. РЕШЕНИЕ УРАВНЕНИЙ КОЛМОГОРОВА МЕТОДОМ ФУРЬЕ
Для решения уравнений Колмогорова могут использоваться все методы, применяемые к уравнениям в частных производных пара болического типа. Наиболее часто решение этих уравнений произ водится методом Фурье, суть которого состоит в разделении пере менных. Применительно ко второму уравнению Колмогорова дан ный метод может быть использован только в том случае, когда коэффициенты сноса и диффузии можно представить в виде
а(т, у) = а(у)с(х) ; Ь(т, у ) = Ь ( у ) с ( т). |
(35.1) |
В частности, если 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) нужно искать таким, чтобы выполнялись соответствующие гранич ные условия для функции %(у). При указанных условиях уравне-