解非线性方程-1.ppt

第四章 解非线性方程和方程组的迭代法 在实际应用中有许多非线性方程的例子，例如： (1) 在光的衍射理论（the theory of diffraction of light)中,我们需要求 x – tan x = 0 的根 (2) 在行星轨道( planetary orbits)的计算中,对任意的a和b，我们需要求 x – a sin x = b 的根 (3) 在数学中，需要求n次多项式 xn+ a1 xn-1+...+an-1 x + an ＝0的根 Remark1：求奇数个根 Find solutions to the equation x3 - 6x2 +10x – 4 = 0 Remark2: 要区别根与奇异点 例3: 已知方程 2x – 7 - lgx＝0，求方程的含根区间，考查用迭代法解此方程的收敛性。 解：在这里我们考查在区间[3.5,4]的迭代法的收敛性 很容易验证：f (3.5)0, f (4)0 将方程变形成等价形式：x＝( lg x + 7 ) / 2 例4: 用一般迭代法求x3- x -1=0 的正实根x* 将方程化为等价方程：x＝2＋ln x k xi 0 3.000000000 1 3.098612289 2 3.130954362 3 3.141337866 4 3.144648781 5 3.145702209 6 3.146037143 7 3.146143611 另一种迭代格式： * 简介（Introduction） 求 f (x) = 0 的根 本章目的：讲述用于实际计算中求 f (x) = 0 的根的近似值的几种常用方法。 方程根的数值计算大致可以分为三个步骤： (1)判断根的存在性； (2)确定根的分布范围(根的隔离)； (3)根的精确化。 根的隔离 1. 分析法：利用对函数 f (x) 的各种性质的分析来确定根的分布范围。 例：试确定f (x) = x3- 6x2 + 9x-1=0各根的分布范围。 隔根区间为(0,1)、(1,3)、(3,4) 2. 逐步搜索法： 先确定方程 f (x) = 0的所有实根所在的区间[a , b]，再按照选定的步长 h =(b –a )/n，取点 xk = a + k h ( k =0,1,2,…., n)，逐步计算函数值 f (xk), 依据函数值异号及实根的个数确定隔根区间. 必要时可以调整步长 h, 总可以把隔根区间全部找出. 代数方程根的模上下界定理： 定理：设代数方程 f (x)= xm+am-1xm-1+…+a1x+a0=0 则: (1) 若μ=max{|am-1|,…,|a1|,|a0|}, 方程根的模小于μ+1； (2) 若 v =1/|a0| max{1, |am-1|,…,|a1|}, 方程根的模大于1/(v+1). 例：求方程 x3-3.2 x2 +1.9 x +0.8 = 0 的隔根区间。 解：设方程的根为α 由 μ = max{ |-3.2| , |1.9| , |0.8| } = 3.2 ; v = (1/0.8) max{1, |-3.2| , |1.9| } = 4 得： 0.2 | α | 4.2 即： -4.2 α -0.2 , 0.2 α 4.2 取 n = 8 , h = 0.5 , 计算 f ( xk) … 0.14 -0.31 0.20 … 1.06 0.28 -2.44 … f (xk) … 2.2 1.7 1.2 … 0.2 -0.2 -0.7 … xk 由上表可知隔根区间为 [-0.7 , -0.2] , [1.2 , 1.7] , [1.7 , 2.2] 3. 图解法： 由函数图像来确定根的大体位置。 §4.1 二分法(对分区间法) (Bisection Method ) 原理：若 f (x) ?C[a, b]单调，且 f (a) · f (b) 0，则 f (x) 在 (a, b) 上有且仅有一实根。 基本思想：通过计算隔根区间的中点，逐步将隔根区间缩小，从而得到方程的近似根数列{ xn }。 a b x1 x2 a1 b2 x* b1 a2 (1) 先将[a , b]等分为两个小区间,分点记为x0=(a + b)/2, 判断根属于哪个小区间,舍去无根区间保留有根区间[a 1, b1].即,若 f (x0)=0, 则x0= x*. 设 f (x0) ≠0, 若 f (a) f (x0)0, 则x*∈ (a , x0),