チェビシェフフィルタの周波数特性
本稿では、チェビシェフフィルタの振幅特性、位相特性、群遅延特性を詳細に解説する。 バターワースフィルタとの比較を通じて、チェビシェフフィルタの利点と欠点を明らかにする。 また、RF用途で重要なVSWR(電圧定在波比)についても説明する。
振幅特性
通過域の等リップル特性
チェビシェフI型フィルタの最大の特徴は、通過域での等リップル特性である。
チェビシェフ多項式 $T_n(\omega/\omega_c)$ は区間 $[-1,1]$ 上で $+1$ と $-1$ の間を $n$ 回交替(等リップル)する。 これに対応して、振幅 $|H|$ は通過域で 0 dB($T_n=0$ の点)と $-r$ dB($|T_n|=1$ の点)の間を等リップルで振動する。
リップルと特性の関係: チェビシェフフィルタのリップルを小さくするほど、通過域は平坦に近づき、遷移は緩やかになる。 すなわち実用上はバターワース的な性質(平坦な通過域・緩やかなロールオフ)に近づく。 ただし、バターワース(最大平坦)とチェビシェフ(等リップル)は異なる多項式族($\omega^n$ と $T_n$)に基づくため、 $\varepsilon \to 0$ の極限で厳密にバターワースへ収束するわけではない点には注意したい。
遷移帯域の急峻さ(ロールオフ特性)
チェビシェフフィルタは、同じ次数のバターワースフィルタより急峻な遷移特性を持つ。 フィルタの減衰がどれだけ急速に増加するかをロールオフ(roll-off)と呼ぶ。
dB/octave と dB/decade
ロールオフの急峻さは、周波数の増加に対する減衰量の増加率で表す。 よく使われる単位が dB/octave と dB/decade である。
- オクターブ(octave)
- 周波数が2倍になること。例: 1 kHz → 2 kHz → 4 kHz はそれぞれ1オクターブ上昇
- ディケード(decade)
- 周波数が10倍になること。例: 1 kHz → 10 kHz → 100 kHz はそれぞれ1ディケード上昇
これらの単位の変換は以下の関係式による:
\begin{equation} 1\text{ decade} = \log_2(10) \text{ octaves} \approx 3.32 \text{ octaves} \end{equation}したがって、20 dB/decade ≈ 6 dB/octave となる(正確には 6.02 dB/octave)。
n次フィルタの漸近的ロールオフ
すべての $n$ 次ローパスフィルタ(バターワース、チェビシェフ、楕円など)は、 十分に高い周波数で以下の漸近的ロールオフを示す:
\begin{equation} \text{ロールオフ} = 20n \text{ dB/decade} = 6n \text{ dB/octave} \end{equation}| 次数 $n$ | ロールオフ (dB/decade) | ロールオフ (dB/octave) |
|---|---|---|
| 1 | 20 | 6 |
| 2 | 40 | 12 |
| 3 | 60 | 18 |
| 4 | 80 | 24 |
| 5 | 100 | 30 |
| 6 | 120 | 36 |
これは「漸近的」な値であり、カットオフ周波数から十分離れた領域($\omega \gg \omega_c$)で成り立つ。
チェビシェフの優位性:遷移帯域での減衰
漸近的ロールオフはすべてのフィルタで同じだが、遷移帯域(カットオフから阻止域までの領域)での振る舞いが異なる。 チェビシェフフィルタは、通過域でのリップルを許容する代わりに、遷移帯域での減衰がより急峻になる。
阻止域での減衰量 $A(\omega) = -20\log_{10}|H(j\omega)|$(正の dB。表の利得 dB とは符号が逆)は、$\omega \gg \omega_c$ で次のように近似できる:
\begin{equation} A(\omega) \approx 20\log_{10}\varepsilon + 6.02\,(n-1) + 20n\log_{10}\left(\dfrac{\omega}{\omega_c}\right) \quad (\omega \gg \omega_c) \end{equation}ここで $\varepsilon$ はリップルパラメータ、$6.02\,(n-1)$ dB は $T_n$ の最高次係数 $2^{\,n-1}$ に由来する項である。 バターワース(同じ $n$ では $A \approx 20n\log_{10}(\omega/\omega_c)$)と比べると、チェビシェフは $20\log_{10}\varepsilon + 6.02\,(n-1)$ dB だけ減衰量が大きい(例: 6 次・1 dB では約 24 dB)。 これが遷移帯域を超えた後の急峻さの差として現れる。
バターワースとの比較
| 周波数 | バターワース 6次 | チェビシェフ 6次 (1dB) |
|---|---|---|
| $\omega = \omega_c$ | -3.0 dB | -1.0 dB |
| $\omega = 1.5\omega_c$ | -21.2 dB | -38.3 dB |
| $\omega = 2\omega_c$ | -36.1 dB | -56.7 dB |
| $\omega = 3\omega_c$ | -57.3 dB | -80.0 dB |
同じ次数で比較すると、チェビシェフフィルタは遷移帯域を超えた後の減衰が約 20 dB 大きい。
位相特性
位相応答
フィルタの位相応答 $\phi(\omega)$ は、伝達関数 $H(j\omega)$ の偏角として定義される:
\begin{equation} \phi(\omega) = \arg H(j\omega) = \tan^{-1}\left(\dfrac{\text{Im}[H(j\omega)]}{\text{Re}[H(j\omega)]}\right) \label{eq:phase} \end{equation}$n$ 次の全極型ローパスフィルタでは、位相は $\omega = 0$ で $0°$ から始まり、 $\omega \to \infty$ で $-n \times 90°$ に漸近する。
チェビシェフの位相特性の問題点
チェビシェフフィルタの位相特性は、バターワースフィルタより非線形性が強い。 特に通過域端(カットオフ周波数付近)で位相の変化が急激になる。
これは、チェビシェフフィルタの極が楕円上に分布し、実軸により近い位置にあるためである。 極が実軸に近いほど、その極に対応する周波数での位相変化が急になる。
位相の線形性が重要な用途: オーディオ信号処理やビデオ信号処理など、波形の形状を保持したい用途では 位相の線形性が重要である。このような用途ではベッセルフィルタや 線形位相FIRフィルタが好まれる。
群遅延特性
群遅延の定義
群遅延(group delay)$\tau_g(\omega)$ は、位相の周波数に対する微分の負値として定義される:
\begin{equation} \tau_g(\omega) = -\dfrac{d\phi(\omega)}{d\omega} \label{eq:group-delay} \end{equation}群遅延は、その周波数成分の信号が通過するのに要する時間を表す。 狭帯域信号の包絡線遅延は、中心周波数での群遅延で近似できる。
チェビシェフの群遅延特性
チェビシェフフィルタの群遅延は、以下の特徴を持つ:
- 通過域でのリップル:振幅特性と同様に、群遅延も通過域で振動する
- カットオフ付近でのピーク:カットオフ周波数付近で群遅延が急激に増加する
- リップルが大きいほど群遅延のピークも大きい
2次セクションの群遅延
極 $p = -\sigma + j\omega_0$($\sigma > 0$)を持つ2次セクションの群遅延は:
\begin{equation} \tau_g(\omega) = \dfrac{2\sigma}{\sigma^2 + (\omega - \omega_0)^2} + \dfrac{2\sigma}{\sigma^2 + (\omega + \omega_0)^2} \label{eq:tau-2nd} \end{equation}極が実軸に近い($\sigma$ が小さい)ほど、$\omega = \omega_0$ 付近での群遅延ピークが大きくなる。 チェビシェフフィルタの極はバターワースより実軸に近いため、群遅延のピークが大きい。
群遅延の数値例
| フィルタ | DC での群遅延 | 最大群遅延 | 変動比 |
|---|---|---|---|
| バターワース 6次 | $6.0 / \omega_c$ | $8.2 / \omega_c$ | 1.37 |
| チェビシェフ 6次 (0.5dB) | $7.5 / \omega_c$ | $16.8 / \omega_c$ | 2.24 |
| チェビシェフ 6次 (1dB) | $8.7 / \omega_c$ | $22.1 / \omega_c$ | 2.54 |
| チェビシェフ 6次 (3dB) | $12.4 / \omega_c$ | $45.3 / \omega_c$ | 3.65 |
リップルが大きいほど群遅延の変動が激しくなる。 3 dB リップルでは、群遅延の最大値が DC の約 3.6 倍に達する。
群遅延の歪みによる波形劣化: 群遅延が周波数によって大きく変化すると、パルス信号・ステップ応答・方形波など高調波成分を多く含む信号で、 通過帯域内の群遅延の不均一性によりオーバーシュートやリンギングが生じ、波形が著しく歪む。 時間領域での波形品質が重要な用途では、リップルを小さく設定して群遅延を平坦化するか、 群遅延が平坦なベッセルフィルタや線形位相FIRフィルタ(ガウシアンフィルタなど)を検討する。
Python による群遅延の計算
import numpy as np
from scipy.signal import cheby1, group_delay
import matplotlib.pyplot as plt
# 6次チェビシェフI型、1dBリップル
n = 6
rp = 1.0
b, a = cheby1(n, rp, 1, analog=False, fs=100)
# 群遅延の計算
w, gd = group_delay((b, a), fs=100)
# プロット
plt.figure(figsize=(10, 4))
plt.plot(w, gd)
plt.xlabel('Frequency [Hz]')
plt.ylabel('Group delay [samples]')
plt.title(f'Chebyshev Type I (n={n}, ripple={rp}dB) Group Delay')
plt.grid(True)
plt.xlim(0, 5)
plt.show()
ポールQと実装上の制約
チェビシェフフィルタの極(ポール)は楕円上に配置され、そのQ値(Quality Factor)は 次数とリップルによって決まる。Q値が高い極は、回路実装において以下の問題を引き起こす。
Q値とは
2次セクションの極 $p = \sigma + j\omega_d$($\sigma = \mathrm{Re}(p) < 0$)に対して、 Q値は固有角周波数 $\omega_0 = |p|$(極の絶対値)を用いて以下で定義される:
\begin{equation} Q = \dfrac{|p|}{2\,|\mathrm{Re}(p)|} = \dfrac{\omega_0}{2|\sigma|}, \qquad \omega_0 = |p| = \sqrt{\sigma^2 + \omega_d^2} \end{equation}Q値が高いほど、極が虚軸に近い(減衰が小さい)ことを意味し、 周波数応答にシャープなピークが現れる。
チェビシェフフィルタのQ値の傾向
チェビシェフフィルタでは、次数が高いほど、またリップルが大きいほど、 最大Q値が急激に増加する。以下は最も高いQ値を持つ2次セクションのQ値である。
| 次数 $n$ | 0.1 dB リップル | 0.5 dB リップル | 1.0 dB リップル | 3.0 dB リップル |
|---|---|---|---|---|
| 2 | 0.96 | 0.86 | 0.96 | 1.30 |
| 4 | 1.80 | 2.94 | 3.56 | 5.58 |
| 6 | 2.79 | 5.56 | 7.99 | 17.4 |
| 8 | 3.84 | 8.41 | 13.0 | 36.4 |
| 10 | 4.91 | 11.4 | 18.4 | 62.9 |
高Q値がもたらす問題
- 部品感度:Q値が高いと、抵抗やコンデンサの値が少し変わっただけで特性が大きく変化する
- ノイズ増幅:高Qセクションは特定の周波数でゲインが高くなり、ノイズを増幅しやすい
- 発振のリスク:オペアンプの帯域やスルーレートが不足すると、発振することがある
- 部品の現実的な制約:Sallen-Key回路では $Q > 10$ で部品比が極端になる
バターワースフィルタとの部品感度比較
チェビシェフフィルタは、同じ次数のバターワースフィルタと比較して部品公差に対する感度が高い。 これは、チェビシェフフィルタの極がより虚軸に近い位置にあり、高いQ値を持つためである。
| フィルタタイプ | 部品公差感度 | 推奨部品精度 |
|---|---|---|
| バターワース | 低〜中 | 5% でも実用可 |
| チェビシェフ(0.5 dB以下) | 中 | 2〜5% |
| チェビシェフ(1 dB以上) | 高 | 1〜2% |
高精度な部品が入手困難な場合や、温度変化が大きい環境では、 チェビシェフよりバターワースを選択する方が安定した特性が得られることがある。
実用上の目安: アクティブフィルタ(Sallen-Key、MFB)では、$Q < 10$ 程度に抑えることが望ましい。 高Qが必要な場合は、状態変数フィルタやバイカッド回路を使うか、 低リップルに設計して次数を増やす方が安定した実装になることが多い。
Q値を下げるための選択肢
- リップルを小さくする:0.1 dB 以下のリップルではQ値が大幅に下がる
- 次数を下げる:必要な減衰特性をぎりぎり満たす最小次数を選ぶ
- バターワースを検討:同じ次数ではQ値が低く、実装が容易
インパルス応答とステップ応答
時間領域での特性
フィルタの時間領域特性は、周波数特性(特に群遅延)と密接に関連する。
- インパルス応答:群遅延の変動が大きいほど、振動(リンギング)が長く続く
- ステップ応答:群遅延のピークが大きいほど、オーバーシュートが大きくなる
オーバーシュートの比較
| フィルタ | ステップ応答のオーバーシュート |
|---|---|
| ベッセル 6次 | 0.8% |
| バターワース 6次 | 14.3% |
| チェビシェフ 6次 (0.5dB) | 18.7% |
| チェビシェフ 6次 (1dB) | 21.0% |
| チェビシェフ 6次 (3dB) | 27.5% |
チェビシェフフィルタは周波数選択性は優れているが、時間領域での過渡応答は劣る。 用途に応じてトレードオフを考慮する必要がある。
VSWR(電圧定在波比)
RF フィルタにおけるVSWR
RF(高周波)やマイクロ波用途では、フィルタのインピーダンス整合が重要である。 VSWR(Voltage Standing Wave Ratio、電圧定在波比)は、 インピーダンス不整合の程度を表す指標である。
\begin{equation} \text{VSWR} = \dfrac{1 + |\Gamma|}{1 - |\Gamma|} \label{eq:vswr} \end{equation}ここで $\Gamma$ は反射係数である。VSWR = 1 が理想(完全整合)、値が大きいほど不整合が大きい。
チェビシェフフィルタとVSWR
両端を定抵抗で終端した無損失 LC ラダーとして実現したチェビシェフフィルタは、通過域で等リップルの反射特性を示す。 これはVSWRの観点からは有利な特性である(この性質はこのプロトタイプを前提とするものであり、アクティブフィルタや任意の実装で一般に成り立つわけではない)。
通過域リップル $r$ [dB] と最大 VSWR の関係:
\begin{equation} \text{VSWR}_{\max} = \dfrac{1 + \sqrt{1 - 10^{-r/10}}}{1 - \sqrt{1 - 10^{-r/10}}} \label{eq:vswr-ripple} \end{equation}| リップル [dB] | 最大 VSWR | 反射損失 [dB] |
|---|---|---|
| 0.010 | 1.10 | 26.4 |
| 0.036 | 1.20 | 20.8 |
| 0.075 | 1.30 | 17.7 |
| 0.122 | 1.40 | 15.6 |
| 0.177 | 1.50 | 14.0 |
| 0.302 | 1.70 | 11.7 |
| 0.512 | 2.00 | 9.5 |
RF設計での推奨: 多くのRFシステムでは VSWR ≤ 1.5 が要求される。 これに対応するリップルは約 0.18 dB である。 より厳しい VSWR ≤ 1.2 が必要な場合は、リップル約 0.04 dB を選択する。
バターワースとの VSWR 比較
バターワースフィルタは通過域で最大平坦特性を持つが、VSWR は通過域端で最大となり、 DC に向かって徐々に改善する。
チェビシェフフィルタは VSWR が通過域全体で等リップルとなるため、 最悪ケースの VSWR を保証しやすいという利点がある。 これが RF フィルタでチェビシェフ特性が好まれる理由の一つである。
VSWRと接続時の安定性
チェビシェフフィルタの等リップル特性は、他の回路との接続において有利に働く。 通過域でのリップルは、言い換えれば通過帯域全体で VSWR が一定の範囲内に収まることを意味する。
これにより、バターワースフィルタのような最大平坦特性を持つフィルタと比較して、 接続された回路との間での信号の反射によるインピーダンス変動(いわゆる「信号のあばれ」)が 予測可能で均一になり、システム全体の設計が容易になる。
RF/マイクロ波設計での利点: 通過帯域内でVSWRが保証されることで、増幅器やミキサーなど後段回路の 入力インピーダンス整合設計が容易になる。特にカスケード接続された 複数のフィルタ間の相互作用を予測しやすい。
II型の特性
振幅特性
チェビシェフII型は、通過域が単調減少で阻止域に等リップルを持つ。 阻止域では減衰量が一定の下限値(設計時に指定する阻止域減衰量)を下回ることはない。
位相・群遅延特性
II型は通過域に零点がないため、I型より位相特性が良好である。 群遅延の変動も I型より小さい傾向がある。
ただし、阻止域に零点(伝達零点)があるため、阻止域での群遅延は複雑な挙動を示す。
Python による特性解析
import numpy as np
from scipy.signal import cheby1, cheby2, butter, freqs, group_delay
import matplotlib.pyplot as plt
def analyze_filter(b, a, label, omega):
"""振幅、位相、群遅延を計算"""
w, H = freqs(b, a, omega)
magnitude_db = 20 * np.log10(np.abs(H))
phase_deg = np.unwrap(np.angle(H)) * 180 / np.pi
# 群遅延(数値微分)
phase_rad = np.unwrap(np.angle(H))
group_delay = -np.gradient(phase_rad, w)
return magnitude_db, phase_deg, group_delay
# 周波数軸
omega = np.logspace(-1, 1, 1000)
# フィルタ設計
n = 6
b_butter, a_butter = butter(n, 1, analog=True)
b_cheby1_05, a_cheby1_05 = cheby1(n, 0.5, 1, analog=True)
b_cheby1_1, a_cheby1_1 = cheby1(n, 1.0, 1, analog=True)
b_cheby2, a_cheby2 = cheby2(n, 40, 1, analog=True)
# 解析
filters = [
(b_butter, a_butter, 'Butterworth'),
(b_cheby1_05, a_cheby1_05, 'Chebyshev I (0.5dB)'),
(b_cheby1_1, a_cheby1_1, 'Chebyshev I (1dB)'),
(b_cheby2, a_cheby2, 'Chebyshev II (40dB)')
]
fig, axes = plt.subplots(3, 1, figsize=(10, 10))
for b, a, label in filters:
mag, phase, gd = analyze_filter(b, a, label, omega)
axes[0].semilogx(omega, mag, label=label)
axes[1].semilogx(omega, phase, label=label)
axes[2].semilogx(omega, gd, label=label)
axes[0].set_ylabel('Magnitude [dB]')
axes[0].set_ylim(-80, 5)
axes[0].legend()
axes[0].grid(True)
axes[1].set_ylabel('Phase [degrees]')
axes[1].legend()
axes[1].grid(True)
axes[2].set_xlabel('Frequency [rad/s]')
axes[2].set_ylabel('Group Delay [s]')
axes[2].set_ylim(0, 30)
axes[2].legend()
axes[2].grid(True)
plt.tight_layout()
plt.show()
フィルタタイプ総合比較
チェビシェフフィルタを他の代表的なフィルタタイプと比較する。 用途に応じた適切なフィルタ選択の参考にされたい。
| 特性 | ベッセル | バターワース | チェビシェフ | 楕円 |
|---|---|---|---|---|
| 群遅延の平坦性 | 最良 | 良 | 悪い | 最悪 |
| 通過域の平坦性 | 悪い | 最良(最大平坦) | リップルあり | リップルあり |
| ロールオフの急峻さ | 最も緩やか | 緩やか | 急峻 | 最も急峻 |
| 位相の線形性 | 最良(線形) | 良 | 非線形 | 最も非線形 |
| ステップ応答の オーバーシュート |
最小(≈1%) | 中程度(≈14%) | 大きい(18〜28%) | 最大 |
| 部品感度 | 低 | 低〜中 | 中〜高 | 高 |
| 同じ減衰仕様での次数 | 最大 | 大 | 中 | 最小 |
用途別推奨フィルタ
- 波形を保持したい(オーディオ、映像、データ通信)
- ベッセルフィルタ、または線形位相FIRフィルタ
- 通過域の平坦性が重要(計測機器、センサー信号処理)
- バターワースフィルタ
- 隣接チャンネル除去、周波数選択性が重要(RF、無線通信)
- チェビシェフフィルタ(I型またはII型)
- 最小次数で最大の遷移帯域急峻さ(部品数制限がある場合)
- 楕円フィルタ
- レーダー信号処理、無線通信システム
- チェビシェフまたは楕円フィルタ(アプリケーション要件による)
設計のトレードオフ: フィルタ設計は常にトレードオフの問題である。 ロールオフを急峻にすれば群遅延特性が悪化し、 群遅延を平坦にすればロールオフが緩やかになる。 要求仕様の中で何を優先するかを明確にしてからフィルタタイプを選択することが重要である。
よくある質問
Q1: チェビシェフフィルタとはどのようなフィルタか
A: チェビシェフフィルタはチェビシェフ多項式を用いた伝達関数を持つアナログフィルタで、タイプI(通過域等リップル・阻止域単調)とタイプII(通過域単調・阻止域等リップル)がある。バターワースより急峻な遷移特性を同次数で実現できる。
Q2: チェビシェフフィルタのリップルと次数はどうトレードオフするか
A: 通過域リップル量 ε を大きくすると、同じ次数でより急峻な遷移特性(高い阻止域減衰量)が得られる。逆にリップルを小さくすると急峻性が下がる。設計仕様(通過域リップル・阻止域減衰量・遷移帯域幅)から必要最小次数を計算式で求めることができる。
Q3: チェビシェフ多項式はフィルタ設計においてなぜ重要か
A: チェビシェフ多項式 T_n(x) は [-1,1] 上で|T_n(x)|<=1 を満たしながら他の任意のn次多項式(最高次係数同じ)より最大値が小さい等偏差性を持つ。この性質を利用してフィルタの周波数応答を等リップル化し、同次数で最も急峻な遷移を実現できる。