Метод модифицированных периодограмм

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

(30)

называют модифицированной периодограммой.

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

Как было отмечено выше, для получения состоятельных оценок спектральной плотности мощности используется усреднение нескольких периодограмм. Среди известных методов наибольшее применение получил метод Уэлча. В этом методе K сегментов данных длины M перекрываются и периодограммы вычисляются по K взвешенным сегментам. Далее периодограммы нормируются на величину U, чтобы компенсировать потери энергии вследствие процедуры взвешивания.

(31)

Таким образом, оценка Уэлча спектральной плотности мощности представляется в следующем виде

(32)

где

(33)

– модифицированная периодограмма, вычисленная по i-ому сегменту.

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

(34)

где – истинная спектральная плотность мощности анализируемого процесса x(n), а

(35)

– дискретное преобразование Фурье оконной функции.

Из (34) следует, что математическое ожидание искомой оценки равно свертке истинной спектральной плотности мощности с квадратом модуля Фурье-преобразования последовательности окна.

Для вычисления дисперсии используется тот факт, что дисперсия среднего арифметического K независимых одинаково распределенных случайных величин равна произведению дисперсии любой из них (они равны между собой) на множитель 1/K. Следовательно,

(36)

или, учитывая то, что

(37)

Таким образом, дисперсия периодограммы Уэлча обратно пропорциональна числу усредняемых периодограмм и стремится к нулю с его ростом. Легко видеть, что увеличение объема выборки N приводит к увеличению M и K. Поэтому при стремлении N к бесконечности как смещение, так и дисперсия оценки Уэлча стремятся к нулю. Отсюда следует, что усредненная периодограмма является асимптотически несмещенной оценкой спектра мощности

При 50 %-ном перекрытии сегментов и использовании треугольного окна дисперсия оценки Уэлча будет определяться выражением

(38)

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

Дата добавления: 2015-07-17