非线性动力学常微分方程组高精度卷积分方法 I.理论和数值方法

Abstract

建立了一种求解非线性动力学常微分方程组初值问题的新方法. 若非线性函数一阶导数存在，则给出解的积分方程表达式，计算得到按规定误差要求的高精度数值解. 引入一般自治或非自治非线性系统的首次近似Jacobi矩阵,不作任何假设重构等价的非线性常微分方程组，简捷而有广泛的适应性，不改变方程的本质，但其主项构成线性化方程组，其它项则代表非线性函数高阶余项而不涉及Taylor级数展开计算，给出该方程组初值问题的Duhamel 卷积分解析表达式，在时间步长内进行数值积分迭代求解,在指定误差内快速收敛，逐步递推获得非线性常微分方程的瞬态响应和全时域高精度数值解. 积分解连续满足微分方程组而不是在离散的步长端点上满足代数方程组，打破了传统用增量法在离散点上建立的代数方程组迭代求解，从而使传统Euler型逐步积分法的各种差分格式算法改变成真正的积分格式算法. 数值计算中给出指数矩阵递增展开式，变矩阵乘法为乘积系数的加法，避免了大量矩阵自乘而大大提高计算效率. 算法验证为无条件稳定，则保证对线性常微分方程而言，计算中舍入误差的传播不会扩散，不出现计算机字长有限而引起舍入误差导致计算不确定性问题. 基于以上理论和数值方法, 计算了线性非线性算例并进行了分析, 验证了本方法简捷而有广泛的适应性，可以有足够的精确性.