デジタルフィルタ設計と窓関数

理想フィルタは実現できない、から始める

デジタルフィルタ設計の出発点は、周波数領域で欲しい応答（理想の通過帯域と阻止帯域）を決め、それを時間領域のインパルス応答へ落とすことです。ところが理想的な「壁のように切り立った」特性は物理的に実現できません。本稿では、なぜ実現不能なのか、そのギャップを 窓関数法・周波数サンプリング法（FIR）と双一次変換（IIR） でどう埋めるのか、そして避けられない 通過帯域リプルと遷移幅のトレードオフ を原理から開封します。信号処理の土台は プログラミング の数値計算とも地続きです。

理想ローパスの周波数応答 H(e^jω) はカットオフ ωc まで 1、それ以上で 0 の矩形です。これを逆離散時間フーリエ変換すると、インパルス応答は次の sinc になります。

h[n] = (ωc / π) * sinc(ωc * n / π)
     = sin(ωc * n) / (π * n)      （n ≠ 0）
h[0] = ωc / π

この h[n] は n が正負の無限まで伸び、しかも非因果的（未来のサンプルを要する）です。実装するには 有限の長さに切り詰め、時間をずらして因果的にする 必要があり、ここで理想からのズレが生まれます。

窓関数法：打ち切りとギブス現象

最も素朴な方法は、無限長 h[n] を -M ≤ n ≤ M（長さ N = 2M+1）で単純に打ち切ることです。これは時間領域で矩形窓を掛けるのに等しく、周波数領域では 理想応答と矩形窓のスペクトル（ディリクレ核）の畳み込み になります。畳み込みの結果、カットオフ付近に振動が乗ります。

そこで矩形以外の、両端でなだらかに 0 へ減衰する窓 w[n] を掛けます。実際に実装する係数は次式です。

h_fir[n] = h[n - M] * w[n]     n = 0 .. N-1   （M だけ右シフトして因果化）

窓のスペクトルは中央の メインローブ と周囲の サイドローブ を持ちます。設計特性はこの2つで決まります。

• メインローブ幅 → 遷移帯域幅（通過から阻止へ移る急峻さ）を支配する。狭いほど急峻。
• サイドローブの高さ → 阻止帯域の減衰量と通過帯域リプルを支配する。低いほど深い阻止が得られる。

窓ごとのトレードオフを読む

窓の選択は「遷移帯を狭くする（メインローブを細く）」ことと「阻止量を深くする（サイドローブを低く）」ことの綱引きです。両立はできません。代表的な窓を比較します。

矩形は遷移帯こそ最も急峻ですが阻止量が浅く、ブラックマンは深い阻止と引き換えに遷移帯が広くなります。同じ阻止量なら、フィルタ長 N を伸ばすほど遷移帯は狭められる ため、実務では「必要な阻止量から窓を選び、必要な遷移幅から長さを決める」順で設計します。

周波数サンプリング法：欲しい応答を直接置く

窓関数法が「時間領域の理想応答を削る」のに対し、周波数サンプリング法 は「周波数領域で欲しい応答を等間隔に指定する」設計です。長さ N のフィルタに対し、N 点で所望の周波数応答 Hd[k]（k = 0 .. N-1）を並べ、逆離散フーリエ変換（IDFT）でインパルス応答を得ます。

h[n] = (1/N) * Σ_k  Hd[k] * exp( j * 2π * k * n / N )    k = 0 .. N-1

サンプリング点 ω = 2πk/N の上では応答が厳密に Hd[k] に一致します。問題は 点と点の間 です。理想的なステップをそのまま指定すると、サンプル間で大きなリプルが生じます。

FIR はインパルス応答が有限で、係数を左右対称（h[n] = h[N-1-n]）にすれば 厳密な直線位相（群遅延が全周波数で一定）を保証できます。これは波形の相対的な形を崩さない大きな利点で、画像や通信で重宝されます。ただし急峻な特性には長い N が要り、演算量とレイテンシが増えます。

IIR と双一次変換：アナログ設計を写す

急峻な特性を短い次数で得たいときは IIR（無限インパルス応答）を使います。定番は、成熟したアナログフィルタ（バターワース、チェビシェフ、楕円）の伝達関数 H(s) を、離散系の H(z) へ写す 双一次変換 です。連続系の s を次の写像で置き換えます。

s = (2/T) * (z - 1) / (z + 1)       T: サンプリング周期

この写像は s 平面の左半面（安定極）を z 平面の単位円内へ1対1で写すため、元のアナログフィルタが安定なら離散化後も必ず安定 です。これは単純な微分の差分近似（前進差分）にはない、双一次変換の決定的な長所です。

FIR と IIR、どちらを選ぶか

同じ仕様なら IIR の方が圧倒的に低次数（少ない乗算・少ないメモリ）で急峻な特性を実現できます。代償は、位相特性が一般に非線形（群遅延が周波数依存）で波形がひずむこと、そしてフィードバックを持つため係数の量子化やオーバーフローで 不安定化しうる ことです。FIR は必ず安定し直線位相を作れる代わりに、同じ急峻さに高い次数が必要です。

まとめ

デジタルフィルタ設計は「実現不能な理想を、制御されたトレードオフで近似する」営みです。要点は、(1) 理想 sinc は無限長・非因果なので打ち切りと因果化が要り、矩形打ち切りはギブス現象で約9%のリプルを残すこと、(2) 窓関数法ではメインローブ幅（遷移帯）とサイドローブ高（阻止量・リプル）が綱引きし、カイザー窓は β でこれを連続調整できること、(3) 周波数サンプリング法は欲しい応答を直接置き、遷移サンプルでリプルを抑えること、(4) IIR は双一次変換でアナログ設計を安定に写せるが、tan による周波数ワーピングをプリワープで補正する必要があること、です。急峻さ・位相・演算量・安定性の要求を天秤にかけ、FIR と IIR を使い分けてください。