求非线性方程近似解的通用程序.pdf

第 21 卷 　第 3 期 　　　　　　　　　　　　　 　　　高 师 理 科 学 刊 　　　　　　　　　　　　　　　Vol. 21 　No. 3 2001 年 　8 月 　　　　　　　　　Journal of Science of Teachers′College and University 　　　　　　　　　　Aug. 　2001 ( ) 　　文章编号 :1007 - 9831 2001 03 - 0013 - 03 求非线性方程近似解的通用程序 1 2 赵喜清 ,徐永春 ( 1 张家口师专 计算机系 ,河北 张家口 075000 ;2 张家口师专 数学系 ,河北 张家口 075000) 摘要 :根据牛顿迭代法思想 ,用 MA TLAB 语言给出求非线性方程近似解的通用程序. 关键词 :MA TLAB ;牛顿迭代法 ;数学实验 中图分类号 :O17514 　　　文献标识码 :A 1 　核心思想 ———牛顿迭代法 ( ) ( ) ( ) 　　设f x = 0 为非线性方程, f 为[ a , b] 上二阶可导函数, 且 f ′x , f ″x 均不为零, 不妨设f ′ 0 , ( ) ( ) f ″ 0 即f 为[ a , b] 上递增凸函数, 又设 f a 0 , f b 0 , 那么由连续函数的介值性定理知, 方程 ( ) ξ ξ f x = 0 在[ a , b] 内存在唯一解 . 牛顿迭代法的基本思想是构造一点列{ X n} , 使lim X n = 且当n 充分 n →∞ ξ ( ) 大时, 可以 X n 作为 的近似值, 即方程 f x = 0 的近似解. ( ) ( ) 　　其中 　X 0 = b 　　X n = X n - 1/ f ″X n - 1 　　 n = 1 , 2 , 3 … ξ ξ ( ) 以 X n 中为 的近似值的误差估计为 | X n - | f X n / m | ( ) 　　其中 m = min { | f ′x | } x ∈[ a , b] 2 　程序设计思路及其算法 2 . 1 　用图形放大法确定迭代函数根所在区间[ a , b] 　　MA TLAB 提供了强大的图形功能, 二维图形用plot ( x , y) 函数实现, 其中plot 命令打开一个称为图形 窗口的窗口, 将坐标轴扩以适应数据, 绘制数据, 然后以 x , y 为横纵坐标绘制连线图. 　　调用上述函数可以非常方便地作出迭代函数 f ( x) ( ) X 轴的交点坐标的横 的图形 曲线