楕円フィルタ

楕円フィルタ（ヴィルヘルム・カウアーにちなんでカウアーフィルタとも呼ばれる、あるいはエゴール・ゾロタレフにちなんでゾロタレフフィルタとも呼ばれる）は、通過帯域と阻止帯域の両方で等化リップル（等リップル）動作を行う信号処理フィルタである。各帯域のリップルの量は独立に調整可能であり、リップルの値が与えられた場合（リップルが等化されているかどうかに関係なく）、同じ次数のフィルタで通過帯域と阻止帯域の間でのゲイン遷移がこれより高速になるものは他にない。あるいは、通過帯域と阻止帯域のリップルを独立に調整する機能をあきらめて、部品のばらつきの影響を最大限受けないフィルタを設計することもできる。

阻止帯域のリップルがゼロに近づくと、フィルタはタイプIチェビシェフフィルタになります。通過帯域のリップルがゼロに近づくと、フィルタはタイプIIチェビシェフフィルタになり、最終的に両方のリップル値がゼロに近づくと、フィルタはバターワースフィルタになります。

角周波数ωの関数としての ローパス楕円フィルタのゲインは次のように表されます。

${\displaystyle G_{n}(\omega )={1 \over {\sqrt {1+\epsilon ^{2}R_{n}^{2}(\xi ,\omega /\omega _{0})}}}}$

ここで、R _nはn次の楕円有理関数（チェビシェフ有理関数とも呼ばれる）であり、

${\displaystyle \omega _{0}}$ カットオフ周波数

${\displaystyle \epsilon }$ 波及要因

${\displaystyle \xi}$ 選択係数である

リップル係数の値は通過帯域リップルを指定し、リップル係数と選択係数の組み合わせは阻止帯域リップルを指定します。

プロパティ

ε  = 0.5、ξ = 1.05の4次楕円ローパスフィルタの周波数応答 。通過帯域の最小ゲインと阻止帯域の最大ゲイン、および正規化周波数1とξの間の遷移領域も示されている。

上記のグラフの遷移領域のクローズアップ。

• 通過帯域において、楕円有理関数は0から1の間で変化します。したがって、通過帯域のゲインは1から の間で変化します。${\displaystyle 1/{\sqrt {1+\epsilon ^{2}}}}$
• ストップバンドでは、楕円有理関数は無限大と次のように定義される識別係数の間で変化します。${\displaystyle L_{n}}$

${\displaystyle L_{n}=R_{n}(\xi ,\xi )\,}$

したがって、ストップバンドのゲインは 0 から の間で変化します。${\displaystyle 1/{\sqrt {1+\epsilon^{2}L_{n}^{2}}}}$

• 楕円有理関数の極限ではチェビシェフ多項式となり、したがってフィルタはリップル係数εを持つチェビシェフI型フィルタとなる。${\displaystyle \xi \rightarrow \infty }$
• バターワースフィルタはチェビシェフフィルタの極限形であるため、 の極限において、フィルタがバターワースフィルタとなることが分かる。${\displaystyle \xi \rightarrow \infty }$ ${\displaystyle \omega _{0}\rightarrow 0}$ ${\displaystyle \epsilon \rightarrow 0}$ ${\displaystyle \epsilon \,R_{n}(\xi ,1/\omega _{0})=1}$
• の極限、および 、およびとなる場合、フィルタはゲイン${\displaystyle \xi \rightarrow \infty }$ ${\displaystyle \epsilon \rightarrow 0}$ ${\displaystyle \omega _{0}\rightarrow 0}$ ${\displaystyle \xi \omega _{0}=1}$ ${\displaystyle \epsilon L_{n}=\alpha }$

${\displaystyle G(\omega )={\frac {1}{\sqrt {1+{\frac {1}{\alpha ^{2}T_{n}^{2}(1/\omega )}}}}}}$

極と零点

複素周波数空間（s = σ + jω）における8次楕円フィルタのゲインの絶対値の対数。ε = 0.5、ξ = 1.05、ω ₀ = 1とする。白い点は極、黒い点は零点である。極と零点は合計16個ある。遷移領域付近で1つの極と零点のように見えるものは、実際には下の拡大図に示すように4つの極と2つの零点である。この図では、黒はゲインが0.0001以下、白はゲインが10以上である。

上の画像の遷移領域の拡大図。4 つの極と 2 つの二重ゼロが示されています。

楕円フィルタのゲインのゼロは、楕円有理関数の記事で導出される楕円有理関数の極と一致します。

