[2018年最新整理]Chap5Sec2数值积分与数值微分2n.pptVIP

杨东武 ydw_1978@126.com 主要内容 高斯求积公式的定义 高斯求积公式的构造 数值微分与误差估计 插值型求积公式的最高代数精度讨论 数值微分 数值微分 数值微分 * 机电工程学院 上节课内容回顾 2. 用复合梯形求积公式计算 ，则至少应将[0,1]分为多少等份才能保证所得积分的近似值有5位有效数字? 1. 已知f (1)=1,f (3)=5,f (5)=-3,用抛物线求积公式求 ≈( )。 3. 龙贝格求积公式所用的型值节点个数至少是 个。 4. 求积分公式 的代数精度。 前面讨论求积公式构造的过程中，节点xk (k=0,1, 2, … ,n)是确定的，而且通常情况下是等距分布的。 对于一个两节点的插值求积公式，待求的求积系数个数为2。一般情况下，如果求积公式具有1次代数精度，则所有求积系数可通过方程式求解.如下例: 试确定一个具有1次代数精度的公式 试确定一个具有1次代数精度的公式 解：因为当 故要公式具有1次代数精度，则必须有 解此方程组得： 梯形公式 两个变量 三个独立方程 四个变量 三个独立方程 还需补充一个方程 如果求积节点位置可以任意选取，求解公式的代数精度是否会有所提高？ 四个变量 四个独立方程 最终得到的求积公式具有3次代数精度 归纳总结 N+1个节点的积分公式最多可达到多少次的代数精度？ 2N+2个未知变量，需要2N+2个独立方程，即满足0次到2N+1次代数精度时所组成的2N+2个方程 具有最高代数精度(2n+1次代数精度)的求积公式 称为Gauss 型求积公式。其节点称为Gauss 点。 高斯求积公式 问题难点: 非线性方程组的求解 多节点的高斯求积公式岂不是很难求得？ 高斯求积公式 思考： 由勒让德多项式构造出来的高斯求积公式对应的求积区间是多少？ 首项系数为1 以此类推，便可以得到2节点~5节点时的高斯型求积公式。 关键步骤1:变量代换,转换求积区间为[-1,1] 因此,一般只讨论在插值节点上的数值微分(求导)