微积分经济类课程侯亚君第6章节多元函数微分学数学实验6演示文档.pptVIP

* 数学实验 6 二元函数的图像、多元 函数的偏导数与极值 首页 上页 下页 返回 结束 1.二元函数图像的绘制 mesh(x,y,z) ： 画网格曲面，这里x，y，z是三个 来绘制二元函数图像． 本实验将介绍如何用MATLAB中的命令mesh，surf 据矩阵，分别表示数据点的横坐标，纵坐标和函数值， 该命令将数据点在空间中描出,并连成网格． 首页 上页 下页 返回 结束 surf(x,y,z) ： 画完整曲面,这里x,y,z是三个数 据矩阵，分别表示数据点的横坐标，纵坐标和函数值， 画出函数 该命令将数据点所表示曲面画出． 可以用help mesh， help surf查阅有关这些命令 的详细信息． 例6-36 的图形,不妨将区 域限制在 解 在命令窗口中输入语句： 首页 上页 下页 返回 结束 clear ; x=-3:0.1:3; y=-3:0.1:3; [X,Y]=meshgrid(x,y); Z=sqrt(X.^2+Y.^2); mesh(X,Y,Z) x的范围为[-3,3] y的范围为[-3,3] 将向量x，y指定的区域转化为矩阵X,Y 产生函数值Z 首页 上页 下页 返回 结束 运行结果： 上图是网格线图，如果要画完整的曲面图,只需将 上述的mesh(X,Y,Z)改为surf(X,Y,Z)． 首页 上页 下页 返回 结束 运行结果为： 要画等高线,需用contour,contour3命令．其中 contour为二维等高线, contour3为三维等高线． 首页 上页 下页 返回 结束 为画出锥面网格线图的三维等高线，输入语句： clear ; x=-3:0.1:3; y=-3:0.1:3; [X,Y]=meshgrid(x,y); Z=sqrt(X.^2+Y.^2); contour3(X,Y,Z,10) xlabel(X-axis),ylabel(Y-axis),zlabel(Z-axis) title(Contour3 of Surface) grid on 画10条等高线 三个坐标轴的标记 标题 画网格线 首页 上页 下页 返回 结束 运行结果： 首页 上页 下页 返回 结束 为画出锥面网格线图的二维等高线，输入语句： clear; x=-3:0.1:3; y=-3:0.1:3; [X,Y]=meshgrid(x,y); Z=sqrt(X.^2+Y.^2); contour (X,Y,Z,10) xlabel(X-axis),ylabel(Y-axis) title(Contour of Surface) grid on 首页 上页 下页 返回 结束 运行结果： 首页 上页 下页 返回 结束 如果要画z=1的等高线，则输入命令： clear; x=-3:0.1:3; y=-3:0.1:3; [X,Y]=meshgrid(x,y); Z=sqrt(X.^2+Y.^2); contour(X,Y,Z,[1 1]) 运行结果： 首页 上页 下页 返回 结束 例6-36中，函数值 可简单算出． 在有些情况下，函数值z=f (x,y)不能简单算出． 这是 因为x和y的值可能是非均匀间隔的甚至是随机分布的， 也可能使用了不同的坐标系，比如非长方形的网． 出现这些情况时,MATLAB中的函数griddata就用 来产生经查值后的均匀间隔数据以作图． 例6-37 二次曲面的方程如下 讨论参数a，b，c 对其形状的影响． 首页 上页 下页 返回 结束 本例的关键在于如何作出三维曲面图形， 分析： 特别注意在给定x，y值求z时，若有开方运算， 一是会出 现虚数， 二是对实数也有正负两个解． 为使虚数不出现 在图中，采用一种技巧，就是将虚数都换成非数(NaN)． 解 在命令窗口中输入语句： a=input(a=); b=input(b=); c=input(c=); d=input(d=); N=input(N=); xgrid=linspace(-abs(a), abs(a),N); ygrid=linspace(-abs(b), abs(b),N); 输入参数，N为网格线数目 建立x网格坐标 建立y网格坐标 首页 上页 下页 返回 结束 [x,y]=meshgrid(xgrid,ygrid); z=c*sqrt(d-y.*y/b^2-x.*x/a^2); u=1; z1=real(z); for k=