高斯积分的基本理论

关志金

已于 2023-11-27 20:06:20 修改

阅读量4.4k

点赞数 2

分类专栏：笔记之数值算法文章标签： 1024程序员节数学

于 2021-10-24 21:32:19 首次发布

本文链接：https://blog.csdn.net/Jinzhi_Guan/article/details/120941476

版权

笔记之数值算法专栏收录该内容

6 篇文章 2 订阅

订阅专栏

高斯积分

高斯积分：一类具有最高的代数精度的插值型求积

数值积分

基本思想

数值积分的特点：将被积函数在某些节点上的函数值加权求和并以该和值作为积分的近似值：
$\int_{a}^{b}f(x) dx \approx \sum\limits_{k=0}^{n} A_{k} f(x_{k}),\ a\leq x_{0}<x_{0}<\cdots<x_{n}\leq b$
从而将求函数积分问题归结为函数值的计算问题，避开了 $N e wt o n - L e ibni z$ 公式需要求原函数的困难。

由积分的几何意义，我们可以得到一些简单的数值积分公式：

梯形公式：
$\int_{a}^{b} f(x) dx \approx \frac{b-a}{2}[f(a)+f(b)]$
中矩形公式：
$\int_{a}^{b} f(x) dx \approx (b-a)f\left(\frac{a+b}{2}\right)$
辛普森(Simpson)公式
$\int_{a}^{b} f(x) dx \approx \frac{b-a}{6}\left[f(a)+4 f\left(\frac{a+b}{2}\right)+f(b)\right]$

代数精度

数值求积公式是近似的，该如何刻画求积公式的精确程度呢？为此引入 “代数精度” 的概念，借助多项式来刻画求积公式的精确程度。

定义： 如果求积公式对所有次数不超过 $m$ 的多项式是精确的，但对 $m + 1$ 次多项式不精确，则称该求积公式具有 $m$ 次代数精度。

可以验证，梯形公式、中矩形公式具有 $1$ 次代数精度，辛普森公式具有 $3$ 次代数精度。自然的，欲使求积公式具有 $m$ 次代数精度，只要令它对于 $\cdots, x^{m}$ 都能精确成立，即

$\begin{cases} \sum\limits_{k=0}^{n} A_{k} = b - a \\ \sum\limits_{k=0}^{n} A_{k}x_{k} = \dfrac{1}{2}(b^{2} - a^{2}) \\ \qquad \cdots \cdots\\ \sum\limits_{k=0}^{n} A_{k}x_{k}^{m} = \dfrac{1}{m+1}(b^{m+1} - a^{m+1}) \\ \end{cases}$
由于其系数矩阵为范德蒙德 $(Van d er m o n d e)$ 矩阵，所以当 $x_{i},i=0,1,\cdots,n$ 为 $n + 1$ 个互异节点时，总存在 $n + 1$ 个相应的求积系数 ${A_{k}\}$ ，使得求积公式 $(1)$ 至少具有 $n$ 次代数精度。(待定系数法求积分公式)

正交多项式

相关理论

定义： 区间 $[a, b]$ 上的非负函数 $\rho(x)$ ，若它满足：

$(1)$ $\int_{a}^{b}|x|^{n} \rho(x) dx$ 对一切非负整数 $n$ 可积且有限；

$(2)$ 假设对某个非负的连续函数 $g (x)$ ， $\int_{a}^{b}\rho(x)g(x)dx = 0$ ，则在区间 $[a, b]$ 上函数 $g (x) = 0$ ，

则称 $\rho(x)$ 为权函数。

定义： 定义在 $[a, b]$ 上的函数系 ${g_{l}(x)\}_{l = 0}^{n}$ ，如果满足：

$(1)$ $g_{l}(x) = \sum\limits_{k=0}^{l}\alpha_{k}x^{k}$ 恰为 $l$ 次多项式，即 $\alpha_{l}\not=0$ ；

$(2)$ $(g_{i},g_{j}) = \begin{cases}0,& i\not= j \\ \int_{a}^{b}\rho(x)g_{i}^{2}(x)dx > 0, &i=j \end{cases}$ ，

则函数系 ${g_{l}(x)\}_{l = 0}^{n}$ 称为 $P_n$ 的带权 $\rho(x)$ 正交基 ( $g_{l}(x)$ 称为 $[a, b]$ 上的带权 $l$ 次正交多项式 )。特别地，若 $(g_{i},g_{j}) = 1, i = 0,1,\cdots,n$ ，则称 ${g_{l}(x)\}_{l = 0}^{n}$ 为 $[a, b]$ 上 $P_n$ 的正规正交基。

