数值计算方法(第2章)

《数值计算方法(第2章)》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第2章非线性方程与方程组的数值解法本章重点介绍求解非线性方程的几种常见和有效的数值方法,同时也对非线性方程组求解,简单介绍一些最基本的解法.无论在理论上,还是在实际应用中,这些数值解法都是对经典的解析方法的突破性开拓和补充,许多问题的求解,在解析方法无能为力时,数值方法则可以借助于计算机出色完成.例2.1.将一半径为r,密度为ρ<1的球浸入水中,求球体没入水中的高度.解:2.1二分法求非线性方程确定方程的有根区间计算根的近似值的根的方法分为两步：首先确定有限区间：依据零点定理。设，且，则方程在区间上至少有一

2、个根。如果在上恒正或恒负，则此根唯一。二分法用二分法(将区间对平分)求解。令若，则为有根区间，否则为有根区间记新的有根区间为，则且对重复上述做法得且设所求的根为，则即取为的近似解求方程f(x)=0的根的二分法算法例题例2设方程kakbkxkf(xk)012…12131.01.01.25…1.3649902351.3651123052.01.51.5…1.3652343751.3652343751.51.251.375…1.3651123051.3651733902.375-1.7970.162…-0.002

3、0.00007等步长扫描法求有根区间用计算机求有根区间：等步长扫描法。设h>0是给定的步长，取，若则扫描成功；否则令，继续上述方法，直到成功。如果则扫描失败。再将h缩小，继续以上步骤。等步长扫描算法算法：（求方程的有根区间）(1)输入；(2)；(3)，若输出失败信息，停机。(4)若。输出，已算出方程的一个根，停机。(5)若。输出为有根区间，停机(6)，转3）注：如果对足够小的步长h扫描失败。说明：在内无根在内有偶重根例题例1设方程解：取h=0.1,扫描得：又即在有唯一根。Δaa+Δa+2Δa+3Δa+4Δh

4、<0<0<0-0.610.344h/2<0+0.009375h/4<0+0.0012h/8<0-0.0505+0.0012h/16<0-0.0775h=0.1,a=1,a=1.3,b=1.4,a=1.3,b=1.35.a=1.3,b=1.325,a=1.3,b=1.3125,a=1.3125,b=1.3252.2一般迭代法2.2.1迭代法及收敛性对于有时可以写成形式如：迭代法及收敛性考察方程。这种方程是隐式方程，因而不能直接求出它的根，但如果给出根的某个猜测值，代入中的右端得到，再以为一个猜测值，代入的右端

5、得反复迭代得迭代法及收敛性若收敛，即则得是的一个根例题例2设方程kxk012341.51.373333331.3652620151.3652300141.365230013该方法比二分法快!迭代法的几何意义交点的横坐标y=x简单迭代法将变为另一种等价形式。选取的某一近似值，则按递推关系产生的迭代序列。这种方法算为简单迭代法。例题例2.2.1试用迭代法求方程在区间(1，2)内的实根。解：由建立迭代关系k=10,1,2,3…….计算结果如下:例题精确到小数点后五位例题但如果由建立迭代公式仍取，则有，显然结果越来

6、越大，是发散序列迭代法的收敛性定理2.2.1（压缩映像原理）设迭代函数在闭区间上满足（1）（2）满足Lipschitz条件即有且。压缩映像原理则在上存在唯一解，且对，由产生的序列收敛于。并且有误差估计式:压缩映像原理证明：不失一般性,不妨设否则为方程的根。首先证明根的存在性令压缩映像原理则，即由条件2）是上的连续函数是上的连续函数。故由零点定理在上至少有一根压缩映像原理再证根的唯一性设有均为方程的根则因为0

7、原理误差估计若满足定理2.2.1条件，则这是事后估计，也就是停机标准。L越小，收敛速度越快。这是事前估计。选取n，预先估计迭代次数。例题例2.2.2证明函数在区间[1，2]上满足迭代收敛条件。证明：例题例题若取迭代函数，不满足压缩映像原理，故不能肯定收敛到方程的根。简单迭代收敛情况的几何解释压缩映像原理推论推论设迭代函数在闭区间上满足定理2.2.2设迭代函数在上连续可微2.2.2Steffensen加速收敛法迭代法收敛的阶定义2.2.1设序列收敛到，若有实数和非零常数C，使得其中，，则称该序列是p阶收敛的，

8、C称为渐进误差常数。迭代法收敛的阶当p=1时，称为线性收敛(01时，称为超线性收敛；当p=2时，称为平方收敛或二次收敛。迭代法收敛的阶定理2.2.2’设是方程的不动点，若为足够小的正数。如果且，则从任意出发，由产生的序列收敛到，当时敛速是线性的。迭代法收敛的阶证明：满足压缩映像原理迭代法收敛的阶敛速是线性的线性收敛到。Steffensen迭代格式由线性收敛知当n充分大时有即Steffensen迭