[理学]第05章数值积分与数值微分.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点，相应的公式为Ga