数值计算的误差

数值计算的误差

误差的来源

1. 模型误差
2. 观测误差
3. 截断误差
4. 舍入误差

其中模型误差和观测误差是不可避免的，不属于数值分析考虑范畴

截断误差是指数学模型（数值计算方法）与精确解直接存在的误差

舍入误差是指受计算机器字长限制导致的误差，由于计算机每个单元存储位数有限，且是按二级制存储在计算机中，没有办法存储所有的数字

误差的分类

绝对误差：
$$e^{\ast }=x^{\ast }-x$$
其中 $x$ 代表准确值，$x^{\ast }$表示近似值

当 $e^{\ast }=x^{\ast }-x\ge 0$ 说明近似值大于准确指是一个强近似值

当 $e^{\ast }=x^{\ast }-x\le 0$ 说明近似值小于准确指是一个若近似值

相对误差：
$$e_r^{\ast }=\frac{e^{\ast }}{x^{\ast }}$$
相对误差更能表示一种误差成都

有效数字

若近似值与绝对值之间的误差，是某一个单位的$\frac{1}{2}$,从第一个非0数字开始到这个位数共有n个位，那么就说近似值具有n位有效数字

计算公式：

$x^{\ast }-x<\frac{1}{2}\ast 10^{m-n}$

其中$x,x^{\ast }$是用科学计数表示的，m为10的次方，n表示有效位数的个数

相对误差与有效数字之间的转换

1. 若$x^{\ast }$ 具有n位有效数字，那么相对误差满足：

   $e_r^{\ast }\le \frac{1}{2a_1}\ast 10^{-n+1}$

2. 如果相对误差满足

   $e_r^{\ast }\le \frac{1}{2a_1+1}\ast 10^{-n+1}$

   那么x至少具有n位有效数字

误差限的计算法则

误差限的四则运算
$$\begin{gathered} \epsilon (x_1^{\ast }\pm x_2^{\ast })=\epsilon (x_1^{\ast })+\epsilon (x_2^{\ast }) \\ \epsilon (x_1^{\ast }\ast x_2^{\ast })=\epsilon (x_1^{\ast })|x_2^{\ast }|+\epsilon (x_2^{\ast })|x_1^{\ast }| \\ \epsilon (x_1^{\ast }/x_2^{\ast })=\frac{\epsilon (x_1^{\ast })|x_2^{\ast }|+\epsilon (x_2^{\ast })|x_1^{\ast }|}{|x_2^{\ast }{|}^2} \end{gathered}$$
函数误差线的计算
$$\begin{gathered} \epsilon (f(x^{\ast }))\approx |f^{\prime }(x^{\ast })|\epsilon (x^{\ast }) \\ \epsilon (f(x_1,x_2,x_3...x_n))\approx \sum _{k=1}^n|\frac{\partial f(x_1^{\ast },x_2^{\ast },...,x_n^{\ast })}{\partial x_k}||\epsilon (x_k^{\ast }) \end{gathered}$$

避免误差危害的若干原则

1. 避免两个相近的数相减
2. 避免大数吃小数，可以适当调整运算顺序
3. 避免除数远小于被除数的除法，防止溢出
4. 简化公式