传统隐式Runge-Kutta方法转换方法.pdfVIP

第45卷第 16期 数学的实践与认识 Vo1．45．NO．16 2015年 8月 MATHEMATICSIN PRACTICEAND THEORY Aug．，2015 传统隐式 Runge-Kutta方法的转换方法 汪芳宗，廖小兵，谢 雄 (三峡大学电气与新能源学院，湖北 宜昌443002) 摘 要：多级隐式Runge—Kutta(RK)方法簇中，除Gauss类方法是 s级2s阶的辛 方法以外，Radau类方法和Lobatto类方法既不是 s级2s阶的方法也不是辛方法． 基于隐式RK方法是一类转换 RK方法这一特征，利用 V一变换和Pa 对角逼近， 提出了构造高阶RK方法的转换定理．依据转换定理，导出了s级 2s阶的Radau 方法和S级2s阶的Lobatto方法．利用V一变换和待定系数法，导出了辛Radau方 法和辛 Lobatto方法．在此基础上，发现并证明了辛Radau方法是 s级 2s阶的方 法． 关键词：隐式RK方法；转换RK方法；V一变换；Pad@逼近；辛方法 1 引言 工程领域中绝大多数随时问演变的物理现象或过程，主要是采用各种复杂的常微分方程 来描述．有关常微分方程数值计算方法的研究，既是数值分析计算领域中较古老的一个分支， 也是众多工程计算领域中经常涉及的一类重要课题． 早期的常微分方程数值计算方法主要是采用数值积分方法．迄今为止，研究人员已提出 了种类繁多的数值积分方法，其中最有代表性的数值积分方法应该是 Runge—Kutta(RK)方法 簇 [1--3]．早在上世纪60年代，为研究多级隐式RK方法的计算精度即阶数，Butcher提出了根 树 (rootedtrees)结构理论，并在此基础上提出了著名的Butcher基本阶定理 [1{3]．上世纪 70 年代，Hairer和Wanner在根树结构理论的基础上，提出了B 级数 (B—series)的概念，并证明 了B一级数的组合定理，应用这种定理可以建立各类RK方法的阶条件 (orderconditions)[ ]． 关于多级隐式 RK方法的线性稳定性分析，研究人员应用指数函数的Pad6逼近理论 [1J3]，证 明了传统的多级隐式RK方法均是A一稳定的． 上述研究结果完整地建立了多级隐式 RK方法的理论基础体系．因此，RK方法发展到 上世纪 80年代末，可以说基本上画上了一个句号．然而在 1984年，冯康先生首次系统地提 出了哈密尔顿 (Hamilton)系统的辛几何算法 [4-6](symplecticgeometrymethod)．与传统的 数值积分方法相比较，辛算法在数值稳定性与长时间跟踪能力方面具有独特的优势．受冯康 先生相关研究结果的启发，Sanz—Serna于 1988年率先提出并证明了RK方法是辛方法的充 要条件 _7J1由此引发了研究辛 RK方法的热潮．随后，孙耿利用 Hairer所提出的 w一变换 (w —transformation)[01导出了 s级 2s一1阶的辛RadauIB和辛 RadauIIB方法 Is-9J同时发 现 LobattoIIIE方法是 8级 2s一2阶的辛方法 _8J．迄今为止，研究人员已将 RK方法辛几何 收稿 日期：2015—01一i0 资助项 目：国家 自然科学基金 51377098) 16期 汪芳宗，等：传统隐式Runge—Kutta方法的转换方法 173 化，建立了完整的辛 RK系列方法． 本文分析、推导了转换 RK方法 [10--11】的特性，在此基础上提出了构造高阶RK方法的 转换定理．依据转换定理，导出了s级 2s阶的Radau方法和 8级 2s阶的Lobatto方法．利用 一 变换 [12--13]和待定系数法，导出了辛Radau方法和辛Lobatto方法．有趣的是，本文利用 V一变换所导出的辛 Radau方法与孙耿利用 w一变换所导出的辛 Radau方法是一致的，所 导出的辛 Lobatto方法即是 LobattoIIIE方法．但从计算结果发现：辛Radau方法与本文利 用Pad6；~角逼近所导出的s级2s阶的Radau方法是一致的．由此可知：Ra