数值积分与数值微分1.pptVIP

数值积分与数值微分 7.5 自适应辛普森求积公式 将Sn 分半，得对分前后 和 得 自适应辛普森积分：根据被积函数在积分区间上的变化陡缓来自动确定每个子区间的积分步长，终止准则 步骤 记 H=b-a 误差 1）对分区间得 若 则 ，结束 2）若 将[a,a+H]左右两个区间分别计算误 差小于 以左半区间[a,a+H/2]为例 通过 否则再分成H/4的子区间。 数值积分与数值微分 3）重复上述过程，得划分 更准确结果 例7 计算 利用自适应辛普森积分计算程序Adapt，得6个子区间 输出项各列依次为：每个子区间 Sn对分前 对分后 计算得到SRmat= 数值积分与数值微分 quad = -1.5490 err = 2.7107e-004 7.6 高斯求积公式 7.6.1高斯型求积公式 插值型求积公式 但限制了代数精度。 【注】N-C公式，插值节点等距。简化了求积公式 数值积分与数值微分 存在n+1个节点和系数，使代数精度高于n+1 ? 最高是多少？ 高斯在1814年得到 【定理3】 的代数精度不超过n+1 证 证明对2n+1次多项式求积公式不能准确成立 取 则 且为2n+2次多项式 另一方面 【定义2】 高斯型求积公式 具有2n+1次代数精度的 插值型求积公式 高斯点 相应的求积节点 数值积分与数值微分 【定理4】 节点 为高斯点 是 的根 其中，令 与次数≤n的多项式带权正交 【注】 （1）高斯点与正交多项式的联系。即高斯点为 [a,b]上带权 的n+1次正交多项式的零点 （2）求 根据代数精度的定义或者利用插值 基函数带权积分公式 7.6.2高斯型求积公式的截断误差及稳定性、收敛性 1、高斯型求积公式的截断误差 数值积分与数值微分 2、高斯型求积公式的稳定性 【定理5】 高斯求积公式的系数 全为正，且 数值稳定性 设 ，利用高斯型求积公式 因此，高斯型求积公式是数值稳定的。 3、高斯型求积公式的收敛性 【定理6】 高斯型求积公式都收敛，即 7.6.3常用的高斯型求积公式 数值积分与数值微分 1、高斯-勒让德公式 [-1，1] 高斯点为勒让德多项式的零点 （1）单点公式 的零点 令它对 准确成立可求出 （2）两点公式 的两个零点 令它对 都准确成立，则 解得 三点高斯-勒让德求积公式 高斯-勒让德求积公式的余项 数值积分与数值微分 【注】 积分区间是[a,b] ，令 化为[-1,1] 例8 计算积分 解 (1)高斯公式 (2)牛顿-柯特斯公式 数值积分与数值微分 两点梯形公式 三点辛普森公式 【注】 高斯两点公式比三点辛普森公式精度还要高， 且少算一个函数值。 2、高斯-切比雪夫公式 [-1，1] 高斯点为切比雪夫多项式的零点 余项 3、高斯-拉盖尔公式 数值积分与数值微分 高斯点和求积系数见表 4、高斯-埃尔米特公式 高斯点和求积系数见表 数值积分与数值微分 【注】 高斯公式代数精度高，带权高斯公式可计算无界函数的广义积分、无穷限的广义积分。为满足精度可采用复化高斯公式。 7.7 二重积分 可推广到多重积分 是积分域上的连续函数 是[a,b]上的连续函数 数值积分与数值微分 记 则 7.7.1复化辛普森求积方法 对I用定步长辛普森公式，g（x）用变步长辛普森公式 其中 采用变步长为 方法计算。构成了求二重积分的复化辛普森公式 。 的辛普森求积 特别： 数值积分与数值微分 例10 计算积分 解 在各节点的函数值见表 数值积分与数值微分 准确值 本题用复化辛普森公式较精确。 7.7.2 二维数值积分的Gauss求积公式 数值积分与数值微分 7.7.3 蒙特卡罗方法 基本思想：积分可看作某个随机变量的期望值。用随 机变量的平均值来近似积分值。 1、均值估计方法 X 的概率密度是 当n很大时， 可以用 的平均值近似 为（0，1）随机数 对于任意区间[a,b],作变量代换化为[0,1]区间 数值积分与数值微分 为（0，1）随机数。 2、重积分的计算 是相互独立的（0，1）随机数。 落在区域 内的 个点记作 不属于 【注】 （1）当 时，需要先作变换 （2）Monte Carlo方法的精度和收敛速度 与问题的维数无关。 （3）缺点是计算量大，结果具有随机性，精度较低 数值积分与数值微分 3、数学建模实例-射击命中概率 数值积分与数值微分 数值积分与数值微分 7.8 数值微分 用函数值的线性组合近似函数在某点的导数值。 用插值多项式的导数近似函数的导数，用差商近 似导数。 方法 7.8.1 差商求导方法 1、一阶求导公