不动点迭代法.docxVIP

中国矿业高校(北京)理学院 数值分析试验报告 试验名称 不动点迭代法求方程的近似根 试验时间 201年3月 20日 组长签名 龙纯鹏 班级 信息与计算 科学(1)班 学号 成果 组员签名 一、试验目的，内容二、相关背景学问介绍三、代码 四、数值结果五、计算结果的分析六、计算中消失的问 题，解决方法及体会 一、试验目的，内容 5、先确定方程x =的一个收敛的有根区间切，然后用不动点迭代法求* = ￡;在此 有根区间的近似根，初值与自己确定，要求根精确到IO并求迭代次数. ,把握二分法，会编写相关代码。 二、相关背景学问介绍 (1)算法原理或计算公式 二分法的基本思路是通过计算隔根区间的中点，逐步将隔根区间缩小，从而可得方程 的近似根数列上｝。％* -々 ￡ (2)程序设计思路 二分法原理： 不妨设方程/(X)= 0在隔根区间［a,b］上f(a) 0,/S) 0.. ⑴将区间出,切二分得中点岁，计算/(%)在中点的函数值/(彳)，若/(4) = 0,则 X*=9心就是方程的根； 2 否则，若/(岁)0,由于/(%)在内不变号，所以隔根区间变为区芋川.同理，若 八4)0,隔根区间变为彳］.将新区间记为［%闵. xO =0.1673 x1 =0.1707 k = 147 xO =0.1707 x1 =0.1674 k = 148 xO =0.1674 x1 = x1 =0.1676 k = xO =0.1676 x1 = k = 158 xO =0.1678 x1 =0.1702 k = 159 xO = 0.1702 x1 =0.1679 k=160xO= 0.1691 x1 =0.1689 k= 250xO= 0.1689 x1 =0.1691 k= 251xO= 0.1691 x1 =0.1689 k= 252xO= 0.1689 x1 =0.1691 k= 253xO= 0.1691 x1 =0.1689 k= 254 xO= 0.1689 x1 =0.1691 k= 255 x0=0.1691 x1 =0.1689 k= 256 x0=0.1689 x1 =0.1691 k= 257 x0=0.1691 x1 =0.1689 k= 258 x0=0.1689 x1 =0.1691 k= 259 x0=0.1691 x1 =0.1689 k= 260 x0=0.1689 x1 =0.1691 k= 261 x0=0.1691 x1 =0.1690 k= 262 x0=0.1690 x1 =0.1691四、数值结果 k 1 2 3 4 5 6 7 8 9 X 0.5000 0.7500 0.6250 0.6875 0.6563 0.6406 0.6484 0.6445 0.6426 k 10 11 12 13 14 15 16 17 X 0.6416 0.6411 0.6414 0.6412 0.6412 0.6412 0.6412 0.6412 f= -5.4870e-006 k 二17 五、计算结果的分析 这个方程运用二分法求得四位有效数字的解为x=0.6414o 图一 六、计算中消失的问题，解决方法及体会.刚开头设计程序，消失死循环，无法得出结果，要经过不断调试。 .不熟识Matlab的用法，导致解决问题效率下降，要学好相关基础学问。 3?由图一可以看出二分在区间很大时，收敛接近根的速度很快，而当区间较小时，靠近要求 精度根速度变得相当缓慢，计算量也很大，对满意高精度的要求比较困难。我们想象此法可 以和别的后半段接近根较快的方法结合，先求出适用于别的方法的初值，应会优化算法。 指导老师: 2013 年3 月12 日 (2)将时4］重复上述步骤⑴.得到一系列隔根区间: [a,……刀卜攵力/^. 并有f（%）*f（a）＜o,x* c（以也），且后一区间的长度都是前一区间长度的一半，所以的长度为a-%（左foo），区间上,幻的长度趋向于零，即这些区间最终收缩于一点 乙 一明显￡就是方程/（x） = o的根。 三、代码(Matlab) clearxO = l e=10A (-5)k=l xl=sqrt(0.2*eAxO)while (abs (xO-xl)e) k=k+l xO=xl xl=sqrt(0.2*eAxO) end xOx0=0.1753 x1 =0.1630k= 100 x0=0.1630x1 =0.1750 k= 101x0=0.1750 x1 =0.1633k= 102 x0=0.1633x1 =0.1747 k= 103xO = 0.1747x1 =0.1636 k= 104x0=0.1636 x1 =0,1744k= 105 x0=0.1744x1 =0.1639 k= 106x0=0.1639 x1 =