计算方法之分线性方法求根.ppt

西华大学计算机科学与工程系 定义：设f(x)为一元连续函数，称方程f(x)=0为函数方程。特别地，当f(x)不是 x 的线性函数时，称对应的函数方程为非线性方程。 零点定理指出： 若f(x)?[a , b]，且满足f(a)f(b)0，则在区间[a , b]上至少有一点?，使f(?)=0。 f(c9)= -0.01605，含根区间为 [a10 , b10]=[ 1.364257813 , 1.365234375] INPUT 端点a , b；误差EPS；迭代次数上界N；定义函数f(x)。 OUTPUT 近似解 c 或迭代失败的信息 Step1 FA= f(a) ; FB= f(b) If（FA?FB＜0）Then k=0 ; Do Step2 Else OUTPUT （‘不能用对分法求解’) STOP End If Step2 While k＜N Do Steps3~6 Step3 c=(a+b)/2 ; FC= f(c) Step4 If（|FC|＜EPS＝Or（(b-a)/2EPS）Then OUTPUT（c , k） //迭代成功 STOP End If Step5 k=k+1 Step6 If（FA?FC0）Then a=c ; FA=FC Else b=c End If //计算ak，bk转向Step3 Step7 OUTPUT（‘次数已达上界，迭代不收敛’ , k） STOP [练习题] 求方程cosx=x的根（1）验证方程在区间[0,?/2]内有单根；（2）令a0=0，b0=?/2，用对分法求解可得一系列的含根区间[a0 , b0]、[a1 , b1]、[a2 , b2]……，试求出a2、b2。 1、由非线性方程f(x)=0，将其变换为等价形式 x=?(x)； 2、构造迭代公式 xk+1=?(xk)； 3、取初值 x0，计算x1=?(x0)、x2=?(x1)、…，即可获得一个数列{xk}，并以此数列的收敛点作为所求根的数值计算方法—简单迭代法，其中?(x)——迭代函数， {xk}—迭代数列，如果该数列收敛，即 ： x7=2.094552215, x8=2.094551593, x9=2.094551498, x10=2.094551484, x11=2.094551482=x12 由于x11=x12，再迭代已无变化，可见? ? x11。 说明：方程f(x)=0 改写为等价形式x= ?(x)，则求f(x)=0的根即求直线y=x与曲线y= ?(x)的交点。由xk+1= ?(xk)的迭代过程示意如下图： II） | ?’(x) |1 1）取初值x0 2）x1= ?(x0) 3）x2= ?(x1)，…… 迭代计算的结果xk+1离直线与曲线的交点(x*, ?(x*)越来越远。即迭代序列不能收敛。 在区间[1 , 2]上显然有|?’(x)|1，而且当x?[1 , 2]时显然也有?(x)?[1 , 2]，由迭代法收敛定理可知，迭代格式 由迭代法收敛定理可知，在区间[1，3]上，迭代式 2、不同迭代法的收敛速度也存在较大差异，为此引入收敛阶的概念。 可知，当迭代过程收敛且?’(x)连续时，有 这一操作总结为如下式表示的过程： （2）过点(x1,y1)再作切线则与x轴交于x2，同理可得： 如此继续下去，即构成了迭代格式（牛顿迭代格式）： 对以上三个条件的说明： （1） 保证了根的存在； 故知?x1x0。同理可证：?xkxk-1…x1x0。这表明迭代序列{xk}单调减小而且有下界?，即收敛于所求根?。 取初值x0=2时显然有f(x0)f”(x0)0，故牛顿迭代收敛。 牛顿迭代法为平方收敛（收敛阶为2）。 证明： 由f(x)在xk处的泰勒展开式 弦截法又称割线法。如下图所示，函数方程f(x)=0有根为?，为求出?： （1）过点(x0,y0)和点(x1,y1)作曲线y=f(x)的割线，该割线 与x轴交于x2；该割线的点斜式方程为 每次所作割线与x轴的交点将逐步逼近所求根?。此种弦割法称为单点弦割法。 一）数据说明： （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 ;