2020第二章误差分析.pptVIP

粗略说来， 收敛性主要是研究方法误差问 题，而稳定性则更关注舍入误差问题 。在一定 条件下，两者又可以是关联的、等价的。 由于计算物理的复杂性，要弄清楚所采用 的方法的收敛性和稳定性，往往相当困难。我 们不能等待解决了收敛性、稳定性之后再去计 算，而是要在计算机实验中进行解决。但另一 方面，我们又必须对计算方法的收敛性、稳定 性问题有所了解，以利于对计算结果进行分析。 误差危害的防止 1. 选用数值稳定的计算公式，控制舍入误差的传 播。 在计算过程中误差不会增长的计算公式称 为数值稳定的，否则就是不稳定的，为了不影响 数值计算结果的精确度与真实性，在实际中我们 应尽量选用数值稳定的计算公式，避免使用数值 不稳定的公式。 第二章 误差分析 误差 数值计算中得到的近似值与精确值之差 本章主要研究误差的来源和分类， 误差的计算和估计方法以及如何减小误 差的危害等问题。 误差的来源和分类 1. 模型误差 抽象 , 简化 例：自由落体运动忽略了空气的阻力。 2. 观测误差 在数学模型中包含的参量一般是通过 实验观测确定的，其精度依赖于仪器和人 的操作。 3. 截断误差 模型的准确解和数值方法求出的近似 解之差称为截断误差，或叫方法误差。 x x 例 : sin x ? x ? ? ? ...... , 3 ! 5 ! 3 5 如果取等式右边前两项为 sinx 的近似值，则 x x sin x ? ( x ? ) ? ? ...... 3 ! 5 ! 3 5 等式的右边就是截断误差。 例：指数函数 e x 可展开成下列幂级数形式 2 n x x x e ? 1 ? x ? ? ... ? ? ... 2 ! n ! 但在实际计算时，不可能计算无穷多项，只能截取 有限项，取 2 n x x S n ( x ) ? 1 ? x ? ? ... ? 2 ! n ! 用 S n ( x ) 作为 e x 的近似值，其截断误差为 n ? 1 x x S n ( x ) ? e ? ? ... n ? 1 ! 4. 舍入误差 计算机对数据四舍五入 例 : 用 3.1415926 作为 π 的近似值 模型误差和观测误差是客观存在的，而截 断误差和舍入误差是由计算方法引起的，因此 后两种误差是我们在计算物理课中主要研究的； 讨论误差在计算过程中的传递和对计算结果的 影响，并找出误差的界，对研究误差的渐进特 性和改进算法的近似程度具有重大的意义。 误差与有效数字 1. 绝对误差 ( error ) 和绝对误差限 设某物理量的精确值为 x * ，用某种方法计算 出的近似值为 x ，则称 e= x *- x 为 x * 的绝对误差， 简称误差。 e 显示出近似值 x 的准确程度，在同一量的不 同近似值中， e 越小， x 的准确度越高。 注意：绝对误差 e 可正可负 一般某一量 x* 的精确值是不知道的，因而 也无法求出，但是往往可以估计出 e 的范围，即 存在 ε ，使得 e ? x ? x ? ? * 成立 （ 2-1 ） ε 称为 x 的绝对误差限，简称误差限或精度； ε 越小，近似值 x 的精度越高。 另外一种写法： x ? x ? ? * 2. 相对误差 ( relative error ) 与相对误差限 相对误差定义为绝对误差与精确值之比，即 * * e x ? x x ? x e r ? x * ? x * ? x （ 2-2 ） 称为近似值的相对误差。 注意 :1. 实际运算中 x* 不可知，此处用 x 代替。 2. 相对误差是无量纲的数，通常用百分比表示， 称为百分误差。 类似于绝对误差的情况，可以估计相对误 差的范围，若存在 ε r ，使得 x * ? x e x ? ? ? r x （ 2-3 ） 称 ε r 为近似值 x 的相对误差限。 例：已知 | x-x* | ≤1 cm ，问测量长度为 100 m 物体与 10 m 物体时的相对误差分别是多少？ 100 m 物体 10 m 物体