由该定义易知，对任给的 $k$ 次多项式 $p_{k}(x)$ 有
$p_{k}(x),g_{l}) = 0,\ l>k$
记 $g_{l}^{*}(x) = g_{l}(x)/\alpha_{l}, l =0,1,\cdots,n$ ，则 ${g_{l}^{*}(x)\}_{l = 0}^{n}$ 构成了 $P_{n}$ 的首项系数为 $1$ 的带权正交多项式基。类似于 $G r am - S c hmi t$ 正交化变换，可以得到带权正交基的递推公式：
$\begin{cases} g_{0}^{*}(x) = 1,\ g_{1}^{*}(x) = x - (xg_{0}^{*},g_{0}^{*})/(g_{0}^{*},g_{0}^{*}),\\ g_{k+1}^{*}(x) = (x-\beta_{k})g_{k}^{*}(x) - \alpha_{k}g_{k-1}^{*}(x),\ k = 1,2,\cdots,n-1, \end{cases}$
其中，
$\beta_{k} = \dfrac{(xg_{k}^{*},g_{k}^{*})}{(g_{k}^{*},g_{k}^{*})}, \ \alpha_{k} = \dfrac{(g_{k}^{*},g_{k}^{*})}{(g_{k-1}^{*},g_{k-1}^{*})}.$
利用公式 $(6)$ 这一性质，可得如下定理：

定理： 带权正交多项式 $g_{l}(x)(l\geq1)$ 有个互异的是根，并且全部位于区间 $[a, b]$ 内。

证明： 首先由 $g_{l},g_{0})=0$ ，则有 $g_{l}^{*}$ 在区间 $[a, b]$ 内有零点。在区间 $[a, b]$ 里，如果有 $m (m < l)$ 个不同的零点 $x_{1},x_{2},\cdots,x_{m}$ ，而且是奇次零点，那么一定有多项式 $g_{l}^{*}(x-x_{1})(x-x_{2})\cdots(x-x_{m})$ 在区间 $[a, b]$ 上保号，这就意味着
$\int_{a}^{b}\rho(x)g_{l}^{*}(x-x_{1})(x-x_{2})\cdots(x-x_{m})dx > 0$
这与公式 $(6)$ 矛盾，所以必然有 $m = l$ 。如果其中有偶次重根，在 $g_{l}(x-x_{1})(x-x_{2})\cdots(x-x_{m})$ 中除去 $x$ 与该根之差，一样可以构造式 $(9)$ ，这与公式 $(6)$ 矛盾，于是 $g_{l}^{*}$ 不存在偶次根，从而再次得出 $m = l$ 。命题得证。

结合递推公式和数学归纳法，容易证明如下结论：

定理： 设 ${g_{l}(x)\}_{l=0}^{n}$ 为带权正交多项式系。对 $l\geq 1$ ，多项式 $g_{l}(x)$ 和 $g_{l+1}(x)$ 的零点必交错，即若 $\xi_{1}<\xi_{2}<\cdots<\xi_{l},\ \eta_{1}<\eta_{2}<\cdots<\eta_{l+1}$ ，分别为 $g_{l}(x)$ 和 $g_{l+1}(x)$ 的零点，则
$a<\eta_{1}<\xi_{1}<\eta_{2}<\xi_{2}<\cdots<\xi_{l}<\eta_{l+1}<b.$

证明： 略

一些例子

勒让德 ( $L e g e n d re$ ) 多项式：

设 $\rho(x) = 1$ ， ${l_{i}(x)\}_{0}^{n}$ 为 $[- 1, 1]$ 上 $P_n$ 的正规正交基，则称 $l_{i}(x)$ 为 $L e g e n d re$ 多项式。
$l_{k}(x) = \frac{1}{2^{k}k!}\cdot\sqrt{\frac{2k+1}{2}}\cdot\frac{d^{k}}{dx^{k}}[(x^{2}-1)^{k}],\ k=0,1,\cdots,n$

$l_{k+1}(x) = \frac{\sqrt{(k+1)(2k+3)}}{k+1}xl_{k}(x)-\frac{k}{k+1}\sqrt{\frac{2k+3}{2k-1}}l_{k-1}(x),\ k=1,2,\cdots$

$C h e b ys h e v$ 多项式：

