定积分的近似计算.pptVIP

定积分的近似计算一、问题的背景和目的二、问题分析三、例题

一、问题的背景和目的定积分计算的基本公式是牛顿-莱布尼兹公式，但当被积函数的原函数不知道时，如何计算？这时就需要利用近似计算。特别是在许多实际应用中，被积函数甚至没有解析表达式，而是一条实验记录曲线，或一组离散的采样值，此时只能用近似方法计算定积分。本讲主要介绍定积分的三种近似计算算法：矩形法、复化梯形法和辛普森公式，及其误差分析。

二、问题的分析我们知道定积分，不论在实际问题中的意义是什么，在数值上都等于(设f(x)0)，直线与x轴所围成的曲边梯形的面积。因此，只要近似地算出相应的曲边梯形的面积，就得到了所给定积分的近似值。近似计算方法的基本思想还是分割、取点、求和这三个步骤。

定积分的定义定积分的近似基本思路

分割：为计算方便，一般采取等分法。把区间[a?b]n等分，即用分点a?x0，x1，x2，???，xn?1，xn?b?把区间[a?b]分成n个长度相等的小区间，每个小区间的长为取点：点可以任意选取，通常的取法有：左端点、右端点，就是左矩形法和右矩形法。

右矩形法

右矩形法如果取xi=xi，则可得近似公式这种方法称为右矩形法，上述公式称为右矩形公式。

梯形法

梯形法这种求定积分近似值的方法为复化梯形法，此公式称为复化梯形公式

辛普森公式

辛普森公式把区间[a?b]2n等分，每个小区间的长为过三点可以确定一条抛物线先计算[-h,h]上，以过三点的抛物线为曲线的曲边梯形的面积S：

辛普森公式同理可以得到区间上以抛物线为曲边的曲边梯形的面积：

辛普森公式特别的，当n=1时上述公式称为辛普森公式或抛物线公式。

误差估计当n=1时，左矩形法和右矩形法的余项的绝对值为两个矩形公式均具有一次代数精度。

误差估计复化梯形法的误差估计为

误差估计辛普森公式的误差估计为

三、例题试用复化梯形公式计算积分将区间[0，1]1000等分，并估计误差。

解：在Matlab中编写程序：functiont=ftrapz(a,b,n)h=(b-a)/n;t=h*(f(a)+f(b))/2;fori=1:1:(n-1)t=t+h*f(a+h*i);endI=t及functionf=f(x)ifx==0f=1;elsef=sin(x)/x;end应用复化梯形法求得.

要估计误差，要求而得复合梯化公式误差