牛顿科特斯公式及其积分应用.pdfVIP

数值计算实习报告 ——牛顿—柯特斯积分方法及其应用 姓名：杨银月 学号：139084154 班级：数131 牛顿—柯特斯积分算法及其应用 一、引言 ●数值积分的必要性 现实生活当中往往会遇到这样的问题：拉着一块物体在一粗糙平面上沿着一直线从一点 移动到另一点所发生的功是多少？一个不规则的平面图形的面积是多少？已知边际成本— 产量函数求在一产量下的总成本…… 通过物理和几何学以及经济学的知识容易知道这类问题是需要计算积分的，根据微积分 b 定理对于积分I a f (x)dx ，只要找到f (x) 的原函数F(x) ，便有牛顿—莱布尼茨公式： b a f (x)dx F(b) F(a) 但是，现实生活中往往只能得到一些离散的点，无法得到连续的函数f (x) ，即便是给 sin x 定了f (x) ，也不一定就是容易找到原函数的 （原函数往往非初等），比方说 (x  0) ， x 2 ex 等，故不能用N-L 公式进行积分运算。即使能求得原函数的积分有时候计算也是非常困 1 难的。例如对于被积函数f (x) 6 ，其原函数 1x 2 1 1 1 1 x x 3 1 F(x) arctan x arctan(x )  ln 2  C 3 6 x 4 3 x x 3 1 计算F(a),F(b) 仍然很困难，因此有必要研究积分的数值计算问题。 ●数值积分的基本思想  由积分中值定理 （如图1）知，在积分区间[a,b] 内存在一点 ，成立 b a f (x)dx (ba)f ( ) 图1 就是说，底为ba而高为 的矩形面积恰等于所求曲边梯形的面积 。问题在于点 的I f ( )  具体位置一般是不知道的，因而难以准确算出f ( ) 的值，称f ( ) 为区间[a,b] 上的平均高 度。这样只要对平均高度f ( ) 提供一种算法，相应地便获得一种数值求积方法。如果我们 用两端点 “高度”f (a) 与f (b) 的算术平均值作为平均高度f ( ) 的近似值，这样导出的求 积公式 b ba a f (x)dx  2 [f (a) f (b)] (11) ab 便是众所周知的梯形公式。而如果改用区间的中点c 的“高度” f