吉林大学数值计算方法第一章数值计算方法与误差精讲.pptVIP

（三）注意简化计算步骤、减少运算次数、避免误差积累 同一个计算问题，如果能减少运算次数，不但可以提高计算速度，而且能减少误差的积累。 简化计算步骤、减少运算次数、避免误差积累的例子 例6 计算多项式 P4(x)=0.0625x4+0.425x3+1.215x2+1.912x+2.1296 的值。 如果先计算各项然后相加，需做十次乘法和四次加 法。 如改用下式计算 (((0.062 5x+0.425)x+1.215)x+1.912)x+2.129 6 则只需做四次乘法和四次加法。 简化计算步骤、减少运算次数、避免误差积累的例子 又如计算 1/(1*2)+1/(2*3)+…+1/(1000*1001) 的值。 若一项一项进行计算，不仅计算次数多，而 且误差积累也很大。若简化成 1－1/1001 进行计 算，则整个计算只要一次求倒数和一次减法。 （四）要避免绝对值小的数作除数 由式 ε(x1/x2)≈d(x1/x2)≈[x2ε(x1)－x1ε(x2)]/ x22 , (x2≠0) 可知，当除数x2接近于零时，商的绝对误差就可能很大。因此 , 在数值计算中要尽量避免绝对值小的数作除数, 避免的方法是把 算式变形或改变计算顺序。 例8 当x接近于0时 (1-cosx)/sinx 的分子、分母都接近0，为避免绝对值小的数作除数，可将原式 化为 (1-cosx)/sinx=sinx/(1+cosx) 例9 当x 很大时，可化 x/[(x+1)0.5－x0.5]=x[(x+1)0.5 + x0.5] 控制误差传播的例子 例10 计算积分 In=∫01 xn ex-1dx，n=0,1, 2, … , 9 利用分部积分法，可得 ? In= xn ex-1| 01 –∫01 ex-1dxn =1– n∫01 xn-1 ex-1dx =1– nIn-1 从而有递推公式 I0= ∫01 ex-1dx= ex-1 | 01 = 1-e-1 ≈0.6321 In= 1– nIn-1 （n=0, 1, 2, … , 9） 计算积分In=∫01 xn ex-1dx的过程 如果直接应用递推公式 I0= ∫01 ex-1dx= ex-1| 01 = 1-e-1 ≈0.6321 In= 1– nIn-1 （n=1, 2, … , 9） 用四位小数计算依次得到： 0.6321, 0.3679, 0.2642, 0.2074, 0.1704 0.1480, 0.1120, 0.2160, -0.7280, 7.5520 由此看到I8为负值、 I9 1，显然与一切0In1 （由于 e-1/(n+1)= min(ex-1)∫01xndxIn (0≤x≤1) max (ex-1)∫01xndx=1/(n+1) ) 矛盾。事实上，从I7开始已经连一位有效数字也没有了 （ I71/8=0.125,而上面算得I7 = 0.2160）。是什么原因 造成这种结果呢？ 计算积分In=∫01 xn ex-1dx 的误差分析 根据 I0=∫01 ex-1dx= ex-1| 01 = 1-e-1 ≈0.6321 In= 1– nIn-1 （n=1, 2, … , 9） 计算In, 假定I0与准确值的误差为ε(I0) 。容易看出误差传递规律 是: ε(In)= –nε(In-1)= … =(-1)nn! ε(I0) （如：I4=1- 4I3=1- 4(1-3I2)=1- 4(1-3(1-2I1