数值分析笔记(一)误差及其传播

基本概念

截断误差

s i n ( x ) = x − x 3 3 ! + x 5 5 ! − x 7 7 ! + . . . sin(x) = x - \frac{x^3}{3!}+\frac{x^5}{5!}-\frac{x^7}{7!}+... sin(x)=x3!x3+5!x57!x7+...

若取 s i n ( x ) ≈ x − x 3 3 ! + x 5 5 ! − x 7 7 ! + . . . + ( − 1 ) n − 1 x 2 n − 1 ( 2 n − 1 ) ! sin(x) \approx x - \frac{x^3}{3!}+\frac{x^5}{5!}-\frac{x^7}{7!}+...+(-1)^{n-1}\frac{x^{2n-1}}{(2n-1)!} sin(x)x3!x3+5!x57!x7+...+(1)n1(2n1)!x2n1,由于截断了无穷级数自第n项起的后段而产生的截断误差。

舍入误差

把一个浮点数在计算机中表示,可能引起的误差,这样的误差叫舍入误差。

收敛性:研究截断误差是否趋近于0。

稳定性:研究原始数据,在逐步计算过程中是否可控。

绝对误差

x ∗ x^* x x x x的近似值,则 e ( x ∗ ) = x ∗ − x e(x^{*})=x^{*}-x e(x)=xx x ∗ x^* x的绝对误差。

∣ e ( x ∗ ) ∣ = ∣ x ∗ − x ∣ ≤ ϵ |e(x^*)|=|x^*-x| \leq \epsilon e(x)=xxϵ,称 ϵ \epsilon ϵ x ∗ x^* x的绝对误差限。

相对误差

e r ( x ∗ ) = x ∗ − x x ∗ e_{r}(x^{*})=\frac{x^{*}-x}{x^{*}} er(x)=xxx

同理, ϵ r ( x ∗ ) = ∣ ϵ ( x ∗ ) ∣ ∣ x ∗ ∣ \epsilon_{r}(x^{*})=\frac{|\epsilon(x^{*})|}{|x^{*}|} ϵr(x)=xϵ(x)为相对误差限。

误差的传播

利用Taylor公式,可得

e ( y ∗ ) = y ∗ − y ≈ ∑ i = 1 n f i ( x 1 ∗ , x 2 ∗ , ⋯   , x n ∗ ) e ( x i ∗ ) e(y^*)=y^*-y \approx \sum_{i=1}^{n} f_{i}\left(x_{1}^{*}, x_{2}^{*}, \cdots, x_{n}^{*}\right) e\left(x_{i}^{*}\right) e(y)=yyi=1nfi(x1,x2,,xn)e(xi)

即函数值的绝对误差近似等于自变量绝对误差的线性组合,其组合系数为相应的偏导数。

对于一元函数,有

e ( y ∗ ) ≈ f ′ ( x ∗ ) e ( x ∗ ) , e\left(y^{*}\right) \approx f^{\prime}\left(x^{*}\right) e\left(x^{*}\right), e(y)f(x)e(x),

e r ( y ∗ ) ≈ f ′ ( x ∗ ) x ∗ y ∗ e r ( x ∗ ) e_{r}\left(y^{*}\right) \approx f^{\prime}\left(x^{*}\right) \frac{x^{*}}{y^{*}} e_{r}\left(x^{*}\right) er(y)f(x)yxer(x)

  • 11
    点赞
  • 13
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值