数值计算(计算方法第一章)精讲.ppt

算法优劣的标准 从截断误差观点看，算法必须是截断误差小，收敛速 速要快。即运算量小，机器用时少。 从舍入误差观点看，舍入误差在计算过程中要能控 制，即算法的数值要稳定。 从实现算法的观点看，算法的逻辑结构不宜太复杂， 便于程序编制和上机实现. 设计算法时应遵循的原则 要具有数值稳定性，即能控制误差的传播。 避免大数吃小数，即两数相加时，防止较小的数加 不到较大的数上。 避免两相近的数相减，以免有效数字的大量丢失。 避免分母很小或乘法因子很大，以免产生溢出。 三、浮点数与有效数字 定点数：小数点的位置固定在个位数后。 机器数：计算机中可表示的数。 为了提高精度，机器数通常是用浮点数表示的。 称为基数 称为尾数或数码 称为阶码 其中基数是正整数，一般取为２，但为照顾习惯和书写方便，通常化为十进制数输入或输出。阶码是整数。 一定型号的计算机，尾数的位数t是固定的，称为计算机 的位数；阶码m也有一定的取值范围： 有4位有效数字 有6位有效数字 定义1.2.3 有8位有效数字 只有4位有效数字！ 由于计算机只能表示有限个数，故通常利用某种舍入规则 如四舍五入，截断误差等 ，将数进行浮点化。因而势必产生舍入误差。 n+m位有效数字 n-m位有效数字 n位有效数字 如何确定有效数字、绝对误差限、相对误差限？ 说明有效数字位数与小数点的位置无关。只有写成规格化 形式后，小数点后的位数才能反映出其有效位数的多少。 因此，根据上述分析，对有效数字有如下结果： 定理1.2.1 例1.2.3 求下列四舍五入近似值的有效数字位数. 3位 3位 4位 4位 3位 5位 补充 例1.2.4 实际上只有1位！ 试求它们的有效数字位数。 解 k 1, n 2, m 2 例1.2.5 从以上分析可见，四舍五入的近似值的数字都是有效数字 而不是四舍五入得到的近似值的数字不一定是有效数字。 k 3, m 4 n 3 k 4, m 5 n 4 定理1.2.2 证明 下面的结果论述了相对误差与有效数字的关系 补充 即 则有 由定理1.2.1可知， 例1.2.6 解 定理1.2.3 该结论可以参照定理1.2.2的证明，请同学们自证 补充 定理说明:有效数字位数越多相对 误差限就越小，反之亦然。 例1.2.7 解 则根据定理1.2.3，相对误差满足 即应取4位有效数字，近似值的误差不超过0.1%. 四、误差的传播 1、数据误差的传播 由多元函数的Taylor展开公式可得， 的绝对误差为： 相对误差为： 称为 f 的条件数，其绝对值的大小可反映函数值对数据的敏感程度 利用上面的误差估计公式，可以得到两个数的和、差、积、商的误差估计 2、舍入误差的传播 因舍入导致的相对误差限仅与计算机的字长有关，通常 称相对误差限 为计算机的相对精度。 即 在计算机中，数需首先转化为机器数，比如浮点数，在 运算器中参与运算后仍需将运算结果转化成浮点数的形 式进行存储。 由上面的讨论可以看出，为了求得满意的计算解，在选 用计算公式和设计算法时，都应注意如下普遍原则： 1 防止大数吃小数 主要由计算机的位数引起 选用算法应遵循的原则 计算机中数的计算特点： 加法先对阶，后运算，再舍入。 乘法先运算，再舍入。 不在计算机数系中的数做四舍五入处理。 计算机在进行运算时，首先要把参加运算的数对阶， 即把两数都写成绝对值小于1而阶码相同的数。 例1.2.8 在四位浮点十进制数的计算机上计算1+ 104 解 1+ 104 0.1000 ?101+ 0.1000 ?105 0.00001 ?105 + 0.1000 ?105 对阶计算 0.10001 ?105 0.1000 ?105 104 作一个有效数字为4位的连加运算 而如果将小数放在前面计算 在作连加时，为防止大数吃小数，应从小到大进行相加， 如此，精度将得到适当改善。当然也可采取别的方法。 例1.2.9 2 作减法时应避免两个相近数相减 两个相近的数相减，会使有效数字的位数严重损失！ 例1.2.10 用四位浮点数计算 解 只有一位有效数字,有效数字大量损失，造成相对 误差扩大。 结果仍然有四位有效数字。 这说明了算法设计的重要性。在算法设计中，若可能出现两个相近数相减，则改变计算公式，如使用三角变换、有理化等等。 例1.2.11 解方程 解 方程的精确解为 而如果在字长为8，基底为10的计算机上利用求根公式 机器吃了 因此在计算机上 上式是解二次方程的数值公式 的值与精确解差别很大！ 3 避免小数作除数和大数作乘数 小数作除数或大数作乘数会产生溢出错误，因而产生大的误差。 在算法设计时，要避免这类情况在计算公式中出现。此时可以 根据一些具体情况, 把某些算式改写成另一种等价的形式，如 分