- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
非方程的数值解法选编
教学目的
1. 掌握解非线性方程(组)的二分法和插值法;
2. 掌握解非线性方程(组)的一般迭代法及有关收敛性的证明与牛顿法;
3. 掌握解非线性方程(组)的牛顿法
4. 了解加速收敛的方法。
教学重点及难点
重点是解非线性方程(组)的牛顿法;
难点是迭代法的收敛性的证明。
; 代数方程求根问题是一个古老的数学问题.早在16世纪就找到了三次,四次方程的求根公式.但直到19世纪才证明了n≥5次的一般代数方程式不能用代数公式求解.因此需要研究用数值方法求得满足一定精度的代数方程式的近似解. ; 求解非线性方程的根,就是求解高次方程或超越方程(含有指数和对数等),因为这类方程没有固定的求根公式。
用 f(x)表示方程左端的函数,则一般的非线性方程可表示为 f (x) = 0.
本节的任务就是上述方程的根或函数的零点。; 通常方程根的数值解法大致分为三个步骤进行
①?判定根的存在性。即方程有没有根?如果有根,有几个根?
② 确定根的分布范围。即将每一个根用区间隔离开来,这个过程实际上是获得方程各根的初始近似值。
③ 根的精确化。将根的初始近似值按某种方法逐步精确化,直到满足预先要求的精度为止. ; 由高等数学知识知, 设f (x)为区间[a,b]上的单值连续, 如果f (a)·f (b)0 , 则[a,b]中至少有一个实根。如果f (x)在[a,b]上还是单调地递增或递减,则仅有一个实根。;(1) 画图法;画图法;对于某些看不清根的函数,可以扩大一下曲线;y;例1 方程f(x)=x3-x-1=0 确定其有根区间
解:用试凑的方法,不难发现 f(0)0 f(2)0
在区间(0,2)内至少有一个实根
设从x=0出发,取h=0.5为步长向右进行根的搜索,列表如下; 用逐步搜索法进行实根隔离的关键是选取步长h要选择适当h ,使之既能把根隔离开来,工作量又不太大。
为获取指定精度要求的初值,可在以上隔离根的基础上采用对分法继续缩小该含根子区间
下面的二分法可以看作是搜索法的一种改进。;(3) 方程求根的二分法(对分法或分半法) (bisection method);例;实根,要求准确到小数点后的第2位。;1.3203
1.3281
1.3242
-; 上述二分法的优点是算法简单,而且在有限区间内,收敛性总能得到保证.值得注意的是,为了求出足够精确的近似解,往往需要计算很多次函数值,是一种收敛较慢的方法,通常用二分法给出根的大致范围,再利用下面将介绍的更有效的方法求解方程.另一方面,二分法只使用于求一元方程的奇数重实根. ;7.3 一元方程的不动点迭代法; 迭代法是一种逐次逼近法。它是求解代数方程,超越方程及方程组的一种基本方法,但存在收敛性及收敛快慢的问题。; 简单迭代法又称为不动点迭代法,基本思想是首先构造不动点方程 x= ?(x),即由方程 f(x)=0变换为等价形式 x= ?(x), 式中?(x)称为迭代函数。然后建立迭代格式:xk+1 = ?(xk)称为不动点迭代格式。;例1*对应的迭代法分别为;表 7-2;数值分析;数值分析;y; 若从任何可取的初值出发都能保证收敛,则称它为大范围收敛。如若为了保证收敛性必须选取初值充分接近于所要求的根,则称它为局部收敛。
通常局部收敛方法比大范围收敛方法收敛得快。因此,一个合理的算法是先用一种大范围收敛方法求得接近于根的近似值(如对分法),再以其作为新的初值使用局部收敛法(如迭代法)。
这里讨论迭代法的收敛性时,均指的是局部收敛性。;7.3.2 局部收敛性;定理3 (迭代法的局部收敛定理);Y;Y; 有时,对于一些不满足定理7.1的条件问题,可以通过转化,化为适合于迭代的形式。这要针对具体情况进行讨论。;7.3.3 迭代法的收敛速度
一种迭代法具有实用价值,首先要求它是收敛的,其次还要求它收敛得比较快。
定义7.2 设迭代过程 收敛于 的根 ,记迭代误差 ,若存在常数p(p≥1)和c(c0),使 ;定理 设迭代过程 , 若 在所求根
的邻域连续且
则迭代过程在 邻域是p阶收敛的。;例5 已知迭代
文档评论(0)