《牛顿科特斯公式及其积分应用》由会员分享，可在线阅读，更多相关《牛顿科特斯公式及其积分应用（15页珍藏版）》请在金锄头文库上搜索。
1、数值计算实习报告数值计算实习报告牛顿牛顿柯特斯积分方法及其应用柯特斯积分方法及其应用姓名：杨银月 学号：139084154 班级：数 131牛顿牛顿柯特斯积分算法及其应用柯特斯积分算法及其应用一、一、引言引言数值积分的必要性数值积分的必要性现实生活当中往往会遇到这样的问题： 拉着一块物体在一粗糙平面上沿着一直线从一点移动到另一点所发生的功是多少？一个不规则的平面图形的面积是多少？已知边际成本产量函数求在一产量下的总成本通过物理和几何学以及经济学的知识容易知道这类问题是需要计算积分的， 根据微积分定理对于积分badxxfI)(，只要找到)(xf的原函数)(xF，便有牛顿莱布尼茨公式：)()()(
2、aFbFdxxfba但是，现实生活中往往只能得到一些离散的点，无法得到连续的函数)(xf，即便是给定了)(xf，也不一定就是容易找到原函数的（原函数往往非初等） ，比方说)0(sinxxx，2xe等，故不能用 N-L 公式进行积分运算。即使能求得原函数的积分有时候计算也是非常困难的。例如对于被积函数611)(xxf，其原函数Cxxxx xxxxF1313ln341)1arctan(61arctan31)(22计算)(),(bFaF仍然很困难，因此有必要研究积分的数值计算问题。数值积分的基本思想数值积分的基本思想由积分中值定理积分中值定理（如图 1）知，在积分区间,ba内存在一点，成立)()()
3、(fabdxxfba图 1就是说， 底为ab而高为)(f的矩形面积恰等于所求曲边梯形的面积I。 问题在于点的具体位置一般是不知道的，因而难以准确算出)(f的值，称)(f为区间,ba上的平均高度。这样只要对平均高度)(f提供一种算法，相应地便获得一种数值求积方法。如果我们用两端点“高度”)(af与)(bf的算术平均值作为平均高度)(f的近似值，这样导出的求积公式)()(2)(bfafabdxxfba) 11 ( 便是众所周知的梯形公式梯形公式。而如果改用区间的中点2bac的“高度”)(cf近似地取代平均高度)(f，则又可导出所谓的“中矩形公式” 。更一般地，可在区间,ba上适当选取某些节点kx，
4、 然后用)(kxf的加权平均得到平均高度)(f的近似值，这样构造出的求积公式具有下列形式： nkkkbaxfAdxxf0)()(，)21 ( 式中kx称为求积节点；kA称为求积系数，亦称伴随节点kx的权。权kA仅仅与节点kx的选取有关，而不依赖于被积函数)(xf的具体形式。这类数值积分通常称为机械求积，其特点是将积分求值问题归结为被积函数值的计算，这就避开了 N-L 公式需要寻求原函数的困难。梯形公式即机械型求积公式1n的情况，而由公式)21 ( ，给定不同的n可以得到不同的求积公式。下面讨论不同n的情况下，求积公式的形式与性质。二、二、N-C 公式及其简单形式（公式及其简单形式（2, 1nn
5、）N-C 公式推导公式推导对于大部分的函数)(xf都是不容易求积的，但是可以用插值法的思想将)(xf用多项式插值表示，因此可以定义插值型的求积公式： nkkkbaxfAdxxf0)()(其中bakkdxxlA)(。等距节点的插值型求积公式称为牛顿柯特斯公式，推导如下：取等距节点：ihaxi，nabh,ni,.,2 , 1，令thax，得： njiinnijbaijjijbaiidtjtinniabhdtjijtdxxxxxdxxlA 00)()!( !) 1)()(，称 njiin idtjtinniabA0)()!( !) 1(为柯特斯系数， 记做)(n iC， 所以有牛顿柯特斯公式： ni
6、in ibaxfCabdxxf0)()()()(由于是多项式积分，柯特斯系数的计算不会遇到困难。当1n时，21)1( 1)1( 0 CC,这时的求积公式就是梯形公式；当2n，这时的柯特斯系数为:61)2)(1(4120)2( 0dtttC；64)2(2120)2( 1dtttC；61) 1(4120)2( 2dtttC.相应的求积公式是如下辛普森公式辛普森公式).()2(4)(6bfbafafabS) 12( 稳定性分析稳定性分析下表是柯特斯公式的系数表（n=1,2，8） （程序程序 21）n)(n iC10.50000.500020.16670.66670.166730.12500.3750
7、0.37500.125040.07780.35560.13330.35560.077850.06600.26040.17360.17360.26040.066060.04880.25710.03210.32380.03210.25710.048870.04350.20700.07660.17300.17300.07660.20700.043580.03490.2077-0.03270.3702-0.16010.3702-0.03270.20770.0349从表中可以看出当8n时，柯特斯系数)(n iC出现负值，于是有10)(0)( nin inin iCC,如果0)()(iin ifxfC，且
8、 )(iifxf，则有 niiin iniiin innfxfCfxfCfIfI0)(0)( )()()()( nin iniiin iCfxfC0)(0)()(.这表明初始数据误差将会引起计算结果误差增大，即计算不稳定，故8n时的牛顿柯特斯公式是不用的。牛顿柯特斯公式通常只用4 , 2 , 1n时的公式，下面只讨论2 , 1n时的误差。 梯形公式与辛普森公式的误差分析梯形公式与辛普森公式的误差分析代数精度代数精度：如果某个求积公式对于次数不超过m的多项式均能准确地成立，但对于1m次多项式就不准确成立，则称该求积公式具有m次代数精度。定理定理：若求积公式的代数精度为m，则求积公式的余项可表示为
9、：)()()()1(0mniiibaKfxfAdxxffR,)22( 其中K为不依赖于)(xf的待定参数，),(ba。不难得到梯形公式的代数精度为 1，辛普森公式的代数精度为 3，所以梯形公式的误差余项可表示为：),(),(bafKfR ，其中：332233)(121)(6121)(2)(3121ababbaababK,所以梯形公式的余项为),(),(12)(3 bafabfR ，)32( 同理可得，辛普森公式的误差余项为 ),(,)2(180)4(4bafababfR，)42( 复合求积公式（梯形复合求积公式（梯形）（程序程序 22）由上述分析知当8n时 N-C 公式不具有稳定性，因此不能通
10、过提高求积公式的阶数来提高求积精度，要换一种思路来解决提高精度的问题。在前面的分析中，求积运算都是在一整个给定的区间上进行的， 既然能够在一个大区间上面进行求积计算， 那么把区间分成若干个小的区间时也能够在每一个区间上进行求积计算的。不妨将区间,ba分成n等分，若我们令分点为ihaxi，步长为nabh，ni,.,2 , 1 , 0，在每一个子区间) 1,.,1 , 0(,1nixxii上采用梯形公式，则得)()()(2)()(101101fRxfxfhdxxfdxxfInniiinixxbaii,记111 01 ( )() ( )2( )( ),22nnniii iihhTf xf xf af
11、 xf b )52( 称之为复合梯形公式复合梯形公式，其余项为： 1013 ),(),(12niiiiinnxxfhTIfR。由于,)(2baCxf，且)(max)(1)(min101010ininiiiniffnf ，所以),(ba使得 10)(1)(niifnf，于是复合梯形公式的余项为：)(122fhabfRn ，)62( 可以看出误差是2h阶，且当,)(2baCxf，有0lim fRn n，所以有 banndxxfT)(lim，因此复合梯形公式是收敛的。事实上，只要设,)(baCxf，则可得到收敛性，因为只要把nT改写为 )()(21110 niiniinxfnabxfnabT当n时，
12、 上式的右端括号内两个和式均收敛到积分dxxfba)(， 所以复合梯形公式收敛。此外，nT的求积系数为正，故复合梯形公式是稳定的。复合求积公式（辛普森复合求积公式（辛普森）（程序程序 23）将区间,ba分为n等分，在每个子区间,1iixx上采用辛普森公式，若记,22/1hxxii则得：)()(4)(6)()(1012/1101fRxfxfxfhdxxfdxxfInniiiinixxbaii,记,)()(2)(4)(6)()(4)(611102/11012/1niiniiniiiinbfxfxfafhxfxfxfhS)72( 称为复合辛普森求积公式复合辛普森求积公式，其余项为：),(, )()2
13、(180110)4(4 iiiniinnxxfhhSIfR，于是当,)(4baCxf时，与复合梯形公式相似有),(),()2(180)4(4bafhabSIfRnn，)82( 可以看出，误差阶为4h，收敛性是显然的，实际上，只要,)(baCxf，则可以得到收敛性，即.)(lim banndxxfS此外，由于nS中的求积系数均为正数，故知复合辛普森公式计算稳定。三、三、两种复化求积方法的步骤两种复化求积方法的步骤复化求积公式（梯形）的解题步骤复化求积公式（梯形）的解题步骤步骤一：步骤一：计算)(xf 的最大值，带入公式).(122fhabfRn ，由公式epsxfnabfhabfR bxa )(
14、max)1(12)()(12023 2确定区间等分次数n；步骤二步骤二：将区间,ba分成n等分，每个小区间分别为,1iixx，1,.,1 , 0ni，步长为nabh;步骤三：步骤三：计算每一个分点处的函数值)(ixf;步骤四：步骤四：将每一个)(ixf带入公式)52( 计算得到复合梯形公式的数值解；复合求积公式（辛普森）的解题步骤复合求积公式（辛普森）的解题步骤步骤一：步骤一：计算)()4(xf的最大值，带入公式).(2880)4(4fhabfRn，由公式epsxfnabfhabfR bxa )(max)1(2880)()(2880)4(45 )4(4确定区间等分次数；步骤二步骤二：将区间,ba分成n等分，每个小区间分别为,1iixx，1,.,1 , 0ni，步长为nabh记hxxii212/1，并求之；步骤三：步骤三：计算每一个分点处的函数值)(ixf和)(2/1ixf;步骤四：步骤四：将每一个)(ixf和)(2/1ixf带入公式)72( ，计算得到复合辛普森公式的数值解；四、两种方法比较四、两种方法比较由前面的理论推导不难发现，对于 N-C 公式进行复合后得到的计算结果精确度更高，而且对于同为 N-
