求解电报方程的自适应重要性抽样蒙特卡罗算法  

Numerical Solution of Telegraph Equation Using Adaptive Importance Sampling Monte Carlo Algorithm

作　　者：洪志敏[1,2] 陈雪[1] 李强[1] 

机构地区：[1]内蒙古工业大学理学院,呼和浩特010051 [2]北京大学数学科学学院,北京100871

出　　处：《应用数学学报》2016年第6期859-870,共12页Acta Mathematicae Applicatae Sinica

基　　金：国家自然科学基金(11461051;11361036);高等学校博士学科点专项科研基金(20131514110005)联合资助课题;内蒙古自然科学基金(2015MS0104)资助项目

摘　　要：提出一种常系数二阶双曲型电报方程的自适应重要性抽样蒙特卡罗算法.通过使用无条件稳定的紧有限差分格式将电报方程离散化为线性代数系统,对得到的线性系统使用具有动态松弛因子的自适应重要性抽样蒙特卡罗算法,加速了蒙特卡罗算法的收敛.一些数值算例的实现证明了提出方法的有效性和适用性.提出的方法容易且适合在计算机上编程实现,所得数值解接近文献提供的精确解.In this paper, the authors proposed a Monte Carlo algorithm based on adaptive importance sampling to find numerical solutions of second order hyperbolic telegraph equations with constant coefficients. Using unconditionally stable compact finite difference formula, we reduce the telegraph equations to a linear algebraic system. The adaptive importance sampling Monte Carlo algorithm with dynamic relaxation factor is applied to this system and accelerating the convergence of the Monte Carlo algorithm. Some numerical ex- amples are considered to demonstrate the validity and applicability of the proposed scheme. It is easy and suitable for computer implementation and obtains numerical solutions close to the exact solutions and available in the literature.

关 键 词：电报方程 紧有限差分 自适应重要性抽样 几乎最优蒙特卡罗算法 

分 类 号：O242[理学—计算数学] O212[理学—数学]