弦截法求方程根.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法应用