Оценка спектральной плотности

Метод обработки сигналов

В статистической обработке сигналов цель оценки спектральной плотности ( SDE ) или просто спектральной оценки состоит в том, чтобы оценить спектральную плотность (также известную как спектральная плотность мощности ) сигнала из последовательности временных выборок сигнала. ^[1] Интуитивно говоря, спектральная плотность характеризует частотный состав сигнала. Одной из целей оценки спектральной плотности является обнаружение любых периодичностей в данных путем наблюдения пиков на частотах, соответствующих этим периодичностям.

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

Обзор

Пример формы голосового сигнала и его частотного спектра

Периодическая форма волны ( треугольная волна ) и ее частотный спектр, показывающий «основную» частоту 220 Гц, за которой следуют кратные (гармоники) 220 Гц.

Спектральная плотность мощности музыкального сегмента оценивается двумя разными методами, для сравнения

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

Спектральный анализ может быть выполнен для всего сигнала. Альтернативно сигнал можно разбить на короткие сегменты (иногда называемые кадрами ), и к этим отдельным сегментам можно применить спектральный анализ. Периодические функции (такие как ) особенно хорошо подходят для этого подразделения. Общие математические методы анализа непериодических функций относятся к категории анализа Фурье . ${\displaystyle \sin(t)}$

Преобразование Фурье функции создает частотный спектр, который содержит всю информацию об исходном сигнале, но в другой форме. Это означает, что исходная функция может быть полностью восстановлена ​​( синтезирована ) обратным преобразованием Фурье . Для идеального восстановления анализатор спектра должен сохранять как амплитуду , так и фазу каждой частотной составляющей. Эти две части информации могут быть представлены как двумерный вектор, как комплексное число или как величина (амплитуда) и фаза в полярных координатах (т. е. как вектор ). Распространенным методом обработки сигналов является рассмотрение квадрата амплитуды или мощности ; в этом случае результирующий график называется спектром мощности .

Из-за обратимости преобразование Фурье называется представлением функции в терминах частоты, а не времени; таким образом, это представление в частотной области . Линейные операции, которые можно выполнить во временной области, имеют аналоги, которые часто легче выполнить в частотной области. Частотный анализ также упрощает понимание и интерпретацию эффектов различных операций во временной области, как линейных, так и нелинейных. Например, только нелинейные или изменяющиеся во времени операции могут создавать новые частоты в частотном спектре.

На практике почти все программное обеспечение и электронные устройства, генерирующие частотные спектры, используют дискретное преобразование Фурье (ДПФ), которое работает с выборками сигнала и обеспечивает математическую аппроксимацию полного интегрального решения. ДПФ почти всегда реализуется с помощью эффективного алгоритма, называемого быстрым преобразованием Фурье (БПФ). Массив компонентов ДПФ с квадратичной величиной представляет собой тип спектра мощности, называемый периодограммой , который широко используется для изучения частотных характеристик бесшумных функций, таких как импульсные характеристики фильтра и оконные функции . Но периодограмма не обеспечивает выигрыша в обработке при применении к шумоподобным сигналам или даже синусоидам с низким отношением сигнал/шум. Другими словами, дисперсия его спектральной оценки на заданной частоте не уменьшается по мере увеличения количества выборок, используемых в расчетах. Это можно смягчить путем усреднения по времени ( метод Уэлча ^[2] ) или по частоте ( сглаживание ). Метод Уэлча широко используется для оценки спектральной плотности (SDE). Однако методы, основанные на периодограмме, вносят небольшие погрешности, которые неприемлемы в некоторых приложениях. Другие альтернативы представлены в следующем разделе.

Техники

Многие другие методы спектральной оценки были разработаны для устранения недостатков базовой периодограммы. Эти методы обычно можно разделить на непараметрические , параметрические и , в последнее время, полупараметрические (также называемые разреженными) методы. ^[3] Непараметрические подходы явно оценивают ковариацию или спектр процесса, не предполагая, что процесс имеет какую-либо конкретную структуру. Некоторые из наиболее распространенных оценщиков, используемых для основных приложений (например, метод Уэлча ), являются непараметрическими оценщиками, тесно связанными с периодограммой. Параметрические подходы, напротив, предполагают, что лежащий в основе стационарный случайный процесс имеет определенную структуру, которую можно описать с помощью небольшого числа параметров (например, с помощью модели авторегрессии или модели скользящего среднего ). В этих подходах задачей является оценка параметров модели, описывающей случайный процесс. При использовании полупараметрических методов основной процесс моделируется с использованием непараметрической структуры с дополнительным предположением, что количество ненулевых компонентов модели невелико (т. е. модель разрежена). Подобные подходы могут также использоваться для восстановления недостающих данных ^[4] , а также реконструкции сигнала .

