MATLAB绘制二元函数图形资料下载.pdf

1、3;%x 的范围为-3，3y=-3:%y 的范围为-3，3X,Y=meshgrid（x,y）;%将向量 x,y 指定的区域转化为矩阵 X,YZ=sqrt（X.2+Y.2）;%产生函数值 Zmesh（X,Y,Z）运行结果为图 5.3如果画等高线，用 contour，contour3 命令。contour 画二维等高线。contour3 画三维等高线。画图 5.3 所示的三维等高线的 MATLAB代码为clear;y=-3:X,Y=meshgrid（x,y）;Z=sqrt（X.2+Y.2）;contour3（X,Y,Z,10）;%画 10 条等高线xlabel（X-axis）,ylabel（Y-a

2、xis）,zlabel（Z-axis）;%三个坐标轴的标记title（Contour3 of Surface）%标题grid on%画网格线运行结果为图 5.4如果画图 5.4 所示的二维等高线，相应的 MATLAB 代码为clear;contour（X,Y,Z,10）;xlabel（X-axis）,ylabel（Y-axis）;title（Contour3 of Surface）grid on运行结果为如果要画 z=1 的等高线，相应的 MATLAB 代码为clear;contour（X,Y,Z,1 1）运行结果为练习练习 2 二次曲面的方程如下222222xyzdabc讨论参数 a,b,c

3、 对其形状的影响。相应的 MATLAB 代码为a=input（a=）;b=input（b=）;c=input（c=）;d=input（d=）;N=input（N=）;%输入参数，N 为网格线数目xgrid=linspace（-abs（a）,abs（a）,N）;%建立 x 网格坐标ygrid=linspace（-abs（b）,abs（b）,N）;%建立 y 网格坐标x,y=meshgrid（xgrid,ygrid）;%确定 NN 个点的 x,y 网格坐标z=c*sqrt（d-y.*y/b2-x.*x/a2）;u=1;%u=1,表示 z 要取正值z1=real（z）;%取 z 的实部 z1for

4、k=2:N-1;%以下 7 行程序的作用是取消 z 中含虚数的点for j=2:N-1if imag（z（k,j）=0 z1（k,j）=0;endif all（imag（z（k-1:k+1,j-1:j+1）=0 z1（k,j）=NaN;endendendsurf（x,y,z1）,hold on%画空间曲面if u=1 z2=-z1;surf（x,y,z2）;%u=1 时加画负半面axis（-abs（a）,abs（a）,-abs（b）,abs（b）,-abs（c）,abs（c）;endxlabel（x）,ylabel（y）,zlabel（z）hold off运行程序，当 a=5,b=4,c=3,

5、d=1,N=50 时结果为当 a=5i,b=4,c=3,d=1,N=15 时结果为当 a=5i,b=4i,c=3,d=0.1,N=10 时结果为【练习与思考】【练习与思考】1.画出空间曲面22221sin10yxyxz在30,30yx范围内的图形，并画出相应的等高线。解：clear;close;u=-30:0.5:30;v=-30:x,y=meshgrid（u,v）;z=10*sin（sqrt（x.2+y.2）./sqrt（1+x.2+y.2）;subplot（1,2,1）;mesh（x,y,z）subplot（1,2,2）;contour（x,y,z,10）2.根据给定的参数方程，绘制下列曲

6、面的图形。a）椭球面vuxsincos3，vuycoscos2，uzsin；u=-4:4;v=-4:U,V=meshgrid（u,v）;X=3.*cos（U）.*sin（V）;Y=2.*cos（U）.*cos（V）;Z=sin（U）;surf（X,Y,Z）;axis equalb）椭圆抛物面vuxsin3，vuycos2，24uz；X=3.*U.*sin（V）;Y=2.*U.*cos（V）;Z=4.*U.2;shading interp;colormap（hot）;axis equalc）单叶双曲面vuxsinsec3，vucossec2，uztan4；u=-2:2;v=-2:X=3*sec（

7、U）.*sin（V）;Y=2*sec（U）.*cos（V）;Z=4*tan（U）;mesh（X,Y,Z）;colormap（jet）;axis equald）双叶抛物面ux，vy，322vuz；X=U;Y=V;Z=（U.2-V.2）/3;axis equale）旋转面vuxsinln，vuycosln，uz；X=log（U）.*sin（V）;Y=log（U）.*cos（V）;Z=U;axis equalf）圆锥面vuxsin，vuycos，uz；X=U.*sin（V）;Y=U.*cos（V）;axis equalg）环面vuxcos）cos4.03（，vuysin）cos4.03（，vzsin

8、4.0；X=（3+0.4*cos（U）.*cos（V）;Y=（3+0.4*cos（U）.*sin（V）;Z=0.4*sin（V）;axis equalh）正螺面vuxsin，vuycos，vz4。Z=4*V;axis equal3.在一丘陵地带测量高程，x和y方向每隔 100 米测一个点，得高程见表 5-2，试拟合一曲面，确定合适的模型，并由此找出最高点和该点的高程。表 5-2 高程数据yx100200300400100200636698697712624630478478300400680662674626598552412334解：x=100 100 100 100 200 200 200

9、 200 300 300 300 300400 400 400 400;y=100 200 300 400 100 200 300 400 100 200 300 400100 200 300 400;z=636 697 624 478 698 712 630 478 680 674 598 412662 626 552 334;xi=100:5:400;yi=100:X,Y=meshgrid（xi,yi）;H=griddata（x,y,z,X,Y,cubic）;surf（X,Y,H）;view（-112,26）;hold on;maxh=vpa（max（max（H）,6）r,c=find（H=single（maxh）;stem3（X（r,c）,Y（r,c）,maxh,fill）