第2章 节 一元线性方程的解发 计算方法 《代码优化》课件.pptVIP

第2章 一元线性方程的解发;§1二分法 ; 等等,看上去形式简单,但却不易求其准确根。为此,只能求方程达到一定精度的近似根。 方程的形式很多,我们主要讨论一元非线性方程,也即 f(x)=0 (2―1) ; 方程(2―1)可以有实根,也可以有复根或者重根等。本章主要讨论它的实根的数值计算问题。 方程根的数值计算大致可分三个步骤进行: (1) 判定根的存在性。 (2)确定根的分布范围,即将每一个根用区间隔离开来。 (3)根的精确化,即根据根的初始近似值按某种方法逐步精确化,直至满足预先要求的精度为止。 ; 设f(x)为定义在某区间上的连续函数,方程(2―1)存在实根。虽然方程(2―1)的根的分布范围一般比较复杂,但我们不难将函数f(x)的定义域分成若干个只含一个实根的区间。 例如考虑方程 x2-2x-1=0 由图2.1所示,该方程的一个负实根在-1和0之间,另一个正实根在2和3之间。 ; 设函数f(x)在区间［a,b］上单调连续,且 f(a)·f(b)＜0 则方程(2―1)在区间(a,b)内有且仅有一个实根x。下面在有根区间(a,b)内介绍二分法的基本思想。 计算f(a)与f(x0),若 ? f(a)·f(x0)＜0 ? 则根x∈(a,x0),令 a1=a,b1=x0 否则x∈(x0,b),令 ? a1=x0,b1=b;图 2 .2 ; 如此逐次往复下去,便得到一系列有根区间 (a,b),(a1,b1),(a2,b2),…,(ak,bk),… 其中 ; 我们把每次二分后的有根区间(ak,bk)的中点 ? ; 则结果xk就是方程(2―1)满足预给精度ε的近似根,也即 ;1.计算步骤 ①输入有根区间的端点a,b及预先给定的精度ε; ②(a+b)/2 x; ③若f(a)f(x)＜0,则xb,转向④;否则xa,转向④。 ④若b-a＜ε,则输出方程满足精度的根x,结束;否则转向②。 ; 2. 计算框图 例1 求方程 ? f(x)=x3-x-1=0 ?在区间(1,1.5)内的根。要求用四位小数计算,精确到10-2。 解 这里 a=1,b=1.5 ??区间(1,1.5)的中点 ; 图 2.3 ; 由于f(1)＜0,f(1.25)＜0,则令 ?a1=1.25, b1=1.5 ?得到新的有根区间(1.25,1.5) ; 表 2―1 ;§2 迭代法 ; 在x=1.5附近的一个根(用六位有效数字计算)。 首先将原方程改写成等价形式 ; x1与x0相差较大,如果改用x1作为近似根代入式(2―5)的右端得 ; 表 2―2 ; 对于一般形式的方程(2―1),首先我们设法将其化为下列等价形式 x=g(x) (2―7) 然后按(2―7)构造迭代公式 从给定的初始近似根x0出发,按迭代公式(2―8)可以得到一个数列 ?x0,x1,x2,…,xk,… 若这个数列｛xk｝有极限,则迭代公式(2―8)是收敛的。此时数列的极限; 就是原方程(2―1)的根。 虽然迭代法的基本思想很简单