《数值分析3牛顿迭代法》由会员分享，可在线阅读，更多相关《数值分析3牛顿迭代法（5页珍藏版）》请在装配图网上搜索。
1、3牛顿迭代法NewtonIteration切线法牛顿迭代法是最著名的方程求根方法。已经通过各种方式把它推广到解其他更为困难的非线性问题。【例如】非线性方程组、非线性积分方程和非线性微分方程。虽然牛顿法对于给定的问题不一定总是最好的方法，但它的简单形式和快的收敛速度常常使得解非线性问题的人优先考虑它。迭代一般理论告诉我们，构造好的迭代函数可使收敛速度提高。然而迭代函数的构造方法又各不相同，方法多样。牛顿法是受几何直观启发，给出构造迭代函数的一条重要途径。牛顿迭代的基本思想：方程f(x)=O的根，几何意义是曲线y=f(x)与ox轴y=0的交点。求曲线与y=0的交点没有普遍的公式，但直接与0x轴的交
2、点容易计算。用直线近似曲线y=f(x)，从而用直线方程的根逐步代替f(x)=0的根。即把非线性方程逐步线性化。方法：设x是f(x)=0的一个近似根，把f(x)在乂处作一阶Taylorkk展开，得到f(x)沁f(x),f(x)(xx)(19)kkk设广(x)MO，由于kf(x)，f(x)(xx)沁f(x)0kkk所以求得解记为x，有k,120)f(x)牛顿迭代公式：x,x-仁k+1kf(X,按牛顿迭代计算称为牛顿迭代法。牛顿法的几何意义：选初值X以后，过P(x,/(x)点，作kkk曲线y=f(x)的切线，其切线方程为(21)则20)20)迭代法的收敛性：如果取g(X)，X-f(X)/f(X)，则
3、有kkX=g(x)，从而牛顿迭代公式就是Xk+1因此就可以由考察g(x)的性质，来讨论迭代法的收敛性及收敛速度。迭代过程的收敛速度是指迭代过程中误差的下降速度。【收敛阶定义设迭代过程Xk+1=g(Xk)收敛于方程x=g(x)的根x*，如果迭代误差ek，Xk-X*，当k时成立e/epc(c0的常数)k+1k24)20)则称该迭代过程是p阶收敛的。特别p=1时称为线性收敛，p1时称超线性收敛，p=2时称为平方收敛。【定理3】|若f(x)在根附近存在连续的二阶导数，x*是f(x)的单根，且初始值x0充分接近x*，则牛顿迭代过程收敛，而且有X-X*k+1f(X*)/2f(X*)-25)证明1)对于f(
4、X),取g(x)，x-f(x)/f(x),则牛顿迭代过程为x，g(x)，注意到k+1kg(x)，f(x)f(x)/f(x)2；g(x)，f(x)f(x)+f(x)f(x)/f(x)2-2f(x)f(x)2/f(x)3由于x*是f(x)=0的单根，即f(x)=0,f(x)0,所以有g(x)，0,g(x)，f(x)/f(x)0(26)由定理2知，迭代过程是局部收敛的。2)将g(x)在x*处进行泰勒展开并代入x，x，有g(x*)2!(xx*)2kkg(x)沁g(x*)+g(x*)(xx*)+kk，g(x*)+-f(x*)/f(x*)(xx*)22k注意到x，g(x)，x*=g(x*)，得到k+1kx
5、k+1x)2x(f(x)/2f(x)-(xk因此有xxk+1|f(x)/2f(x)-|xx定理证毕【说明】e/e2tc(c，|f(x*)/2f(x*)|0)，即牛顿迭代过程在k+1kX*附近具有平方收敛速度。牛顿迭代法的优点:快速收敛性,算法简单、容易实现缺点：初值x0必须选在X*附近，否则，可能不收敛【例】用牛顿法解下面方程在x=0.5附近的根，要求精确到S，10-5。xex一10解牛顿迭代公式为x一e-xkxx一k+1k1,xk取初值x0.5，迭代计算，得到0x0.5710204,x0.5671555,x0.5671432,x0.56714321234【注】牛顿法的收敛速度非常快附：牛顿法
6、的算法函数程序：(存至work目录中)functiony=newton(x0)x1=x0-f1(x0);n=1;while(norm(x1-x0)=1.0e-6)&(n=1000)x0=x1;x1=x0-f1(x0);n=n+1;endvpa(x1,7),n%输出n,方程的近似解functiony=f1(x)y=(x-exp(-x)/(1+x);endnewton(0.5)观察初值的选取对解的影响【例】用牛顿法求方程f(x)x3x1=0的根。首先选取初值：(转到指令窗)x=-10:0.01:10;y=x.A3-x-1;plot(x,y,r,x,0*x)解牛顿迭代公式为x3x1X二x-kkk+1k3x21kfunctiony=f1(x)y=(xA3-x-1)/(3*xA2-1);end分别取初值x0二0.6和1.3,比较所用迭代次数.newton(0.6)newton(1.3)
数值分析3牛顿迭代法
