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

ЕГОРОВ В.А.- Обработка траекторий случайных процессов

.pdf
Скачиваний:
93
Добавлен:
08.07.2017
Размер:
1.25 Mб
Скачать

Определение 2.7. Ковариационной функцией случайного процесса X(t) называются функция K(s, t) EX (s)X (t). ▼

Определение 2.8. Стационарным или строго стационарным случайным процессом называется процесс, распределение которого не зависит от сдвига по времени. В терминах конечномерных распределений это свойство процесса сводится к выполнению условия

P

(C) P

h, ..., tn h

(C)

t1, ..., tn

t1

 

для всех допустимых наборов t1, ..., tn, h T и для всех измеримых C из про-

странства S n .

Определение 2.9. Стационарным в широком смысле называется случай-

ный процесс с постоянным математическим

ожиданием, для которого

 

 

 

 

 

 

EXt h X

s h EXt Xs , т. е. K(t h, s h) K(t, s),

t, h T . ▼

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

Для стационарных в широком смысле случайных процессов всегда справедливо равенство K(t, s) K(t s, 0), т. е. ковариационная функция является по существу функцией одной переменной. В этом случае функция K1(t) K (t, 0) также называется ковариационной функцией стационарного в широком смысле процесса X (t) и для нее обычно используется прежнее обозначение K (t) .

Определение 2.10. Случайный процесс называется нормальным (гауссовским), если все его конечномерные распределения нормальны. ▼

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

Если имеется случайный процесс с K(s, t) EX (s) X (t), m(t) EX (t), то можно построить нормальный случайный процесс с теми же K(s, t) и m(t) . Пусть X (t) – случайный процесс.

Для каждого t определим два множества случайных величин F t и F t , где множество F t состоит из ограниченных случайных величин, являющихся функциями от величин {X s, s t}, а множество F t состоит из огра-

11

E(B t

ниченных случайных величин, являющихся функциями от величин

{X s, s t}.

Определение 2.11. Марковское свойство случайного процесса X (t) состоит в том, что при любом выборе момента времени t и любом выборе случайных величин B t F t и B t F t выполняется равенство

B t | X (t)) E(B t | X (t)) E(B t | X (t)).

Это определение можно записать в следующей форме: при любом t и любом выборе случайных величин B t F t и B t F t выполняется равен-

ство E(B t | F t ) E(B t | X (t)).

Замечание. Марковское свойство случайного процесса качественно можно выразить в следующей форме: для любого момента времени t при известном настоящем X (t) прошедшее и будущее (определяемые в терминах значений процесса) независимы.

Определение 2.12. Случайный процесс, обладающий марковским свойством, называется марковским. ▼

Теорема 2.1. Пусть X (t), t T R1, вещественный нормальный случайный процесс. Он является марковским тогда и только тогда, когда для

любых a, b, c T ,

a b c, выполняется равенство

 

K (a, c)

K (a, b)K (b, c)

,

 

K (b, b)

 

 

 

где K – ковариационная функция процесса X (t) .

Доказательство необходимости. Используя свойства условных математических ожиданий из 1.6, получаем

K (a, c) EX (a) X (c) EE( X (a) X (c) | F b)

E[ X (a) E( X (c) | F b )] E[ X (a) E( X (c) | X (b))]

 

E[ X (c) X (b)]

 

 

E[ X (a) X (b)] E[ X (c) X (b)]

 

E X (a)

X (b)

 

.

EX 2 (b)

 

 

 

 

EX 2 (b)

 

 

 

 

 

 

Доказательство достаточности в теореме 2.1 не приводится, поскольку оно значительно сложнее. ▼

Упражнение. Доказать, что если X (t) является стационарным гауссовским процессом, то

12

при T Z

K (n) K ρ|n|

для некоторых | ρ | 1, K

0

0 ,

 

0

 

 

при T R1

K (u) K0 e |uc| для некоторых c 0, K0 0.

3. СЛУЧАЙНЫЕ ФУНКЦИИ ВТОРОГО ПОРЯДКА (ДИФФЕРЕНЦИРОВАНИЕ)

