6实验八主要命令.docVIP

PAGE 1 PAGE 80 实验八 多元函数积分学 一 实验目的 掌握用MATLAB计算二重积分与三重积分的方法；深入理解曲线积分、曲面积分的概念和计算方法. 提高应用重积分和曲线、曲面积分解决各种应用问题能力. 二 学习MATLAB命令 重积分命令int 命令int也可用于计算重积分，例如要计算, 输入 syms x y; int(int(x*y^2, y, 0,x), x,0,1) 命令int用于计算三重积分时的使用格式类似. 由此可见，用MATLAB计算重积分，关键是确定各个积分限. 2. 二元函数的数值积分 函数dblquad的功能是求矩形区域上二元函数的数值积分.其格式为 q = dblquad(fun,xlower,xupper,ymin,ymax,tol) % 用指定的精度tol代替缺省精度10-6，再进行计算. dblquad(inline(sin(exp(x.*y))),0,1,0,1) %可以把inline函数去掉也可以 dblquad(inline(sqrt(max(1-(x.^2+y.^2),0))),-1,1,-1,1) dot（vecf,dr） % 计算向量vecf与dr点积.例8会用到 qx=[x,y,z]; dx=diff(qx) % 向量求导.例9会用到 f=subs(f,’z’,z) % 把f中的符号z用z的新表达式代换.例10会用到 mj=eval(mj) % 计算mj的值. 三 实验内容 1. 计算重积分 例1 计算, 其中D为由,,所围成的有界区域. 先做出区域D的图，确定积分限.先对x后对y积分，输入 x=0:0.01:2; y1=2-x; y2=x.^2; y3=2; plot(x,y1,r-,x,y2,b-,x,y3,k-) syms x y; int(int(x*y^2, x, 2-y, sqrt(y)), y, 1,2) 输出为193/120. 预备知识：直角坐标重积分转换为极坐标的几种形式 一、极坐标—二重积分 选取极点?为直角坐标系的原点、极轴为轴，由直角坐标与极坐标的关系 有 二重积分的极坐标形式为 公式表明，要把一个二重积分变为极坐标系下的二重积分，只要把被积函数中的、分别换成、，面积元素换成即可。 计算, 其中D为. ， theta=0:0.1:2*pi; x=cos(theta); y=sin(theta); plot(x,y,r-) clear; syms r theta; f=exp(-(r^2))*r ; int(int(f,r,0,1),theta,0,2*pi) 输出为 -exp(-1)*pi+pi 请在直角坐标系下做做上面这个题目？？ 计算三重积分, 其中由曲面与 围成. 先作出区域的图形. 输入 [x,y]=meshgrid(-1:0.05:1); z= sqrt(2-x.^2-y.^2); surf(x,y,z) hold on z=sqrt(x.^2+y.^2); surf(x,y,z) 输出了区域的图形(图8.1) 图8.1 参照图形, 可以用手工确定积分限. 如果用直角坐标, 则输入 clear; syms x y z f=x^2+y^2+z; int(int(int(f, z, sqrt(x^2+y^2),sqrt(2-x^2-y^2)), y,-sqrt(1-x^2), sqrt(1-x^2)), x,-1,1) %y在交出来的那个圆上，那个圆的方程可以这样得到：由题目给的这两个方程 %和得，所以, %x只能从-1到1，而y的范围就是-sqrt(1-x^2)≤y≤sqrt(1-x^2). %执行上面语句后得不到明确结果，所以需要用下面的极坐标形式来做 执行后未得到明确结果. 改用柱坐标和球坐标计算. 预备知识： 二、柱面坐标—三重积分 令，于是 得到 三、球坐标—三重积分 令， 于是dv ==r 2 sinj drdqdj 得到 用柱坐标计算时输入 clear; syms r theta z f=(r^2+z)*r; %因为 int(int(int(f,z,r,sqrt(2-r^2)),r,0,1),theta ,0,2*pi) 输出为 -5/6*pi+16/15*2^(1/2)*pi 用球坐标计算时输入 clear; syms r phi theta f=(r^2*sin(phi)^2+r*cos(phi))*r^2*sin(phi); int(int(int(f, r, 0,sqrt(2)), phi,0,pi/4), theta,0,2*pi) simple(int(in