工程数学绪论及几个关键问题.ppt

例6: 以下近似数是经四舍五入得到的 x1=0.12, x2=0.120, x3=1.73, x4=0.00073, x5=73000 x6=0.73*103 问它们分别有几位有效数字? 解： x1=0.12有2位有效数字, x2=0.120, 有3位有效数字。 x3=1.73有3位有效数字, x4=0.00073有2位有效数字。 x5=73000有5位有效数字， x6=0.73*103有2位有效数字。 例7：设计算机数系为F(10,t,L,U),将实数 x=±0.a1a2…atat+1…×10c,(a1≠0),用四舍五入法表为机 器数fl(x);求其有效数字、绝对误差限、相对误差限。 解： 相对误差为 绝对误差为： 机器数的相对误差与x无关，只与字长t有关。 定义 eps=0.5×10-t+1为机器精度。 fl(x)有t位有效数字。 注： （1）在数值计算中，尽可能多地保留近似数 的有效数字，有效数字越多， 相对误差越小， 计算结果越精确。 （2）对IEEE国际通用标准双精度数系 机器精度 Eps=2.2204e-16 最小实数 Realmin=2.2251e-308 最大实数 Realmin=1.7977e+308 (1) 加减法先对阶,后运算,再舍入 fl(x) ± fl (y)= fl (fl(x) ± fl (y)) 例9:在F(10,4,-33,33)的计算机上计算1+ 104 解: 1+ 104 =0.1000 ?101+ 0.1000 ?105 = 0.00001 ?105 + 0.1000 ?105 (对阶，靠高阶) = 0.10001 ?105 = 0.1000 ?105 = 104 (2)乘除法先运算,再舍入 fl (x) × fl (y)= fl (fl(x) × fl (y))， fl (x) ÷ fl (y)= fl (fl(x) ÷ fl (y)) (3) 不在计算机数系中的数做四舍五入处理 注：浮点数加法 不满足结合律 4、计算机中数的计算特点: 四、数值运算中的误差估计1、数值运算的绝对误差和相对误差 四、数值运算中的误差估计2、和、差、积、商的误差估计式 四、数值运算中的误差估计2、和、差、积、商的误差估计式 (1) 防止相近的两数相减 例1:当x较大时,计算 例2 计算 解：当x很小时，分子出现相近数相减，将以上算式变形 四、数值运算中的误差估计3.数值计算中值得注意的问题 (2) 防止大数吃小数 当两个绝对值相差很大的数进行加法或减法运算时,绝对值小的数有可能被绝对值大的数吃掉从而引起计算结果很不可靠. 在上式中，重新排序计算 上式= 0.2+0.4+0.4+ 23456=1+ 23456= 0.00001 ×105+ 0.23456×105=23457 例3:在F(10,5,-119,119)中，计算 23456+0.2+0.4+0.4 上式=0.23456×105+0.000002 ×105+ 0.000004 ×105+ 0.000004 ×105 =23456 (3) 防止接近零的数做除数 分母接近零的数会产生溢出错误,因而产生大的误差,此时可以用数学公式化简后再做. 四、数值运算中的误差估计3.数值计算中值得注意的问题 (4) 注意计算步骤的简化,减小运算次数 简化计算步骤是提高程序执行速度的关键，它不仅可以节省时间，还能减少舍入误差。 误差的定性分析法：即研究算法的数值稳定性。 误差的定量分析法： 定义：一个算法如果输入数据有扰动（即有误差），而计算过程中舍入误差不增长，则称此算法是数值稳定的，否则称此算法为不稳定的。 一、算法的数值稳定性 第三节 效果评价（稳定性） 此算法是不稳定的 无条件稳定与条件稳定 对任何输入数据都是稳定的算法称为无条件稳定。 对某些数据稳定，而对另一些数据不稳定的算法称为条件稳定。 例3 在F(10,4,-19,19)数系中，求解二次方程： 解法1 按求根公式 解得，x1=0.3199×103, x2=0.1000 ×10 ; 解法2 解得x1=0.3199 ×103