楕円フィルタのゲインの極は、タイプIチェビシェフフィルタのゲインの極の導出と非常によく似た方法で導出できます。簡単のため、カットオフ周波数は1であると仮定します。楕円フィルタのゲインの極は、ゲインの分母の零点になります。複素周波数を用いると、これは次の式を意味します。 ${\displaystyle (\omega _{pm})}$ ${\displaystyle s=\sigma +j\omega }$

${\displaystyle 1+\epsilon ^{2}R_{n}^{2}(-js,\xi )=0\,}$

cd() をJacobi の楕円余弦関数として定義し、楕円有理関数の定義を使用すると次のようになります。 ${\displaystyle -js=\mathrm {cd} (w,1/\xi )}$

${\displaystyle 1+\epsilon ^{2}\mathrm {cd} ^{2}\left({\frac {nwK_{n}}{K}},{\frac {1}{L_{n}}}\right)=0\,}$

ここで、および。wを解くと${\displaystyle K=K(1/\xi )}$ ${\displaystyle K_{n}=K(1/L_{n})}$

${\displaystyle w={\frac {K}{nK_{n}}}\mathrm {cd} ^{-1}\left({\frac {\pm j}{\epsilon }},{\frac {1}{L_{n}}}\right)+{\frac {mK}{n}}}$

ここで、逆 cd() 関数の複数の値は、整数インデックスmを使用して明示的に指定されます。

楕円ゲイン関数の極は次のようになります。

${\displaystyle s_{pm}=i\,\mathrm {cd} (w,1/\xi )\,}$

チェビシェフ多項式の場合と同様に、これは明示的に複素形式で表現することができる（Lutovac & et al. 2001、§12.8）。

${\displaystyle s_{pm}={\frac {a+jb}{c}}}$

${\displaystyle a=-\zeta _{n}{\sqrt {1-\zeta _{n}^{2}}}{\sqrt {1-x_{m}^{2}}}{\sqrt {1-x_{m}^{2}/\xi ^{2}}}}$

${\displaystyle b=x_{m}{\sqrt {1-\zeta _{n}^{2}(1-1/\xi ^{2})}}}$

${\displaystyle c=1-\zeta _{n}^{2}+x_{i}^{2}\zeta _{n}^{2}/\xi ^{2}}$

ここでは の関数であり、と は楕円有理関数の零点である。はすべてのnに対してヤコビ楕円関数で表現できる。あるいは、特に 1、2、3 次では代数的に表現できる。 1 次と 2 次では ${\displaystyle \zeta _{n}}$ ${\displaystyle n,\,\epsilon }$ ${\displaystyle \xi }$ ${\displaystyle x_{m}}$ ${\displaystyle \zeta _{n}}$

${\displaystyle \zeta _{1}={\frac {1}{\sqrt {1+\epsilon ^{2}}}}}$

${\displaystyle \zeta _{2}={\frac {2}{(1+t){\sqrt {1+\epsilon ^{2}}}+{\sqrt {(1-t)^{2}+\epsilon ^{2}(1+t)^{2}}}}}}$

どこ

${\displaystyle t={\sqrt {1-1/\xi ^{2}}}}$

の代数表現はかなり複雑です（Lutovac & et al. (2001、§ 12.8.1)を参照）。 ${\displaystyle \zeta _{3}}$

楕円有理関数のネスト特性を利用して、次のような高階表現を構築することができます。 ${\displaystyle \zeta _{n}}$

${\displaystyle \zeta _{m\cdot n}(\xi ,\epsilon )=\zeta _{m}\left(\xi ,{\sqrt {{\frac {1}{\zeta _{n}^{2}(L_{m},\epsilon )}}-1}}\right)}$

どこ。 ${\displaystyle L_{m}=R_{m}(\xi ,\xi )}$

設計上の考慮事項

ξ = 1.1の8次楕円フィルタの各極の正規化されたQ係数を、リップル係数ε の関数として示したものです。複素共役極ペアと正負極ペアは同じQ係数を持つため、各曲線は4つの極を表しています。（青とシアンの曲線はほぼ一致しています。）すべての極のQ係数は、ε _Qmin  = 1 /  √ L _n  = 0.02323 で同時に最小化されます。

Lutovac & et al.を参照してください。 (2001、§ 12.11、13.14)。

