数值计算方法--数值微分与数值积分-part-.pptVIP

五. Gauss型求积公式 求积公式的最大代数精度 Gauss点与正交多项式的关系 从简单问题得到的启示（续） 推广到一般n节点的情况 定理7.4 Gauss积分点的充要条件 定理7.4 【证明】(续-1) 定理7.4 【证明】(续-2) Gauss点正好是正交多项式的根 如何求Gauss积分系数Ai Gauss求积公式的误差及稳定性 定理7.5 证明（续） Gauss积分是数值稳定性的 Gauss积分数值稳定性的证明 (续) 几种常见的Gauss 型求积公式 Gauss-Legendre 求积公式中的系数 Gauss-Legendre 求积系数的证明 Gauss-Legendre求积系数的证明(续-1) Gauss-Legendre求积系数的证明(续-2) Gauss-Legendre 积分的截断误差 任意区间上Gauss-Legendre 积分公式 Gauss-Legendre 积分的例子 (续) 高斯-切比雪夫求积公式 Gauss-Chebyshev 积分系数 Gauss-Chebyshev积分系数(续-1) Gauss-Chebyshev 积分系数 (续-2) Gauss-Chebyshev 求积公式的误差 高斯-拉盖尔求积公式 一般形式的Gauss-Laguerre积分公式 Gauss-Laguerre 求积公式系数表 高斯-厄米特求积公式 Gauss-Hermite 求积公式系数表 关于高斯型积分的几点总结 六. 振荡函数的积分 例子（续） 求振荡函数积分的方法 具体表达式（续） 例 分段积分的误差 第7章 数值微分与数值积分 5 4 3 2 n 0.019 953 24 0.393 619 32 0.945 308 72 0.081 312 84 0.804 914 09 0.295 408 98 1.181 635 90 0.866 226 93 4,5 2,3 1 3,4 1,2 2,3 1 1,2 k 1.181 488 63 ±2.020 182 87 0.986 581 00 ±0.958 572 46 0.945 308 72 0.000 000 00 1.240 225 82 ±1.650 680 12 1.059 964 48 ±0.524 647 62 1.323 931 18 ±1.224 744 87 1.181 635 90 0.000 000 00 1.461 141 18 ±0.707 106 78 将Gauss-Hermite 积分点及积分系数计算列表, 部分数据如下: 第7章 数值微分与数值积分 把所求积分化为区间[-1,1], [0,∞] 或 [-∞,+∞] 上的形式, 根据被积函数确定用何种求积公式; Gauss型积分的计算步骤 Gauss 型积分的优缺点及应用 选定Gauss积分的精度 (即积分点数n), 查表获得Gauss点{xk} 及系数 {Ak}; 按积分公式求积分的近似值. 优点: 在积分点数确定的情况下, 它的代数精度是最高的; 或者说, 在获得相同的代数精度下, 它计算函数值的点数是最少的; 缺点: Gauss点{xk} 及系数 {Ak}须查表, 无法根据误差公式确定积分点数; 当增加节点时, 须重新查表计算所有的函数值; 应用: 在大规模力学计算中, 如有限元计算中, 经常要用到 Gauss积分, 应用十分广泛. 第7章 数值微分与数值积分 例子. 在工程计算中, 经常会遇到这类积分: 其中f(x)为非振荡函数, ω为较大的正数. 如果采用Newton-Cotes积分或Gauss积分, 会引起较大的误差. 误差: 问题的提出 若用n=10 的复化梯形公式 精确解: 所以复化梯形公式连1位有效数字都无法保证. 计算积分: 第7章 数值微分与数值积分 若采用Gauss-Legendre 积分 用5个Gauss点 (n=5), 代数精度 2n-1=9 误差估计: 同样, 精度无法得到保证. 引起问题的原因 当ω较大时, 是高度振荡的函数, 很难用多项式精度来衡量. 第7章 数值微分与数值积分 基本思想 将[a,b]等分成n 个小区间, 在每个小区间上用低阶的多项式代替函数f(x), 然后求得精确的积分结果. 在[xi,xi+1]上对f(x)作线性插值. 具体表达式 记, 或者, 计算积分 计算得 (推导过程略), 第7章 数值微分与数值积分 计算积分 或者, 其中, 计算得 (推导过程略), ……………………… (1) 其中 同(1) 式. * 第7章 数值微分与数值积分 目的 求积公式: 当节点数n固定时, 选取适当的节点{xk}及系数{Ak}, 使其具有最高的代数精度. 为权函数