将反射系数转换为预测滤波器多项式
一个= rc2poly (k)
[a, efinal] = rc2poly (k, r0)
一个= rc2poly (k)
转换反射系数k
对应于预测滤波器多项式的晶格结构一个
,(1)
=1
.输出一个
行向量是长度吗长度(k)
+ 1。
[a, efinal] = rc2poly (k, r0)
返回最终的预测错误efinal
基于零滞后自相关,r0
.
rc2poly
计算输出一个
使用莱文森的递归[1].这个函数
设置输出向量一个
到第一个元素k
.
的其余元素循环k
.
对于每个循环迭代我
,A = [A + A (i-1:-1:1)*k(i) k(i)]
.
实现了A = [1 A]
.
Kay, Steven M。现代谱估计.恩格尔伍德悬崖,新泽西州:Prentice-Hall, 1988。