第四章 数值积分与数值微分

第四章 数值积分与数值微分

4.1 Newton-cotes求积

将有限区间$[a,b]$分成$n$等份，取等距基点：$a\le x_0<x_1<\cdot \cdot \cdot <x_n\le b$，则可以先根据这些点求得$f(x)$的插值函数，然后对插值函数进行求导（这样会方便一点的）。

当n=1时：
$$\begin{array}{rl}f(x)& =\frac{x-x_1}{x_0-x_1}{f}_0+\frac{x-x_0}{x_1-x_0}{f}_1\\ I& =\frac{h}{2}(f_0+f_1)\\ E& =-\frac{1}{12}{h}^3{f}^{\prime \prime }\end{array}$$
当n = 2时：
$$\begin{array}{rl}f(x)& =\frac{(x-x_1)(x-x_2)}{(x_0-x_1)(x_0-x_2)}{f}_0+\frac{(x-x_0)(x-x_2)}{(x_1-x_0)(x_1-x_2)}{f}_1+\frac{(x-x_0)(x-x_1)}{(x_2-x_0)(x_2-x_1)f_2}\\ I& =\frac{h}{6}(f(x_0)+4f(x_1)+f(x_2))\\ E& =-\frac{1}{90}{h}^5{f}^{(4)}(\delta )\end{array}$$
当n= 4时：
$$\begin{array}{rl}I& =\frac{h}{90}[7f(x_0)+32f(x_1)+12f(x_2)+32f(x_3)+7f(x_4)]\\ E& =-\frac{8}{945}{h}^7{f}^{(6)}(\delta )\end{array}$$

4.2 复合积分

在区间$[a,b]$中的小区间里使用积分公式进行积分，然后得到结果

4.3 变步长的求积公式

选定了某种积分方法之后，如何选择适当的$n$，使得计算结果达到预先选定的精度？

针对于n=2的梯形公式
$$\begin{array}{rl}R(f,T_n)& =-\frac{b-a}{12}{h}^2{f}^{\prime \prime }({\eta }_n)\\ R(f,T_{2n})& =-\frac{b-a}{12}(\frac{h}{2}{)}^2{f}^{\prime \prime }({\eta }_{2n})\\ T_{2}n-T_n& =-\frac{b-a}{12}(\frac{h}{2}{)}^2[4f^{\prime \prime }({\eta }_n)-f^{\prime \prime }({\eta }_{2n})]\\ 若充分大& ,f^{\prime }(\eta )在(a,b)上连续且f^{\prime \prime }({\eta }_n)\approx f^{\prime \prime }({\eta }_{2n})则\\ T_{2n}-T_n& =3(I-T_n)\\ I-T_{2n}& =\frac{1}{3}(T_{2n}-T_n)\end{array}$$
同理对于n=3的simpson公式有
$$I-T_{2n}=\frac{1}{15}(T_{2n}-T_n)$$
n=4的cotes公式有
$$I-T_{2n}=\frac{1}{63}(T_{2n}-T_n)$$

4.4 龙贝格公式

4.5 Gauss-Legendre求积公式

定义：

求积公式具有$2n+1$次的代数精度，则称该求积公式为Gauss型求积公式。

节点的来源：$Gauss-Legendre$公式的根

$Guass-Legendre$公式：
$$\begin{array}{rl}{P}_0(x)& =1\\ P_1(x)& =x\\ P_2(x)& =\frac{1}{2}(3x^2-1)\\ & \dots \dots \\ P_n(x)& =\frac{2n-1}{n}x{P}_{n-1}(x)-\frac{n-1}{n}{P}_{n-2}(x)\end{array}$$