一种新的常微分方程初值问题数值解法.pdfVIP

20 11 8 ( ) Aug. ，20 11 年 月 重庆文理学院学报 自然科学版 30 4 Journal of Chongqing University of Arts and Sciences (Natural Science Edition) Vol. 30 No. 4 第 卷 第 期 一种新的常微分方程初值问题数值解法 钟 巍 ( ， 7 10024 ) 西北核技术研究所 陕西 西安 ［ ］ ， 摘 要 有别于传统的常微分方程数值解法 从一个新的角度出发提供一套新的求解常微分 ; ( )、 Vander Pol 方程初值问题的数值计算方法 通过对一阶常微分方程 组 二阶常微分方程和 ， 0 ． 5 ． 方程的求解验证了方法的正确性 数值解与解析解相对误差小于 % ［ ］ ( ); ;Vander Pol ; ; 关键词 一阶常微分方程 组 显式二阶常微分方程 方程 初值问题 数值解法 ［ ］O24 1． 81 ［ ］A ［ ］1673 － 80 12 (20 11)04 － 00 15 － 04 中图分类号 文献标志码 文章编号 y (x )= y ， 很多科学技术和工程问题用微分方程的形 0 0 并假定初值问题的解全为符号相同 ， ， ， ， ． 式建立数学模型 其中不乏常微分方程 所以 常 的值 本文中全为正值 ． 本文提供如下方法求解该初值问题: 微分方程的数值求解很有实际意义 随着计算机 ， dy f (x ，y ) 计算能力的增强以及数值计算方法的发展 常微 1) (1) = dx ， k 将方程 变形为 并令 分方程的数值求解方法越来越多，比较成熟的有 y y Euler 、 Euler 、 、Runge － Kutta f (x ，y ) 法 后退 法 梯形方法 = dx ，k “ ”; 为 判定因子 ［1］ y 、 、 ， ． 方法 投影法 单步法和多步法 等等 这些方