rc2poly
反射係数の予測フィルター多項式への変換
構文
a = rc2poly(k)
[a,efinal] = rc2poly(k,r0)
説明
a = rc2poly(k)
では、a(1)
= 1
を使用して、ラティス構造に対応する反射係数 k
が予測フィルター多項式 a
に変換されます。出力 a
は、長さ length(k)
+ 1 の行ベクトルです。
[a,efinal] = rc2poly(k,r0)
では、ゼロラグ自己相関 r0
に基づいて最終予測誤差 efinal
が返されます。
例
アルゴリズム
rc2poly
では、レビンソン再帰法 (参考文献[1]) を使用して出力 a
が計算されます。その関数は以下を実行します。
出力ベクトル
a
を、k
の最初の要素に設定します。k
の残りの要素について、ループを使用して繰り返し計算します。ループ
i
の反復ごとに、a = [a + a(i-1:-1:1)*k(i) k(i)]
が計算されます。a = [1 a]
が実装されます。
参考文献
[1] Kay, Steven M. Modern Spectral Estimation. Englewood Cliffs, NJ: Prentice-Hall, 1988.
拡張機能
バージョン履歴
R2006a より前に導入