楕円フィルタは一般に、通過帯域リップル、阻止帯域リップル、およびカットオフの鋭さを特定の値に設定することによって指定されます。これらは、フィルタの動作を記述する有理楕円関数の最小次数nを指定します。具体的には、 ω _p未満の周波数を通過させ、 ω _sを超える周波数を阻止するフィルタを考え、補助次数を定義します。通過帯域ゲインを[1- G _p、 1]に、阻止帯域ゲインを[0、 1- G _s ]に制限するには、次数n が次式を満たしている必要があります^{[ 1 ]} 。この式は、完全な楕円積分関数 Kの近似から生じます。正確な定式化は^{[ 2 ]}です。もう 1 つの設計上の考慮事項は、フィルタの構築に使用する電子部品の値に対するゲイン関数の感度です。この感度は、フィルタ伝達関数の極の品質係数 ( Q 係数) に反比例します。極の Q 係数は次のように定義されます。 $${\displaystyle u={\frac {1}{2}}\cdot {\frac {1-{\sqrt[{4}]{1-\left({\frac {\omega _{p}}{\omega _{s}}}\right)^{2}}}}{1+{\sqrt[{4}]{1-\left({\frac {\omega _{p}}{\omega _{s}}}\right)^{2}}}}}}$$ $${\displaystyle n\geq {\frac {\log {\left(16{\frac {G_{s}}{G_{p}}}\right)}}{-\log(u+2u^{5}+15u^{9}+150u^{13})}}}$$ $${\displaystyle n\geq {\frac {K\left({\sqrt {1-\tau _{1}^{2}}}\right)K(\tau _{2})}{K(\tau _{1})K\left({\sqrt {1-\tau _{2}^{2}}}\right)}}}$$ $${\displaystyle {\begin{aligned}\tau _{1}&={\sqrt {\frac {G_{p}}{G_{s}}}}\\\tau _{2}&={\frac {\omega _{p}}{\omega _{s}}}\end{aligned}}}$$

${\displaystyle Q=-{\frac {|s_{pm}|}{2\mathrm {Re} (s_{pm})}}=-{\frac {1}{2\cos(\arg(s_{pm}))}}}$

これは、極がゲイン関数に与える影響の尺度です。楕円フィルタの場合、与えられた次数において、リップル係数と選択係数の間に、伝達関数内のすべての極のQ値を同時に最小化する関係が存在します。

${\displaystyle \epsilon _{Qmin}={\frac {1}{\sqrt {L_{n}(\xi )}}}}$

この結果、部品の変動に対する不感性が最大限に高まるフィルタが得られますが、通過帯域と阻止帯域のリップルを個別に指定する機能は失われます。このようなフィルタでは、次数が増加するにつれて、両帯域のリップルが減少し、遮断率が増加します。フィルタ帯域で特定の最小リップルと特定の遮断率を達成するために最小 Q 楕円フィルタを使用する場合は、必要な次数は一般に、最小 Q 制限がない場合に必要な次数よりも大きくなります。ゲインの絶対値の図は、極が楕円ではなく円内に配置されていることを除けば、前のセクションの図と非常によく似ています。極は等間隔ではなく、ω 軸上にゼロが存在します。これは、極が等間隔の円内にゼロなしで配置されている Butterworth フィルタとは異なります。

他の線形フィルタとの比較

以下は、同じ数の係数で得られた他の一般的な種類のフィルターの横に楕円フィルターを示した画像です。

画像から明らかなように、楕円フィルターは他のすべてのフィルターよりもシャープですが、帯域幅全体にリップルが現れます。

チェビシェフ伝達零点からの構築

楕円フィルタの阻止帯域は、本質的には透過零点を持つチェビシェフフィルタであり、透過零点は等リップル阻止帯域を形成するように配置されている。これを踏まえると、分子にチェビシェフ反射零点を含み、分母に透過零点を含まないチェビシェフフィルタの特性方程式を、分子に楕円反射零点を含み、分母に楕円透過零点を含む楕円フィルタに変換することが可能である。これは、チェビシェフ反射零点のスケール逆数から透過零点を反復的に作成し、次に透過零点から等リップルチェビシェフ通過帯域を再構築し、反復によって に対する有意な変化が生じなくなるまで繰り返すことによって行われる。^{[ 3 ]} ここで使用されるスケーリング係数 は、阻止帯域と通過帯域のカットオフ周波数の比であり、「選択係数」の逆数としても知られている。^{[ 1 ]} 楕円設計は一般に阻止帯域減衰要件から指定されるため、は上記の最小次数nを確立する式から導くことができる。 ${\displaystyle K(s)}$ ${\displaystyle K(s)}$ ${\displaystyle K(s)}$ ${\displaystyle \Omega _{c}}$ ${\displaystyle \Omega _{c}}$

比は、上記の最小位数nの問題をnから逆に解いて求めることで導出できる。^{[ 1 ]}${\displaystyle \omega _{s}/\omega _{p}}$ ${\displaystyle \Omega _{c}}$ ${\displaystyle \Omega _{c}}$