Ниже приводится неполный список методов оценки спектральной плотности:

• Непараметрические методы , для которых выборки сигналов могут быть неравномерно разнесены во времени ( записи могут быть неполными )
  • Спектральный анализ методом наименьших квадратов , основанный на подборе методом наименьших квадратов известных частот.
  • Периодограмма Ломба – Скаргла , приближение спектрального анализа методом наименьших квадратов
  • Неравномерное дискретное преобразование Фурье

• Непараметрические методы , для которых выборки сигналов должны быть равномерно разнесены во времени ( записи должны быть полными ):
  • Периодограмма , квадрат модуля дискретного преобразования Фурье
  • Метод Бартлетта представляет собой среднее значение периодограмм, взятых из нескольких сегментов сигнала, чтобы уменьшить дисперсию оценки спектральной плотности.
  • Метод Уэлча — оконная версия метода Бартлетта, в которой используются перекрывающиеся сегменты.
  • Мультиконусность — это метод, основанный на периодограмме, который использует несколько конусов или окон для формирования независимых оценок спектральной плотности с целью уменьшения дисперсии оценки спектральной плотности.
  • Анализ сингулярного спектра - это непараметрический метод, который использует разложение по сингулярным значениям ковариационной матрицы для оценки спектральной плотности.
  • Кратковременное преобразование Фурье
  • Критический фильтр — это непараметрический метод, основанный на теории информационного поля , который может работать с шумом, неполными данными и инструментальными функциями отклика.
• Параметрические методы (неполный список):
  • Оценка модели авторегрессии (AR), которая предполагает, что n- й образец коррелирует с предыдущими p- выборками.
  • Оценка модели скользящего среднего (MA), которая предполагает, что n -я выборка коррелирует с шумовыми условиями в предыдущих p выборках.
  • Оценка авторегрессионного скользящего среднего (ARMA), которая обобщает модели AR и MA.
  • Классификация множественных сигналов (МУЗЫКА) — популярный метод сверхразрешения .
  • Спектральная оценка максимальной энтропии — это всеполюсный метод, полезный для SDE, когда ожидаются сингулярные спектральные особенности, такие как острые пики.
• Полупараметрические методы (неполный список):
  • Оценка SParse Iterative Covariance Based Estimation (SPICE), ^[3] и более обобщенная SPICE. ^[5] ${\displaystyle (r,q)}$
  • Оценка итеративного адаптивного подхода (IAA). ^[6]
  • Lasso , аналогичный спектральному анализу методом наименьших квадратов , но со штрафом за разрежение. ^[7]

Параметрическая оценка

При параметрической спектральной оценке предполагается, что сигнал моделируется стационарным процессом , который имеет функцию спектральной плотности (SDF) , которая является функцией частоты и параметров . ^[8] Тогда проблема оценки становится проблемой оценки этих параметров. ${\displaystyle S(f;a_{1},\ldots ,a_{p})}$ ${\displaystyle f}$ ${\displaystyle p}$ ${\displaystyle a_{1},\ldots ,a_{p}}$

Наиболее распространенная форма параметрической оценки SDF использует в качестве модели авторегрессионную модель порядка . ^{[8] : 392} Последовательность сигналов, подчиняющаяся процессу с нулевым средним, удовлетворяет уравнению ${\displaystyle {\text{AR}}(p)}$ ${\displaystyle p}$ ${\displaystyle \{Y_{t}\}}$ ${\displaystyle {\text{AR}}(p)}$

${\displaystyle Y_{t}=\phi _{1}Y_{t-1}+\phi _{2}Y_{t-2}+\cdots +\phi _{p}Y_{t-p}+\epsilon _{t},}$

где — фиксированные коэффициенты и — процесс белого шума с нулевым средним и инновационной дисперсией . SDF для этого процесса ${\displaystyle \phi _{1},\ldots ,\phi _{p}}$ ${\displaystyle \epsilon _{t}}$ ${\displaystyle \sigma _{p}^{2}}$

