数值分析221方程求根(二分法).ppt

* 1.数值计算中的基本原则 2.算法的数值稳定概念 3.方程求根问题引例 4.二分法算法及其应用 《科学计算方法》2 ? ? ? ? 数值计算中的基本原则 (1)避免绝对值小的数做除数; (2)避免两相近数相减; (3)防止大数“吃”小数现象 a = 109，b = 9，在8位浮点数系统中做加法 a + b =1.0000000 ×109+ 0.000000009 ×109 由于只保留8位有效数，处于第九、十位的数09被舍去,实际操作是: 将 a 的数据作为加法计算的最终结果. 2/16 (4)尽量减少计算工作量(乘、除法次数) 例 计算 P(x) = 1+2x+3x2+4x3 + 5x4 的值 秦九韶算法 P(x)=1+x(2+x(3+x(4+5x))) 应用: 2进制数转换为10进制数算法 (1 1 1 0 1 1 1 0)2 = 27+26 +25 +0 +23 +22 +2 +0 =((((((1·2+1)2+1)2+0)2+1)2+1)2+1)2+0 =238 3/16 例1 计算 ( n =0,1,···, 20 ) 4/16 (分部积分法) 初值： I0 = 1 – e –1 ≈0.63212055882856 n=20时，S20= -30.19239488558378 递推公式: In = 1 – nIn-1 (I0 = 1- e-1) S0=1-exp(-1); S(1)=1-S0; for n=2:20 S(n)=1-n*S(n-1) end 实际计算: Sn=1-nSn-1,S0=0.63212055882856 |e(S0)|=|S0 –I0|0.5·10-14 error 5/16 由In=1 - nIn-1 Sn- In=- n(Sn-1 - In-1) 则e(Sn)= –ne(Sn-1)=······= (n!)(–1)ne(S0) 新算法: In-1 = (1 - In)/n S(30)=1/31 for n=30:-1:2 S(n-1)=(1-S(n))/n; end S0=1-S(1),S(1:21) 初值误差在算法执行过程中不断增大,这种算法称为数值不稳定算法。 6/16 Sn-1-In-1= –(Sn - In)/n 和Sn=1-nSn-1得 Sn-1 = (1 - Sn)/n 在算法执行过程中,舍入误差对计算结果影响不大的一类算法被称为数值稳定算法;否则称为不稳定算法. 初始误差在算法执行过程中不断减小,这种算法称为数值稳定算法。 |e(S20)|=|S20-I20|=|(1-S21/21)-(1-I21/21)| =|S21-I21| /21 =·······=|S30-I30|/(21·22·23·····30) 7/16 r d 例2.水中浮球问题 有一半径r =10 cm的球体,密度? =0.638.球体浸入水中后,浸入水中的深度d 是多少？ 根据阿基米德定律,物体排开水的质量就是水对物体的浮力。 整理得: d 3 – 3 r d 2 + 4 r 3? = 0 注意：上式已知物体密度? =0.638, r = 10,水的密度为1. 由? =0.638, r = 10.代入,得d 3 – 30 d 2 + 2552 = 0 令 f (x) = x 3 – 30 x 2 + 2552 ,函数图形如下所示 求解方程 f(x)=0,即 是求函数 f(x)的零 点. f(x) 的零点所 在区间为:[0, 20] 9/16 roots([l -30 0 2552]) ans = 26.3146 11.8615 -8.1761 第一步:对根进行隔离,找出隔根区间,或在隔根区间内确定一个解的近似值x0; 设f(x) = 0的根为 x*,通过迭代计算,产生序列: x0 ? x1 ? x2 ? ··· ? xn·········