第8章非线性方程求解学案.pptVIP

第八章 非线性方程的数值解法 第八章 非线性方程的数值解法 §1 二分法 §2 迭代法 2.1 迭代格式 2.2收敛性条件 2.3迭代法的收敛阶 §3 牛顿迭代法 3.1 迭代格式 3.2 迭代法的收敛阶 §4 弦割法 §1 二分法 §2、一般迭代解法 一、迭代格式的构造 二、收敛性条件 三、迭代法的收敛阶 §3、牛顿迭代法 一、迭代格式 二、收敛性与收敛阶 §4 弦截法（割线法） 一、双点弦截法 二、单点弦截法 第八章 总 结 二 分 法：根的存在唯一性、计算公式及精度控制 一般迭代法：迭代法的构造、收敛性和局部收敛性定 理、收敛阶、精度控制 牛顿迭代法：迭代公式、局部收性与收敛阶、几何意义 弦 截 法：双点弦截公式与单点弦截公式、几何意 义、收敛性与收敛阶 算法的程序： 各个算法的程序设计 第八章 练习题 令 故 解得 如果 并称（1）为方程求根的牛顿迭代格式。 则 （1） 而 f (x*)=0 方程求根的牛顿迭代法为又称为切线法，可以通过其 几何意义明确这一称呼，如下图所示： x0 x y o a b x* x1 x2 1). 在解的附近任取一点 x0 , 在曲线上过点(x0 , f(x0))作切线 y= f(x0)+f’(x0)(x-x0) (x0 , f(x0)) 令 y=0 ,得到切线与 x 轴的交点 2). 再过点(x1 , f(x1))作切线 y= f(x1)+f’(x1)(x-x1) (x1 , f(x1)) 令 y=0 ,得到切线与 x 轴的交点 以此类推，最后得到的近似解 x0 , x1 , x2 , … 越来越靠近 x*。 对于牛顿迭代格法 迭代函数为 导数为 在精确解 x* 处 , 由 f(x*)=0 得到 | g’ (x*)|=0 =L 1 由局部收敛性知：牛顿法局部收敛（单根附近至少二阶）。 得到： 关于收敛阶，由牛顿迭代格式 再由Taylor展式得到： 两式相减得： 即： 当 k?∞ 时 xk ? x* ,ξ? x* 这时 如果 则 Newton 迭代法具有二阶收敛速度（重根线性收敛） 。 例8.5 用牛顿迭代法求方程 xex-1=0在 x0=0.5 附近的根。 解:已知方程在[0.5 , 0.6]内有唯一实根, 令f(x) = xex-1 则 f ’(x)= (x+1)ex , 这样可以得到Newton迭代格式： 或者 取初值进行计算： x0=0.5 x1=0x2=0x3=0精度为10-5 。如果用一般迭代法 xk+1= e- xk , k=0,1,2, …， 要达到同样精度，则需要迭代 22 次。 x22= 0 例8.6 求 的值(a 0). 例如，对于 解：利用迭顿迭代法，令 得计算格式 则有 xn=a 再令 f(x)=xn-a , 得到 f’(x)=nxn-1, 代入牛顿迭代格式 将 n=2 代入上式得到 ： 如果分别取a=2,a=3,a=5, 并要求精度为ε=10-5，计算结果如下： a =2,ε=10-5 a =3,ε=10-5 a =5,ε=10-5 x0 =1x1 =1x2 =1x3 =1x4 =1x5 =1x0 =1x1 =2x2 =1x3 =1x4 =1x5 =1x0 =1x1 =3x2 =2x3 =2x4 =2x5 =2x6 =2简化牛顿法（平行弦法：线性收敛）；拟牛顿法 对于牛顿迭代法 当 f ’(x) 不存在时，可以用 作近似代替,得到 并称其为双点弦截法。 其几何几何解释为： 在曲线上任取两点(x0 , f(x0))、 (x1 , f(x1))，作割线，方程为 多点迭代法 令 y=0, 解得： 再过两点(x1 , f(x1))、 (x2 , f(x2))作割线，得割线方程： 再令 y=0, 解得： 以此类推，最后得到的近似解 x0 , x1 , x2 , … 越来越靠近 x*。