Преобразование коэффициентов отражения в полином фильтра прогнозирования
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)].
Орудия a = [1 a].
[1] Кей, Стивен М. Современная спектральная оценка. Энглвуд Клиффс, Нью-Джерси: Прентис-Холл, 1988.