第一节误差的基本技巧.ppt

于是I7* , I8* 与精确值已经面目全非。 n 精确值 In 近似值In* n 精确值 In 近似值In* 0 1 2 3 4 0.63212? 0.36787? 0.26424? 0.20727? 0.17089? 0.6321 0.3678 0.2642 0.2074 0.1704 5 6 7 8 9 0.14553? 0.12680? 0.11238? 0.10093? 0.09161? 0.1408 0.1120 0.2180 -0.7280 7.5520 算法一 In =1-nIn-1 , 代入得下表 由于计算I0有误差 不计中间再产生的舍入误差 |e(In* )|= n! |e(I0* ) | 到 I8 时 |e(I8* )|= 8!ε = 40320ε 误差扩大了4万倍, 因而该算法是不稳定的。 e(In* ) = - n e( In-1* ) 分析： In=1-nIn-1 , 可以估计出 故 算法二 如果递推式改为 In-1 =(1-In )/n 则 In-1* =(1- In* )/n 误差