数值积分与微分.pptVIP

1、牛顿-莱布尼兹公式 另外若给出的函数f(x)是数据表，也不好求函数的积分。 计算定积分的方法： 但是求函数f(x)的原函数F(x)不一定比计算积分容易， 例如函数 找不到用初等函数表示的原函数。 一、数值求积的基本思想 实验4 数值积分与微分 主讲人：魏志强 数学实验课件 2、积分中值定理 但是点 的具体位置一般不知道，故难以准确算出 的值。 有两种近似方法： 得到 梯形公式 另一种是矩形法， 左矩形公式 一种是梯形法，用 代替 中矩形公式 右矩形公式 简称矩形公式 二、 数值积分基本概念 定义1：公式 叫做数 值求积公式。 其中xk称为求积节点， Ak称为求积系数,亦称伴随节点xk的权。 定义2：代数精度 若求积公式（1）对于次数不超过m的多项式均能准确成立，但对于m+1次的多项式就不能准确成立，则称此求积公式的代数精度为 m 。 如何判断代数精度？ 三、插值型的求积公式 设给定一组节点: 且已知函数在这些节点处的函数值f(xi)(i=0,1,…n),由第二章知可以作插值函数Ln(x),由于Ln(x)为多项式，所以其积分可以很容易求得： 记 则上式= 定义：称公式 其中 为求积系数。 为插值型的求积公式。 其余项为： 对等分的情况 ( k=0,1…,n ) 代入插值求积有 称为牛顿-柯特斯求积公式,Ck（n)称为柯特斯系数 引进记号 ( k=0,1…,n ) 则 n = 1: 梯形公式 代数精度 = 1 n = 2: 辛普森公式 代数精度 = 3 余项 余项 n = 4: 柯特斯公式 代数精度 = 5 余项 变步长求积方法 复化辛普森公式及其误差 将积分区间[a,b]划分为n等分,记子区间 的中点为 在每个小区间上应用辛卜生 公式，则有 记 称为复化辛普森公式 ，余项 龙贝格求积公式 梯形法的递推化 为了提高求积精度，可在复化求积的基础上将积分小区间[xk , xk+1]二分一次，增加了一个分点x k+1/2= (x k+x k+1) /2,记h=(b-a)/n，则： 称之为梯形法的递推化 龙贝格算法 由复化梯形公式的余项知： 有 即： 所以说明复化梯形公式二分前后两次计算值的线性组合就为复化辛普森求积公式。 同理对辛普生公式进行二分处理，前后两次计算值的误差进行比较， 复化柯特斯公式 有 即： 重复同样的操作，可进一步得到龙贝格（Romberg)公式： 高斯求积公式？ 四、 数值积分的实现方法 1．变步长辛普森法 基于变步长辛普森法，MATLAB给出了quad函数来求定积分。该函数的调用格式为： [I,n]=quad(fname,a,b,tol,trace) 其中fname是被积函数名。 a和b分别是定积分的下限和上限。tol用来控制积分精度，缺省时取tol=0.001。trace控制是否展现积分过程，若取非0则展现积分过程，取0则不展现，缺省时取trace=0。返回参数I即定积分值，n为被积函数的调用次数。 例1 求定积分。 (1) 建立被积函数文件fesin.m。 function f=fesin(x) f=exp(-0.5*x).*sin(x+pi/6); (2) 调用数值积分函数quad求定积分。 [S,n]=quad(fesin,0,3*pi) S = 0.9008 n = 77 To Matlab 2．牛顿－柯特斯法 MATLAB给出quadl函数来求定积分。 该函数的调用格式为： [I,n]=quadl(fname,a,b,tol,trace) 其中参数的含义和quad函数相似，只是tol的缺省值取10^(-6)。该函数可以更精确地求出定积分的值，且一般情况下函数调用的步数明显小于quad函数，从而保证能以更高的效率求出所需的定积分值。 如：[S,n]=quadl(fesin,0,3*pi) 例2 求定积分。 (1) 被积函数文件fx1.m。 function f=fx1(x) f=x.*sin(x)./(1+cos(x).*cos(x)); (2) 调用函数quadl求定积分。 I=quadl(fx1,0,pi) I = 2.4674 调用函数quad求定积分(可以文件存储）： format long; fx=inline(exp(-x)); [I,n]=quad(fx,1,2.5,1e-1