资源描述:
《数学建模集训小题目1-19题源程序》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、1•求下列积分的数值解:(•+8dx1%.Vx2-3x4-2附录:主程序:%辛普森公式zl=quad('1./x.*(x.2-3.*x+2)."(1/3)',2,inf)f1二inline(T./x.*(x."2-3.*x+2).八(1/3)');%inlineobject2.己知/(Z,h)=et+hcos(/+/?)+(/+/i)2sin(r+/z),g(h)=,画ill展[-10,10]时,g(防的图形。主程序:symst,symsh;f=exp(t+h)*cos(t+h)+(t+h)A2*sin(t+h);int(f,t,0,10)ans=1/2*exp(10
2、+h)*cos(10+h)+l/2*exp(10+h)*sin(l0+h)-98*cos(10+h)-20*cos(10+h)*h-cos(l0+h)*hA2+20*sin(10+h)+2*sin(10+h)*h-l/2*exp(h)*cos(h)-l/2*exp(h)*sin(h)+cos(h)*hA2-2*cos(h)-2*sin(h)*hezplot('1/2*exp(10+h)*cos(10+h)+l/2*exp(10+h)*sin(l0+h)-98*cos(l0+h)-20*cos(l0+h)*h-cos(10+h)*hA2+20*sin(10+h)+2*s
3、in(10+h)*h-l/2*exp(h)*cos(h)-l/2*exp(h)*sin(h)+cos(h)*hA2-2*cos(h)-2*sin(h)*h[-10,10])3.画岀F+(y_5)2=16绕兀轴一周所围成的图形,并求所产生的旋转体的体积。主程序:alpha=[0:0.1:2*pi]';beat=0:0.1:2*pi;x=4*cos(alpha)*ones(size(beat));y=(5+4*sin(alpha))*cos(beat);z=(5+4*sin(alpha))*sin(beat);surf(xyz)图2求体积dv=JJJ16_(尹-5)2d
4、ydz1Wy2+z2W9先计算在第I卦限的体积计算fJ16_(y_5),dy1dz2V9-y2Jclearsymszmym=sqrt(16-(y-5)A2)*(sqrt(9-yA2)-sqrt(l-yA2));int(m,l,9)ans=-118/3*i-19*3A(l/2)*EllipticK(1/2*2A(l/2))-35/3*3A(1/2)*EllipticPi(1/3,1/2*2A(1/2))+50*3A(1/2)^EllipticE(1/2*2A(1/2))-58/3*i*3A(l/2)*EllipticK(1/2*2A(1/2))+35/6*i*3A(l/
5、2)*EllipticPi(3/4,l/2*2A(l/2))+50*i*3A(l/2)*EllipticE(l/2*2A(l/2))+75/2*log(5)-75/2*log(-3+4*i)V=8(118/3*i-l9*3A(l/2)*EllipticK(l/2*2A(l/2))-35/3*3A(l/2)*EllipticPi(l/3,l/2*2A(l⑵)+50*3A(1/2)*EllipticE(l/2*2A(l/2))-58/3*i*3A(1/2)*EllipticK(l/2*2A(1/2))+35/6*i*3A(1/2)*EllipticPi(3/4,l/2*2
6、A(l/2))+50*i*3A(l/2)*EllipticE(1/2*2A(1/2))+75/2*log(5)-75/2*log(-3+4*i))4••画出下列曲面的图形⑴旋转单叶双曲面守-寻“主程序:t=-pi/4:0.1:pi/4;r=0:0.1:2*pi;[r,t]二meshgrid(r,t);x=3*sec(t).*sin(r);y=3*sec(t).*cos(r);z=2*tan(t);surf(x,y,z)(2)马鞍面z=xy;x=-2:0.1:2;y二-2:0.1:2;[xx,yy]二meshgrid(x,y);zz=xx.*yy;surf(xx,yy,
7、zz)5.画出隐函数sin+cosy=1的图形。程序:ezplot(,sin(x)+cos(y)-1')图表如下:6.(1)主程序symsx;diff(log((x+2)/(l-x)),3)(2)程序:a=[0,0,5,2,4];da=diff(a)da=05・327.⑴程序:[x,y]=solve(,xA2+y-6,,,yA2+x-67x,y,)(2)[x,y]=solve(,exp(x+sin(y))+cos(y)-107tan(x)+10*log(y)-5')&主程序:x=sym(,x,);diff(xA3+6*xA2+8*x-1)ans=3*