二阶非线性振动方程的两种混合数值.pptVIP

二阶非线性振动方程的两种混合数值格式 吴奇 ，化存才 云南师范大学数学学院，昆明, 650092； artlass103@， cuncai-hua@ 二阶非线性振动方程的两种混合数值格式 引言 混合数值格式的构造思想 “中点数值积分/差商”混合的数值格式 “中点数值积分/中心差商”混合的数值格式 数值结果 结论 引言 二阶非线性振动方程 渐近方法: KBM 数值方法: Euler, Runge-Kutta, 有限元法, 线性多步法 , 对称线性多步法 混合数值格式的构造思想 初值问题 算子方程形式 混合数值格式的构造思想 将二次求导运算进行分离，不同次的求导运算分别采用不同的数值微积分近似公式. “中点数值积分/差商”混合的数值格式 一阶微分方程组 混合数值格式:第一个方程用中点数值积分， 第二个方程用向前差商公式 隐式格式 显式格式 “中点数值积分/差商”混合的数值格式 隐式格式的局部误差 隐式格式的整体误差 显式格式的整体误差: 同隐式格式 “中点数值积分/中心差商”混合的数值格式 一阶微分方程组 混合数值格式:第一个方程用中点数值积分， 第二个方程用中心差商公式 两步显式格式: “中点数值积分/中心差商”混合的数值格式 格式的局部误差 格式的整体误差 数值结果 显式格式 Van der Pol振动方程 结果与二级二阶Runge-Kutta法相同 结论 通过二阶微分算子运算的分离,将二阶非线性振动方程化为一阶方程组,对两个方程的导数分别运用不同的数值积分和差商公式，得到了解二阶非线性振动微分方程的两种混合数值格式，不仅格式易于显式化，而且局部截断误差优于Euler方法，整体误差分析和算例结果都表明：混合数值方法与Runge-Kutta法相同 。 利用适当的数值方法进行巧妙地混合组合，就可能得到许多关于非线性振动方程的高性能算法。 参考文献 李家春,周显初. 数学物理中的渐近方法[M]，科学出版社，2002 余德浩,汤华中.微分方程数值解法[M]. 科学出版社，2003 P．亨利西著，包雪松等译. 常微分方程离散变量方法[M].科学出版社，1985 J.D.Lambert I.A.Watson. Symmetric multi-step methods for periodic initial value problem [J]. J.Inst.Math.Appl.,18 (1976),189-202 李世航，肖箭. 关于差分不等式的估计和计算[J]. 大学数学，2005,21(2):98-101