几类力梯度辛算法的研究.pdfVIP

摘要 摘要 辛算法由于具有严格保持H锄ilton系统的辛结构、不存在人工耗散等特 点，已经成为数值方法的研究热点。 通常二阶以上的辛算法很难避免使用负步长，但力梯度辛算法却可以全为 正步长。对于可分解为动能T(p)部分和势能V(q)部分的哈密顿函数，当要求每 个算法的四阶误差项系数平方和最小时我们构造了两个全为正系数的非对称的 三阶最优化力梯度辛算法。在数值稳定性和拓扑结构方面，相比非优化方法， 三阶最优化力梯度辛算法并没有很大的优势，但是在计算经典力学模型的能量 误差时，精度优势就会得到明显的体现。同时，本文证明了全为正步长的三阶 最优化力梯度辛算法对于求解某些时间不可逆的量子系统是有利的，特别是在 应用辛．打靶法求解定态薛定谔方程能量本征值问题时，新算法表现出明显的精 度优势。与最优化非力梯度辛算法相比，甚至是与常用的四阶非力梯度辛算法 相比，新算法仍能够保持其优越性。其后，将得到的三阶最优化力梯度辛算法 以对称组合的方式得到两个四阶辛算法。这种对称组合格式在求解摄动Kepler 混沌问题的能量精度和一维定态薛定锷方程的能量本征值精度方面比 Forest．Ruth四阶标准辛算法要好得多，甚至还要明显优越于