Ch07非线性方程求根-2010秋.ppt_第1页
Ch07非线性方程求根-2010秋.ppt_第2页
Ch07非线性方程求根-2010秋.ppt_第3页
Ch07非线性方程求根-2010秋.ppt_第4页
Ch07非线性方程求根-2010秋.ppt_第5页
已阅读5页,还剩15页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

1、第7章 非线性方程求根,1根的存在性。方程有没有根?如果有根,有几个根?,2这些根大致在哪里?如何把根隔离开来?,3根的精确化,定理1:设函数 f (x) 在区间a, b上连续,如果 f (a) f (b) 0,则方程 f (x) = 0 在a, b内至少 有一实根x*。,x*,f(x),1画出 f(x) 的略图,从而看出曲线与x 轴交点的位置。,2从左端点x = a出发,按某个预先选定的步长h 一步一步地向右跨,每跨一步都检验每步起点x0 和终点x0 + h的函数值,若,那么所求的根x*必在x0与x0+h之间,这里可取x0 或x0+h作为根的初始近似。,例1:考察方程,x1,x2,a,b,或

2、,不能保证 x 的精度,x*,2,1 二 分 法,执行步骤,1计算f (x)在有解区间a, b端点处的值,f (a),f (b)。,2计算f (x)在区间中点处的值f (x1)。,3判断若f (x1) = 0,则x1即是根,否则检验:,(1)若f (x1)与f (a)异号,则知解位于区间a, x1, b1=x1, a1=a;,(2)若f (x1)与f (a)同号,则知解位于区间x1, b, a1=x1, b1=b。,反复执行步骤2、3,便可得到一系列有根区间: (a, b), (a1, b1), , (ak, bk), ,简单; 对f (x) 要求不高(只要连续即可) .,无法求复根及偶重根

3、收敛慢,定义f (x),f (a) f (b)0,f (a) f (b)=0,f (a) =0,打印b, k,打印a, k,结束,是,是,是,否,否,否,m=(a+b)/2,|a-b|,f(a)f(b)0,打印m, k,a=m,b=m,结束,k=K+1,是,是,否,否,输入,k = 0,例2: 求方程,2 迭 代 法,1简单迭代法,x1 = 0.4771 x2 = 0.3939 x6 = 0.3758 x7 =0.3758,2迭代过程的收敛性,f (x) = 0,x = g (x),迭代格式,定理2:如果 (x) Ca, b满足下列条件 (1)当x a, b时,(x) a, b (2)当任意x

4、 a, b ,存在0 L 1,使 则方程x = (x)在a, b上有唯一的根x*,且对任意初值x0 a, b时,迭代序列xk+1= (xk) (k = 0, 1, )收敛于x*。,(7.1),3迭代法的结束条件,例4:求方程 在0, 0.5内的根, 精确到10-5。,x1 = (0.25) = 0.3385416,(7.2),x2 = (x1) = 0.3462668,x3 = (x2) =0.3471725,x4 = (x3) =0.3472814,x5 = (x4) =0.3472945,x6 = (x5) =0.3472961,x7 = (x6) =0.3472963,4.迭代过程的收敛

5、速度,设由某方法确定的序列xk收敛于方程的根x*, 如果存在正实数p,使得,(C为非零常数),定义:,则称序列xk收敛于x*的收敛速度是p阶的,或称该方法 具有p 阶敛速。当p = 1时,称该方法为线性(一次)收敛; 当p = 2时,称方法为平方(二次)收敛;当1 p 2时, 称方法为超线性收敛。,3 牛 顿 法,1、牛顿法的迭代公式,x*,x0,x1,x2,原理:将非线性方程线性化 Taylor 展开,取 x0 x*,将 f (x)在 x0 做一阶Taylor展开:,, 在 x0 和 x 之间。,将 (x* x0)2 看成高阶小量,则有:,只要 f C1,每一步迭代都有f ( xk ) 0,

6、 而 且 ,,,则 x*就是 f (x)的根。,(7.3),2、牛顿法的收敛性,定理3: 设f (x)在a, b上满足下列条件: (1)f (a) f (b) 0 则由(7.3)确定的牛顿迭代序列xk收敛于f (x) 在a, b上的唯一根x*。,注:Newton法的收敛性依赖于x0 的选取。,x*,1)若 ,则当 时,取x* xk+1,计算过程结束;当 时,则把xk+1作为新的xk值,并重复回到(3)。,3牛顿下山法,牛顿下山法计算步骤可归纳如下:,(1)选取初始近似值x0;,(2)取下山因子 = 1;,(3)计算,(4)计算f (xk+1),并比较 与 的大小,分以下二种情况,2)若 ,则当且,取x* xk,计算

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论