北京交通大学(数字分析研究生课程)2求根课件11.docVIP

第2章 非线性方程的求根方法 本章探讨函数方程求根的常用数值方 法的构造和原理，主要介绍非线性方程求根 方法的有关知识和方法. 重点论述二分法、简单迭代法、牛顿迭 代法及其变形的原理、构造、收敛性等内容。 2.1实际案例 如下关于角度a的方程 cos a ? (1 - 0.5 tan a)14 - 0.32 = 0 来自某种门的气压控制问题，工程师要求出 满足如上方程的oc的值以设计出一种自动控 2. 2问题的描述与基本概念 定义2.1设几兀丿为一元连续函数，称 f(x)=O (2.1) 为函数方程； 当幷r丿不是兀的线性函数时，称对应的函 数方程为非线性方程。 非线性方程中，当函数几兀丿是多项式函 数时，称为代数方程，否则称为超越方程。 ?n次代数方程 + an + an iX711 + ? ? ? + d^X + 61 q 在非线性方程中，绝大部分是没有求根 公式的，因此寻找求近似根的方法是非常重 要!。 ?求根问题的本质 根的存在性、根的范围和根的精确化根的精确化是方程求根问题的核心。 ?数值分析中求根方法分类 两类：区间法，迭代法 定义2.2设数列｛兀｝收敛于兀*，若存在 正数p和C,满足 .\xk^~x\-r lim c i g - 天 I 则称仪｝的收敛阶为P或方法具有阶敛速。 ?当p = l且Cvl时，称方法线性收敛； ?当#=2时，称方法平方收敛； ? pl时称方法超线性收敛。、 收敛阶越大，收敛越快，方法越好！ 3二分法 二分法也称对分区间法、对分法等，是 最简单的求根方法，属于区间法求根类型。 基本思想 利用连续函数零点定理，将含根区间逐 次减半缩小，构造点列｛?｝来逼近根兀*。 1.构造原理 假设fM = 0在区间[°,切中只有一个根，且满足 /(a)/(b) 0,则二分法求根数列｛“ ｝的构造过程为： 记/。=[讪,取/()的中点竝=0.5(a + b),计算/(%)； 判别于(兀)的值 若/un)-o,则扌=无0,终止； 若 f(xQ)f(a)Of 则X* e[fl,x0],取Ii = [a,X0]； 若/(xo)/(tz)O,则x*g[x0,Z?],取h = [x0,b]; ③记= 胡］,取/i的中点兀1 =0?5（坷+勺） ④判别召-兀＜￡是否成立（￡为给定的精度） 若成立，取兀*2兀1,终止； 否则用，代替/0,转①。 按上述步骤求根的方法称为二分法。 （构造过程一般是重复的,其描述要写出其重复体） ?求根数列比｝描述 若记第k次二分区间处理得到的含根区间为 ?=[%仇],则有二分法对应的求根数列算式为 无=0.5（%+々），“0,1,2,…。 2扮析 因为 x* g 因为 x* g Ik ,k =0,1,…, 兀的误差满足 且 Z0O/lO/2O...,所以有 ~xk 右（切-绞）=連（如一1 -％1）=???=季（b-a）（2.2） 于是当k—g时，由 右（b-d）TO得到无 说明由二分法产生的数列｛忑｝总是收敛于根T的。 ?计算次数控制 给定精度￡0后，要卜”-无|￡成立，取氐满足 (b-a)￡即可，解出匕有 , \n(b - a) - \n , \n(b - a) - \n ￡ t k 1 In 2 (2.3)(事先估计) 这样就可以保证进行氐次二分计算得到的心 就是满 足精度要求的近似根。 式(2.3)确定的E往往偏大，主要用于理论估计。 由二分法的构造，有 故总成立有 故总成立有 （2.4）（事后估计） （2.4）（事后估计） 因而当比- | § ￡时就有F -檢|,此时心是满足 精度的近似根。 式（2.4）确定的E往往较小，主要用于实际控制。 例2?1用二分法求方程2x3-5x-1 = 0在区间[1,2] 内的根，绝对误差沐lot 解 令/(x) = 2x3-5x-1,贝IJ/(x) = 6x2-5, 因为/(D0,几2)0,且广⑴在[1, 2]内不变号。 可知方程fM = 0在[1, 2]内有唯一根。 由二分法算法有 Io = [1,2], x0 = 0.5(1 + 2) = 1.5 由 /(x0)/(l)0,得 /]=[兀0,2] = [1.5,2]，西=0.5(1.5 4-2) = 1.75 由 /(^0)/(^)0,得 Z2 = [1.59x1] = [1.5,1.75], x2 = 0.5(1.5 +1.75) = 1.625 类似计算，可得 /4 = [1.625,1.6875] 9x4 = 1.65625 ； 化=[1.6525,1.6875]，丘=1 ?671875—兀41 = °?° 15625 102 ； /6 = [1.671875, L6875], x6 = 1.679688 , x6-x5 = 0.007813 = 0.7813x 10-2