第4章 数值积分与数值微分.pptVIP

?Gauss求积公式的稳定性Gauss型求积公式（*）总是稳定的。证明：只需证明：因为Gauss型求积公式（*）对所有不超过2n+1次的多项式都精确成立：取是n次的Lagrange插值基函数第62页，共93页，星期日，2025年，2月5日?Gauss求积公式的收敛性则Gauss型求积公式（*）是收敛的。设证明：由Weierstrass定理知对存在m次多项式满足下证当时第63页，共93页，星期日，2025年，2月5日第64页，共93页，星期日，2025年，2月5日三、Gauss求积公式的构造根据前面的讨论，只需要取n+1次正交多项式的n+1个零点为求积节点，构造的求积公式即为Gauss求积公式?区间的转化问题任意区间经过下列变换可变为区间下面仅以Legendre多项式和Chebyshev多项式为例第65页，共93页，星期日，2025年，2月5日?Gauss-Legendre求积公式其中求积节点是n+1次Legendre多项式的零点求积系数可通过求解方程组得到，或者利用下式第66页，共93页，星期日，2025年，2月5日复合Simpson公式复合Simpson公式的几何意义小抛物面积之和近似-第30页，共93页，星期日，2025年，2月5日复合Simpson公式的余项设由介值定理余项估计式第31页，共93页，星期日，2025年，2月5日复合Simpson公式的收敛性类似地可以得到复合Cotes公式第32页，共93页，星期日，2025年，2月5日例2：分别利用复合梯形公式、复合Simpson公式计算积分的近似值，要求按复合Simpson公式计算时误差不超过。解：首先来确定步长复合Simpson公式的余项：其中第33页，共93页，星期日，2025年，2月5日本题的求法：由归纳法知第34页，共93页，星期日，2025年，2月5日解不等式得将区间8等分，分别采用复合Simpson、梯形公式01/81/43/810.9973980.9896880.9767271/25/86/87/810.9588510.9361560.9088580.8771930.841471第35页，共93页，星期日，2025年，2月5日复合梯形公式(n=8)复合Simpson公式(n=4)第36页，共93页，星期日，2025年，2月5日（1）使用复合梯形公式、Simpson公式，首先要确定步长；（2）而步长要根据余项确定，这就涉及到高阶导数的估计；（3）高阶导数的估计一般比较困难，且估计值往往偏大；（4）计算机上实现起来不方便，通常采用“事后估计法”。三、积分步长的自动选取：注意事项：?基本思想：将积分区间逐次分半?终止法则：前后两次近似值的误差小于已知精度第37页，共93页，星期日，2025年，2月5日?具体过程（以复合梯形公式为例）1、首先将区间n等分：2、再将区间2n等分，即步长减半：第38页，共93页，星期日，2025年，2月5日上述条件满足，程序终止；否则，继续分半计算。3、终止条件：由复合梯形公式的余项知变化不大时由此得到近似关系式误差控制条件第39页，共93页，星期日，2025年，2月5日收敛速度慢?对于复合Simpson公式、Cotes公式可以类似得到不足?对于复合梯形公式第40页，共93页，星期日，2025年，2月5日加速收敛应用步长逐次减半得到的复合梯形值、复合Simpson值、复合Cotes值与精确值的比较第41页，共93页，星期日，2025年，2月5日§4.4Romberg积分法/*RombergIntegrationMethod*/?Romberg积分思想由上节分析知，用复合梯形公式计算积分值的误差大约为：令由复合梯形公式知第42页，共93页，星期日，2025年，2月5日第43页，共93页，星期日，2025年，2月5日梯形加速公式：利用复合梯形公式前后两次积分近似值和，按照上式作出的线性组合得到了具有更高精度的积分值。上述公式说明：Romberg积分公式正是由此思想产生第44页，共93页