チェビシェフフィルタ

チェビシェフフィルタは、バターワースフィルタよりも急峻なロールオフ特性を持つアナログまたはデジタルフィルタであり、通過帯域リップル(タイプI)または阻止帯域リップル(タイプII)のいずれかを持ちます。チェビシェフフィルタは、フィルタの動作周波数範囲全体にわたって、理想的なフィルタ特性と実際のフィルタ特性との誤差を最小化する特性を持ちますが[1] [2] 、この特性は周波数応答にリップルを伴います。このタイプのフィルタは、その数学的特性がチェビシェフ多項式に由来することから、パフヌティ・チェビシェフにちなんで名付けられました。タイプIチェビシェフフィルタは通常「チェビシェフフィルタ」と呼ばれ、タイプIIフィルタは通常「逆チェビシェフフィルタ」と呼ばれます[3] 。チェビシェフフィルタに固有の通過帯域リップルのため、特定の用途では、通過帯域ではより滑らかな応答を示し、阻止帯域ではより不規則な応答を示すフィルタが好まれます[4] 。

I型チェビシェフフィルタ(チェビシェフフィルタ)

4次I型チェビシェフローパスフィルタの周波数応答は、

タイプIチェビシェフフィルタは、最も一般的なチェビシェフフィルタです。次ローパスフィルタ角周波数の関数としてのゲイン(または振幅)応答 は、における伝達関数の絶対値に等しくなります。

ここで、 はリップル係数、カットオフ周波数、は 番目の次数チェビシェフ多項式です

通過帯域は等リップル挙動を示し、リップル係数 によってリップルが決定されます。通過帯域では、チェビシェフ多項式は -1 と 1 の間を交互に変化するため、フィルタゲインは で最大値、 で最小値の間を交互に変化します

リップル係数εは、デシベル単位の通過帯域リップルδと次の関係があります。

カットオフ周波数ではゲインは再び値を持ちますが、周波数が上昇するにつれて 阻止帯域まで低下し続けます。この挙動は右の図に示されています。カットオフ周波数を-3dBに定義するという一般的な方法は、チェビシェフフィルタでは通常適用されません。代わりに、ゲインが最終的にリップルの値まで低下する点をカットオフとします。

3 dB 周波数は次の関係にあります

チェビシェフ フィルタの次数は、アナログ エレクトロニクスを使用してフィルタを実現するために必要なリアクタンス部品 (たとえば、インダクタ)の数に等しくなります

阻止帯域にリップルを許容し、複素平面の - 軸上に零点を許容することで、さらに急峻なロールオフ特性を得ることができます。これにより、零点およびその近傍においてほぼ無限大の抑制効果が得られます(ただし、部品のQ値、寄生成分、および関連要因によって制限されます)。一方、阻止帯域全体の抑制効果は減少します。この結果は楕円フィルタ、またはカウアーフィルタと呼ばれます。

極と零点

複素周波数空間s  =  σ  +  における8次チェビシェフI型フィルタのゲインの絶対値の対数(ε = 0.1、ω = 0.2 )です。白い点は極であり、σ方向で0.3836...、ω方向で1.071...の半軸を持つ楕円上に配置されています。伝達関数の極は、左半平面にある極です。黒はゲインが0.05以下、白はゲインが20以上に対応します。

簡略化のため、カットオフ周波数は1であると仮定します。チェビシェフフィルタのゲイン関数の極は、ゲイン関数の分母の零点です。複素周波数を用いると、これらは以下の場合に発生します。

チェビシェフ多項式の三角関数の定義と使用により次の式が得られます。

解決する

ここで、逆余弦関数の複数の値は、整数インデックスを用いて明示的に表されます。チェビシェフゲイン関数の極は次のようになります。

三角関数と双曲線関数の性質を利用すると、これを明示的に複素形式で書くことができます。

どこで そして

これは、 の媒介変数方程式として見ることができ、を中心とする空間の楕円上に極があり、実半軸の長さが で、虚半軸の長さが であることを示しています。

伝達関数

上記の式はゲインの極を表す。各複素極に対して、複素共役極が1つ存在し、各共役対に対して、その負の極がさらに2つ存在する。伝達関数は安定でなければならない。つまり、その極はゲインの極のうち負の実部を持つ極であり、したがって複素周波数空間の左半平面に位置する。伝達関数は次のように表される。

ここでは、上記の式から得られた、実数項の前に負の符号が付いたゲインの極のみを示します。

群遅延

ε = 0.5 の 5 次 I 型チェビシェフ フィルタのゲインとグループ遅延。

群遅延は、角周波数に対する位相の微分として定義されます

ε=0.5の5次I型チェビシェフフィルタのゲインと群遅延が左のグラフに示されています。阻止帯域にはリップルがありません。しかし、通過帯域における群遅延のリップルは、異なる周波数成分に異なる遅延があることを示しており、これが通過帯域におけるゲインのリップルと相まって波形の歪みを引き起こします。