${\displaystyle {\begin{aligned}n&={\text{ number of poles (order of the filter)}}\\\alpha _{p}{\text{ and }}\alpha _{s}&={\text{pass band and stop band attenuation, respectively}}\\\omega _{p}{\text{ and }}\omega _{s}&={\text{pass band and stop band frequencies, respectively}}\\D&={\frac {G_{s}}{G_{p}}}={\frac {10^{\alpha _{s}/10}-1}{10^{\alpha _{p}/10}-1}}\\q&=(16D)^{(-1/n)}\\0&=-q+u+2u^{5}+15u^{9}+150u^{13}\\u&={\text{ real root of above equation}}\\k&={\text{ the selectivity factor }}={\sqrt {1-{\bigg (}{\frac {1-2u}{1+2u}}{\bigg )}^{4}}}\\\Omega _{c}&={\frac {\omega _{s}}{\omega _{p}}}={\frac {1}{k}}={\frac {1}{\sqrt {1-{\bigg (}{\frac {1-2u}{1+2u}}{\bigg )}^{4}}}}\\\end{aligned}}}$

減衰要件から計算される特性多項式は、古典的な変換によって伝達関数多項式に変換される。ここで、および は通過帯域リップルである。^{[ 3 ] [ 4 ]}${\displaystyle K(s)}$ ${\displaystyle \Omega _{c}}$ ${\displaystyle G(s)}$ ${\displaystyle G(s)={\sqrt {1/(1+\varepsilon ^{2}K(s)K(-s))}}{\bigg |}_{\text{LHP poles}}}$ ${\displaystyle \varepsilon ^{2}=10^{Ap/10.}-1.}$ ${\displaystyle A_{p}}$

簡単な例

0 ～ 1 rad/秒で通過帯域リップルが 1 dB、少なくとも 1.25 rad/秒から 40 dB の阻止帯域リップルを持つ楕円フィルタを設計します。 ${\displaystyle \infty }$

上記の計算をceil()関数を適用する前の n の値に適用すると、n は 4.83721900 となり、ceil()関数を適用することで次の整数 5 に切り上げられます。つまり、指定された設計要件を満たすには 5 極の楕円フィルタが必要であることを意味します。上記の計算を、正確に 40 dB の減衰を持つ阻止帯域を設計するために必要な n に適用すると、n は 1.2186824 となります。 ${\displaystyle \Omega _{c}}$ ${\displaystyle \Omega _{c}}$

多項式スケール逆変換関数は、各根sを に変換することによって実行できます。これは、図に示すように、多項式を反転し、 でスケール変換することで簡単に実行できます。 ${\displaystyle \Omega _{c}/s}$ ${\displaystyle \Omega _{c}}$

${\displaystyle {\begin{aligned}&as^{n}+bs^{n-1}\dots cs^{2}+ds^{1}+es^{0}\Longrightarrow ({\frac {e}{\Omega _{c}^{n}}})s^{n}+({\frac {d}{\Omega _{c}^{n-1}}})s^{n-1}\dots ({\frac {c}{\Omega _{c}^{2}}})s^{2}+({\frac {b}{\Omega _{c}^{1}}})s^{1}+({\frac {a}{\Omega _{c}^{0}}})s^{0}\\\end{aligned}}}$

楕円の設計手順は次のようになります。^{[ 3 ]}

1. 1 dB の通過帯域リップルを持つチェビシェフ フィルターを設計します。
2. 反射ゼロをすべて反転して透過ゼロを作成する${\displaystyle \Omega _{c}}$
3. チェビシェフ伝送零点法で概説されているプロセスを使用して、伝送零点から等リップル通過帯域を作成する。
4. 通過帯域と阻止帯域の両方で顕著な変化がなくなるまで、手順2と3を繰り返します。通常、15～25回の反復で係数の差は1.e-15程度になります。

手順を説明するために、以下のK(s)方程式は標準的なチェビシェフK(s)から開始し、その後、プロセスを反復します。最初の3回の反復では目に見える違いが見られます。18回の反復に達すると、K(s)の差は無視できるようになります。K(s)係数の変化が設計精度要件を満たすほど十分に小さくなった時点で、反復を中止できます。以下のK(s)反復はすべて正規化されているため、 となっていますが、必要に応じてこの手順を最後の反復まで延期できます。 ${\displaystyle |K(j)|=1}$

