第3.3讲 数值积分与MATLAB求解.ppt

将积分区间[a,b] 划分为n等分, 步长 求积节点为 为了计 算系数Ak, 由于 , 所以 作变量代换 当 时,有 ,于是可得 求积系数 ( k=0,1…,n ) 代入插值求积公式(4)有 称为牛顿-柯特斯求积公式,Ck称为柯特斯系数 引进记号 ( k=0,1…,n ) 则 求积系数 容易验证 ∵ ∴ 显然, Ck是不依赖于积分区间[a,b]以及被积函数f(x)的常数,只要给出n,就可以算出柯特斯系数。 当n=2时 当n=1时 几个低阶求积公式 在牛顿-柯特斯求积公式中n=1,2,4时，就分别 得到下面的梯形公式、辛卜生公式和柯特斯公式。(1)?? 梯形公式 当n=1时，牛顿-柯特斯公式就是梯形公式 定理2 （梯形公式的误差）设f(x)在[a,b]上具有连续的二阶导数，则梯形公式的误差（余项）为 （2） 辛卜生公式 当n=2时，牛顿-柯特斯公式就是辛卜生公式（或 称抛物线公式） 定理3（辛卜生公式的误差）设在[a,b]上具有连续的四阶导数，则辛卜生求积公式的误差为： 定理证明从略。 （3） 柯特斯公式 当n=4时，牛顿-柯特斯公式为 定理4（柯特斯公式的误差）设在[a,b]上具有连续的6阶导数，则柯特斯求积公式的误差为 定理的证明从略。 例11 分别用梯形公式、辛卜生公式和柯特斯 公式计算定积分 的近似值 (计算结果取5位有效数字) (1) 用梯形公式计算 (2) 用辛卜生公式 (3) 用柯特斯公式计算，系数为 积分的准确值为 可见，三个求积公式的精度逐渐提高。 例12 用辛卜生公式和柯特斯公式计算定积分 的近似值,并估计其误差(计算结果取5位小数) 解: 辛卜生公式 由于 由辛卜生公式余项 知其误差为 例12 用辛卜生公式和柯特斯公式计算定积分 的近似值,并估计其误差(计算结果取5位小数) 解:柯特斯公式 知其误差为 例12 用辛卜生公式和柯特斯公式计算定积分 的近似值,并估计其误差(计算结果取5位小数) 该定积分的准确值 ,这个例子告诉我们，对于同一个积分，当n≥2时，公式却是精确的，这是由于辛卜生公式具有三次代数精度，柯特斯公式具有五次代数精度，它们对被积函数为三次多项式当然是精确成立的。 牛顿—柯特斯MATLAB应用 牛顿－柯特斯法 基于牛顿－柯特斯法，MATLAB给出了quad8函数来求定积分。 该函数的调用格式为： [I,n]=quad8(fname,a,b,tol,trace) 其中参数的含义和quad函数相似，只是tol的缺省值取10-6。 该函数可以更精确地求出定积分的值，且一般情况下函数调用的步数明显小于quad函数，从而保证能以更高的效率求出所需的定积分值。 例2 求定积分。 (1) 被积函数文件fx.m。 function f=fx(x) f=x.*sin(x)./(1+cos(x).*cos(x)); (2) 调用函数quad8求定积分。 I=quad8(fx,0,pi) I = 2.4674 牛顿－柯特斯法 例3 分别用quad函数和quad8函数求定积分的近似值，并在相同的积分精度下，比较函数的调用次数。 调用函数quad求定积分： format long; fx=inline(exp(-x)); [I,n]=quad(fx,1,2.5,1e-10) I = 0.28579444254766 n = 65 牛顿－柯特斯法 调用函数quad8求定积分： format long; fx=inline(exp(-x)); [I,n]=quad8(fx,1,2.5,1e-10) I = 0.28579444254754 n = 33 牛顿－柯特斯法 离散函数求积分 被积函数由一个表格定义 在MATLAB中，对由表格形式定义的函数关系的求定积分问题用trapz(X,Y)函数。 其中向量X,Y定义函数关系Y=f(X)。 例4 用trapz函数计算定积分。 命令如下： X=1:0.01:2.5; Y=exp(-X); %生成函数关系数据向量 trapz(X,Y) ans = 0.28579682416393 二重、三重定积分数值求解 二重