数值与计算方法chapter1-2.pptVIP

* 机器上用很小的数作除数时，可能会因为数据溢出而造成停机； 很小的数作除数，若除数有少量误差，会造成结果出现很大误差。 计算时既要避免两个相近的数相减，更要避免用相减之后的差作除数。 0.0011 ≈2471.1 1.6.2 避免绝对值小的数作除数 * 1.6 误差危害的防止 避免相近数相减 避免用绝对值很小的数作除数 防止大数吃掉小数 简化计算步骤，减少运算次数 使用数值稳定的算法 * 1.6.3 防止大数吃掉小数 a = 109，b = 9，设想在8位浮点数系中相加 a + b =01010+ 0.10-1 = 01010+0.0000000009 ×1010 = 01010 由于只保留8位有效数，数据09被舍去,实际加法操作 a + b计算结果是 将 a 的值作为计算结果赋给 a+ b. 多个数相加，应按绝对值从小到大的顺序依次进行，以避免被大数吃掉。 * 例1-13：计算0.4994+1000+0.0006000+0.4090, 保留4位有效数字。 解： 0.4994+1000≈1000 1000+0.0006000≈1000 1000+0.4090 ≈ 1000 1.6.3 防止大数吃掉小数 * 改变顺序： 0.4994+0.0006000 ≈ 0.5000 0.5000+0.4090 ≈ 0.9090 0.9090+1000 ≈ 1001 1.6.3 防止大数吃掉小数 * 1.6 误差危害的防止 避免相近数相减 避免用绝对值很小的数作除数 防止大数吃掉小数 简化计算步骤，减少运算次数 使用数值稳定的算法 * 1.6.4 减少运算次数 例1-14： 计算x31的值。 方法一：将x的值逐个相乘，共需30次乘法。 方法二： x31= x x2 x4 x8 x16 而 x2 = x x ，x4 = x2 x2，x8 = x4 x4 ，x16 = x8 x8 共8次乘法。 * 例1-3： 对给定的x，求多项式 方法一： 的值。 直接计算 的每一项，再逐项求和 。 乘法次数： 1.6.4 减少运算次数 * 方法二： 乘法次数： ，加法次数： …… 1.6.4 减少运算次数 练习：请写出此算法的递推公式。 * 秦九韶算法 1.6.4 减少运算次数 编程计算： 递推公式 * 求多项式值的秦九韶算法 P(x)=a0 xn + a1xn-1 + ······ + an-1 x+ an 输入 x；a0，a1，…，an b← an；u←1 k 从 n到 1循环 u←x×u b←b+ak-1× u 输出数据b ；结束 输入 x；a0，a1，…，an b ← a0 k 从 1到 n循环 b←ak+ x×b 输出数据b ；结束 秦九韶算法 * 求多项式值的秦九韶算法 a0 a1 a2 …… an-1 an x=x0 b0x0 b1x0 …… bn-2x0 bn-1x0 b0 b1 b2 …… bn-1 bn 例1-15： 求 f(x)=2+x-x2+3x4 在x0=2 的值。 3 0 -1 1 2 x0=2 6 12 22 46 3 6 11 23 48 =f(x0) =f(2) * 1.6 误差危害的防止 避免相近数相减 避免用绝对值很小的数作除数 防止大数吃掉小数 简化计算步骤，减少运算次数 使用数值稳定的算法 1.6.5 使用数值稳定的算法 例1-16: 建立积分 的递推关系式。 解：因为 递推关系式： 初始值： (1) 1.6.5 使用数值稳定的算法 计算I0时的舍入误差设为e0，且求得的I0的近似值为 (1)-(2)式，并不断递推，可得 e0对第n步的影响是扩大了5n倍，算法不稳定 ；n较大时，误差会非常大。 (2) 1.6.5 使用数值稳定的算法 改变算法，可以使初始误差减小，算法稳定。 如果第n＋1步的误差en＋1与第n步的误差en满足 则称此计算公式是绝对稳定的，否则称为不是绝对稳定的。 * 小结 数据误差的影响 误差危害的防止方法 习题一 4（1），8，10，11 说明： 第4题：要求同时计算绝对误差限和相对误差限（利用第1题的数据进行计算）。 第8题，不必计算出y100的值，只需估计误差。 第11题：要求：同时采用手工计算，和编程计算，并给出二者的计算结果 * 第一章 绪论（二） 有效数字 数据误差的影响 误差危害的防止 小结 * 1.5.3