matlab实验三 定积分的近似计算.doc

５３实验三定积分的近似计算一、问题背景与实验目的利用牛顿莱布尼兹公式虽然可以精确地计算定积分的值，但它仅适用于被积函数的原函数能用初等函数表达出来的情形．如果这点办不到或者不容易办到，这就有必要考虑近似计算的方法．在定积分的很多应用问题中，被积函数甚至没有解析表达式，可能只是一条实验记录曲线，或者是一组离散的采样值，这时只能应用近似方法去计算相应的定积分．本实验将主要研究定积分的三种近似计算算法矩形法、梯形法、抛物线法．对于定积分的近似数值计算，MATLAB有专门函数可用．二、相关函数(命令)及简介1．SUMA求数组A的和．2．ATLONG长格式，即屏幕显示15位有效数字．(注由于本实验要比较近似解法和精确求解间的误差，需要更高的精度)．3．DOUBLE若输入的是字符则转化为相应的ASCII码；若输入的是整型数值则转化为相应的实型数值．4．QUAD抛物线法求数值积分．格式QUADFUN,A,B，注意此处的FUN是函数，并且为数值形式的，所以使用、/、等运算时要在其前加上小数点，即、/、等．例QQUAD 1/X32X5 ,0,25．TRAPZ梯形法求数值积分．格式TRAPZX,Y其中X为带有步长的积分区间；Y为数值形式的运算(相当于上面介绍的函数FUN)例计算0SINDXX0PI/100PIYSINXTRAPZX,Y6．DBLQUAD抛物线法求二重数值积分．格式DBLQUADFUN,XMIN,XMAX,YMIN,YMAX，FUN可以用INLINE定义，也可以通过某个函数文件的句柄传递．例1Q1DBLQUADINLINE YSINX ,PI,2PI,0,PI顺便计算下面的Q2，通过计算，比较Q1与Q2结果(或加上手工验算)，找出积分变量X、Y的上下限的函数代入方法．Q2DBLQUADINLINE YSINX ,0,PI,PI,2PI例2Q3DBLQUADINTEGRND,PI,2PI,0,PI这时必须存在一个函数文件INTEGRNDM５４FUNCTIONZINTEGRNDX,YZYSINX7．FPRINTF(文件地址，格式，写入的变量)把数据写入指定文件．例X011YXEXPXFID EXPTXT , W 打开文件FPRINTFFID, 62F128F\N ,Y写入FCLOSEFID关闭文件8．SYMS变量1变量2定义变量为符号．9．SYM 表达式 将表达式定义为符号．解释MATLAB中的符号运算事实上是借用了MAPLE的软件包，所以当在MATLAB中要对符号进行运算时，必须先把要用到的变量定义为符号．10．INTF,V,A,B求F关于V积分，积分区间由A到B．11．SUBSF， X ，A将A的值赋给符号表达式F中的X，并计算出值．若简单地使用SUBSF，则将F的所有符号变量用可能的数值代入，并计算出值．三、实验内容1．矩形法根据定积分的定义，每一个积分和都可以看作是定积分的一个近似值，即1DNBIAIFXFX在几何意义上，这是用一系列小矩形面积近似小曲边梯形的结果，所以把这个近似计算方法称为矩形法．不过，只有当积分区间被分割得很细时，矩形法才有一定的精确度．针对不同的取法，计算结果会有不同，我们以为例(取)，I120DX10N(1)左点法对等分区间，BXINABXAXNI10在区间上取左端点，即取，,1II1II078789399673078，1201DNIIXFX理论值，此时计算的相对误差1204X78936078403178５５(2)右点法同(1)中划分区间，在区间上取右端点，即取，,1IIXIIX078289399673078，1201DNIIXFX理论值，此时计算的相对误差1204X78293607840318(3)中点法同(1)中划分区间，在区间上取中点，即取1,IIX，2IIIX078540024673078，1201DNIIXFX理论值，此时计算的相对误差1204X67852463078425310如果在分割的每个小区间上采用一次或二次多项式来近似代替被积函数，那么可以期望得到比矩形法效果好得多的近似计算公式．下面介绍的梯形法和抛物线法就是这一指导思想的产物．2．梯形法等分区间，BXINABXAXNI10NA相应函数值为()．NY,10IFI,10,曲线上相应的点为XFY()NP,10NIYXII,,将曲线的每一段弧用过点，的弦(线性函数)来代替，这使得I1IIIP1每个上的曲边梯形成为真正的梯形，其面积为,1IIX５６，．XYII21NI,21于是各个小梯形面积之和就是曲边梯形面积的近似值，，11D2NNBIIIIAIYFXXY即，011BNAAFYN称此式为梯形公式．仍用的近似计算为例，取，120DX078539399673078，0112NYYBAN理论值，此时计算的相对误差1204X6078539607845301很显然，这个误差要比简单的矩形左点法和右点法的计算误差小得多．3．抛物线法由梯形法求近似值，当为凹曲线时，它就偏小；当为凸曲XFYXFY线时，它就偏大．若每段改用与它凸性相接近的抛物线来近似时，就可减少上述缺点，这就是抛物线法．将积分区间作等分，分点依次为,BAN2，，BXINABXXNI2102NAX2对应函数值为()，NY210,IFI,10,曲线上相应点为()．NP210,NIYXII2,,现把区间上的曲线段用通过三点，，,20XFY,0YXP,1YX的抛物线,2YP12XPXY５７来近似代替，然后求函数从到的定积分1XP02201DXP20DX3022030XXX462020222X由于，代入上式整理后得201X201DXP461220202XXX42102YX6210YNAB同样也有42DXP4322NX6212NNYYAB将这个积分相加即得原来所要计算的定积分的近似值N，222111DD4IINNBXIIIAIIFP即02132124246BNNNAFXYYYY这就是抛物线法公式，也称为辛卜生(SIMPSON)公式．仍用的近似计算为例，取，120DX002132124246NNNBAYYYY078539816339745，理论值，此时计算的相对误差120D4X160785391674528704直接应用MATLAB命令计算结果５８(1)数值计算120DX方法1INT 1/1X2 , X ,0,1符号求积分方法2QUAD 1/1X2 ,0,1(抛物线法求数值积分)方法3X000011Y1/1X2TRAPZX,Y(梯形法求数值积分)(2)数值计算2120DXY方法1INTINT XY2 , Y ,1,1, X ,0,2(符号求积分)方法2DBLQUADINLINE XY2 ,