数值分析 第一章-误差.ppt

  1. 1、本文档共52页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
有 所以 于是 若取 由此得到 该方法对 的舍入误差 传递情况是 递推可得 所以 误差传递逐步缩小,实际计算求得 它是 的具有8位有效数字的近似值。 如果第 步的误差 与第 步的误差 满足 则称计算公式是绝对稳定的。 可以看出: 可能比 大得多。 当 与 相近, 就可能很大,从而导 致计算结果的有效数字位数的减少。 三、避免两接近的近似数相减 由两近似值差 的相对误差关系式 例3 、试计算 (准确值为0.000000033921915…) 解、 可见,计算结果不可靠。若用 得到一个精度很高的近似值。 四、避免“大数除以小数” 由二元函数的误差传播规律式知 可知,当 相对 小时, 会很大。 由于计算机采用浮点制,在数值运算中,如果数据的数量级相差很大,如不注意运算次序,就可能出现大数“吃掉”小数的现象。 五、 防止大数“吃掉”小数 例4 、求一元二次方程实根 有两个互异实根 如果用求根公式 采用八位十进制在计算机上计算, , 其中 明显错误。 表明 与1两数相加,大数 “吃掉”了小数1。 于是所得两根为 求出 ,再利用韦达定理 若利用 求出 ,就可以保证一元二次方程 两个根都是可靠的。 一、向量的范数 ,其模为 我们用其度量向量的“大小”。 §3 向量范数与矩阵范数 2、(齐次性).任意 ,有 3、(三角不等式). 将向量模的概念加以推广,便得到向量范数概念。 实质上,向量的模 满足三个条件: 1、(非负性). 当且仅当 定义1 定义映射 ,若满足条件: ② ③ 当且仅当 则称其为 上的一种范数。 1—范数 2—范数 范数 最常用的如下三种向量范数: 下面只证 是向量范数: 证明:(1)由定义有 (3)易见, (2) 由Cauchy-Schwarz不等式 向量的2—范数也称为Euclid范数。 其实,向量的1—范数、2—范数、 范数都是p—范数 的特例。 易证: 中范 定理(范数的相容性). 设 为 数,则存在 使得 例5、 设 ,求 解:由定义 定义2 定义映射 ,若满足: 1、 当且仅当 2、 3、 则称其为 上一种的矩阵范数。 二、矩阵的范数 在实际应用问题中,需要矩阵空间中范数与向量空间中范数满足下面不等式(范数的相容性): 称其为由向量范数诱导的矩阵空间的算子范数。 对于 ,令 定理1 在 中一种给定向量范数 为矩阵 的范数,记为 。 (1)当 时, ;当A≠0时,必有 有 (2)对任一数 定理证明 证毕。 (3)对任意的矩阵 此结果显然也适用于Y=0的情形。 证、 证明向量范数与其诱导的矩阵空间的算子范数相容。 1―范数(列范数) 2―范数(谱范数) ∞―范数(行范数) 其中, 为 的特征值。 常用的三种矩阵算子范数为: 例6、 设 ,求 。 解、 定理 2. 设 ,若存在 中算子范数 使得 ,则 非奇异,且有 证明. 若 为奇异矩阵,则存在 使得 显然矛盾。 练习:已知矩阵A和向量X,求 其中, 数 值 分 析 第一章 误差 主讲老师: 汪达成 数值分析是计算数学核心课程,它研究高等数学和线性代数中基本数学问题的数值解法,以及在求解过程中出现的收敛性、数值稳定性和误差估计等问题。 §1 误差 数值解:满足一定精度的近似解。 精度:我们用误差或近似数的有效数字刻划。 一、 误差的产生和分类 实际 问题 验证结果 编程求解 选择算法 收集数据 建立模型 模型误差 观测误差 截断误差 舍入误差 (2)观测误差:数学模型中一些物理量的观测值,如:电压、温度、长度等,不可避免会带来误差,称为观测误差。 (1)模型误差:由实际问题转化为数学问题,在建立数学模型时,数学模型与实际问题之间出现的误差称为模型误差。 (3)截断误差:计算机在求解数学模型时选用数值计算方法,由此方法产生的误差称为截断误差。 (4)舍入误差:由于计算机字长有限,只能对有限位进行运算,因而往往进行四舍五入,这样产生的

文档评论(0)

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

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

1亿VIP精品文档

相关文档