【2017年整理】计算方法大作业.docxVIP

计算方法大作业 学院：电子工程 姓名： 班级： 学号： 大作业选题： 分析方程求根问题中牛顿法的性能，包括收敛性等，并用该方法求解一个问题，给出过程和结果。 一、牛顿迭代法介绍： 用迭代法求方程的根时，首先要构造一个迭代函数，迭代函数构造的好坏，不仅影响收敛速度，而且有可能使迭代序列发散，构造迭代函数的一条重要途径，是用近似方程代替原方程去求根，因此如果能将非线性方程用线性方程来近似代替，那么求近似根问题就容易得到解决，而且十分方便。牛顿法就是把非线性方程线性化的一种方法。 二、牛顿迭代法原理 设已知方程的近似根,则在附近可用一阶泰勒多项式近似代替.因此, 方程可近似地表示为.用表示的根,它与的根差异不大. 设,由于满足解得 重复这一过程,得到迭代公式 这就是著名的牛顿迭代公式,它相应的不动点方程为 . 用牛顿迭代公式求方程根的方法称为牛顿迭代法，简称牛顿法。 三、牛顿迭代法的几何解析 在处作曲线的切线，切线方程为。令，可得切线与轴的交点坐标，这就是牛顿法的迭代公式。因此，牛顿法又称“切线法”，其几何意义即为点处的切线方程。 四、牛顿迭代法的收敛性 计算可得,设是的单根,有,则 , 故在附近,有.根据不动点原理知牛顿迭代法对单根收敛.同理可知当是的重根时也收敛，则可分析出牛顿法不论对单根还是重根均是局部收敛的，只要初值足够靠近，牛顿迭代序列均收敛于。 五、牛顿迭代法的收敛速度 定理(牛顿法收敛定理) 设在区间上有二阶连续导数,且满,在上不变号, 在上不等于0,令 有.则对任意,牛顿迭代格式收敛于在中的唯一实根,并且: (1) (2) (3) ,牛顿迭代法为2阶收敛. 六、牛顿迭代法的优缺点及比较： 迭代法是求解方程近似根的一个重要方法，也是计算方法中的一种基本方法，它的算法简单，适用于求方程或方程组近似根的一种常用的算法设计方法。牛顿迭代法最大的优点是在方程的单根附近具有平方收敛，而且该法还可以用来求方程的重根、复根。牛顿法方法简单，每次迭代都是简单地重复运算，易于编制程序，与求解线性方程的精确法相比，迭代法对于字长位数较少的计算机更为适用，他可以用增加迭代次数来弥补字长位数少的不足，初值可以任意取，因而中间结果的偶然错误不影响最后结果的获得，但迭代速度稍慢 七、牛顿迭代求根的方法： 设方程为，用某种数学方法导出等价的形式x=g(x),按下面步骤执行： 选一个方程的近似根，赋给变量x0； 将x0的值保存于变量x1，然后计算g(x1)，并将结果存于变量x0； 当x1与x0的差的绝对值还小于指定的精度要求时，重复步骤（2）的计算。若方程有根，并且用上述方法计算出的近似根序列收敛，则按上述方式求得的x0就认为是方程的根。 八、用牛顿迭代法求解一个问题： 题：用牛顿法求下列方程的根：x*e^x-1=0; 相应的MATLAB代码为： x0=1.0; N=100; k=0; eps=5e-6; delta=1e-6; while(1) x1=x0-fc1(x0)/fc2(x0); k=k+1; if kN disp(New method failed) break end if(abs(x1-x0)delta || abs(fc1(x1))delta) break; end x0=x1; end fprintf(%f,x0) fprintf( %f , abs(fc1(x1)) ) %原函数 function fc1=fc1(x) fc1=x*exp(x)-1; end %导数 function fc2=fc2(x) fc2=(x+1)*exp(x); end 相应的结果： 原方程的根为：x=0.5671435302 牛顿迭代法由于设置delta=1e-6，所以算出的误差e1.0*10^-5; 注意：牛顿迭代方法能够有效的基本条件是：迭代公式必须是收敛的。