Каноническое разложение случайных процессов.
Каноническим разложением случайного процесса Х(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