第1章--1(矩阵与数值分析).ppt

现将上述估计式应用到四则运算. 特别地,当 时, (1)加法 两个近似数相加，其运算结果的 精度不比原始数据的任何一个精度高。 (2)减法 应避免两个相近的数相减。 两个近似数相减，其运算结果的 精度不比原始数据的任何一个精度高。 应避免大数作乘数。 (3)乘法 (4)除法 应避免小数作除数。 练习 已知 均为有效数字， 求 的相对误差界。 解：取 得 由已知， 从而 如果 ，则 ，用上述公式计算时 总之，两者其中之一必将会损失有效数字。 一元二次方程 例 有两个根，其求根公式为 , 如果 ，则有 如果 ，则有 一般解二次方程 其中 是b 的符号函数。 （设 均不为零）， 应取 方程 的根为 若取三位有效数字计算，有 ，则由习惯的公式 ，有三位有效数字。 ，只有一位有效数字。 ． 例 而 ，如果改用公式 ，计算有 的精确值是 具有二位有效数字。 时发生了两个相近数相减，造成 其原因为在计算 有效数字损失。 用某一种数值方法求一个问题的数值解，如果在方法的计算过程中舍入误差在一定条件下能够得到控制（或者说舍入误差的增长不影响产生可靠的结果），则称该算法是数值稳定的；否则，即出现与数值稳定相反的情况，则称之为数值不稳定的. 1.2.5 数值方法的稳定性和避免误差危害的基本原则 1.数值方法的稳定性 由于 则递归算法如下： 1. 2. 计算积分 例题 解： 计算出 由 计算出 由 设 的近似值为 ，然后按方法1计算 的近似值 ，如果最初计算时误差为 递推过程的舍入误差不记，并记 ，则有 由此可见，用该方法计算 时， ，那么计算 时产生的舍入误差 放大了 倍，因此，该方法是数值不稳定的。 按方法2计算时， 记初始误差为 ，则有 生的舍入误差为 时产 若计算 方法2不会放大舍入误差。因此，该方法是数值稳定的。 为了用数值方法求得数值问题满意的近似解，在数值运算中应注意下面两个基本原则。 2、避免误差危害的基本原则 (I)避免有效数字的损失 （3）避免小数做除数或大数做乘数。 在四则运算中为避免有效数字的损失，应注意以下事项： （1）在做加法运算时，应防止“大数吃小数”； （2）避免两个相近数相减； 在五位十进制的计算机上计算 解 计算机作加减法时，先将所相加数阶码对齐，根据 因其中 的舍入结果为0，所以上式的计算结果是 。这种现象被称为“大数吃小数”。 后一种方法的结果是正确的，前一种方法的舍入误差影响太大。 例 字长舍入，再加减。 规范化和阶码对齐后的数表示为 相加，再和63015相加，即 次序，先把1000个 1000个 ，那么上式用 如果用63015依次加各个 如果改变运算 例如在八位十进制计算机上，计算 为9位尾数左移变成机器零，这便说明用小数做除数或用 3.712 与 在计算机上做和时， 大数做乘数时，容易产生大的舍入误差，应尽量避免。 3.712由于阶码升 就是所求的值。 如果直接逐项求和计算，需要大约 次乘法运算, 即 例如，多项式求值运算，设 若取 ， (II)减少运算次数 次 次 则有递推公式： 总的计算量需进行 次乘法。 若将公式变成如下递推公式，即令 若令 则有递推公式： 总的计算量为 次乘法。 就是所求的的值。 39 10 5 0 1 -3 1 -1 5 0 1 -3 1 -1 令X=2 2 78 10 2 20 2 42 21 2 2 79 158 157 以上计算过程称之为 秦九韶算法 利用 计算 ， 要计算十万项的和，计算量很大，另 一方面舍入误差的积累也十分严重。 取 只须计算前9项的和，截断误差便小于 ，若要 精确到 如果改用级数 例 秦九韶是一位既重视理论又重视实践，既善于继承又勇于创新的数学家．他所 提出的大衍求一术和正负开方术及其名著《数书九章》，是中国数学史上光彩夺目 的一页，对后世数学发展产生了广泛的影响．美国著名科学史家G．萨顿(Sarton， 1884－1956)说过，秦九韶是“他那个民族，他那个时代，并且确实也是所有时代最 伟大的数学家之一”． 秦九韶的数学成就及对世界数学的贡献主要表现在以下方面： 1、秦九韶的《数书九章》是一部划时代的巨著 2、秦九韶的“大衍求一术”，领先高斯554年，被康托尔称为“最幸运的天才” 3、秦九韶的任意次方程的数值解领先英国人霍纳（