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

物理计算的保真与代数动力学算法物理计算的保真与代数动力学算法 物理计算的保真与代数动力学算法物理计算的保真与代数动力学算法 II. 截断近似与代数动力学算法截断近似与代数动力学算法 截断近似与代数动力学算法截断近似与代数动力学算法 王顺金王顺金 张华张华 王顺金王顺金 张华张华 四川大学四川大学 物理系物理系, 成都成都 610064 四川大学四川大学 物理系物理系 成都成都 摘要摘要 摘要摘要 代数动力学算法的目的是要给出常微分方程的数值解. 代数动力学算法的目的是要给出常微分方程的数值解 在常微分方程的代数动力学在常微分方程的代数动力学 代数动力学算法的目的是要给出常微分方程的数值解代数动力学算法的目的是要给出常微分方程的数值解 在常微分方程的代数动力学在常微分方程的代数动力学 精确解的基础上精确解的基础上, 对泰勒级数表示的局域收敛的精确解做有限项截断对泰勒级数表示的局域收敛的精确解做有限项截断, 从而设从而设 精确解的基础上精确解的基础上 对泰勒级数表示的局域收敛的精确解做有限项截断对泰勒级数表示的局域收敛的精确解做有限项截断 从而设从而设 在 N 阶截断近似下阶截断近似下, 估计了相空间轨道, 运运 计出新的常微分方程的代数动力学算法计出新的常微分方程的代数动力学算法. 在 估计了相空间轨道 计出新的常微分方程的代数动力学算法计出新的常微分方程的代数动力学算法 在在 阶截断近似下阶截断近似下 估计了估计了相空间轨道相空间轨道 运运 动学变量的代数关系动学变量的代数关系(包括哈密顿系统的辛几何结构包括哈密顿系统的辛几何结构)和代数和代数-几何不变量以及动力学守恒量几何不变量以及动力学守恒量 动学变量的代数关系动学变量的代数关系包括哈密顿系统的辛几何结构包括哈密顿系统的辛几何结构和代数和代数几何不变量以及动力学守恒量几何不变量以及动力学守恒量 如哈密顿系统的能量) 的的精度 对哈密顿系统, 提出了提出了改进辛几何结构精度或 (如哈密顿系统的能量 精度. 对哈密顿系统 改进辛几何结构精度或改进能量精度改进能量精度 如哈密顿系统的能量如哈密顿系统的能量 的的精度精度 对哈密顿系统对哈密顿系统 提出了提出了改进辛几何结构精度改进辛几何结构精度或或改进能量精度改进能量精度 的计算方案. 的计算方案 的计算方案的计算方案 [1] 前文 ( 以后称文(I) ) 详细地介绍了常微分方程的代数动力学解法, 给出了自治的常 微分方程的用泰勒级数表示的局域收敛的精确的解析解. 本文的目的是要解决常微分方程 的数值求解问题. 这是本文要介绍的代数动力学算法的任务. 代数动力学算法的目的是要给出常微分方程的数值解. 它是基于常微分方程的代数动 力学精确解. 对泰勒级数表示的局域收敛的精确解做有限项截断, 从而设计出新的常微分方 程的代数动力学算法. 为此, 我们需要前一篇文章的与本文有关的主要结果. [1] I. 常微分方程的代数动力学精确解常微分方程的代数动力学精确解 常微分方程的代数动力学精确解常微分方程的代数动力学精确解 1-1. 常微分方程的结构和常微分方程的结构和代数动力学精确解代数动力学精确解 常微分方程的结构和常微分方程的结构和代数动力学精确解代数动力学精确解 [2,3] 多变量常微分方程及其结构 : n 个变量 多变量常微分方程及其结构 p 个参数的常微分方程可写为