第三章.数值分析

数值积分概述

研究对象：的数值计算方法，定积分是和的极限，数值积分就是将定积分的计算用和式近似，可表为

其中为求积系数，为求积节点，系数与被积函数无关

定义1.若求积公式对时精确成立，而对不精确成立。则称该求积公式具有次代数精确度。

特别的，我们得到插值型求积公式的代数精度至少是次。

等距节点插值型求积公式称为牛顿-柯特斯公式，表示为

对牛顿-柯特斯公式通常只用低阶公式，为提高精度可在区间上等分为个小区间。

定理1 牛顿-柯特斯求积公式的误差当为奇数,，则

当为偶数,假定，则

通常会使用复合梯形和复合辛普森公式。
误差分析可通过小区间叠加得到。

收敛性是指求积公式无限细分时可以收敛到积分，稳定性是指任给只要插值点的误差小于就有求积误差小于.

定理2 若求积公式中系数，则此求积公式是稳定的。

当时求积系数出现负值，因此公式的稳定性不能保证。

高斯求积方法

考虑带权积分
其中为上的权函数,希望找到形如

的求积公式，使它具有最高的代数精确度。由于系数及节点对共有个待定参数，根据代数精确度定义，要使公式具有次代数精确度，则方程个数，由此可得。也就是只要使当选择及可使上式达到次代数精确度。另一方面，若令,
则

而，这表明求积公式不管系数如何选择，对的代数多项式不可能精确成立。

定义4 具有最高代数精确度的插值求积公式的节点,称为高斯(Gauss)点，相应求积公式称为高斯型求积公式。
定理3 插值型求积公式的求积节点是高斯点的充分必要条件是，在上以这组节点为根的多项式与任何次数的多项式带权正交，即

定理4 高斯型求积公式的系数

高斯-勒让德求积公式
即取

高斯-切比雪夫求积公式
区间上权函数的正交多项式为切比雪夫多项式，此时的高斯型求积公式称为高斯-切比雪夫求积公式。表示为

其中为的零点，余项为
。
公式由于，计算简单，又因被积函数包含因子，
故可处理含此类因子的奇异积分。

固定部分节点的高斯求积公式
如果在求积公式中有个节点固定，个节点待定，则求积公式的代数精确度为次。最常用的情形是对区间为，权值函数时要求区间端点-1或1固定，或两个端点均为固定。前者被称为高斯-拉道(Gauss-Radau)公式，后者称为高斯-罗巴托(Gauss-Lobatto)公式。

自适应求积方法

利用了相邻次细分之间积的误差和总误差的关系。

奇异积分与振荡函数积分的计算

反常积分通常是指被积函数在有限区间上无界的积分。下面介绍几种处理方法。

1. 变量置换法
   通过变量置换将奇点消除，使反常积分转换为正常积分。例如计算
   
   为充分光滑的函数。
   在为奇点，若令，则得
   
   为正常积分。

2. 区间截断法
   设,点为奇点，则可将分为
   .

3. 用高斯型积分
   
   把原式拆出分母。

对于没有现成求积公式可用得奇异积分也可用建立高斯型求积的方法处理。

无穷区间积分

1. 变量置换
   用变换或可将无穷区间变换为，用可将区间变换为，如变换后被积函数有界则可用正常积分方法计算，若变换后仍为反常积分则可用反常积分方法处理。

2. 无穷区间截断
   将无穷区间截去“尾巴”转化为有线区间，若选，使

3. 无穷区间上的高斯型求积公式
   无穷区间，权函数的正交多项式为拉盖尔多项式，由此构造的高斯-拉盖尔求积公式为
   

震荡函数积分

如傅里叶积分当较大时就是震荡函数积分，。对这类积分用通常数值积分方法计算效果都不好，因此需要用特殊方法处理。下面介绍几种方法。

1. 在零点之间积分
   设被积函数振荡部分在上的零点为
   
   则将分解为子区间并将各子区间积分之和。在每个子区间端点上被积函数之值为零，故可采用高斯-罗巴托求积公式，这样不必增加计算量就可得到较高精度。

2. Filon方法
   把它分割2N各子区间，然后每两个子区间上用的二次插值函数近似.

二重积分计算方法

利用不同的单积分求积公式就可得到不同的重积分求积法。

高斯求积公式

积分方程数值解法

考虑积分方程

其中为积分方程的核，为自由项，为参数，均为已知，为未知函数。求积分方程的解的数值方法就是在区间的某些点上求值的近似，使误差满足精度要求。通常可用数值积分方法把方程离散化。