第一章误差2017年秋.pptVIP

* 数值分析是做什么用的？ 数值 分析 输入复杂问题或运算 ? ? ? ? 计算机 近似解 ? 研究对象 那些在理论上有解而又无法手工计算的数学问题 ? 例 解300阶的线性方程组 求6阶矩阵的全部特征值 主要内容 ? 数值代数 近似求解线性方程组 (直接解法, 迭代解法) 矩阵特征值的计算 ? 数值逼近： 插值法，函数逼近 ? 数值微分与数值积分 ? 微分方程近似求解 ? 非线性方程求解 常微分方程数值解法 偏微分方程数值解法 第一章 误差 §2 绝对误差、相对误差和有效数字 §3 数值计算中误差的传播 §4 数值计算中应注意的几个问题 §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*=1.24是由精确值 x 经过四舍五入得到的近似值, 求x*的绝对误差限和相对误差限. 由已知可得: 所以 ? =0.005, 解 ? 一般地, 凡是由准确值经过四舍五入得到的近似值, 其绝对误差限等于该近似值末位的半个单位. 有 位有效数字，精确到小数点后第 位 ? 有效数字 ? 若近似值 x*满足 则称 x*准确到小数点后第n位. 并把从第一个非零数字到这一位的所有数字均称为有效数字. 例： 问： 有几位有效数字？ 解： 4 3 数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 位有效数字. ? 基本运算中(? ? ? ? )的误差估计 问 §3 数值计算中误差的传播 如 例 计算 A=f (x1, x2). 如果x1, x2的近似值为 x1*, x2*, 则A的近似值为 A*=f (x1*, x2*), 用多元函数微分近似公式可以得到 ? 绝对误差 e 运算可近似看成微分运算. 由此可以得到基本运