注文の変更も

偶数次チェビシェフフィルタは、通常はインダクタ、コンデンサ、伝送線路などの受動素子を用いて、両側に等しい値の終端を施して実装されますが、従来のチェビシェフ伝達関数では結合コイルを使用せずに実装することはできません。これは、特に高周波においては望ましくなく、実現不可能な場合もあります。これは、偶数次チェビシェフ反射零点を物理的に収容できないことに起因し、散乱行列のS12 値が における S12 値を超えることになります。通過帯域 S12 に対応するために終端の 1 つを増減させてフィルタを設計することが実現不可能な場合は、通過帯域の等リップル応答を維持しながら、最も低い偶数次反射零点を に移動するようにチェビシェフ伝達関数を変更する必要があります。[5]

必要な修正は、チェビシェフ伝達関数の各極を、最低周波数の反射零点を零点にマッピングし、残りの極を等リップル通過帯域を維持するために必要な位置にマッピングすることです。最低周波数の反射零点は、チェビシェフノード、から見つけることができます。完全なチェビシェフ極マッピング関数を以下に示します。[5]


どこ:

nはフィルタの次数です(偶数である必要があります)

Pは伝統的なチェビシェフ伝達関数の極である

P' は、修正偶数次伝達関数のマッピングされた極です。

「左半平面」は、負の実数値を含む平方根を使用することを示します。

完了すると、等リップル受動回路網を実装する場合、S12の反射零点散乱行列値が1、S11の反射零点散乱行列値が0となる置換等リップル伝達関数が作成されます。下の図は、等リップル通過帯域周波数応答を維持しながら、最低周波数反射零点を有限周波数から0に移動することで、偶数次の等リップル受動回路網をサポートするように修正された8次チェビシェフフィルタを示しています。

偶数次修正チェビシェフの図
偶数次修正チェビシェフの図

カウアートポロジーにおけるLC素子の値の公式は、偶数次修正チェビシェフ伝達関数には適用できず、使用できません。したがって、LC値は、反射係数から導出できるインピーダンス関数の従来の連分数から計算する必要があります。反射係数は、伝達関数から導出できます。

最小注文数

必要最小限の要素数でチェビシェフフィルタを設計するには、チェビシェフフィルタの最小次数は次のように計算できます。[6] これらの式は標準的なローパスチェビシェフフィルタのみを考慮しています。次数の変更や有限の阻止帯域透過零点によっても、これらの式では考慮されていない誤差が生じます。

どこ:

通過帯域リップル周波数と最大リップル減衰量(dB)である。

およびは阻止帯域周波数とその周波数における減衰量(dB)である。

最小の極数、つまりフィルタの次数です。

ceil []は次の整数への切り上げ関数です。

カットオフ減衰の設定

チェビシェフフィルタの通過帯域カットオフ減衰量は、通常、上記の計算によって設定される通過帯域リップル減衰量と同じである。しかし、ダイプレクサやトリプレクサ[5]などの多くの用途では、必要な反射を得るために-3.0103 dBのカットオフ減衰量が必要となる。その他の特殊な用途では、様々な理由から、カットオフ減衰量に異なる特定の値が必要となる場合がある。したがって、チェビシェフフィルタの通過帯域カットオフ減衰量を、通過帯域リップル減衰量とは独立して、-1 dB、-10 dBなどに設定できる手段があると便利である。カットオフ減衰量は、伝達関数の極を周波数スケーリングすることによって設定できる。

スケーリング係数は、定義するチェビシェフフィルタ関数 (および を含む)の直接的な代数操作によって決定できます。チェビシェフ関数の一般定義が必要です。これは、チェビシェフ多項式方程式、および逆チェビシェフ関数から導出できます。 の値に対して実数値を維持するために複素双曲恒等式を用いて方程式を およびと書き直すことができます

上記の方程式と参照に簡単な代数を使用すると、各チェビシェフ極をスケーリングする式は次のようになります。

どこ:

希望するカットオフ減衰を設定するために再配置されたポールです。

楕円上にあるリップルカットオフポールです。

通過帯域減衰リップル(dB 単位、0.05 dB、1 dB など)です。

カットオフ周波数における所望の通過帯域減衰量(dB 単位、1 dB、3 dB、10 dB など)

極の数(フィルタの次数)です。

通過帯域リップル減衰を通過帯域カットオフ減衰として使用して上記の式を簡単に検証すると、この場合、極調整は予想どおり 1.0 になることがわかります。

偶数次修正カットオフ減衰調整

受動等終端フィルタの通過帯域リップルを偶数次用に修正したチェビシェフフィルタを設計する場合、減衰周波数の計算には、計算された減衰周波数に対して偶数次調整演算を実行することにより、偶数次調整を含める必要があります。これにより、周波数を実変数(この場合は )として扱うことができるため、偶数次調整演算が若干簡素化されます


どこ:

希望するカットオフ減衰を設定するために再配置されたポールです。

