数值代数精品课件-数值代数---第二章第三节.pptVIP

§2.3 基本运算的舍入误差分析 计算机浮点数 IEEE Floating Point Arithmetic 课外知识阅读 IEEE Floating Point Arithmetic Institute for Electrical and Electronices Eigineer. IEEE is pronounced “I triple E” IEEE Floating Point Arithmetic -Published in 1985 /cs/faculty/overton/book Three very important requirements: Consistent representation of floating point number by all machines; Correctly rounded floating operations, using various models; Consistent treatment of exceptional situations such as division by zero. 浮点数与舍入误差 非零的二进制数一般表示 例 考虑 此浮点数集合中最小元素为1/8，最大的元素为3.5。该集合F共有正数20个，负数20个及零元素。由于正、负数是关于零点对称的。 机器精度=1/4。 结论：浮点数分布的不均匀性 程序介绍：第二章密度与误差分析 \floatgui Matlab-使用IEEE双精度二进制 使用64个位存储一个数 浮点数分布的不均匀有影响吗？ 浮点数在接近其下界m时比较稠密，而接近上届M时浮点数则相对比较稀疏。 浮点数接近M时，由于这时浮点数的间隔大，所以近似的绝对误差则会比较大； 只要被近似的数m和M之间，相对误差在各处是一样的（是什么？源程序：epstry;）。 数值计算的焦点问题更多地还是相对误差问题，因此浮点数分布的不均匀并不影响算法的讨论。 计算机有限精度带来的问题 Format long X=4/3-1 Y=3*x Z=1-y 精确计算：z=0 数值结果：z=2.220446049e-16 （舍入得到的） 舍入误差是如何影响数值计算的？ 例 考虑用Matlab程序 基本运算的舍入误差 a和b是两个浮点数，考虑两者之间的加、减、乘、除运算。 何谓上溢？何谓下溢？ 定理2.3.2（公理） 基本浮点运算的相对误差界 浮点数运算满足的法则 内积运算的误差界-例2.3.1 符号介绍“矩阵绝对值” 矩阵大小的定义 矩阵基本运算的舍入误差 矩阵与数乘的舍入误差 矩阵加法 矩阵乘积 上述均为向前误差分析 何谓向后误差分析？具体表达形式。何谓向后稳定的？ 人物介绍：J.Wilkinson 研究的主要问题是线性代数方程组和矩阵的特征值问题的数值方法，特别是他的向后误差分析方法，奠定了数值分析和数值计算的理论基础 其传世之作：The Algebraic Eigenvalue Problem.Clarendon Press.1965 Jim Wilkinson作为数值分析和数值计算的开拓者和奠基人是当之无愧的。 Wilkinson-- 数学软件的开拓者 英国的NAG(Numerical Algorithm Group)1970年正式成立，是世界上最早的数学软件项目之一。当时NAG的大部分线性代数方程组和矩阵特征值计算方面的程序，都是由Wilkinson主持编写的 * * 实数x→计算机对应的浮点数fl(x) 方式一、舍入法 方式二、截断法 满足上述条件的数的集合称为浮点数。 s称为阶码且满足 这里 NaN-Not a number m计算机表示的最小的实数； M计算机表示的最大的实数。 源程序：realmintry;realmaxtry 所给出的浮点数的集合。 7/8 6/8 5/8 4/8 十进制尾数 0.111 0.110 0.101 0.100 二进制尾数 7/8 6/8 5/8 4/8 1 0 7/2 6/2 5/2 4/2 4 2 2 ? ? 2^s 7/4 6/4 5/4 (1/4机器精度) 4/4 1 7/16 6/16 5/16 4/16 -1 7/32 6/32 5/32 4/32 -2 浮点数 指数s 实数x→浮点数fl(x) (舍入或截断) 尾数符号 1位；尾数52位； 阶码（包括符号）11位； 机器精度 Plot(x,y) 源程序：第二章密度与误差