_第一章_误差资料.ppt

* * 数值计算方法 数值计算方法课程主要讨论如何构造求数学模型近似解的算法，讨论算法的数学原理、误差和复杂性，配合程序设计进行计算试验并分析试验结果。 与纯数学的理论方法不同，用数值计算方法所求出的结果一般不是解的精确值或者准确的解析表达式，而是所求真解的某些近似值或近似曲线。 数值分析的任务 实际问题 数学模型 数值计算方法 程序设计 上机计算 数值结果 数值计算方法讲述的基本内容 如何把数学模型归结为数值问题 如何制定快速的算法 如何估计一个给定算法的精度 分析误差在计算过程中的积累和传播 如何构造精度更高的算法 如何使算法较少的占用存储量 如何分析算法的优缺点 * 教材 《数值计算方法》，徐 涛 编著，吉林科学技术出版社 * 主要内容 ? 数值代数 近似求解线性方程组 (直接解法, 迭代解法) 矩阵特征值的计算 ? 数值逼近： 插值法，函数逼近 ? 数值微分与数值积分 ? 微分方程近似求解 ? 非线性方程求解 常微分方程数值解法 偏微分方程数值解法 * 第一章 误差 §2 绝对误差、相对误差和有效数字 §3 数值计算中误差的传播 §1 误差的来源 ? * §1 误差的来源 ? 误差按来源可分为： ? 模型误差 ? 观测误差 ? 截断误差 ? 舍入误差 * 实际计算时只能截取有限项代数和计算，如取前5项有： 这里产生误差 (记作R5 )截断误差 例如，利用 ln(x+1) 的Taylor公式计算 ln2, * ? 舍入误差 由于计算机只能对有限位数进行 原则保留有限位，这时产生的误差称为舍入误差。 等都要按舍入 运算，在运算中像 在数值计算方法中，均假定数学模型是准确的，因而不考虑模型误差和观测误差，只讨论截断误差和舍入误差对计算结果的影响. * ? 设x* 是准确值x 的一个近似值，记 e=x?x* 称 e为近似值 x* 的绝对误差，简称误差. 绝对误差一般很难准确计算, 但可以估计上界. ? 绝对误差 e 0 不唯一，当然 e 越小越具有参考价值. 则称 ?为近似值 x* 的绝对误差限，简称误差限. ? 若? 满足 §2 绝对误差、相对误差和有效数字 * 例 用毫米刻度的米尺测量一长度 x, 如读出的长度是 x*=765 mm, 由于误差限是 0.5 mm, 故准确值 ? 精确值x , 近似值 x* 和误差限 ? 之间满足： 通常记为 * ? 绝对误差有时并不能完全地反映近似值的好坏，如测量 100 m 和 10 m 两个长度，若它们的绝对误差都是 1 cm，显然前者的测量结果比后者的准确. ? 因此，决定一个量的近似值的精确度，除了要看绝对误差外，还必须考虑该量本身的大小. * 称 er 为近似值 x* 的相对误差. ? 记 由于 x 未知，实际使用时总是将 x* 的相对误差取为 ? 相对误差 称为近似值x*的相对误差限. ? * 有 位有效数字，精确到小数点后第 位 ? 有效数字 ? 若近似值 x*满足 则称 x*准确到小数点后第n位. 并把从第一个非零数字到这一位的所有数字均称为有效数字. 例： 问： 有几位有效数字？ 解： 4 3 * 例 已知下列近似值的绝对误差限都是0.005, 问 它们具有几位有效数字? a=12.175, b=－0.10, c=0.1, d=0.0032 由于0.005＝0.5×10－2, 解 所以 a 有4位有效数字1, 2, 1,7; b 有2位有效数字1, 0; c 有1位有效数字1; d 没有有效数字. * 数x*总可以写成如下形式 x* 作为x的近似值, 具有n位有效数字当且仅当 其中m是整数, ai是0到9中的一个数字, 由此可见, 近似值的有效数字越多, 其绝对误差越小. ? 有效数字的另一等价定义 * 故取 n=6，即取 6 位有效数字. 此时 x*=1.41421. 解 则近似值x*可写为 由于 令 例 为了使 的近似值的绝对误差小于10－5，问应取几位有效数字？ * ? 相对误差限与有效数字之间的关系. ? 有效数字 ? 相对误差限 已知 x* =?0.a1a2…an×10m有 n 位有效数字，则其相对误差限为 * ? 相对误差限 ? 有效数字 已知 x* 的相对误差限可写为 则 可见 x* 至少有 n 位有效数字. * 例 计算 A=f (x1, x2). 如果x1, x2的近似值为 x1*, x2*, 则A的近似值为 A*=f (x1*, x2*), 用多元函数微分近似公式可以得到 ? 绝对