Студопедия
Случайная страница | ТОМ-1 | ТОМ-2 | ТОМ-3
АвтомобилиАстрономияБиологияГеографияДом и садДругие языкиДругоеИнформатика
ИсторияКультураЛитератураЛогикаМатематикаМедицинаМеталлургияМеханика
ОбразованиеОхрана трудаПедагогикаПолитикаПравоПсихологияРелигияРиторика
СоциологияСпортСтроительствоТехнологияТуризмФизикаФилософияФинансы
ХимияЧерчениеЭкологияЭкономикаЭлектроника

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

Читайте также:
  1. G. Методические подходы к сбору материала
  2. I. Методический блок
  3. I. Методы исследования в акушерстве. Организация системы акушерской и перинатальной помощи.
  4. I. Общие методические требования и положения
  5. I. Организационно-методический раздел
  6. I.9.1.Хемилюминесцентный метод анализа активных форм кислорода
  7. I.Организационно-методический раздел

 

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

(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; просмотров: 204 | Нарушение авторских прав


Читайте в этой же книге: Цифровой спектральный анализ | Взвешивание. Свойства весовых функций | Задание |
<== предыдущая страница | следующая страница ==>
Паразитная амплитудная модуляция спектра| Использование программы SPTool (Signal Processing Toolbox) системы Matlab для задач спектрального оценивания

mybiblioteka.su - 2015-2024 год. (0.007 сек.)