3.1. Критерии непрерывности и дифференцируемости случайных процессов

Пусть X (t), t R1, − случайный процесс с конечными вторыми моментами. Введем обозначения K(s, t) EX (s)X (t), m(t) EX (t). Пусть H – гильбертово пространство случайных величин, заданных на одном вероятностном пространстве и обладающих конечными вторыми моментами. В про-

странстве H

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

X ,Y EXY

 

,

||X||2 E |X|2 . Сходимость случайных величин в дальнейшем

 

понимается в среднеквадратичном, т. е. по норме пространства H . Определение 3.1. Случайный процесс X (t) называется непрерывным в

точке t , если lim

 

 

 

X (t h) X (t)

 

 

 

0 . ▼

 

 

 

 

h0

 

 

 

 

 

 

 

 

Определение 3.2. Случайный процесс Y (t) называется среднеквадратической производной случайного процесса X (t) в точке t , если

lim

 

X (t h) X (t)

Y (t)

 

0.

 

 

 

h

h0

 

 

 

 

 

 

 

 

 

 

 

Если производная существует для любого t (a, b), то процесс X (t) называеся дифференцируемым на интервале (a, b) . ▼

Для среднеквадратичных производных используется обычное обозначе-

ние Y (t) X (t) dX (t) . dt

Теорема 3.1 (критерий существования предела у процесса X (t) ). Сле-

дующие утверждения равносильны: а) существует предел lim X (t) , б) суще-

t t0

ствует предел lim K (s,t) .

s,t t0

Доказательство. Докажем сначала, что из а) следует б).

13

Скалярное произведение в любом гильбертовом пространстве непре-

рывно

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

 

 

совокупности

 

 

 

 

 

 

 

переменных,

 

поскольку при

H

H

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X X0

, Y Y0 справедливы соотношения:

 

 

 

 

X ,Y X0,Y0

 

 

 

 

 

X X0, Y X0, Y Y0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X X0

 

 

Y

 

 

 

 

 

 

 

X0

 

 

 

 

 

 

 

Y Y0

 

 

 

0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Следовательно,

из условия {X (t) Y} следуют равенства lim K(t, s)

lim X (t), X (s)

 

Y , Y

K(Y , Y ),

т. е. справедливо соотношение б).

Докажем, что из б) следует а). Пусть s, t t0 . Тогда

 

 

 

 

||X (t) X (s)||2

X (t) X (s), X (t) X (s)

=

 

 

K (t, t) K (t, s) K (s, t) K (s, s) 0.

 

 

Из этого соотношения и полноты гильбертова пространства H следует существование предела процесса X (t) . При t t0.

Из теоремы 3.1 непосредственно следует следующая теорема.

Теорема 3.2 (критерий непрерывности). Следующие утверждения рав-

носильны: а) lim X (t) X (t0) , б)

lim

K (s, t) K (t0, t0) .

 

 

 

 

 

t t0

 

s, t t0

 

 

 

 

 

 

 

 

 

 

 

Теорема 3.3 (критерий дифференцируемости процесса X (t) ).

 

 

а) Производная X (t0 ) существует тогда и только тогда, когда существу-

ет предел

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

lim

1

K (t h, t k) K (t h, t

 

) K (t

 

, t

 

k) K (t

 

, t

 

) .

 

0

0

0

0

0

 

0

0

0

 

 

 

 

 

 

h, k 0 hk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

б) Чтобы процесс был непрерывно дифференцируем на интервале (α, β) , необходимо и достаточно, чтобы существовала непрерывная при t, s (α, β)

смешанная вторая производная 2K (t, s) .

t s

Доказательство. а) По случайному процессу X (t) и точке t0 построим случайный процесс Y (h) 1h ( X (t0 h) X (t0)), h 0. Существование произ-

водной случайного процесса X (t) в точке t0 равносильно существованию предела для процесса Y (h) при h 0 . В силу теоремы 3.1 существование

14

предела

