双一次変換

双一次変換（そういちじへんかん Bilinear transform 双一次Z変換 タスティン変換、台形差分法 Trapezoidal methodとも呼ばれる）は、デジタル信号処理において、連続時間領域における線型時不変 (LTI)フィルタの伝達関数formula_1(アナログフィルタと呼ばれる)を離散時間領域における線形シフト不変フィルタの伝達関数formula_2(スイッチトキャパシタで構成されるアナログフィルタも離散時間フィルタだが、デジタルフィルタと呼ばれる)に変換するのによく用いられる等角写像のひとつである。
この変換では、s平面上のformula_3、formula_4をz平面上のformula_5の単位円に写像する。
双一次変換は元のフィルタの安定性を保存し、連続時間フィルタformula_6の周波数応答のすべての点を離散時間フィルタformula_7の周波数応答の対応する点に1対1に写像する。ただし以下の周波数歪みの項で述べるようにもとの周波数とはすこし異なる周波数へ写像される。この歪みは低周波域ではほとんど目立たないが、ナイキスト周波数に近づくほど顕著になる。
単位インパルスによってサンプリングされた離散時間信号にラプラス変換を行うと結果は正確にZ変換として表されるが、双一次変換は以下の式に示すようにz平面からs平面への正確な写像を行う自然対数関数の一次近似である。
（ただしformula_8は離散時間フィルタのサンプル時間でサンプリング周波数の逆数）により近似できる。上の式をformula_9について解くか、formula_10を同様に近似すると、この変換の逆写像とその双一次近似は
となる。双一次変換とは、この一次近似を用い連続時間の伝達関数formula_1において
双一次変換はメビウス変換と呼ばれる等角写像の特殊な場合であり、以下のように定義される。
安定性と最小位相の保持.
連続時間フィルタはその伝達関数の極がs平面上の左半平面にあるとき、安定である。離散時間フィルタはその伝達関数の極がz平面上の単位円内にあるとき安定である。双一次変換はs平面の左半平面をz平面上の単位円内に写像する。このため連続時間領域で設計されたフィルタが安定ならそれを変換した離散時間領域のフィルタも安定である。
同様に、連続時間フィルタが最小位相ならその伝達関数の零点はs平面上の左半平面にあるため、離散時間フィルタの伝達関数の零点はz平面上の単位円内に位置し、最小位相となる。
周波数歪み.
連続時間フィルタの周波数応答は、伝達関数formula_1でformula_13とすれば求められる。同様に離散時間フィルタの周波数応答は伝達関数formula_2でとすれば求められる。双一次変換で設計された離散時間フィルタに実際にformula_15が入力されたとき、その周波数応答が連続時間フィルタのどの周波数formula_16に対応するかというと、
となり、これは離散時間フィルタにおいてz平面内の単位円上のすべての点、が連続時間フィルタのs平面上のformula_17軸、formula_18に写像されることを示している。よって双一次変換での離散時間周波数から連続時間周波数への写像は
となり、その逆は
となる。
離散時間フィルタは周波数formula_19において連続時間フィルタの周波数formula_20での振る舞いと同じ振る舞いをする。特にゲインと位相に関して、離散時間フィルタは周波数formula_19において連続時間フィルタの周波数formula_20でのゲイン・位相に等しくなる。その特性の現れる周波数はわずかに異なるが、連続時間フィルタの周波数応答のすべての特徴が離散時間フィルタに現れることを意味する。低周波域においてはformula_23となる。
連続時間周波数の範囲
は周波数区間
に写像される。
連続時間フィルタの周波数formula_24なら対応する離散時間フィルタの周波数formula_25となり、連続時間フィルタの周波数formula_26なら対応する離散時間フィルタの周波数formula_27となる。
さらにformula_16 と formula_15の関係は非線形であり、これは周波数歪みと呼ばれている。
連続時間フィルタの仕様として与えられている周波数（遮断周波数や中心周波数）を、この
formula_30
によってあらかじめ補正して設計することもでき、これはプリワーピングと呼ばれる。
この周波数歪みによる主な利点は、インパルス不変法で見られるような周波数応答のエイリアシングが発生しないことである。