高等工程数学（张韵华，汪琥庭，宋立功）—— 第二篇：数值计算

第八章：数值积分和数值微分

1. 用复化梯形公式计算${\int }_0^1\frac{2}{1+x^2}dx$,使其截断误差不超过0.01.

   $f(x)=\frac{2}{1+x^2}$

   $\frac{d}{dx}(\frac{2}{1+x^2})=-\frac{4x}{(1+x^2{)}^2}$

   $\frac{d}{dx}(-\frac{4x}{(1+x^2{)}^2})=-\frac{4(-3x^2+1)}{(1+x^2{)}^3}$

   $f^{\prime \prime }(\eta )=-\frac{4(-3x^2+1)}{(1+x^2{)}^3}$

   $\text{ 记 }{M}_2={\max }_{a⩽x⩽b}\left|f^{\prime \prime }(x)\right|,\text{ 则有 }$，$M_2=4$，(a=0，b=1)

   $\frac{(b-a{)}^3}{12n^2}{M}_2<\epsilon$

   求得$6⩽n$，

   记 $n$ 等分的复化梯形公式为 $T_n(f)$ 或 $T(h),$ 有

$$T_n(f)=h\left(\frac{1}{2}f(a)+\sum _{i=1}^{n-1}f(a+ih)+\frac{1}{2}f(b)\right)$$

​ $n=6$，$h=\frac{1}{6}$，$$\begin{gathered} T_6(f)=\frac{1}{6}(\frac{1}{2}f(0)+f(0+1\ast h)+f(0+2\ast h)+f(0+3\ast h) \\ +f(0+4\ast h)+f(0+5\ast h)+\frac{1}{2}f(1)) \end{gathered}$$，

​ $f(x)=\frac{2}{1+x^2}$

​ $T_6(f)=1.56848$.

​ 准确值：${\int }_0^1\frac{2}{1+x^2}dx=\frac{\pi }{2}(\text{ Decimal: }1.57079\dots )$

2. 用逐次分半的复化梯形公式计算${\int }_0^1\frac{2}{1+x^2}dx$,使其截断误差不超过0.01.

复化求积公式是提高精度的一种有效方法，但在使用复化求积公式之前，必须根据复化求积公式的余项进行先验估计，以确定节点数目，从而确定合适的等分步长。因为余项表达式中包含了被积函数的导数，而估计各阶导数的最大值往往是很困难的，且估计的误差上界往往偏大。

所以实际中，常常使用“事后估计误差”的方法，通过区间的逐次分半，在步长逐次分半的过程中，反复利用复化求积公式进行计算，查看相继两次计算结果的差值是否达到要求，直到所求得的积分值满足精度要求。

用T(2n)做为积分Ｉ的近似值时，其误差大约为

计算时只需检验是否满足？若不满足，则再把区间分半进行计算，直到满足要求为止 。这样就实现了步长的自动选取。

​ $T(1)=\frac{1}{2}[f(0)+f(1)]=1.5$

​ $n=1$时

​ $T(2)=\frac{1}{2}T(n)+\frac{1-0}{2\ast 1}(f(0+\frac{1-0}{2\ast 1}(2\ast 1-1))=1.55$

​ $T(2)=\frac{1}{2^2}(f(0)+2f(\frac{1}{2})+f(1))=1.55$

​ $n=2$时

​ $T(4)=\frac{1}{2}T(2)+\frac{1}{2\ast 2}(f(\frac{1}{4})+f(\frac{3}{4}))$

​ $T(4)=\frac{1}{2^3}(f(0)+2f(\frac{1}{4})+2f(\frac{2}{4})+2f(\frac{3}{4})+f(1))=\frac{1}{8}(2+1+\frac{16}{5}+\frac{64}{17}+\frac{64}{25})=1.5656$.

​ $T(8)=\frac{1}{2^4}(f(0)+2f(\frac{1}{8})+2f(\frac{2}{8})+2f(\frac{3}{8})+2f(\frac{4}{8})+2f(\frac{5}{8})+2f(\frac{6}{8})+2f(\frac{7}{8})+f(1))=1.5695$.

​ 误差：$|(T(8)-T(4))|=0.0039$.