数值计算(计算方法第一章).pptVIP

算法优劣的标准 从截断误差观点看，算法必须是截断误差小，收敛速 速要快。即运算量小，机器用时少。 从舍入误差观点看，舍入误差在计算过程中要能控 制，即算法的数值要稳定。 从实现算法的观点看，算法的逻辑结构不宜太复杂， 便于程序编制和上机实现. 设计算法时应遵循的原则 要具有数值稳定性，即能控制误差的传播。 避免大数吃小数，即两数相加时，防止较小的数加 不到较大的数上。 避免两相近的数相减，以免有效数字的大量丢失。 避免分母很小或乘法因子很大，以免产生溢出。 三、浮点数与有效数字 定点数：小数点的位置固定在个位数后。 机器数：计算机中可表示的数。 为了提高精度，机器数通常是用浮点数表示的。 称为基数 称为尾数或数码 称为阶码 其中基数是正整数，一般取为２，但为照顾习惯和书写方便，通常化为十进制数输入或输出。阶码是整数。 一定型号的计算机，尾数的位数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位的连加运算 而如果将小数放在前面计算 在作连加时，为防止大数吃小数，应从小到大进行相加， 如此，精度将