AnyMath 文档

rc2波利

反射系数转换成预测滤波器多项式。

库::`工程师`

语法

函数调用

* [参数:a]=rc2poly(<参数:k>>) -返回向量 [参数:a] 矢量得到的预测滤波器的系数 [参数:k] 晶格结构的反射系数。

* [参数:a],[参数:eFinal]=rc2poly([参数:k],[参数:r0]) -还返回最终预测的误差 [参数:eFinal] 基于零延迟自相关 [参数:r0].

争论

输入参数

# *k* — 反射系数

+ 向量资料

Details

指定为矢量的反射系数的列表。

数据类型

漂浮物32, 漂浮64</无翻译> 支持复数::是

# *r0* — 零延迟自相关

+ 标量,标量

Details

零延迟自相关,定义为标量。

数据类型

漂浮物32, 漂浮64</无翻译> 支持复数::是

输出参数

# *一个* — 预测滤波器的系数

+ 向量资料

Details

作为长度的向量串返回的预测滤波器的系数 ,在哪里 -矢量元素的数量 [参数:k]. 函数返回向量的第一列 a 平等的 1,因此, a[1]=1.

# *eFinal* — 最终预测误差的功率

+ 标量,标量

Details

最终预测的误差功率,作为标量返回。

例子:

预测滤波器形式的等效表示

Details

考虑由一组反射系数定义的晶格IIR滤波器。 让我们以预测滤波器的形式找到它的等效表示。

import EngeeDSP.Functions: rc2poly

k = [0.3090 0.9800 0.0031 0.0082 -0.0082]

a = rc2poly(k)
1×6 Matrix{Float64}:
 1.0  0.614816  0.989881  2.42604e-5  0.00315796  -0.0082

算法

功能 *rc2poly* 计算输出向量 [参数:a] 使用Levinson递归[1]。 功能:

  1. 设置输出向量 a 等于第一个元素 [参数:k].

  2. 循环遍历其余元素。 k.

    对于每个迭代 i 周期 a=[a + a[i−1:-1:1]*k[i];k[i]].

  3. 工具/工具 a=[1a].

文学作品

  1. Kay,Steven M. _Modern Spectral Estimation。_恩格尔伍悬崖,NJ:普伦蒂斯-霍尔,1988.