c语言用六种方法求定积分.docxVIP

C 语言实验报告 求定积分 班级 10 信息与计算科学一班 姓名 戴良伟 学号 21 描述问题 利用①左矩形公式，②中矩形公式，③右矩形公式 ，④梯形公式，⑤ simpson 公式，⑥Gauss 积分公式求解定积分。 分析问题 定积分 定积分的定义 定积分 就是求函 数 f ?x? 在 区 间 ?a, b?中图线下包 围 的面积。 即 y ? 0, x ? a, x ? b, y ? f ?x?所包围的面积。这个图形称为曲边梯形，特例是曲边梯形。如下图： （图 1） 设一元函数 y ? f ?x?,在区间 ?a, b?内有定义。将区间 ?a, b?分成 n 个小区 间?a, x ?,?x , x ?,?x , x ?......?x ,b?。设?x ? x ? x ，取区间 ?x 中曲线上任意一 0 ? 0 ? 1 1 2 i i i i?1 i 点记做 f ? i ，作和式：  ?n ?lim ? ? ? n ? ??? i?1  f ?? i?xi?? ?? ? 若记 λ 为这些小区间中的最长者。当 ? ? 0 时，若此和式的极限存在， 则称这个和式是函数 f ?x? 在区间 ?a, b?上的定积分。 记作： ? b f ?x ?dx a 其中称 a 为积分下限， b 为积分上限， f ?x?为被积函数， f ?x?dx 为被积式，∫ 为积分号。 之所以称其为定积分，是因为它积分后得出的值是确定的，是一个数， 而不是一个函数。 定积分的几何意义 [1] 它是介于 x 轴、函数 f(x) 的图形及两条直线 x=a ，x=b 之间的各个部分面积的代数和。在 x 轴上方的面积取正号；在 x 轴下方的面积取负号。如图 言实现定积分计算的算法 利用复合梯形公式实现定积分的计算 假设被积函数为 f ?x?，积分区间为?a, b?，把区间?a, b?等分成n 个小区间， 各个区间的长度为h ，即h ? ?b ? a?/ n ，称之为“步长”。根据定积分的定义及几 何意义，定积分就是求函数 f ?x?在区间 ?a, b?中图线下包围的面积。将积分区间 n 等分，各子区间的面积近似等于梯形的面积，面积的计算运用梯形公 式求解，再累加各区间的面积，所得的和近似等于被积函数的积分值， n 越大，所得结果越精确。以上就是利用复合梯形公式实现定积分的计算的算法思想。 复合梯形公式： ?h ? ? 2T ? f 2 n ? ?a? ? 2?n?1 i?1 f ?x i ?? f ?b?? [2] ?? ? 具体算法如下： 1：输入积分区间的端点值 a 和 b ；  算法一 2：输入区间的等分个数 n （要求 n 尽可能大，以保证程序运行结果有较高的精确度）； 3：计算步长 h ? ?b ? a?/ n ； 4：对累加和赋初值 T ? ? f a 5：计算累加和 ? f ?/ 2 ； b T ? ?n?1 i?1 f ?x ? i 6：算出积分值 T n ? T ? h ； 7：输出积分近似值T n ，完毕。 1.2.2 利用 Smpson 公式实现定积分的计算 假设被积函数为 f ?x?，积分区间为?a, b?，把区间?a, b?等分成n 个小区间， 各个区间的长度为h 。在复合梯形公式的基础上，构造出一种加速计算积分的方法。作为一种外推算法， 它在不增加计算量的前提下提高了误差的精度。 具体算法如下： 算法二 1：输入积分上限 b 和下限 a ； 2：输入区间的等分个数 n （要求 n 尽可能大，以保证程序运行结果有较高的精确度）； 3：利用辛甫生公式： S ?n?? ?4 ?T ?2n?? T ?n??/ 3[2]， 实现对定积分的求解（其中 T ?2n?，T ?n?均为梯形公式计算所得的结果，由此可见辛甫生公式是以梯形公式为基础的）； 4：算出积分值 S n ； 5：输出积分近似值S n ，完毕。 利用 Guass 公式实现定积分计算 Guass 型求积公式是构造高精度差值积分的最好方法之一。他是通过让节点和积分系数待定让函数f(x)以此取i=0,1,2. n 次多项式使其尽可能多的能够精 确成立来求出积分节点和积分系数。高斯积分的代数精度是 2n-1，而且是最高的。通常运用的是 -1---+1 的积分节点和积分系数，其他积分域是通过变换x=(b-a)t/2 +(a+b)/2 变换到-1 到 1 之间积分。 算法三 1：输入积分上限 b 和下限 a ； 2：利用 Guass 公式，求定积分 4：算出积分值 S ； n 5：输出积分近似值S n ，完毕。 程序的编写 程序一（左矩形公式） 源程序 #include #include void main() {double f(double x); /