Построения спектральной плотности мощности.

[Map]

Есть какой-то сигнал, а точнее просто массив точек. Нужно оценить спектральную плотность мощности и построить полосовой фильтр по срезке спектра.
Для оценки СПМ нужно использовать авторегрессионную модель порядка p. Индентификация p+1 параметра a выполняется путём решения р+1 уравнения Юла-Уолкера, решение которой выполняется по алгоритму Левинсона-Дурбина.
Я так понял что сначало нужно найти корреляционную матрицу сигнала, который должен не иметь тренда. Потом подставить эту матрицу, которая имеет цеплицев вид, в уравнения и, путём решения методом ЛД, должны получить параметры АР-модели и её порядок. Позже строим спектр и делаем срез.
Всё ли так? Потому что это новая для меня сфера и тут я мало что понимаю.