[論文レビュー] A generalization of Osgood's test and a comparison criterion for integral equations with noise

Marcos J. Ceballos-Lira, Jorge E. Macías‐Díaz+1 more|arXiv (Cornell University)|Dec 8, 2010

Probabilistic and Robust Engineering Design参考文献 13被引用数 5

ひとこと要約

この論文は、確率的ノイズを含む積分方程式における有限時刻爆発のためのOsgoodのテストを一般化し、比較基準を確立し、爆発時間の確率的境界を導出する。この枠組みを材料内のき裂伝播のモデルに適用し、特にブラウン運動の絶対値としてのノイズが、破壊時刻の予測に顕著な影響を及ぼすことを示し、有限時刻における爆発の確率的推定値を明示的に提供する。

ABSTRACT

In this work, we prove a generalization of Osgood's test for the explosion of the solutions of initial-value problems. We also establish a comparison criterion for the solution of integral equations with noise, and provide estimations of the time of explosion of problems arising in the investigation of crack failures where the noise is the absolute value of the Brownian motion

研究の動機と目的

• 常微分方程式におけるOsgoodの古典的テストを、確率的ノイズを含む積分方程式へと拡張すること。
• 爆発挙動の相対的分析を可能にするために、ノイズを含む積分方程式の解の間の比較基準を確立すること。
• ノイズがブラウン運動の絶対値としてモデル化される場合の爆発時刻の確率的推定値を導出すること。
• 剛体材料におけるき裂破壊ダイナミクスなどの物理的問題への理論的枠組みの応用。
• こうしたノイズ付きシステムの爆発時刻の上界および下界を提供し、工学的文脈における予測精度を向上させること。

提案手法

• ノイズを含む一般積分方程式を導入：$ X_t = x_0 + \int_0^t a(s)b(X_s)ds + g(t) $、ここで $ g(t) $ はノイズを表す。
• 変数変換 $ Y_t = X_t - g(t) $ を用いてノイズなし系に変換し、$ Y_t = x_0 + \int_0^t a(s)b(Y_s + g(s))ds $ を得る。
• Osgoodの基準を一般化するための補助関数 $ A(t) = \int_0^t a(s)ds $ および $ B_r(x) = \int_{x_0 - r}^x \frac{ds}{b(s)} $ を定義する。
• 関数 $ b $ と $ a $ の単調性および成長度に応じて、異なる系の解を関連付ける比較原理を適用する。
• $ B $ および $ A $ の逆関数を用いて、条件 $ B(\infty) < A(\infty) $ の下で爆発時刻を $ T_e = A^{-1}(B(\infty)) $ と表現する。
• 標準正規分布の尾確率を用いて、ブラウン運動の絶対値の分布を応用し、爆発時刻の確率的境界を導出する。

実験結果

リサーチクエスチョン

• RQ1 ノイズを含む積分方程式の解が有限時刻に爆発する条件は何か、そしてそれがどのように解析的に特徴付けられるか？
• RQ2 特にブラウン運動の絶対値としてのノイズが、き裂伝播モデルにおける爆発までの時間にどのように影響を与えるか？
• RQ3 異なるノイズ付き積分方程式の解の間で、爆発挙動を予測可能な比較基準を確立できるか？
• RQ4 こうしたシステムにおける爆発時刻の最もきつい上界および下界は何か、そしてそれらがノイズ関数および駆動関数にどのように依存するか？
• RQ5 与えられた時間間隔内に爆発が発生する確率は何か、そしてそれが確率過程を用いてどのように推定できるか？

主な発見

• 一般化されたOsgoodのテストにより、爆発が有限時刻に発生するのは $ B(\infty) < A(\infty) $ であるときであり、爆発時刻は $ T_e = A^{-1}(B(\infty)) $ と与えられる。
• ノイズが $ |\widehat{W}_T| $ の場合、時刻 $ t $ までに爆発する確率は $ P(T_e \leq t) \leq 1 - \Phi\left( \frac{(\alpha a_0 t)^{-1/\alpha} - x_0}{\sqrt{T}} \right) $ で抑えられ、ここで $ T = (\alpha a_0 x_0^\alpha)^{-1} $ である。
• 改良された境界が得られ、$ P(T_e \leq t) \leq 1 - \Phi\left( \frac{x_0}{\sqrt{t}} \left( (\alpha a_0 x_0^\alpha t)^{-1/(1+\alpha)} - 1 \right) \right) $ となる。これは $ \alpha a_0 x_0^\alpha t > 1 $ のとき、以前の推定よりもきつい。
• 比較基準により、$ b $ が増加的で $ a $ が非減少的であれば、初期値が大きい、またはノイズレベルが高い解は、後から爆発する。
• 特定の例として $ b(x) = x^{1+\alpha} $ を取り上げると、95%の信頼度で爆発が発生する臨界時刻は $ t = \frac{1}{\alpha a_0} \left[ x_0 + \sqrt{T} \Phi^{-1}(0.95) \right]^{-\alpha} $ となる。
• 本論文は、ノイズの存在が $ b $ の関数的形に応じて爆発を遅らせたり早めたりする可能性があることを示し、その効果を厳密に確率論的枠組みで定量化する手法を提供する。