${\displaystyle {\begin{aligned}{\text{iteration 0: }}K(s)&={\frac {16s^{5}+20s^{3}+5s}{1}}\\{\text{iteration 1: }}K(s)&={\frac {9.2965947s^{5}+12.999133s^{3}+4.0025668s}{0.14167325s^{4}+0.84164496s^{2}+1}}\\{\text{iteration 2: }}K(s)&={\frac {8.6496472s^{5}+12.270597s^{3}+3.8746611s}{0.19518773s^{4}+0.94147634s^{2}+1}}\\&\vdots \\{\text{iteration 17: }}K(s)&={\frac {8.550086786383502s^{5}+12.157269873073034s^{3}+3.854163602012615s}{0.2043607336740334s^{4}+0.9573802183509494s^{2}+1}}\\{\text{iteration 18: }}K(s)&={\frac {8.550086786383422s^{5}+12.157269873072942s^{3}+3.854163602012599s}{0.2043607336740334s^{4}+0.9573802183509494s^{2}+1}}\\\end{aligned}}}$

伝達関数を求めるには、次の手順に従います。^{[ 3 ]}${\displaystyle G(s)}$

${\displaystyle {\begin{aligned}\varepsilon ^{2}&=10^{1dB/10.}-1.=.25892541\\G(s)&={\sqrt {G(s)G(-s)}}{\bigg |}_{\text{LHP poles}}={\sqrt {\frac {1}{1+\varepsilon ^{2}K(s)K(-s)}}}{\bigg |}_{\text{LHP poles}}={\sqrt {\frac {K(s)_{den}K(-s)_{den}}{K(s)_{den}K(-s)_{den}+\varepsilon ^{2}K(s)_{num}K(-s)_{num}}}}{\bigg |}_{\text{LHP poles}}\\&={\frac {0.20436073s^{4}+0.95738022s^{2}+1}{{\sqrt {-18.928479s^{10}-53.78661s^{8}-54.942632s^{6}-22.939175s^{4}-1.931467+1}}{\bigg |}_{\text{LHP poles}}}}\\\end{aligned}}}$

左半平面からを求めるには、分子と分母を因数分解し、求根アルゴリズムを用いて根を求めます。分母の右半平面にあるすべての根と、分子の重複根の半分を破棄し、残りの根で再構築します。^{[ 3 ] [ 4 ]} 一般に、 では1 に正規化します。 ${\displaystyle G(s)}$ ${\displaystyle G(s)}$ ${\displaystyle |G(s)|}$ ${\displaystyle s=0}$

${\displaystyle {\begin{aligned}&G(s)={\frac {0.20436073s^{4}+0.95738022s^{2}+1}{4.3506872s^{5}+4.0174213s^{4}+8.0362343s^{3}+4.9129149s^{2}+3.4288915s+1}}\\\end{aligned}}}$

例が正しいことを確認するために、通過帯域リップルが1 dB、カットオフ周波数が1 rad/sec、ストップバンド減衰が1.21868 rad/secから40 dBの場合のプロットを以下に示し ます。${\displaystyle G(s)}$ ${\displaystyle G(s)}$ ${\displaystyle j\omega }$

5極楕円シミュレーション

注文の変更も

偶数次楕円フィルタは、通常はインダクタ、コンデンサ、伝送線路などの受動素子と、両側に等しい値の終端で実装されますが、結合コイルを使用せずに従来の楕円伝達関数で実装することはできず、これは望ましくなかったり実現不可能であったりする場合があります。これは、偶数次チェビシェフ反射零点と透過零点を物理的に収容できないためです。この零点によって散乱行列S12 値が で S12 値を超え、 で有限の S12 値が存在することになります。通過帯域 S12 に合わせるために終端の 1 つを増減してフィルタを設計することが実現不可能な場合は、楕円伝達関数を変更して、通過帯域と阻止帯域の等リップル応答を維持しながら、最低偶数次反射零点を に、最高偶数次透過零点を に移動する必要があります。 ^{[ 5 ]}${\displaystyle \omega =0}$ ${\displaystyle \omega =\infty }$ ${\displaystyle \omega =0}$ ${\displaystyle \omega =\infty }$

必要な修正は、楕円伝達関数の各極と零点を、最も低い周波数の反射零点を零点に、最も高い周波数の透過零点を に、そして残りの極と零点を等リプル通過帯域と阻止帯域を維持するために必要な値にマッピングすることです。最も低い周波数の反射零点は分子を因数分解することで見つけることができ、最も高い周波数の透過零点は分母を因数分解することで見つけることができます。 ${\displaystyle \infty }$ ${\displaystyle K(s)}$ ${\displaystyle K(s)}$

反射零点を変換する場合、次の式を のすべての極と零点に適用します。^{[ 5 ]} 理論上は変換操作は または のいずれかで実行できますが、反射零点は から抽出する必要があるため、 で変換操作を実行する方が通常は効率的です。 ${\displaystyle K(s)}$ ${\displaystyle K(s)}$ ${\displaystyle G(s)}$ ${\displaystyle K(s)}$ ${\displaystyle K(s)}$

