实验名称数值积分1汇编.doc

探索实验7 数值积分法 实验目的 了解求积公式及代数精度概念，理解并掌握求定积分的求积公式的算法构造和计算，学习用计算机求定积分的一些科学计算方法和简单的编程技术和能用程序实现这些算法。 二、概念与结论 1. 求积公式： 计算定积分的如下形式的近似公式： 称为求积公式。 2．代数精度： 若求积公式 对一切不高于m次的 多项都准确成立，而对于m+1次多项式等号不成立，则称此求积公式 的代数精度为m。 代数精度越高，求积公式越好。 3.求积余项： 4.Newton-Cotes求积公式的代数精度 n点Newton-Cotes求积公式的代数精度至少可以达到n-1,且当n为奇数时，可以达到n。 5.Richardson外推定理： 设函数F1(h)逼近量F*的余项为： F*-F1(h)=a1h p1+a2h p2+····+a k p k+··· 式中p kp k-1···p2p10, F*和a i (i=1,2, ···)都是与h无关的常数,且k(1时，a k (0，则由： 定义的函数F2(h)也逼近F*，且有 F*-F2(h)= b2h p2+····+b k p k+··· 6. 关于复合梯形公式的展开定理 设f(x)在[a,b]区间上无穷次可微，则有如下展开式： T(h)=I+a1h2+a2h4+a3h6+…+amh2m+… 式中T(h)是函数f(x)在[a,b]区间上的复化梯形值Tn, 三、程序中Mathematica语句解释： 1. 随机函数 Random[] 随机给出闭区间[0,1]内的一个实数 Random[Real, xmax] 随机给出闭区间[0,xmax]内的一个实数 Random[Real, {xmin, xmax}] 随机给出闭区间[xmin,xmax]内的一个实数 Random[Integer] 随机给出整数0或1 Random[Integer, {xmin, xmax}] 随机给出xmin到xmax之间的一个整数 Random[Complex] 随机给出单位正方形内的一个复数 2．{a1,a2,…,an} 表示由元素a1,a2,…,an组成的一个表，元素可以是任何内容。 如：{1，3，4，5}，{1,x,{2,3},x+y}，{{1，3}，{1，2，3}，{3，2，4}}等 3．list[[k]] 表listk个元素 4．list[[i,j]] 表list第i元第j元,此时list中的第i元 四、方法与程序 在实际问题中，往往会遇到被积函数f(x)的原函数无法用初等函数来表示，或函数只能用表格表示，或有的虽然能用初等函数表示，但过分复杂，所以这些情形都需要去建立定积分的近似计算公式来做积分计算。数值积分是进行定积分计算的一种方法，它可以解决不能用定积分基本公式计算的所有定积分问题。数值积分涉及很多计算公式，这里主要介绍Newton-Cotes求积公式、复合求积公式、Romberg求积方法和Monte-Carlo方法的构造过程和算法程序。 n点 Newton—Cotes求积公式 n点 Newton—Cotes求积公式又称为等距节点求积公式，它是利用被积函数f(x)在积分区间[a,b]的n个等分节点上的函数值构造的插值函数((x)代替f(x)做定积分计算所构造求积公式。这个求积公式是通常做定积分近似计算的梯形公式和抛物线公式的推广，主要在理论上用的多些。 n点 Newton—Cotes求积公式的构造过程： 将积分区间[a,b] 分为n-1等分，其中n个节点 xi=a+(i-1}h, i=1,2,…,n，h=(b(a)/(n-1)，然后用f(x)在这n个节点上建立插值于f(x)的n-1次代数多项式Pn-1(x)，引入变换 x=a+th, 0(t(n-1 则有 带入定积分，有： Ck(n)称为Cotes(柯特斯)系数, 则得到n点 Newton—Cotes求积公式: n点 Newton—Cotes求积公式的求积余项为 当n=2时，2点的 Newton—Cotes求积公式就是如下梯形公式： 梯形求积公式求积余项为 当n=3时，3点的Newton—Cotes求积公式就是如下抛物线（Simpson）公式： Simpson求积公式求积余项为 如果想得到其