${\displaystyle S(f;\phi _{1},\ldots ,\phi _{p},\sigma _{p}^{2})={\frac {\sigma _{p}^{2}\Delta t}{\left|1-\sum _{k=1}^{p}\phi _{k}e^{-2i\pi fk\Delta t}\right|^{2}}}\qquad |f|<f_{N},}$

с интервалом времени дискретизации и частотой Найквиста . ${\displaystyle \Delta t}$ ${\displaystyle f_{N}}$

Существует ряд подходов к оценке параметров процесса и, следовательно, спектральной плотности: ^{[8] : 452-453.} ${\displaystyle \phi _{1},\ldots ,\phi _{p},\sigma _{p}^{2}}$ ${\displaystyle {\text{AR}}(p)}$

• Оценщики Юла -Уокера находятся путем рекурсивного решения уравнений Юла-Уокера для процесса ${\displaystyle {\text{AR}}(p)}$
• Оценки Бурга находятся путем рассмотрения уравнений Юла – Уокера как разновидности обычной задачи наименьших квадратов. Оценщики Бурга обычно считаются превосходящими оценки Юла – Уокера. ^{[8] : 452} Бург связал их с оценкой спектра максимальной энтропии . ^[9]
• Оценщики методом наименьших квадратов вперед-назад рассматривают этот процесс как задачу регрессии и решают эту проблему, используя метод вперед-назад. Они конкурируют с оценщиками Бурга. ${\displaystyle {\text{AR}}(p)}$
• Оценщики максимального правдоподобия оценивают параметры, используя подход максимального правдоподобия . Это предполагает нелинейную оптимизацию и является более сложным, чем первые три.

Альтернативные параметрические методы включают подгонку к модели скользящего среднего (MA) и полной авторегрессионной модели скользящего среднего (ARMA).

Оценка частоты

Оценка частоты — это процесс оценки частоты , амплитуды и фазового сдвига сигнала при наличии шума с учетом предположений о количестве компонентов. ^[10] Это контрастирует с общими методами, описанными выше, которые не делают предварительных предположений о компонентах.

Один тон

Если нужно оценить только частоту одного самого громкого чистого тона сигнала , можно использовать алгоритм определения высоты тона .

Если доминирующая частота меняется со временем, то проблемой становится оценка мгновенной частоты , определенной в представлении время-частота . К методам мгновенной оценки частоты относятся методы, основанные на распределении Вигнера-Вилля и функциях неоднозначности более высокого порядка . ^[11]

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

Несколько тонов

Типичная модель сигнала состоит из суммы комплексных экспонент при наличии белого шума . ${\displaystyle x(n)}$ ${\displaystyle p}$ ${\displaystyle w(n)}$

${\displaystyle x(n)=\sum _{i=1}^{p}A_{i}e^{jn\omega _{i}}+w(n)}$.

Спектральная плотность мощности состоит из импульсных функций в дополнение к функции спектральной плотности, обусловленной шумом. ${\displaystyle x(n)}$ ${\displaystyle p}$

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

метод Писаренко

${\displaystyle {\hat {P}}_{\text{PHD}}\left(e^{j\omega }\right)={\frac {1}{\left|\mathbf {e} ^{H}\mathbf {v} _{\text{min}}\right|^{2}}}}$

МУЗЫКА

${\displaystyle {\hat {P}}_{\text{MU}}\left(e^{j\omega }\right)={\frac {1}{\sum _{i=p+1}^{M}\left|\mathbf {e} ^{H}\mathbf {v} _{i}\right|^{2}}}}$,

Метод собственных векторов

${\displaystyle {\hat {P}}_{\text{EV}}\left(e^{j\omega }\right)={\frac {1}{\sum _{i=p+1}^{M}{\frac {1}{\lambda _{i}}}\left|\mathbf {e} ^{H}\mathbf {v} _{i}\right|^{2}}}}$

Метод минимальной нормы

${\displaystyle {\hat {P}}_{\text{MN}}\left(e^{j\omega }\right)={\frac {1}{\left|\mathbf {e} ^{H}\mathbf {a} \right|^{2}}};\ \mathbf {a} =\lambda \mathbf {P} _{n}\mathbf {u} _{1}}$

Пример расчета

Предположим , от до — временной ряд (дискретное время) с нулевым средним значением. Предположим, что это сумма конечного числа периодических составляющих (все частоты положительны): ${\displaystyle x_{n}}$ ${\displaystyle n=0}$ ${\displaystyle N-1}$

