资源描述：
《数值分析3牛顿迭代法》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。
1、§3牛颉連代法NewtonIteration切线法牛顿迭代法是最著名的方程求根方法。已经通过各种方式把它推广到解其他更为困难的非线性问题。【例如】非线性方程组、非线性积分方程和非线性微分方程。虽然牛顿法对于给定的问题不一定总是最好的方法，但它的简单形式和快的收敛速度常常使得解非线性问题的人优先考虑它。迭代一般理论告诉我们，构造好的迭代函数可使收敛速度提高。然而迭代函数的构造方法又各不相同，方法多样。牛顿法是受几何直观启发，给出构造迭代函数的一条重要途径。牛顿迭代的基本思想：方程f(x)=0的根，几何意义是曲线y=f(x)与ox轴y=0的交点。求曲线与y=0的交点没有普遍的公
2、戎，但直接与Ox轴的交点容易计算。用直线近似曲线y=f(x)，从而用直线方程的根逐步代替f(x)=0的根。即把非线性方程逐步线性化。方法：设xk是f(x)=0的一个近似根，把f(x)在xk处作一阶Taylor展开，得到.,(又)=/(%々)+fxk)(x-xj(19)设/'(X)判,由于+fxkx-Xk)^f(x)=0所以求得解记为+
3、，有牛顿迭代公式:按牛顿迭代计算称为牛顿迭代法。牛顿法的几何意义：选初值Xk以后，过/XA，/(A))点，作(21)则曲线y=f(x)的切线，其切线方程为/(X)=f(xk)+fxk){x-xk)^+
4、=-/⑷//'(々)切线与ox
5、轴的交点，hxk(22)牛顿迭代法也称为切线法迭代法的收敛性:如果取=X=g(x)，从而牛顿迭代公式就是人+1=
6、若f(x)在根附近存在连续的二阶导数,x*是f(x)的单根，且初始值X。充分接
7、近x*，则牛顿迭代过程收敛，而且有卜1-x^fx^l2fx^xk-Z
8、2(25)证明1)对于f(x),取g(x)=x-/(x)//(x),则牛顿迭代过程为注意到gx)=f(x)fx)/[fx)]2;gx)=/(x)r(x)+/(x)，(x)//(x)2-2，(x)r(ww由于X*是f(x)=O的单根，即/(x*)=0,/'(Z)弇0,所以有/(Z)=0,=fx)//(Z)矣0(26)由定理2知，迭代过程是局部收敛的。2)将g(X)在X*处进行泰勒展开并代入•¥=',有+-x*)+gg)(〜•-%*)2=g(x*)+1•[/"(X*)//’(x*)](xA
9、,-X*)2注意到+
10、=g(X*),X*=g(X*),得到A^+,-X*=(,(/
11、)/2/(Z)).(x々-Z)2因此有
12、^+1-Zp
13、定理证毕【说明】q+1/eA244c=
14、/7%*)/2/z(x*)p0)，即牛顿迭代过程在x*附近具有平方收敛速度。牛顿迭代法的优点:快速收敛性，算法简单、容易实现缺点：初值&必须选在X*附近，否则，可能不收敛【例】用牛顿法解下面方程在x=0.5附近的根，要求精确到g二1(T5。解牛顿迭代公式为從A-1=0A+l=心_X,=0.5710204,又2=0.5671555,&=0.5671432,又4=0.5671432【注】牛顿法的收敛速
15、度非常快附：牛顿法的算法函数程序：(存至work目录中)functiony=newton(x0)xl=x0-fl(xO);n=l；while(norm(xl-xO)>=1.Oe-6)&(n<=1000)x0=xl;xl=x0-f1(xO);n=n+l;endvpa(xl,7),n%输出n,方程的近似解functiony=f1(x)y=(x-exp(-x))/(1+x);endnewton(0.5)观察初值的选取对解的影响【例】用牛顿法求方程/(幻=?一*-1=0的f首先选取初值：(转到指令窗)x=-10:0.01:10;y=x.A3-x-l;plot(x,y,▼r1,x,0
16、*x)二xkxl~xk~lfunctiony=f1(x)y=(xA3_x_l)/(3*xA2-l);end分别取初值=0.6和1.3,比较所用迭代次数.newton(0.6)newton(1.3)
