[理学]08实验八 多元函数积分学.pptVIP

图8.7 为了用线积分计算曲线的弧长, 必须把曲线用参数方程表示出来。因为空间曲线的投影曲线的方程为 , 它可以化成 , , 再代入锥面方程 , 得 ( )。因为空间曲线的弧长的计算公式是 , 因此输入: clear; syms t real x=cos(t)^2; y=cos(t)*sin(t); z=cos(t); qx=[x,y,z]; dx=diff(qx); f=sqrt(dot(dx,dx)); f=simplify(f); f=inline(f); quad(f,-pi/2,pi/2) 输出为: ans= 3.8202 8.2.4 计算曲面积分 【例10】计算曲面积分 , 其中 为锥面 被柱面 所截得的有限部分。 因为面积微元 , 根据曲面积分化作二重积分的计算公式, 并采用极坐标计算重积分。注意到投影曲线 的极坐标方程为 , 。输入: clear; syms x y z r t f=x*y+y*z+z*x; z=sqrt(x^2+y^2); f=subs(f,z,z); mj=sqrt(1+diff(z,x)^2+diff(z,y)^2); x=r*cos(t); y=r*sin(t); f=eval(f); mj=eval(mj); f1=f*mj*r; int(int(f1,r,0,2*cos(t)),t,-pi/2,pi/2) 输出为: ans= (64*2^(1/2))/15 【例11】计算曲面积分 , 其中 球面 的外侧。 可以利用两类曲面积分的关系, 化作对曲面面积的曲面积分 。这里 , 。因为球坐标的体积元素 , 注意到在球面 上 , 取 后得到面积元素的表示式: ( )。把对面积的曲面积分直接化作对 , 的二重积分。 输入: clear; syms x y z r s t a real A=[x^3,y^3,z^3]; n=[x,y,z]./a; x=a*sin(s)*cos(t); y=a*sin(s)*sin(t); z=a*cos(s); A=eval(A); n=eval(n); ds=a^2*sin(s); f1=dot(A,n); int(int(f1*ds,t,0,2*pi),s,0,pi) 输出为: ans= (12*pi*a^5)/5 如果用高斯公式计算, 则化为三重积分 , 其中 为 。采用球坐标计算, 输入: clear; syms x y z r s t a f=3*(x^2+y^2+z^2); x=r*sin(s)*cos(t); y=r*sin(s)*sin(t); z=r*cos(s); f=eval(f); int(int(int(f*r^2*sin(s),t,0,2*pi),s,0,pi),r,0,a) 输出结果相同。 MATLAB 高等数学实验 实验八 多元函数积分学 实验目的 掌握用MATLAB计算二重积分与三重积分的方法。深入理解曲线积分、曲面积分的概念和计算方法。提高应用重积分和曲线、曲面积分解决各种应用问题的能力。 8.1 学习MATLAB命令 8.1.1 重积分命令 命令int也可用于计算重积分, 例如要计算 输入: syms x y; int(int(x*y^2,y,0,x),x,0,1) 得到: ans= 1/15 用于计算三重积分时, 命令int的使用格式与此类似。 由此可见, 用MATLAB计算重积分, 关键是确定各个积分限。 8.1.2 二元函数的数值积分 函数dblquad的功能是求矩形区域上二元函数的数值积分。其格式为: q=dblquad(fun,xlower,xupper,ymin,y