基于符号计算的对角隐式辛龙格-库塔算法结构研究.pdfVIP

第39卷第9期赤峰学院学报渊自然科学版冤Vol.39No.9

2023年9月JournalofChifengUniversity(NaturalScienceEdition)Sep.2023

基于符号计算的对角隐式辛龙格-库塔算法结构研究

余英袁卢圳

渊凯里学院袁贵州凯里556011冤

摘要院龙格-库塔算法是数值计算非线性微分方程常用的方法之一遥本文通过在计算代数系统

Sagemath中建立软件包dis_rk.sage袁探讨了对角隐式辛龙格-库塔算法的系数所构成的仿射簇具有的结构

性质遥将所得到的研究结果与已有结果进行比较袁验证所设计软件包的可行性与正确性袁同时也给出了在

软件包下进行的数学实验所得到的新结论遥

关键词院对角隐式辛龙格-库塔算法曰仿射簇曰积分不变量曰Sagemath

中图分类号院O242.1文献标识码院A文章编号院1673-260X渊2023冤09-0031-04

1引言我们不期望的遥换句话说袁可以认为这样的数值解存

随着科技的进步袁非线性现象在生活和生产领在巨大的误差遥早在1928年袁Courant,Friedrichs

域有着广泛的应用遥很明显袁非线性现象较线性现和lewy就在构造算法时渗透了保持能量守恒的思

象更能反映实际模型遥鉴于目前许多非线性微分方想袁只不过当时的这个保持的量是一个正定的量袁还

程的解析解很难求得袁因而学者们转向数值研究这没有能量守恒这种提法遥在国内袁据了解袁学者郭本

[4-6]

类模型遥对非线性微分方程的数值计算目前没有统瑜是比较早构造具有这种结构算法的学者之一遥

[7]

一的算法袁龙格-库塔算法(Runge-Kuttamethod)是再后来袁冯康先生提出了辛几何算法思想袁并利用

非线性方程数值计算最常用的一种方法遥在对非线生成函数与辛方法的联系袁证明了任意的高阶辛算

性方程的数值模拟中袁由于隐式的龙格-库塔算法法是存在的这个结论遥该思想在国内引起了极大的

在每一步的计算中袁都要求一个非线性方程组的反响袁直接推动了保结构算法的研究遥构造保结构

解遥这不仅增加了计算成本袁而且求非线性方程组算法的原则是问题的模型在离散化后其基本特征

[1][8-11]

的解至今没有一致且有效的算法袁因而给计算带应当得到最大限度地保持遥随后袁保结构算法的

来不便遥在这一点上袁显式的龙格库塔法就比较受构造成为世界前沿的问题之一遥近年来袁保结构算

欢迎遥它被设计成专门的程序袁并在各种计算代数法已在量子力学尧流体力学尧结构力学尧水动力学袁

系统中都有实现它的软件包遥电磁学尧地球物理学等科学和工程领域得到了广泛

近三十年来袁构造守恒的算法受到了广大学者的应用遥

的高度重视