для процесса

 

Y (h) равносильно

существованию

 

предела

lim

KY (h, k) , где KY (h, k) – ковариационная функция процесса Y (h) :

h, k t0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K (h, k)

1

K (t h, t k) K (t h, t

 

) K (t

 

, t

 

k) K (t

 

, t

 

) .

 

 

0

0

0

0

0

 

 

Y

 

 

 

 

hk

 

 

0

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

б) Если процесс X (t) непрерывно дифференцируем на указанном ин-

тервале,

то выражение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

K (t h, s k) K (t, s k) K (t h, s) K (t, s)

 

 

 

 

 

 

 

 

 

 

hk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X (t h) X (t)

,

X (s k) X (s)

 

 

 

 

 

 

 

 

(3.1)

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

сходится к X '(t),

X

'(s) при h, k 0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

С другой стороны, устремив в выражении (3.1) h к нулю, получим пре-

 

 

1

 

K (t, s k)

 

K (t, s)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

дел

 

 

 

 

 

t

 

 

t

. Частные производные ковариационной функ-

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ции непрерывны, поэтому,

устремив k к нулю, получим предельную функ-

цию

2K (t, s)

. Следовательно, функция

2K (t, s)

существует.

 

 

 

 

 

 

t

s

t

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

2K (t, s) , тогда левая часть (3.1), являясь разностным аналогом этой частной

t s

производной имеет предел при h, k 0. Дифференцируемость процесса X (t) теперь непосредственно следует из соотношения (3.1). ▼

Следствие. Если производная X (t), t (α, β), непрерывна, то

