• 数值积分中的梯形法则



    \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*}

  • 相关阅读:
    godaddy 亚太机房 更换 美国机房 全过程(图)
    博客园设置访问密码
    GoDaddy Linux主机支持机房的更换
    今天电信宽代终于装上光纤了,升级或安装光纤需购光猫,可以自购。我来扫盲一下
    我来科普一下为毛很多人升级了20M的电信光纤宽带反而感觉速度更卡了
    百度浏览器使用率统计
    hdu 1281
    C#基于SMTP协议和SOCKET通信,实现邮件内容和附件的发送,并可隐藏收件人
    如何处理标注打架
    提高你的Java代码质量吧:谨慎包装类型的比较
  • 原文地址:https://www.cnblogs.com/yeluqing/p/3827934.html
Copyright © 2020-2023  润新知