非线性方程求根--补充知识.ppt

非线性方程的数值解 §1 二分法 §2 迭代法 2.1 迭代格式 2.2收敛性条件 2.3迭代法的收敛阶 §3 牛顿迭代法 3.1 迭代格式 3.2 迭代法的收敛阶 §4 弦割法 解析解与数值解 由推理的方法计算--- 解析解(精确解)。 近似计算（数值计算）---数值解(近似解）。 常借助于计算机。 对于工科而言，各门后续课程和未来的工程实践中遇到最大量的将是数值计算问题。 §1 二分法 §2、迭代解法 一、迭代格式的构造 二、收敛性条件 三、迭代法的收敛阶 §3、牛顿迭代法 一、迭代格式 二、收敛性与收敛阶 §4 弦截法（割线法） 一、双点弦截法 二、单点弦截法 总 结 二 分 法：根的存在唯一性、计算公式及精度控制 一般迭代法：迭代法的构造、收敛性和局部收敛性定 理、收敛阶、精度控制 牛顿迭代法：迭代公式、局部收性与收敛阶、几何意义 弦 截 法：双点弦截公式与单点弦截公式、几何意 义、收敛性与收敛阶 算法的程序： 各个算法的程序设计 得到： 关于收敛阶，由牛顿迭代格式 再由Taylor展式得到： 两式相减得： 即： 当 k?∞ 时 xk ? x* ,ξ? x* 这时 如果 则说明 Newton 迭代法具有二阶收敛速度。 例6 用牛顿迭代法求方程 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 次。 例7 求 的值(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) 不存在时，可以用 作近似代替,得到 并称其为双点弦截法。 其几何几何解释为： 令 y=0, 解得： x0 x y o a b x* x1 x2 (x0 , f(x0)) (x1 , f(x1)) x3 在曲线上任取两点(x0 , f(x0))、 (x1 , f(x1))，作割线，方程为 再过两点(x1 , f(x1))、 (x2 , f(x2))作割线， 得割线方程： 再令 y=0, 解得： 以此类推，最后得到的近似解 x0 , x1 , x2 , … 越来越靠近 x*。 x0 x y o a b x* x1 x2 (x0 , f(x0)) (x1 , f(x1)) x3 在曲线上过两点(x0 , f(x0))、 (x1 , f(x1))作割线，方程为 令 y=0, 解得： 再过两点(x0 , f(x0))、 (x2 , f(x2))作割线 得割线方程： 再令 y=0, 解得： x4 以此类推，可以得到单点弦截公式： * * x y o a b * x* 这种方程往往无法求得其精确解，只能通过数值方法求其近似解。这里我们将介绍两种数值方法： 非线性方程求根是我们经常碰到的问题，例如： (1). 二分法; (2). 迭代法：一般迭代法、牛顿迭代法、弦截法. 对于 f(x)=0