|
二、MATLAB实例应用 例2、试求出曲面积分，其中是椭圆球面的的部分，取椭球面外侧为正侧.[2] copyright paper51.com [求解]：画出积分图用下面语句画出椭球 paper51.com >>[X,Y,Z]=ELLIPSOID(0,0,0,1,2,3,40); copyright paper51.com >> surf(X,Y,Z) paper51.com
http://www.paper51.com 图四 椭球 paper51.com 引入参数方程且，，我们就可以用下面语句求出所需的曲面积分. copyright paper51.com >>syms u v;syms a b c positive; http://www.paper51.com x=a*sin(u)*cos(v);y=b*sin(u)*sin(v);z=c*cos(u);A=diff(y,u)*diff(z,v)-diff(z,u)*diff(y,v);I =int(int(x^3*A,u,0,pi/2),v,0,2*pi) paper51.com I = 2/5*pi*a^3*c*b 内容来自www.paper51.com 例3、计算[3] 内容来自论文无忧网 www.paper51.com 要用传统的算法去求繁琐复杂，根本达不到所需的精确度.下面我们用MATLAB语言来求其精确值. paper51.com >> syms x y;s1=int(int(y*sin(x)+x*cos(y),x,pi,2*pi),y,0,pi) 内容来自www.paper51.com s1 =-pi^2 内容来自论文无忧网 www.paper51.com
>> s2=vpa(s1) 内容来自www.paper51.com s2 =-9.8696044010893586188344909998761 内容来自www.paper51.com
内容来自论文无忧网 www.paper51.com |