重庆大学数值分析课件第4章数值积分与数值微分.pptxVIP

第四章 数值积分与数值微分;§1 数值积分概述;?;二. 求积公式的代数精（确）度;?;特别地,若取n = m，即取m+1个节点 x0 Lxm ， 则可通过给定的m+1个节点得到上述含m+1个 未知数、m+1个方程的方程组:;?;定理： 在区间[a,b]上，对于给定的m+1个互;例：确定求积公式;;如果;则;§2 Newton Cotes 公式;由;?;16;西北工业大学 理学院 欧阳洁;(n) k;当n=2时，可计算得;2 Newton-Cotes公式的截断误差;n?1;?;(2) 几种低阶求积公式的误差;证：由;② 若f(x)在[a,b]有四阶连续导数，辛浦生公 式有下列误差估计：;证：注意辛浦生公式的代数精度至少为3。 构造次数不高于三次的多项式 H3 (x)，使之 满足下列插值条件： H3 (a) ? f (a) H3 (b) ? f (b);;③ 若f(x)在[a,b]有六阶连续导数，柯特斯 公式有下列误差估计：;三. 复化求积公式;1 复化型求积公式 (1) 复化梯形公式;?;若 f ? (x) 在[a,b]连续，设m为 f ? (x) 的最小 值，M为 f ? (x) 的最大值，则;(2) 复化辛浦生公式;;其中;;(3) 复化柯特斯公式;例：用积分;39;;四．区间逐次分半求积法;对复化梯形公式，当区间分为n等份时，由;上式可改写为;;;2 区间逐次分半的复化梯形公式 目的：利用前一次计算结果以减少重复;即;§3 Romberg求积法;西北工业大学 理学院 欧阳洁 49;类似可以证明：;序列{T(n)}，{S(n)}，{C(n)}和{R(n)} 分别称为梯形序列， Simpson 序列， Cotes序列和Romberg序列。 利用公式(1)、(2)、(3)构成计算积 分的数值算法称为Romberg算法。;二. Romberg算法的实现;§4 Gauss型求积公式;一．带权积分的插值型求积公式;的;二．高斯公式和高斯点;由;有;两条结论： ① 高斯型求积公式一定是插值型求积公 式，其系数由高斯点唯一确定。 引理：求积公式至少具有n次代数精度 的充要条件是它是插值型的。高斯型求积公 式的代数精度为2n+1，故是插值型求积公 式。 ② 高斯型求积公式是代数精度最高的求 积公式。;?;3 高斯点的选取;f ( x) ? p ( x)? n ?1 ( x) 的次数不超过2n+1。故;?;注意 f (xk ) ? p(xk )?n?1 (xk ) ? q(xk ) ? q(xk );当 fn?1 (x)为正交多项式系中的n+1次多项 式，取;确定。;4 高斯型求积公式的收敛性和稳定性;5 高斯型求积公式的余项 定理 设f(x)在[a,b]内具有2n+2阶导数，则高斯型 求积公式的余项为：;?;三. 几种特殊的高斯型求积公式;1 高斯-勒让德求积公式;;常用的高斯型求积公式;四．一般高斯型求积公式的构造;解：上面的求积公式为高斯型求积公式，故对 f (x) ? 1, x, x2 , x3 精确成立的。故有;;x2 ? 0.289949;五．奇异积分与振荡积分的计算;?;（３）用高斯型积分公式;２ 无穷区间积分的计算 无穷区间上的积分可类似于有限区间上 无界函数的积分进行处理。 （1）变量代换法 通过变量代换将无穷区间转化为有限区 间，课本给出几个常用变换。 （2）区间截断法 将无穷区间截去“尾巴”转化为有限区间。 如：对于预先给定的误差限?，若存在M0，使得;（３）用无穷区间上的高斯型积分公式 对于下述无穷区间的积分：;积分;?;西北工业大学 理学院 欧阳洁 85;§5 数值微分;为讨论方便，一般假定所给节点为等距分布。;得;2 三点公式;将两边对t求导，得;;（1）;证明(3)：;西北工业大学 理学院 欧阳洁