求解薛定谔–泊松方程组的时间分裂紧致差分格式

本文通过应用高阶紧致差分格式,时间分裂法与Crank-Nicolson法等方法求解非线性薛定谔–泊松方程组。基于快速Sine变换,我们设计了一种求解离散系统的快速算法。我们用该算法分别求解一维、二维、三维的非线性薛定谔–泊松方程组。我们列举具体的数值例子,使用MATLAB软件编写出算法程序,并且通过程序计算近似误差与画出近似解的图像。数值计算结果证实了该算法在空间方向具有谱精度,也证实了该算法的高效性与稳定性。