简化牛顿法.docVIP

第七章 非线性方程解法 ⒊Newton法 ①????方法 构造近似解的一个思路是‘以直代曲’,‘以线性函数代非线性函数’. Newton法是在初始近似x0作Taylor展开: f(x)= f(x0)+f′(x0)(x-x0)+f″(η)(x-x0)2/2 略去余项,取其零点为新近似x1, x1=x0-f(x0)/f′(x0) 一般地,对n=0,1,2,…计算 x n+1=xn-f(xn)/f′(xn) 直到∣x n+1-xn∣≤ε.这就是Newton法. 几何上也就是由(x0,f(x0))作f(x)的切线,交x轴于 x1.因此,Newton法也叫切线法. 同例2,取x0=1.5,用Newton迭代解方程. 解：取迭代初值x0=1.5，迭代公式x(k+1)=x(k)-f(x(k))/y(x(k)); 计算结果见下表. n xn 0 1.50000000000000 1 1.34782608695652 2 1.32520039895091 3 1.32471817399905 4 1.32471795724479 5 1.32471795724475 6 1.32471795724475 n xn 0 1.50000000000000 1 1.41666666666667 2 1.41421568627451 3 1.41421356237469 4 1.41421356237309 5 1.41421356237309 例６：用Newton迭代解方程x2-2=0的根. 解：取迭代初值x0=1.5，迭代公式xn+1=(xn+2/xn)/2。计算结果见下表. ②???? 收敛性 Newton迭代法在f(x)=0单根ξ附近是二阶收敛的,并且有 ∣x n+1-ξ∣≈∣f″(ξ)/(2f′(ξ))∣∣x n-ξ∣2 由 0= f(ξ)=f(x n)+f′(x n)(ξ-x n)+f″(η)(ξ-xn)2/2 得 x n+1-ξ= f″(η)/(2f′(x n))(ξ-xn)2 据此可推出在ξ附近取x0,误差越来越小,是二阶收敛的. 注意到Newton迭代法也可视为不动点迭代法 φ(x )=x-f(x)/f′(x) 从而应用其结果,亦得二阶收敛。在这里单根的假设是必要的.例如,在求(x-1)2=0的二重根时.Newton迭代 x n+1=x n-(x n-1)/2=(x n+1)/2 x n+1-1=(x n-1)/2 是线性收敛的. ③???? 变形 Newton迭代法有多种变形. 简化Newton法. 为减少计算导数的化费,可只求f′(x0)以后所有导数 值就用它,不另求.这相当于第一次作切线,以后作其平行线.当然,这样收敛要慢些.还可以取折衷方案,隔几步计算导数. 简化牛顿实验⑶ 上机题目:简化牛顿实验Newton法 实验目的:编制求单变量非线性方程组的程序. 实验要求: ①上机前充分准备，复习有关内容，写出计算步骤，查对程序； ②完成实验后写出完整的实验报告，内容应该包括：所用的算法语言， 算法步骤陈述，变量说明，程序清单，输出计算结果，结果分析等； ③用编好的程序在Ｍatlab环境中执行简化Newton法.来解方程； 计算步骤：①准备 选定初始近似值x，计算f=f(x),。 ②迭代 按公式 x迭代一次，得新的近似值x, 计算. ③控制 如果 x满足或，则终止迭代，以x作为所求的根；否则转步骤4，此处是允许误差，而 时。 其中c是取绝对误差或相对误差的控制常数，一般可取c=1. ④修改 如果迭代次数达到预先指定的次数N，或者=0，则方程失败；否则以（x，）代替（）转步骤2继续迭代。 例题：用简化牛顿方法解方程 x-x-1=0，取迭代初值 x=1.5, d=10. Matlab程序: function x=snewtonfa(f,x0,d,max) y=diff(f); %取导数 y=inline(y); %定义y f=inline(f); % 定义f disp(k x ); %以指定格式输出k,x. for k=1:max x(k+1)=x(k)-f(x(k))/y(x(k)); % 计算公式 disp(sprintf(%d %f ,k,x(k))); %输出计算结果 if ab