高斯积分公式第十一次课技术总结.ppt

二:用复化梯形公式 令h=1/8=0.125 三：用复化抛物线公式(simpson) 令h=1/8=0.125 例 题 3 四、Gauss公式 令x=(t+1)/2, 用2个节点的Gauss公式 用3个节点的Gauss公式 例 题 3 此例题的精确值为0.9460831... 由例题的各种算法可知： 对Newton-cotes公式，当n=1时只有1位有效数字，当n=2时有3位有效数字，当n=5时有7位有效数字。 对复化梯形公式有2位有效数字，对复化抛物线公式有6位有效数字。 用复化梯形公式，对积分区间[0，1]二分了11次用2049个函数值，才可得到7位准确数字。 用Gauss公式仅用了3个函数值，就得到结果。 各种方法的比较 1：梯形求积公式和抛物线求积公式是低精度的方法。复化梯形公式和抛物线求积公式，精度较高，计算较简，使用非常广泛。 2：Gauss型求积，它的节点是不规则的，所以当节点增加时，前面的计算的函数值不能被后面利用。计算过程比较麻烦，但精度高，是其他方法所不能比的。 小结 主要内容： 1:高斯型积分公式 2:高斯-勒让德求积公式 3:高斯型积分公式余项 本 节 小 结 * 第3节高斯型求积公式 一、 复化梯形公式 二、 复化辛普森求积公式 复习回顾 复化梯形公式： 求和可得 复化梯形公式 复化辛普森公式 例1 解 NIntegrate[Sin[1/x],{x,1,1.5}] 答案： 0.360811 估计误差: 注意 若使用梯形公式 解 例2 把积分区间[1，1.5]10等分，使用复化辛普森公式 估计误差 这说明使用复化梯形公式计算量比复化辛普森公式大得多 用误差估计公式不仅可以计算所求近似值的误差，还可由给定的精度估计应取多大的步长。 例3 若用复化求积公式计算积分 的近似值，要求计算结果有四位有效数字，n应取多大？ 解 所以若用复化梯形公式，n应等于41才达到精度 若用复化辛普森公式，则 高斯求积公式 数值求积公式提高精度的方法 1、增加节点的个数。 如复化梯形公式，复化辛普森公式 2、适当选取 插值型求积公式 代数精度至少是？ 代数精度至多是？ 2n+1 n 定义： 高斯求积公式 面临的问题： 哪些节点是高斯点？各节点 对应的系数是多少？ 分析： 此方程组为非线性方程组，求解十分困难。 定理： 插值型求积公式 证明： 高斯求积定理 必要性： 推出： 必要性： 充分性： 推出： 充分性： 高斯求积定理 结论： 　 　(1) 　 其中高斯点为Legendre多项式的零点 　对于一般有限区间[a,b]，用线性变换x=(a+b)/2+(b-a)t/2使它变成为[-1,1]。 Gauss - Legendre 求积公式 计算相应的系数，就可得到高斯求积公式： 高斯求积公式 具有三次代数精度 具有五次代数精度 高斯求积公式 解： 用三个节点的高斯—勒让德公式 用Gauss型求积公式计算积分近似值时，Gauss点与求积系数都是预先给出的，见p122表 例1 用高斯—勒让德求积公式计算 使其具有五次代数精度。 例1 例 题 1 解：作变换 若用n=2的Gauss-Legendre公式，则 求积分 的近似值。 例2 例 题 2 若用n=3的Gauss-Legendre公式，则 例 题 2 定理： 来源于埃尔米特插值的余项 高斯求积公式的误差 证明：… 定理： 各种做法比较如下： 1、Newton-Cotes公式 分别用不同方法计算 积分,并做比较 例 题 3 当n=1时，即用梯形公式， 当n=2时, 即用Simpson公式, 当n=3时, 当n=4时, 当n=5时, *