数值积分的牛顿——科茨求积.doc

数值积分的牛顿——科茨求积 [摘要]：在实际生活中我们常遇到数值积分的求积问题，虽然我们也学过求数值积分的一些方法，但是由于用插值多项式近似表达函数f(x)时存在截断误差，即有插值余项，因此插值型求积公式也有相应的余项。存在求函数f(x)在区间[a,b]上的定积分以及在给定点上的值的数值方法，为了克服求的原函数可能遇到的困难和便于计算，我们利用牛顿——科茨来计算。 其中还推导它的两种特殊形式——梯形求积公式和辛普森求积公式，并对这三种求积公式（梯形公式、辛普森公式和柯茨公式）进行了分析和比较。现在要对数值积分进行求积需要运用matlab对梯形求积公式、辛普森求积公式和牛顿—柯茨公式进行编程实现，程序简洁、直观、求解速度快并且方法实用性强。 [关键字]：插值积分、梯形求积公式、辛普森求积公式、牛顿——科茨公式 梯形求积公式 梯形求积公式即使当n=1时，过a,b两点，做直线： 用代替，得 用梯形面积近似替代曲面梯形的面积，所以（1）式叫做梯形求积公式。 辛普森求积公式 辛普森求积公式即是当n=2时，把区间2等分即是过a、b和 三点，做抛物线： 用代替，则可求得 式（2）就叫做辛普森（Simpson）公式。从几何意义上来看，因为辛普森公式是用抛物线围成的曲边梯形来近似代替所围成的曲边梯形面积，所以辛普森公式也叫做抛物线求积公式。 牛顿-科次（Newton-Cotes）公式 牛顿-科次（Newton-Cotes）公式即把区间[a,b]n等分，其分点为，过这n+1节点，可以构造一个n次差值多项式： 其中，用代替被积函数则有 公式（3）叫做牛顿-科次（Newton-Cotes）公式，使用牛顿-科次（Newton-Cotes）公式的关键是计算系数，用变量替换，于是 这时是依赖于函数和区间[a,b]的常数，可以事先计算出来，叫做牛顿-科茨系数。 利用式（3）和式（5）得到牛顿-科次系数后，便可以写出相应的牛顿-科次公式。当n=1时，牛顿-科次公式为 即之前讨论过的梯形求积公式，当n=2时，牛顿-科次公式为 即辛普森公式。 所以我们可以看出梯形求积公式和辛普森求积公式是牛顿-科次公式的特例。当n=4时，牛顿-科次公式为 其中xi=a+k*h(k=0,1,2,3,4)、，式（8）也称为科次公式。 例1 试分别用梯形求积公式、辛普森求积公式和科次求积公式计算定积分。并用Mat lab编写程序，求解积分要求给出实验结果。 解：由梯形求积公式可得 由辛普森求积公式得到 利用科次求积公式，由n=4可得到x0=0，x1=0.2500，x2=0.5000，x3=0.7500，x4=1.0000则 原积分的准确值为3.1421，可见三个求积公式得到的数值解与准确值之间的误差是逐渐减少的。 用Mat lab编写程序如下： disp(已知 y=4/(1+x.^2), ) disp( 用辛甫生公式、梯形公式和柯次公式分别求积分值PI=int(y,x,0,1)) %辛甫生公式 f(x)在[a,b]上的积分为：(b-a)/6*( f(a)+f(b)+4*f( (a+b)/2) ); %梯形公式 f(x)在[a,b]上的积分为：(b-a)/2*( f(a)+f(b) ); %柯次公式 f(x) 在[a,b]上的积分为：(b-a)/90*[7*f(x0)+32*f(x1)+12*f(x2)+32*f(x3)+ 7*f(x4)] fprintf(\n方法一： 辛甫生公式计算\n) b=1; a=0; x0=0;x1=0.2500;x2=0.5000;x3=0.7500;x4=1.0000; PI=(b-a)/6*( f(a)+f(b)+4*f( (a+b)/2) ) fprintf(\n方法二： 梯形公式计算\n) PI=(b-a)/2*( f(a)+f(b) ) fprintf(\n方法三： 柯次公式计算\n) PI=(b-a)/90*[7*f(x0)+32*f(x1)+12*f(x2)+32*f(x3)+ 7*f(x4)] 对应程序：（保存为文件名f.m文件） function y=f(x) y=4/(1+x.^2); 实验结果： 方法一： 辛普森公式计算 PI = 3．1333 方法二： 梯形公式计算 PI = 3 方法三： 柯次公式计算 P