博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
数值积分中的梯形法则
阅读量:6036 次
发布时间:2019-06-20

本文共 2644 字,大约阅读时间需要 8 分钟。

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

 

转载于:https://www.cnblogs.com/yeluqing/archive/2012/12/30/3827934.html

你可能感兴趣的文章
centos7 安装wps 后 演示无法启动
查看>>
git简单命令
查看>>
LAMP编译部署
查看>>
XenDesktop7.6安装部署入门教程
查看>>
HashMap的工作原理及HashMap和Hashtable的区别
查看>>
GregorianCalendar日历程序
查看>>
Sublime 中运行 Shell 、Python、Lua、Groovy...等各种脚本
查看>>
【Java集合源码剖析】ArrayList源码剖析
查看>>
linux的目录结构
查看>>
这次逻辑通了,
查看>>
HTMLHelper
查看>>
快速构建Windows 8风格应用29-捕获图片与视频
查看>>
OC语言Block和协议
查看>>
使用xpath时出现noDefClass的错误(找不到某个类)
查看>>
.Net规则引擎介绍 - REngine
查看>>
CSS3 transforms 3D翻开
查看>>
利用传入的Type类型来调用范型方法的解决方案
查看>>
Top命令内存占用剖析
查看>>
转 网络IO模型:同步IO和异步IO,阻塞IO和非阻塞IO
查看>>
求带分数(蓝桥杯)
查看>>