matlab的多元函数_李志强.pptVIP

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);2 空间图形的绘制; M文件： t=0:0.01:60 =pi/60; =2; =3; x= *sin( )*cos( .*t).*t; y= *sin( )*sin( .*t).*t; z= *t.*cos( ); comet3(x,y,z); ; ; 空间曲线图的绘制 基本命令： 1 plot3——根据向量x,y,z绘制空间曲线图 例1：绘制参数方程组 的空间曲线图。 M文件： t=1:0.01:10*pi; x=sin(t); y=cos(t); z=t; plot3(x,y,z); grid on;%加网格线 ; 2 box on ——显示边框 t=1:0.01:10*pi; x=sin(t); y=cos(t); z=t; plot3(x,y,z,r-o); box on; %加边框 ; 空间网线图的绘制 基本命令：mesh（x,y,z）%绘制由二元函数z=f(x,y) 所确定的网格线图。 例2：绘制函数z=x^2+y^2在定义域 的 空间网线图。 x=linspace(0,4,30); （x=0:0.2:4;） y=linspace(0,3,20); ( y=0:0.15:3;) [X,Y]=meshgrid(x,y); %生成数据点阵 Z=X.^2+Y.^2; mesh(X,Y,Z); title(‘网线图’)； box on;%(加边框) ; 空间曲面图的绘制 基本命令：surf（x,y,z）%绘制由二元函数 z=f(x,y)所确定的空间曲面图。 例3：x=linspace(0,4,30); （x=0:0.2:4;） y=linspace(0,3,20); ( y=0:0.15:3;) [X,Y]=meshgrid(x,y); %生成数据点阵 Z=X.^2+Y.^2; surf(X,Y,Z); title(‘曲面图’)； box on； ; 例4：绘制空间曲面图： x=0:0.15:2*pi; y=0:0.15:2*pi; z=sin(y)*cos(x); surf(x,y,z) title(3-D surf); ;3 重积分的计算;2、长方体上三重积分数值计算函数triplequad 设三重积分 ，其中积分区域 是长方体域，则三重 积分函数triplequad的调用格式如下： 1）q = triplequad (fun,a,b,c,d,e,f)：输入量fun是被积函数 ； 2）q = triplequad (fun,a,b,c,d,tol)：同上，但限定精度tol； 3）q = triplequad (fun,a,b,c,d,tol, @quadl)：指定用quadl代替quad计算三重积分； 4）q = triplequad (fun,a,b,c,d,tol, @myquadf)：指定使用用户自己编写得求三重积分的Matlab程序myquadf.m代替quad计算三重积分.;解：M文件： syms x y; f=inter1(x,y); f=exp(-x.^2/2).*sin(x.^2+y); q = vpa(dblquad(@inter1,-2,2,-1,1),30) q = vpa(dblquad(@in