请问如何用matlab画出一个三维的正方体？

Posted 2023-05-10

请问如何用matlab画出一个三维的正方体？

最好能有解释

非常感谢！
我在后面加了这样一段：
boxplot3(1,2,3,5,5,5)
hold on
boxplot3(2.5,3.5,8,2,2,2)
hold on
boxplot3(3,4,10,1,1,1)
得到一个复杂一点的图形，不知怎么样才能给三个正方体染上不同的颜色呢？
如果更复杂一点，在他们的上表面加上一些大小不等的矩形、圆形、三角形等图形，请问该怎样写程序呢？

参考技术A 可以先编写个M文件,保存;

function boxplot3(x0,y0,z0,Lx,Ly,Lz)
%(x0,y0,z0)是第一个顶点的位置; (Lx,Ly,Lz)是长方体的长宽高.
x=[x0 x0 x0 x0 x0+Lx x0+Lx x0+Lx x0+Lx];
y=[y0 y0 y0+Ly y0+Ly y0 y0 y0+Ly y0+Ly];
z=[z0 z0+Lz z0+Lz z0 z0 z0+Lz z0+Lz z0];
index=zeros(6,5);
index(1,:)=[1 2 3 4 1];
index(2,:)=[5 6 7 8 5];
index(3,:)=[1 2 6 5 1];
index(4,:)=[4 3 7 8 4];
index(5,:)=[2 6 7 3 2];
index(6,:)=[1 5 8 4 1];
for k=1:6
plot3(x(index(k,:)),y(index(k,:)),z(index(k,:)))
hold on
end
保存为 boxplot3.m 文件后
在窗口中,调用

boxplot3(1,2,3,5,5,5)

就可得到边长为5的正方体.. (1,2,3)为第一个顶点的位置

补充程序:
编写两个M文件,并保存
(1)
function [xx,yy,zz,l]=planarsurface(p0,p1,p2)
v=p1-p0;
w=p2-p0;
s=0:0.2:1;
l=length(s);
[s,t]=meshgrid(s,s);
xx=p0(1)+s*v(1)+t*w(1);
yy=p0(2)+s*v(2)+t*w(2);
zz=p0(3)+s*v(3)+t*w(3);

(2)function boxsurface(p0,l)
[x,y,z]=planarsurface(p0,p0+[0 0 l(3)],p0+[0 l(2) 0]);
surf(x,y,z)
hold on
[x,y,z]=planarsurface(p0+[l(1) 0 0],p0+[l(1) 0 l(3)],p0+[l(1) l(2) 0]);
surf(x,y,z)
[x,y,z]=planarsurface(p0,p0+[0 0 l(3)],p0+[l(1) 0 0]);
surf(x,y,z)
[x,y,z]=planarsurface(p0+[0 l(2) 0],p0+[0 l(2) l(3)],p0+[l(1) l(2) 0]);
surf(x,y,z)
[x,y,z]=planarsurface(p0,p0+[l(1) 0 0],p0+[0 l(2) 0]);
surf(x,y,z)
[x,y,z]=planarsurface(p0+[0 0 l(3)],p0+[l(1) 0 l(3)],p0+[0 l(2) l(3)]);
surf(x,y,z)
axis equal
axis off

运行: boxsurface([1 1 1],[2 2 2])
%[1 1 1]表示初始位置.[2 2 2]: 边长本回答被提问者采纳

如何用Matlab快速画出带有3D渲染效果的复杂曲面

　  Matlab是一个很常用的理工科数学软件，我们平常会用它来画一些平面函数或者简单的曲线或者简单的3D平面图，但是通常很少用到它的稍微高级一点的画图功能。这里介绍一些高级渲染功能和画图技巧，先看结果。

写论文中可能会经常遇到这样的事情，想画一些高级点的漂亮的图不知道怎么画。当然我们通常可能会用solidworks, 3dsmax等等软件画这种3D的图，但是Matlab的好处是可以比较容易生成函数控制的复杂3D曲面形状，就比如上面的鸡蛋盒子一样的晶格图形，看似简单，用其他的软件还真的不太容易画。还有比如像下面这样的波浪图。想要画出这样的图，除了通常的一些画图命令外，还需要掌握一些小的tips，下面来逐一介绍。

 

1、 光照设置

以第二个图为例，它的Matlab代码如下

x=-6*pi:pi/100:6*pi;
y=x;
[X,Y]=meshgrid(x,y);
Z=5*(exp(-(X.^2+Y.^2)/100)+exp(-(X.^2+Y.^2)/4)).*cos(sqrt(X.^2+Y.^2));
figure;surf(X,Y,Z);shading flat; axis off;
view(15,45);axis([min(x),max(x),min(y),max(y),-6,12]);colormap(‘summer‘);

%光照渲染

light(‘position‘,[0,0,20],‘style‘,‘local‘,‘color‘,‘w‘);lighting phong;
light(‘position‘,[0,0,20],‘style‘,‘local‘,‘color‘,‘w‘);lighting phong;

前面四句是定义二元函数，meshgrid函数是将两个一维数组转化为二维数组以便定义二元函数，此时X,Y,Z都是二维矩阵。

figure;surf(X,Y,Z);是画出3D的表面图，但是这个图通常是带有网格的，而且由于定义的数据间隔小，网格也会很密，整张图看起来一片漆黑。因此需要用shading flat命令来去掉这些网格，此时显示的结果如下图所示，可以看到没什么光泽，不是很好看。上面的view命令是用来设置观察的角度的，axis([min(x),max(x),min(y),max(y),-6,12])设置图显示的区域大小，axis off意思是不显示坐标轴，colormap用来选择Matlab自带的一些色调，比如summer、pink、gray、winter等等。

 