偶数次通過帯域用に修正されたリップル カットオフ ポールです。

通過帯域減衰リップル(dB 単位、0.05 dB、1 dB など)です。

カットオフ周波数における所望の通過帯域減衰量(dB 単位、1 dB、3 dB、10 dB など)

極の数(フィルタの次数)です。

最小の偶数次チェビシェフノードである

タイプIIチェビシェフフィルタ(逆チェビシェフフィルタ)

5次II型チェビシェフローパスフィルタの周波数応答は、

逆チェビシェフフィルタとも呼ばれるタイプIIチェビシェフフィルタは、タイプIほど急激に減衰せず、部品点数も増えるため、あまり一般的ではありません。通過帯域にはリップルがありませんが、阻止帯域には等リップルがあります。ゲインは以下のとおりです。

阻止帯域では、チェビシェフ多項式は-1と1の間で振動し、ゲインは0と1の間で振動する。

そして、この最大値に達する最小の周波数がカットオフ周波数である。したがって、パラメータεはデシベル単位の阻止帯域減衰量γと以下の関係にある。

阻止帯域減衰量が5 dBの場合、ε = 0.6801、10 dBの場合、ε = 0.3333です。周波数f 0 = ω 0 /2 πはカットオフ周波数です。3 dB周波数f Hはf 0と以下の関係があります

極と零点

複素周波数空間(s = σ + jω)における8次チェビシェフII型フィルタのゲインの絶対値の対数。ε = 0.1 、。白い点は極、黒い点は零点である。16個の極がすべて示されている。各零点は多重度2で、12個の零点が示されている。そのうち4個は図の外側、正のω軸上に2個、負のω軸上に2個配置されている。伝達関数の極は左半平面上の極であり、伝達関数の零点は多重度1の零点である。黒はゲインが0.05以下、白はゲインが20以上である。

カットオフ周波数が1であると仮定すると、チェビシェフフィルタのゲインの極はゲインの分母のゼロになります。

タイプ II チェビシェフ フィルタのゲインの極は、タイプ I フィルタの極の逆になります。

ここで、タイプIIチェビシェフフィルタの零点はゲインの分子の零点である。

したがって、タイプ II チェビシェフ フィルタの零点は、チェビシェフ多項式の零点の逆数になります。

のために

伝達関数

伝達関数はゲイン関数の左半平面の極によって与えられ、同じゼロを持ちますが、これらのゼロは二重ゼロではなく単一ゼロです。

群遅延

ε = 0.1 の 5 次タイプ II チェビシェフ フィルタのゲインとグループ遅延。

ε=0.1の5次タイプIIチェビシェフフィルタのゲインと群遅延を左のグラフにプロットしています。阻止域ではゲインにリップルがありますが、通過域ではリップルがないことがわかります。

注文の変更も

チェビシェフ偶数次フィルタと同様に、標準チェビシェフ II 偶数次フィルタは結合コイルを使用せずに等終端受動素子で実装することはできず、これは望ましくなかったり実現不可能である可能性があります。チェビシェフ II の場合、これは阻止帯域における S12 の有限減衰によるものです。[5] ただし、チェビシェフ II 阻止帯域の等リプル関数を維持しながら、最高周波数の有限伝送零点を無限大に変換することで、偶数次チェビシェフ II フィルタを変更できます。この変換を行うには、標準チェビシェフ関数の代わりに偶数次修正チェビシェフ関数を使用して、偶数次修正チェビシェフ II 伝達関数を作成するために必要なチェビシェフ II 極を定義します。零点は、偶数次修正チェビシェフ多項式の根、つまり偶数次修正チェビシェフ ノードを使用して作成されます。

下の図は、等リップル停止帯域周波数応答を維持しながら、最高周波数伝送ゼロを有限周波数からに移動することにより、偶数次の均等に終端された受動ネットワークをサポートするように変更された 8 次逆チェビシェフ フィルタを示しています。

偶数次修正逆チェビシェフの図
偶数次修正逆チェビシェフの図

最小注文数

逆チェビシェフフィルタを必要最小限の要素数で設計するには、逆チェビシェフフィルタの最小次数は以下のように計算できます。[7] これらの式は標準的なローパス逆チェビシェフフィルタのみを考慮しています。次数を変更しても、これらの式では考慮されていない誤差が生じます。この式はチェビシェフフィルタの最小次数に使用されている式と同じですが、変数の定義が若干異なります。

どこ:

通過帯域周波数とその周波数における減衰量(dB)である。

およびは阻止帯域周波数と最小阻止帯域減衰量(dB)である。

最小の極数、つまりフィルタの次数です。

ceil []は次の整数への切り上げ関数です。

カットオフ減衰の設定

前述の標準的なカットオフ減衰量は、通過帯域リップル減衰量と同じです。しかし、チェビシェフフィルタと同様に、カットオフ減衰量を任意の値に設定することが有用であり、その理由も同じです。チェビシェフIIフィルタのカットオフ減衰量の設定は、チェビシェフフィルタの場合と同じです。ただし、算術減衰量とリップル成分が式の中で反転され、極と零点が結果に乗算されます(チェビシェフフィルタの場合は除算)。

