论文部分内容阅读
本文主要研究Stratonovich意义的乘性噪声驱动的随机Schr(o)dinger方程数值方法.基于随机变分原理,提出无穷维随机Hamilton系统的一般形式,并应用于随机Schr(o)dinger方程,给出其无穷维随机Hamilton系统的表达式,证明了它的解流保持无穷维随机辛性.对于随机Schr(o)dinger方程,电荷守恒定律仍然成立,但是能量却不再守恒,而是满足随机演化公式[6].利用该公式及随机分析的工具,证明能量具有一致有界性.该性质蕴含着d-维随机Schr(o)dinger方程解的H1(O)-稳定性.特别地,对于一维问题,证明了解的H2(O)-稳定性.同时,对于d-维问题(d=1,2,3),我们证明了解在L2(O)-范数下的H(o)lder连续性,特别地,在一维情形,解在H1(O)-范数下也具有H(o)lder连续性. 为构造保持无穷维随机辛性的数值方法,给出一类在时间方向上半离散随机Schr(o)dinger方程的随机Runge-Kutta方法,并得到该方法保持无穷维随机辛性的条件,证明了满足该条件的随机Runge-Kutta方法保持离散的电荷守恒定律.给出了随机Schr(o)dinger方程时间半离散格式的基本收敛性定理,该定理表明一个格式的全局均方阶只依赖于局部误差和算子半群的性质,它是全局均方收敛阶的一般性判别准则.基于该定理,以中点格式为例分析了它的均方收敛阶,验证了定理的有效性. 因随机Schr(o)dinger方程的能量具有一致有界性,为使数值方法对应的离散能量继承该有界性,提出一类θ-方法,并利用该方法在时间方向上离散随机Schr(o)dinger方程.为了产生足够的数值耗散性以控制离散噪声的影响,这里θ在区间[1/2+O(√(τ)),1]中取值.证明了对于一般的d-维问题,这类方法对应的离散能量具有一致上界,该上界只依赖于初始值ψ0、终端时间T和问题的区域.离散能量的一致有界性蕴含着数值解的H1(O)-稳定性,基于此利用紧致理论得到数值解收敛到原问题温和解的结论.为了更进一步刻画该方法的收敛性,即得到其收敛阶,对于一维问题,我们利用解的稳定性、H(o)lder连续性、Sobolev嵌入定理以及随机分析的工具证明了当初始值ψ0∈L8(Ω;H10∩H2(O))时,数值方法的收敛阶为1/2;当初始值ψo∈L6(Ω;H△∩H4(O))时,数值方法的最优收敛阶为1.