[論文レビュー] On the numerical solution of second order differential equations in the high-frequency regime

James Bremer|arXiv (Cornell University)|Sep 21, 2014

Numerical methods for differential equations参考文献 21被引用数 4

ひとこと要約

本稿では、非周期的位相関数を活用することで、高周波数の2階線形常微分方程式を周波数に依存しない高速なアルゴリズムで解く手法を提示する。Kummerの微分方程式の対数形を用いて問題を再定式化し、高価な周期的挙動の解消を回避することで、解の位相関数を計算し、λに依存しない実行時間で解を得られる。Bessel関数、ルジャンドル関数、および回転楕円体波関数に対して、高い精度と効率で検証された。

ABSTRACT

We describe an algorithm for the numerical solution of second order linear differential equations in the highly-oscillatory regime. It is founded on the recent observation that the solutions of equations of this type can be accurately represented using nonoscillatory phase functions. Unlike standard solvers for ordinary differential equations, the running time of our algorithm is independent of the frequency of oscillation of the solutions. We illustrate the performance of the method with several numerical experiments.

研究の動機と目的

• 高周波数領域における標準的ODEソルバーの非効率性、特に解の周期的挙動に起因するλに比例する計算コストの増大を是正すること。
• 特殊関数を定義するような、広範な2階ODEに適用可能な汎用アルゴリズムの開発。
• Kummerの微分方程式の対数形から導出される非周期的位相関数を用いることで、λに依存しない解の精度を達成すること。
• Bessel関数、ルジャンドル関数、および回転楕円体波関数などの特殊関数を、高精度かつ最小限の計算負荷で高速に評価可能にすること。

提案手法

• ODE y'' + λ²q(t)y = 0 の解が u(t) = cos(α(t))/|α'(t)|¹ᐟ² および v(t) = sin(α(t))/|α'(t)|¹ᐟ² で表される非周期的位相関数 α(t) を構築する。
• r'(t) = λ exp(r(t)/2) を用いてKummerの3階非線形ODEをより取り扱いやすい対数形に変換し、r'' - ¼(r')² + 4λ²(exp(r) - q) = 0 を得る。
• この対数形を、特異点付近ではグレーデッドメッシュを用い、部分区間では高次数のチェビシェフスペクトルコロケーション法で数値的に解く。
• 各部分区間で高次数のチェビシェフ展開を用いて解を表現することで、任意の点での効率的評価を可能にする。
• 位相関数を用いることで、直接的に高速な周期的挙動を解消するのを避け、計算コストがλに依存しないことを保証する。
• Bessel関数、ルジャンドル関数、および回転楕円体波関数の式に対して検証を行い、Osipov-Rokhlin法など既存のソルバーと比較した。

実験結果

リサーチクエスチョン

• RQ1 周波数λに依存しない、高周波数2階ODEに一般に適用可能な数値的手法を開発できるか？
• RQ2 特殊関数を定義するような広範なODEクラスに対し、非周期的位相関数が存在し、効率的に計算可能か？
• RQ3 Kummerの微分方程式の対数形を高精度に数値的に解くことで、解の高速評価が可能か？
• RQ4 回転楕円体波関数に対して、Osipov-Rokhlin法と比較して本手法の性能と精度はどの程度か？
• RQ5 高次のルジャンドル関数に対して、本手法が機械精度に近い精度を達成できるか？

主な発見

• 本アルゴリズムは、μがqに依存するがλに依存しない (μλ)⁻¹exp(−μλ) のオーダーの解の精度を達成し、λが大きくても高精度を維持する。
• 本アルゴリズムの実行時間はλに依存せず、標準的ソルバーとは異なりλに比例して増加しない。
• ν = π×10⁵ のルジャンドル関数に対して、グレーデッドメッシュ上での16次チェビシェフ展開を用い、最大絶対誤差が10⁻¹⁴未満で解を計算した。
• 回転楕円体波関数の評価において、Osipov-Rokhlin法と同等の精度を維持しながら、評価速度で本手法が優れた性能を示した。
• 初等関数で閉形式で表せないBessel方程式の非周期的位相関数も、正確に計算可能であり、本手法の一般性を示している。
• 数値実験により、端点に特異係数を有する多様なODEに対しても、本手法が高精度かつ高効率を維持することが確認された。