数值积分及其应用.PDFVIP

第一讲 数值积分及其应用 —— 自适应数值积分 1 误差分析 定理：设I 是定积分精确值，T 是由梯形法计算出来的近似 n 2 值，若f (x) ∈C [a, b] ，则存在ξ∈(a, b) ，使得 ( ) 2 b −a b −a h I −Tn =− f (ξ) h 12 n 定理：设I 是定积分精确值，Sn 是由抛物线法计算出来的近 4 似值，若f (x) ∈C [a, b] ，则存在ξ∈(a, b) ，使得 4 (b −a) h   (4) b −a I −Sn =−   f (ξ) h 180 2   n 注：抛物线法事实上使用了2n+1 个节点 2 自适应数值积分 取等分点的缺点 当被积函数在部分区域变化较剧烈，而其他部分变化较平缓 时，采用等分点会增加工作量，效率低下。 自适应数值积分 变化剧烈的地方取较小步长，在变化平缓的地方取较大步长， 使得在满足计算精度的前提下工作量尽可能小。 作业：试给出梯形法的递推计算公式：Tk 1 2 3 4 k （n=2 , 2 , 2 , 2 , ..., 2 , ... ，即对积分区间不断对分） 3 自适应梯形法 算法：基于梯形法的自适应数值积分 递归算法 (1) 取当前积分区间(x , x )=(a, b) ； 0 1 (2) 用梯形公式计算函数在当前积分区间上的积分近似值； (3) 判断近似值的误差，若满足要求，则该值即为该区间上的 计算结果，该区间的计算过程结束。 (4) 否则，将区间二等分：令x =(x + x )/2 ，分别将子区间 c 0 1 (x , x ) 和(x , x ) 设为当前积分区间，重复步骤(2)-(4) ，直到 0 c c 1 它们的计算误差满足要求，然后将子区间(x , x ) 和(x , x ) 0 c