双一次変換のソースを表示

'''双一次変換'''（そういちじへんかん Bilinear transform 双一次Z変換 タスティン変換、台形差分法 Trapezoidal methodとも呼ばれる）は、[[デジタル信号処理]]において、[[連続]][[時間領域]]における[[線型性|線型]][[時不変]] ([[線形時不変システム|LTI]])フィルタの[[伝達関数法|伝達関数]]<math> H_a(s) \ </math>([[アナログフィルタ]]と呼ばれる)を[[離散信号|離散]]時間領域における線形シフト不変フィルタの伝達関数<math> H_d(z) \ </math>([[スイッチトキャパシタ]]で構成されるアナログフィルタも離散時間フィルタだが、[[デジタルフィルタ]]と呼ばれる)に変換するのによく用いられる[[等角写像]]のひとつである。
この変換では、[[s平面]]上の<math> Re[s]=0 \ </math>、<math> j \omega \ </math>を[[z平面]]上の<math> |z| = 1 \ </math>の[[単位円]]に写像する。
双一次変換は元のフィルタの安定性を保存し、連続時間フィルタ<math> H_a(j \omega_a) \ </math>の[[周波数応答]]のすべての点を離散時間フィルタ<math> H_d(e^{j \omega_d T}) \ </math>の周波数応答の対応する点に1対1に写像する。ただし以下の周波数歪みの項で述べるようにもとの周波数とはすこし異なる周波数へ写像される。この歪みは低周波域ではほとんど目立たないが、[[ナイキスト周波数]]に近づくほど顕著になる。

単位インパルスによってサンプリングされた離散時間信号に[[ラプラス変換]]を行うと結果は正確に[[Z変換]]として表されるが、双一次変換は以下の式に示すように[[z平面]]から[[s平面]]への正確な写像を行う自然対数関数の一次近似である。

{{Indent|<math>
\begin{align}
z &= e^{sT}   \\
  &= \frac{e^{sT/2}}{e^{-sT/2}} \\
  &\approx \frac{1 + s T / 2}{1 - s T / 2}
\end{align}
</math>}}
（ただし<math> T \ </math>は離散時間フィルタのサンプル時間で[[サンプリング周波数]]の逆数）により近似できる。上の式を<math> s \ </math>について解くか、<math> s = (1/T) \ln(z) \  \ </math>を同様に近似すると、この変換の逆写像とその双一次近似は
{{Indent|<math>
\begin{align}
s &= \frac{1}{T} \ln(z)  \\
  &= \frac{2}{T} \tanh^{-1} \frac{z-1}{z+1} \\
  &= \frac{2}{T} \left[\frac{z-1}{z+1} + \frac{1}{3} \left( \frac{z-1}{z+1} \right)^3  + \frac{1}{5} \left( \frac{z-1}{z+1} \right)^5  + \frac{1}{7} \left( \frac{z-1}{z+1} \right)^7 + \cdots \right] \\
  &\approx  \frac{2}{T} \frac{z - 1}{z + 1} \\
  &\approx  \frac{2}{T} \frac{1 - z^{-1}}{1 + z^{-1}}
\end{align}
</math>}}
となる。双一次変換とは、この一次近似を用い連続時間の伝達関数<math> H_a(s) \ </math>において
{{Indent|<math>s \leftarrow \frac{2}{T} \frac{z - 1}{z + 1}</math>とし<br />
<math>H_d(z) = H_a(s) \bigg|_{s = \frac{2}{T} \frac{z - 1}{z + 1}}= H_a \left( \frac{2}{T} \frac{z-1}{z+1} \right)\ </math>とするものである。}}

双一次変換は[[メビウス変換]]と呼ばれる[[等角写像]]の特殊な場合であり、以下のように定義される。
{{Indent|<math>z^{\prime} = \frac{a z + b}{c z + d}</math>}}

== 安定性と最小位相の保持 ==
連続時間フィルタはその伝達関数の極が[[s平面]]上の左半平面にあるとき、[[BIBO安定|安定]]である。離散時間フィルタはその伝達関数の極が[[z平面]]上の[[単位円]]内にあるとき安定である。双一次変換はs平面の左半平面をz平面上の単位円内に写像する。このため連続時間領域で設計されたフィルタが安定ならそれを変換した離散時間領域のフィルタも安定である。

同様に、連続時間フィルタが[[最小位相]]ならその伝達関数の[[零点]]はs平面上の左半平面にあるため、離散時間フィルタの伝達関数の零点はz平面上の単位円内に位置し、最小位相となる。

