定积分的运.doc

实验四 定积分及应用 实验的目的 1、掌握利用Matlab进行积分运算； 2、掌握积分在计算面积、体积等问题中的应用； 3、掌握各种积分指令的区别与特点。 实验的基本理论与方法 1、定积分定义： 函数在区间上的定积分定义为：设函数在上有界，在区间上任取个分点：，把分成个小区间， 。这些分点构成对区间的一个分割，用T表示。小区间的长度为。记，称为分割T的模。在区间上取点，做函数值与小区间长度的乘积，并作和。当时，和总趋于确定的极限，这时这个极限为函数在区间上的定积分，记作。即 。 2、定积分的应用 ①计算平面图形的面积：由连续曲线，直线及x轴所围成的曲边梯形面积为； ②计算旋转体的体积：由连续曲线，直线及x轴所围成的曲边梯形绕x轴旋转一周所成立体的体积为； ③计算平面曲线的弧长：设曲线弧由直线坐标方程给出，其中在上具有一阶连续导数，则曲线弧长； 设曲线弧由参数方程给出，其中在上具有连续导数，则曲线弧长； 设曲线弧由极坐标方程给出，其中在上具有连续导数，则曲线弧长。 ④定积分的近似计算： 用分点将区间分成个长度相等的小区间，每个小区间长度为。 矩形法公式： 梯形法公式： 抛物线法公式： 实验使用的函数与命令 1、符号积分命令int int(fun)：求函数fun的不定积分； int(fun,var)：求函数fun关于变量var的不定积分； int(fun, var, a，b,)：求函数fun的在[a,b]间的定积分或广义积分； 2、数值计算定积分quad， quad1，trapz ①quad(fun,a,b,tol)：fun为被积的函数名，a, b微积分上下限，tol为精度，若缺省，其缺省值为1.0e-6； ②quad(fun,a,b,tol,trace)：参数fun,a,b,tol用法与上面相同，而输入第五个非零参数trace，是对积分过程通过被积函数上的图像进行跟踪。 对于quad使用自适应步长Simpson法，而quadl的调用格式与quad一致，但它使用Lobbato算法，其精度和速度要比quad的精度高。 ③Z = trapz(X,Y)：用于进行梯形积分，精度低，适用于数值函数和光滑性不好的函数。 其中X表示积分区间的离散化变量，Y与X同维的向量，表示被积函数，Z返回积分的近似值。 quad， quad1，trapz都不能用于广义积分，此外由于数值积分的特点，对一些假奇异积分也不能直接求解。 3、Inline 实验指导 例1 用符号积分命令int计算下列各题 1）不定积分与；2）定积分；3）广义积分 解：1）由于该题比较简单，我们直接在命令窗口输入Matlab代码 clear;syms x y z; int(sin(x*y+z),z) 得到的计算结果 ans = -cos(x*y+z) 再继续输入 int(x^2*sin(x)) 得到的计算结果 ans = -x^2*cos(x)+2*cos(x)+2*x*sin(x) 在以上的结果基础上，可以利用diff命令验证计算结果的正确性。 2）在命令窗口输入Matlab代码 syms x; int(1/(3+2*x+x^2),x,0,1) 得到定积分的值 ans = 1/2*atan(2^(1/2))*2^(1/2)-1/2*atan(1/2*2^(1/2))*2^(1/2) 3）在2）的基础上继续输入 int(1/(3+2*x+x^2),x,-inf,inf) 得到广义积分的值 ans = 1/2*pi*2^(1/2) 例2 计算数值积分。 解：我们先用trapz(X,Y)来求，在命令窗口输入Matlab代码 x=-1:0.1:2; y=x.^3;%取积分步长为0.1 trapz(x,y) 输出结果为 ans = 3.7575 我们现在把步长取到0.01 x=-1:0.01:2; y=x.^3; trapz(x,y) 输出结果为 ans = 3.7501 通过上面的运算结果可以看到，利用不同的步长进行计算，可以看到步长和精度之间的关系。如果取步长为0.001则输出结果为3.7500与步长取到0.01值接近。 可以利用符号积分命令计算 syms x; int(x^3,x,-1,2) ans = 15/4 下面我们利用quad，quad1来计算，首先要建立名为jifen.m的文件。 function y=jifen(x) y=x.^3; 然后进行在命令窗口进行调用并输出结果 quad(jifen,-1,2) ans = 3.7500： 利用quad1同样得到上面结果。 例3 计算数值积分。并用符号积分指令求解，观察输出结果。 解：利用trapz计算数值积分 x=0:0.1:1; y=sin(x.^2)./(1+x); tr