Преобразуйте полином предсказательного фильтра в отражательные коэффициенты
k = poly2rc(a)
[k,r0] = poly2rc(a,efinal)
k = poly2rc(a) преобразует полином предсказательного фильтра a к отражательным коэффициентам соответствующей структуры решетки. a может быть действительным или комплексным, и a(1) не может быть 0. Если a(1) не равно 1, poly2rc нормирует полином предсказательного фильтра на a(1)K вектор-строка из размера length(a)-1.
[k,r0] = poly2rc(a,efinal) возвращает автокорреляцию нулевой задержки, r0, на основе итоговой ошибки предсказания, efinal.
Если abs(k(i)) == 1 для любого i, нахождение отражательных коэффициентов является плохо обусловленной проблемой. poly2rc возвращает некоторый NaNs и предоставляет предупреждающее сообщение в тех случаях.
Простой, быстрый способ проверять, если a имеет все его корни в модульном кругу, должен проверять если каждый из элементов k имеет величину меньше чем 1.
stable = all(abs(poly2rc(a))<1)
poly2rc реализации это рекурсивное отношение:
Это отношение основано на рекурсии Левинсона [1]. Реализовывать его, poly2rc циклы через a в обратном порядке после отбрасывания его первого элемента. Для каждой итерации цикла i, функция:
Наборы k(i) равняйтесь a(i)
Применяет второе отношение выше к элементам 1 через i из векторного a.
a = (a-k(i)*fliplr(a))/(1-k(i)^2);
[1] Кей, Стивен М. Современная спектральная оценка. Englewood Cliffs, NJ: Prentice Hall, 1988.