误差分析
误差的来源与分类
- 观测误差
- 模型误差
- 截断误差:近似计算出现的误差
- 舍入误差:由于计算级只能进行有限位的运算而产生的误差
- 绝对误差: e = x ∗ − x e=x^{*}-x e=x∗−x 其中e为绝对误差
如果 ϵ \epsilon ϵ满足:
∣ e ∣ ≤ ϵ |e|\leq\epsilon ∣e∣≤ϵ, 则 ϵ \epsilon ϵ为 x x x的绝对误差限 - 相对误差:
e r = x ∗ − x x = e x e_r = \frac{x^* - x}{x} = \frac{e}{x} er=xx∗−x=xe, e r e_r er为相对误差 - 相对误差限:
ϵ r = ϵ ∣ x ∣ \epsilon_r=\frac{\epsilon}{|x|} ϵr=∣x∣ϵ, ϵ r \epsilon_r ϵr为相对误差限,同时 ∣ e r ∣ ≤ ϵ r |e_r|\le\epsilon_r ∣er∣≤ϵr
有效数字
- 有效数字与绝对误差限的关系:
∣ x ∗ − x ∣ ≤ 1 2 × 1 0 m − n |x^*-x|\leq\frac{1}{2}\times10^{m-n} ∣x∗−x∣≤21×10m−n,其中 n n n为有效数字的位数, m m m为 x ∗ x^* x∗化为标准浮点数形式下的10的 m m m次幂。
注:精确值的有效数字有无穷多位! - 有效数字和相对误差限的关系(其中 a 1 a_1 a1为 x x x化为标准浮点数下小数点后第一个非零的数):
-
已知有效位数 n n n:
ϵ r ≤ 1 2 a 1 × 1 0 − n + 1 \epsilon_r\le\frac{1}{2a_1}\times10^{-n+1} ϵr≤2a11×10−n+1 -
已知 x x x的相对误差限:
ϵ r ≤ 1 2 ( a 1 + 1 ) × 1 0 − n + 1 \epsilon_r\le\frac{1}{2(a_1+1)}\times10^{-n+1} ϵ
-