复化抛物线积分公式.pdfVIP

摘要 求函数在给定区间上的定积分，在微积分学中已给出了许多计算方法， 但是，在实际问题计算中，往往仅给出函数在一些离散点的值，它的解析 表达式没有明显的给出，或者，虽然给出解析表达式，但却很难求得其原 函数。这时我们可以通过数值方法求出函数积分的近似值。 在用近似值代替真实值时，遇到的问题就是近似值的代数精度是否足 够。当代数精度不足够时，很显然提高插值函数的次数是一种方法，但是 考虑到数值计算的稳定性，当次数过高时，会出现龙格现象，用增大 n 的 方法来提高数值积代数精度是不可取的。因此，提出类似于分段插值，为 了减少数值积分的误差，可以把积分区间分成若干个小区间，在每个小区 间上采用低阶数值积分公式，然后把这些小区间上的数值积分结果加起来 作为函数在整个区间上的近似值，这个就是复化数值积分的思想。 本实验针对在每个小区间上利用抛物线积分公式，即阶数为 2，进行 实验。 关键词：龙格现象 复化数值积分 代数精度 复化抛物线积分公式 1、实验目的 1) 通过本次实验体会并学习复化抛物线积分公式的优点。 2) 通过对复化抛物线积分公式进行编程实现，提高自己的编程能力。 3) 用实验报告的形式展现，提高自己在写论文方面的能力。 2、算法流程 已知定积分的抛物线积分公式及其误差为 ‒ + ∫ () ≈ () + 4( ) + () 6 [ 2 ] { 2 [] =‒ ( ‒ )54 (), ∈ (,) 2880 根据数学知识，我们知道积分区间可划分，且不改变积分值，即如下所示： ∫ () = ∑∫ () a = 1 ‒ 1 针对上式，在每一个小区间上利用抛物线积分公式有 ℎ ℎ5 ( ) ( ) 4 ∫ () ≈ + 4 1 + ‒ (ℰ ) ‒ 1 6[ k ‒ 1 ( ‒ 2) k ] 2880 得到 ∫ () = ∑ ∫ () a = 1 ‒ 1 ℎ ℎ5 = ∑ ( ) + 4 1