バターワースフィルタ

**バターワースフィルタ**（）は、フィルタ回路設計の一種。通過帯域が数学的に可能な限り平坦な周波数特性となるよう設計されている。

バターワースフィルタは1930年、イギリスの技術者 が論文  で発表した。

また、特定のフィルタ回路構成を指す用語ではなく、フィルタの応答特性を指す用語であるため、**バターワースフィルタ特性**（あるいは**バターワース特性**）と呼ぶ場合もある。

## 概要
バターワースフィルタの周波数応答は通過帯域では最大限平坦であり（リップルがない）、除去帯域に向かってゼロに近づいていく。対数目盛のボーデ図で見ると、応答曲線は線形に負の無限大に近づいていく。一次フィルタの場合、応答曲線の傾斜は -6dB/octave または -20dB/decade となる（回路構成に関わらず、一次のバターワースフィルタは全てこの特性を示す）。二次バターワースフィルタの場合、応答曲線の傾斜は -12dB/octave、三次の場合 -18dB/octave となる。バターワースフィルタは、ωに対して振幅が単調に変化する。バターワースフィルタは高次になっても特性曲線が同じ形状（ただし、傾斜はきつくなる）だが、他のフィルタ（ベッセル、チェビシェフ、楕円など）は高次になると曲線の形状が変わる。

他のフィルタに比べるとバターワースフィルタによる減衰は緩やかであるため、特定の除去帯域仕様を実装するには高次な実装を必要とする。しかし、通過帯域は他のフィルタより線形な位相応答を示す。

## 簡単な例
[300px|right|thumb|三次ローパスフィルタ（Cauer形）。このフィルタが遮断周波数 &omega;<sub>c</sub>=1 のバターワースフィルタとなるのは、（例えば）C<sub>2</sub>=4/3 [ファラド|F、R4=1 Ω、L1=3/2 H、L3=1/2 H の場合である。]]

バターワースフィルタの簡単な例として三次ローパスフィルタを右図に示す。C_2=4/3F、R_4=1Ω、L_1=3/2H、L_3=1/2H とする。s=\sigma+j\omega は複素周波数とする。コンデンサ *C* のインピーダンスを *1/Cs*、コイル *L* のインピーダンスを *Ls* としたとき、この回路の伝達関数は以下のようになる。

:H(s)=\frac{V_{out}(s)}{V_{in}(s)}=\frac{1}{1+2s+2s^2+s^3}

周波数応答の大きさ（利得） G(\omega) は以下の式で得られる。

G^2(\omega)=|H(j\omega)|^2=\frac{1}{1+\omega^6}\,

また、位相は以下の式で得られる。

:\Phi(\omega)=\arg(H(j\omega))\,

256px|left|thumb|&omega;<sub>c</sub>=1 の三次バターワースフィルタの利得（緑）と群遅延（赤）

群遅延は、角周波数についての位相の微分と定義され、異なる複数の周波数間の位相差による信号の歪みの尺度である。このフィルタの利得と遅延をプロットしたものを左図に示す。利得曲線を見ると、通過帯域にも除去帯域にもリップルがないことがわかる。

伝達関数 *H(s)* の絶対値の対数を複素平面にプロットしたものが右図である。複素平面の左半分に3つの極がある。これらは単位円上にあり、実数軸を中心として対称に位置する。利得関数は右半分に3つの極を持ち、全体として単位円が完成する。

この回路のコイルとコンデンサを入れ替えると、ハイパス・バターワースフィルタとなる。うまく計算した値のコイルとコンデンサを並列接続したものをそれぞれの位置に入れると、バンドパス・バターワースフィルタになる。

## 伝達関数
thumb|350px|一次から五次までのバターワース特性のローパスフィルタの利得。''n''次のとき、傾斜は 20''n'' dB/decade になっている。

*n*次バターワース・ローパス・フィルタの利得 G(\omega) は、伝達関数 *H(s)* から次のように得られる。

:G^2(\omega)=\left |H(j\omega)\right|^2 = \frac {G_0^2}{1+\left(\frac{\omega}{\omega_c}\right)^{2n}} 

ここで
* n はフィルタの次数
* ωc は遮断周波数（約 -3dB となる周波数）
* G_0 はDC利得（ゼロ周波数での利得）
である。

*n* が無限大に近づくと、利得は矩形関数となり、ωc 以下の周波数は利得 G_0 で通過し、ωc 以上の周波数は抑止される。*n* が小さいほど、遮断は緩やかになる。

s=\sigma+j\omega の伝達関数 *H(s)* を決定することを考える。*s = jω* のときの *H(s)H(-s)* を計算すると |*H(jω)*|2 と同じになるため、次が得られる。

:H(s)H(-s) = \frac {G_0^2}{1+\left (\frac{-s^2}{\omega_c^2}\right)^n}

この式の極は半径 ωc の円上に等間隔で現れる。伝達関数自体は複素平面 *s* 上の実数が負の側の極で決定される。*k*番目の極は次の式で決定される。

