Type of signal processing filter
バターワースの1930年の論文の周波数応答プロット。 [1] バターワース フィルタは、 通過帯域 において可能な限り平坦な 周波数応答 を持つように設計された 信号処理フィルタ の一種です。 最大平坦振幅フィルタ とも呼ばれます。1930年にイギリスの技術者であり物理学者でもあった スティーブン・バターワース が、論文「フィルタ増幅器の理論について」で 初めて記述しました。 [1]
原著論文 バターワースは、「不可能」と思われていた非常に複雑な数学的問題を解くことで名声を博していました。当時、 フィルタ設計には、 当時使用されていた理論 の限界により、相当な設計経験が必要でした 。このフィルタは、発表後30年以上もの間、広く使用されることはありませんでした。バターワースは次のように述べています。
「理想的な電気フィルタは、不要な周波数を完全に除去するだけでなく、必要な周波数に対して均一な感度も持つ必要があります」。
このような理想的なフィルタは実現不可能ですが、バターワースは、適切な値のフィルタ要素の数を増やすことで、より近似値が得られることを示しました。当時のフィルタは通過帯域に大きなリップルを生成し、部品値の選択は大きく影響していました。バターワースは、周波数の関数としてゲイン(つまり、 周波数応答 の大きさ )が
次のように表さ れる ローパスフィルタ を設計できることを示しました。
G ( ω ) = 1 1 + ω 2 n , {\displaystyle G(\omega )={\frac {1}{\sqrt {1+{\omega }^{2n}}}},} ここで 、 は ラジアン/秒で表した 角周波数 、 はフィルタの 極 数で、パッシブフィルタのリアクタンス素子の数に等しい。 カットオフ周波数 ( 電力の半分の点 、およそ -3 dB または電圧ゲイン 1/ √2 ≈ 0.7071)は、𝜔 = 1 ラジアン/秒に正規化されている。 Butterworth は論文では偶数極のフィルタのみを扱っているが、奇数次フィルタは偶数次フィルタの出力に単極フィルタを追加することで作成できる。彼は高次フィルタを、真空管アンプで区切られた 2 極フィルタから構築した。 2、4、6、8、10 極フィルタの周波数応答のプロットは、彼の元のグラフで A、B、C、D、E として示されている。 ω {\displaystyle \omega } n {\displaystyle n}
バターワースは2極フィルタと4極フィルタの方程式を解き、後者は 真空管 アンプで分離することでカスケード接続が可能になり、 インダクタ 損失にもかかわらず高次フィルタを構成できることを示しました。1930年当時、 モリパーマロイ のような低損失コア材料は まだ発見されておらず、空芯オーディオインダクタはむしろ損失が大きいものでした。バターワースは、フィルタの部品定数を調整することでインダクタの巻線抵抗を補正できることを発見しました。
彼は直径1.25インチ、長さ3インチのプラグイン端子付きコイルフォームを使用しました。巻かれたコイルフォーム内には、関連するコンデンサと抵抗器が内蔵されていました。コイルはプレート負荷抵抗器の一部を形成していました。真空管1本につき2極が使用され、後続の真空管のグリッドとはRC結合されていました。
バターワースは、基本的なローパス フィルタを変更して、 ローパス 、 ハイパス 、 バンドパス 、 バンドストップの 機能を持たせることができることも示しました。
概要 一次ローパスフィルタの ボード 線図 バターワースフィルタの周波数応答は、 通過帯域で最大限に平坦(つまり、 リップルがない)で、 阻止帯域 でゼロに向かってロールオフします。 [2]対数 ボード線図 で見ると 、応答は負の無限大に向かって直線的に傾斜しています。1 次フィルタの応答は、 オクターブあたり −6 dB(10 倍 あたり−20 dB )でロールオフします(すべての 1 次ローパスフィルタは同じ正規化された周波数応答を持ちます)。2 次フィルタはオクターブあたり −12 dB で減少し、3 次フィルタは −18 dB で減少します。バターワースフィルタは、 通過帯域や阻止帯域で非単調なリップルを持つ他のフィルタタイプとは異なり、 で単調に変化する振幅関数を持ちます。 ω {\displaystyle \omega }
チェビシェフ I 型/II 型フィルタや 楕円フィルタ と比較すると 、バターワース フィルタは ロールオフ が遅いため、特定の ストップバンド 仕様を実装するにはより高い次数が必要になりますが、バターワース フィルタは、チェビシェフ I 型/II 型フィルタや楕円フィルタが実現できるものよりも、通過帯域においてより線形な位相応答を実現します。
例 右の図に示す 3 次ローパス バターワース フィルタ設計の伝達関数は次のようになります。
V o ( s ) V i ( s ) = R 4 s 3 ( L 1 C 2 L 3 ) + s 2 ( L 1 C 2 R 4 ) + s ( L 1 + L 3 ) + R 4 {\displaystyle {\frac {V_{o}(s)}{V_{i}(s)}}={\frac {R_{4}}{s^{3}(L_{1}C_{2}L_{3})+s^{2}(L_{1}C_{2}R_{4})+s(L_{1}+L_{3})+R_{4}}}} 3次ローパスフィルタ( カウアートポロジー )。例えば、 =4/3 F、 =1 Ω、 =3/2 H 、 =1/ 2 Hのとき、このフィルタは カットオフ周波数 =1のバターワースフィルタになります。 ω c {\displaystyle \omega _{c}} C 2 {\displaystyle C_{2}} R 4 {\displaystyle R_{4}} L 1 {\displaystyle L_{1}} L 3 {\displaystyle L_{3}} バターワースフィルタの簡単な例として、右の図に示す3次ローパス設計があります。ここで、 = 4/3 F、 = 1 Ω、 = 3/2 H、 = 1/2 Hです。 [3] コンデンサの インピーダンス を 、インダクタのインピーダンスを (複素周波数)とすると、回路方程式 から この デバイスの 伝達関数 が得られます。 C 2 {\displaystyle C_{2}} R 4 {\displaystyle R_{4}} L 1 {\displaystyle L_{1}} L 3 {\displaystyle L_{3}} C {\displaystyle C} 1 / ( C s ) {\displaystyle 1/(Cs)} L {\displaystyle L} L s {\displaystyle Ls} s = σ + j ω {\displaystyle s=\sigma +j\omega }
H ( s ) = V o ( s ) V i ( s ) = 1 1 + 2 s + 2 s 2 + s 3 . {\displaystyle H(s)={\frac {V_{o}(s)}{V_{i}(s)}}={\frac {1}{1+2s+2s^{2}+s^{3}}}.} 周波数応答(ゲイン)の大きさは 次のように表される。 G ( ω ) {\displaystyle G(\omega )}
G ( ω ) = | H ( j ω ) | = 1 1 + ω 6 , {\displaystyle G(\omega )=|H(j\omega )|={\frac {1}{\sqrt {1+\omega ^{6}}}},} から得た
G 2 ( ω ) = | H ( j ω ) | 2 = H ( j ω ) ⋅ H ∗ ( j ω ) = 1 1 + ω 6 , {\displaystyle G^{2}(\omega )=|H(j\omega )|^{2}=H(j\omega )\cdot H^{*}(j\omega )={\frac {1}{1+\omega ^{6}}},} そして 位相 は次のように与えられる。
Φ ( ω ) = arg ( H ( j ω ) ) . {\displaystyle \Phi (\omega )=\arg(H(j\omega )).\!} 3次バターワースフィルタの ゲインと 群遅延は、 ω c = 1 {\displaystyle \omega _{c}=1} 群遅延 は 、位相シフトの角周波数に対する負の微分として定義され、異なる周波数における位相差によって生じる信号歪みの尺度となります。このフィルタのゲインと遅延は、左側のグラフに示されています。通過帯域と阻止帯域のどちらにおいても、ゲイン曲線にリップルはありません。
右側の2番目のグラフは、 伝達関数の絶対値の対数を複素周波数空間にプロットしたものです。この関数は、複素周波数平面の左半分にある3つの極によって定義されます。 H ( s ) {\displaystyle H(s)}
3次バターワースフィルタ(=1)の 複素周波数空間 における 伝達関数の対数密度プロット。3 つの 極は 左半平面上の単位半径の円上に位置している。 H ( s ) {\displaystyle H(s)} ω c {\displaystyle \omega _{c}} これらは、実軸を中心に対称な 半径1の円 上に配置されます 。ゲイン関数は、円を完成させるために右半平面上にさらに3つの極を持ちます。 s {\displaystyle s}
各インダクタをコンデンサに置き換え、各コンデンサをインダクタに置き換えることで、ハイパス バターワース フィルタが得られます。
バンドパス・バターワースフィルタは、各インダクタに直列にコンデンサ、並列にインダクタを配置して共振回路を形成することで実現されます。新しい部品の値は、対象となる周波数において古い部品と共振するように選択する必要があります。
バンドストップ・バターワースフィルタは、各インダクタに並列にコンデンサ、各コンデンサに直列にインダクタを配置して共振回路を形成することで実現されます。新しい各部品の値は、除去したい周波数において古い部品と共振するように選択する必要があります。
伝達関数 カットオフ周波数 が である1次から5次のバターワースローパスフィルタのゲインのプロット 。傾きは20 dB/decadeであることに注意する。ここで はフィルタ次数である。 ω c = 1 {\displaystyle \omega _{c}=1} n {\displaystyle n} n {\displaystyle n} すべてのフィルタと同様に、典型的な プロトタイプ はローパス フィルタです。ローパス フィルタはハイパス フィルタに変更することも、他のフィルタと直列に接続して バンドパス フィルタや バンドストップ フィルタ、さらにこれらの高次バージョンを形成することもできます。
次バターワースローパスフィルタ の ゲインは、伝達関数で 次のよう
に表される。 G ( ω ) {\displaystyle G(\omega )} n {\displaystyle n} H ( s ) {\displaystyle H(s)}
G 2 ( ω ) = | H ( j ω ) | 2 = G 0 2 1 + ( ω ω c ) 2 n {\displaystyle G^{2}(\omega )=\left|H(j\omega )\right|^{2}={\frac {{G_{0}}^{2}}{1+\left({\frac {\omega }{\omega _{c}}}\right)^{2n}}}} ここで 、はフィルタの次数、 は カットオフ周波数 (およそ -3 dB 周波数)、は DC ゲイン(ゼロ周波数でのゲイン)です。 n {\displaystyle n} ω c {\displaystyle \omega _{c}} G 0 {\displaystyle G_{0}}
が無限大に近づくにつれて、ゲインは矩形関数となり、それ以下の周波数は ゲイン で通過し 、それを超える周波数は 抑制されること がわかります。 の値が小さいほど 、カットオフは緩やかになります。 n {\displaystyle n} ω c {\displaystyle \omega _{c}} G 0 {\displaystyle G_{0}} ω c {\displaystyle \omega _{c}} n {\displaystyle n}
伝達関数を( ラプラス変換 より)で 決定したい 。 および は、 におけるラプラス変換の一般的な性質として 、であるため 、 とすると、 となる 。 H ( s ) {\displaystyle H(s)} s = σ + j ω {\displaystyle s=\sigma +j\omega } | H ( s ) | 2 = H ( s ) H ( s ) ¯ {\displaystyle \left|H(s)\right|^{2}=H(s){\overline {H(s)}}} s = j ω {\displaystyle s=j\omega } H ( − j ω ) = H ( j ω ) ¯ {\displaystyle H(-j\omega )={\overline {H(j\omega )}}} H ( s ) {\displaystyle H(s)}
H ( s ) H ( − s ) = G 0 2 1 + ( − s 2 ω c 2 ) n , {\displaystyle H(s)H(-s)={\frac {{G_{0}}^{2}}{1+\left({\frac {-s^{2}}{\omega _{c}^{2}}}\right)^{n}}},} すると、 バターワースフィルタの周波数応答が得られます。 s = j ω {\displaystyle s=j\omega }
この式の極は、半径 の円周上の等間隔の点に現れ 、 負の実軸の周りに対称である。 したがって、安定性のために、伝達関数 は、 の負の実半平面内の極のみを含むように選択される 。第 極は次のように指定される
。 n {\displaystyle n} ω c {\displaystyle \omega _{c}} H ( s ) {\displaystyle H(s)} s {\displaystyle s} k {\displaystyle k}
− s k 2 ω c 2 = ( − 1 ) 1 n = e j ( 2 k − 1 ) π n k = 1 , 2 , 3 , … , n {\displaystyle -{\frac {s_{k}^{2}}{\omega _{c}^{2}}}=(-1)^{\frac {1}{n}}=e^{\frac {j(2k-1)\pi }{n}}\qquad k=1,2,3,\ldots ,n} そしてそれゆえ
s k = ω c e j ( 2 k + n − 1 ) π 2 n k = 1 , 2 , 3 , … , n . {\displaystyle s_{k}=\omega _{c}e^{\frac {j(2k+n-1)\pi }{2n}}\qquad k=1,2,3,\ldots ,n.} 伝達関数(またはシステム関数)は、これらの極を用いて次のように表される。
H ( s ) = G 0 ∏ k = 1 n ω c s − s k = G 0 ∏ k = 1 n ω c s − ω c e j ( 2 k + n − 1 ) π 2 n {\displaystyle H(s)=G_{0}\prod _{k=1}^{n}{\frac {\omega _{c}}{s-s_{k}}}=G_{0}\prod _{k=1}^{n}{\frac {\omega _{c}}{s-\omega _{c}e^{\frac {j(2k+n-1)\pi }{2n}}}}} 。 ここで は 数列 演算子の積です。分母は のバターワース多項式です 。 ∏ {\displaystyle \textstyle {\prod }} s {\displaystyle s}
正規化バターワース多項式 バターワース多項式は上記のように複素形式で表記されることもあるが、通常はや の ような複素共役な極対を乗じて実係数で表記される 。多項式は と設定することで正規化される 。正規化されたバターワース多項式は、一般的な積の形をとる。 s 1 {\displaystyle s_{1}} s n {\displaystyle s_{n}} ω c = 1 {\displaystyle \omega _{c}=1}
B n ( s ) = ∏ k = 1 n 2 [ s 2 − 2 s cos ( 2 k + n − 1 2 n π ) + 1 ] n = even {\displaystyle B_{n}(s)=\prod _{k=1}^{\frac {n}{2}}\left[s^{2}-2s\cos \left({\frac {2k+n-1}{2n}}\,\pi \right)+1\right]\qquad n={\text{even}}} B n ( s ) = ( s + 1 ) ∏ k = 1 n − 1 2 [ s 2 − 2 s cos ( 2 k + n − 1 2 n π ) + 1 ] n = odd . {\displaystyle 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]\qquad n={\text{odd}}.} 次の表に、1 次から 10 次までのバターワース多項式の因数を示します (小数点以下 6 桁まで)。
n バターワース多項式の因数 B n ( s ) {\displaystyle B_{n}(s)} 1 ( s + 1 ) {\displaystyle (s+1)} 2 ( s 2 + 1.414214 s + 1 ) {\displaystyle (s^{2}+1.414214s+1)} 3 ( s + 1 ) ( s 2 + s + 1 ) {\displaystyle (s+1)(s^{2}+s+1)} 4 ( s 2 + 0.765367 s + 1 ) ( s 2 + 1.847759 s + 1 ) {\displaystyle (s^{2}+0.765367s+1)(s^{2}+1.847759s+1)} 5 ( s + 1 ) ( s 2 + 0.618034 s + 1 ) ( s 2 + 1.618034 s + 1 ) {\displaystyle (s+1)(s^{2}+0.618034s+1)(s^{2}+1.618034s+1)} 6 ( s 2 + 0.517638 s + 1 ) ( s 2 + 1.414214 s + 1 ) ( s 2 + 1.931852 s + 1 ) {\displaystyle (s^{2}+0.517638s+1)(s^{2}+1.414214s+1)(s^{2}+1.931852s+1)} 7 ( s + 1 ) ( s 2 + 0.445042 s + 1 ) ( s 2 + 1.246980 s + 1 ) ( s 2 + 1.801938 s + 1 ) {\displaystyle (s+1)(s^{2}+0.445042s+1)(s^{2}+1.246980s+1)(s^{2}+1.801938s+1)} 8 ( s 2 + 0.390181 s + 1 ) ( s 2 + 1.111140 s + 1 ) ( s 2 + 1.662939 s + 1 ) ( s 2 + 1.961571 s + 1 ) {\displaystyle (s^{2}+0.390181s+1)(s^{2}+1.111140s+1)(s^{2}+1.662939s+1)(s^{2}+1.961571s+1)} 9 ( s + 1 ) ( s 2 + 0.347296 s + 1 ) ( s 2 + s + 1 ) ( s 2 + 1.532089 s + 1 ) ( s 2 + 1.879385 s + 1 ) {\displaystyle (s+1)(s^{2}+0.347296s+1)(s^{2}+s+1)(s^{2}+1.532089s+1)(s^{2}+1.879385s+1)} 10 ( s 2 + 0.312869 s + 1 ) ( s 2 + 0.907981 s + 1 ) ( s 2 + 1.414214 s + 1 ) ( s 2 + 1.782013 s + 1 ) ( s 2 + 1.975377 s + 1 ) {\displaystyle (s^{2}+0.312869s+1)(s^{2}+0.907981s+1)(s^{2}+1.414214s+1)(s^{2}+1.782013s+1)(s^{2}+1.975377s+1)}
1 次から 6 次までのバターワース多項式の因数を次の表に示します (正確)。
n バターワース多項式の因数 B n ( s ) {\displaystyle B_{n}(s)} 1 ( s + 1 ) {\displaystyle (s+1)} 2 ( s 2 + 2 s + 1 ) {\displaystyle (s^{2}+{\sqrt {2}}s+1)} 3 ( s + 1 ) ( s 2 + s + 1 ) {\displaystyle (s+1)(s^{2}+s+1)} 4 ( s 2 + 2 − 2 s + 1 ) ( s 2 + 2 + 2 s + 1 ) {\displaystyle (s^{2}+{\sqrt {2-{\sqrt {2}}}}s+1)(s^{2}+{\sqrt {2+{\sqrt {2}}}}s+1)} 5 ( s + 1 ) ( s 2 + φ − 1 s + 1 ) ( s 2 + φ s + 1 ) {\displaystyle (s+1)(s^{2}+\varphi ^{-1}s+1)(s^{2}+\varphi s+1)} 6 ( s 2 + 2 − 3 s + 1 ) ( s 2 + 2 s + 1 ) ( s 2 + 2 + 3 s + 1 ) {\displaystyle (s^{2}+{\sqrt {2-{\sqrt {3}}}}s+1)(s^{2}+{\sqrt {2}}s+1)(s^{2}+{\sqrt {2+{\sqrt {3}}}}s+1)}
ここで、ギリシャ文字の ファイ ( または )は 黄金比 を表す。これは 無理数であり、 二次方程式 の解は [4] [5] である。 φ {\displaystyle \varphi } ϕ {\displaystyle \phi } x 2 − x − 1 = 0 , {\displaystyle x^{2}-x-1=0,}
φ = 1 + 5 2 = 1.618033988749... {\displaystyle \varphi ={\frac {1+{\sqrt {5}}}{2}}=1.618033988749...} ( OEIS :A001622 ) 番目 のバターワース多項式は、 n {\displaystyle n}
B n ( s ) = ∑ k = 0 n a k s k , {\displaystyle B_{n}(s)=\sum _{k=0}^{n}a_{k}s^{k}\,,} その係数は 再帰式 [6] [7]で与えられる。 a k {\displaystyle a_{k}}
a k + 1 a k = cos ( k γ ) sin ( ( k + 1 ) γ ) {\displaystyle {\frac {a_{k+1}}{a_{k}}}={\frac {\cos(k\gamma )}{\sin((k+1)\gamma )}}} そして積の式によって
a k = ∏ μ = 1 k cos ( ( μ − 1 ) γ ) sin ( μ γ ) , {\displaystyle a_{k}=\prod _{\mu =1}^{k}{\frac {\cos((\mu -1)\gamma )}{\sin(\mu \gamma )}}\,,} どこ
a 0 = 1 and γ = π 2 n . {\displaystyle a_{0}=1\qquad {\text{and}}\qquad \gamma ={\frac {\pi }{2n}}\,.} さらに、 最初の10個のバターワース多項式の 丸められた係数 は次のとおりです。 a k = a n − k {\displaystyle a_{k}=a_{n-k}} a k {\displaystyle a_{k}} B n ( s ) {\displaystyle B_{n}(s)}
バターワース係数( 小数点4桁まで) a k {\displaystyle a_{k}} n {\displaystyle n} a 0 {\displaystyle a_{0}} a 1 {\displaystyle a_{1}} a 2 {\displaystyle a_{2}} a 3 {\displaystyle a_{3}} a 4 {\displaystyle a_{4}} a 5 {\displaystyle a_{5}} a 6 {\displaystyle a_{6}} a 7 {\displaystyle a_{7}} a 8 {\displaystyle a_{8}} a 9 {\displaystyle a_{9}} a 10 {\displaystyle a_{10}} 1 {\displaystyle 1} 1 {\displaystyle 1} 1 {\displaystyle 1} 2 {\displaystyle 2} 1 {\displaystyle 1} 1.4142 {\displaystyle 1.4142} 1 {\displaystyle 1} 3 {\displaystyle 3} 1 {\displaystyle 1} 2 {\displaystyle 2} 2 {\displaystyle 2} 1 {\displaystyle 1} 4 {\displaystyle 4} 1 {\displaystyle 1} 2.6131 {\displaystyle 2.6131} 3.4142 {\displaystyle 3.4142} 2.6131 {\displaystyle 2.6131} 1 {\displaystyle 1} 5 {\displaystyle 5} 1 {\displaystyle 1} 3.2361 {\displaystyle 3.2361} 5.2361 {\displaystyle 5.2361} 5.2361 {\displaystyle 5.2361} 3.2361 {\displaystyle 3.2361} 1 {\displaystyle 1} 6 {\displaystyle 6} 1 {\displaystyle 1} 3.8637 {\displaystyle 3.8637} 7.4641 {\displaystyle 7.4641} 9.1416 {\displaystyle 9.1416} 7.4641 {\displaystyle 7.4641} 3.8637 {\displaystyle 3.8637} 1 {\displaystyle 1} 7 {\displaystyle 7} 1 {\displaystyle 1} 4.4940 {\displaystyle 4.4940} 10.0978 {\displaystyle 10.0978} 14.5918 {\displaystyle 14.5918} 14.5918 {\displaystyle 14.5918} 10.0978 {\displaystyle 10.0978} 4.4940 {\displaystyle 4.4940} 1 {\displaystyle 1} 8 {\displaystyle 8} 1 {\displaystyle 1} 5.1258 {\displaystyle 5.1258} 13.1371 {\displaystyle 13.1371} 21.8462 {\displaystyle 21.8462} 25.6884 {\displaystyle 25.6884} 21.8462 {\displaystyle 21.8462} 13.1371 {\displaystyle 13.1371} 5.1258 {\displaystyle 5.1258} 1 {\displaystyle 1} 9 {\displaystyle 9} 1 {\displaystyle 1} 5.7588 {\displaystyle 5.7588} 16.5817 {\displaystyle 16.5817} 31.1634 {\displaystyle 31.1634} 41.9864 {\displaystyle 41.9864} 41.9864 {\displaystyle 41.9864} 31.1634 {\displaystyle 31.1634} 16.5817 {\displaystyle 16.5817} 5.7588 {\displaystyle 5.7588} 1 {\displaystyle 1} 10 {\displaystyle 10} 1 {\displaystyle 1} 6.3925 {\displaystyle 6.3925} 20.4317 {\displaystyle 20.4317} 42.8021 {\displaystyle 42.8021} 64.8824 {\displaystyle 64.8824} 74.2334 {\displaystyle 74.2334} 64.8824 {\displaystyle 64.8824} 42.8021 {\displaystyle 42.8021} 20.4317 {\displaystyle 20.4317} 6.3925 {\displaystyle 6.3925} 1 {\displaystyle 1}
正規化されたバターワース多項式は、以下のように 任意のローパスフィルタのカットオフ周波数に対する伝達関数を決定するために使用できる。 ω c {\displaystyle \omega _{c}}
H ( s ) = G 0 B n ( a ) {\displaystyle H(s)={\frac {G_{0}}{B_{n}(a)}}} 、 どこ a = s ω c . {\displaystyle a={\frac {s}{\omega _{c}}}.} 他のバンドフォームへの変換も可能です。 プロトタイプ フィルター を参照してください。
最大限の平坦性 とを仮定する と 、ゲインの周波数に関する微分は次のように表される。 ω c = 1 {\displaystyle \omega _{c}=1} G 0 = 1 {\displaystyle G_{0}=1}
d G d ω = − n G 3 ω 2 n − 1 {\displaystyle {\frac {dG}{d\omega }}=-nG^{3}\omega ^{2n-1}} ゲインは常に正なので、 これは すべて 単調 減少する。したがって、バターワースフィルタのゲイン関数にはリップルがない。ゲインの級数展開は次のように表される。 ω {\displaystyle \omega } G {\displaystyle G}
G ( ω ) = 1 − 1 2 ω 2 n + 3 8 ω 4 n + … {\displaystyle G(\omega )=1-{\frac {1}{2}}\omega ^{2n}+{\frac {3}{8}}\omega ^{4n}+\ldots } 言い換えれば、ゲインの2 次導関数を除くすべての導関数は においてゼロとなり 、「最大平坦性」が実現されます。単調性という要件が通過帯域のみに限定され、阻止帯域におけるリップルが許容される場合、逆 チェビシェフフィルタ のような、通過帯域において「最大平坦性」を持つバターワースフィルタよりも平坦な、同じ次数のフィルタを設計することが可能です。 n {\displaystyle n} ω = 0 {\displaystyle \omega =0}
高周波ロールオフ 再び と仮定すると 、 が大きい場合のゲインの対数の傾き は ω c = 1 {\displaystyle \omega _{c}=1} ω {\displaystyle \omega }
lim ω → ∞ d log ( G ) d log ( ω ) = − n . {\displaystyle \lim _{\omega \rightarrow \infty }{\frac {d\log(G)}{d\log(\omega )}}=-n.} したがって、デシベル では 、高周波ロールオフは 20 dB/ディケード、または 6 dB/オクターブになります (電力は電圧ゲインの 2 乗に比例するため、係数 20 が使用されます。20 対数規則 を参照してください)。 n {\displaystyle n} n {\displaystyle n}
最小注文数 必要最小限の要素数でバターワースフィルタを設計するには、バターワースフィルタの最小次数は次のように計算できます。 [8]
n = ⌈ log ( 10 α s / 10 − 1 10 α p / 10 − 1 ) 2 log ( ω s / ω p ) ⌉ {\displaystyle n=\left\lceil {\frac {\log {{\bigr (}{\frac {10^{\alpha _{s}/10}-1}{10^{\alpha _{p}/10}-1}}}{\bigr )}}{2\log {(\omega _{s}/\omega _{p})}}}\right\rceil } どこ:
ω p {\displaystyle \omega _{p}} および は通過帯域周波数と、その周波数における減衰量(dB)です。 α p {\displaystyle \alpha _{p}} ω s {\displaystyle \omega _{s}} および はストップバンド周波数と、その周波数での減衰量(dB)です。 α s {\displaystyle \alpha _{s}} n {\displaystyle n} 最小の極数、つまりフィルタの次数です。 ⌈ ⋅ ⌉ {\displaystyle \lceil \cdot \rceil } 天井関数 を表します 。
非標準カットオフ減衰 バターワースフィルタのカットオフ減衰量は通常 -3.01 dB と定義されます。カットオフ周波数で異なる減衰量を使用したい場合は、各極に以下の係数を適用します。この場合、極は円周上に残りますが、半径は1ではなくなります。 [8] カットオフ減衰量式は、ページ上部に記載されているバターワース定義式を代数的に操作することで導出できます。 [9]
p A = p 1 × ( 10 α / 10 − 1 ) − 1 / 2 n For 0 ≤ α < ∞ {\displaystyle {\begin{aligned}p_{A}=p_{1}\times (10^{\alpha /10}-1)^{{-1}/{2n}}&\qquad {\text{For 0}}\leq \alpha <\infty \end{aligned}}} どこ:
p A {\displaystyle p_{A}} 希望するカットオフ減衰を設定するために再配置されたポールです。 p 1 {\displaystyle p_{1}} は単位円上にある -3.01 dB カットオフ ポールです。 α {\displaystyle \alpha } カットオフ周波数における希望の減衰量(dB 単位、1 dB、10 dB など)です。 n {\displaystyle n} 極の数、つまりフィルターの次数です。
フィルタの実装と設計 線形アナログフィルタを実装するためには、いくつかの異なる フィルタトポロジーが 利用可能です。パッシブ実現で最もよく使われるトポロジーはCauerトポロジーであり、アクティブ実現で最もよく使われるトポロジーはSallen-Keyトポロジーです。
カウアー位相幾何学 カウアー位相 を用いたバターワースフィルタ カウ アートポロジーは 、受動部品(シャントコンデンサと直列インダクタ)を用いて線形アナログフィルタを実現する。与えられた伝達関数を持つバターワースフィルタは、カウアー1型を用いて実現できる。k 番目の要素は [10] で与えられる。
C k = 2 sin [ ( 2 k − 1 ) 2 n π ] k = odd {\displaystyle C_{k}=2\sin \left[{\frac {(2k-1)}{2n}}\pi \right]\qquad k={\text{odd}}} L k = 2 sin [ ( 2 k − 1 ) 2 n π ] k = even . {\displaystyle L_{k}=2\sin \left[{\frac {(2k-1)}{2n}}\pi \right]\qquad k={\text{even}}.} フィルタは必要に応じて直列インダクタから始めることができます。その場合、 L k は k の 奇数、 C k は kの偶数となります。これらの式は、 L k と C k の 両方をg k と することで組み合わせることができます 。つまり、 g k はイミタンスを s で割った 値です 。
g k = 2 sin [ ( 2 k − 1 ) 2 n π ] k = 1 , 2 , 3 , … , n . {\displaystyle g_{k}=2\sin \left[{\frac {(2k-1)}{2n}}\pi \right]\qquad k=1,2,3,\ldots ,n.} これらの式は、ω c = 1 の両端終端フィルタ(つまり、信号源インピーダンスと負荷インピーダンスがともに1)に適用されます。この プロトタイプフィルタは、他のインピーダンス値と周波数に対してもスケーリング可能です。片端終端フィルタ(つまり、理想的な電圧源または電流源で駆動されるフィルタ)の場合、素子値は [3] で与えられます。
g j = a j a j − 1 c j − 1 g j − 1 j = 2 , 3 , … , n {\displaystyle g_{j}={\frac {a_{j}a_{j-1}}{c_{j-1}g_{j-1}}}\qquad j=2,3,\ldots ,n} どこ
g 1 = a 1 {\displaystyle g_{1}=a_{1}} そして
a j = sin [ ( 2 j − 1 ) 2 n π ] j = 1 , 2 , 3 , … , n {\displaystyle a_{j}=\sin \left[{\frac {(2j-1)}{2n}}\pi \right]\qquad j=1,2,3,\ldots ,n} c j = cos 2 [ j 2 n π ] j = 1 , 2 , 3 , … , n . {\displaystyle c_{j}=\cos ^{2}\left[{\frac {j}{2n}}\pi \right]\qquad j=1,2,3,\ldots ,n.} 電圧駆動フィルタは直列素子から、電流駆動フィルタはシャント素子から始める必要があります。これらの形式は、ダイプレクサ や マルチプレクサ の設計に役立ちます 。 [3]
サレン・キー位相 サレン・キー位相 Sallen -Keyトポロジーは 、能動部品と受動部品(非反転バッファ、通常は オペアンプ 、抵抗、コンデンサ)を用いて線形アナログフィルタを実装します。各Sallen-Key段は共役極対を実装し、フィルタ全体はすべての段を直列にカスケード接続することで実装されます。実極が存在する場合(が奇数の場合)、通常は RC回路 として別途実装し 、能動段とカスケード接続する必要があります。 n {\displaystyle n}
右に示す2次サレン・キー回路の伝達関数は次のように表される。
H ( s ) = V out ( s ) V in ( s ) = 1 1 + C 2 ( R 1 + R 2 ) s + C 1 C 2 R 1 R 2 s 2 . {\displaystyle H(s)={\frac {V_{\text{out}}(s)}{V_{\text{in}}(s)}}={\frac {1}{1+C_{2}(R_{1}+R_{2})s+C_{1}C_{2}R_{1}R_{2}s^{2}}}.} 分母がバターワース多項式の2次項のいずれかであることを望みます。 と仮定すると 、これは次のことを意味します。 ω c = 1 {\displaystyle \omega _{c}=1}
C 1 C 2 R 1 R 2 = 1 {\displaystyle C_{1}C_{2}R_{1}R_{2}=1\,} そして
C 2 ( R 1 + R 2 ) = − 2 cos ( 2 k + n − 1 2 n π ) . {\displaystyle C_{2}(R_{1}+R_{2})=-2\cos \left({\frac {2k+n-1}{2n}}\pi \right).} これにより、自由に選択できる 2 つの未定義のコンポーネント値が残ります。
1つのオペアンプ のみを使用した3次および4次のサレン・キートポロジーのバターワースローパスフィルタは、 ヒュールスマンによって説明されています。 [11] [12] また、さらに高次の単一アンプバターワースフィルタは、ジュリシッチらによって提供されています。 [13]
デジタル実装 バターワースフィルタなどのデジタルフィルタの実装は、多くの場合、 双一次変換 法または 整合Z変換法 に基づいています。これらは、アナログフィルタ設計を離散化する2つの異なる手法です。バターワースフィルタなどの全極フィルタの場合、整合Z変換法は インパルス不変法と同等です。高次になると、デジタルフィルタは量子化誤差の影響を受けやすくなるため、多くの場合、カスケード接続された 双二次セクション と、奇数次の場合は1次または3次セクションを1つ追加して 計算されます。
他の線形フィルタとの比較 バターワース フィルタのプロパティは次のとおりです。
これは、離散時間バターワースフィルタのゲインを他の一般的なフィルタタイプと並べて示した画像です。これらのフィルタはすべて5次です。
バターワース フィルターは、 チェビシェフ フィルター や 楕円フィルター よりも、カットオフ周波数付近でゆっくりとロールオフしますが、リップルはありません。
参照
参考文献 ^ ab Butterworth, S. (1930). 「フィルタ増幅器の理論について」 (PDF) . 『実験的無線と無線技術者』 . 7 : 536– 541. ^ Bianchi, Giovanni; Sorrentino, Roberto (2007). 電子フィルタのシミュレーションと設計. McGraw-Hill Professional. pp. 17– 20. ISBN 978-0-07-149467-0 。 ^ abc Matthaei, George L.; Young, Leo; Jones, EMT (1964). マイクロ波フィルタ、インピーダンス整合回路、および結合構造 . McGraw-Hill. pp. 104– 107, 105, and 974. LCCN 64007937. ^ Weisstein, Eric W. 「黄金比」. mathworld.wolfram.com . 2020年8月10日 閲覧 。 ^ OEIS : A001622 ^ ボッセ、G. (1951)。 「Siebketten ohne Dämpfungsschwankungen im Durchlaßbereich (Potenzketten)」。 周波数 。 5 (10): 279–284 。 Bibcode :1951Freq....5..279B。 土井 :10.1515/FREQ.1951.5.10.279。 S2CID 124123311。 ^ Weinberg, Louis (1962). ネットワーク分析と統合. マラバー, フロリダ: Robert E. Krieger Publishing Company, Inc. (1975年出版). pp. 494– 496. hdl :2027/mdp.39015000986086. ISBN 0-88275-321-5 . 2022年6月18日 閲覧 。 ^ ab Paarmann, Larry D. (2001). アナログフィルタの設計と解析:信号処理の観点. 米国マサチューセッツ州ノーウェル:Kluwer Academic Publishers. pp. 117, 118. ISBN 0-7923-7373-1 。 ^ バターワースフィルタ#原著論文 ^ US 1849656、ウィリアム・R・ベネット、「伝送ネットワーク」、1932年3月15日公開 ^ Huelsman, LP (1971年5月). 「等価容量アクティブ RC ネットワークによる3次ローパス・バターワース特性の実現」. Electronics Letters . 7 (10): 271– 272. Bibcode :1971ElL.....7..271H. doi :10.1049/el:19710185. ^ Huelsman, LP (1974年12月). 「等価容量アクティブ RC ネットワークによる4次ローパス・バターワース特性の実現」. Proceedings of the IEEE . 62 (12): 1709. doi :10.1109/PROC.1974.9689. ^ ジュリシッチ、ドラジェン;モシッツ、ジョージ S.ミジャット、ネブン (2008)。 「テーブルを使用した低感度、単一アンプ、アクティブ RC 全極フィルター」。 オートマティカ 。 49 ( 3–4 ): 159–173 .