地球物理计算方法 第四章方程求根的迭代法
- 格式:pdf
- 大小:815.07 KB
- 文档页数:25
计算方法4方程求根的迭代法四方程求根的迭代法是一种用于解决非线性方程的数值方法。
在计算方法中,非线性方程指的是形如f(x)=0的方程,其中f(x)包含x的非线性项。
在实际中,非线性方程的求解是非常常见的问题,因此有很多不同的迭代法可以用于解决这些问题。
以牛顿迭代法为例,它是一种基于线性近似的迭代方法。
该方法的基本思想是将非线性方程转化为线性方程,通过不断迭代来逼近方程的根。
具体而言,牛顿迭代法的步骤如下:1.选择初始估计值x0作为方程的根,并计算f(x0)的值。
2.计算f(x)的导数f'(x),并计算方程的线性近似式x-x0=-f(x0)/f'(x0)。
3.计算下一个近似值x1,即x1=x0-f(x0)/f'(x0)。
4.判断,x1-x0,是否小于给定的收敛条件,如果是则停止迭代,否则转到步骤55.将x1作为新的近似值x0,转到步骤2牛顿迭代法具有快速收敛的特点,尤其适用于具有单根的方程。
然而,该方法也存在一些限制,如在计算f'(x)时需要知道方程的导数,当方程的导数不易计算时,该方法可能不适用。
除了牛顿迭代法,还有其他一些常用的四方程迭代方法,如割线法、弦截法等。
每种方法都有其特点和适用范围,选择合适的方法对于求根问题的解决至关重要。
总结起来,四方程求根的迭代法是一种用于解决非线性方程的数值方法。
牛顿迭代法是其中一种常用的方法,通过不断迭代来逼近方程的根。
根据方程的特点和计算条件,选择合适的迭代方法是解决求根问题的关键。
希望以上的介绍可以帮助您更好地理解和应用这一方法。
Newton迭代法第四章一元方程求根/非线性方程组数值解法初步4.3 Newton 迭代法 1. Newton 迭代法解一元非线性方程组0)(=x f (4.3.1)的Newton 迭代法是不动点迭代法的一种特殊形式。
可从不同途径导出Newton 迭代公式,这里采用T aylor 展开。
设方程0)(=x f 的根*x 的一个近似值0x ,将)(x f 在0x 附近展开得20000)(!2)(''))((')()(0x x f x x x f x f x f -+-+==ξ 或表示为 200000)()('2)('')(')(x x x f f x f x f x x ---=ξ (4.3.2)其中设0)('0≠x f ,''f 存在、连续,而ξ在x 与0x 之间。
忽略上式最后一项*x 的一个新近似值)(')(0001x f x f x x -= 把1x 代替上式右端的0x ,并设0)('1≠x f ,于是又得新近似值)(')(1112x f x f x x -=如此继续,可知当),2,1,0(0)(' =≠k x f k 可得),2,1,0()(')(1 =-=+k x f x f x x k k k k (4.3.3)这就是著名的Newton (牛顿)迭代公式。
在迭代序列收敛的情况下,取一定精度的迭代值kx 作为方程0)(=x f 的根*x 的近似值,这就是解方程组0)(=x f 的Newton 迭代法。
显然,它以在*x 附近函数0)(=x f 线性化为基础,并以),2,1,0(0)(' =≠k x f k 为前提。
例 4.3.1 用Newton 迭代法求下列方程的近似根:1-x xe =0解令 1)(-=x xe x f ,则x x xe e x f +=)(',于是迭代公式为),1,0(11 =+--=+k ex e e x x x kkkxk x x k k k 整理得),1,0(11 =+--=-+k x e x x x kxk k k k取5.00=x 请同学们自己动手完成2.Newton 迭代法的收敛性Newton 迭代公式作为不动点迭代,其迭代函数为 )(')()(x f x f x x -=? 从而有222)]('[)('')()]('[)('')()]('[1)('x f x f x f x f x f x f x f x =--=? 可见,如果在方程0)(=x f 的根*x 的某个邻域内0)('0≠x f (从而有0)('*≠x f ,即*x 是单根的情况),''f 存在并连续(从而有界),则只要x 足够靠近*x ,(从而|)(|x f 足够靠近0),就有1|)('|<≤L x,于是根据定理4.2.1的推论,Newton 迭代公式收敛于*x ,并且0)(*=x f 导致 0)('*=x ?,于是又根据收敛阶的判定定理4.2.4 ,可知Newton 迭代公式在单根附近至少是2阶的。
方程求根的迭代法一、考核知识点:区间二分法,弦位法(单点弦法、双点弦法)、切线法、一般迭代法,收敛性。
二、考核要求:1.熟练掌握用区间二分法求方程近似根的方法。
2.熟练掌握用单点弦法、双点弦法求方程近似根的方法。
了解其收敛性。
3.熟练掌握用切线性求方程近似根的方法。
了解其收敛性。
4.掌握用一般迭代法求方程的方法近似根的方法。
了解其收敛性。
三、重、难点分析例1 证明计算)0(>a a 的切线法迭代公式为:,1,0),(211=+=+n x a x x nn n 并用它求2的近似值(求出1x 即可)解 (1)因计算a 等于求02=-a x 正根,a x x f -=2)(,x x f 2)(=' 代入切线法迭代公式得)(21221nn n n n n x a x x x x x +=-=+ ,1,0=n (2) 设2)(2-=x x f ,因,0121)1(2<-=-=f 025.1)5.1(2>-=f 所以 []5.1,12*∈=x在[]5.1,1上 02)(>='x x f 02)(>=''x f由 0)()(0≥''x f x f ,选5.10=x用上面导出的迭代公式计算得 4167.11217)2(21001≈=+=x x x例2用单点弦法求方程 0153=+-x x 的最小正根(计算出1x ) 解:由于0375.1)5.0(,01)0(<-=>=f f 则]5.0,0[*∈x 在[0,0.5],,06)(,053)(2≥=''<-='x x f x x f 由,0)()(≥''x f c f 取5.0,00==x c 则单点弦法迭代公式 ,1,0)15(51),15(151032331=+--+=+--+---=+n x x x x x x x x x x x n n n n n n n n n n n 计算得 21.075.4375.15.01≈-=x 例3 用双点弦法,一般迭代法求0243=-+x x 的最小正根(求出2x 即可)。
§4.1 引 言绪论中讲到方程求根得二分法,但二分法收敛速度慢,有必要掌握新的方法。
§4.1.1迭代法的思想迭代法是一种逐次逼近法,使用某个固定公式(迭代公式)反复校正,逐步精确,直到满足精度。
迭代法求根分两步: 1) 猜测初值 2)迭代如求解初值问题00')(),,(y x y y x f y ==用梯形公式111[(,)(,)2n n n n n n h y y f x y f x y +++≈++ (1)看作关于1+n y 的函数方程,按欧拉公式提供猜测值),()0(1n n n n y x hf y y +=+,代入(1)得)],(),([2)0(11)1(1+++++=n n n n n n y x f y x f h y y若)1(1+n y 仍不满足要求,则将它代入(1)式,继续得到校正值)2(1+n y ,写成迭代公式)],(),([2)(11)1(1k n n n n n k n y x f y x f h y y ++++++= (2)一般地,为了求一元非线性方程0)(=x f 的根,可以先将其转换为如下的等价形式()x x ϕ= (3)式(3)中连续函数()x ϕ称为迭代函数,其右端含未知数,不能直接求解。
先用根的某个猜测值0x 代入(3),构造迭代公式:()k k x x ϕ=+1。
如果迭代值k x 有极限,则称迭代收敛,极限值k k x x ∞→=lim *就是方程(3)的根。
几何意义P127图4-1为使迭代法有效,必须保证它的收敛行,()x ϕ满足什么条件,才能保证收敛?以最简单的线性迭代()d kx x +=ϕ,可以看出收敛的充分必要条件()1'<=k x ϕ。
几何意义P127图4-2,3,4,5。
§4.1.3 压缩映像原理设*x 是方程()x x ϕ=的根,则由微分中值定理))(()()(*'*1*k k k x xx x x x-=-=-+εϕϕϕ,如果存在10<≤L ,使得],[b a x ∈有()k k x x L x x L x -≤-⇒≤+*1*'ϕ,则迭代误差0e L e kk ≤,由于10<≤L ,故0→k e ,即迭代收敛。
方程求根的迭代法--------------牛顿法和弦截法实验目的:实际问题中碰到的函数f(x)是各种各样的,有的表达是很复杂,这时求解函数方程f(x)=0的根就会变得很困难,而工程应用中,对计算的结果只要保证在某个误差范围之内就足够了,这就要求我们设计一种方法能够求解复杂函数方程的根,迭代法就是这样一种求解复杂函数方程的根的方法。
又由于通常对某些函数方程的根要求比较精确,误差要控制在一定的范围之内,所以计算过程比较复杂,这就要求这种算法能够便于利用计算机编程实现,牛顿法和弦截法为我们提供了一种既能够利用迭代法求解复杂方程的根,又能够便于利用计算机实现,从而方便我们求解,节省我们时间的求解方法。
本实验的目的就在于熟练掌握利用牛顿法和弦截法编程求解方程根,并且比较牛顿法和弦截法的收敛速度,比较两者的不同之处,进而加深对牛顿法和弦截法这两种方法的数学原理的理解。
实验编程实现方程f=@(x)x-exp(-x)及f=@(x)x*exp(x)-1分别利用牛顿法和弦截法求解,并且比较牛顿法和弦截法各自迭代多少次才能达到要求的某一精度,并且理解二者的差别和各自的优缺点。
本实验中要注意在牛顿法的迭代过程中,分母上会出现f(x)的导数项,所以在编程过程中,一定要在程序的一开始就先判断f'(x)是否等于0,如果等于0,则直接跳出,只有在其不为0的情况下才继续执行程序。
另外某些函数可能会发散或者是迭代过程收敛的非常慢,这时用牛顿法就是不合适的,所以这就要求我们另取其他方法。
所以在用牛顿法求解方程根的过程中,要设置一个最大的迭代次数,以免造成电脑资源不足,系统崩溃。
在弦截法的迭代过程中,除了与牛顿法相同的一些注意事项之外,还要注意在计算之前必须要先提供两个开始的值x0,x1。
实验原理: 迭代法的设计思想:迭代法是一种逐次逼近法,这种方法使用某个固定公式,即所谓的迭代公式,反复校正根的近似值,使之逐步精确化,直到得出满足精度要求的结果。