计算机数值方法第3章.pptVIP

* 如此继续下去，即构成了迭代格式（牛顿迭代格式）： 对于函数方程f(x)=0，且满足下面的条件： （1） f(x)在区间[a,b]上一、二阶导数均存在 且f ’(x)?0，f ”(x)不变号； （2） f(a)f(b)0； （3） 给定初值x0?[a,b]，有 则采用牛顿迭代法所得的迭代序列xk将收敛于函数方程f(x)=0在区间[a,b]上的唯一根?。 * 对以上三个条件的说明： （1） 保证了根的存在； （2） 保证了函数f(x)的单调性及根的唯一性； （3）三个条件则共同保证了序列{xk}的收敛性。 下面只证明在满足上述条件时的收敛性： 为方便叙述，不妨设： 即函数曲线y=f(x)单调增加、上凹。此时，因取初值x0使 则由 得： 从而有： * 故知?x1x0。同理可证：?xkxk-1…x1x0。这表明迭代序列{xk}单调减小而且有下界?，即收敛于所求根?。 [例] 用牛顿迭代法求x3+2x2+10x-20=0在x=1附近的根，计算结果误差小于10-5。 解：令f(x)= x3+2x2+10x-20，考虑区间[0，2]上f(0)f(2)0且f’(x)=3x2+4x+100，f”(x)=6x+40，可知f(x)=0在[0，2]上有唯一的根，按牛顿迭代法构造迭代式： * 取初值x0=2时显然有f(x0)f”(x0)0，故牛顿迭代收敛。 [例] 用牛顿迭代法求 解： 令 则方程f(x)=x2-3=0的根即为所求，取含根区间[1 , 2]时有 所以方程有唯一的根。构造牛顿迭代格式： 取初值x0=2，显然有 故迭代数列收敛于函数方程的根。 * 牛顿迭代法为平方收敛（收敛阶为2）。 证明： 由f(x)在xk处的泰勒展开式 将x=x*代入上式，并由f(x*)=0可得： 故有： 因而： 所以由收敛阶的定义可知，牛顿迭代法具有2阶收敛。 * 弦截法又称割线法。如下图所示，函数方程f(x)=0有根为?，为求出?： 弦 截 法 y=f(x) x1 x0 x2 x3 x4 ? * （1）过点(x0,y0)和点(x1,y1)作曲线y=f(x)的割线，该割线 与x轴交于x2；该割线的点斜式方程为 令y=0，得该割线与x轴的交点x2 : （2）过点(x0,y0)和(x2,y2)再作割线则与x轴交于x3，可得： 如此继续下去，即构成了迭代格式： * 每次所作割线与x轴的交点将逐步逼近所求根?。此种弦割法称为单点弦割法。 双点弦截法：如下图所示，过点(x0,y0)和点(x1,y1)作曲线y=f(x)的割线，该割线与x轴交于x2；过点(x1,y1)和(x2,y2)再作割线则与x轴交于x3，……。 o y=f(x) x y x1 x0 x2 x3 x4 双点弦截法迭代格式为： * 一）数据说明： （1）精度控制量EPS，最大迭代次数MAXREPT； （2）迭代初值x0，x1； （3）x_k0、x_k1、x_k2，用于进行迭代计算； （4）自定义函数f(x) 二）操作说明： step1 输入迭代初值x0和x1； step2 x_k0=x0 ; x_k1=x1 ; x_k2 ; k=0 ; p=1.0; step3 While pEPS And k=MAXREPT Do step4 x_k2=x_k1-f(x_k1)*(x_k1-x_k0)/(f(x_k1)-f(x_k0)); step5 p= fabs(x_k2-x_k1); step6 x_k0=x_k1 ; x_k1=x_k2; step7 k=k+1; step8 输出x_k2即为所求根。 * 一、问题描述 以二阶非线性方程组为例： 非 线 性 方 程 组 求 解 其中至少有一个方程为非线性方程，求其根，即一组数据(x* , y*) 二、基本思想（牛顿迭代法） 将非线性方程组作向量化处理，即： 令 ，则原方程组为 所求根即： ，使 * 如此则可用迭代法求根，即用u0,u1,……,uk,uk+1序列逐次逼近根u*。 三、构造方法 将f1(x , y),f2(x , y)在(xk , yk)处作泰勒展开，并取其线性部分，得到方程组： 令 ，则有 * 将u=uk+1代入上式得牛顿迭代格式： ，即 令 ，得 上式左右两边同乘 ，得 从初值u0开始，逐次计算 ，及 ，直到 * 取初值 四、实例 求解非线性方程组： 解： * 由 按此种方法继续做下去，直到…… * 定义：设f(x)为一元连续函数，称方程f(x)=0为函数方程。特别地，当f(x)不是 x 的线性函数时，称对应的函数方程为非线性方程。 在非线性方程中，当f(x)为多项式函数时，称为代数方程，