- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
计算方法7.1-7.3综述
第七章 非线性方程的解法;*; 在科学研究和工程设计中, 经常会遇到的一大类问题是非线性方程; 当f(x)不是x的线性函数时,称对应的函数方程为非线性方程。如果f(x)是多项式函数,则称为代数方程,否则称为超越方程(三角方程,指数、对数方程等)。一般称n次多项式构成的方程 ; ;本章介绍方程的迭代解法,它既可以用来求解代数方程,也可以用来解超越方程,并且仅限于求方程的实根。
运用迭代法求解方程的根应解决以下两个问题:
确定根的初值;
将进一步精确化到所需要的精度。; 为了确定根的初值,首先必须圈定根所在的范围,
称为圈定根或根的隔离。
在上述基础上,采取适当的数值方法确定具有一定
精度要求的初值。
对于代数方程,其根的个数(实或复的)与其次数
相同。至于超越方程,其根可能是一个、几个或无
解,并没有什么固定的圈根方法
求方程根的问题,就几何上讲,是求曲线 y=f (x)
与x轴交点的横坐标。; 由高等数学知识知, 设f(x)为区间[a,b]上的单值连续, 如果f(a)·f(b)0 , 则[a,b]中至少有一个实根。如果f(x)在[a,b]上还是单调地递增或递减,则仅有一个实根。;(1) 画图法;对于某些看不清根的函数,可以扩大一下曲线;y;运用零点定理可以得到如下逐步搜索法:;例 求方程 x3-3.2x2+1.9x+0.8=0的隔根区间。;例1 方程f(x)=x3-x-1=0 确定其有根区间
解:用试凑的方法,不难发现 f(0)0 f(2)0
在区间(0,2)内至少有一个实根
设从x=0出发,取h=0.5为步长向右进行根的
搜索,列表如下; 用逐步搜索法进行实根隔离的关键是选取步长h
要选择适当h ,使之既能把根隔离开来,工作量
又不太大。
为获取指定精度要求的初值,可在以上隔离根的
基础上采用对分法继续缩小该含根子区间
二分法可以看作是搜索法的一种改进; 对分法又称二分区间法,是求解方程(1)的近似根的一种常用的简单方法。
设函数f(x)在闭区间[a,b]上连续,且f(a)f(b)0,根据连续函数的性质可知, f(x)=0在(a,b)内必有实根,称区间[a,b]为有根区间。为明确起见,假定方程f(x)=0在区间[a,b]内有惟一实根x*。
二分法的基本思想: 首先确定有根区间,将区间二等分, 通过判断f(x)的符号, 逐步将有根区间缩小, 直至有根区间足够地小, 便可求出满足精度要求的近似根。;① 取有根区间[a,b]之中点, 将它分为两半,分点 ,这样就可缩小有根区间;; ② 对压缩了的有根区间 重复①的步骤,即取中点
,将区间 再分为两半,然后再确定有
根区间 ,其长度是 的二分之一;
③ 如此反复下去,若不出现 ,即可得出一系列有
根区间序列:
上述每个区间都是前一个区间的一半,因此 的长
度:; 每次二分后,取有根区间 的中点
作为根的近似值,得到一个近似根的序列
该序列以根x*为极限
只要二分足够多次(即k足够大),便有
这里ε为给定精度,由于 ,则 ;;例1 求方程f(x)=x3-x-1=0在区间[1.0,1.5]内 的一
个实根, 使误差不超过0.5×10-2。
例2 证明方程 在区间[2,3]内有一个根,
使用二分法求误差不超过0.5×10-3 的根要二
分多少次?
证明 令 ; 误差限为 只要取k满足 ; 二分法的优点是不管有根区间[a,b]多大,总能求出满足精度要求的根,且对函数f(x)的要求不高,只要连续即可,计算简单;它的局限性是只能用于求函数的实根,不能用于求复根及重根,它的收敛速度与比值为1/2的等比级数相同,不算太快,因此一般在求方程近似根时,不太单独使用,常用它来为其它方法求方程近似根提供好的初始值。 ;;*;*;*;*;*;迭代法的几何意义 ;*;*;*;2.迭代法的加速;当 范围不大时,设 变化不大,其估计值为L,则有 ;令;例5. 用加速迭代法求方程 在0.5附近的根。
解: 因为在 附近
取L=-0.6,建立如下迭代公式;(2)埃特金(Aitken)加速法
在加权法中, 估计L的值有时不太方便
文档评论(0)