E( X (t) X (s)E( X (t) X (s)

E( X (t) X (s) E( X (t) X (s)

K (t, s)

K (t, s)

t

K (t, s)

 

 

 

s

 

 

 

.

 

2

 

 

K (t, s)

 

 

 

 

t s

 

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

15

При h 0 справедливы соотношения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

EX

[ X (t h) X (t)]X

(s)

 

K (t h, s) K (t, s)

 

K (t, s)

,

 

h

 

 

 

 

 

t

 

 

 

 

 

 

h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

[ X (t h) X (t)]X

(s)

 

 

 

 

 

 

EX

E`X (t) X

(s),

 

 

 

 

 

h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

которые доказывают следствие. ▼ Упражнение. Используя приведенное доказательство, выразить через

ковариационную функцию случайного процесса X (t) взаимную ковариационную функцию EX (t) X (s) случайных процессов X (t) и X (t) .

3.2. Свойства дифференцируемых функций

Дифференцируемые в среднем квадратичном процессы обладают некоторыми свойствами обычных дифференцируемых функций. В качестве примера рассмотрим следующую теорему.

Теорема 3.4. Пусть процесс X (t) непрерывен при t [α, β] и X (t) 0, t [α, β]. Тогда X (t) с вероятностью единица является постоянной случайной величиной на этом интервале.

Доказательство. Пусть H – гильбертово пространство случайных вели-

чин. Тогда для любого Y из H определим функцию

f (t) X (t),

Y

. Функ-

ция f (t) − непрерывная комплекснозначная функция на [α, β] и

f

 

(t) 0 ,

поэтому f (t) f (α), что равносильно соотношению

X (t) X (α), Y

0. По-

следнее соотношение справедливо для любого Y из H , поэтому X (t) X (α) с вероятностью единица. ▼

Случайные процессы со значениями в H не могут быть исследованы с помощью среднеквадратичных производных так, как исследуются с помощью производных обычные функции. Это связано с тем, что для процессов со значениями в H не определены отношения «больше» и «меньше». Как показывает приведенное доказательство теоремы 3.4, с помощью среднеквадра-

тических производных можно исследовать функционалы типа

f (t)

X (t), Y

или

g(t) X (t), Y (t) , используя соотношения

 

 

 

 

f (t)

X (t),Y и

 

 

 

 

 

 

g (t)

X '(t),Y (t) X (t), Y (t) .

 

 

 

 

16

Определение 3.3. Нормальный случайный процесс W (t), t 0, называется винеровским процессом, если W (0) 0 и ковариационная функция процесса равна K(s, t) σ2 min (s, t).▼

Упражнение. Доказать, что винеровский процесс имеет независимые приращения, а его траектории непрерывны и нигде не дифференцируемы в среднем квадратичном.

4. СЛУЧАЙНЫЕ ФУНКЦИИ ВТОРОГО ПОРЯДКА (ИНТЕГРИРОВАНИЕ)

4.1. Определение интеграла от случайного процесса

Пусть

X X (t) –

случайный процесс, определенный на

T R

или

T R1, П –

некоторое

разбиение

интервала [α, β] : α t

t

... t

n

β,

 

 

 

 

0

1

 

 

δ(П) max(ti ti 1), 0 i n 1

– диаметр

разбиения;

ξi [ti , ti 1],

i n 1.

 

 

 

 

 

 

 

 

 

n 1

 

 

 

 

 

 

 

Сумма SП X (ξi ) (ti 1 ti ) называется интегральной суммой.

 

 

 

i0

 

 

 

 

 

 

 

Определение 4.1. Среднеквадратическим интегралом называется сред-

неквадратичный предел

lim SП , не зависящий от выбора разбиений и то-

 

 

δ(П)0

 

 

 

 

 

 

β

чек ξi . Этот интеграл обозначается как обычный интеграл X (t) dt . ▼

α

Выше определен интеграл от случайного процесса по неслучайной мере Лебега. Позже будет определен другой интеграл от неслучайной функции по случайному процессу.

Определение 4.2. Случайный процесс X (t) называется интегрируемым на интервале [α, β] в среднем квадратичном, если для него существует среднеквадратический интеграл по этому интервалу. ▼

4.2. Свойства интеграла от случайного процесса

Теорема 4.1 (критерий сходимости Коши для интегральных сумм). Ин-

β

 

 

 

 

 

 

 

 

 

теграл X (t)dt существует тогда и только тогда, когда

lim

 

 

 

SП SΠ

 

 

 

0 ,

 

 

 

 

 

 

 

 

 

1

 

 

 

 

α

 

 

 

 

 

 

 

 

 

17

где предел берется по всем таким разбиениям , 1 интервала [α,β], для ко-

торых δ(S ) 0, δ(S 1 ) 0.

Доказательство теоремы 4.1 полностью повторяет доказательство аналогичной теоремы анализа для интеграла Римана. ▼

Теорема 4.2 (критерий интегрируемости случайного процесса в среднем квадратичном в терминах ковариационной функции).

β

Интеграл X (t) dt существует тогда и только тогда, когда существует конеч-

α

β β

ный двойной интеграл K (s, t) ds dt .

α α

Доказательство. Предположим, что интеграл существует. Пусть

β

S , S 1 две интегральные суммы для интеграла X (t) dt по разбиениям, схо-

α

дящимся к нулю. Не умаляя общности, можно считать, что все разбиения для интегральных сумм одинаковы. Тогда

 

|| S

S ||2

S S , S S

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

1

1

 

 

 

 

 

 

 

 

 

 

 

 

S , S

S

, S

 

S , S

 

S , S

 

.

 

 

 

(4.1)

 

 

 

 

 

1

 

 

1

 

1

 

 

 

 

1

 

 

 

 

 

Рассмотрим одно из слагаемых в (4.1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n1

 

 

 

 

n1

 

 

 

 

 

 

 

 

 

 

S

1

, S

 

 

X

i,1

) (t

t ), X

j

) (t

j

1

t

j

)

 

 

 

 

i0

 

 

i 1

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j0

 

 

 

 

 

 

 

 

 

 

n1n1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

EX i,1) X

j ) (ti 1 ti ) (t j 1 t j ).

 

 

 

 

 

 

 

 

 

 

i1 j1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Правая часть этого равенства является интегральной суммой для двой-

 

 

 

β β

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ного интеграла

K (s, t) ds dt, поэтому она

сходится

 

к этому

интегралу.

α α

Аналогичные рассуждения можно провести относительно всех скалярных произведениях в правой части равенства соотношений (4.1). Следовательно, все четыре скалярных произведения имеют один и тот же предел, равный

18

β β

K (s, t) ds dt, поэтому левая часть соотношения (4.1) сходится к нулю. В

α α

силу критерия сходимости Коши интеграл X (t)dt существует.

β

Наоборот, предположим, что интеграл X (t) dt существует. Из непре-

α

рывности скалярного произведения для любого гильбертова пространства

 

 

 

 

 

 

 

β

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

следует равенство

lim (S , S )

 

X (t) dt

 

 

 

.

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

α

 

 

 

 

 

 

 

 

 

 

 

С другой стороны, выражение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n1 n1

 

 

 

 

 

 

 

 

 

 

 

 

 

S

1

, S

 

 

K

, ξ

j, 1

) (t

t ) (t

j1, 1

t

j, 1

)

 

 

 

i0 j0

i

 

 

 

i1

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

β β

 

 

 

является интегральной суммой Римана для интеграла K (s, t) ds dt , поэто-

α α

 

 

β

 

2

 

 

 

 

 

 

 

му этот интеграл существует и равен

 

X (t) dt

 

.▼

 

 

α

 

 

Замечания.

1. В ходе доказательства теоремы 4.2 было установлено, что при выполнении ее условий справедливо равенство

 

2

 

 

X (t) dt

K(s, t) ds dt.

 

 

2. Как само определение интеграла, так и его свойства во многом аналогичны определению и свойствам обычного интеграла Римана. В частности, если процесс X (t) непрерывен на интервале [α, β], то он интегрируем этом интервале.

t

3. Интеграл с переменным верхним пределом Y (t) X (s) ds от ста-

α

ционарного случайного процесса X (t) может быть нестационарным процессом. Действительно, пусть X (t) − вещественный стационарный

19

случайный процесс с нулевым математическим ожиданием. Ковариационная функция процесса Y (t) равна

t s

K (u v) du dv

t s v

 

 

 

 

 

K (w) dw dv.

(4.2)

 

 

 

α α

 

α

α

 

 

D(Y (t)) не зависит от α , поэтому вычисляя дисперсию по формуле (4.2), по-

t

лучим равенство D(Y (t)) 2 (t τ) K (τ) dτ . Отсюда следует, что D(Y (t)) все-

0

гда зависит от t.

4.3. Задачи на интегрирование случайных процессов

Задача 4.1. Определить вероятность того, что производная Y(t) X (t) от нормального стационарного случайного процесса X(t) в фиксирован-

 

 

 

 

 

 

 

 

 

 

 

ный момент времени

будет иметь

 

значение

больше

 

5, если mX 0,

K X (τ) σ2e α|τ|

 

 

α

 

 

, σ2X 4,

 

 

 

 

cos (βτ)

sin (β |τ|)

 

α 1.1,

β 2.1.

 

 

 

 

β

 

 

 

 

 

 

Решение. Распределение нормального процесса в фиксированный момент времени определяется математическим ожиданием и дисперсией. Вычислим эти характеристики. Очевидно, mY 0,

 

d 2

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2

 

2

 

α|τ|

 

α

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

DY

 

 

K X (τ)

τ 0

σ X

 

 

β

 

)e

 

 

cos (βτ)

 

sin (β |τ|)

 

dτ2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

β

τ 0

 

σ2 2 β2 ) 20.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Искомая вероятность равна

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

x

2

 

 

 

1

 

 

 

 

 

P

 

 

 

 

exp

 

dx

(1 Φ(0.5))

0.3085.▼

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

40

2

 

2 5 2π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

Задача 4.2. Известно, что Y (t) X (τ) dτ . Определить дисперсию про-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

цесса Y (t) при t 20 , если

K

X

(τ) σ2 e α|τ|(1 α |τ|), σ2

10, α 0.51.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

 

 

 

 

X

 

 

 

Решение. При t 20

20