实验4__积分计算[精选].pptVIP

实验4 积分计算 实验目的： 学习 Matlab 命令 计算不定积分 定积分的概念 不定积分和广义积分 二重积分的计算 曲线积分 例3 . （2）Simpson公式 ------f 在[a，b]区间上积分的Simpson公式 。 用通过三点(a, f (a)), ((a+b)/2, f((a+b)/2)), (b, f(b))的抛物线围成的曲边梯形的面积来代替由 f 围成的曲边梯形的面积，由此获得积分的近似值. ------将区间[a，b] 分为n等分，步长为h=(b?a)/n ，分点为xk=a + kh (k=0,...,n)。在每个小区间[xk , xk+1]上采用辛普森公式，得到小区间上 f 积分的近似值。对每个小区间上积分的近似值求和，得到f 在[a，b]区间积分的近似值，这就是复合辛普森公式： (3) adaptive Simpson公式 ------通常被积函数f 在整个区间[a，b] 上变化不是很均匀的，如在某点附近函数变化非常急剧，而在其余地方上的变化比较平缓。这种情况用等距剖分小区间的复合求积公式不很合适。为了使积分达到预定的精度又要节省工作量，则可以在函数变化急剧的部分增多节点，即子区间分得细，而在函数变化平缓的地方减少节点，即子区间分得大，这个方法称为自适应积分法(adaptive quadrature methods) ------对于自适应辛普森方法，采用逐次将区间二等分的方法。为写法统一，将区间[a，b]记为[a，a + h]，其中h= b?a为区间的长度，称为0级区间。在区间[a, a + h]上采用辛普森公式，把结果记作 将区间分成两个相等的子区间[a, a+h/2]和[a+h/2, a+h], 这两个子区间称为1级子区间，其长度为h/2. 在每个子区间上 采用辛普森公式计算积分，然后相加并令 再将1级子区间中的一个或所有的两个二等分，所得的子区间称为2级子区间，其长度为h/22, ......。如此继续下去，最后将区间 [a, a+h]分成n个子区间[ai, ai+1] (i=0,1,...,n ? 1)，这样有 子区间的长度一般是不同的，如果子区间 [ai, ai+1]是r级，则其长度为 实际上，区间的划分是根据被积函数的变化而定的，函数变化平缓的地方，子区间大，函数变化急剧的地方，子区间小。 设Sa, a+h表示积分I(f )的近似值，那么有 ------自适应辛普森方法的matlab的命令为 z=quad(f, a, b, tol, trace) f(x) 为被积函数 a 为积分下限 b 为积分上限 tol 为计算精度，缺省为0.001 trace非0时，以动态点图的形式实现积分的整个过程。 ------注意：调用quad函数时，先要建立一个描述被积函数f(x) 的函数文件或语句函数。 例12. 用自适应辛普森方法求积分 解：使用语句函数（内联函数） （高精度近似值为1 clear; g=inline(exp(-x.^2)); z=quad(g,-1,1) z=1例12. 用自适应辛普森方法求积分 （高精度近似值为1 clear; g=inline(exp(-x.^2)); z=quad(g,-1,1,1e-6) z=1例12. 用自适应辛普森方法求积分 解：建立函数文件 （高精度近似值为1 function g=li12(x) g=exp(-x.^2); z=quad(li12,-1,1) z=1练习：用自适应辛普森方法法计算积分 6 二重积分 解：输入指令： syms x y; int(int(x*y,y,2*x,x^2+1),x,0,1) ans = 1/12 syms x y; int(int(sin(pi*(x^2+y^2)),y,-sqrt(1-x^2),sqrt(1-x^2)),x,-1,1) 解：积分区域可用不等式表成： 二重积分可化为二次积分 * * 1. 通过实验加深理解积分理论中分割、近似、求和、取极限的思想方法； 2. 学习并掌握matlab求不定积分、定积分、二重积分、曲线积分的方法； 3. 学习matlab命令sum，symsum，int。 实验内容： 1. 学习Matlab命令 1). 求和命令 sum 调用格式： sum(x) 给出向量x的各元素的累加和。若x为矩阵，则是一个元素为每列列和的行向量。 例1 x=[1,2,3,4,5,6,7