calcu-1计算方法要点.ppt

* * * * * * * * * * * * * * * * * * * * * * 绪论—机器数系 数的符点表示 X=(±0.a1a2…at) βp=s * βp 尾数部：s=±0.a1a2…at=±(a1/β+a2/β2…at/βt) 定位部： βp 基数：β，一般=2，8，10，16 计算机的字长：t 阶码（指数）：p，有固定的下限L和上限U L, U, t由计算机的硬件所决定的常数 规格化符点数：a1≠0，即β-1≤|s|1。 绪论—机器数系 机器数系 计算机中浮点数所组成的集合加上“机器零”，记为F（机器数系） F被4个参数描述：基数β、字长t、阶码范围[L, U] 绪论—机器数系 机器数系 例：一个二进制的2位字长的计算机，即β=2，t=2，其阶码p的范围：[-1，1]，则该计算机能表示的数只有13个数： p=-1时，±(0.10*2-1)2=±(0.25)10, ±(0.11*2-1)2=±(0.375)10 p=0时，±(0.10*20)2=±(0.5)10, ±(0.11*20)2=±(0.75)10 p=1时，±(0.10*21)2=±(1)10, ±(0.11*21)2=±(1.5)10 机器零0 -1 -0.5 -0.25 0.25 0.5 1 -1.5 -0.75 -0.375 0 0.375 0.75 1.5 X 绪论—机器数系 机器数系 X=(±0.a1a2…at) βp=s * βp 计算机中浮点数所组成的集合加上“机器零”，记为F（机器数系） F是一个离散的、分布不均匀的、可数的有限集合 实数是一个稠密的、连续的、不可数的无限集合 绪论—机器数系 机器数系 集合F的数据个数：2*(β-1)*βt-1*(U-L+1)+1 而且这些数不等间隔地分布在数轴上 例：当β=10，t=4，-L=U=99时，计算机的机器数系F仅含有3582001个数 绝对值最小的非零数：±0.1000*10-99 最大的数：0.9999*1099 最小的数：-0.9999*1099 溢出：如计算结果超出了上述范围，则为溢出 绪论—机器数系 机器数（fl(x)）：一个实数x进入计算机后，成为计算机里的数，称为x的机器数。 一般fl(x)只是x的一个近似值。 例：对于x=0.6，在上述二进制2位字长计算机的机器数系F中找不到这个数，取与x最靠近的数0.5作为x的近似值，即x的机器数fl(0.6)=(0.10*20)2=0.5 目前计算机分为截断机和舍入机两种 截断机：fl(x)取x的前t位数字 舍入机：fl(x)按四舍五入原则取x的前t位数字 绪论—机器数系 机器数的误差限和相对误差限 设x=(±0.a1a2…atat+1…) βp, a1≠0 对于舍入机： 当|at+1|≥β/2时，fl(x)= x=(±0.a1a2…at+1) βp 当|at+1|β/2时，fl(x)= x=(±0.a1a2…at) βp 则：|x-fl(x)|≤1/2* β-t βp 绪论—机器数系 机器数的误差限和相对误差限 设x=(±0.a1a2…atat+1…) βp, a1≠0 对于截断机： |x-fl(x)|≤β-t βp 所以： 第一章 绪论 研究对象和特点 误差来源及基本概念 机器数系 误差危害的防止 绪论—误差危害的防止 误差分析在数值运算中是一个非常重要的问题，因为每一步运算都有可能产生误差，而一个工程或科学计算问题往往计算成千上万次，如果每步都分析误差，是不可能也是没必要的。 这里提出误差分析的几个原则，可以鉴别计算结果的可靠性，并且防止误差危害现象的产生。 绪论—防止误差危害的原则 使用数值稳定的计算公式 数值稳定的算法：计算过程中舍入误差对计算结果影响不大的算法。 绝对稳定的计算公式：如果第n+1步的误差en+1与第n步的误差en满足|en+1/en|≤1，则称此计算公式是绝对稳定的。 p8-9，自己看 实习：p185-187 例：计算 ? 公式一： 注意此公式精确成立 记为 则初始误差 ? ?? ? ! ! ! 怎么 回事 ?! 考察第n步的误差 我们有责任改变。 造成这种情况的是不稳定的算法 /* unstable algorithm */ 迅速积累，误差呈递增走势。 可见初始的小扰动 ? 公式二： 注意此公式与公式一在理论上等价 方法：先估计一个IN ,再反推要求的In ( n N )。 可取 取 是不是运气很好? 考察反推一步的误差： 以此类推，对 n N 有： 误差逐步递减, 这样的算法称为稳定的算法 /* stable algorithm */