实验四+不定积分、定积分及其应用.docVIP

实验四 不定积分、定积分及其应用 【实验类型】验证性 【实验学时】2学时 【实验目的】 1．掌握用MATLAB求函数方法； 2．； 【实验内容】【实验方法或步骤】 一、实验使用的MATLAB函数 的不定积分； 2．int( f(x), x , a , b)；求在上的定积分； 3．int( f(x) , x , -inf, inf )；计算广义积分； 4．solve(eqn1,eqn2,...,eqnN,var1,var2,...,varN)；求解n元方程组； 二、实验指导 。 syms x; int(exp(x)*cos(2*x),x) ans = 1/5*exp(x)*cos(2*x)+2/5*exp(x)*sin(2*x) 例2 计算不定积分。 syms x; int(1/(x^4*sqrt(1+x^2))) ans = -1/3/x^3*(1+x^2)^(1/2)+2/3/x*(1+x^2)^(1/2) 例3 以几何图形方式演示、理解定积分概念，并计算近似值。 先将区间任意分割成n份，为保证分割加细时，各小区间的长度趋于0，在取分点时，让相邻两分点的距离小于，分点取为（为随机数），在每一区间上任取一点（为随机数）作积分和进行计算，程序如下： function juxs(fname,a,b,n) % 定积分概念演示，随机分割、 随机取近似，并求近似值 xi(1)=a; xi(n+1)=b; for i=1:n-1 xi(i+1)=a+(i+rand(1))*(b-a)/n; end I=0; hold on; for i=1:length(xi)-1 sxi=xi(i)+rand(1)*(xi(i+1)-xi(i)); syi=feval(fname,sxi); I=I+syi*(xi(i+1)-xi(i)); xii=[xi(i) xi(i) xi(i+1) xi(i+1) xi(i)]; yii=[0 syi syi 0 0]; fill(xii,yii,c); end x=a:(b-a)/100:b; y=feval(fname,x); plot(x,y,r,markersize,20); hold off; fprintf(n=%6d, I=%12.5f\n,n,I); function y=fex(x) y=x.^2+1; 以积分为例，调用上述函数，观察如下： （1）几何上 我们知道，当时，定积分的值表示由所围成的曲边梯形的面积，从图形3.1上可以看出，用小矩形面积和逼近曲边梯形面积的过程。 值得注意的是，虽然每次运行后的图形可能有所差异（相同的参数下），但总的趋势是，分点个数越多，小矩形的面积之和越逼近曲边梯形的面积，即积分和越逼近积分值。 （2） 数值上 当对区间逐步进行细分时，反复调用上述程序，可得一系列积分近似值（运行结果可能有差异），可以看到，随着区间数的增大，近似值越来越接近精确值（精确值为28/3）。 n= 20, I= 9.12818 n= 40, I= 9.38262 n= 160, I= 9.34459 n= 320, I= 9.33352 n= 640, I= 9.33158 n= 1280, I= 9.33324 n= 2480, I= 9.33364 例4 计算定积分。 syms x; int(1/(5+3*sin(x)),x,0,2*pi) ans = 1/2*pi 例5 计算定积分。 syms x; int(x^2*sqrt(a^2-x^2),x,0,a) ans = 1/16*a^5*(1/a^2)^(1/2)*pi 例6 计算广义积分。 syms x; int(1/(x^2+2*x+2),x,-inf,inf) ans = 1/4*pi 例7 判别广义积分的敛散性。 syms x; int(1/sqrt(1-x^2),x,0,1) ans = 1/2*pi 所以原积分收敛。 例8 求由抛物线和直线所围图形的面积。 首先画出函数图形，如图3.2所示。 x=0:0.1:9; plot(x,-x+4,’b’,x,sqrt(2*x),r,x,-sqrt(2*x),r) 求解方程组，得到两曲线交点 [x,y]=solve(y^2-2*x=0,y+x-4=0)； x = [8 2] , y = [-4 2] 以 y为积分变量求面积 int(-y+4-y^2/2,y,-4,2) ans = 18 例9 求由圆和双纽线所围