${\displaystyle R_{i}'={\sqrt {\frac {R_{i}^{2}+\omega _{LO}^{2}}{1-\omega _{LO}^{2}}}}}$

どこ：

${\displaystyle R_{i}}$ 元の楕円関数の零点または極

${\displaystyle R_{i}'}$ 修正偶数次伝達関数のマッピングされた零点または極です。

${\displaystyle \omega _{LO}}$ 通過帯域内の最も低い周波数の反射ゼロです。

の虚数部の符号は、元の の符号によって決まります。 ${\displaystyle R_{i}'}$ ${\displaystyle R_{i}}$

透過零点を変換する場合、次の式を のすべての極と零点に適用します。^{[ 5 ]} 理論的には、変換操作は または のどちらでも実行できますが、反射零点を から抽出する必要がある場合は、 で変換操作を実行する方が効率的です。 ${\displaystyle K(s)}$ ${\displaystyle K(s)}$ ${\displaystyle G(s)}$ ${\displaystyle K(s)}$ ${\displaystyle K(s)}$

${\displaystyle R_{i}'={\sqrt {\frac {(\omega _{HI}^{2}-1)R_{i}^{2}}{\omega _{HI}^{2}+R_{i}^{2}}}}}$

どこ：

${\displaystyle R_{i}}$ 元の楕円関数の零点または極

${\displaystyle R_{i}'}$ 修正偶数次伝達関数のマッピングされた零点または極です。

${\displaystyle \omega _{HI}}$ 通過帯域内の最も高い周波数の伝送ゼロです。

の虚数部の符号は、元の の符号によって決まります。 を操作する場合、左半平面の要件を満たすために、 の実数部の符号は負でなければなりません。${\displaystyle R_{i}'}$ ${\displaystyle R_{i}}$ ${\displaystyle G(s)}$ ${\displaystyle R_{i}'}$

すべてのアプリケーションにおいて、通過帯域と阻止帯域の両方の変換が必要であることに注意することが重要です。例えば、パッシブネットワークダイプレクサは、偶数次の阻止帯域の変換のみを必要とし、通過帯域を変換せずに偶数次の阻止帯域のみを変換することで、より効率的に動作します。^{[ 5 ]}

が完了すると、における S12 の散乱行列値が 1および 0 である等リプル伝達関数が作成されます。これは、受動的な等終端ネットワークで実装できます。 ${\displaystyle G(s)}$ ${\displaystyle \omega =0}$ ${\displaystyle \omega =\infty }$

下の図は、等リプル通過帯域と阻止帯域の周波数応答を維持しながら、最低周波数反射ゼロを有限周波数から 0 に、最高周波数透過ゼロを に移動することにより、偶数次の均等終端受動ネットワークをサポートするように変更された 8 次楕円フィルタを示しています。 ${\displaystyle \infty }$

偶数次修正楕円の図

上記の楕円関数構築の段落における および 次数の計算は、修正されていない楕円関数のみを対象としています。偶数次数の変更は通過帯域または阻止帯域の減衰には影響しませんが、次数と計算には小さな誤差が生じることが予想されます。したがって、通過帯域および阻止帯域の減衰を維持したい場合は、すべての反復処理が完了した後に偶数次数の変更を適用することが重要です。要件に基づいて偶数次数変更楕円関数を作成した場合、実際の は設計値よりもわずかに大きくなります。同様に、次数n の計算では、実際に必要な次数よりも小さい値になる場合があります。 ${\displaystyle \Omega _{c}}$ ${\displaystyle \Omega _{c}}$ ${\displaystyle K(s)}$ ${\displaystyle \Omega _{c}}$ ${\displaystyle \Omega _{c}}$ ${\displaystyle \Omega _{c}}$

砂時計の実装

砂時計フィルタは、反射零点が1 rad/秒の3.01 dB正規化カットオフ減衰周波数の周りの透過零点の逆数である特殊なフィルタであり、フィルタのすべての極が単位円上に存在する。^{[ 6 ]} 楕円砂時計フィルタの実装は、通過帯域がより平坦であるという点で逆チェビシェフフィルタよりも優れており、グループ遅延がカットオフ周波数でそれほど急峻ではないという点で従来の楕円フィルタよりも優れています。

7極砂時計型レシプロカルS11およびS12周波数応答

合成プロセス

砂時計フィルタを合成する最も簡単な方法は、指定された設計阻止帯域減衰量Asと、散乱パラメータというロスレス2ポートネットワーク要件を満たす計算された通過帯域減衰量を持つ楕円フィルタを設計することです。^{[ 7 ]} よく知られている振幅dBから算術変換、とともに代数操作により、次の通過帯域減衰計算要件が得られます。 ${\displaystyle |S_{11}|^{2}+|S_{12}|^{2}=1}$ ${\displaystyle (S_{ij})_{dB}=20log_{10}(|S_{ij}|_{arith})}$

${\displaystyle A_{p}=-10\log _{10}{(1.-10^{(-A_{s}/10)})}}$

上記で定義したA _pは、まだ未知の3.01 dBのカットオフ周波数付近で、相互反射および透過零点を生成します。通過帯域周波数が1 rad/secの楕円フィルタを設計するには、3.01 dBの減衰周波数を決定し、その周波数を用いて楕円設計多項式を逆スケーリングする必要があります。その結果、正規化周波数1 rad/secで減衰量が3.01 dBとなる多項式が得られます。 ニュートン法、または解探索アルゴリズムを用いて方程式を直接解くことで、3.01 dBの減衰周波数を決定できます。

ニュートン法による周波数スケーリング

が砂時計伝達関数で3.01 dBの周波数を見つけ、が3 dBの周波数を見つける場合、以下の手順で見つけることができます。${\displaystyle G(s)}$ ${\displaystyle \omega _{c}}$ ${\displaystyle \omega _{c}}$

1. がまだ利用できない場合は、 を掛けてを取得します。${\displaystyle G(s)G(-s)}$ ${\displaystyle G(s)}$ ${\displaystyle G(-s)}$ ${\displaystyle G(s)G(-s)}$
2. が で割り切れる場合、 のすべての項を反転します。つまり、、、 などとなります。修正された関数は と呼ばれ 、この修正により、多項式とその導関数を評価する際に、複素数ではなく実数を使用できるようになります。これで、複素数の代わりに実数が使用できるようになります。${\displaystyle s^{n}}$ ${\displaystyle (n+2)}$ ${\displaystyle 4}$ ${\displaystyle s^{2}}$ ${\displaystyle s^{6}}$ ${\displaystyle s^{10}}$ ${\displaystyle G_{2}(s)G_{2}(-s)}$ ${\displaystyle \omega _{a}}$ ${\displaystyle j\omega _{a}}$
3. を用いて、 dB単位の所望の減衰量 を二乗算術ゲイン値 に変換します。例えば、3.010 dBは0.5に、1 dBは0.79432823に変換されます。${\displaystyle A_{dB}}$ ${\displaystyle B_{arith}^{2}}$ ${\displaystyle B_{arith}^{2}=10^{A_{dB}/10}}$
4. 実数値 を用いてニュートン法で修正 を計算します。常に絶対値を取ります。${\displaystyle |G_{2}(s)G_{2}(-s)|}$ ${\displaystyle \omega _{a}}$
5. 実数値 に関して修正された導関数を計算します。導関数の絶対値は取らないでください。${\displaystyle G_{2}(\omega _{a})G_{2}(-\omega _{a})}$ ${\displaystyle \omega _{a}}$

手順 1) から 4) が完了すると、ニュートン法を含む式は次のように表すことができます。

${\displaystyle \omega _{a}=\omega _{a}-([G_{2}(\omega _{a})G_{2}(-\omega _{a})|-B^{2})/(d[G_{2}(\omega _{a})G_{2}(-\omega _{a})]/d\omega _{a})}$

複雑な演算を必要とせず、実数値を用いて を計算できます。信頼性を高めるため、反復計算の初期段階で の動きが負にならないように制限する必要があります。収束が完了すると、を に使用して、元の伝達関数の分母をスケーリングできます。修正された の減衰量は、1 rad/sec でほぼ正確な目標値になります。適切に実行すれば、小次フィルタと非常に大次フィルタの両方において、広範囲にわたる目標減衰値を設定するのに数回の反復計算で済みます。 ${\displaystyle \omega _{a}}$ ${\displaystyle \omega _{a}}$ ${\displaystyle \omega _{a}}$ ${\displaystyle \omega _{c}}$ ${\displaystyle G(s)}$ ${\displaystyle G(s)}$

根探索による周波数スケーリング

には位相情報が含まれていないため、伝達関数を直接因数分解しても有用な結果は得られません。しかし、伝達関数を で乗じて の奇数乗をすべて消去することで修正できます。これにより、 はすべての周波数において実数となり、その後、目的の注目度の2乗となる周波数を求めることができます。 ${\displaystyle |G(j\omega _{a})|}$ ${\displaystyle G(-s)}$ ${\displaystyle G(j\omega a)}$ ${\displaystyle G(j\omega a)}$

1. がまだ利用できない場合は、 を掛けてを取得します。${\displaystyle G(s)G(-s)}$ ${\displaystyle G(s)}$ ${\displaystyle G(-s)}$ ${\displaystyle G(s)G(-s)}$
2. を用いて、 dB単位の所望の減衰量 を二乗算術ゲイン値 に変換します。例えば、3.010 dBは0.5に、1 dBは0.79432823に変換されます。${\displaystyle A_{dB}}$ ${\displaystyle B_{arith}^{2}}$ ${\displaystyle B_{arith}^{2}=10^{A_{dB}/10}}$
3. 探す${\displaystyle P(S)=G_{num}(S)G_{num}(-S)-B_{arith}^{2}G_{den}(S)G_{den}(-S)}$
4. 根を求めるアルゴリズムを使用して P(S) の根を求めます。
5. 上記の根の集合から、すべての次数フィルタの正の虚数根と、偶数次フィルタの正の実数根を選択します。${\displaystyle \omega _{c}}$

伝達関数のスケーリング

が決定された場合、砂時計伝達関数多項式は次のようにスケーリングできます。 ${\displaystyle \omega _{c}}$

${\displaystyle {\begin{aligned}G(s)_{\text{orig}}&={\frac {N_{nn}s^{nn}+\dots +N_{2}s^{2}+N_{1}s+N_{0}}{D_{nn}s^{nd}+\dots +D_{2}s^{2}+D_{1}s+D_{0}}}{\text{ (original unscaled transfer function polynomials)}}\\G(s)_{\text{ scaled}}&={\frac {N_{nn}(\omega _{c}s)^{nn}+\dots +N_{2}(\omega _{c}s)^{2}+N_{1}\omega _{c}s+N_{0}}{D_{nn}(\omega _{c}s)^{nd}+\dots +D_{2}(\omega _{c}s)^{2}+D_{1}\omega _{c}s+D_{0}}}{\text{ (}}3.01{\text{ dB at 1 rad/sec scaled transfer function polynomials)}}\\\omega _{c}&=|G(s)|{\text{ 3.01 dB attenuation frequency }}\\nn,nd&={\text{ order of numerator and denominator, respectively}}\\N,D&={\text{ coefficients of numerator and denominator, respectively}}\\\end{aligned}}}$

注文の変更も

偶数次砂時計フィルタは、他の楕円フィルタと同様に、等終端受動回路に関して同様の制限があります。楕円フィルタの問題を解決するのと同じ偶数次修正が、砂時計フィルタの問題も解決します。

参考文献

1. Rorabaugh, C. Britton (1993年1月1日).デジタルフィルタ設計者ハンドブック（復刻版）. ブルーリッジサミット, ペンシルバニア州, 米国: Tab Books, Division of McGraw-Hill, Inc. pp. 93-95. ISBN 978-0830644315。
2. Paarmann, Larry D. (2001).アナログフィルタの設計と解析：信号処理の観点から. 米国マサチューセッツ州ノーウェル：Kluwer Academic Publishers. pp. 182～198. ISBN 0-7923-7373-1。
3. Dr. Byron Bennett's filter design lecture notes, 1985, Montana State University Archived 2024-04-23 at the Wayback Machine , EE Department Archived 2023-03-28 at the Wayback Machine , Bozeman , Montana, US
4. Sedra, Adel S.; Brackett, Peter O. (1978). 『フィルタ理論と設計：能動フィルタと受動フィルタ』Beaverton, Oegon, US: Matrix Publishers, Inc. pp.  45– 73. ISBN 978-0916460143。
5. ザール、ルドルフ (1979 年 1 月)。フィルター設計ハンドブック(英語およびドイツ語) (第 1 版)。ミュンヘン、ドイツ: Allgemeine Elektricitais-Gesellschaft。pp . 25、26、56–61、116、117。ISBN​ 3-87087-070-2。
6. Bennett, Byron J. (1988年12月). 「新しいフィルタ合成技術 - 砂時計」 . IEEE Transactions on Circuits and Systems . 35 (12): 1469– 1477. doi : 10.1109/31.9910 – IEEE経由.
7. Matthaei, George L.; Young, Leo; Jones, EMT (1984).マイクロ波フィルタ、インピーダンス整合回路、および結合構造. 610 Washington Street, Dedham, Massachusetts, US: Artech House, Inc. (1985年出版). p. 44. ISBN 0-89006-099-1。

• ダニエルズ、リチャード・W. (1974). 『電子フィルタ設計のための近似法』 ニューヨーク：マグロウヒル. ISBN 0-07-015308-6。
• Lutovac, Miroslav D.; Tosic, Dejan V.; Evans, Brian L. (2001). MATLABとMathematicaを用いた信号処理のためのフィルタ設計. ニュージャージー州, 米国: Prentice Hall. ISBN 0-201-36130-2。