rc2poly

反射係数の予測フィルタ，多項式への変換

構文

A = rc2poly(k) [a,efinal] = rc2poly(k,r0)

A = rc2poly(k)では，(1)＝1を使用して，ラティス構造に対応する反射係数kが予測フィルタ，多項式一个に変換されます。出力一个は，長さ长度(k)+ 1の行ベクトルです。

[a,efinal] = rc2poly(k,r0)では，ゼロラグ自己相関r0に基づいて最終予測誤差efinalが返されます。

ラ@ @ブスクリプトを開く

一連の反射係数によって与えられるラティスiirフィルタ，を考えます。それと等価の予測フィルタ，の表現を求めます。

K = [0.3090 0.9800 0.0031 0.0082 -0.0082];A = rc2poly(k)

一个=1×61.0000 0.6148 0.9899 0.0000 0.0032 -0.0082

rc2polyでは，レビンソン再帰法(参考文献[1])を使用して出力一个が計算されます。その関数は以下を実行します。

出力ベクトル一个を,kの最初の要素に設定します。
kの残りの要素にいて，ルプを使用して繰り返し計算します。
ルプ我の反復ごとに，A = [A + A (i-1:-1:1)*k(i) k(i)]が計算されます。
A = [1 A]が実装されます。

[1] Kay, Steven M.现代光谱估计。恩格尔伍德悬崖，新泽西州:Prentice-Hall, 1988。

使用上の注意および制限:

R2006aより前に導入