粒子自旋4阶辛格式计算误差研究

摘要

利用辛算法求解粒子自旋问题的薛定谔方程,得到波函数的数值解.研究了4阶辛差分格式计算结果的误差,并与2阶辛差分格式的结果进行了比较.利用4阶辛格式计算波函数实部和虚部结果的精密度比2阶格式高出7个数量级,即绝对误差低7个数量级,但二者演变规律基本相同,即绝对误差随着时间的推演均周期性地在正数和负数之间来回变动,变化方式类似于正弦、余弦函数,其振幅不断增大.4阶辛格式结果误差的变化图形较2阶辛格式略微滞后.2阶的绝对误差随时间的变化恰好与波函数本身的时间变化率成正比,即波函数绝对误差与其时间变化率的比值随时间的变化呈严格的直线图像,而4阶辛格式结果没有这样的关系.但若考虑到4阶绝对误差在时间上的滞后,也能够变换出类似的直线关系.