数学建模培训 课后十九小题详细详解.docVIP

1.M文件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 求出Q = 1.4389 Q = 1.4396 Q = 1.4393 因此，可以将1.4396作为此定积分的值。syms t h f=exp(t+h)*cos(t+h)+(t+h)^2*sin(t+h); g=int(f,t,1,10); ezplot(g) xlim([-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); 4.（1）程序： x=@(h,t)3*sec(h)*cos(t); y=@(h,t)3*sec(h)*sin(t); z=@(h,t)2*tan(h); ezmesh(x,y,z) （2）ezsurf(@(x,y)x*y) 5. ezplot(sin(x)+cos(y)=1) 6．（1）程序： syms x; diff(log((x+2)/(1-x)),3) 结果： ans = (6/(1-x)^3+6*(x+2)/(1-x)^4)/(x+2)*(1-x)-2*(2/(1-x)^2+2*(x+2)/(1-x)^3)/(x+2)^2*(1-x)-2*(2/(1-x)^2+2*(x+2)/(1-x)^3)/(x+2)+2*(1/(1-x)+(x+2)/(1-x)^2)/(x+2)^3*(1-x)+2*(1/(1-x)+(x+2)/(1-x)^2)/(x+2)^2 （2）程序： a=[0,0.5,2,4]; cha=diff(a) 结果： cha = 0.5000 1.5000 2.0000 7.（1） 程序：[x,y]=solve(x^2+y-6,y^2+x-6) 结果：x = 2 -3 1/2+1/2*21^(1/2) 1/2-1/2*21^(1/2) y = 2 -3 1/2-1/2*21^(1/2) 1/2+1/2*21^(1/2) （2） 程序： [x,y]=solve(exp(x+sin(y))+cos(y)-10,tan(x)+10*log(y)-5) 结果： x = 1.8015294326165098294767777230651 y = 2.5235903731067595761070805492882 8x=sym(x); diff(x^3+6*x^2+8*x-1) ans = -26 -43 -8 11 -12 78 -26 -44 -7 13 -14 78 -75 4 solve(3*x^2+12*x+8=0) ans = -2+2/3*3^(1/2) -2-2/3*3^(1/2) ezplot(x^3+6*x^2+8*x-1,[-200,200]) 9. 目标函数：minz=0*x1+0.1*x2+0.6*x3+0.7*x4+0.8*x5; 目标方程：x1+2*x2=100; 3*x3+2*x4+x5=100; 结果： Feasible solution found. Total solver iterations: 0 Variable Value MINZ 80.00000 X1 100.0000 X2 0.000000