数学分析实验..docVIP

多变量函数的极限 若想求出二元函数的极限： ，则可以嵌套使用函数： 使用公式为： 如果不是确定的值，例如是，则上述极限的求取顺序不能交换。 例1.1: 试求出二元函数的极限值： Matlab 执行代码如下： syms x y a; f=exp(-1/(y^2+x^2))*sin(x)^2/x^2*(1+1/y^2)^(x+a^2*y^2); L=limit(limit(f,x,1/sqrt(y)),y,inf) 输出结果: L = exp(a^2) 多元函数的偏导数 例2.1 试求出二元函数的偏导数，并用图形表示。 Matlab 代码如下： syms x y z=(x^2-2*x)*exp(-x^2-y^2-x*y); zx=simple(diff(z,x)) zy=diff(z,y) 执行结果： zx = -exp(-x^2-y^2-x*y)*(-2*x+2+2*x^3+x^2*y-4*x^2-2*x*y) zy = (x^2-2*x)*(-2*y-x)*exp(-x^2-y^2-x*y) 在区域区域内生成网格，则可以分别得出原函数及偏导数的数值解： [x,y]=meshgrid(-3:0.2:3,-2:0.2:2);%产生网格数据 z=(x.^2-2*x).*exp(-x.^2-y.^2-x.*y); surf(x,y,z) %绘制曲面 如图1 contour(x,y,z,30),hold on %绘制等高线 zx=-exp(-x.^2-y.^2-x.*y).*(-2*x+2+2*x.^3+x.^2.*y-4*x.^2-2*x.*y); zy=(x.^2-2.*x).*(-2*y-x).*exp(-x.^2-y.^2-x.*y);%求偏导的数值解 quiver(x,y,zx,zy)%绘制引力线（向量场） 图1 图2 偏导数等值线图及引力线图 多元函数的Jacobic矩阵 调用格式： 其中x是自变量构成的变量，y是各函数构成的向量 例3.1 设 （教材P156页例3实例验证） ，试求其Jacobic矩阵 解：Matlab执行代码如下： syms r theta phi x=r*sin(theta)*cos(phi); y=r*sin(theta)*sin(phi); z=r*cos(theta); J=jacobian([x y z],[r theta phi]) simple(det(J))%????Jabcobian???ó??ó|μ?DDáDê?μ??μ?￡ 运行结果： J = [ sin(theta)*cos(phi), r*cos(theta)*cos(phi), -r*sin(theta)*sin(phi)] [ sin(theta)*sin(phi), r*cos(theta)*sin(phi), r*sin(theta)*cos(phi)] [ cos(theta), -r*sin(theta), 0] ans = sin(theta)*r^2 隐函数的偏导数 已知隐函数的数学表达式为，则可通过下面的公式求出 ,用Matlab语言可如下计算： 例4.1 二元函数 解：程序代码： syms x y; f=(x^2-2*x)*exp(-x^2-y^2-x*y); pretty(-simple(diff(f,x)/diff(f,y))) 运算结果： ans = (-(2*x-2)*exp(-x^2-y^2-x*y)-(x^2-2*x)*(-2*x-y)*exp(-x^2-y^2-x*y))/(x^2-2*x)/(-2*y-x)/exp(-x^2-y^2-x*y) pretty(-simple(diff(f,x)/diff(f,y))) 3 2 2 -2 x + 2 + 2 x + x y - 4 x - 2 x y - ------------------------------------- x (x - 2) (2 y + x) 多变量函数的T aylor幂级数展开式 Matlab符号工具箱并未提