${\displaystyle {\begin{aligned}x_{n}&=\sum _{k}A_{k}\sin(2\pi \nu _{k}n+\phi _{k})\\&=\sum _{k}A_{k}\left(\sin(\phi _{k})\cos(2\pi \nu _{k}n)+\cos(\phi _{k})\sin(2\pi \nu _{k}n)\right)\\&=\sum _{k}\left(\overbrace {a_{k}} ^{A_{k}\sin(\phi _{k})}\cos(2\pi \nu _{k}n)+\overbrace {b_{k}} ^{A_{k}\cos(\phi _{k})}\sin(2\pi \nu _{k}n)\right)\end{aligned}}}$

Дисперсия для функции нулевого среднего, как указано выше, определяется выражением ${\displaystyle x_{n}}$

${\displaystyle {\frac {1}{N}}\sum _{n=0}^{N-1}x_{n}^{2}.}$

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

Теперь, для простоты, предположим, что сигнал распространяется бесконечно во времени, поэтому мы переходим к пределу: если средняя мощность ограничена, что почти всегда имеет место в действительности, то существует следующий предел, который представляет собой дисперсию данных. ${\displaystyle N\to \infty .}$

${\displaystyle \lim _{N\to \infty }{\frac {1}{N}}\sum _{n=0}^{N-1}x_{n}^{2}.}$

Опять же, для простоты, перейдем к непрерывному времени и предположим, что сигнал бесконечно распространяется во времени в обоих направлениях. Тогда эти две формулы станут

${\displaystyle x(t)=\sum _{k}A_{k}\sin(2\pi \nu _{k}t+\phi _{k})}$

и

${\displaystyle \lim _{T\to \infty }{\frac {1}{2T}}\int _{-T}^{T}x(t)^{2}dt.}$

Среднеквадратичное значение равно , поэтому дисперсия равна Следовательно , вклад в среднюю мощность, исходящий от компонента с частотой , равен Все эти вклады в сумме составляют среднюю мощность ${\displaystyle \sin }$ ${\displaystyle 1/{\sqrt {2}}}$ ${\displaystyle A_{k}\sin(2\pi \nu _{k}t+\phi _{k})}$ ${\displaystyle {\tfrac {1}{2}}A_{k}^{2}.}$ ${\displaystyle x(t)}$ ${\displaystyle \nu _{k}}$ ${\displaystyle {\tfrac {1}{2}}A_{k}^{2}.}$ ${\displaystyle x(t).}$

Тогда мощность как функция частоты равна и ее статистическая кумулятивная функция распределения будет равна ${\displaystyle {\tfrac {1}{2}}A_{k}^{2},}$ ${\displaystyle S(\nu )}$

${\displaystyle S(\nu )=\sum _{k:\nu _{k}<\nu }{\frac {1}{2}}A_{k}^{2}.}$

${\displaystyle S}$— ступенчатая функция , монотонно неубывающая. Его скачки происходят на частотах периодических компонент , а значение каждого скачка является степенью или дисперсией этого компонента. ${\displaystyle x}$

Дисперсия — это ковариация данных относительно самих себя. Если мы теперь рассмотрим те же данные , но с задержкой , мы можем взять ковариацию с и определить ее как автокорреляционную функцию сигнала (или данных) : ${\displaystyle \tau }$ ${\displaystyle x(t)}$ ${\displaystyle x(t+\tau )}$ ${\displaystyle c}$ ${\displaystyle x}$

${\displaystyle c(\tau )=\lim _{T\to \infty }{\frac {1}{2T}}\int _{-T}^{T}x(t)x(t+\tau )dt.}$

Если он существует, то это четная функция: Если средняя мощность ограничена, то существует повсюду, конечна и ограничена тем, чем является средняя степень или дисперсия данных. ${\displaystyle \tau .}$ ${\displaystyle c}$ ${\displaystyle c(0),}$

Можно показать, что можно разложить на периодические составляющие с теми же периодами, что и : ${\displaystyle c}$ ${\displaystyle x}$

${\displaystyle c(\tau )=\sum _{k}{\frac {1}{2}}A_{k}^{2}\cos(2\pi \nu _{k}\tau ).}$

Фактически это спектральное разложение по различным частотам, связанное с распределением мощности по частотам: амплитуда частотной составляющей является ее вкладом в среднюю мощность сигнала. ${\displaystyle c}$ ${\displaystyle x}$ ${\displaystyle c}$

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

Смотрите также

