- 1、本文档共34页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
于是I7* , I8* 与精确值已经面目全非。 n 精确值 In 近似值In* n 精确值 In 近似值In* 0 1 2 3 4 0.63212? 0.36787? 0.26424? 0.20727? 0.17089? 0.6321 0.3678 0.2642 0.2074 0.1704 5 6 7 8 9 0.14553? 0.12680? 0.11238? 0.10093? 0.09161? 0.1408 0.1120 0.2180 -0.7280 7.5520 算法一 In =1-nIn-1 , 代入得下表 由于计算I0有误差 不计中间再产生的舍入误差 |e(In* )|= n! |e(I0* ) | 到 I8 时 |e(I8* )|= 8!ε = 40320ε 误差扩大了4万倍, 因而该算法是不稳定的。 e(In* ) = - n e( In-1* ) 分析: In=1-nIn-1 , 可以估计出 故 算法二 如果递推式改为 In-1 =(1-In )/n 则 In-1* =(1- In* )/n 误差 e (In-1* )= - e (In* )/n n 精确值 In 近似值In* n 精确值 In 近似值In* 9 8 7 6 5 0.09161? 0.10093? 0.11238? 0.12680? 0.14553? 0.1000 0.1000 0.1125 0.1268 0.1455 4 3 2 1 0 0.17089? 0.20272? 0.26424? 0.36787? 0.63212? 0.1709 0.2073 0.2642 0.3679 0.6321 若计算 I9 有误差 ε , 由e (In-1* )= - e (In* )/n ,其传播到 I0所引起的误差仅为 故该算法是稳定的。 由 I9* = 0.1000, In-1* =(1- In* ) /n,列表如下 * * 数学学院 信息与计算科学系 第一节 误 差 定义:数值计算方法是计算数学的一个分支, 又称数值分析或计算方法, 它是研究用数字计算机求解各种数学问题的数值方法及其理论的一门学科, 是程序设计和对数值结果进行分析的依据和基础。 应用计算机解决科学技术和工程问题的步骤: (1) 提出实际问题 (2) 建立数学模型 (3) 选用数值计算方法 (4) 编程上机计算得出数据结果。 一、误差的来源 1. 模型误差: 在建立数学模型过程中, 不可能将所有因素均考虑, 必然要进行必要的简化, 这就带来了与实际问题的误差。 2.?观测误差: 测量已知参数时,数据带来的误差。 3.?截断误差: 在设计算法时,近似处理带来的误差。 函数 用泰勒多项式 近似代替时,有误差 其中 在 与 之间。这种误差就是截断误差。 例如: 4.?舍入误差: 计算机的字长是有限的, 每一步运算 均需四舍五入, 由此产出的误差。 例如:用3.14159近似代替 , 产生的误差 就是舍入误差。 二、浮点数 任何一个浮点数均可表示为 其中,r叫做这个数的基,p是阶,是一个整数,取正数,负数或零。w是尾数,由t位小数构成, 若 ,则该浮点数为规格化浮点数。 对于一个特定的机器,尾数的位数t是固定的,也称其精度有t个r进位数字。 二、误差的基本概念 1. 误差和误差限 设 是准确值x 的一个近似值,称 为近似值 的绝对误差, 简称误差. 又简记 . 误差是无法计算的 (因为准确值 x 不知道), 但可以估计出它的一个上界。即 ,称 是近似值 的绝对误差限, 简称误差限. 误差是有量纲的,可正可负。 2. 相对误差和相对误差限 实际计算中, 由于准确值 x 总是未知的, 且由于 称 为近似值 的相对误差。简记为 。 相对误差是无量纲的, 也可正可负, 它的绝对值的上界称为该近似值的相对误差限, 记作 简记为 即 是 的平方项级, 故当 较小时, 常取 三、有效数字
文档评论(0)