理工数学实验答案matlab.pdfVIP

第五章 5 ． syms x y u v dux duy dvx dvy f=x^2+y^2+u^2+v^2-R^2; g=x+y+u+v; fx=diff(f,x);fy=diff(f,y);fu=diff(f,u);fv=diff(f,v); gx=diff(g,x);gy=diff(g,y);gu=diff(g,u);gv=diff(g,v); ffx=fu*dux+fv*dvx+fx; ggx=gu*dux+gv*dvx+gx; [dux,dvx]=solve(ffx,ggx,dux,dvx) dux -(x-v)/(u-v) dvx (-u+x)/(u-v) ffy=fu*duy+fv*dvy+fy; ggy=gu*duy+gv*dvy+gy; [duy,dvy]=solve(ffy,ggy,duy,dvy) duy -(y-v)/(u-v) dvy (-u+y)/(u-v) 6 ． syms x y z f=x+y+z-exp(-(x+y+z)); fx=diff(f,x);fz=diff(f,z); dzx=-fx/fz; fxy=diff(dzx,y) fxy exp(-x-y-z)/(1+exp(-x-y-z))+(-1-exp(-x-y-z))/(1+exp(-x-y-z))^2*exp(-x-y-z) 7 ． syms x y z=x^2-x*y-2*y^2; zx=diff(z,x);zy=diff(z,y); fzx=inline(zx); fzy=inline(zy); a=pi/6;b=pi/3; f1=fzx(1,2)*cos(a)+fzy(1,2)*cos(b) f1 -4.5000 8 ． syms t x=sin(t) y=cos(t) z=t; dx=diff(x,t);dy=diff(y,t);dz=diff(z,t); x1=inline(dx); x2=inline(dy); x3=inline(dz); x10=x1(1) x10=x1(pi/4) x10 0.7071 x20=x2(pi/4) x20 -0.7071 x30=x3(pi/4) x30 1 9. syms x y f=x^2+y^2; fx=diff(f,x); fy=diff(f,y); x=1; y=2; fx0=subs(fx) fx0 2 fy0=subs(fy) fy0 4 10. syms x y f=x^3-y^3+3*x^2+3*y^2-9*x; fx=diff(f,x);fy=diff(f,y); [x0 y0]=solve(fx,fy) x0 1 -3 1 -3 y0 0 0 2 2 fxx=diff(diff(f,x),x) fxx 6*x+6 fxy=diff(diff(f,x),y); fyy=diff(diff(f,y),y); delta=inline(fxy^2-fxx*fyy); delta(x0,y0) ans -72 72 72 -72 x=1;y=0; fmin=subs(f) fmin -5 x=-3;y=2; fmax=subs(f) fmax 31 10． syms x y z lamda a L=x*y*z+lamda*(x^2/2+y^2/2+z^2-a*a); Lx=diff(L,x); Ly=diff(L,y); Lz=diff(L,z); Llamda=diff(L,lamda); [lamda x y z]=solve(Lx,Ly,Lz,Llamda) lamda 0 0 0 0 0 0 1/3*3^(1/2)*a -1/3*3^(1/2)*a 1/3*3^(1/2)*a -1/3*3^(1/2)*a -1/3*3^(1/2)*a 1/3*3^(1/2)*a -1/3*3^(1/2)*a 1/3*3^(1/2)*a x 0 0 2^(1/2)*a -2^(1/2)*a 0 0 -1/3*6^(1/2)*a -1/3*6^(1/2)*a 1/3*6^(1/2)*a 1/3*6^(1/2)*a 1/3*6^(1/2)*a 1/3*6^(1/2)*a -1/3*6^(1/2)*a