具有指数型激励的非齐次自治系统的一种高精度精细算法.pdf

万方数据 万方数据 万方数据 具有指数型激励的非齐次自治系统的一种高精度精细算法 摘 要 自钟万勰院士1994年提出齐次线性自治动力系统的精细算法 HPD 以来，这一计算力学、工程应用与计算数学的学术交叉点迅速发 展，已成为学术热点。本文基于已有的研究成果，围绕多项式展开技 术，主要是Chebyshev正交多项式展开技术将精细算法引向深入,开展 了含指数型长效精细算法的系列性研究。在实际问题中，经常碰到含 指数激励项的非齐次方程，因此研究这类方程的精细算法十分必要。 本文的创新工作主要有以下四个方面： （1）新的外推方法及初始N值确定定理 本文从精细算法推导过程入手，在分析其误差来源的基础上，通 过对不同N值对应的传递矩阵H做适当的线性组合来研究一种新的精 细算法的外推，并且在一阶外推的基础上，给出了二阶外推的表达式 以及更高阶外推的外推方法。该外推方法不同于传统的外推方法：不 x 是对结果 进行线性组合，而是对传递矩阵H本身进行外推线性组合， 因此称为内置型外推方法，记为NW-HPD。算例表明，NW-HPD的计算精 度相比于HPD有显著提高。 另外在分析误差来源时，使误差项造成的相对误差小于给定的相 对误差最大值，给出了初始N值确定定理，算例表明，用该定理选取 的N值恰好能够满足给定的精度要求，这样的N值是十分理想的。 （2）带指数激励的第二类Chebyshev正交多项式展开技术 万方数据 本文设计的算法记为HHPD-SC，该算法把含有指数型激励的右端 函数按照第二类Chebyshev正交多项式经改进变换为含有指数项的特 殊多项式展开。文中给出了构造HHPD-SC的具体方法，并对该方法的 精度及计算量给出了分析。理论及算例表明，HHPD-SC相对于龙格库 塔方法的计算精度及计算量，优势是十分理想的，且HHPD-SC属于长 效精细算法，即步长一定时，其传递矩阵保持不变，可多次重复使用。 （3）TB-HHPD-SC算法 在实际应用中，时常见到测量型右端激励，本文针对含有指数项 的测量型右端激励，设计了拓展型可并行计算的长效精细算法 TB-HHPD-SC,理论与算例表明：本文的TB-HHPD-SC算法十分有效。 （4）PS-HHPD-SC算法 在系统矩阵A为实Hamilton矩阵时，考虑到Hamilton 系统辛几何 算法的优良特性，巧妙地将Hamilton 系统辛几何算法与HHPD-SC算法 结合，构造了PS-HHPD-SC算法，算例表明，PS-HHPD-SC算法的计算结 果更为精确。 关键词：Chebyshev正交多项式，外推精细算法，齐次扩容精细算法， 辛算法，Hamilton 系统 注：本文中的定义、定理、命题、算例、算法均指本文作者的 工作。 万方数据 A High-accuracy High Precision Direct Algorithm to Solve Non- homogeneous Autonomy System with Exponent Stimulus. ABSTRACT Since Zhong Wanxie, a well-known academician proposed the High Precision Direct(HPD) integration scheme in 1994, this scheme, which is