:-\frac{s_k^2}{\omega_c^2} = (-1)^{\frac{1}{n}} = e^{\frac{j(2k-1)\pi}{n}}
\qquad\mathrm{k = 1,2,3, \ldots, n}

従って、次が得られる。

:s_k = \omega_c e^{\frac{j(2k+n-1)\pi}{2n}}\qquad\mathrm{k = 1,2,3, \ldots, n}

伝達関数はこれらの極を使って次のようにも表せる。

:H(s)=\frac{G_0}{\prod_{k=1}^n (s-s_k)/\omega_c}

分母は、*s* におけるバターワース多項式である。

### 正規化バターワース多項式
バターワース多項式は上記のように複素数形式でも書けるが、複素共役な極同士を掛け合わせることで実数形式で書くこともできる。この多項式は \omega_c=1 と設定することで正規化される。正規化バターワース多項式の一般形式は以下のようになる。

:B_n(s)=\prod_{k=1}^{\frac{n}{2}} \left[s^2-2s\cos\left(\frac{2k+n-1}{2n}\,\pi\right)+1\right] （n が偶数の場合）
:B_n(s)=(s+1)\prod_{k=1}^{\frac{n-1}{2}} \left[s^2-2s\cos\left(\frac{2k+n-1}{2n}\,\pi\right)+1\right] （n が奇数の場合）

小数点以下第4位まで表すと、以下のようになる。




### 最大平坦性
\omega_c=1 および G_0=1 としたとき、周波数毎の利得の導関数は以下の式となる。

:\frac{dG}{d\omega}=-nG^3\omega^{2n-1}

利得 *G* は常に正なので、全ての \omega について単調減少している。つまり、バターワースフィルタの利得関数にはリップルがない。さらに、利得を級数展開すると次のようになる。

:G(\omega)=1 - \frac{1}{2}\omega^{2n}+\frac{3}{8}\omega^{4n}+\ldots

言い換えれば、利得の導関数は2*n*次導関数を超えるまでゼロであり、それにより最大平坦性を生じる。

### 高周波ロールオフ
再度 \omega_c=1 としたとき、ω が大きいときの利得の対数の傾斜は以下のようになる。

:\lim_{\omega\rightarrow\infty}\frac{d\log(G)}{d\log(\omega)}=-n

デシベルで表すと、高周波ロールオフは 20*n* dB/decade または 6*n* dB/octave となる（電力は電圧利得の二乗に比例するため、20 という係数が使われる）。

## フィルタ設計
線形のアナログフィルタ回路の実装には、様々なトポロジーが存在する。あるトポロジーの回路は構成は同じだが、個々の部品の特性値が異なる。

### Cauer形

right|450px|Cauer形のバターワースフィルタ

Cauer形は受動部品だけで線形アナログフィルタを構成する。バターワースフィルタの伝達関数はCauer形の回路で実装できる。右図のk番目の部品の特性値は以下のようになる。

:C_k = 2 \sin \left [\frac {(2k-1)}{2n} \pi \right ] （k は奇数）

:L_k = 2 \sin \left [\frac {(2k-1)}{2n} \pi \right ] （k は偶数）

### Sallen-Key形
right|450px|Sallen-Key形（単にトポロジーを示しているだけで、個々の値はバターワースフィルタとは限らない）

Sallen-Key形は能動部品（オペアンプ）も使って線形アナログフィルタを実装するトポロジーである。各段のSallen-Key形回路で複素共役の2つの極を実装する。全体としては、Sallen-Key形の回路をカスケード接続してフィルタを構成する。*n* が奇数の場合実数の極ができるが、それは別途実装する必要があり、一般にRC回路で構成して、それをオペアンプの回路とカスケード接続する。

Sallen-Key形の伝達関数は次の通りである。

:H(s)=\frac{1}{1+C_2(R_1+R_2)s+C_1C_2R_1R_2s^2}

この分母がバターワース多項式の二次項の1つになるよう設定すればよい。\omega_c=1 とすると、

:C_1C_2R_1R_2=1\,

かつ

:C_2(R_1+R_2)=2\cos\left(\frac{2k+n-1}{2n} \pi\right)

となるように設定する。この場合、2つの部品の値が定まらないが、好きなように選べばよい。

### デジタルでの実装
バターワースフィルタをデジタルフィルタとして実装する場合、双一次変換やZ変換を使ってアナログフィルタを離散化することが多い。高次の場合、量子化誤差の影響が出やすくなる。そのためをカスケード接続したものとして計算することが多い（奇数次の場合は、一次フィルタもカスケード接続する）。

## 他の線形フィルタとの比較
下図は、離散時間バターワースフィルタと他のフィルタの利得を示したものである。いずれも五次のフィルタである。

500px|center

これらは同じ次数（この場合は五次）である。五次のフィルタであるとは、decade（周波数が10倍になる区間）当たり20dB×5 すなわち 100dB のロールオフとなることを意味する。バターワースフィルタは他のフィルタに比べて遮断周波数付近でのロールオフが緩やかだが、リップルが見られない。

## 参考文献

*