双一次変換シミュレーター — アナログからデジタルIIRへ

オーディオ機器のイコライザー：ミキサーやプラグインのEQ、トーンコントロール、クロスオーバーネットワークは、ほとんどがIIRフィルタで実装されます。設計者はまずアナログのバターワースやベル型シェルビングフィルタを s 平面で設計し、双一次変換でデジタル係数に落とします。オーディオはサンプリング周波数が 44.1 kHz や 48 kHz と固定なので、プリワープでカットオフ点を正確に合わせることが音作りの精度に直結します。

センサー信号のノイズ除去：加速度センサー、ひずみゲージ、温度センサーなどの計測データには高周波ノイズが乗ります。マイコンやDSP上で1次〜2次のIIRローパスを走らせて平滑化するのが定番で、本ツールの差分方程式 y[n] = b0·x[n] + b1·x[n−1] − a1·y[n−1] がそのまま組み込みコードになります。FIRフィルタより少ない係数・少ない演算で同等の遮断特性が得られるため、リソースの限られた組み込み機器で重宝します。

制御システムの離散化：PID制御器やリードラグ補償器をアナログで設計し、デジタルコントローラに実装する際にも双一次変換（タスティン変換）が使われます。連続時間の伝達関数をそのままサンプリング周期 T で離散化でき、安定性が保存されるため、アナログ設計の安定余裕がデジタル実装でも崩れません。サーボモーターやインバータのファームウェアで広く使われています。

通信・生体信号処理：無線受信機のチャネル選択フィルタ、心電図（ECG）の基線変動除去、脳波（EEG）の帯域分割など、リアルタイム性が求められる分野でもIIRフィルタが使われます。アナログプロトタイプフィルタ（楕円フィルタなど）を双一次変換で離散化することで、急峻な遮断特性を少ない次数で実現できます。

まず最大の落とし穴が、「プリワープを忘れてカットオフがズレる」ことです。アナログのカットオフ ωc をそのまま K = ωc·T/2 に入れて係数を計算すると、できあがったデジタルフィルタの −3 dB 点は目標の fc より低くなります。低い周波数・高いサンプリング周波数ではズレが小さく見過ごされがちですが、カットオフがナイキスト周波数に近いほどズレは無視できません。本ツールの「周波数ワーピング率」は、プリワープ後 ωc と元の 2π·fc の比です。この値が 1 から離れるほど、プリワープしないと設計どおりにならない、という危険信号です。

次に、「双一次変換は振幅も位相も完全に保存する」と思い込むこと。双一次変換が保存するのは「振幅応答の形」と「安定性」であって、周波数軸そのものは非線形にワーピングされます。つまりアナログとデジタルで、ある周波数での利得は一致させられても、周波数軸の目盛りが伸び縮みしています。特に位相応答（群遅延）はアナログプロトタイプと同じにはなりません。線形位相が必要な用途では、IIRではなくFIRフィルタを検討すべきです。双一次変換はあくまで「振幅特性とエイリアシング回避」を優先する手法だと理解してください。

最後に、「係数を計算して終わり」ではないという点。本ツールは理想的な実数演算で b0・b1・a1 を求めますが、実際の組み込み実装では係数を有限ビットで量子化します。特に固定小数点DSPでは、係数の丸め誤差で極が単位円の外へわずかにはみ出し、安定なはずのフィルタが発振することがあります。高次フィルタを1つの伝達関数で実装せず、2次セクション（バイクワッド）の縦続接続に分解するのは、この量子化感度を下げるためです。双一次変換で得た係数は「設計値」であり、実装時には量子化と数値精度の検証が別途必要です。