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)では,最終予測誤差efinalに基づいてゼロラグ自己相関r0が返されます。

例

すべて折りたたむ

予測フィルター多項式からの反射係数の決定

ライブスクリプトを開く

予測フィルター多項式一个と最終予測誤差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では,以下の再帰関係が実装されます。

$\begin{matrix} k （ n ）＝ {一个}_{n} （ n ） \\ {一个}_{n - 1} （米）＝ \frac{{一个}_{n} （米） - k （ n ） {一个}_{n} （ n - 米）}{1 - k {（ n ）}^{2}} ，米＝ 1 ， 2 ， \dots ， n - 1 \end{matrix}$

この関係は,レビンソン再帰法(参考文献［1］)に基づいています。この関係を実装するために,poly2rcは一个の最初の要素を空行列とし,その後の要素を反転してループ状に上の計算を実行します。ループ我を反復するごとに,この関数では以下のことが実行されます。

k(我)を(我)と設定します。
上記の2番目の関係をベクトル一个の1 ~我の要素に適用します。
```
= (a - k(我)* fliplr (a)) / (1 k(我)^ 2);
```

参考文献

凯，史蒂文·m·现代光谱估计。Englewood Cliffs，新泽西:Prentice-Hall, 1988年。

拡張機能

C / c++コード生成
MATLAB®编码器™を使用してCおよびc++コードを生成します。

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

ツールボックス関数のコード生成に対する可変サイズの制限(MATLAB编码器)を参照してください。

参考

R2006aより前に導入

poly2rc

構文

説明

例

予測フィルター多項式からの反射係数の決定

制限

ヒント

アルゴリズム

参考文献

拡張機能

C / c++コード生成
MATLAB®编码器™を使用してCおよびc++コードを生成します。

参考

信号处理工具箱ドキュメンテーション

サポート

MATLABによる信号処理向けディープラーニング

poly2rc

構文

説明

例

予測フィルター多項式からの反射係数の決定

制限

ヒント

アルゴリズム

参考文献

拡張機能

C / c++コード生成MATLAB®编码器™を使用してCおよびc++コードを生成します。

参考

信号处理工具箱ドキュメンテーション

サポート

MATLABによる信号処理向けディープラーニング

C / c++コード生成
MATLAB®编码器™を使用してCおよびc++コードを生成します。