参考资料:
数值分析教程 刘长安 西北工业大学出版社
数值计算方法 黄云清 科学出版社
数值分析简明教程 王能超 高等教育出版社 第二版
误差来源
误差来源主要分为以下几种:模型误差、观测误差、舍入误差以及截断误差。我们主要讨论一下截断误差以及舍入误差。
以下我们的x表示某一个β进制数为x=0.a1a2...as×βcx=0.a_{1}a_{2}...a_{s} \times \beta ^{c}x=0.a1a2...as×βc
截断误差
模型的准确解与数值方法的准确解之差就称为截断误差。什么意思呢?举个简单点的例子就是:
已知x>0x\gt0x>0,求exe^{x}ex时,有泰勒展开
ex=1+x+12x2+16x3+...e^{x}=1+x+\frac1 2x^{2}+\frac1 6x^{3}+...ex=1+x+21x2+61x3+...
取部分和
e(x)=1+x+12x2+16x3e(x)=1+x+\frac1 2x^{2}+\frac1 6x^{3}e(x)=1+x+21x2+61x3
作为exe^{x}ex的近似值,则误差项就为
R(x)=ex−e(x)=124e−ξx4,其中0<ξ<x.R(x)R(x)=e^{x}-e(x)=\frac1 {24}e^{-\xi}x^{4},其中0\lt\xi\lt x.R(x)R(x)=ex−e(x)=241e−ξx4,其中0<ξ<x.R(x)即为e(x)e(x)e(x)作为e−xe^{-x}e−x的近似时所产生的误差,这样的误差就是我们所谓的截断误差 .
用数学的语言描述的话就是将fl(x)fl(x)fl(x)(这是上次留下的一个坑,等会会补上)取成满足∣fl(x)∣≤∣x∣|fl(x)|\le|x|∣fl(x)∣≤∣x∣的误差最小者,即:
∣x−fl(x)∣=minf∈F,∣f∣≤∣x∣∣x−f∣|x-fl(x)|=\min \limits_{f\in F,|f|\leq|x|}|x-f|∣x