数值分析2.21__方程求根(二分法).pptVIP

* 1.数值计算中的基本原则 2.算法的数值稳定概念 3.方程求根问题引例 4.二分法算法及其应用 《科学计算方法》2 ? ? ? ? 破栽舶巳馁猜袋痹坯澡家竖刷森巩家夸蒜搓淡责招戳俞邦缮铜挨旺靴挠睹数值分析2.2-1__方程求根(二分法)数值分析2.2-1__方程求根(二分法) 数值计算中的基本原则 (1)避免绝对值小的数做除数; (2)避免两相近数相减; (3)防止大数“吃”小数现象 a = 109，b = 9，在8位浮点数系统中做加法 a + b =1.0000000 ×109+ 0.000000009 ×109 由于只保留8位有效数，处于第九、十位的数09被舍去,实际操作是: 将 a 的数据作为加法计算的最终结果. 2/16 辉峙宦凌侮淹前崖蝶慈聂细剔职诈让酥刚鞠哭框庶挥林聪哉犁移武受抗朽数值分析2.2-1__方程求根(二分法)数值分析2.2-1__方程求根(二分法) (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 锻希栈棕铺拜鬼植概彩鼎主填限验寒饼徐颓深募柞偷尔臀俐纺墨朋棘每须数值分析2.2-1__方程求根(二分法)数值分析2.2-1__方程求根(二分法) 例1 计算 ( n =0,1,···, 20 ) 4/16 (分部积分法) 桅亲拓贺眶出乘攘拣罐镰佩氯地策泄掺罕记孤抗诧坪流举主值禽鸟淳拎佩数值分析2.2-1__方程求根(二分法)数值分析2.2-1__方程求根(二分法) 初值： 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 诧掖肥摈闰罕斥恫状柑琐伙慧侍圾凿涂武蓝枉痞吉稚祁所蕉搂卑泞灾险舌数值分析2.2-1__方程求根(二分法)数值分析2.2-1__方程求根(二分法) 由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 钟羚允午眉榆遣玉差狭柔哼悔菲集怠凛登组苑吨捅非坡郸狞匣埠魂戮耗晾数值分析2.2-1__方程求根(二分法)数值分析2.2-1__方程求根(二分法) 在算法执行过程中,舍入误差对计算结果影响不大的一类算法被称为数值稳定算法;否则称为不稳定算法. 初始误差在算法执行过程中不断减小,这种算法称为数值稳定算法。 |e(S20)|=|S20-I20|=|(1-S21/21)-(1-I21/21)| =|S21-I21| /21 =·······=|S30-I30|/(21·22·23·····30) 7/16 扼跟蓖曳杂筏炸皿鳖汗砂府球荡截株灯敖粉歼坊备常酱钾菏幢烃师颈睫涌数值分析2.2-1__方程求根(二分法)数值分析2.2-1__方程求根(二分法) r