[信息与通信]华南师范大学2012年数学建模集训小题目.doc

2011年数学建模集训小题目 1.求下列积分的数值解 解：syms x f=1/(x*(x*x-3*x+2)^(1/3)),F=int(f,x,2,inf),y=double(F) f = 1/x/(x^2-3*x+2)^(1/3) F = int(1/x/(x^2-3*x+2)^(1/3),x = 2 .. Inf) y = 1.4396 2.已知，，画出时，的图形 解：syms t, h; f=exp(t+h)*cos(t+h)+(t+h)^2*sin(t+h); int(f,t,0,10) ans = -1/2*exp(h)*cos(h)-1/2*exp(h)*sin(h)+cos(h)*h^2-2*cos(h)-2*sin(h)*h+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 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]=meshgrid(-9:0.01:9); x=sqrt(16-(sqrt(y.^2+z.^2)-5).^2); x(16-(sqrt(y.^2+z.^2)-5).^20)=NaN; mesh(x,y,z);hold on;mesh(-x,y,z); 4.画出下列曲面的图形 （1）旋转单叶双曲面； 解：[x,y]=meshgrid(-10:0.1:10); z=2*sqrt((x.^2+y.^2)/9-1); z((x.^2+y.^2)/9-10)=NaN; mesh(x,y,z) （2）马鞍面； x=-10:0.1:10; y=-10:0.1:10; [X,Y]=meshgrid(x,y); Z=X.*Y; mesh(Z) 5.画出隐函数的图形。 解： 6.（1）求函数的三阶导数； 解：运用MATLAB syms x dx y=log((x+2)/(1-x)) y = log((x+2)/(1-x)) yx3=diff(y,x,3); y3=simple(yx3) y3 = -18*(x^2+x+1)/(x+2)^3/(-1+x)^3 （2）求向量的一阶向前差分。 a=[0 0.5 2 4] a = 0 0.5000 2.0000 4.0000 diff(a) ans = 0.5000 1.5000 2.0000 7.求解非线性方程组 （1） 解：syms x y;eq1=x*x+y-6; eq2=y*y+x-6; [x,y]=solve(eq1,eq2) 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） 解：syms x y; eq1=exp(x+sin(x))+cos(y)-10; eq2=tan(x)+10*log(y)-5; [x,y]=solve(eq1,eq2) x = 1.2981159100626495218126621342789 y = 1.1530435960638058760233858415577 8.求函数的极值点，并画出函数图形. syms x diff(x^3+6*x^2+8*x-1) ans =3*x^2+12*x+8 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,[-100,100]) 函数的极值点，x1=-2+2/3*3^(1/2), x2=-2-2/3*3^(1/2) 9.某单位需要加工制作100套钢架，每套用长为2.9m，2.1m和1m的圆钢各一根。已知原料长6.9m，问应如何下料，使用的原材料最省。 解：目标函数：minz=0*x1+0.1*x2+0.6