[论文解读] 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 $。
• 定义辅助函数 $ A(t) = \int_0^t a(s)ds $ 和 $ B_r(x) = \int_{x_0 - r}^x \frac{ds}{b(s)} $，以推广Osgood的判别准则。
• 应用比较原理，基于 $ b $ 和 $ a $ 的单调性与增长性，关联不同系统解之间的关系。
• 利用 $ B $ 和 $ A $ 的反函数，将爆破时间表达为 $ T_e = A^{-1}(B(\infty)) $，前提是 $ B(\infty) < A(\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 $ 的函数形式，并提供了一个严格的概率框架以量化该效应。