第7章-非线性方程和方程组的数值解法.pptxVIP

第七章 非线性方程迭代解法 ;求方程 几何意义;§1 二分法;终止法则？; 二分法算法 给定区间[a,b] ，求f(x)=0 在该区间上的根x. 输入: a和b; 容许误差 TOL; 最大对分次数 Nmax. 输出: 近似根 x. Step 1 Set k = 1; Step 2 Compute x=(a+b)/2; Step 3 While ( k ? Nmax) do steps 4-6 Step 4 If f(x)=0 or |b-a| 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 =(a+b)/2 ;Go To Step 3 ; Step 7 Output the solution of equation: x; STOP.;3、;解：;k;①简单; ② 对f (x) 要求不高(只要连续即可) .;§2 不动点迭代法及其收敛性;;例2：;取;Lipschitz条件成立的充分条件;证明：① φ(x) 在[a, b]上存在不动点？;L 越 收敛越快; 算法: 不动点迭代 给定初始近似值 x0 ，求x = φ(x) 的解. 输入: 初始近似值 x0; 容许误差 TOL; 最大迭代次数 Nmax. 输出: 近似解 x 或失败信息. Step 1 Set i = 1; Step 2 While ( i ? Nmax) do steps 3-6 Step 3 Set x = φ(x0); /* 计算 xi */ Step 4 If | x ? x0 | TOL then Output (x); /*成功*/ STOP; Step 5 Set i ++; Step 6 Set x0 = x ; /* 更新 x0 */ Step 7 Output (The method failed after Nmax iterations); /*不成功 */ STOP.;二、局部收敛性 /*Local Convergence*/;;例3：已知方程 在1.5附近有根，把方程写成三 种不同的等价形式(1) 对应迭代格式 ； (2) 对应迭代格式 ；(3) 对应 迭代格式 ;判断迭代格式在 的收敛性，选 一种收敛格式计算，精确到小数点后第二位。 ;（收敛阶/*the order of Convergence*/);;必要性;k;一、使用两个迭代值的组合方法：;;又如：;二、Steffensen(斯蒂芬森)加速迭代法：（三个迭代值组合）;一般地，;若令 ;例5：;取;;反之，设;;代入上述结果：;;取;将 (x* ? x0)2 看成高阶小量，则有：;无开方运算，又无除法运算。;;;有根;由中值定理， 使得;重复以上过程，可得（归纳法）;注：Newton’s Method 收敛性依赖于x0 的选取。;;;? 求复根 /* Finding Complex Roots */ —— Newton 公式中的自变量可以是复数; §5 弦割法与抛物线法 /* Secant Method and Parabola Method */; 局部收敛性 设 表示区间 ， x*为方程 f (x) =0的根， 函数f (x)在 中有足够阶连续导数， 且 满足 ;Corollary(推论);;? Muller方法的具体实现:; （局部收敛性） 设 ， 在x*的某邻域内连续，则存 在 x* 的一个邻域 ，当 时,由抛物线法产