目次

キーワード

概要

バターワースフィルタ(Butterworth filter)は、 通過域で最大平坦な振幅特性を示すローパスフィルタです。 位相特性も線形に近いという特徴があります。

バターワースフィルタの周波数特性
バターワースフィルタの周波数特性

周波数特性

バターワースフィルタは以下のような周波数特性を持ちます。

|H(ω)|2
1
1 + ω2n

この特性はω>1で急激に減衰していきます。 ωの変わりにω/ωcを代入することでカットオフ周波数ωcを任意に設定できます。 また、この伝達関数は1階から2n-1階までの全ての導関数がω = 0において0であるという性質を持っています。 この性質を最大平坦と呼びます。

図2に、例として、3次、5次、9次のバターワースフィルタの振幅特性を示します。

3~9次のバターワースフィルタの周波数特性
3~9次のバターワースフィルタの周波数特性

アナログプロトタイプの設計

極配置

|H(ω)|2の分母は2n次の多項式なので、 |H(ω)|22n個の極を持ちます。 この2n個の極は、 ω2n = -1より、以下のようになります。

ω = exp(i
2k + 1
2n
π)

この2n個の中から安定なn個の極を選び、 バターワースフィルタの極 s にします。

s = - cos
n - 2k - 1
2n
π ± i sin
n - 2k - 1
2n
π
ただし、k0~(n-1) / 2までの整数です。

アナログプロトタイプ伝達関数

決定された極配置から、バターワースフィルタの伝達関数H(s)は以下のようになります。

H(s)
n
2
- 1
k=0
1
s2 + αk s + 1
   (nが偶数のとき)
H(s)
1
s + 1
n-1
2
k=0
1
s2 + αk s + 1
   (nが奇数のとき)

ただし、αk は以下の式で表される値です。

αk = 2 cos(
n - 2 k - 1
2 n
π)

設計仕様

透過域/阻止域の周波数/リプル (Ap, rs, ωp, ωs) を仕様として与えたとき、 仕様を満たす最小の次数 N を求める方法について説明します。

バターワースフィルタの設計仕様
バターワースフィルタの設計仕様

このような次数 N を求めるためには、 図3に示すように、

Ap2
1
1 + ωp2 N
,   rp2
1
1 + ωs2 N

という条件を満たすような最小の N を求めます。 これらの式から、

N ≧
log(
1
Ap2
- 1 )
2 log ωp
,   N ≧
log(
1
rs2
- 1 )
2 log ωs

という式が得られるので、 これらの最大値を選ぶことにより、

N ≧ Max(
log(
1
Ap2
- 1 )
2 log ωp
,

log(
1
rs2
- 1 )
2 log ωs

)

という式が得られます。 この不等式を満たすような最小の N を選びます。

ディジタルフィルタ設計

アナログプロトタイプの設計」で得られた式を、 カットオフ周波数ωcで双1次変換することで以下の式が得られます。

c = cos ωc , s = sin ωc
H(z)
n
2
-1
k=0
(
1-c
2
)( 1 + 2 z-1 + z-2)
(1+αk, n) - 2 c z-1(1-αk, n) z-2
,    (nが偶数のとき)
H(z)
s (1 + z-1)
(s + c + 1)(s - c - 1) z-1
×
n-1
2
k=0
(
1-c
2
)( 1 + 2 z-1 + z-2)
(1+αk, n) - 2 c z-1(1-αk, n) z-2
   (nが奇数のとき)

ただし、αk, n は以下の式で表される値です。

αk, n = s cos(
n - 2k + 1
2n
π)

ハイパスフィルタ

周波数変換」で説明したように、 ローパスフィルタのアナログ伝達関数の変数 ss →

1
s
に入れ替えると、 ハイパスフィルタになります。 したがって、 ハイパスバターワースフィルタのアナログプロトタイプ伝達関数は以下のようになります。

H(s)
n
2
- 1
k=0
s2
s2 + αk s + 1
   (nが偶数のとき)
H(s)
s
s + 1
n-1
2
k=0
s2
s2 + αk s + 1
   (nが奇数のとき)

この式の分母は、ローパス伝達関数の分母と全く同じものになっています。 (バターワースフィルタの伝達関数は、分母の定数項と2次の係数が同じなのでこうなる。) その結果、ディジタルフィルタ伝達関数も、以下のように、 ローパス伝達関数と分母が全く同じものになります。

H(z)
n
2
-1
k=0
(
1+c
2
)( 1 - 2 z-1 + z-2)
(1+αk, n) - 2 c z-1(1-αk, n) z-2
   (nが偶数のとき)
H(z)
(1c)(1 z-1)
(s + c + 1)(s - c - 1) z-1
×
n-1
2
k=0
(
1c
2
)( 1 2 z-1 + z-2)
(1+αk, n) - 2 c z-1(1-αk, n) z-2
   (nが奇数のとき)

この式を、ローパス伝達関数と見比べてみると、 分母が同じなだけでなく、 分子も似たものになっています。 式中で強調表示している部分の符号が異なっているだけです。 このことから、同じカットオフ周波数のローパス・ハイパスフィルタは、 その処理の大部分を共通化することができます。

更新履歴

ブログ