偶数次修正カットオフ減衰調整

チェビシェフ偶数次修正カットオフ減衰に使用されたのと同じ、極と零点に対する偶数次調整は、チェビシェフ II の場合にも使用できます。極と零点の両方が結果に掛けられます。

実装

カウアー位相

パッシブLCチェビシェフ・ローパスフィルタは、カウアートポロジを用いて実現できる。次チェビシェフプロトタイプフィルタのインダクタンスまたはコンデンサの値は、以下の式から計算できる。[8]

G 1、 G kはコンデンサまたはインダクタ要素の値です。3 dB周波数f Hは次のように計算されます。

係数AγβA k、およびB kは次の式から計算できます。

ここで、通過帯域リップルはデシベル単位です。数値は正確な値から切り上げられます

カウアー位相を用いたローパスフィルタ

計算されたG k値は、右図に示すようにシャントコンデンサと直列インダクタに変換できます。また、直列コンデンサとシャントインダクタに変換することもできます。例えば、

  • C 1 シャント= G 1L 2 シリーズ= G 2、...

または

  • L 1 シリーズ= G 1C 2 シャント= G 2、...

G 1がシャントコンデンサまたは直列インダクタの場合、G 0 はそれぞれ入力抵抗またはコンダクタンスに対応することに注意してください。G n+1とG nについても同様の関係が成り立ちます。結果として得られる回路は正規化ローパスフィルタです。周波数変換インピーダンススケーリングを用いることで、この正規化ローパスフィルタは、任意のカットオフ周波数または帯域幅を持つハイパスフィルタバンドパスフィルタバンドストップフィルタに変換できます

デジタル

ほとんどのアナログフィルタと同様に、チェビシェフフィルタは双線形変換によってデジタル(離散時間)再帰型に変換できます。ただし、デジタルフィルタは帯域幅が有限であるため、変換後のチェビシェフフィルタの応答形状は歪んでしまいます。代わりに、応答形状を歪ませないマッチングZ変換法を使用することもできます。

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

次の図は、チェビシェフ フィルターと、同じ数の係数 (5 次) で得られる他の一般的なフィルター タイプを並べて示しています。

チェビシェフ フィルタはバターワース フィルタよりも鋭いです。楕円フィルタほど鋭くはありませんが、帯域幅全体にわたってリップルが少なくなります。

チェビシェフフィルタの高度なトピック

チェビシェフフィルタの設計柔軟性は、このセクションで説明するより高度な設計手法によって強化される可能性があります。 透過零点を阻止帯域に挿入して、特定の不要な周波数を中和したり、カットオフ減衰量を増やしたり、軸外に挿入してより望ましい群遅延を実現したりできます。周波数非対称の設計要件をより効率的に満たすために、通過帯域の両側に異なる数の極を持つ非対称チェビシェフバンドパスフィルタを作成できます。チェビシェフフィルタの特徴である等リップル通過帯域は、通過帯域の一部のみを等リップルにするという設計要件をより効率的に満たすために、通過帯域の一定の割合に制限される場合があります。[9]

チェビシェフ伝達零点

チェビシェフフィルタは、等リップル通過帯域を維持しながら、阻止帯域に任意の位置に配置された有限の透過零点を持つように設計できます。軸に沿った阻止帯域零点は、通常、不要な周波数を除去するために使用されます。実軸に沿った阻止帯域零点、または複素平面における4つの阻止帯域零点は、群遅延をより望ましい形状に修正するために使用できます。透過零点の設計では、特性多項式K(S)を用いて透過零点と反射零点を配置し、これらの零点を用いて伝達関数[10]を作成します

K(S)の計算は、以下の観察された等式に基づいています。[10]

すべての、虚共役ペア、四重共役ペア 、または実数の反対符号ペア に対して。

通過帯域()内で振幅が常に 1 であるとすると、有理項と無理項は 0 から 1 の間で変化する必要があります。したがって、有理項のみを使用して特性関数を作成すると、通過帯域内で等リプル応答が予想され、すべての で特性極(伝送零点)が予想されます

上記の式を使用したK(S)の設計プロセスは以下のとおりです。

実数虚数のペアには正の解を用います。複素数4つ組には、正の実数と共役虚数の解を用います

必要に応じて、 となるように正規化する必要があります

「有理項のみ」とは、積の有理数部を保持し、無理数部を破棄することを意味します。有理数部は、多項式演算を手動で実行することで得られるほか、以下のショートカット(多項式演算から導出され、二項係数を使用する解)を使用して得ることができます。このアルゴリズムは、二項係数が事前に計算された値の参照テーブルから実装されている場合、非常に効率的です。

