[論文レビュー] An elementary proof of asymptotic behavior of solutions of $u"=Vu$

Giorgio Metafune, Motohiro Sobajima|arXiv (Cornell University)|May 22, 2014

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

ひとこと要約

本稿では、ルベーグ積分や優越収束定理といった高度な解析的手法に依存せずに、グロンウォールの不等式を用いて2階常微分方程式 $ u'' = V(x)u $ の解の漸近的挙動を初等的に証明する。主な貢献は、ポテンシャル $ V(x) $ に対して最小限の可積分性仮定のもとで、特にベッセル関数や修正ベッセル関数のようなケースにおいて、解の長期的挙動を簡潔かつ自己完結的に導出することにある。

ABSTRACT

We provide an elementary proof of the asymptotic behavior of solutions of second order differential equations.

研究の動機と目的

• 高階の積分論に依存せずに、$ u'' = V(x)u $ の解の漸近的挙動を自己完結的かつ初等的に証明すること。
• 逐次近似の代わりにグロンウォールの不等式を用いることで、$ f = \pm 1 $ の基本的ケースにおいて、古典的な結果を簡略化すること。
• 無限遠点および孤立特異点近傍における挙動の両方の分析を拡張し、特異ポテンシャルを有する方程式に対しても適用可能にすること。
• ベッセル関数や修正ベッセル関数といった特殊関数に本結果を適用し、統一的な枠組みで既知の漸近展開を回復すること。
• 一般形 $ u'' + g(x)u' = V(x)u $ が変数変換によって標準形に還元され、その方法の適用性が保たれることを示すこと。

提案手法

• 変数変換法を用い、$ f \equiv 1 $ および $ f \equiv -1 $ の場合に、指数関数的または振動的解の摂動として解を表現すること。
• グロンウォールの補題を用いて摂動問題における解の成長を制御し、反復近似法の代わりにすること。
• $ u = e^{\int g} v $ の変換により、一般の2階常微分方程式を $ u'' = V(x)u $ の形に還元すること。
• リウヴィル変換を用い、$ \psi_{f,g} \in L^1(0,\infty) $ の条件下で、一般形 $ u'' = (f + g)u $ をモデルケース $ f \equiv \pm 1 $ に還元すること。
• 変数 $ x \mapsto x^{-1} $ の逆変換を用いて $ x = 0 $ 近傍の挙動を無限遠点での挙動に変換し、解析すること。
• 関数 $ z(x) = e^{-\zeta x}u(x) $ の推定を通じて漸近的同値性を確立し、$ x \to \infty $ のとき $ z(x) $ および $ z'(x) $ が定数に収束することを示すこと。

実験結果

リサーチクエスチョン

• RQ1 ルベーグ積分や優越収束定理を避けて、$ u'' = V(x)u $ の解の漸近的挙動を最小限の解析的道具でどのように導出できるか。
• RQ2 $ f(x) = \pm 1 $ かつ $ g \in L^1(0,\infty) $ のとき、$ V(x) = f(x) + g(x) $ に対する解の正確な漸近的挙動は何か。また、特殊関数とどのように関係するか。
• RQ3 孤立特異点（例：$ x = 0 $）近傍の解を解析するために、変数変換をどのように適応できるか。
• RQ4 この枠組みから得られるベッセル関数および修正ベッセル関数の漸近展開は何か。また、既知の結果と一致するか。
• RQ5 1階項を含む一般の2階常微分方程式が、漸近解析を保ったまま標準形 $ u'' = V(x)u $ に還元可能かどうか。

主な発見

• $ u'' = (1 + g(x))u $ で $ g \in L^1(0,\infty) $ のとき、$ e^{-x}u_1(x) \to 1 $、$ e^x u_2(x) \to 1 $（$ x \to \infty $）を満たす解 $ u_1 $, $ u_2 $ が存在し、対応する微分係数はそれぞれ 1 および -1 に収束する。
• $ u'' = (-1 + g(x))u $ の場合、解は $ e^{-ix}u_1(x) \to 1 $、$ e^{ix}u_2(x) \to 1 $ を満たし、微分係数はそれぞれ $ i $ および $ -i $ に収束する。
• $ x = 0 $ 近傍で、$ V(x) = f(x) + g(x) $ で $ f(x) \sim x^{-2} $ のとき、解は $ r^{-\nu}I_\nu(r) \to c_3 \neq 0 $、$ r^{\nu}K_\nu(r) \to c_4 \neq 0 $（$ r \downarrow 0 $）として振る舞い、$ \nu \neq 0 $ のとき成り立つ。
• $ \nu = 0 $ のとき、修正ベッセル関数 $ I_0(r) $ は $ r \to 0 $ で有界であり、$ |\log r|^{-1}K_0(r) \to c_6 \neq 0 $ を満たす。
• ベッセル関数に関しては、$ r \downarrow 0 $ のとき $ r^{-\nu}J\_\nu(r) \to c_1 \neq 0 $、$ r^{\nu}Y_\nu(r) \to c_2 \neq 0 $ であり、$ \nu = 0 $ のとき $ \log r \cdot J_0(r) \to c_3 \neq 0 $ が成り立つ。
• $ r \to \infty $ のとき、$ \sqrt{r}J_\nu(r) \sim c_5 \cos(r + \theta_1) $、$ \sqrt{r}Y_\nu(r) \sim c_6 \cos(r + \theta_2) $ であり、$ c_5, c_6 \neq 0 $ かつ $ \theta_1 \neq \theta_2 $ である。