Matlab与数学建模综合演习2.doc

求的所有根。（先画图后求解）syms x; f=exp(x)-3*x^2; ezplot(f) solve(exp(x)-3*x^2) ans = -2*lambertw(-1/6*3^(1/2)) -2*lambertw(-1,-1/6*3^(1/2)) -2*lambertw(1/6*3^(1/2)) 求下列方程的根。 1) solve(x^5+5*x+1) ans = 1.1044655068824455162575638841973+1.0598296691525201166749456468980*i -1.0044974557968355184823910746206+1.0609465064060406435760940804509*i -.19993610217121999555034561915339 -1.0044974557968355184823910746206-1.0609465064060406435760940804509*i 1.1044655068824455162575638841973-1.0598296691525201166749456468980*i 2） solve(x*sin(x)-1/2) ans = -.74084095509549062101093540994313 f=x*sin(x)-1/2; ezplot(f) 3） f=sin(x)*cos(x)-x^2; solve(sin(x)*cos(x)-x^2) ans = 0. ezplot(f) 求解下列各题： 1） syms x; limit((x-sin(x))/x^3,x,0) ans = 1/6 2） syms x; diff(exp(x)*cos(x),x,10) ans = -32*exp(x)*sin(x) 3） syms x; f=int(f,x,0,1/2); vpa(f,17) ans = .54498710418362220 4）syms x; int(x^4/(25+4*x^2),x) ans = 1/12*x^3-25/16*x+125/32*atan(2/5*x) 5） syms x; f=sqrt(x+1); taylor(f,9) ans = 1+1/2*x-1/8*x^2+1/16*x^3-5/128*x^4+7/256*x^5-21/1024*x^6+33/2048*x^7-429/32768*x^8 6） syms x; f=diff(exp(sin(1/x)),x,3) f = cos(1/x)/x^6*exp(sin(1/x))+6*sin(1/x)/x^5*exp(sin(1/x))+3*sin(1/x)/x^6*cos(1/x)*exp(sin(1/x))-6*cos(1/x)/x^4*exp(sin(1/x))-6*cos(1/x)^2/x^5*exp(sin(1/x))-cos(1/x)^3/x^6*exp(sin(1/x)) limit(f,x,2) ans = -23/64*cos(1/2)*exp(sin(1/2))+3/64*exp(sin(1/2))*sin(1/2)*cos(1/2)-3/16*exp(sin(1/2))*cos(1/2)^2+3/16*exp(sin(1/2))*sin(1/2)-1/64*exp(sin(1/2))*cos(1/2)^3 1)求矩阵 的逆矩阵 及特征值和特征向量。 a=[-2 1 1;0 2 0;-4 1 3]; inv(a) ans = -1.5000 0.5000 0.5000 0 0.5000 0 -2.0000 0.5000 1.0000 2)求点(1,1,4)到直线l: (x-3)/-1 =y/0=(z+1)/2的距离。 已知分别在下列条件下画出的图形： 、 画 下列函数的图形：（1） t=[0:0.01:20]; u=[0:0.01:2]; [T,U]=meshgrid(t,u); x=U.*sin(T); y=U.*cos(T); z=T/4; mesh(x,y,z) (2) （2）x=[0:0.01:3]; y=[0:0.01:3]; [X,Y]=meshgrid(x,y); z=sin(X.*Y); mesh (x,y,z) （3）.t=[0:0.1:2*pi]; u=[0:0.1:2*pi]; [T,U]=meshgrid(t,u); x=sin(T*(3+cos(U))); y=cos(T*(3+cos(U))); z=sin