MATLAB三维图形的绘制

bai_lan_ya 2024-06-28 13:05:04 阅读 80

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 all

x=-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 all

x=-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 all

f=@(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 all

f=@(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 all

f=@(x,y,z)x.^2/9-y.^2/4-z.^2-1;

fimplicit3(f)

5.3 抛物柱面

示例:y^2=x

clc,clear,close all

fsurf(@(y,z)y.^2)

5.4 椭圆锥面

示例:

clc,clear,close all

subplot(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 all

subplot(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 all

fsurf(@(x,y)x.*y)

5.7 椭圆柱面

clc,clear,close all

subplot(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 all

r=@(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 all

fmesh(@(x,y) erf(x)+cos (y),[-5,0,-5,5]), hold on

fmesh (@ (x,y)sin(x)+cos(y),[0,5,-5, 5]), hold off

7,四维数据曲面图

 



声明

本文内容仅代表作者观点,或转载于其他网站,本站不以此文作为商业用途
如有涉及侵权,请联系本站进行删除
转载本站原创文章,请注明来源及作者。