Преобразуйте отражательные коэффициенты в полином фильтра прогноза
a = rc2poly(k)
[a,efinal] = rc2poly(k,r0)
a = rc2poly(k) преобразовывает отражательные коэффициенты, k, соответствующий структуре решетки к прогнозу, фильтрует полиномиальный a с a(1) = 1. Вывод a является вектором - строкой из длины length(k) + 1.
[a,efinal] = rc2poly(k,r0) возвращает итоговую ошибку прогноза efinal на основе автокорреляции нулевой задержки, r0.
rc2poly вычисляет вывод a с помощью рекурсии Левинсона [1]. Функция
Устанавливает выходной вектор a на первый элемент k.
Циклы через остающиеся элементы k.
Для каждой итерации цикла i, a = [a + a(i-1:-1:1)*k(i) k(i)].
Implements a = [1 a].
[1] Кей, Стивен М. Современная спектральная оценка. Englewood Cliffs, NJ: Prentice Hall, 1988.