すべての M の値が 1 に設定されている場合、は標準のチェビシェフ方程式になります。これは、すべての伝達零点が であるため予想されます。偶数次の有限伝達零点チェビシェフ フィルタには、等終端受動ネットワークを使用して構築できないという、全極の場合と同じ制限があります。同じ偶数次の修正を偶数次の特性多項式 に加えることで、等終端受動ネットワークの実装が可能になります。ただし、偶数次の修正により、有限伝達零点もわずかに動きます。この動きは、最低のチェビシェフ ノードを使用して、偶数次の修正の逆で伝達零点を提案することで、大幅に緩和される可能性があります

単純な伝送ゼロの例

1 dB の通過帯域、2 rad/sec の透過ゼロ、次の透過ゼロを持つ 3 極チェビシェフ フィルタを設計します


伝達関数を求めるには、次の手順に従います。[10] [11]


左半平面から求めるには、分子と分母を因数分解して根を求めます。分母の右半平面にあるすべての根と、分子の重複根の半分を捨て、残りの根で再構築します。一般に、 で1 に正規化します

例が正しいことを確認するために、通過帯域リップルが 1 dB、カットオフ周波数が 1 rad/秒、ストップ帯域ゼロが 2 rad/秒である のプロットを以下に示します

2 rad/sec でのチェビシェフ伝達零点

非対称バンドパスフィルタ

チェビシェフ帯域通過フィルタは、上記チェビシェフ透過零点方程式のより一般化された形[10]と下記を用いて、零点と無限大に所望の数の透過零点を配置することにより、幾何学的に非対称な周波数応答を持つように設計することができる。下記の方程式は、1から の周波数正規化された通過帯域を考慮している。零点における透過零点の数が における透過零点の数と同じでない場合、フィルタは幾何学的に非対称となる。有限の透過零点が幾何学的中心周波数(この場合は )の周りに対称的に配置されていない場合も、フィルタは非対称となる。非対称方程式が使用可能な結果を​​生成するためには、フィルタが正味偶数次、すなわちすべての極の合計が偶数でなければならないという制約がある。実数および複素数の4つ組透過零点もこの技術を使用して作成することができ、ローパスフィルタの場合と同様に群遅延応答を修正するのに有用である。非対称チェビシェフバンドパスフィルタを作成するための特性方程式の導出を以下に示します。

必要に応じて、 となるように正規化する必要があります

単純な非対称の例

1~2 rad/秒の1dB通過帯域リップル、に1つの透過零点、0に3つの透過零点を持つ非対称チェビシェフフィルタを設計します。上記の式に数値を適用すると、特性多項式は次のように計算できます。

無理数部分を破棄し、s=j で 1 に正規化すると、

ローパスの場合と同じ手順でからを求め、定数を使って大きさをスケーリングします。[10] [11]

左半平面の極から分母を再構成する際には、反射零点が0dBで発生するように振幅を設定する必要があります。そのためには、通過帯域コーナー周波数において = -1dB、かつとなるようにスケーリングする必要があります。これが完了すると、設計された非対称チェビシェフフィルタの最終的な伝達関数は以下のようになります。

s=j および s=2j で評価すると、どちらの場合も -1dB の値が得られ、例が正しく合成されたことが保証されます。周波数応答は以下のとおりです。 に対してチェビシェフの 1dB 等リップル通過帯域応答、通過帯域端でのカットオフ減衰量が -1dB、 に向かって -60dB / デケード減衰、 に向かって-20dB / デケード減衰、そして通過帯域端付近でチェビシェフ型の急勾配が見られます。

非対称チェビシェフを示すシミュレーション

通過帯域リップルの抑制

標準的なローパスチェビシェフフィルタ設計では、0 rad/secから周波数正規化値1 rad/secまでの等リップル通過帯域が形成されます。しかし、設計要件によっては、低周波数域で等リップル通過帯域を必要としない場合もあります。この用途に標準的な完全等リップルチェビシェフフィルタを適用すると、過剰設計のフィルタとなってしまいます。等リップルを通過帯域の一定割合に制限することで、より効率的な設計が可能になり、フィルタのサイズが縮小され、部品点数を1~2点削減できる可能性があります。これは、基板スペース効率を最大化し、量産品の製造コストを最小限に抑えるのに役立ちます。[9]

狭窄された通過帯域リップルは、この記事で上で説明した手法を使用し、0 次非対称ハイパス側(0 に透過零点がない)と狭窄されたリップル周波数に設定された を持つ非対称チェビシェフ バンド パス フィルタを設計することによって実現できます。ローパス側の次数は、奇数次フィルタの場合は N-1、偶数次修正フィルタの場合は N-2、標準偶数次フィルタの場合は N です。この結果、 における S12 は 1 未満になりますが、これは偶数次標準チェビシェフ設計の典型であるため、標準偶数次チェビシェフ設計の場合、このステップでプロセスは完了です。奇数次設計の場合は に 1 つの反射零点を挿入し、偶数次修正設計の場合は に 2 つの反射零点を挿入する必要があります。追加された反射零点により、通過帯域に顕著な誤差が生じ、問題になる可能性があります。この誤差は、連立方程式のニュートン法を使用して有限反射零点の位置を変更することで、迅速かつ正確に除去できます

