数值分析---第三章 非线性方程数值解.pptVIP

非线性方程数值解 二分法 例题 迭代法 2.收敛定理 迭代法的收敛阶 二、迭代加速公式 求解非线性方程组的迭代法 例题 2.Newton迭代法的几何意义 与 y =0求交点，解出 x ，并令 x =xn+1 , 则 过点(xn，f(xn))的切线方程为 x x0 x1 x2 x* y O y=f(x) Newton切线法 3. Newton迭代法收敛定理 （1） 连续，且分别不变号； 则 Newton 迭代法产生的数列{xn},收敛到根x*;并且当x*是根时， Newton 迭代格式是二阶收敛的。 在 有根x*，且 在 设 ,使 （2） 取初值 为例证明（其它情况类似） 证： 以 将f(x*)在xn处Taylor展开 说明数列{xn}有下界x*, 又 故{xn}单调递减。 所以{xn}收敛。设 于是， ， 所以，此时Newton法至少二阶收敛. Newton迭代序列收敛于方程的根 例 解：设 取 ，则由（2.1） 用 Newton 迭代法求 Newton 迭代法的特点是： 1.对初值 x0 的选取要求较高。一般的， Newton 迭代法只有局部收敛性，当初值 x0 在收敛区间里时，收敛速度很快(平方收敛)。但初值 x0离方程根x*较远时，不能保证Newton 迭代法收敛。 x0 x1 x* 2. Newton 迭代法求单根时，收敛速度很快(平方收敛)。但如果方程根x*是重根，则收敛速度较慢，且重数越高速度越慢。 则至少能保持平方收敛。 但可以证明：当x*是m重根时,用下面的迭代格式： 3. Newton下山法： 由于当初值 x0离方程根x*较远时，不能保证Newton 迭代法收敛，但一旦xk进入收敛区间，则收敛速度很快。为使xk尽快进入收敛区间，常采用Newto下山法： ?称为下山因子 具体做法如下： 2.取下山因子?=1 ； 1.选取初值 x0； 3.计算xk+1: 4.计算f(xk+1)并比较| f(xk+1) |与| f(xk) |的大小： (1)若| f(xk+1) || f(xk) |, 则 I)当| xk+1 ? xk |?时，取x*? xk+1 ，结束； 2)若| f(xk+1) || f(xk) |, 则 取?=?/2，返回3； II)当| xk+1 ? xk |??时，将xk+1作为新的xk值继续计算； 开始 输入x0,e i=1 f ? (x0)=0? y0=f(x0) y’0=f ’(x0) u=1 x1=x0-u*y0/y’0 y1=f(x1) |y1||y0| ? |x1-x0|e ? 输出奇异标志 Y Y N N u=u/2 i=i+1 x0=x1 N Y 输出x1 结束 Newto下山法框图 将Newton切线法中的切线斜率f ?(x)用弦的斜率替换： x x0 x1 x2 x* y O y=f(x) 已知 弦截法: 一.简单迭代法 求解方程组： 将其写成矩阵方程形式：F(X)=0 , 其中 * 理学院 University of Shanghai for Science and Technology College of Science 上海理工大学 基本概念 方程的根：使方程 f(x)＝0 成立的x*称为方程的根，也叫做函数f(x)的零点。 如果f(x)为多项式函数，则称方程f(x)＝0为代数方程；否则称为超越方程。 如果f(x)可分解为f(x)=(x- x*)m g(x)，其中g(x*)?0，则x*称为f(x)=0的m重根。 m＝1时x*称为单根。 设函数f(x)在区间[a , b]上连续、单调，并且 f(a) f(b)0,则方程f(x)＝0在(a , b)内有且只有一个解。 二分法的基本思想是将含有根的区间二分，通过判断函数值的符号确定根所在的新区间，逐步将含有根的区间缩小，从而获得满足精度要求的根的近似值。 具体做法：设f(a) f(b) 0，则 存在x* ?(a , b)为方程的根 ， 1）取x0=(a + b)/2，并判定f(a) f(x0) 的符号， 若f(a) f(x0) 0，取a1 = x0 ， b1 = b; 若f(a) f(x0) 0，取a1 =a ， b1 = x0 于是x* ?(a1 , b1) ; 2）再取x1=(a1 + b1)/2，并判定f(a1) f(x1) 的符号， 若f(a1) f(x1) 0，取a2 = x1 ， b2 = b1; 若f(a1)