Использование оценки корреляционной функции

Предположим, что мы исследуем эргодический процесс и знаем одну реализацию на интервале от 0 до некоторого . Выше было показано, что по этим данным можно построить оценку корреляционной функции. Если бы мы знали полностью непрерывную корреляционную функцию , для оценки спектральной плотности можно было бы использовать преобразование Фурье: .

В реальности известны лишь значения в отдельных точках, поэтому последнюю формулу нужно перевести в дискретный вид, заменив интеграл на конечную сумму:

. (2)

Этот метод оценки спектральной плотности называют методом Блэкмана-Тьюки.

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

Окна

Чтобы исправить ситуацию, нужно как-то «сгладить» незнание корреляционной функции при больших и сделать оценку спектральной плотности более надежной. Для этого используются так называемые «окна» – четные функции, на которые умножается корреляционная функция перед тем, как применить к ней преобразование Фурье. Одно из простейших «окон» – окно Хэмминга:

.

На рисунке слева показано окно Хэмминга, а справа – исходная оценка корреляционной функции и результат применения к ней окна Хэмминга (красная линия):

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

Для оценки спектральной плотности с учетом окна применяют формулу, аналогичную (2):

.

(3)

Не стоит печалиться по поводу того, что окно вносит дополнительное искажение. Так или иначе, «окно» используется всегда. Фактически, усекая корреляционную функцию, мы применяем прямоугольное окно:

.

На следующем рисунке показаны оценки спектра сигнала, полученные при использовании прямоугольного окна (, синяя линия) и окна Хэмминга (, красная линия).

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