ニュートン法の応用

ニュートン法で反射ゼロを配置するには、次の 3 つの情報が必要です。

  1. 制限されたリップル周波数よりも高い周波数に存在する各通過帯域リップル最小値の位置。
  2. 正規化された振幅の値、すなわち、狭窄周波数および狭窄周波数を超える各極小値における振幅の値。今後この関数を参照する場合は、またはと表記する。
  3. 狭窄周波数および狭窄周波数を超える各最小値における偏導関数ヤコビ行列。各反射ゼロに関して。

チェビシェフ特性方程式 では、すべての反射零点が軸上にあり、すべての透過零点が軸上または軸対称(受動素子の実装に必要)にあるため、通過帯域リップルの最小値の位置は、 の導関数の分子を因数分解し、根を求めるアルゴリズム を用いることで求めることができますこの多項式の根は、通過帯域の最小値周波数となります。 は標準的な多項式導関数の定義から得られ、 は です

偏微分は を用いてデジタル的に計算できますが、連続偏微分は一般的に精度が高く収束時間が短いため、推奨されます。反射零点に関する の連続偏微分を得るには、常にを作用させる連続的な式を得る必要があります。これは、以下に示すように、 をその共役根の対の関数として表すことで実現できます。

ここで、は有限の反射零点と透過零点を含み、 と は反射零点と透過零点の共役対の数を表し、 と は反射零点と透過零点の共役対です。奇数項は、奇数次チェビシェフフィルタに発生する、0 における単一の反射零点を表します。4つ組の透過零点を使用する場合は、4つ組の項に対応するように式を修正する必要があることに注意してください。上記の式では 、必ず であることが分かります。

チェビシェフ通過帯域を形成するには反射零点の移動のみが必要なので、偏微分表現は項についてのみ行えばよく項は定数として扱われます。各 に対する偏微分表現の決定を容易にするために、上記の式は以下のように書き直すことができます。

ここで、特定の反射ゼロ共役対を指定します。

この式の に関する微分は、標準的な微分規則に従って簡単に計算できます。定数 は、関数の完全性を保つために項を除算する必要があります。これを行う最も簡単な方法は、前に移動した項の逆数を乗じることです。微分可能な式は次のように書き直すことができます。

偏微分は、標準的な微分手順を適用して簡約化することで決定できます。結果は以下の通りです。

関連する周波数は狭窄点の周波数と の根だけなので、ヤコビ行列は次のように構築できます。

ここで、 は狭窄限界周波数、 は残りの通過帯域最小値の根の大きさ、 は反射零点です。

フィルタのカットオフ減衰がリップル振幅と同じであると仮定すると、 の値はにおいて常に1なので、解ベクトルの要素はすべて1となり、ニュートン法を解くための反復方程式は次のようになる。

収束は、すべてのフィルタとフィルタの合計がアプリケーションに対して十分に小さい値、通常は1.e-05~1.e-16の範囲にあるときに達成されます。より大きなフィルタの場合、収束初期における過度な変動を防ぐために各フィルタのサイズを制限する必要がある場合があります。また、収束中に各フィルタの値が制限されたリップル範囲内に収まるように、各フィルタのサイズを制限する必要がある場合もあります。

狭帯域通過帯域の例

通過帯域の 55% に制限された 1 dB 等リップル通過帯域を持つ 7 極チェビシェフ フィルタを設計します。

ステップ 1:上記の非対称合成プロセス (コーナー周波数 = 0.45 を使用) を使用して、 に 6 つのローパス極、 0 つのハイパス極を持つ、 .45 から 1 までの非対称周波数応答の特性多項式を設計します

狭窄波紋ステップ1
ステップ1:
7極55%狭帯域リップル通過帯域
1dB等リップル通過帯域
線形周波数スケール


ステップ2:ステップ1のフィルタに反射ゼロを1つ挿入します。(偶数次修正フィルタの場合は、反射ゼロを2つ追加する必要があります)

狭窄波紋ステップ2
ステップ2:
7極55%狭帯域リップル通過帯域
1dB等リップル通過帯域
線形周波数スケール

ステップ3:通過帯域から の根の正の実数値または虚数値を計算してゼロ導関数周波数を決定し、 の最小根を 0.45 の狭窄周波数に置き換えます

計算された反復
10.450.646707850.89924235
20.450.680100030.9147864
30.450.67105970.91089712
40.450.669699720.91042253
50.450.669677630.9104163
60.450.669677620.9104163


ステップ 4 :各制限された微分ゼロ点におけるの値を決定します。

計算された反復
10.450.64035786  0.89703503
21.38865451.16380331.0148793
31.0451081.01337210.99991225
41.00072891.00010940.99998768
51.000000211
6111

