第5篇数值积分.ppt

第5章 数值积分 5.1 引言 利用牛顿—莱布尼兹（Newton—Leibniz）公式 (5.1) 解决函数 在 上的积分问题在理论和应用上都有重大的意义。然而，在实际问题中，往往会遇到一些困难。有些形式上较简单的函数，其原函数 不易求出或不能用初等函数表示成有限形式；有些被积函数的原函数过于复杂；而有些函数的函数值是由实验、观测等方法得出，并没有给出具体的解析表达式。这些情形说明公式(5.1)在应用上是有局限性的，因此研究定积分的数值计算问题就显得十分必要。 本章主要介绍一些常用的数值积分方法，包括梯形积分法、辛卜生积分法、变步长积分法、牛顿—柯特斯积分法、高斯积分法、龙贝格积分法以及高振荡函数的积分法。 5.2 梯形积分法 5.2.1 梯形积分法的基本思想 梯形积分法的基本思想：在积分区间 上，根据给定的插值条件 和 ，构造一个一次二项式 ，并以 的积分值近似地代替 。从几何角度而言，是以梯形面积近似地代替曲边梯形的面积。 5.2.2 梯形求积公式 依据梯形积分法的基本思想，将区间 分成 个 相等的小区间,则每个小区间的长度为 ,对每个小区间均实施如下的梯形求积： 将这些小梯形的求积值加起来，可以得到如下梯形求积公式： 其中， 5.2.3 实现梯形积分法的基本步骤 (1) 输入区间 的端点 值以及分割数 ； (2) 将区间 等分成 个小区间，每一个小区间的 长度 ； (3) 计算每一个等分点的函数值 (4) 计算 (5) 输出 的值； (6) 结束。 例5.1 使用梯形求积公式求下列定积分的值。 main() { float a,b,h,s,f[N+1]; clrscr(); printf(input a,b=); scanf(%f,%f,a,b); h=(b-a)/(float)N; gedianzhi(f,a,h); s=trapeze(f,h); printf(s=%f\n,s); } 程序运行结果： input a,b=0,1↙ s=3.140942 辛卜生积分法的基本思想：在积分区间 上，根据给定的插值条件 、 和 ，构造—个二次插值求积多项式 ，并以 的积分值近似地代替 。从几何角度而言，是用过三点的抛物线面积近似地代替积分的曲边面积。 5.3.2 辛卜生求积公式 依据辛卜生积分法的基本思想，将区间 分成 ( 必须是偶数）个相等的小区间，则每个小区间的长度为 ，在小区间 均实施如下的辛卜生求积： 将这些求积值加起来，可以得到如下辛卜生求积公式： 其中: 5.3.3 实现辛卜生积分法的基本步骤 (1) 输入区间 的端点的 值以及分割数 ； (2)将区间 等分成 个小区间，每一个小区间的长度 ； (3) 计算每一个等分点的函数值 ； (4) 计算: （计算奇数项的函数值之和） （计算偶数项的函数值之和） (5) 计算 ； (6) 输出 的值； (7) 结束。 例5.2 使用辛卜生求积公式求下列定积分的值。 float simpson(float y[],float h)