Параметрические методы

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

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

Методы AR имеют тенденцию соответственно описывать спектры данных, которые являются “остроконечными”, то есть, данные, PSD которых является большим на определенных частотах. Данные во многих практических применениях (таких как речь) имеют тенденцию иметь “остроконечные спектры” так, чтобы модели AR были часто полезны. Кроме того, модели AR приводят к системе линейных уравнений, которая относительно проста решить.

Методы AR Signal Processing Toolbox™ для спектральной оценки включают:

Все методы AR приводят к оценке PSD, данной

$\hat{P} (f) = \frac{1}{F_{s}} \frac{ε_{p}}{{| 1 - \sum_{k = 1}^{p} {\hat{a}}_{p} (k) e^{- j 2 π k f / F_{s}} |}^{2}} .$

Различные методы AR оценивают параметры немного по-другому, приводя к различным оценкам PSD. Следующая таблица предоставляет сводные данные различных методов AR.

Методы AR

	Город	Ковариация	Измененная ковариация	Уокер Рождества
Характеристики	Не применяет окно к данным	Не применяет окно к данным	Не применяет окно к данным	Применяет окно к данным
Характеристики	Минимизирует прямые и обратные ошибки прогноза в смысле наименьших квадратов, с коэффициентами AR, ограниченными удовлетворить рекурсию L-D	Минимизирует прямую ошибку прогноза в смысле наименьших квадратов	Минимизирует прямые и обратные ошибки прогноза в смысле наименьших квадратов	Минимизирует прямую ошибку прогноза в смысле наименьших квадратов (также названный “Метод автокорреляции”)
Преимущества	Высокое разрешение для коротких записей данных	Лучшее разрешение, чем Y-W для коротких записей данных (более точные оценки)	Высокое разрешение для коротких записей данных	Выполняет, а также другие методы для больших записей данных
	Всегда производит стабильную модель	Способный извлечь частоты от данных, состоящих из p или более чистых синусоид	Способный извлечь частоты от данных, состоящих из p или более чистых синусоид	Всегда производит стабильную модель
	Всегда производит стабильную модель		Не переносит спектральное разделение строки	Всегда производит стабильную модель
Недостатки	Пиковые местоположения, очень зависящие от начальной фазы	Может произвести нестабильные модели	Может произвести нестабильные модели	Выполняет относительно плохо для коротких записей данных
	Может перенести спектральное разделение строки для синусоид в шуме, или когда порядок является очень большим	Смещение частоты для оценок синусоид в шуме	Пиковые местоположения, немного зависящие от начальной фазы	Смещение частоты для оценок синусоид в шуме
	Смещение частоты для оценок синусоид в шуме	Смещение частоты для оценок синусоид в шуме	Незначительное смещение частоты для оценок синусоид в шуме
Условия для неособенности		Порядок должен быть меньше чем или равен половине входного формата кадра	Порядок должен быть меньше чем или равен 2/3 входной формат кадра	Из-за смещенной оценки матрица автокорреляции гарантируется положительно-определенному, следовательно несингулярному

Метод AR Уокера Рождества

Скрипт Open Live Script

Метод Уокера Рождества АРА спектральной оценки вычисляет параметры AR путем решения следующей линейной системы, которые дают уравнения Уокера Рождества в матричной форме:

$[\begin{matrix} r (0) & r (1) & \dots & r (p - 1) \\ r (1) & r (0) & \dots & r (p - 2) \\ ⋮ & ⋮ & ⋱ & ⋮ \\ r (p - 1) & r (p - 2) & \dots & r (0) \end{matrix}] [\begin{matrix} a (1) \\ a (2) \\ ⋮ \\ a (p) \end{matrix}] = [\begin{matrix} r (1) \\ r (2) \\ ⋮ \\ r (p) \end{matrix}]$ .

На практике смещенная оценка автокорреляции используется для неизвестной истинной автокорреляции. Метод Уокера Рождества АРА приводит к тем же результатам как максимальное энтропийное средство оценки.

Использование смещенной оценки автокорреляционной функции гарантирует, что матрица автокорреляции выше положительна определенный. Следовательно, матрица является обратимой, и решение, как гарантируют, будет существовать. Кроме того, параметры AR, таким образом вычисляемые всегда, приводят к стабильной модели все-полюса. Уравнения Уокера Рождества могут быть решены эффективно с помощью алгоритма Левинсона, который использует в своих интересах Эрмитовую структуру Теплица матрицы автокорреляции.

Функция тулбокса pyulear реализует метод Уокера Рождества АРА. Например, сравните спектр речевого сигнала использование метода валлийцев и метода Уокера Рождества АРА. Первоначально вычислите и постройте валлийскую периодограмму.

load mtlb
pwelch(mtlb,hamming(256),128,1024,Fs)

Спектр Уокера Рождества АРА более сглажен, чем периодограмма из-за простой базовой модели все-полюса.

order = 14;
pyulear(mtlb,order,1024,Fs)

Метод города

Скрипт Open Live Script

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

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

Точность метода Города ниже для старших моделей, длинных записей данных и высоких отношений сигнал-шум (который может вызвать разделение строки или генерацию постороннего peaks в оценке спектра). Спектральная оценка плотности, вычисленная методом Города, также восприимчива к сдвигам частоты (относительно истинной частоты) следующий из начальной фазы шумных синусоидальных сигналов. Этот эффект увеличен при анализе коротких последовательностей данных.

Функция тулбокса pburg реализует метод Города. Сравните оценки спектра речевого сигнала, сгенерированного и методом Города и методом Уокера Рождества АРА. Первоначально вычислите и постройте оценку Города.

load mtlb
order = 14;
pburg(mtlb(1:512),order,1024,Fs)

Оценка Уокера Рождества очень похожа, если сигнал достаточно длинен.

pyulear(mtlb(1:512),order,1024,Fs)

Сравните спектр сигнала с шумом, вычисленного с помощью метода Города и валлийского метода. Создайте двухкомпонентный синусоидальный сигнал с частотами 140 Гц и 150 Гц, встроенных в белый Гауссов шум отклонения 0,1 ². Второй компонент имеет дважды амплитуду первого компонента. Сигнал выбирается на уровне 1 кГц в течение 1 секунды. Первоначально вычислите и постройте валлийскую оценку спектра.

fs = 1000;
t = (0:fs)/fs;
A = [1 2];
f = [140;150];
xn = A*cos(2*pi*f*t) + 0.1*randn(size(t));

pwelch(xn,hamming(256),128,1024,fs)

Вычислите и постройте оценку Города с помощью модели порядка 14.

pburg(xn,14,1024,fs)

Ковариация и измененные методы ковариации

Скрипт Open Live Script

Метод ковариации для AR спектральная оценка основан на минимизации прямой ошибки прогноза. Измененный метод ковариации основан на минимизации прямых и обратных ошибок прогноза. Функции тулбокса pcov и pmcov реализуют соответствующие методы.

Сравните спектр речевого сигнала, сгенерированного и методом ковариации и измененным методом ковариации. Сначала вычислите и постройте оценку метода ковариации.

load mtlb
pcov(mtlb(1:64),14,1024,Fs)

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

pmcov(mtlb(1:64),14,1024,Fs)

Документация