ステップ 5 :各周波数でのターゲット値を減算して、線形方程式のBベクトルを作成します。この場合、カットオフ減衰量は、この特定の例では通過帯域リップル減衰量に等しいため、ターゲット値はすべて 1 になります。 カットオフ周波数は です

計算された線形方程式のベクトル反復
1-0.55-0.35964214-0.10296497
20.388654450.16380330.014879269
30.0451080430.013372137-8.7751135e-05
47.2893112e-041.0943442e-04-1.2324941e-05
51.7276985e-075.2176787e-09-2.6640391e-09
61.8873791e-141.5765167e-14-2.553513e-15

ステップ6:各反射零点に関して、各 の偏微分のヤコビ行列を決定する

反復処理
反復1
9.13452413.500252317.567498
-0.35964214-3.121026425.682621
-0.10296497-0.422311545.32731
反復2
18.97830811.68478467.247144
-5.569348515.01497457.94421
-0.46259286-4.758309563.000455
反復3
15.7242518.575108348.268068
-4.830957312.86004248.094251
-0.45645647-4.345539159.167024
反復4
15.3426668.187135546.007638
-4.751692112.65538547.240959
-0.45514037-4.304696358.87818
反復5
15.3370798.180871645.971655
-4.750678912.65328347.233095
-0.45510227-4.304239158.875318
反復6
15.3370788.180870245.971647
-4.750678712.65328347.233094
-0.45510225-4.304239158.875317

ステップ 7 :ステップ 5 のBベクトルを使用して線形方程式を解き、反射ゼロの移動を取得します。

1-0.033937389-0.040973291-0.0054977233.02680
20.0101591030.0104363530.001099011.00723149
30.00181702710.0013144721.090765e-04.00108019
43.4653892E-051.6843291E-051.2899974E-061.75957e-05
59.0033707E-092.9081531E-092.3695501E-104.04949e-08
60000

ステップ 8:上記で計算した値を過去の反射ゼロ位置の反復から差し引いて、新しい反射ゼロ位置を計算します。

10.539825090.816376410.97841993
20.529665990.805940060.97732092
30.527848960.804625590.97721185
40.527814310.804608740.97721056
50.52781430.804608740.97721056
60.52781430.804608740.97721056

適用収束基準を満たすまで、手順3から8を繰り返します。この例では、収束基準は1.e-12としています。収束基準を満たすと、最終的な反射零点の位置、+/-j0.5278143、+/-J0.80460874、+/-J0.97721056、および0から最終的な が構築されます。振幅を となるように正規化すると、構築された は以下のようになります。

合成プロセスは、ステップ3の各周波数における簡単なチェックを行うことで検証できます。これにより、これらの周波数で1dBの減衰が確保され、カットオフ減衰も1dBであることが保証されます。以下の計算結果の要約は、この例の合成プロセスを検証するものです。

検証の概要
-1 dB
-1 dB
-1 dB
-1 dB

順方向伝達関数の最終的な振幅周波数応答を以下に示します。

狭窄したリップル最終
最終ステップ:
7極55%狭帯域リップル通過帯域
1dB等リップル通過帯域
線形周波数スケール

チェビシェフIIストップバンドリップル狭窄

標準的なローパス逆チェビシェフフィルタ設計では、正規化された値1 rad/secから始まる等リップル阻止帯域が生成される。しかし、設計要件によっては、高周波域で等リップル通過帯域を必要としない場合もある。この用途に標準的な完全等リップル逆チェビシェフフィルタを適用すると、過剰設計のフィルタとなってしまう。等リップルを阻止帯域の一定割合に制限することで、より効率的な設計が可能になり、フィルタのサイズが縮小され、部品点数を1~2個削減できる可能性がある。これは、基板スペース効率を最大化し、量産品の生産コストを最小化するために役立つ。[9]

狭窄リップルを持つ逆チェビシェフフィルタは、標準的な逆チェビシェフと全く同じプロセスで合成されます。狭窄リップルチェビシェフは、反転した (ここではdB単位のストップバンド減衰量)を用いて設計され、設計された狭窄リップルチェビシェフフィルタの極と零点は反転され、カットオフ減衰量は に設定されます。標準的なチェビシェフ方程式は狭窄リップル設計には適用できないため、カットオフ減衰量は楕円砂時計設計で説明したプロセスを使用して設定する必要があります。

以下は、3dB カットオフ減衰を備えた 7 極狭窄リップル逆チェビシェフ フィルタの|S11| および |S12|散乱パラメータです。

逆チェビシェフ狭窄リップル
7極逆チェビシェフ狭窄阻止帯域リップル

非標準カットオフ減衰と透過ゼロ

上記の狭窄リップルの例は、カットオフ減衰量を通過帯域リップル減衰量と等しくし、任意の透過零点を省略し、偶数次への変更を必要としない奇数次を用いることで、意図的に単純化されています。ただし、ステップ5で目標値をカットオフ周波数 に存在する必要な 1 からオフセットするように計算し透過零点を含む微分定数の一部として分母を含め、ステップ2で元の に反射零点を1つではなく2つ挿入することで、非標準的なカットオフ減衰量にも対応できます

