2.3多元函数微积分学教程.ppt

2.3 多元函数微积分学;2.3.1 多元函数的微分;;syms x y dx dy f df f=sin(x*y)+(cos(x^3+y^2))^2; fx=diff(f,x) %求f x fy=diff(f,y) %求f y df=fx*dx+fy*dy %求全微分 f2x2=diff(fx,x) %求f’’ xx f2xy=diff(fx,y) %求f’’ xy f3xyx=diff(f2xy,x) %求f’’’ xyx;M文件： syms x y z; f=x^2+y^2+z^2-4*z; fx=diff(f,x);fz=diff(f,z);fy=diff(f,y); zx=collect(-simple(fx/fz)); z2xx=collect(diff(zx,x)+diff(zx,z)*zx);解：求其Jacobi矩阵，建立Matlab的M文件程序： syms u v; % x=exp(u)*cos(v);y=exp(u)*sin(v); x=exp(u)*cos(v);y=exp(u)*sin(v); J=jacobian([x,y],[u,v]);2.3.2 极值问题; ;解：首先求解偏导数 syms x y; f=x^3-y^3+3*x^2+3*y^2-9*x; diff(f,x),diff(f,y);第四步，分别判别四点是否为极值，建立M文件， 自动判断四点的极值情况： xx=[1 -3 1 -3]; %驻点横坐标 yy=[0 0 2 2]; %驻点纵坐标 for i=1:4 D=(6*xx(i)+6)*(-6*yy(i)+6) if D0 if (6*xx(i)+6)0 x=xx(i) y=yy(i) disp(为极大值点；) disp(极大值为) fmax=x^3-y^3+3*x^2+3*y^2-9*x elseif (6*xx(i)+6)0 x=xx(i) y=yy(i) ;disp(为极小值点；) disp(极小值为) fmin=x^3-y^3+3*x^2+3*y^2-9*x end end if D0 x=xx(i) y=yy(i) disp(该点不是极值点；) end if D==0 x=xx(i) y=yy(i) disp(无法确定！) end end ;syms x y z r k; L=8*x*y*z+k*(x^2+y^2+z^2-r^2); diff(L,x),diff(L,y),diff(L,z),diff(L,k) 再求解驻点坐标： [k,x,y,z]=solve(8*y*z+2*k*x=0,8*x*z+2*k*y=0, 8*x*y+2*k*z=0,x^2+y^2+z^2-1=0);2.3.3 空间曲面的绘制;2）mesh(x,y,z)：用空间中的两组相交的平行平面上的网状线的方式表示曲面，即绘制网格图. 3) meshc(x,y,z)：绘制网格图，并带有等高线. 4) meshz(x,y,z)：屏蔽的网格图. 5) surf(x, y, z)：用空间中网状线与网格中填充色彩的方式表示曲面，即绘制表面图. 6) surfc(x, y, z)：绘制表面图，并带有等高线. 7) surfl (x, y, z)：绘制表面图，并带有阴影.; ;3、二元函数简捷绘图指令 1) 绘制网格图函数ezmesh ① ezmesh(f)：生成二元函数f(x, y) 网格图，其中f为一个字符串. ② ezmesh(f,domain)：绘制二元函数f(x, y) 在指定邻域上的网格图，该邻域可以是4×1 向量[xmin, xmax, ymin, ymax] 或 2×1 向量[min, max] (其中, min x max, min y max). ③ ezmesh(x, y, z)：绘制在默认邻域-2 s 2, -2 t 2上的参数方程x = x(s,t)，y = y(s,t)和z = z(s,t)确定的曲面. ④ ezmesh(x, y, z, [smin, smax, tmin, tmax]) or ezmesh(x, y, z, [min, m