网站大量收购独家精品文档,联系QQ:2885784924

计算方法第一章误差资料.ppt

  1. 1、本文档共67页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
(5) 要选用数值稳定的计算公式 对于一个数学问题的求解往往有多种数值方法在选择数值方法时,要注意所用的数值方法不应将计算过程中难以避免的误差放大的较快,造成计算结果完全失真。 例13 计算积分 并估计误差 解 容易得到递推公式 即 为 则准确的理论递推式 实际运算的递推式 两式相减有 这就是说,若 与 的误差为 = - ,即 ,则误差的递推规律为 于是 计算 时的误差被扩大了 倍,显然算法是数值不稳定的。 如果将递推公式 变换一种形式 准确的理论递推式 实际运算的递推式 从而有 即 于是有 则这个算法的误差传递规律为 即每计算一步的误差的绝对值是上一步的十分之一,误差的传播逐步缩小,得到很好的控制,这个算法是数值稳定的 本章小结 误差在数值计算中是不可避免的,误差的传播和积累直接影响到计算结果的精度。在研究算法的同时,必须注重误差分析,使建立起来的算法科学有效。 按照误差产生的来源可分为模型误差、观测误差,截断误差、和舍入误差等。 误差的表示法有绝对误差和相对误差两种。 在表示一个近似数时,要用到有效数字的概念,这在数值计算中非常有用,有效数字是由绝对误差决定的 通常用函数的泰勒展开对误差进行估计 * * test * * test * * * * * * * 1.3 .3 有效数字 例5. 3.142作为π的近似值时有几位有效数字 解: 3.141592…= 0.3141592…× 3.142 = 0.3142× p = 1 |π-3.142 |=|0.3141592…× -0.3142×? | < 0.000041× < 0.0005= × p –n =1–n =-3 所以 n =4,具有4位有效数字 例6. 当取3.141作为?的近似值时 ??-3.141?=?0.3141592…?101 -0.3141?101? ≤0.0000592 ?101 0.005=1/2 ? 10-2 p-n=1-n=-2 所以n=3具有3位有效数字 推论 如果近似数x*误差限是某一位的半个单位, 由该位到x*的第一位非零数字一共有n位 x*就有n位有效数字,也就是说准确到该位 再如3.1416作为?的近似值时 ? ?-3.1416 ?=? 0.3141592…?101-0.31416?101 ? ≤0?101 ≤ 000005 0.5 ? 10-4 p-n=1-n=-4 所以 n=5 x*= 3.1416有5位有效数字 例7 若 是 的具有六位 有效数字的近似值,求其误差限。 解 由于 具有六位有效数字,即 近似代替 时,准确到末位,于是误差限为 若 是 的具有五位有 效数字的近似值,求其误差限。 解 由于 具有五位有效数字, 即 近似代替 时,准确到末位,于是误差限为 关于有效数字说明 ① 用四舍五入取准确值的前n位x*作为近似值,则 x*必有n位有效数字。如3.142作为 ?的近似值 有4位有效数字,而3.141为3位有效数字 ② 有效数字相同的两个近似数,绝对误差不一定 相同。例如,设x1*=12345,设x2*=12.345,两者 均有5位有效数字但绝对误差不一样 ?x- x1* ?=?x- 12345? ≤ 0.5= 1/2 ? 100 ?x- x2* ?=?x- 12.345?≤0.0005=1/2?10-3 ③ 把任何数乘以10p(p=0,?1,…)不影响有效位数 ④ 准确值具有无穷多位有效数字,如三角形面积 S=1/2ah=0.5ah 因为0.5是真值,没有误差 ?*=0,因此n??,准确值具有无穷位有效数字 1.3.4 有效数字与相对误差 定理1.1 若近似数 具有 n 位有效数字,则其相对误差限 证明 由定义知 是具有n位有效 数字的近似值,因此 相对误差限 而

文档评论(0)

四娘 + 关注
实名认证
内容提供者

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

1亿VIP精品文档

相关文档