阻止帯域透過零点を含める場合、 の根には の阻止帯域極大値が含まれることを覚えておくことが重要です。これらの根は、計算で使用される通過帯域極小値には含めないでください。

は のカットオフ減衰を設定するために使用できるため、ステップ 5 の目標値は 1 に関して作成できます。ステップ 5 の目標値は、上記の式から取得できる の式を使用して計算できます。

%constriction = 55、次数 = 8、単一透過零点 = 1.1、通過帯域リップル減衰量 = 0.043648054(ロスレスネットワークの関係式に基づくとS12 = 20dB減衰に相当[12])、通過帯域カットオフ減衰量 = 20dBのフィルタ設計を考えてみましょう。

ステップ 5 の目標値は .01010101 で、計算する値は 99 です。完了すると、特性多項式、、および順方向伝達関数、は以下のようになります。

検証は、以下に示すように、狭窄周波数、カットオフ周波数、その間の残りの通過帯域最小周波数、および透過ゼロ周波数の散乱パラメータ (それぞれおよび )を計算することから構成されます。

8極非標準カットオフ減衰および伝送ゼロ検証の概要
-0.043648054 dB-20dB
-0.043648054 dB-20dB
-0.043648054 dB-20dB
-20 dB-0.043648054 dB
-0 dB


最終的な振幅周波数応答を以下に示します。

有限の透過零点と非標準のカットオフ減衰を持つ8極の狭窄リップルチェビシェフ
最終ステップ:
8極55%狭帯域リップル通過帯域
20dB S11等リップル通過帯域
1.1 rad/sec での有限透過率ゼロ
非標準S12カットオフ減衰20dB
幾何学的頻度スケール

参照

参考文献

  1. ^ ダニエルズ、リチャード・W. (1974). 『電子フィルタ設計のための近似法』 ニューヨーク: マグロウヒル. ISBN 0-07-015308-6
  2. ^ ルトヴァック、ミロスラフ D.;ルトヴァック、D.トシッチ、デヤン V.エヴァンス、ブライアン・ローレンス (2001)。 MATLAB と Mathematica を使用した信号処理用のフィルター設計。プレンティス・ホール。ISBN 9780201361308
  3. ^ Weinberg, Louis; Slepian, Paul (1960年6月). 「TakahasiのTchebycheffおよびButterworth Ladder Networksに関する結果」. IRE Transactions on Circuit Theory . 7 (2): 88– 101. doi :10.1109/TCT.1960.1086643.
  4. ^ ウィリアムズ、アーサー・B.、テイラーズ、フレッド・J. (1988). 『電子フィルタ設計ハンドブック』 ニューヨーク:マグロウヒル. ISBN 0-07-070434-1
  5. ^ abcd ザール、ルドルフ (1979 年 1 月)。フィルター設計ハンドブック (英語およびドイツ語) (第 1 版)。ミュンヘン、ドイツ: Allgemeine Elektricitais-Gesellschaft。pp . 25、26、56–61、116、117。ISBN 3-87087-070-2
  6. ^ Paarmann, Larry D. (2001). アナログフィルタの設計と解析:信号処理の観点から. 米国マサチューセッツ州ノーウェル:Kluwer Academic Publishers. pp. 137, 138. ISBN 0-7923-7373-1
  7. ^ Paarmann, Larry D. (2001). アナログフィルタの設計と解析:信号処理の観点から. 米国マサチューセッツ州ノーウェル:Kluwer Academic Publishers. pp. 161, 162. ISBN 0-7923-7373-1
  8. ^ Matthaei, George L.; Young, Leo; Jones, EMT (1980). 『マイクロ波フィルタ、インピーダンス整合回路、および結合構造』 ノーウッド、マサチューセッツ州: Artech House. ISBN 0-89-006099-1
  9. ^ abc Pelz, Dieter (2005). 「狭帯域等リップル通過帯域を備えたマイクロ波ローパスフィルタ」(PDF) . AMW . 13 (7): 28~34 – APPLIED MICROWAVE & WIRELESS経由.
  10. ^ abcde バイロン・ベネット博士のフィルタ設計講義ノート、1985年、モンタナ州立大学電気工学部、 ボーズマン、モンタナ州、米国
  11. ^ ab Sedra, Adel S.; Brackett, Peter O. (1978). 『フィルタ理論と設計:能動フィルタと受動フィルタ』Beaverton, Oegon, US: Matrix Publishers, Inc. pp.  45– 73. ISBN 978-0916460143
  12. ^ 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{{cite book}}: CS1 maint: location (link)
  • ウィキメディア・コモンズのチェビシェフフィルタ関連メディア
Retrieved from "https://en.wikipedia.org/w/index.php?title=Chebyshev_filter&oldid=1321673877"