主要内容

rc2poly.

将反射系数转换为预测滤波器多项式

句法

a = rc2poly(k)
[A,eFinal] = RC2Poly(k,R0)

描述

a = rc2poly(k)转换反射系数K.对应于预测滤波器多项式的晶格结构一种, 和A(1)=1。输出一种是长度的行矢量长度(k)+ 1。

[A,eFinal] = RC2Poly(k,R0)返回最终预测错误efinal.基于零滞后自相关,R0.

例子

全部收缩

考虑由一组反射系数给出的格子IIR过滤器。找到其等效预测过滤器表示。

K = [0.3090 0.9800 0.0031 0.0082 -0.0082];a = rc2poly(k)
A =1×6.1.0000 0.6148 0.9899 0.0000 0.0032 -0.0082

算法

rc2poly.计算输出一种使用Levinson的递归[1]。功能

  1. 设置输出矢量一种到了第一个元素K.

  2. 循环通过剩余的元素K.

    对于每个循环迭代一世a = [a + a(i-1:-1:1)* k(i)k(i)]

  3. 实施a = [1 a]

参考

[1] Kay,Steven M.现代光谱估计。Englewood Cliffs,NJ:Prentice-Hall,1988年。

扩展能力

在R2006A之前介绍