zoukankan      html  css  js  c++  java
  • 数值积分中的梯形法则


    \begin{align*}
    \int_{x_0}^{x_1}f(x)dx=\frac{h}{2}[f(x_0)+f(x_1)]-\frac{h^3}{12}f''(\xi)
    \end{align*},其中$h=x_1-x_0$,$x_0=a,x_1=b$.且$f(x)$在$(a,b)$二阶可导,且在$[a,b]$上一阶导函数连续.

    证明:这是梯形法则.我们对$f(x)$进行牛顿插值.设立插值点
    \begin{align*}
    x_0,x_1
    \end{align*}
    则牛顿插值多项式为
    \begin{align*}
    f(x_0)+(x-x_0)f[x_0,x_1]
    \end{align*}
    插值余项是
    \begin{align*}
    (x-x_0)(x-x_1)\frac{f''(\xi(x))}{2!}
    \end{align*}

    \begin{align*}
    f(x)=f(x_0)+(x-x_0)f[x_0,x_1]+(x-x_0)(x-x_1) \frac{f''(\xi(x))}{2!}
    \end{align*}
    于是
    \begin{align*}
    \int_a^bf(x)dx&=\int_a^b\left(f(x_0)+(x-x_0)f[x_0,x_1]+(x-x_0)(x-x_1)
    \frac{f''(\xi(x))}{2!}\right)dx\\
    &=hf(x_0)+(\frac{1}{2}h^2f[x_0,x_1])+\int_a^b(x-x_0)(x-x_1) \frac{f''(\xi(x))}{2!}dx
    \end{align*}

    \begin{align*}
    hf(x_0)+\frac{1}{2}h^2 \frac{f(x_1)-f(x_0)}{h}=\frac{1}{2}h(f(x_1)+f(x_0))
    \end{align*}
    下面我证明
    \begin{align*}
    \int_a^b(x-a)(x-b) \frac{f''(\xi(x))}{2!}dx=-\frac{h^3}{12}f''(\xi)
    \end{align*}


    易知
    \begin{align*}
    \int_a^b(x-a)(x-b)dx&=[\frac{1}{3}b^3-\frac{1}{2}(a+b)b^2+ab^2]-[\frac{1}{3}a^3-\frac{1}{2}(a+b)a^2+a^2b]\\&=(\frac{1}{2}ab^2-\frac{1}{6}b^3)-(\frac{1}{2}ba^2-\frac{1}{6}a^3)\\&=\frac{1}{2}ab(b-a)-\frac{1}{6}(b-a)(b^2+a^2+ab)\\&=\frac{1}{6}(b-a)(3ab-ab-a^2-b^2)\\&=-\frac{1}{6}(b-a)(b-a)^2\\&=-\frac{1}{6}h^3
    \end{align*}
    可见梯形法则得证(注意这里隐含地用到了加权型的积分中值定理).

    注1:另一种做法如下:

    \begin{align*}
    \int_{x_0}^{x_1}f(x)dx=\frac{h}{2}[f(x_0)+f(x_1)]-\frac{h^3}{12}f''(\xi) \end{align*}
    其中$h=x_1-x_0,x_0=a,x_1=b$.且$f(x)$在$(a,b)$二阶可导,且在$[a,b]$上
    一阶导函数连续.

    前面主体部分是简单的,关键是余项部分.我们要证明
    \begin{align*} \int_a^b(x-a)(x-b)
    f[x,x_0,x_1]dx=-\frac{h^3}{12}f''(\xi) \end{align*}

    \begin{align*}
    g(t)=\int_{x_0}^t(x-a)(x-b)f[x,x_0,x_1]dx
    \end{align*}
    则$g(x_0)=0$,根据微积分第一基本定理,
    \begin{align*}
    g'(x_0)=g'(x_1)=0
    \end{align*}
    我们设立插值点
    \begin{align*}
    x_0,x_0',x_1,x_1'
    \end{align*}
    经过这四个插值点的牛顿插值公式为

    \begin{align*}
    g(x_0)+(x-x_0)g[x_0,x_0']+(x-x_0)(x-x_0')g[x_0,x_0',x_1]+(x-x_0)(x-x_0')(x-x_1)g[x_0,x_0',x_1,x_1']
    \end{align*}


    \begin{align*}
    x_0'\to x_0,x_1'\to x_1
    \end{align*}
    可得相应的Hermite插值多项式为

    \begin{align*}
    q(x)= g(x_0)+(x-x_0)g[x_0,x_0]+(x-x_0)^2g[x_0,x_0,x_1]+(x-x_0)^2(x-x_1)g[x_0,x_0,x_1,x_1]
    \end{align*}


    \begin{align*}
    k(x)=g(x)-q(x)
    \end{align*}
    可得
    \begin{align*}
    k'(x_0)=k'(x_1)=k(x_0)=k(x_1)=0
    \end{align*}
    使用若干次Rolle定理可得
    \begin{align*}
    k'''(\xi)=0
    \end{align*}

    \begin{align*}
    g'''(\xi)=3!g[x_0,x_0,x_1,x_1]
    \end{align*}

    \begin{align*}
    g'''(\xi)=f''(\xi)
    \end{align*}
    因此
    \begin{align*}
    f''(\xi)=3!\frac{g[x_0,x_0,x_1]-g[x_0,x_1,x_1]}{x_0-x_1}
    \end{align*},

    可见,
    \begin{align*}
    f''(\xi)=\frac{12g(x_1)}{-h^3}
    \end{align*}
    可见,
    \begin{align*}
    g(x_1)=\frac{-h^3}{12}f''(\xi)
    \end{align*}

  • 相关阅读:
    gridcontrol中使用右健菜单popupMenu1
    无线路由信号增强办法
    sql server2008企业版和标准版
    使用apache设置绑定多个域名或网站
    VirtualBox 虚拟机复制
    Mysql 性能优化7【重要】sql语句的优化 慢查询
    Mysql 性能优化6【重要】 索引优化
    Mysql 性能优化7【重要】sql语句的优化 浅谈MySQL中优化sql语句查询常用的30种方法(转)
    mysql 高可用架构
    Mysql 复制工作原理
  • 原文地址:https://www.cnblogs.com/yeluqing/p/3827553.html
Copyright © 2011-2022 走看看