§5 在数值计算时应注意的几个问题.pptVIP

  1. 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
  2. 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  3. 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  4. 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  5. 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  6. 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  7. 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
§5 在数值计算时应注意的几个问题

§5 在数值计算时应注意的几个问题 顺序:基本运算顺序 例1 2、对于不同的算法,要注意收敛速度,讲效率 二、 数值计算中要构造和使用数值稳定的计算方法 2、防止两接近的数相减 3、设法控制误差的传播 用递推关系进行计算时必须注意误差的积累. 逆向递推公式: 三、计算过程中应十分小心处理病态的数学问题 小结 * 一、减少运算次数 —— 不仅能提高计算精度,而且能减少误差的积累 1、对同一种算法(计算方式),要选用计算量少的运算次序 例如 精度(数值稳定性). 运算方案能否控制误差的传播和积累以保证计算结果有足够的 一般标准 运算次数的多少(计算效率); 运算过程是否规律(易编程); 需要记录的中间结果的多少(储存量); 精度(数值稳定性). 运算方案能否控制误差的传播和积累以保证计算结果有足够的 运算次数的多少(计算效率); (b) 利用秦九韶算法: (a) 直接计算每一项再求和: 解: 例2 解: (a) 作矩阵和向量的乘法:略. (b) 作向量的内积和加法: 计算次数少 过程规律 乘除法: kn 乘除法: kn2 例3 计算 ln2 的近似值,要求误差小于10??. 解: ① 计算量太大; ②各项的舍入误差会损失和的有效数字 (b) 用级数 来计算 用前 9 项(即取 m = 8)计算就能达到精度要求: 舍入误差 (a) 用级数 来计算 结果不可靠,计算失败 否则,则称这个算法是数值不稳定的。 算法是数值稳定的 ——计算结果受计算过程中舍入误差影响较小时。 1、注意计算机数系运算特点 有理数的有限数集,即浮点集 例4 讨论在计算机数系中分别用公式 解: 无误差时,必相等; 有舍入误差时,可能不相等, 6.799 -1.556 6.799 -1.557 6.800 -1.556 8.355 -8.355 5.243 5.243 准确值 m? m? b a 6.799 -1.556 解: 在实数集上, 取4位有效数字近似计算: 非严格单调序列 且极限也不等于a 4位有效数字舍入运算: 1234+0.4+0.3+0.2+0.1=1234 0.4+0.3+0.2+0.1+1234=1235 若出现“溢出”应立即中断 应避免出现“大数吃小数” ①事先预防 ②事后解决 例7 解: 可防溢出 精确运算: 例6 例5 例8 求下列方程的根 解: 用 8 位浮点数 (有效数字)计算 用 4 位浮数点(有效数字)计算 两接近数相减 损失了有效数字 数值不稳定的方法 仍用4 位浮点数计算 数值稳定 的方法 减法本身完全正确 误差传播的研究十分重要 是因为求 的误差(并不大), 进行减法后导致不应忽视的后果 ||δ|| ? || x||时, 计算结果 的误差较小 准确 逆向递推公式 例9 当 n = 0,1,2,?,8 时,求积分 的近似值. 见本章 §2 例1 解: 错误的原因? 将 的误差扩大到5倍。 :递推公式 引起 A 的很大误差 例10 设有方程组 将系数舍入成 2 位浮点数, 则方程为 例11 病态问题一般要用高精度(双精度)计算或解病态问题的方法解决 以上两问题称为病态的问题,病态是问题本身固有的 线性齐次方程 数值计算中应注意的问题 考虑: P14 5. 本节例 10 编程: *

文档评论(0)

skvdnd51 + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档