[理学]Chap1误差分析.ppt

第一章 绪论与误差分析 数值分析 数 值 分 析 硕士研究生学位课 第一章 绪论与误差分析 §1 计算数学研究的对象和内容 §2 误差的来源和分类 §3 误差的表示 §4 误差的传播 §5 算法设计的若干原则 本章内容安排 §1 计算数学研究的对象和内容 一、计算数学的产生与发展 二、计算数学研究的对象和任务 三、数值分析的学习内容 四、学习要求 §2 误差的来源和分类 §3 误差的表示 §4　误差的传播 §5 算法设计的若干原则 习 题 一 即绝对误差为： 此时，得相对误差为： 例如：对于一元函数 y= f(x) 其中， 例 1-6 测得直角三角形的斜边 c 及一直角边 a 的近似值为c*=75cm,a*=32cm，而且测量误差为 如果计算边 a 对应的角 A 时会产生多大的误差？ 解：由a=c sinA 得到 则由绝对误差估计式： 及 于是： 从而： 即，由于对边的测量产生的误差，影响到角的计算将产生9分的误差。 例 1-7 周长为 10cm 的圆，在计算面积时，欲使其误差不 超过0.1 cm2， 问测量半径时误差应控制在什么范围以内？ 解：首先给出面积的计算公式 其绝对误差为 于是 应有 e(s)= s‘ (r) e(r)= 2πr e(r) 在计算机上进行数值运算时，由于计算机的字长有限，只能保留有限位有效数字，因而每一步计算都可能产生误差，比如计算舍入误差。在反复多次的计算过程中，将产生误差的传播和积累。当误差积累过大时，会导致计算结果失真。因而，为减少舍入误差的影响，设计算法时应遵循如下一些原则。 一、避免两个相近的数相减 在数值计算中，两个相近的数相减会使有效数字受到损失，有效数位减少。例如 都有四位有效数字， 但 x-y=0.005 却仅有一位有效数字。 x=5.143 , y=5.138 事实上，如果 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 的绝对误差可能很大。 此外，当商过大时，或者其数值超出计算机表示的范围而引发“溢出”现象，或者作为一个大数它将吃掉参与运算的一些小数。 的值。如果采用逐项计算然后相加的算法： 四.注意简化计算程序，减少计算次数 同一个问题的计算，可以有不同的计算方法。若方法选取得当能减少计算次数，则不仅可提高计算速度，也可减少误差积累。例如，对给定的，计算多项式 所需的乘法次数为 加法次数为 n 次。 如果把 pn(x) 改写为 采用如下算法：秦九韶算法 这时，只有 n 次乘法，加法次数为 n 次。 五、选用数值稳定性好的算法 利用分部积分法可得到