论文部分内容阅读
本文提出了一种高阶多辛分裂方法,在解决一类确定性和随机哈密尔顿系统方程时,具有很重要的作用。它在应用过程中包含了高阶差分算法、多辛算法和分裂算法,是一种很有效的复合差分方法。本文把这类方法应用到了确定性非线性薛定谔方程和随机非线性薛定谔方程的数值计算中,并通过对孤立波和碰撞波的数值模拟,从中可以发现此类方法在解决随机问题和确定性问题时体现出的优越性。简单来讲,首先针对随机非线性薛定谔方程和确定性的非线性薛定谔方程,我们分别构造了相对应的哈密尔顿多辛结构形式;接着在原有问题的基础上,通过使用二步分裂方法或者三步分裂方法,把它分裂成了一个线性子问题和一个或者两个非线性子问题,并且得到了相对应的二步多辛分裂格式和三步多辛分裂格式;然后在面对线性子问题时,应用了四阶差分方法,而面对非线性子问题时,直接给出了其相对应的数值显式表达式;最后对得到的数值差分结构,我们给出了满足随机非线性薛定谔方程的离散电荷守恒律和能量递推公式,以及满足确定性情况下的离散电荷守恒律和离散能量守恒律。 本文的最后部分给出了相对应的数值模拟实验,证明了该算法的收敛精度是四阶的,且分别观察了不同随机扰动项系数对孤立波和碰撞波的传播影响情况。此外,通过数值实验得到的数据,用图展现了满足这两种波的电荷守恒律和能量递推公式,较符合理论结果。