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