第7章 非线性方程和方程组的数值解法(一).ppt

这样导出的迭代公式(5.2)可以看做牛顿公式 中的导数　　 用差商 取代的结果. 　　(5.2)式有明显的几何意义： 设曲线y=f(x)上横坐标为xk-1和xk的点分别为Pk-1和Pk, 则差商 表示弦 的斜率, 弦 的方程为 O x* xk+1 xk Pk xk-1 y x Pk-1 因此，按(5.2)式求得xk+1实际上是两点弦线 与x轴交点的横坐标(令y=0解出x即可).这种算法因此而形象地称为割线 (弦截)法. 　　割线法与牛顿法 (切线法)都是线性化分法，但两者有本质的区别.牛顿法在计算xk+1时只用到前一步的值xk，而割线法要用到前面两步的结果xk-1,xk，因此使用这种方法必须先给出两个开始值x0, x1. 　　定理6 假设f(x)在根x*的邻域内△: |x-x*|≤δ具有二阶连续导数，且对任意x?△有f?(x)≠0，所取的初值x0, x1?△，那么当邻域△充分小时，弦截法(5.2)将按阶 收敛到x*. 这里p是方程λ2-λ-1=0的正根. 　　因为(5.2)式用到前两点xk-1和xk的值，故此方法又称为双点割线法. 每步只用一个新点xk的值，此方法称为单点割线法. 　　如果把(5.2)式中的xk-1改为x0，即迭代公式为 　　例题 用牛顿迭代法和割线法求方程 f(x)=x4+2x2–x–3=0， 在区间(1, 1.5)内之根(误差为10-9). 解 取x0=1.5，用牛顿法, 可得 x6=1.12412303030 取x0=1.5, x1=1，用双点割线法，迭代6次得到同样的结果，而采用单点割线法，则迭代18次得 x18=1.124123029. x y x* xk y=f(x) xk+1 Pk Pk+1 xk+2 y=f(x) x0 x1 x2 x* O x* xk+1 xk Pk xk-1 y x Pk-1 O x* xk+1 xk y=P2(x) xk-2 y x y=f(x) xk-1 y=f(x) 牛顿法 简化牛顿法 割线法 抛物线法 2 抛物线法 　　设已知方程f(x)=0的三个近似根xk,xk-1,xk-2，我们以这三点为节点构造二次插值多项式p2(x)，并适当选取p2(x)的一个零点xk+1作为新的近似根，这样确定的迭代过程称为抛物线法，亦称为密勒(Müller)法. 在几何图形上, 这种方法的基本思想是用抛物线y=p2(x)与x轴的交点xk+1作为所求根x*的近似位置. O x* xk+1 xk y=P2(x) xk-2 y x y=f(x) xk-1 　　抛物线法的几何意义见下面图形. 　　现在推导抛物线法的计算公式. 插值多项式 有两个零点 式中 　　为了在(5.3)式定出一个值xk+1，我们需要讨论根式前正负号的取舍问题. 　　在xk, xk-1, xk-2三个近似值中，自然假定xk更接近所求的根x*，这时，为了保证精度，我们选(5.3)式中接近xk的一个值作为新的近似根xk+1. 为此，只要取根式前的符号与ω的符号相同. 例 用抛物线法求解方程f(x)=xex-1=0. 解 取x0=0.5, x1=0.6, x2=0.56532开始，计算得 f(x0)=-0.175639, f(x1)=0.093271, f(x2)=-0.005031. f[x1,x0]=2.68910, f[x2,x1]=2.83373, f[x2,x1,x0]=2.21418. 故 代入(5.3)式求得 　　事实上, 在一定条件下可以证明, 对于抛物线法，迭代误差有下列渐近关系式 由此式可见抛物线法也是超线性收敛的，其收敛的阶是p=1.840(是方程λ3-λ2-λ-1=0的根)，收敛速度比割线法更接近于牛顿法. 7.6 解非线性方程组的牛顿迭代法 　　考察方程组 　　其中f1,?,fn均为(x1,?,xn)的多元函数. 若用向量记号记x=(x1,?,xn)T∈Rn, F=(f1,?,fn)T, (6.1)就可写成 F(x)=0. (6.2) 当n≥2，且 f1,?,fn 中至少有一个是自变量 x1,?,xn 的非线性函数，则称方程组(6.1)为非线性方程组. 非线性方程组求根问题是前面介绍的方程(即n=2)求根的直接推广，