34)Определение характеристик эргодической стационарной функции по одной реализации.

Ответ: Рассмотрим стационарную случайную функцию , обладающую эргодическим свойством, и предположим, что в нашем распоряжении имеется всего одна реализация этой случайной функции, но зато на достаточно большом участке времени Т. Для эргодической стационарной случайной функции одна реализация достаточно большой продолжительности практически эквивалентна (в смысле объема сведений о случайной функции) множеству реализаций той же общей продолжительности; характеристики случайной функции могут быть приближенно определены не как средние по множеству наблюдений, а как средние по времени t. В частности, при достаточно большом Т математическое ожидание т_х может быть приближенно вычислено по формуле (6.7) Аналогично может быть приближенно найдена корреляционная функция при любом . Действительно, корреляционная функция, по определению, представляет собой не что иное, как математическое ожидание случайной функции :

. (6.8) Это математическое ожидание также, очевидно, может быть приближенно вычислено как среднее по времени. Фиксируем некоторое значение и вычислим указанным способом корреляционную функцию k_х( ). Для этого удобно предварительно «центрировать» данную реализацию х(t), т. е. вычесть из нее математическое ожидание (6.7):

. (6.9)

Вычислим при заданном математическое ожидание случайной функции как среднее по времени. При этом, очевидно, нам придется учитывать не весь участок времени от 0 до Т, а несколько меньший, т. к. второй сомножитель известен нам не для всех t, а только для тех, для которых . Вычисляя среднее по времени указанным выше способом, получим: . (6.10) Вычислив интеграл (6.10) для ряда значений , можно приближенно воспроизвести по точкам весь ход корреляционной функции. На практике обычно интегралы (6.7) и (6.10) заменяют конечными суммами. Покажем как это делается. Разобьем интервал записи случайной функции на п равных частей длиной и обозначим середины полученных участков (рис. 6.9).

Рис. 6.9

Представим интеграл (6.7) как сумму интегралов по элементарным участкам и на каждом из них вынесем функцию из-под знака интеграла средним значением, соответствующим центру интервала . Получим приближенно:

или . (6.11)

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

и вычислим интеграл (6.10), деля интервал интегрирования

на равных участков длиной и вынося на каждом из них функцию за знак интеграла средним значением. Получим:

или окончательно

. (6.12) Вычисление корреляционной функции по формуле (6.12) производят для последовательно вплоть до таких значений т, при которых корреляционная функция становится практически равной нулю или начинает совершать небольшие нерегулярные колебания около нуля. Общий ход функции воспроизводится по отдельным точкам (рис. 6.10).

Рис. 6.10

Для того чтобы математическое ожидание т_х и корреляционная функция были определены с удовлетворительной точностью нужно, чтобы число точек п было достаточно велико (порядка сотни, а в некоторых случаях даже нескольких сотен). Выбор длины элементарного участка определяется характером изменения случайной функции. Если случайная функция изменяется сравнительно плавно, участок можно выбирать большим, чем когда она совершает резкие и частые колебания. Чем более высокочастотный состав имеют колебания, образующие случайную функцию, тем чаще нужно располагать опорные точки при обработке. Ориентировочно можно рекомендовать выбирать элементарный участок так, чтобы на полный период самой высокочастотной гармоники в составе случайной функции приходилось порядка 5–10 опорных точек. Часто выбор опорных точек вообще не зависит от обрабатывающего, а диктуется темпом работы записывающей аппаратуры. В этом случае следует вести обработку непосредственно полученного из опыта материала, не пытаясь вставить между наблюденными значениями промежуточные, т. к. это не может повысить точности результата, а излишне осложнит обработку.