• Многомерная спектральная оценка
• Периодограмма
• СигСпец
• Спектрограмма
• Частотно-временной анализ
• Частотно-временное представление
• Вероятность Уиттла
• Спектральное распределение мощности

Рекомендации

1. П. Стойка и Р. Мозес, Спектральный анализ сигналов, Прентис Холл, 2005.
2. Уэлч, П.Д. (1967), «Использование быстрого преобразования Фурье для оценки спектров мощности: метод, основанный на усреднении по времени по коротким модифицированным периодограммам», IEEE Transactions on Audio and Electroacoustics , AU-15 (2): 70 –73, Bibcode : 1967ITAE...15...70W, doi : 10.1109/TAU.1967.1161901, S2CID  13900622
3. Стойка, Петре; Бабу, Прабху; Ли, Цзянь (январь 2011 г.). «Новый метод оценки разреженных параметров в разделимых моделях и его использование для спектрального анализа данных с нерегулярной выборкой». Транзакции IEEE по обработке сигналов . 59 (1): 35–47. Бибкод : 2011ИТСП...59...35С. дои :10.1109/TSP.2010.2086452. ISSN  1053-587X. S2CID  15936187.
4. Стойка, Петре; Ли, Цзянь; Линг, Джун; Ченг, Юбо (апрель 2009 г.). «Восстановление отсутствующих данных с помощью непараметрического итеративного адаптивного подхода». 2009 Международная конференция IEEE по акустике, речи и обработке сигналов . IEEE. стр. 3369–3372. дои : 10.1109/icassp.2009.4960347. ISBN 978-1-4244-2353-8.
5. Свард, Йохан; Адалбьорнссон, Стефан Инги; Якобссон, Андреас (март 2017 г.). «Обобщение разреженной итеративной оценки на основе ковариации». Международная конференция IEEE по акустике, речи и обработке сигналов (ICASSP) 2017 г. IEEE. стр. 3954–3958. дои : 10.1109/icassp.2017.7952898. ISBN 978-1-5090-4117-6. S2CID  5640068.
6. Ярдиби, Тарик; Ли, Цзянь; Стойка, Петре; Сюэ, Мин; Баггерер, Артур Б. (январь 2010 г.). «Локализация источника и зондирование: непараметрический итеративный адаптивный подход, основанный на взвешенных наименьших квадратах». Транзакции IEEE по аэрокосмическим и электронным системам . 46 (1): 425–443. Бибкод : 2010ITAES..46..425Y. дои : 10.1109/TAES.2010.5417172. hdl : 1721.1/59588 . ISSN  0018-9251. S2CID  18834345.
7. Панахи, Ашкан; Виберг, Матс (февраль 2011 г.). «О разрешающей способности метода оценки DOA на основе LASSO». 2011 Международный семинар ITG по интеллектуальным антеннам. IEEE. стр. 1–5. дои : 10.1109/wsa.2011.5741938. ISBN 978-1-61284-075-8. S2CID  7013162.
8. Персиваль, Дональд Б.; Уолден, Эндрю Т. (1992). Спектральный анализ для физических приложений . Издательство Кембриджского университета. ISBN 9780521435413.
9. Бург, JP (1967) «Спектральный анализ максимальной энтропии», Материалы 37-го собрания Общества геофизиков-разведчиков , Оклахома-Сити, Оклахома.
10. Хейс, Монсон Х., Статистическая цифровая обработка сигналов и моделирование , John Wiley & Sons, Inc., 1996. ISBN 0-471-59431-8 . 
11. Лерга, Джонатан. «Обзор методов мгновенной оценки частоты сигнала» (PDF) . Университет Риеки . Проверено 22 марта 2014 г.

дальнейшее чтение

• Порат, Б. (1994). Цифровая обработка случайных сигналов: теория и методы . Прентис Холл. ISBN 978-0-13-063751-2.
• Пристли, МБ (1991). Спектральный анализ и временные ряды . Академическая пресса. ISBN 978-0-12-564922-3.

• Стойка, П.; Моисей, Р. (2005). Спектральный анализ сигналов . Прентис Холл. ISBN 978-0-13-113956-5.

• Томсон, ди-джей (1982). «Оценка спектра и гармонический анализ». Труды IEEE . 70 (9): 1055–1096. Бибкод : 1982IEEP..70.1055T. CiteSeerX  10.1.1.471.1278 . дои : 10.1109/PROC.1982.12433. S2CID  290772.