振荡微分系统保结构算法的研究.pdf

摘 要 摘摘摘 要要要 微分方程保结构算法的研究始于上世纪 年代 其基本思想是 在处理振荡微分系统 解的逼近时 不仅要求数值解从定量 的角度能够逼近原有系统的精确解 还期望数值解从 定性的角度模拟精确解的性态 尽可能地保持微分系统本身的物理结构和几何特征 保结 构算法以其独特的优势 被越来越多的科研工作者认可 井开启了关于高精度保结构算法 的研究 本文的主题是振荡微分方程的保结构算法 基于指数型求积公式和矩阵常数变易公 式 针对一阶和二阶非线性高振荡的微分方程 分别构造井分析了一类显式伪两步指数型 算法和一类显式两步 算法 论文内容 由以下三部分构成 第一部分 介绍振荡微分系统保结构算法的发展背景以及研究现状 第二部分 针对一阶高振荡的微分方程 结合显式伪两步 方法和指数型 方法 在满足阶条件 的前提下 构造了一类显式伪两步指数型 算法 分析了算法的收敛阶 数值实验的结果不仅验证了算法的效能 还展示 出该算法具有较好的长时间行为 第三部分 针对二阶高振荡的微分方程 结合矩阵常数变易公式和 插 值多项式 构造了一类可达任意高阶的显式两步 时间积分算法 井分析、讨论了该算法的局部截断误差、非线性稳定性和收敛性 将数值算法分别应用于 非线性波方程 包括 方程、 方程等 对比文献 中己有 的数值算 法 验证了算法的高效性 关关关键键键同同同 振荡微分方程 保结构算法 指数型求积公式 矩阵常数变易公式 误差分析 目 录 目目目 录录录 摘要 第一章 引言 第二章 一类显式伪两步指数型 方法 定义及阶条件 全局误差 方法的推导 数值实验 第三章 一类显式两步 时间积分方法 矩阵值函数及其性质 定义及局部截断误差 非线性稳定性 收敛性 数值实验 第四章 总结