代数动力学算法(3b-4).pdfVIP

物理计算的保真与代数动力学算法物理计算的保真与代数动力学算法 物理计算的保真与代数动力学算法物理计算的保真与代数动力学算法 III.代数动力学算法与辛算法和龙格代数动力学算法与辛算法和龙格-库塔算法的关系库塔算法的关系 代数动力学算法与辛算法和龙格代数动力学算法与辛算法和龙格 库塔算法的关系库塔算法的关系 王顺金王顺金 张华张华 王顺金王顺金 张华张华 四川大学四川大学 物理系物理系, 成都成都 610064 四川大学四川大学 物理系物理系 成都成都 摘要摘要 摘要摘要 从代数动力学算法和李群 李代数的观点, 分析分析 比较了代数动力学算法与辛算法和龙 从代数动力学算法和李群、李代数的观点 、比较了代数动力学算法与辛算法和龙 从代数动力学算法和李群从代数动力学算法和李群 李代数的观点李代数的观点 分析分析 比较了代数动力学算法与辛算法和龙比较了代数动力学算法与辛算法和龙 ˆ ˆ 格格-库塔算法及其库塔算法及其关系 从这三种算法对于微分生成元 和对于精确的时间演化群元 关系. 从这三种算法对于微分生成元L 和对于精确的时间演化群元U(t) 格格 库塔算法及其库塔算法及其关系关系 从这三种算法对于微分生成元从这三种算法对于微分生成元 和对于精确的时间演化群元和对于精确的时间演化群元 的偏离情况, 讨论了讨论了三种算法的精度 分析表明: 辛算法为了保持辛几何结构付出了动力学 的偏离情况 三种算法的精度. 分析表明 辛算法为了保持辛几何结构付出了动力学 的偏离情况的偏离情况 讨论了讨论了三种算法的精度三种算法的精度 分析表明分析表明 辛算法为了保持辛几何结构付出了动力学辛算法为了保持辛几何结构付出了动力学 ˆ 失真的高失真的高昂代价昂代价 库塔算法由于破坏微分生成元 的完整性而产生人为耗散 , 龙格龙格-库塔算法由于破坏微分生成元L 的完整性而产生人为耗散, 代数动代数动 失真的高失真的高昂代价昂代价 龙格龙格库塔算法由于破坏微分生成元库塔算法由于破坏微分生成元 的完整性而产生人为耗散的完整性而产生人为耗散 代数动代数动 力学算法具有更高的精度并且能在很高的精度下兼顾运动学代数力学算法具有更高的精度并且能在很高的精度下兼顾运动学代数-几何保真和动几何保真和动 力学算法具有更高的精度并且能在很高的精度下兼顾运动学代数力学算法具有更高的精度并且能在很高的精度下兼顾运动学代数 几何保真和动几何保真和动 力学保真二者力学保真二者. 力学保真二者力学保真二者 [1 −2] [3 −9] [5 −6] 代数动力学算法 是独立于辛算法 和龙格-库塔算法 的第三种算法. 这三种 算法都致力于常微分方程[10 −12] 的数值求解. 研究表明, 它们之间既有密切的关系. 又有实 质性的区别. 本文正是要研究这种联系和区别. [1 −2] I. 代数动力学算法代数动力学算法 代数动力学算法代数动力学算法 [1] 1-1 常微分方程的代数动力学精确解常微分方程的代数动力学精确解 常微分方程的代数动力学精确解常微分方程