第七章非性方程解法.docVIP

新疆大学 - PAGE 8 - 数学与系统科学院 第七章 非线性方程解法 2. Newton-下山法. 每次迭代在改变量前加一因子以保证收敛: x n+1=xn-λn f(xn)/f′(xn) 这儿λn在0,1间,可用各种方法搜索,例如用分半法取1,1/2,1/4,…试探,使 下山条件∣f(xn+1)∣∣f(xn)∣成立为止。 牛顿下山实验⑷ 上机题目:牛顿下山法上机实验 实验目的:编制求单变量非线性方程组的程序. 实验要求: ①上机前充分准备，复习有关内容，写出计算步骤，查对程序； ②完成实验后写出完整的实验报告，内容应该包括：所用的算法语言，算法步 骤陈述，变量说明，程序清单，输出计算结果，结果分析等等；③用编好的程 序在Ｍatlab环境中执行。利用Newton下山法来解方程； 计算步骤：①准备 选定初始近似值x SKIPIF 1 0 ，计算f SKIPIF 1 0 =f(x SKIPIF 1 0 ), SKIPIF 1 0 SKIPIF 1 0 。 ②迭代 按公式 x SKIPIF 1 0 迭代一次，得新的近似值x SKIPIF 1 0 , 计算 SKIPIF 1 0 SKIPIF 1 0 .最初 SKIPIF 1 0 =1，之后奖 SKIPIF 1 0 减半进行试算直到下山条件 ∣f(xn+1)∣∣f(xn)∣成立。 ③控制 如果满足∣f(xn+1)∣∣f(xn)∣则终止迭代，；否则转步骤4； ④修改 如果迭代次数达到预先指定的次数N，或者 SKIPIF 1 0 =0，则方程失败；否则以（x SKIPIF 1 0 ， SKIPIF 1 0 ）代替（ SKIPIF 1 0 ）转步骤2继续迭代。 算法例题：用牛顿下山方法解方程 x SKIPIF 1 0 -x-1=0，取迭代初值 x SKIPIF 1 0 =1.5, d=10 SKIPIF 1 0 . Newton下山法的Matlab程序： function x=newton_xiashanfa(f,x0,d,max) y=diff(f); %取导数 y=inline(y); %定义y f=inline(f); %定义f x(1)=x0; l=1; disp(k l x ); %以指定格式输出k,x. for k=1:max x(k+1)=x(k)-l*f(x(k))/y(x(k)); % 计算公式 if abs(f(x(k+1)))abs(f(x(k))) %判断是否小于d if abs(x(k+1)-x(k))d break end else if abs(f(x(k+1)))abs(f(x(k))) l=l/2; end w=k; end for k=1:w disp(sprintf(%d %f %10f ,k,l,x(k))) %输出算结果 end x=x(k+1); 运行结果： x=newton_xiashanfa(x^3-x-1,1.5,10^(-6),3) max =3 k x f(x(K+1)) f(x(k)) 3 1.325213 0.662708 0.463111 x = 1.5000 1.3478 1.3252 1.3247 Newton下山法的流程图:： 流程图解释：1） 输入 SKIPIF 1 0 ；2）把1 赋给 SKIPIF 1 0 ;3）把 SKIPIF 1 0 赋给 SKIPIF 1 0 ；4）判断 SKIPIF 1 0 ，若 SKIPIF 1 0 那么表明下山成功，停止；若 SKIPIF 1 0 ，那么到步骤5）；5）判断 SKIPIF 1 0 ，若