最后要介绍的就是light命令了，在Matlab command 窗口输入 help light命令就可以查看light函数的所有性质，读者可以自行查询，这里只介绍它的功能。light(‘position‘,[0,0,20],‘style‘,‘local‘,‘color‘,‘w‘);lighting phong; 这段命令中，‘position‘是用来设置光源的位置，‘local‘表示光源是局域的光源而不是从无穷远处射来的平行光。lighting phong是设置光照的方式。用两次同样的light命令是提高光照的亮度，使图片看起来色泽更鲜亮。

 

2、着色设置

上面第一张的图中的球可以在极坐标下写出来，这里期待能够实现更复杂的着色控制，比如像下图这样的，所以代码稍微复杂些。

这个图的代码如下：

theta=-pi/2:pi/50:pi/2;
phi=0:pi/50:2*pi;
[Phi,Theta]=meshgrid(phi,theta);
x=R*cos(Theta).*cos(Phi);
y=R*cos(Theta).*sin(Phi);
z=R*sin(Theta);

%着色控制
color=zeros(length(theta),length(phi),3);
color(:,:,1) = cos(Theta)*cos(Theta0).*cos(Phi-Phi0)+sin(Theta)*sin(Theta0); % red
color(:,:,2) = sin(Theta)*sin(Theta0).*cos(Phi-Phi0)+cos(Theta)*cos(Theta0); % green
color(:,:,3) = sin(Theta)*sin(Theta0).*cos(Phi-Phi0); % blue
figure;
surf(x,y,z,color);shading flat;axis off;
view(15,45);axis([-1 1 -1 1 -1 1]);colormap(‘summer‘);
light(‘position‘,[0,0,20],‘style‘,‘local‘,‘color‘,‘w‘);lighting phong;
light(‘position‘,[-20,0,0],‘style‘,‘local‘,‘color‘,‘w‘);lighting phong;

前面到x,y,z的定义和上面一样，都是定义二元函数。这里唯一的不同点在于color变量的使用，我们可以自定义每个数据点的红、绿、蓝三个颜色的数值，就可以组合出任意的颜色。这样就几乎可以把这个图涂上任意的想要的颜色。曲面形状拓展一下，还可以画一个铁蒺藜的图。

 

3、材质设定

 Material命令是设置3D 表面的材料的，具体点说主要是反射光的性质，它的命令格式为 material([ka kd ks])；分别设置环境光、漫反射、镜面反射的反射程度，ka,kd,ks的值均在0到1之间。组合设置一下可以得到更细腻的3D表面的质感。

 

 从上图可以看到设置material后物体表面对光线的反射会更柔和点。

 

4、透明度设置

 

透明度的设置可以直接通过surf命令surf(x,y,z,color,‘FaceAlpha‘,0.3);来设置，这里的‘FaceAlpha‘就是设置面的透明度，数值在0到1之间，0表示完全透明，1表示完全不透明。

 

5、动画控制

最后是动画的控制，Matlab可以直接生成动画，动画的每一帧都是一个图，因此每个图都需要画出来。可以用getframe得到每一帧，用movie函数播放动画。下面有个完整的例子。

x=-6*pi:pi/20:6*pi;
y=x;
[X,Y]=meshgrid(x,y);
loops=200;
for j=1:loops
    Z=5*cos(j/20*pi-sqrt(X.^2+Y.^2)).*(exp(-(X.^2+Y.^2)/100)+exp(-(X.^2+Y.^2)/4));
    surf(X,Y,Z);axis off;shading flat;
    view(15,45);axis([min(x),max(x),min(y),max(y),-6,12]);colormap(‘summer‘);
    light(‘position‘,[0,0,20],‘style‘,‘local‘,‘color‘,‘w‘);lighting phong;
    light(‘position‘,[0,0,20],‘style‘,‘local‘,‘color‘,‘w‘);lighting phong;
    drawnow;
    F(j)=getframe;
end
movie(F,1,20); % 播放1次，每秒钟20张图

% 将动画存为一个视频文件

aviobj = VideoWriter(‘water wave.avi‘);
aviobj.FrameRate = 20;
open(aviobj);
writeVideo(aviobj,F);
close(aviobj);

动画的每一帧是存在F这个变量里的，然后用movie(F,1,20)来播放动画。首先需要用VideoWriter函数建立一个新的视频文件，这个视频文件在Matlab里的操作是一个对象aviobj，我们可以 设置Framerate为20，即每秒钟播放20帧。open(aviobj)和close(aviobj)是打开和关闭文件的操作，writeVideo(aviobj,F)是把F中的动画存在aviobj里，也就是存在water wave.avi文件里，matlab默认生成avi文件。之后打开water wave.avi文件就可以看到生成的动画了。

 

6、总结

 本文介绍了用Matlab画3D 物体和3D表面的一些方法。Matlab的优势是可以非常方便的用函数去构建复杂的3D物体和表面，使得可以非常简单地表达复杂的3D曲面，通过一些材质和光线渲染可以画出比较漂亮的图或者动画。其缺点是这样的图的数据量相对较大，所画的曲面不太容易导出为3D 模型，因此Matlab画这种图的主要作用可能是在科研上，而非画模型上。当然一种解决方案是，直接用Matlab写一个程序直接生成3D模型的数据文件，或许可以实现把Matlab的编程优势与3D建模软件的优势结合在一起，实现更复杂的3D建模。