最新のリリースでは,このページがまだ翻訳されていません。このページの最新版は英语でご覧になれます。

yulewalk

再帰型デジタルフィルターの设计

构文

并[b,A] = yulewalk(N,F,M)

说明

yulewalkでは,指定した周波数応答に最小二乘近似を使用して,再帰型IIRデジタルフィルターが设计されます。

并[b,A] = yulewalk(N,F,M)では,ベクトルFで与えられる周波数 - 振幅特性を近似するñ次のIIRフィルターのn + 1个个の系数を含んだ行ベクトルb一种が返されます。

  • Fは,0〜1の范囲で设定する周波数点のベクトルで,1はサンプリング周波数の1/2(ナイキスト周波数)に対応します。Fの最初の点は0で最后の点は1でなければなりません。中间の点はすべて増加する顺でなければなりません。周波数点を重复することは可能で,これは,周波数応答のステップに対応します。

  • は,ベクトルFで指定された点での希望する振幅応答を含むベクトルです。

  • Fは同じ长さでなければなりません。

  • 图(F,M)では,フィルターの形状が表示されます。

出力されるフィルター系数はžの降顺になっています。

ž 一种 ž = b 1 + b 2 ž - 1 + + b ñ + 1 ž - ñ 一种 1 + 一种 2 ž - 1 + + 一种 ñ + 1 ž - ñ

周波数応答を设定するときには,通过帯域から阻止帯域への急激な迁移を避けてください。最良のフィルター设计を得るには,迁移领域の勾配をいくつか试さなければならない场合もあります。

すべて折りたたむ

正规化されたカットオフ周波数0.6をもつ8次のローパスフィルターを设计します。その周波数応答をプロットし,それに,対応する理想的フィルターの応答を重ね合わせます。

F = [0 0.6 0.6 1];M = [1 1 0 0];并[b,A] = yulewalk(8,F,M);[H,W] = freqz(B,A,128);积(W / PI,ABS(h)中,F,M,' - ')xlabel'弧度频率(\欧米加/ \ PI)',ylabel大小传说(“尤拉 - 沃克”'理想'),传说boxoff

アルゴリズム

yulewalkでは,时间领域で最小二乘近似が行われます。指定した周波数応答の逆フーリエ変换によって计算された相关系数を使用し,拡张ユール·ウォーカー方程式によって分母系数が计算されます。分子を计算する场合,yulewalkでは以下のステップが実行されます。

  1. パワー周波数応答の加法分解に対応する分子多项式が计算されます。

  2. 分子多项式および分母多项式に対応する周波数応答が计算されます。

  3. スペクトル因数分解法を使用して,フィルターのインパルス応答が求められます。

  4. このインパルス応答に対する最小二乘近似により,分子多项式が求められます。

参考文献

[1]弗里德兰德,B。,和波阿斯波拉特。“修正的ARMA谱估计的尤拉 - 沃克法”。IEEE®交易对航空电子系统。卷。AES-20,第2号,1984年,第158-173。

拡张机能

R2006aより前に导入