数学建模小题目及答案.docVIP

求下列积分的数值解： function y = myfun(x) y = 1./(x.*(x.^2 - 3*x + 2 ).^(1/3)); warning off all Q = quad(@myfun,2,100000) Q = quad(@myfun,2 Q = quad(@myfun,2,1000000000000000) warning on 当上限为1000001000000000时， 定积分的值为x=1.4389,1.4396,1.4396。 因此，可以将1.4396作为此定积分的值。 已知，，画出时，的图形。 syms t,syms h; f=exp(t+h)*cos(t+h)+(t+h)^2*sin(t+h); int(f,t,0,10) ans = 1/2*exp(10+h)*cos(10+h)+1/2*exp(10+h)*sin(10+h)-98*cos(10+h)-20*cos(10+h)*h-cos(10+h)*h^2+20*sin(10+h)+2*sin(10+h)*h-1/2*exp(h)*cos(h)-1/2*exp(h)*sin(h)+cos(h)*h^2-2*cos(h)-2*sin(h)*h ezplot(1/2*exp(10+h)*cos(10+h)+1/2*exp(10+h)*sin(10+h)-98*cos(10+h)-20*cos(10+h)*h-cos(10+h)*h^2+20*sin(10+h)+2*sin(10+h)*h-1/2*exp(h)*cos(h)-1/2*exp(h)*sin(h)+cos(h)*h^2-2*cos(h)-2*sin(h)*h,[-10,10]) 3.画出绕轴一周所围成的图形，并求所产生的旋转体的体积。 主程序： [y,z]=cylinder(1:0.2:9,100); mesh(sqrt(16-(sqrt(y.^2+z.^2)-5).^2),y,z); hold on; mesh(-sqrt(16-(sqrt(y.^2+z.^2)-5).^2),y,z); 求体积dv=dydz 先计算在第Ⅰ卦限的体积 1≤y+z≤9 计算dydz =- clear syms z m y m=sqrt(16-(y-5)^2)*(sqrt(9-y^2)-sqrt(1-y^2)); int(m,1,9) ans = -118/3*i-19*3^(1/2)*EllipticK(1/2*2^(1/2))-35/3*3^(1/2)*EllipticPi(1/3,1/2*2^(1/2))+50*3^(1/2)*EllipticE(1/2*2^(1/2))-58/3*i*3^(1/2)*EllipticK(1/2*2^(1/2))+35/6*i*3^(1/2)*EllipticPi(3/4,1/2*2^(1/2))+50*i*3^(1/2)*EllipticE(1/2*2^(1/2))+75/2*log(5)-75/2*log(-3+4*i) V=8(118/3*i-19*3^(1/2)*EllipticK(1/2*2^(1/2))-35/3*3^(1/2)*EllipticPi(1/3,1/2*2^(1/2))+50*3^(1/2)*EllipticE(1/2*2^(1/2))-58/3*i*3^(1/2)*EllipticK(1/2*2^(1/2))+35/6*i*3^(1/2)*EllipticPi(3/4,1/2*2^(1/2))+50*i*3^(1/2)*EllipticE(1/2*2^(1/2))+75/2*log(5)-75/2*log(-3+4*i)) 4.画出下列曲面的图形 （1）旋转单叶双曲面； x=@(s,t)3.*sec(s).*cos(t); y=@(s,t)3.*sec(s).*sin(t); z=@(s,t)2.*tan(s); ezmesh(x,y,z) 或者 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）马鞍面； x=-2*pi:0.2:2*pi; [x,y]=meshgrid(x); z=x.*y; surf(x,y,z); xlabel(x-axis),ylabel(y-axis),zlabel(z-axis); title(surf)； 或者 x=-2:0.1:2; y=-2:0.1:2; [xx,yy]=meshgrid(x,y); zz=xx.*yy;