资源描述：
《matlab曲面绘图》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、泰州学院Matlab课堂练习十曲面绘图【实验目的】1．了解二元函数图形的制作。2．空间曲面等高线的制作。3．学习掌握MATLAB软件有关的命令。【实验内容】画出函数的图形,并画出其等高线。【实验准备】1．曲线绘图的MATLAB命令MATLAB中主要用mesh,surf命令绘制二元函数图形。mesh(x,y,z)画网格曲面，这里x,y,z是三个数据矩阵，分别表示数据点的横坐标，纵坐标和函数值，该命令将数据点在空间中描出,并连成网格。surf(x,y,z)画完整曲面,这里x,y,z是三个数据矩阵，分别表示数据点的横坐标
2、，纵坐标和函数值，该命令将数据点所表示曲面画出。可以用helpmesh,helpsurf查阅有关这些命令的详细信息【实验方法与步骤】练习1画出函数的图形,不妨将区域限制在。用MATLAB作图的程序代码为：clearx=-3:0.1:3;y=-3:0.1:3;[X,Y]=meshgrid(x,y);Z=exp(X.^2+Y.^2);mesh(X,Y,Z)21泰州学院Matlab课堂练习十结果如图5.1。图5.1是网格线图,如果要画完整的曲面图,只需将上述的MATLAB代码mesh(X,Y,Z)改为surf(X,Y,Z
3、),结果如图5.2图5.1锥面图5.2锥面要画等高线,需用contour,contour3命令.其中contour为二维等高线,contour3为三维等高线,如画图5.1的三维等高线,MATLAB代码为：clearx=-3:0.1:3;y=-3:0.1:3;21泰州学院Matlab课堂练习十[X,Y]=meshgrid(x,y);Z=sqrt(X.^2+Y.^2);contour3(X,Y,Z,10)xlable('X-axis'),ylable('Y-axis'),zlable('Z-axis')%z字符型加‘’
4、title('三维等高线')gridon21泰州学院Matlab课堂练习十结果如图5.3.图5.3等高线如画图5.1的二维等高线,MATLAB代码为：结果如图5.4.图5.4等高线如果要画的等高线，则用命令clearx=-3:0.1:3;y=-3:0.1:3;[X,Y]=meshgrid(x,y);Z=sqrt(X.^2+Y.^2);21泰州学院Matlab课堂练习十contour(X,Y,Z,[1,1])结果如图5.5。图5.5等高线练习1中，函数值可简单算出。在有些情况下，函数值不能简单算出。这是因为x和y的值
5、可能是非均匀间隔的甚至是随机分布的，也可能使用了不同的坐标系，比如非长方形的网。出现这些情况时,MATLAB中的函数griddata就用来产生经查值后的均匀间隔数据以作图。21泰州学院Matlab课堂练习十练习２二次曲面的方程如下讨论参数对其形状的影响。本练习的关键在于如何作出三维曲面图形，特别注意在给定值求时，若有开方运算，一是会出现虚数，二是对实数也有正负两个解。为了使虚数不出现在绘图中，采用了一种技巧，就是将虚数都换成非数(NaN).MATLAB代码为：a=input('a=');a=b=input('b='
6、);b=c=input('c=');c=d=input('d=');N=input('N=');xgrid=linspace(-abs(a),abs(a),N);ygrid=linspace(-abs(b),abs(b),N);[x,y]=meshgrid(xgrid,ygrid);z=c*sqrt(d-y.*y/b^2-x.*x/a^2)>>fork=2:N-1%以下7行程序的作用是取消z中含虚数的点>>forj=2:N-1ifimag(z(k,j))~=0z1(k,j)=0;end21泰州学院Matlab课堂练
7、习十ifall(imag(z([k-1:k+1],[j-1:j+1])))~=0za(k,j)=NaN;endendend>>surf(x,y,z1),holdon%画空间曲面>>ifu==1z2=-z1;surf(x,y,z2);%u=1时加画负半面axis([-abs(a),abs(a),-abs(b),abs(b),-abs(c),abs(c)]);end>>xlabel('x'),ylabel('y'),zlabel('z')>>holdoff运行程序，当时的结果见图5.6,当时的结果见图5.7,当时的结果
8、见图5.8,图5.6椭球面21泰州学院Matlab课堂练习十图5.7双曲面图5.8椭球双曲面练习3列出求空间两任意曲面的交线的程序。[x,y]meshgrid(-2:0.1:2,-2:0.1:2);z1=x.^2+y.^2z2=2*x-3*ymesh(x,y,z1);hold;mesh(x,y,z2);r0=(abs(z1-z2)<=0.1);zz=ro.*