== 周波数歪み ==
連続時間フィルタの周波数応答は、[[伝達関数]]<math> H_a(s) \ </math>で<math>s = j \omega \ </math>とすれば求められる。同様に離散時間フィルタの周波数応答は伝達関数<math> H_d(z) \ </math>で{{font|<math>z = e^{ j \omega T} \ </math>}}とすれば求められる。双一次変換で設計された離散時間フィルタに実際に<math> \omega \ </math>が入力されたとき、その周波数応答が連続時間フィルタのどの周波数<math> \omega_a \ </math>に対応するかというと、
{{Indent|<math>H_d(z) = H_a \left( \frac{2}{T} \frac{z-1}{z+1}\right) \ </math>}}
{|
|-
|<math>H_d(e^{ j \omega T}) \ </math>
|<math>= H_a \left( \frac{2}{T} \frac{e^{ j \omega T} - 1}{e^{ j \omega T} + 1}\right) \ </math>
|-
|
|<math>= H_a \left( \frac{2}{T} \cdot \frac{e^{j \omega T/2} \left(e^{j \omega T/2} - e^{-j \omega T/2}\right)}{e^{j \omega T/2} \left(e^{j \omega T/2} + e^{-j \omega T/2 }\right)}\right) \ </math>
|-
|
|<math>= H_a \left( \frac{2}{T} \cdot \frac{\left(e^{j \omega T/2} - e^{-j \omega T/2}\right)}{\left(e^{j \omega T/2} + e^{-j \omega T/2 }\right)}\right) \ </math>
|-
|
|<math>= H_a \left(j \frac{2}{T} \cdot \frac{ \left(e^{j \omega T/2} - e^{-j \omega T/2}\right) /(2j)}{\left(e^{j \omega T/2} + e^{-j \omega T/2 }\right) / 2}\right) \ </math>
|-
|
|<math>= H_a \left(j \frac{2}{T} \cdot \frac{ \sin(\omega T/2) }{ \cos(\omega T/2) }\right) \ </math>
|-
|
|<math>= H_a \left(j \frac{2}{T} \cdot \tan \left( \omega \frac{T}{2} \right) \right) \ </math>
|-
|
|<math>= H_a \left(j \omega_a \right) \ </math>
|}

となり、これは離散時間フィルタにおいて[[z平面]]内の単位円上のすべての点、{{font|<math>z = e^{ j \omega T} \ </math>|css=vertical-align:+30%;}}が連続時間フィルタの[[s平面]]上の<math>j \omega \ </math>軸、<math>s = j \omega_a \ </math>に写像されることを示している。よって双一次変換での離散時間周波数から連続時間周波数への写像は

{{Indent|<math> \omega_a = \frac{2}{T} \tan \left( \omega \frac{T}{2} \right) </math>}}
となり、その逆は
{{Indent|<math> \omega = \frac{2}{T} \arctan \left( \omega_a \frac{T}{2} \right) </math>}}
となる。

離散時間フィルタは周波数<math>\omega \ </math>において連続時間フィルタの周波数<math> (2/T) \tan(\omega T/2) \ </math>での振る舞いと同じ振る舞いをする。特にゲインと位相に関して、離散時間フィルタは周波数<math>\omega \ </math>において連続時間フィルタの周波数<math> (2/T) \tan(\omega T/2) \ </math>でのゲイン・位相に等しくなる。その特性の現れる周波数はわずかに異なるが、連続時間フィルタの周波数応答のすべての特徴が離散時間フィルタに現れることを意味する。低周波域においては<math>\omega \approx \omega_a \ </math>となる。

連続時間周波数の範囲
{{Indent|<math> -\infty < \omega_a < +\infty \ </math>}}
は周波数区間
{{Indent|<math> -\frac{\pi}{T} < \omega < +\frac{\pi}{T}. \ </math>}}
に写像される。
連続時間フィルタの周波数<math> \omega_a = 0 \ </math>なら対応する離散時間フィルタの周波数<math> \omega = 0 \ </math>となり、連続時間フィルタの周波数<math> \omega_a = \pm \infty \ </math>なら対応する離散時間フィルタの周波数<math> \omega = \pm \pi / T\ </math>となる。

さらに<math> \omega_a \ </math> と <math> \omega \ </math>の関係は非線形であり、これは'''''周波数歪み'''''と呼ばれている。
連続時間フィルタの仕様として与えられている周波数（[[遮断周波数]]や[[中心周波数]]）を、この
<math> \omega_a = \frac{2}{T} \tan \left( \omega \frac{T}{2} \right) \ </math>
によってあらかじめ補正して設計することもでき、これは'''''プリワーピング'''''と呼ばれる。
この周波数歪みによる主な利点は、[[インパルス不変法]]で見られるような周波数応答の[[エイリアシング]]が発生しないことである。

== 関連項目 ==
* [[デジタル信号処理]]
* [[数値積分]]
* [[台形則]]
* [[前進差分法]]
* [[後退差分法]]

{{デジタル信号処理}}

{{DEFAULTSORT:そういちしへんかん}}
[[Category:デジタル回路]]
[[Category:信号処理]]