数值计算方法第五章 数值积分与数值微分.pptVIP

三、Gauss 求积公式的构造根据前面的讨论，只需要取n+1 次正交多项式的n+1 个零点为求积节点，构造的求积公式即为Gauss 求积公式区间的转化问题任意区间经过下列变换可变为区间下面仅以Legendre 多项式和Chebyshev 多项式为例Gauss-Legendre 求积公式其中求积节点是n+1 次Legendre 多项式的零点求积系数可通过求解方程组得到，或者利用下式习题四（5）复化Simpson 公式的收敛性类似地可以得到复化Cotes 公式例2：分别利用复化梯形公式、复化Simpson 公式计算积分的近似值，要求按复化Simpson 公式计算时误差不超过。解：首先来确定步长复化Simpson 公式的余项：其中本题的求法：由归纳法知解不等式得将区间8 等分，分别采用复化Simpson 、梯形公式0.976727 0.989688 0.997398 1 3/8 1/4 1/8 0 0.841471 0.877193 0.908858 0.936156 0.958851 1 7/8 6/8 5/8 1/2 复化梯形公式(n=8) 复化Simpson 公式(n=4) （1）使用复化梯形公式、Simpson 公式，首先要确定步长；（2）而步长要根据余项确定，这就涉及到高阶导数的估计；（3）高阶导数的估计一般比较困难，且估计值往往偏大；（4）计算机上实现起来不方便，通常采用“事后估计法”。三、积分步长的自动选取：注意事项：基本思想：将积分区间逐次分半终止法则：前后两次近似值的误差小于已知精度具体过程（以复化梯形公式为例）1、首先将区间n 等分：2、再将区间2n 等分，即步长减半：上述条件满足，程序终止；否则，继续分半计算。3、终止条件：由复化梯形公式的余项知变化不大时由此得到近似关系式误差控制条件收敛速度慢对于复化Simpson 公式、Cotes 公式可以类似得到不足对于复化梯形公式加速收敛应用步长逐次减半得到的复化梯形值、复化Simpson 值、复化Cotes 值与精确值的比较§5.4 Romberg 积分法/*Romberg Integration Method */ Romberg 积分思想由上节分析知，用复化梯形公式计算积分值的误差大约为：令由复化梯形公式知梯形加速公式：利用复化梯形公式前后两次积分近似值和，按照上式作出的线性组合得到了具有更高精度的积分值。上述公式说明：Romberg 积分公式正是由此思想产生Romberg 值序列Simpson 加速公式：Cotes 加速公式：类似于梯形加速公式的处理方法，得到：通过上述3个积分值序列求积分近似值的方法，称之为Romberg 积分法。4个积分值序列：梯形值序列Simpson 值序列Romberg 值序列Cotes 值序列Romberg 积分法的一般公式其中Romberg 积分表例3：利用Romberg 积分法式计算积分要求精确到小数点后面7位。解：根据Romberg 积分法计算得具体结果见下表§5.5 Gauss 求积公式/*Gauss Quadrature Formula */ 一、Gauss 积分问题的提法前述Newton—Cotes 求积公式中求积节点是取等距节点，求积系数计算方便，但代数精度要受到限制；为了提高代数精度，需要适当选择求积节点: ①当求积节点个数确定后，不管这些求积节点如何选取，求积公式的代数精度最高能达到多少？②具有最高代数精度的求积公式中求积节点如何选取？积分公式的一般形式：只需证明：对于上述插值型求积公式，存在一个2n+2 次多项式，使得求积公式不能精确成立。形如的插值型求积公式的代数精度最高不超过2n+1 次。证明：其中令因为而故求积公式不能精确成立下面讨论一般积分形式：其中为权函数构造积分公式（*）具有2n+1 次代数精度。其中求积节点求积系数仅与求积节点有关如果一组节点，使得上述插值型求积公式具有2n+1 次代数精度，则称该组节点为Gauss 点，相应的公式为Gauss 型求积公式。求积公式精确成立：由代数精度定义，当时，2n+2 个未知数，2n+2 个方程的非线性方程组二、Gauss 求积公式的性质Gauss 求积公式存在的条件插值型求积公式（*）的节点是Gauss 点的充要条件是以这些节点为零点的多项式与任何不超过n次的多项式带权正交: 证明：必要性设则因为是Gauss 点充分性对于其中即求积公式（*）对一切不超过2n+1 次的多项式精确成立所以节点是Gauss 点上述定理表明：上带权的n+1 次正交多项式的零点就是求积公式（*）的Gauss 点Gauss 求积公式中求积系数的求法由代数精度定义，得到n+1 阶线性方程组：设已知Gauss 点或者Gauss 求积公式的余项证明：设是满足下列条件的Hermite