[論文レビュー] Some asymptotics for the Bessel functions with an explicit error term

Ilia Krasikov|arXiv (Cornell University)|Jul 11, 2011

Mathematical functions and polynomials参考文献 16被引用数 5

ひとこと要約

本稿では、Sonin関数を用いた事前境界を介して、誤差項を明示的に評価可能な方法により、ベッセル関数およびエアリー関数の漸近展開を導出する新規な手法を開発する。標準的なベッセル関数の漸近展開における誤差の上界が、明示的な定数を用いて $\nu^2$ で抑えられることを証明し、オリェンコの懸案していた問題を解決する。また、両関数の零点に対する非常に高精度な近似式を、鋭い誤差制御のもとで提供する。

ABSTRACT

We show how one can obtain an asymptotic expression for some special functions satisfying a second order differential equation with a very explicit error term starting from appropriate upper bounds. We will work out the details for the Bessel function $J_ν(x)$ and the Airy function $Ai(x)$ and find a sharp approximation for their zeros. We also answer the question raised by Olenko by showing that $$c_1 | ν^2-1/4\,| < \sup_{x \ge 0} x^{3/2}|J_ν(x)-\sqrt{\frac{2}{πx}} \, \cos (x-\frac{πν}{2}-\fracπ{4}\,)|

研究の動機と目的

• 特殊関数の漸近展開を、明示的かつ一様な誤差項を伴って導出する体系的な手法の開発。
• 標準的なベッセル関数の漸近展開における誤差の $\nu^2$-バウンドの鋭さに関して、オリェンコが提起した未解決問題の解決。
• エアリー関数 $Ai(-x)$ の正の零点および、それらを拡張してベッセル関数 $J_\nu(x)$ の正の零点に対する、高精度で明示的な近似式の提供。
• ベッセル関数およびエアリー関数の振動領域、遷移領域、単調領域にわたり、誤差を厳密に制御した手法の拡張。

提案手法

• Sonin関数を用いて、ベッセル関数およびエアリー関数の絶対値に対する事前上界を導出し、誤差推定を可能にする。
• 2階常微分方程式を、小さな誤差項 $\epsilon(x)$ を持つ摂動された調和振動子方程式に変換するWKB型変換を適用する。
• 変数係数法を用いて誤差の積分表現を導出し、$\int \epsilon(t) b(t) dt$ に比例する上界を得る。
• 解の上界を用いて、元の関数のノルムに関する漸近近似の誤差を推定する。
• この手法を、振動領域および遷移領域におけるベッセル関数 $J_\nu(x)$ およびエアリー関数 $Ai(-x)$ に適用し、明示的な誤差項を伴うより精緻な漸近展開を改善する。
• 位相補正と振幅補正を一致させることで、$Ai(-x)$ の零点に対するより鋭い近似式を導出し、明示的な誤差項を伴う非常に高精度な零点公式を導く。

実験結果

リサーチクエスチョン

• RQ1 特殊関数の漸近展開を、明示的かつ一様な誤差項を伴って導出する一般的手法を開発できるか。
• RQ2 標準的なベッセル関数の漸近展開 $J_\nu(x) \sim \sqrt{\frac{2}{\pi x}} \cos(x - \frac{\pi\nu}{2} - \frac{\pi}{4})$ における誤差が、明示的な定数 $c_1, c_2$ を用いて $c_1 \nu^2$ と $c_2 \nu^2$ の間に抑えられるか。
• RQ3 この手法により、エアリー関数 $Ai(-x)$ の正の零点に対するより鋭い近似式が得られるか、かつ明示的な誤差項を伴うか。
• RQ4 振動領域、遷移領域、単調領域にわたり、漸近展開の誤差項を一様に制御できるか。

主な発見

• 本稿では、明示的な数値定数 $c_1$ および $c_2$ を用いて、$c_1 \nu^2 < \sup_{x \geq 0} x^{3/2} \left| J_\nu(x) - \sqrt{\frac{2}{\pi x}} \cos\left(x - \frac{\pi\nu}{2} - \frac{\pi}{4}\right) \right| < c_2 \nu^2$ を証明し、オリェンコの問いに答えを提示する。
• 新たな漸近近似式として、$J_\nu(x)$ について $J_\nu(x) = \sqrt{\frac{2}{\pi}} \frac{\cos\left(x - \frac{\mu}{2x} - \omega_\nu\right)}{(x^2 + \mu)^{1/4}} + \theta \frac{25\mu}{24\sqrt{2\pi} x^3 (x^2 + \mu)^{1/4}}$ を導出し、$|\nu| \leq \frac{1}{2}$ に対して有効であり、誤差項が明示的に評価可能である。
• エアリー関数 $Ai(-x)$ の $s$ 番目の正の零点に対する明示的近似式として、$a_s = 16^{-2/3} \left(m + \sqrt{m^2 + 40}\right)^{2/3} + \theta \frac{1280\pi}{9m^3 (m^2 + 40)^{1/6}}$ を提示し、$m = (12s - 3)\pi$ である。誤差項は $O(m^{-3})$ のオーダーである。
• 本手法により、エアリー関数に対するより鋭い近似式が得られる：$Ai(-x) = \frac{2\sqrt{x} \cos\left(\frac{2}{3}x^{3/2} - \frac{5}{48}x^{-3/2} - \frac{\pi}{4}\right)}{\sqrt{\pi} (16x^3 + 5)^{1/4}} + \theta \frac{10\sqrt{3}}{\sqrt{\pi} x^{1/4} (16x^3 + 5)^{3/2}}$ であり、小さな明示的誤差項を伴う。
• エアリー関数の零点近似式は非常に高い精度を達成する：最初の零点において、誤差は 0.00122 未満である。これにより、本手法の有効性が示される。