几类力梯度辛算法研究-计算数学专业论文.docxVIP

.jtb . 摘要摘要 摘要 摘要 辛算法由于具有严格保持H锄ilton系统的辛结构、不存在人工耗散等特 点，已经成为数值方法的研究热点。 通常二阶以上的辛算法很难避免使用负步长，但力梯度辛算法却可以全为 正步长。对于可分解为动能T(p)部分和势能V(q)部分的哈密顿函数，当要求每 个算法的四阶误差项系数平方和最小时我们构造了两个全为正系数的非对称的 三阶最优化力梯度辛算法。在数值稳定性和拓扑结构方面，相比非优化方法， 三阶最优化力梯度辛算法并没有很大的优势，但是在计算经典力学模型的能量 误差时，精度优势就会得到明显的体现。同时，本文证明了全为正步长的三阶 最优化力梯度辛算法对于求解某些时间不可逆的量子系统是有利的，特别是在 应用辛．打靶法求解定态薛定谔方程能量本征值问题时，新算法表现出明显的精 度优势。与最优化非力梯度辛算法相比，甚至是与常用的四阶非力梯度辛算法 相比，新算法仍能够保持其优越性。其后，将得到的三阶最优化力梯度辛算法 以对称组合的方式得到两个四阶辛算法。这种对称组合格式在求解摄动Kepler 混沌问题的能量精度和一维定态薛定锷方程的能量本征值精度方面比 Forest．Ruth四阶标准辛算法要好得多，甚至还要明显优越于已有的四阶最优化 力梯度辛算法，因此它们是所谓的四阶最优化力梯度辛算法。 针对一些特殊的哈密顿系统我们构造了两个四阶力梯度辛算法。两个新算 法的每一积分过程中都只包含三个子步，并且在形式上类似于常用的非力梯度 二阶辛算法。首先将它们应用于谐振子模型、数学摆模型和晶体仍模型中，数 值模拟显示：在计算能量误差时，新算法的精度远高于二阶非力梯度辛算法， 明显优越于四阶非力梯度算法，甚至与四阶最优化力梯度辛算法相当。其次， 本文证明了新算法还适用于求解这些特殊哈密顿系统的变分方程。我们利用新 算法求解晶体仍模型的变分方程，继而应用变分法计算某些混沌判别指标，如 LyapuIlov指标和快速Lyapunov指标等。最后，本文应用新算法考察了单一动 力学参数的变化对轨道有序和混沌转化以及两初始参数的变化对有序区域和混 沌区域转变的影响。 Ⅱ 摘要关键词：辛算法；辛打靶法；天体力学：含时薛定谔方程；能量本征值；数值 摘要 关键词：辛算法；辛打靶法；天体力学：含时薛定谔方程；能量本征值；数值 稳定性；二分法；拓扑结构；混沌；Lyapunov指标 In AbstractABSTRACT Abstract ABSTRACT A symplectic irItegratorwim the coIlseⅣation of the s)，111plectic stmcture of a H锄iltoIlian system alld no secular ch锄ge in me eneFgy e玎ors，has becorne a hot research topic． It ghould be emphaSized that beyond second order硒y scheIne does not avoid using some negative time coe伍cients．For me彻：tural Splitting of a H锄iltoniall syStem irno kinetic ene唱y arld potemial ene唱y，we cons臼uct t、)旧new optimal t11ird order force-掣adient symplectic algoritllms，in each of which the no蛐of f10unh-order帆ation enDrs is miIlimized．They are bo也110t eXplicitly superior to meir no-optimal co珊terpans in t11e mlmerical Stabili够arld tlle topology s咖cture．presen，ing，but they are in the accuracy of ene唱y on classical problems 觚d iIl one of tlle enellgy eigenValues for one-dimensional time-inde∞ndent ＆砌舀姆厂 equatioIls．IIl panicul碣tlley are much better maIl me optimal third-order non—gradient symplectic metllod．ney also haVe锄adVamage over me founh-order non·gradient symplectic imegm