Каноническое разложение случайных процессов.

 

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

¥

Х(t) = mx (t) + å u k× j k (t) , ( 21 )

k=1

где m x (t) - математическое ожидание с.п.; u1, u2, ... , uk, ... - некоррелированные, центрированные случайные величины с дисперсиями D1, D2, ... , Dk, ... ;

j1 (t), j2 (t), ... , jk (t), ... - неслучайные функции аргумента t . Число членов в разложении ( 21 ) может быть как конечным так и бесконечным.

Такое представление случайных процессов было впервые предложено в конце 50 - х годов В.С. Пугачевым и оно позволяет довольно просто проводить анализ преобразований (линейных , нелинейных) случайных процессов.

Пусть с.п. Х(t) задан представлением ( 21 ). Найдем его характеристики:

- математическое ожидание Х(t)

¥

М [ Х(t) ] = m x (t) + å M( uk )× jk (t) = m x (t), т.к. M(uk) = 0; (22)

k = 1

- корреляционная функция Х(t)

 

Кх(t1, t2) = М[Х(t1) Х(t2) ] = M[(Х(t1) - mx(t1)) (Х(t2) - mx(t2))] =

¥ ¥ ¥ ¥

= M[ å uk×jk (t1) × å un×jn (t2) ] = M [ å å uk un jk(t1) jn(t2) ] =

k=1 n = 1 k = 1 n = 1

¥ ¥ ¥ ¥

= å å M(uk ×un) jk(t1) jn(t2) = å å jk(t1) jn(t2) Dk , (23)

k=1 n=1 k =1 n =1

здесь M ( uk un ) =0 при k ¹ n и M (uk un) = Dk при k = n.

Выражение (23) называется каноническим разложением корреляционной функции с.п. Х(t). В частности при t1 = t2 = t получаем

¥

Kx (t) = å jk2 (t) Dk . ( 25 )

k=1

Это выражение называется каноническим разложением дисперсии с.п. Х(t).

 

Операции анализа над случайными функциями

 

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

Не вдаваясь в математические аспекты операций дифференцирования и интегрирования случайных функций остановимся только на том, как числовые характеристики случайного процесса Х( t ) связаны с числовыми характеристиками его производной У( t) = dХ( t) / dt и интеграла

t

Z(t) = ò X(s) ds.

0

Если для случайного процесса Х( t ) его числовые характеристики будут

M Х (t) = m х (t) , К Х (t1, t 2) и D Х (t), то характеристики с.п. У( t) будут связаны с характеристиками с.п. X (t) следующим образом

 

M [ У (t) ] = d m x (t)/ dt,

K У (t 1, t2) = ¶2 Kx (t1, t2)/ ¶t1 ¶t2,

D[ У (t) ] = K У (t , t).

Числовые характеристики случайного процесса Z(t) соотносятся с характеристиками с.п. X (t) таким образом

t

M [ Z (t) ] = ò m x (s) ds,

0

t1 t2

K z (t 1, t 2) = ò ò K x (s 1, s 2) ds 1 ds 2 ,

0 0

 

D[ Z (t) ] = K Z (t , t).

 

 

Пример 4. Случайный процесс задан элементарной функцией У = Х× t + a, где

Х - равномерно распределенная на отрезке [ c,d ] случайная

величина; а - неслучайная величина. Найти характеристики

производной с.п.

 

Решение: На основании определений найдем М [y(t)] = mx t + a ;

Ky (t1, t2 ) = sx2 t1× t2. Далее по формуле (25,26) найдем

my ¢(t) = mx = const;

 

K y ¢ (t1, t2) = ¶2 (sx2 t1 t2) /¶t1 ¶t2= sx2 = const ;

 

D [y¢(t)] = lim K y ¢ (t1,t2) = sx2 .

t1®t2

 

Пример 5 . Случайный процесс X(t) задан каноническим разложением

X(t) = 1 + u×t + v×t2 с характеристиками Du = 3; Dv = 1.

Найти математическое ожидание и ковариационную функцию

производной данного процесса.

Решение. По определению канонического разложения имеем Mx(t) = 1.

Обозначим Z(t) = dX(t)/dt. Тогда M [ dX(t)/ dt ] = dMx(t) / dt = 0 ,

 

Кz (t1,t2) = ¶ 2 Kx (t1,t2) /¶t1 ¶t2. Найдем Кх (t1,t2) = M [ X0(t1) X0(t2)] =

 

= { X0(t) = ut + vt2 - центрированный с.п. }= M[(ut + vt2)(ut + vt2)] =

 

= t1t2M(u2) + (t1t22 + t12t2)M(uv) = (t1t2)2M(v2) = 3t1t2 + (t1t2)2 ,

 

т.к. M(u2) = Du = 3 ; M(v2) = Dv = 1 ; M(uv) = 0.

¶ ¶

Тогда Kz(t1,t2) = ¶2 [3t1t2 + (t1t2)2 ] / ¶t1 ¶t2 = 3 + 4t1t2 .

 

 

Пример 6. В условиях примера 5 найти математическое ожидание и дисперсию

t

процесса Z(t) = ò X(s) ds.

0

 

 

 

t t

Решение. По определению (28) найдем М[Z(t)] = ò Mz(s)ds = ò 1× ds = t ;

0 0

t1 t2

Kz(t1,t2) = ò ds1 ò (3s1s2 + s12 s22) ds2 = 3 (t1t2)2 /4 + (t1t2)3 / 9.

0 0

Дисперсию с.п. Z(t) найдем как

 

D[Y(t)] = lim Ky(t1,t2) = 3 t4 /4 + t6 /9.

t1®t2