主要内容

poly2rc

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

折叠所有页面

语法

k = poly2rc (a)
[k, r0] = poly2rc (a, efinal)

描述

k = poly2rc (a)转换预测滤波器多项式一个为相应晶格结构的反射系数。一个可以是真实的，也可以是复杂的，还有(1)不能是0。如果(1)不等于1， poly2rc通过对预测滤波器多项式进行归一化(1)．k行向量的大小是多少长度(a) 1．

[k, r0] = poly2rc (a, efinal)返回零滞后自相关，r0，根据最终的预测误差，efinal．

例子

全部折叠

打开生活的脚本

给定一个预测滤波器多项式，一个，以及最后一个预测错误，efinal，确定相应晶格结构的反射系数和零滞后自相关。

A = [1.0000 0.6149 0.9899 0.0000 0.0031 -0.0082];efinal = 0.2;[k, r0] = poly2rc (a, efinal)
k =5×10.3090 0.9801 0.0031 0.0081 -0.0082
r0 = 5.6032

限制

如果abs (k (i)) = = 1对于任何我，求反射系数是一个病态问题。poly2rc返回一些南，并在这些情况下提供警告消息。

提示

一个简单快捷的检查方法一个是否所有的根都在单位圆内是为了检查是否每个元素k大小小于1。

稳定= (abs (poly2rc (a)) < 1)

算法

poly2rc实现这个递归关系:

k （ n ） ＝ 一个 n （ n ） 一个 n − 1 （ 米 ） ＝ 一个 n （ 米 ） − k （ n ） 一个 n （ n − 米 ） 1 − k （ n ） 2 ， 米 ＝ 1 ， 2 ， ⋯ ， n − 1

这种关系是基于莱文森的递归[1]．要实现它,poly2rc遍历一个在丢弃第一个元素后，以相反的顺序排列。对于每个循环迭代我功能:

  1. 集k(我)等于(我)

  2. 将上面的第二个关系应用到元素1到元素1我向量的一个．

    = (a - k(我)* fliplr (a)) / (1 k(我)^ 2);

参考文献

Kay, Steven M。现代谱估计．Englewood Cliffs, NJ: Prentice-Hall, 1988。

扩展功能

另请参阅

|||||

之前介绍过的R2006a

信号处理工具箱文档

万博1manbetx

基于MATLAB的信号处理深度学习

基于MATLAB的信号处理深度学习

下载白皮书