1,三维曲线
plot3(x,y,z)通过描点连线画出曲图,这里x,y,z都是n维向量,分别表示该曲线上点集的横坐标,纵坐标,竖坐标。
t=0:pi/50:10*pi;subplot(121),plot3(sin(t),cos(t),t);subplot(122),fplot3(@(t)sin(t),@(t)cos(t),@(t)t,[0,10*pi])
2,网格图
命令mesh(x,y,z)画网格曲图。这里x,y,z是三个同维数的数据矩阵,分别表示数据点的横坐标,纵坐标,竖坐标,命令mesh(x,y,z)将该数据点在空间中描出u,并且连成网格
示例:绘制二元函数z=sin(xy)/xy
clc,clear,close allx=-5:0.2:5;[x,y]=meshgrid(x);%生成网格数据z=(sin(x./y)+eps)./(x.*y+eps);%为避免0/0,分子分母都加eps,变成浮点精度型subplot(121),mesh(x,y,z)subplot(122),fmesh(@(x,y)sin(x.*y)./(x.*y))
3,表面图
命令surf(x,y,z)画三维,这里x,y,z是三个同维数的数据矩阵,分别表示数据点的横坐标,纵坐标,竖坐标
示例:绘制二元函数z=sin(xy)/xy
clc,clear,close allx=-5:0.2:5;[x,y]=meshgrid(x);z=(sin(x./y)+eps)./(x.*y+eps);%为避免0/0,分子分母都加eps,变成浮点精度型subplot(121),surf(x,y,z)subplot(122),fsurf(@(x,y)sin(x.*y)./(x.*y))
4,旋转曲图
示例:
使用命令fmesh画图,matlab程序如下:
clc,clear,close allf=@(x,y,z)x.^2+(sqrt(y.^2+z.^2)-5).^2-16;subplot(121),fimplicit3(f,[-4,4,-9,9,-9,9])x=@(u,v)4*cos(u);y=@(u,v)(5+4*sin(u)).*cos(v);z=@(u,v)(5+4*sin(u)).*sin(v);subplot(122),fsurf(x,y,z)
5,二次曲面图形
matlab中使用绘图命令fmesh或fsurf,画显函数或参数方程表示的二次曲面很方便
或者直接用fimlicit3画隐函数表示的二次曲面
5.1 单页双曲面
示例:
clc,clear,close allf=@(x,y,z)x.^2/4+y.^2/10-z.^2/8-1;subplot(121),fimplicit3(f,[-20,20,-20,20,-15,15])x=@(u,v)2*cosh(v).*cos(u);y=@(u,v)sqrt(10)*cosh(v).*sin(u);z=@(u,v)2*sqrt(2)*sin(v);subplot(122),fmesh(x,y,z,[0,2*pi,-pi,pi]);
5.2 双叶双曲面
示例:
clc,clear,close allf=@(x,y,z)x.^2/9-y.^2/4-z.^2-1;fimplicit3(f)
5.3 抛物柱面
示例:y^2=x
clc,clear,close allfsurf(@(y,z)y.^2)
5.4 椭圆锥面
示例:
clc,clear,close allsubplot(121)fimplicit3(@(x,y,z)x.^2/9+y.^2/4-z.^2,[-6,6,-4,4,-2,2])subplot(122),x=@(s,t)3*tan(s).*cos(t);y=@(s,t)2*tan(s).*sin(t);z=@(s,t)tan(s);fsurf(x,y,z,[-1,1,0,2*pi])
5.5 椭圆球
clc,clear,close allsubplot(121),fimplicit3(@(x,y,z)x.^2/9+y.^2/9+z.^2/6-1)subplot(122),ellipsoid(0,0,0,3,2,sqrt(6))
5.6 马鞍面
示例:z=xy
clc,clear,close allfsurf(@(x,y)x.*y)
5.7 椭圆柱面
clc,clear,close allsubplot(121),fimplicit3(@(x,y,z)x.^2/9+y.^2/4-1)x=@(u,v)3*cos(u);y=@(u,v)2*sin(u);z=@(u,v)v;subplot(122),fsurf(x,y,z)
6,其他曲面图形
clc,clear,close allr=@(s,t)2+sin(7*s+5*t);x=@(s,t)r(s,t).*cos(s).*sin(t);y=@(s,t)r(s,t).*sin(s).*cos(t);z=@(s,t)r(s,t).*cos(t);fmesh(x,y,z,[0,2*pi,0,pi]),alpha(0.8)
clc,clear,close allfmesh(@(x,y) erf(x)+cos (y),[-5,0,-5,5]), hold onfmesh (@ (x,y)sin(x)+cos(y),[0,5,-5, 5]), hold off
7,四维数据曲面图