技術要約：$\sqrt{f}$ 特性を有する分数次フィルタ

問題提起

本論文は、伝達関数が $s^{\alpha}$（$\alpha = \pm 1/2$ に対応し、それぞれ $\sqrt{f}$ および $1/\sqrt{f}$ の周波数特性に対応）に比例する分数次フィルタの設計と解析について扱っている。これらのフィルタは、プリエンファシス／デエンファシス、$1/f$（フリッカー）ノイズのパワースペクトルのモデリング、分数次制御、および生体組織のような媒体の分散特性のシミュレーションなどのアプリケーションにおいて不可欠である。

Oustaloup 再帰的フィルタや Charef 法といった既存の近似手法は広く用いられているが、これらは位相または振幅レスポンスを平坦化するために、対数スケール上に極と零点を配置することに依存している。本論文は、純粋に数学的な近似に基づくのではなく、物理法則に根ざした代替的な導出方法を追求しており、マクスウェル方程式の離散化を通じて、所望の分数次特性を自然に発現させる挙動モデルの導出を目指している。

手法

1. マクスウェル方程式からの物理的導出

核心となる手法は、線形、等方性、かつ均質な導電性ロッドに対するマクスウェル方程式から始まる。変位電流が無視できる（$J \gg \partial D/\partial t$）という仮定の下、長さを $l$、幅を $b$、厚さを $d$ とし、$l \gg b, d$ である矩形ロッドに対する準定常的なアプローチを用いることで、これらの方程式は電場（$E_z$）および磁場（$H_x$）を記述する偏微分方程式（PDE）の系へと簡略化される。

本論文は、強スキンテフェクト領域（スキン深さが幾何学的形状に対して十分に小さい領域）において、導電性材料のスキン深さ（$\Delta x$）が $1/\sqrt{f}$ の法則に従うという既知の物理現象を利用している。目的は、これらの PDE を離散化することで、この挙動を模倣するフィルタを作成することである。

2. 空間離散化と回路モデリング

PDE は有限差分（FD）スキームを用いて離散化される。

• 一様格子： 最初は等間隔の格子が検討され、これにより常微分代数方程式（DAE）の系が導かれる。この系は、インダクタと抵抗（LR）、またはキャパシタとコンダクタ（GC）からなる受動的挙動回路モデル（ラダーネットワーク）へとマッピングされる。
• 非一様格子： 一様格子では、高周波におけるスキンエフェクトを捉えるために非現実的な数の点が必要となることを認識し、本論文では非一様格子を提案している。格子の間隔は、スキンエフェクトが最も顕著な表面（$y = d/2$）で最小となる。隣接する間隔の比はパラメータ $c$ によって定義される。
• 回路実現： 離散化された方程式は、等価な受動回路として実現される。本論文では、LR モデルが自然ではあるものの、修正ノード解析（MNA）における未知数の数を削減できるため、GC モデル（キャパシタとコンダクタを使用）が SPICE ライクなシミュレータにはより適していると述べている。

3. デジタル変換

デジタル等価物を得るために、アナログフィルタは双線形変換（台形則）を用いて変換され、後退微分公式（BDF）と比較される。

• 双線形変換： 虚軸を単位円に写像し、BDF メソッドで見られるような数値的な減衰（数値的ダンピング）を導入しないため、こちらが好ましい。また、高次の BDF メソッド（$p \ge 2$）と比較して、より低いフィルタ次数をもたらす。
• 実装： 結果として得られる差分方程式の効率的な実装戦略について、特定の再帰ステップにおいて乗算器を回避するための補助ベクトルの導入を含めて議論している。

主な貢献

1. 分数次フィルタの物理的導出： 本論文は、導電性ロッドのマクスウェル方程式の離散化から $\sqrt{f}$ および $1/\sqrt{f}$ フィルタを直接導出し、極と零点の配置というヒューリスティックな手法ではなく、分数次挙動の物理的根拠を提供している。
2. 非一様格子の最適化： 特定の間隔比（$c$）を用いた非一様格子を使用することで、精度を維持しながらシステム次数を大幅に削減できることを示している。本論文では、リンギング現象を避けるために、隣接する格子間隔の比はオイラー数よりも小さくなければならないと特定しており、$c=2$ から $c=2.5$ が最適であるとしている。
3. 挙動回路モデル： 導出されたモデルは受動的なラダー回路（LR および GC 両方）であり、オーム損失によって極が左半平面内に留まるため、本質的に安定している。
4. デジタル実装： 双線形変換を用いたデジタル実現への明確な経路を提供しており、安定かつ効率的なデジタルフィルタ構造を実現している。

結果

• Oustaloup との比較： フィルタ次数 $K=11$ におけるシミュレーション結果によれば、本論文で導出されたスキン深さモデルは、Oustaloup 再帰的フィルタと比較して、近似の妥当な周波数範囲が広く、位相精度も高い。
• 周波数範囲： $K=8$ 点、 $c=2$ の非一用格子を用いた場合、モデルは期待される $\sqrt{f}$ 挙動（インピーダンスの二重対数プロットにおける 0.5 の線形傾斜）を、約 4 デケードの周波数にわたって示す。理論的な推定ではより広い範囲が示唆されていたが、シミュレーションでは実用的な妥当範囲はおよそ $f_1/10$ であることが示された。
• パラメータ感度： 格子間隔比 $c$ の値が 2 からオイラー数（$\approx 2.718$）の間にあるとき、システム次数、妥当な範囲、および正確性の間で最良のトレードオフが得られることを確認した。$c$ がオイラー数を超える値になると、望ましくないリンギング現象が発生する。
• デジタル性能： 双線形変換を介して実装されたデジタル対応モデルは、正規化角周波数の約 4 デケードにわたって平方根の挙動を維持しており、変換の有効性を裏付けている。

意義と主張

本論文は、提案手法が、特に妥当な範囲と位相精度に関して、Oustaloup のような既存の再帰的近似に対して優れた代替案を提供すると主張している。導電体の物理的なスキンエフェクトに根ざしていることで、得られるモデルは数学的に健全であるだけでなく、受動ネットワークとして物理的に解釈可能である。

著者は、非一様離散化戦略を用いることで、一様離散化よりも大幅に少ない格子点（したがって計算コストの低減）でスキンエフェクトを正確にモデル化できることを強調している。双線形変換の選択は、安定性を維持し、数値的減衰を回避できるという点において正当化されており、ハードウェア実装に適している。

最後に、本論文は、この研究を分数次フィルタリング（$\alpha = \pm 1/2$）の特定事例として位置づけ、偏微分方程式を離散化するという基礎的な手法が、将来的に $\pm 1/2$ を超える任意の次数 $\alpha$ の分数次フィルタにも一般化できる可能性を示唆している。