非线性方程求根试验.doc

非线性方程求根实验 学号：P0725005 姓名：乔海亮 TOC \h \z \t 数值实验 1级,1,数值实验 2级,2 1. 问题描述 2 2. 问题求解 2 2.1. 不动点迭代法 2 2.2. Steffensen迭代法 2 2.3. Newton法 3 3. 数值结果 4 3.1. 4 3.2. 4 3.3. Steffensen迭代法，结合（１） 4 3.4. Steffensen迭代法，结合（２） 5 3.5. Newton法 5 4. 结论 5 5. 代码和程序说明 5 5.1. 开发环境 5 5.2. 主要变量和函数说明 5 5.3. 程序说明 6 问题描述 求方程在附近的根（精确解是），精度为，使用下列方法求解 （１） （２） （３）Steffensen迭代法，结合（１） （４）Steffensen迭代法，结合（２） （５）Newton法 问题求解 不动点迭代法 设一元函数f(x)是连续的，为了求一元非线性方程 的实根，先将它转化成等价形式 其中是一个连续函数。然后构造迭代公式 对于给定的初始值，若由此迭代生成的序列有极限 则有 从而，也是方程的根。 Steffensen迭代法 对于不动点迭代法，当k充分大时，有 即 已知，则，，上式改写成 这样的迭代法称为Steffensen迭代法。 Newton法 设x *是方程f (x ) = 0的根，又x0 为x * 附近的一个值 ，将f (x ) 在x0附近做泰勒展开 其中在和之间，令，则 去掉的二次项，有： 即 以x1代替x0重复以上的过程，继续下去得： 以此产生的序列{Xn}得到f(x)=0的近似解，称为Newton法，又叫切线法。 数值结果 n=0 x= 1.000000000000 n=1 x= 1.538461538462 n=2 x= 1.295019157088 n=3 x= 1.401825309449 n=4 x= 1.354209390404 n=5 x= 1.375298092487 n=6 x= 1.365929788171 n=7 x= 1.370086003402 n=8 x= 1.368241023613 n=9 x= 1.369059812007 n=10 x= 1.368696397556 n=11 x= 1.368857688629 n=12 x= 1.368786102578 n=13 x= 1.368817874396 n=14 x= 1.368803773144 n=15 x= 1.368810031675 n=16 x= 1.368807253961 n=17 x= 1.368808486789 n=18 x= 1.368807939625 n=19 x= 1.368808182472 n=20 x= 1.368808074689 n=21 x= 1.368808122526 n=22 x= 1.368808101295 n=23 x= 1.368808110718 发散，无解。 Steffensen迭代法，结合（１） n=0 x= 1.000000000000 n=1 x= 1.370813882687 n=2 x= 1.368808169945 n=3 x= 1.368808107821 n=4 x= 1.368808107821 Steffensen迭代法，结合（２） n=0 x= 1.000000000000 n=1 x= 1.333492139114 n=2 x= 1.368415439108 n=3 x= 1.368808058309 n=4 x= 1.368808107821 n=5 x= 1.368808107821 Newton法 n=0 x= 1.000000000000 n=1 x= 1.437500000000 n=2 x= 1.363150785855 n=3 x= 1.369366687777 n=4 x= 1.368753732168 n=5 x= 1.368813408528 n=6 x= 1.368807591163 n=7 x= 1.368808158181 n=8 x= 1.368808102913 n=9 x= 1.368808108300 结论 不动点迭代法是线性收敛的，收敛速度（迭代23次）较慢，且对的选取有限制。 Steffensen迭代至少是平方以上收敛，（1）仅仅4次迭代就达到精度要求。本来发散的（2），在此也仅用5次迭代就达到精度要求。 Newton迭代法最实用，它在单根附近具有2阶收敛。但没有Steffensen迭代快。 代码和程序说明 开发环