Эргодические процессы. Строго корректно характеристики случайных процессов оцениваются путем усреднения по ансамблю реализаций в определенные моменты времени (по сечениям процессов). Но большинство стационарных случайных процессов обладает эргодическим свойством. Сущность его заключается в том, что по одной достаточно длинной реализации процесса можно судить обо всех его статистических свойствах так же, как по любому количеству реализаций. Другими словами, закон распределения случайных величин в таком процессе может быть одним и тем же как по сечению для ансамбля реализаций, так и по координате развития. Такие процессы получили название эргодических (ergodic). Для эргодических процессов имеет место:
mX(t) = M{x(t)} =
x(t) dt, (9.1.16)DХ(t) = M{x(t) - mХ(t)]2} =
(x(t) - mХ(t))2 dt, (9.1.17)RX(t) = M{x(t)x(t+t)} =
x(t)x(t+t) dt. (9.1.18)Эргодичность является очень важным свойством случайных стационарных, и только стационарных процессов. Математическое ожидание эргодического случайного процесса равно постоянной составляющей любой его реализации, а дисперсия является мощностью его флюктуационной составляющей. Так как определение функций производится по ограниченным статистическим данным одной реализации и является только определенным приближением к соответствующим фактическим функциям процессов, целесообразно называть эти функции статистическими.
Свойства эргодичности могут проявляться только по отношению к двум первым моментам случайного процесса, что вполне достаточно для использования соответствующих методик исследования процессов. Практическая проверка эргодичности процесса обычно производится проверкой выполнения условия Слуцкого:
K(t) dt = 0. (9.1.19)Если ковариационная функция процесса стремится к нулю при возрастании значения аргумента (t), то процесс относится к числу эргодических, по крайней мере, относительно моментов первого и второго порядков.
Пример. Случайная функция задана выражением Z(t)=X(t)+Y, где X(t) - стационарная эргодичная функция, Y- случайная величина, некоррелированная с X(t). Эргодична ли функция Z(t)?
mz(t) = mz(x)+my, Kz(t) = Kx(t)+Dy.
Функция Z(t) стационарна, но не эргодична, так как при t Þ ¥ имеет место Kz(t) Þ Dy.
По формулам (9.1.16-9.1.18) можно вычислить моменты и для детерминированных процессов. Например, для периодической функции f(t)=a sin wt автоковариационная функция описывается выражением:
K(t) = (a2/2) cos wt.
Соответственно, для произвольной периодической функции, представленной рядом Фурье (разложенной по рядам Фурье):
K(t) = (1/2)
an2 cos wnt.Таким образом, автоковариационная функция периодической функции также является периодической функцией от аргумента t - величины временного сдвига.
9.2. Функции спектральной плотности [2,25,26].
Каноническое разложение случайных функций. Введем понятие простейшей случайной функции, которая определяется выражением:
X(t) = X×j(t), (9.2.1)
где Х - обычная случайная величина, j(t) - произвольная неслучайная функция. Математическое ожидание простейшей случайной функции:
mx(t) = M{Xj(t)}= j(t)×M{X}= j(t)×mx, (9.2.2)
где mx - математическое ожидание случайной величины Х. При mx = 0 математическое ожидание mx(t) также равно нулю для всех t и функция (9.2.1) в этом случае называется элементарной случайной функцией. Ковариационная функция элементарной случайной функции определится выражением:
Kx(t1,t2) = M{X(t1)X(t2)}= j(t1)j(t2)×M{X2}= j(t1)j(t2)×Dx. (9.2.3)
где Dx - дисперсия случайной величины Х.
Центрированную случайную функцию 0X(t) можно представить суммой взаимно некоррелированных элементарных случайных функций:
0X(t) =
Xi×ji(t), (9.2.4)Из взаимной некоррелированности элементарных случайных функций следует взаимная некоррелированность величин Xi. Математическое ожидание и ковариационная функция случайной функции 0X(t):
M{0X(t)}= M{
Xi×ji(t)}= 0.Kx(t1,t2) = M{0X(t1) 0X(t2)}= M{
Xi×ji(t1)Xj×jj(t2)}= ji(t1)jj(t2)M{XiXj}.В силу взаимной некоррелированности парных значений XiXj имеет место M{XiXj}= 0 при i ¹ j, и все члены суммы в последнем выражении равны нулю, за исключением значений при i = j, для которых M{XiXj}= M{Xi2}= Di. Отсюда:
Kx(t1,t2) =
ji(t1)ji(t2)Di. (9.2.5)Произвольная нецентрированная случайная функция соответственно может быть представлена в виде
X(t) = mx(t) + 0X(t) = mx(t) +
Xi×ji(t), (9.2.6)с математическим ожиданием mx(t) и с той же самой ковариационной функцией (9.2.5) в силу свойств ковариационных функций, где 0X(t) - флюктуационная составляющая случайной функции X(t). Выражение (9.2.6) и является каноническим разложением функции X(t). Случайные величины Xi называются коэффициентами разложения, функции ji - координатными функциями разложения. При t1 = t2 из (9.2.5) получаем функцию дисперсии случайной функции X(t):
Dx(t) =
[ji(t)]2×Di. (9.2.7)Таким образом, зная каноническое разложение (9.2.6) функции X(t), можно сразу определить каноническое разложение (9.2.5) ее ковариационной функции, и наоборот. Канонические разложения удобны для выполнения различных операций над случайными функциями. Это объясняется тем, что в разложении зависимость функции от аргумента t выражается через неслучайные функции ji(t), а соответственно операции над функцией X(t) сводятся к соответствующим операциям математического анализа над координатными функциями ji(t).
В качестве координатных функций разложения, как и при анализе детерминированных сигналов, обычно используются гармонические синус-косинусные функции, а в общем случае комплексные экспоненциальные функции exp(jwt). С учетом последнего предварительно рассмотрим особенности представления случайных функций в комплексной форме.
Комплексные случайные функции. В общем случае случайный процесс может описываться комплексной случайной функцией:
Z(t) = X(t) + jY(t), (9.2.8)
где X(t) и Y(t) - действительные случайные функции. Соответственно, математическое ожидание комплексной функции:
mz(t) = mx(t)+j×my(t). (9.2.9)
Заметим, что комплексное представление случайных функций не более чем удобная для анализа математическая форма их отображения, которая, с использованием выражений Эйлера, всегда может быть переведена в форму вещественных функций. Функции дисперсии, корреляции и ковариации должны представлять собой однозначные и неслучайные вещественные характеристики случайных процессов и функций, независимо от формы их математического представления. Это условие будет выполняться при использовании в выражениях моментов второго порядка операций умножения комплексных функций с комплексно сопряженными функциями. Так, выражение для вычисления корреляционной функции имеет следующий вид:
Rz(t1,t2) = M{Z(t1)×Z*(t2)}= M{[X(t1)+jY(t1)][X(t2)-jY(t2)]}=
= M{X(t1)X(t2)+Y(t1)Y(t2)+j×[Y(t1)X(t2)-X(t1)Y(t2)]} =
= Rx(t1,t2) + Ry(t1,t2) + j×[Ryx(t1,t2) - Rxy(t1,t2)]. (9.2.10)
Если действительные и мнимые части комплексной функции некоррелированы, то Ryx = Rxy = 0 и последний член выражения (9.2.10) также равен нулю.
Аналогичное выражение имеет место и для ковариационной функции. При t1 = t2 = t для функции дисперсии комплексной случайной величины имеем:
Dz(t) = M{|Z(t)-mz(t)|2} = Dx(t) + Dy(t), (9.2.11)
Все приведенные выражения в общем случае могут использоваться для любых комплексных случайных функций с любым физическим смыслом переменной t.
Финитное преобразование Фурье случайных функций. По аналогии с функциями детерминированных сигналов, отдельно взятая на интервале 0-Т реализация xk(t) стационарного случайного процесса 0X(t) может быть представлена в виде ряда Фурье: