[論文レビュー] Stochastic solutions of the wave equation

Sourav Chatterjee|arXiv (Cornell University)|Jun 10, 2013

Stochastic processes and financial applications参考文献 3被引用数 3

ひとこと要約

本稿では、ブラウン運動、標準コーシー分布に従う確率変数、および最初の到達時刻を用いて、任意の次元および有界領域における波動方程式の非自明な解のクラスを確率的表現として提示する。主な貢献は、このような解が波動方程式 $\partial_t^2 u = \Delta u$ を満たすことを証明することであるが、この手法はすべての解をカバーするものではなく、初期データとの直接的な関連性も欠けている。

ABSTRACT

Unlike the heat equation or the Laplace equation, solutions of the wave equation on general domains have no known stochastic representation. This short note gives a simple solution to this well known problem in arbitrary dimensions. The proposed representation has several shortcomings, one of which is that it does not cover all solutions. Still, it is proof that a large class of nontrivial solutions of the wave equation in general dimensions and domains may indeed be represented stochastically.

研究の動機と目的

• 一般領域における波動方程式の解に対する確率的表現を求める長年の未解決問題に取り組む。
• ラプラス方程式および熱方程式に対して知られている確率的表現を、双曲型の波動方程式へと拡張する。
• ブラウン運動と確率変数を用いた確率的手法により、広範な非自明な解のクラスを生成できることを示す。
• 境界データを関数 $f(t,x)$ に符号化した構成的技法を提供する。

提案手法

• 解 $u(t,x)$ を期待値 $\mathbb{E}_x[f(tX + \sqrt{\tau}Z, B_\tau)]$ として定義する。ここで $X$ は標準コーシー分布、$Z$ は標準正規分布に従い、$\tau$ はブラウン運動 $B_t$ が境界 $\partial D$ に最初に到達する時刻である。
• 拡張された領域 $D' = \mathbb{R} \times D$ における強マルコフ性および調和関数論を用いて、条件付き期待値 $v(t,x) = \mathbb{E}_x[f(t + \sqrt{\tau}Z, B_\tau)]$ が $\mathbb{R} \times D$ 上で調和的であることを示し、したがって $\partial_t^2 v = -\Delta v$ を満たすことを示す。
• 導関数の一様有界性を保証するため、$a > 0$ に対して $v_a(t,x) = e^{-a t^2} v(t,x)$ を導入し、優越収束定理および積分の微分可能性を保証する。
• $a \to 0$ の極限を取ることで、元の解 $u(t,x) = \mathbb{E}[v(tX, x)]$ を回復し、二階微分の収束を用いて $\partial_t^2 u = \Delta u$ が成り立つことを示す。
• 任意の停止定理およびモーメント母関数の恒等式（例：$\mathbb{E}[e^{i t X}] = e^{-|t|}$）を用いて、具体的な例（例：$D = (-1,1)$ 上で $u(t,x) = e^{x - t}$）を計算する。
• 導関数と極限の交換を、$v$ および $v_a$ の導関数の一様有界性により正当化し、$u$ が $C^2$ 級であることを保証する。

実験結果

リサーチクエスチョン

• RQ1 任意の次元および一般の有界領域において、波動方程式の非自明な解に対して確率的表現を構成可能か？
• RQ2 提案手法は、ブラウン運動と独立な確率変数のみを用いて、$\partial_t^2 u = \Delta u$ を満たす解を生成可能か？
• RQ3 解は初期条件 $u(0,x)$ および $\partial_t u(0,x)$ とどのように関係するか？ また、それらに適合させるために調整可能か？
• RQ4 すべての有界解がこの形で表現可能か、それとも生成される解のクラスに根本的な制限があるか？
• RQ5 コーシー分布および正規分布に従う確率変数が、期待値を通じて波動方程式の構造を実現するために果たす役割は何か？

主な発見

• 解 $u(t,x) = \mathbb{E}_x[f(tX + \sqrt{\tau}Z, B_\tau)]$ は $t > 0$ および $x \in D$ において $t$ と $x$ に関して $C^2$ 級であり、$\mathbb{R}_+ \times D$ で $\partial_t^2 u = \Delta u$ を満たす。
• 関数 $f: \mathbb{R} \times \partial D \to \mathbb{R}$ が任意の有界可測関数であるため、境界上の多様な時間的変化を表現可能であり、広範な非自明な解のクラスが生成可能である。
• 明示的な例を計算した：$D = (-1,1)$、$f(t,x) = e^x \cos t$ のとき、解は $u(t,x) = e^{x - t}$ であり、$\mathbb{R}_+ \times (-1,1)$ で波動方程式を満たす。
• この表現はすべての解をカバーするものではない。例えば $e^{x + t}$ は解ではあるが、無限大に発散するため、与えられた規定のもとでは表現不可能である。
• 初期データ $u(0,x)$ および $\partial_t u(0,x)$ との直接的な関連性が欠如しており、$f$ の時間依存性が一般に逆写像不能であるため、そのような関係を明示的に取り入れることは困難である。
• 証明は $v_a(t,x) = e^{-a t^2} v(t,x)$ を用いた正則化、導関数の一様有界性、および優越収束定理に基づく $a \to 0$ の極限の正当化に依存しており、極限においても波動方程式が保たれることを保証する。