Chap1.误差分析.pptVIP

数 值 分 析 硕士研究生学位课 第一章 绪论与误差分析 §1 绪论：数值分析的研究内容 §2 误差的来源和分类 §3 误差的表示 §4 误差的传播 §5 算法设计的若干原则 §1 绪论：数值分析的研究内容 §2 误差的来源和分类 §3 误差的表示 §4　误差的传播 §5 算法设计的若干原则 习 题 一 理论分析：如果 x、y 的近似值分别为 x*、y* ，则两数的差为：z= x-y, z* = x* - y*. 可见，当 x*与y* 非常接近时， x* - y* 作为 x-y 的近似值其相对误差有可能很大。 如何克服该问题？ 当 x 接近零时，可有 当 x 0 很大时，可有 　　如果找不到适当方法，可考虑在计算机上采用双倍字长计算，以增加有效数字，提高精度。 在数值计算中，如果遇到两个近似的数相减运算，可考虑能否改变一下算法以避免两数相减。例如： 当 x1 x2 接近时，可有 例如，在八位十进制计算机上计算 A0.1+0. 9 二. 防止大数“吃掉”小数 参加计算的数，有时数量级相差很大，如果不注意采取相应措施，在它们的加、减法运算中，绝对值很小的数往往被绝对值较大的数“吃掉”，不能发挥其作用，造成计算结果失真。 此时，按照加法浮点运算的对阶规则，应有 由于计算机只能存放八位十进制数，上式中后两个数在计算机上变成“机器零”，计算结果为 A=0108+0.000000001×108+0.000000009×108 A= 0108 = 即相对小数0.1和0.9已被大掉,计算结果失真。 如何克服该问题？ 一般情况下：当一组数进行相加运算时，应按照由小到大的次序进行相加。 如果改变计算次序，现将两个小数相加得到整数1，再进行整数加法运算，就可以比避免上述现象。此时 A=(0.1+0.9) +=1 +三. 绝对值太小的数不宜作除数 在计算过程中，用绝对值很小的数作除数会使商的数量级增加。假设 x 、 y 的近似值分别是 x *、 y* ，则 的近似值是 可见，当|y|很小时， z 的绝对误差可能很大。 此外，当商过大时，或者其数值超出计算机表示的范围引发“溢出”现象，或者作为一个大数它将吃掉参与运算的一些小数。 的值。如果采用逐项计算然后相加的算法： 四.注意简化计算程序，减少计算次数 同一个问题的计算，可以有不同的计算方法。若方法选取得当能减少计算次数，则不仅可提高计算速度，也可减少误差积累。例如，对给定的x，计算多项式 所需的乘法次数为 加法次数为 n 次。 如果把 pn(x) 改写为 采用如下算法：秦九韶算法 这时，只有 n 次乘法，加法次数为 n 次。 五、选用数值稳定性好的算法 利用分部积分法可得到递推公式 如果求得初始值的近似值 一种数值算法，如果其计算舍入误差积累是可控制的，则称其为数值稳定的，反之称为数值不稳定的。数值不稳定的算法没有实用价值。考虑积分计算 则可求得所有积分的近似值： 并且得到误差为： 如果按照此算法计算，当n=20时，误差将会很大，因此该方法需要改进。将积分递推式 改写如下： 再估计第一项 根据 可得 取 这时的误差为： 于是得到改变后的算法 当n=20，19，…，1，0 时，误差将会越来越小。因此改进后的算法是一个稳定的算法。 1-1 下列各数都是经过四舍五入得到的近似值。试分别指出它们的绝对误差限，相对误差限和有效数字的位数。 a=0.0315, b=0.3015, c=31.50, d=5000 1-2 下列近似值的绝对误差限都是0.005, a=-1.00031 , b=0.042 , c=-0.00032 试指出它们有几位有效数字。 1-3 为了使 的近似值的相对误差小于0.01%，试问应取几位有效数字？ 1-4 求方程 x2-56x+1=0 的两个根，使它们至少具有四位有效数字 1-6 设 ,假定 g 是精确的，而对时间 t 的测 量有 ± 0.1s 的误差。证明：当t 增大时，S 的绝对误差增大而相对误差减小. 1-5 若取