区间 $[- 1, 1]$ 上带权 $\rho(x) = (1-x^{2})^{-1/2}$ 的 $j$ 次正交多项式： $T_{j}(x) = \cos(j\arccos x),\ j=0,1,\cdots$

第二类 $C h e b ys h e v$ 多项式：

区间 $[- 1, 1]$ 上带权 $\rho(x) = (1-x^{2})^{1/2}$ 的 $n$ 次正交多项式： $U_{n}(x) = \dfrac{\sin((1+n)\arccos x)}{\sqrt{1-x^{2}}},\ n=0,1,\cdots$

拉盖尔 ( $L a gu erre$ ) 多项式：

区间 $[0,+\infty)$ 上带权 $\rho(x) = e^{-x}$ 的 $n$ 次正交多项式： $L_n(x) = e^{x}\dfrac{d^{n}}{dx^{n}}(x^{n}\cdot e^{-x}),\ n=0,1,\cdots$

埃尔米特 ( $Her mi t e$ ) 多项式：

区间 $(-\infty,+\infty)$ 上带权 $\rho(x) = e^{-x^{2}}$ $n$ 次正交多项式： $H_n(x) = (-1)^{n}e^{x^{2}}\dfrac{d^n}{dx^{n}}(e^{-x^{2}}),\ n=0,1,\cdots$

Gauss型求积公式

考虑求积公式
$\int_{a}^{b}\rho(x)f(x)dx\approx \sum\limits_{k=1}^{n}A_{k}f(x_k)$
其中 $x_{i}\not=x_{j}$ (当 $i\not=j$ 时)。若 $A_{k}$ 满足：

$A_{k} = \int_{a}^{b}\rho(x)l_{k}(x)dx , \ l_{k}(x) = \prod_{\large j = 1\atop j\not=k}^{n}\dfrac{x-x_{j}}{x_{k} - x_{j}}$
则称式 $(13)$ 为插值型的。

容易证明如下定理：

定理： 形如式 $(13)$ 的求积公式至少具有 $n - 1$ 次代数精度的充分必要条件是它是插值型的。

证明：略

定义： 形如式 $(13)$ 的求积公式具有 $2 n - 1$ 次代数精度，则称它为 $n$ 点 $G a u ss$ 型求积公式，积分节点 ${x_{k}\}$ 称为 $G a u ss$ 点。

对于一般的 $n$ 直接法求 $G a u ss$ 积分公式是比较困难的，由上述定理可得， $n$ 点 $G a u ss$ 型求积公式 $(13)$ 必定是插值型的，因此只要找到 $G a u ss$ 点 ${x_{k} \}$ ，求积系数 ${A_{k} \}$ 就可以按式 $(14)$ 算得。下面从研究 Gauss 点的基本特性入手来解决 Gauss 公式的构造问题。

定理： 互异节点 $x_{k}, k=1,2,\cdots,n$ 是 $n$ 点 $G a u ss$ 求积公式的 $G a u ss$ 点的充要条件是 $\omega_{n}(x) = \prod\limits_{k=1}^{n}(x-x_{k})$ 与任何次数不超过 $n - 1$ 的多项式关于权函数 $\rho(x)$ 正交，即成立
$\int_{a}^{b}\rho(x)\omega_{n}(x)\cdot x^{j}dx = 0,\ j = 0,1,2,\cdots,n-1.$

证明：略

定理： 对于任何正整数 $n$ ， $n$ 点 $G a u ss$ 求积公式是存在的，其 $G a u ss$ 点就是求积区间上关于权函数正交多项式系 $\{\varphi_{m}(x)\}$ 中 $\varphi_{n}(x)$ 的零点。

证明：略

省略部分下次再补充。

关志金

关注

2
点赞
踩
26

收藏

觉得还不错? 一键收藏
0
评论
高斯积分的基本理论

高斯积分数值积分基本思想代数精度正交多项式相关理论一些例子Gauss型求积公式高斯积分：一类具有最高的代数精度的插值型求积数值积分基本思想数值积分的特点：将被记函数在某些节点上的函数值最为加权求和并以该和值作为积分的近似值：∫abf(x)dx≈∑k=0nAkf(xk), a≤x0<x0<⋯<xn≤b\int_{a}^{b}f(x) dx \approx \sum\limits_{k=0}^{n} A_{k} f(x_{k}),\ a\leq x_{0}<x_{0
复制链接

扫一扫