文档

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返回一些南S，并在这些情况下提供警告消息。

提示

一个简单快捷的方法来检查是否一个它在单位圆内的所有根都是检查是否每个元素的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到元素我向量的一个．

    A = (A -k(i)*fliplr(A))/(1-k(i)^2);

参考文献

[1]凯，史蒂文M。现代光谱估计．恩格尔伍德悬崖，新泽西州:Prentice-Hall, 1988。

扩展功能

另请参阅

|||||

R2006a之前介绍

信号处理工具箱文档
金宝app

试试MATLAB、Sim金宝appulink和其他产品下载188bet金宝搏

现在就去审判吧