才可能有计算精度越高，结果的误差越小-OSCARLab.ppt

* 例1. 例2：著名的错误计算案例 正确答案是：-0.8273960599… 例3. 例4. 例5：著名的错误计算案例 例1. 精度不同, 则会得到不同数量级的数; 精度相同, 但是舍入策略有异, 则“微小”的表示误差可能会导致不一样的结果, 甚至是符号的不同. 例1. 对例 2 的分析 z与x均是37位的整数. 与例1类似, z+x要么为0，要么是一个大的整数. 要么“相减相消”； 要么“相减不相消”. 例2：著名的错误计算案例 正确答案是：-0.8273960599… 先“相减不相消”后“大数吃小数”. 第3个：可能先“大数吃小数”，后“相减相消”. 先“相减相消”, 后“表示误差” 前两个：先“相减不相消”后“大数吃小数”. 前两个：先“相减不相消”后“大数吃小数”. 第3个: 先“相减相消”, 后“表示误差” 例3. 并且 精度越高， 误差越大. 例3. 例4. 因为sin(x)在2100的错数是31，所以有效数字全部错误. 并且 精度越高， 误差越大. 例5：著名的错误计算案例 因为f(x)在x0的错数是8，所以后面的8位数字是错的. 若采用双精度计算, 结果中也只有约一半的有效数字正确. 五、总结 对于一个算术表达式， “错数”是个“坎儿” 若计算精度小于等于它，则1位正确的有效数字也得不到 只有大于它，才可能有：计算精度越高，结果的误差越小 浮点运算中的应用 若其相关错数大于16, 则即使采用双精度计算, 结果中也很难得到 1 位正确的有效数字, 有时甚至连符号都不对; 若其小于0, 则不论是单精度计算, 还是双精度计算, 表示误差均不太容易影响结果的有效数字的正确性。 谢 谢! 请批评指正！ 单精度(24位)、双精度(53位)、扩展精度(64位)下标准的错误答案如下： * 影响航天嵌入式软件可信的十个核心问题分析 计算误差造成的危害 1961 年, 美国麻省理工学院气象学家洛伦兹在仿真天气预报时, 将0.506127 舍入到0.506 , 所得计算结果大相径庭! 这种“差之毫厘, 谬以千里”的现象导致他不得不发出感叹：南美洲亚马逊河流域热带雨林中的一只蝴蝶, 偶尔煽动几下翅膀, 可能在两周后引起美国德克萨斯的一场龙卷风(即“蝴蝶效应”) 1967 年, 统计学家James Longley 通过手工和几个回归软件包对1947 年至1962 年共16 年的“国民生产总值(GNP)”、“人口”等数据进行了“就业”的回归分析. 结果发现软件包给出的答案错得离谱(drastically incorrect answers) 1982 年, 温哥华证券交易所推出一项股票指数, 其初值设定为1000.000. 在经济并无衰退的22 个月后, 指数跌到了520. 正确指数是1098.892. 软件在计算时, 多次从小数点后第4 位开始截断, 只保留了3 位小数 1987 年, 英国政府发现由于软件的舍人误差导致对过去的21 个月的通胀低估了0.1%. 这意味着与通胀挂钩的养老金测算系统的计算有误, 因此不得不重新计算养老金, 并紧急给超过9 百万的客户做补偿, 总金额达上亿英镑 计算误差造成的危害 1991 年2 月25 日, 在海湾战争中, 美军爱国者导弹防御系统未能拦截一枚来犯的萨达姆导弹, 造成美军士兵28 名死亡, 98 人受伤. 此次拦截失败的原因在于实数的不精确的二进制表示所造成的舍入错误累积所致. 爱国者导弹控制系统的系统时钟采用24 位寄存器存储时钟值0.1 秒, 这种编码方式造成0.000000095 秒的误差. 由于这种微小误差的存在, 导致系统长时运行的累积误差达0.34 秒, 进而使得拦截失败 1991 年8 月23 日, 挪威的海上油气平台Sleipner 在建设的最后阶段沉没了. 原因是由于工程师在使用有限元软件NASTRAN 计算时的精度不够, 导致平台的压力被低估了47%. 最终损失达7 亿美元 1994 年10 月, 一位数论研究者发现Pentium 处理器存在除法错误. 比如4195835/3145727, 只能精确到小数点后3 位小数, 第4 位小数是错的. 此次事件导致Intel 公司召回其CPU, 损失4 亿7500 万美元 误差可控计算原理及软件Isreal 误差可控计算原理及软件Isreal 误差可控计算原理及软件Isreal 特点: 保留小数位数 不论计算模型是“良态的”，还是“病态的 ”，结果始终是正确的 误差可控计算原理及软件Isreal 下载网址: /download.htm 误差可控计算原理及软件Isreal * 华东师范大学 20