§6 Mathematica求定积分以及相关应用问题.doc

6 Mathematica求定积分以及相关应用问题6.1 用Mathematica求定积分1 定积分的运算在不定积分中加入积分的上下限便成为定积分(definite integral)。Mathematica的定积分命令和不定积分的命令相同，但必须指定积分变量的上下限。(1) Integratef,x,下限,上限(2) 例6.1 计算定积分。解 Out1=4-2ArcTan2和不定积分一样,除了我们指定的积分变量之外,其它所有符号都被作常数处理.例6.2 计算定积分。解 2 数值积分如果Mathematica无法解出积分的符号表达式或者定积分的结果过于冗长而失去意义时，我们就可以用数值积分求解。数值积分只能进行定积分的运算，即必须指定上、下限。用Mathematica求解数值积分有两种形式：(1) NIntegratef,x,a,b 从到 ，做的数值积分。(2) N 求定积分表达式的数值例6.3 求定积分。解 用Integrate命令无法求的定积分，用NIntegrate命令即可求得其数值积分。In1:=NIntegrateSinSinx,x,0,Pi/3Out1=0.466185求定积分表达式的数值，也能得到与上式相同的结果。Out2=0.466185例6.4 求定积分的近似值。解 被积函数的原函数不能被等函数表示，我们可以计算它的数值积分。In3:=NIntegrateExp-x2,x,0,1Out3=0.7468243 近似值积分 用Mathematica计算定积分的近似值还有矩形法、梯形法和抛物线法用分点将区间分成个长度相等的小区间，每个小区间长度为 矩形法公式：梯形法公式： 抛物线法公式：例6.5 分别用矩形法、梯形法、抛物线法计算定积分。解 为了便于比较，首先计算积分的精确值：In1:=Clearx;yx_=x2;Integrateyx,x,0,1Out1=(1) 矩形法In2:=Cleary,x,s1,n,b,a;n=20;a=0;b=1;yx_:=x2;s1=(b-a)/n*Sumya+i(b-a)/n,i,0,n-1/N;s2=(b-a)/n*Sumya+i(b-a)/n,i,1,n/N;Print“s1=”,s1” s2=”,s2Out2=s1=0.30875 s2=0.35875 (2) 梯形法In3:=Cleary,x,a,b,ss3,s3;yx_:=x2;n=20;a=0;b=1;ss3=Sumya+i*(b-a)/n,i,1,n-1;s3=(ya/2+yb/2+ss3)*(b-a)/n /N;Print“s3=”,s3Out3=0.33375(3) 抛物线法In4:=Cleary,x,a,b,s3;yx_:=x2;n=20;a=0;b=1;m=10;ss1=Sum(1+(-1)i)*ya+i*(b-a)/n,i,1,n-1;(*ss1=2y2+2y4+2yn-2*)ss2=Sum(1-(-1)i)*ya+i*(b-a)/n,i,1,n-1;(*ss2=2y1+2y3+2yn-1*)s4=N(ya+yb+ss1+2ss2)*(b-a)/3/n,20;Prints4=”,s4Out4=0.33333333333333333333由上述结果可知：抛物线法近似程度最好，矩形法近似程度最差。6.2 用Mathematica计算相关定积分应用问题在解有关定积分应用问题时会用到的Mathematica函数有以下几种：1、 Solve方程1，方程2，变量1，变量2：求解二元方程组。2、 Plotfx,x,a,b：画一元函数图形。3、 ParametricPlotxt,yt,t,t1,t2：二维参数作图。4、 Integratefx,x,a,b：计算定积分。5、 Showf1,f2：将函数组合显示。1 利用定积分计算平面图形的面积有连续曲线 ，直线及轴所围成的曲边梯形的面积为 例6.6 求由抛物线和直线所围成图形的面积。解 首先画出函数图形，如图6-1所示In1:=PlotSqrt2x,-Sqrt2x,-x+4,x,0,9Out1=-Graphics- 图6-1然后求出两条曲线的交点：In2:=Solvey2-2=0,y+x-4=0,x,yOut2=x2,y2,x8,y-4再以y为积分变量求面积：In3:=s=Integrate-y+4-y2/2,y,-4,2Out3=18例6.7 求由圆所围图形的面积。解 首先求出两条曲线的交点：In4:=Solver-3Cost=0,r-1-Cost=0,r,tOut4=然后画出两曲线所围成的图形，如图6-2所示In5:=f1=ParametricPlot3CostCost,3CostSint,t,0,2Pi;f2= ParametricPlot(1+Cost)Cost,(1+Cost)Sint,t,0,2PiShowf1,f2,AspectRatioAutomaticOut5= -Graphics- 图 6-2再利用定积分计算面积 In6:=s1=Integrate1+Cost,t,0,Pi/3Out6=In7:=s2=Integrate3Cost,t,Pi/3,Pi/2Out7=In8:=s=2*(s1+s2)Out8=2 利用定积分计算平面曲线的弧长 设曲线弧由参数方程：,给出，其中在上具有连续导数，则曲线弧的弧长为 例6.8 求曲线 上相应于从到的一段弧长。解 首先画出曲线的图形，如图6-3所示In1:=ParametricPlotArcTant,(1/2)*Log1+t2,t,-2,2,AspectRatioAutomaticOut1= -Graphics- 图6-3再利用定积分计算曲线的弧长：In2:=dx=DArcTant,tIn3:=dy=D(1/2)Log1+t2,tIn4:=s=IntegrateSqrtdy2+dx2,t,0,1/NOut4=0.8813743 利用定积分计算旋转体的体积由连续曲线,直线及轴所围成的曲边梯形绕轴旋转一周所成立体的体积为例6.9 将星形线所围成的图形绕轴旋转一周，计算所得旋转体的体积。 解 星形线的参数方程为 取a=1，画出星形线的图形，如图6-4所示In1:=ParametricPlot(Cost3),(Sint3),t,0,2Pi,AspectRatioAutomaticOut1= -Graphics- 图6-4 利用Mathematica计算旋转体的体积：In2:=xt_:=a*Cost3;yt_:=a*Sint3;dx=Dxt,t;V=2*IntegratePi*(yt)2*dx,t,0,Pi/2练习5.61. 用Mathematica求解下列定积分：(1);