数值计算方法课件 第二章 非线性方程的数值解法VIP

反之，设型洛比塔法则其次证明Steffensen 迭代是二阶收敛的由Taylor 公式：在处展开则上式中用替换即证代入上述结果：注意：对本来就是P(1) 阶收敛的方法，改用Stefensen 迭代方法优点不多。取计算结果如下：法2 原迭代次数29 法3 原来不收敛法1 原来不收敛返回§4 牛顿法/* Newton - Raphson Method */ 一、牛顿迭代公式的推导1、待定参数法不动点迭代的关键是构造满足收敛条件的迭代函数一种自然的选择是令为了加速不动点迭代的收敛过程，应尽可能使迭代函数在处有更多阶导数等于零（定理2.5 ）。令现设取满足因此，选取迭代函数Newton –Raphson 迭代格式称之为牛顿—拉夫森方法，简称牛顿法原理：将非线性方程线性化取x0 x* ，将f (x) 在x0 做一阶Taylor 展开: ，在x0 和x 之间2、Taylor 展开法/* Taylor’s expansion Method */ 将(x* x0)2 看成高阶小量，则有：x y x* x0 只要f C1 ，每一步迭代都有而且，则x* 就是f 的根。与x轴交点的横坐标无开方运算，又无除法运算。例1：写出求的Newton 迭代格式；写出求的Newton 迭代格式,要求公式中既解：等价于求方程的正根解法一：等价于求方程的正根解法二：等价于求方程的正根Th2.7 （局部收敛性）设x* 为方程f (x) =0 的根，在包含x* 的某个开区间内连续，且，则存在x* 的邻域，使得任取初值，由Newton’s Method 产生的序列以不低于二阶的收敛速度收敛于x* ，且证明：Newton’s Method 事实上是一种特殊的不动点迭代其中，则收敛由Taylor 展开：在单根/*simple root */ 附近收敛快只要，则令可得结论。Th2.5 有根根唯一产生的序列单调有界保证收敛证明：因为f C2[a, b] ，由（1）和（2）知f (x) 在[a, b] 内有唯一根下面由条件（1）、（2）分4种情况讨论：仅证明第一种情况，其它情况类似讨论Th2.8 （收敛的充分条件）设f (x) =0 且f C2[a, b] ，若(1) f (a) f (b) 0 ；(2) 在整个[a, b] 上不变号且；(3) 选取x0 [a, b] 使得；则Newton’s Method 产生的序列{ xk } 收敛于方程的根, 且由中值定理，使得因此即在上单调递增由另一方面，由Taylor 展开得介于、之间* 第二章非线性方程的数值解法/* Numerical Solutions of Nonlinear Equations*/ 本章主要内容：1 、二分法2 、不动点迭代的构造及其收敛性判定（重点）3 、Newton 和Steffensen 迭代（重点）4 、割线法5 、非线性方程组的迭代解法历史背景代数方程的求根问题是一个古老的数学问题。理论上，次代数方程在复数域内一定有个根(考虑重数)。早在16 世纪就找到了三次、四次方程的求根公式，但直到19 世纪才证明大于等于5次的一般代数方程式不能用代数公式求解，而对于超越方程就复杂的多，如果有解，其解可能是一个或几个，也可能是无穷多个。一般也不存在根的解析表达式。因此需要研究数值方法求得满足一定精度要求的根的近似解。求方程几何意义基本定理如果函数在上连续，且则至少有一个数使得，若同时的一阶导数在内存在且保持定号，即( 或) 则这样的在内唯一。a b x* §1 二分法/* Bisection Method */ 原理：若f C[a, b] ，且f (a) · f (b) 0 ，则f 在(a, b) 上至少有一实根。基本思想：逐步将区间分半，通过判别区间端点函数值的符号，进一步搜索有根区间，将有根区间缩小到充分小，从而求出满足给定精度的根的近似值。以此类推终止法则？a b x1 x2 a b When to stop? 或不能保证x 的精度x* 2 x x* 二分法算法给定区间[a,b] ，求f(x)=0 在该区间上的根x. 输入: a 和b; 容许误差TOL; 最大对分次数Nmax. 输出: 近似根x. Step 1 Set k = 1; Step 2 Compute x=f((a+b)/2); Step 3 While ( k Nmax) do steps 4-6 Step 4 If |x| TOL , STOP; Output the solution x. Step 5 If x*f(a)0 , Set b=x; Else Set a=x; Step 6 Set k=k+1; Compute x=f((a+b)/2);Go To Step 3 ; Step 7 Output the soluti