フィードバック構造
FIR (有限インパルス応答) フィルタは、各出力サンプルを現在と過去の 入力 のみの重み付け和として計算します。フィードバックはありません。インパルス応答の期間は有限です。
IIR (無限インパルス応答) フィルタは、以前の 出力 をフィードバックして計算に戻します:
y_n = Σ b_k · x_{n−k} − Σ a_k · y_{n−k}
再帰項 −Σ a_k · y_{n−k} はフィードバックを作成します。入力の単一インパルスはフィードバックループの周りで無期限にエコーします (安定している場合は幾何学的に減衰)。
なぜフィードバックを使うのか?
IIRフィルタはFIRフィルタよりもはるかに少ない係数で鋭い周波数選択性を達成できます。2極のIIRは、同じストップバンド減衰に必要な50係数のFIRを近似できます。
代償: 潜在的な不安定性。H(z)の極が安定性を決定します。すべての極は単位円の内側に厳密に存在する必要があります。
Hammingのフィードバックシャワーストーリー
Hammingはフィードバック不安定を説明するために鮮やかな個人的なストーリーを使用しました。
彼は疲れているときに自分の方向を見失わないようにするために、何度も同じホテルの部屋に泊まりました。配管工はシャワーに大径の温水パイプを設置していました。これらはノブを調整してから水温の変化を感じるまでの間に有意な 遅延 を作成しました。
毎朝、Hammingは同じパターンに従いました: 水が冷たすぎる → ホットを上げる → まだ冷たい → さらに上げる → 急に沸騰 → ジャンプして出る → 下げる → 繰り返す。
フィードバック経路の遅延は、彼の補正が常に過度になることを意味していました。彼は何度繰り返しても遅延に適応することはできませんでした。
エンジニアリング教訓: 不安定性はフィードバック経路の過度なゲイン、またはフィードバック経路の過度な遅延から生じます。どちらも同じハンティング動作として現れます。フィルタの用語では: 単位円上またはその外側の極はちょうどこの振動的または発散応答を生成します。
不安定性の特性化
Hammingは同じシャワー不安定性を2つの方法で分析できることを観察しました:
1. 彼の反応が強すぎました (補正動作の過度なゲイン)。
2. 彼の検出は遅延しすぎていました (システムが落ち着く前に浴槽に入るのが急いでいた)。
どちらの説明も同じ数学的結果を生じます: フィードバックループの極が単位円の外側に移動しました。
4つの古典的なファミリー
アナログフィルタ理論は、4つの古典的な設計ファミリーの周りで発展しました。各々は異なるトレードオフを表しています。これらのファミリーは双一次変換またはインパルス不変性を介して離散時間に変換されます。
Butterworth (最大限フラット)
パスバンド応答: |H(jω)|² = 1 / (1 + (ω/ω_c)^{2N})。単調に減少。パスバンドやストップバンドにリップルはありません。極はs平面 (またはZ平面で変換された円) の半径 ω_c の円上に位置します。与えられた次数Nに対する最もフラットなパスバンド。
Chebyshev Type I
パスバンドの等リップル、ストップバンドの単調。与えられた次数Nとリップルレベルに対して、Butterworthより鋭いカットオフを達成します。極はs平面の楕円上に位置します。
Chebyshev Type II
ストップバンドの等リップル、パスバンドの単調。周波数領域でType Iの鏡像。
楕円 (Cauer)
パスバンドとストップバンドの両方に等リップル。与えられた次数Nとリップルレベルに対して、パスバンドからストップバンドへの最も鋭い遷移を達成します。楕円関数を使用して極とゼロを最適に配置します。Hamming: 名前は楕円関数が導出に使用されるという事実から来ています。
基本的なトレードオフ
4つのファミリーはすべて同じ基本的なトレードオフをさまざまな方法で達成します: より高い次数Nはより鋭い遷移を与えます。リップルを許可する (Chebyshev、楕円) は同じNに対してより鋭い遷移を達成します。楕円は任意の与えられたNとリップル仕様に対する絶対的に最も鋭い遷移を達成します。
フィルタファミリーの選択
ファミリー間の選択は、アプリケーションが何を許容するかによります。
専門家の主張への質問
Hammingは、特定の専門家がIIR (再帰) フィルタに特定の特性を持つと主張していることを思い出しました。彼は自分自身に、これが本当に真実であるかどうか尋ね、反例を見つけました。
彼のポイント: 専門家は学校で吸収した請求を運ぶことが多く、現在の問題の文脈でそれらを再検討することはありません。あなたがあなたに言われていることが本当に真実であるかどうかを自分自身に尋ねるなら、それがどうであるか、または非常に発展した分野でさえ偽りに達しているのは驚くべきです。
反例は、通常設計するフィルタの種類ではありませんでしたが、請求を浅いことを証明しました。普遍的な請求を反証するために1つの反例で十分です。
実際のIIR設計
Hammingは、IIRフィルタ理論の多くを独立して開発していたことに気づきました。異なる問題を解決している間: 数値常微分方程式の安定した補正者の公式を導出します。
補正者の公式の形: y_n = Σ a_k · y_{n−k} + Σ b_k · f(y_{n−k})
フィードバックはy項 (線形フィードバック) と f(y) 項 (微分方程式を通じた非線形フィードバック) の両方に表示されます。IIRフィルタの安定性は、数値ODE積分器のより一般的な安定性問題の特定のケースです。
ドメイン全体でのフィードバックの接続
同じ数学的構造 — フィードバック、極、安定性の境界 — デジタルフィルタ、数値ODEソルバー、制御システム、生物学的リズム、& 経済モデルに表示されます。
各ドメインで: フィードバックループは前の状態から新しい状態を計算します。安定性には、フィードバックが無限に摂動を増幅しないようにする必要があります。
Z平面の単位円安定性の境界は、対応しています: ラプラスs平面の虚軸 (連続時間)、線形反復のスペクトル半径条件 ρ(A) < 1、& 非線形システムのリアプノフ指数条件 λ < 0。