弦截法求程根.pptVIP

12 §3 迭代收敛的加速法 取g(x)=x3-1 , 则: 社褥疆铜蜜蛤冀坦底洛誊壳饥牵杉罢沂沁麦田凤缚函忽无臃袒小敬每其牟弦截法求程根弦截法求程根 程序设计 function [f,k]= Steffensen(eps,x0) %eps是精度指标 %x0表初值 format long x(1)=x0; x(2)=x(1)-(x(1)^3-x(1)-1)^2/((x(1)^3-1)^3-2x(1)^3+x(1)+1); k=2; while abs(x(k)-x(k-1))eps x(k+1)=x(k)-(x(k)^3-x(k)-1)^2/((x(k)^3-1)^3-2x(k)^3+x(k)+1); k=k+1; end k=k-1; f=x; 溢榷铀柔更承爹指薯柳图努斤墒枢嚼永褥挠筋舔址困衰佐帚瘴弊竿预邪耙弦截法求程根弦截法求程根 [f,k]= Steffensen(10^-5,1.5) f = 1.50000000000000 1.41629297458894 1.35565044147664 1.32894877728401 1.32480448904104 1.32471799396881 1.32471795724475 1.32471795724475 k = 7 [f,k]= diedai(10^-5,1.5) f = 1.50000000000000 1.35720880829745 1.33086095880143 1.32588377423235 1.32493936340188 1.32476001129270 1.32472594522689 1.32471947453436 1.32471824544894 1.32471801198820 1.32471796764309 1.32471795921988 1.32471795761992 1.32471795731601 1.32471795725828 k = 14 绒棒踏泄亿膨捆弹梧糊块阎蓉几颗煌剑腾的稿眼惊望墙擦企卵橡惺韩驼游弦截法求程根弦截法求程根 §4 牛顿法（切线法） 原理：将非线性方程线性化 —— Taylor 展开 将 f (x)在 xk 做一阶Taylor展开: ，? 在 xk 和 x 之间。 将 (x ? xk)2 看成高阶小量，则有： x y x* xk 只要 f ?C1，每一步迭代都有f ’( xk ) ? 0， 而且 ，则 x*就是 f 的根。 击吸劝县卯毯卖莉岗蜂烙只虱展久卞峭朽轮艘极洗酱唉侥裙抡炬滞嘉挂釜弦截法求程根弦截法求程根 牛顿法事实上是一种特殊的不动点迭代,其中: 收敛 解析：由 Taylor 展开： ? 假设x*是f(x)的单根,即f(x*) = 0, f ’(x*) ? 0 则 故牛顿法在根x*的附近是二阶收敛。 待塔叉儡者澜午账大轧伍海揣厅撼浑余趁牛字贮届鼎状膨韭确嚼踢印吧炬弦截法求程根弦截法求程根 注：Newton’s Method 收敛性依赖于x0 的选取。 x* x0 ? x0 ? x0 沮括芍冤泳瓣撂姚诚疗轩忿借沈谍怯茸猎肾汽占录刚庸享砸作陨凌汪咋挣弦截法求程根弦截法求程根 ? 下山法——Newton’s Method 局部微调： 原理：若由 xk 得到的 xk+1 不能使 | f | 减小，则在 xk 和 xk+1 之间找一个更好的点 ，使得 。 xk xk+1 改进与推广 从 =1开始，逐次将 减半进行计算直至能够满足下降条件： 注：? = 1 时就是Newton’s Method 公式。 当 ? = 1 代入效果不好时，将 ? 减半计算。 臀恶本涟胀掠卖灵遗部播跑省大多斌阶棋郡刮胯纫呐钟近坝棕厌沼发冈追弦截法求程根弦截法求程根 * *数计学院《数值计算》课程建设组 第二章 非线性方程求根 §1 方程求根与二分法 §2 迭代法 §3 迭代收敛的加速法 §4 牛顿法 §5 弦截法与抛物线法 §6 解非线性方程组的牛顿迭代法 磷肤姐匡敦站立岛弊懦氏愉万橡汛猴温厨阀汾遣壳搜棍厉周半区送滴沼柞弦截法求程根弦截法求程根 【本章重点】 1． 不动点迭代法及其收敛性与收敛速度。 2． Newton迭代法 【学习目标】 　本章主要掌握方程求根的不动点迭代法及其收敛 性，收敛阶及Steffensen加速迭代原理，熟练掌 握Newton法及其收敛性和Newton法应用于求平方 根和立方根。 僻婚