[2018年最新整理]CHAP6数值积分与数值微分_1_4.pptVIP

第6章 数值积分与数值微分 刘东毅 天津大学理学院数学系 第6章 数值积分与数值微分 主要目的： 讨论数值积分的基本理论与方法 代数精度的概念 数值稳定性 插值型数值积分的思想 复化求积方法的思想 变步长的求积方法 Guass 求积公式 讨论求数值微分的各种方法 主要内容: 数值积分公式及其代数精度 插值型数值积分公式 与 Newton-Cotes 公式 复化求积法 变步长的梯形公式 与 Romberg 算法 Guass 求积公式 数值微分 6.1 数值积分公式及其代数精度 2.数值积分公式的代数精度 例6.1.1. 确定下列数值积分公式 6.2 插值型数值积分公式与 Newton-Cotes 公式 插值型数值积分公式 设 f (x) 在插值节点 a ≤ x0 x1 . . . xn ≤ b 处的函数值为f (xk) (k = 0 , 1 , . . . , n ) ,作 n 次 Lagrange 插值多项式 定义6.2.1 对于数值积分公式 (6.1.2) 2. Newton-Cotes公式 n = 1~8 的Cotes系数 Newton-Cotes公式的余项 Newton-Cotes公式数值稳定性 什么是数值稳定性？ 令 利用 Newton-Cotes 公式求解, 由此引起的计算结果的绝对误差为 6.3 复化求积法 1.几个常用的复化求积法 复化 Cotes 公式 例6.3.1 利用复化梯形公式和simpson公式计算 ，为使误差不超过10-5，需要将各分为多少等份，由此可得出什么结论？ 例6.3.2 利用复化求积公式计算积分 6.4 变步长的梯形公式与Romberg算法 1.变步长的梯形公式 复化求积公式称为定步长的求积公式，它对提高精度是行之有效的。但对于给定的精度，要确定一个合适的步长往往难以办到。因此实际上一般常采用变步长的求积公式。即让步长逐次折半的过程中，反复使用复化求积公式进行计算，直到相邻两次计算结果之差的绝对值小于允许精度ε的要求时终止计算，这种方法称为变步长的求积方法。 例如: 对于积分 变步长的梯形公式的算法 2. Richardson外推算法 3. Romberg 算法 当 m = 3 时 , Romberg 算法的计算过程 如此继续下去, 计算结果如下表 如果 不成立，则 h = h/2 = 1/2 ， 继续计算 。 k k T 2k T 2k 0 1 2 3 4 5 6 7 8 9 10 0.920 735 5 0.939 793 3 0.944 513 5 0.945 690 9 0.945 985 0 0.946 059 6 0.946 076 9 0.946 081 5 0.946 082 7 0.946 083 0 0.946 083 1 从上表可看出, 将积分区间对分了10次, 求得 I 的近似值为0.9460831 (积分精确值为0.9460831. . . ), 可见收敛速度比较缓慢。 若用一个步长为 h 的函数 I1( h ) 去逼近问题 I , 设其 截断误差可表示为 为了提高逼近的精度，选取 q 为满足 的 正数, 将上式(1)中的 h 换为 qh , 则有 其中 是与h无关的常数,并且 ， (1) 由 (1) 可知I1(h)逼近I的误差为 。 (2) (2) 式减上式 , 得 式(1)两端同乘以 得 (1) (2) 则 I2(h) 逼近I 误差降为 令 其中 是与h无关的常数,则有 ， 如此继续。 一般地, 选取 q 为满足 的正数, 由此得到序列 则 Im+1 ( h ) 逼近 I 的误差由下面的定理给出。 定理 6.4.1 设 I1 ( h ) 逼近 I 的截断误差由下式给出 则 Im+1 ( h )逼近 I 的截断误差为 其中 是与 h 无关的常数。 这种利用序列{Im+1(h)} 逐步加速去逼近 I 的方法 称为Richardson外推算法 Richardson外推公式 Romberg 算法是利用变步长的梯形求积序列 外推加速来逼近积分真值的算法. 考虑积分 由复化梯形公式有 现在将 Tn 记为 T1( h ), 即 设 f (x) 在区间 [a , b] 上任意次可微, 根据