MATLAB绘制平行六面体

如果给出一个平行六面体(甚至其他多面体)的各个顶点坐标,如何画出这个平行六面体。

在网上找了找方法,可以参考这篇博客 matlab中patch函数详解。然后我具体查看了 Multifaceted Patches 帮助,记录下来以备后查。

绘制主要就是利用 patch 函数,patch 的一般调用格式为patch(X,Y,C)(二维)和patch(X,Y,Z,C)(三维)。这里主要看三维情况下。X,Y,Z分别是各个点的 X 坐标,Y 坐标和 Z 坐标,C 是颜色参数。假设我们有一个立方体如下图所示,

那么对应于这个立方体的XYZ参数就是下图所示,

可以看出,XYZ都是一个4*6的矩阵,3个矩阵第 n 列合起来表示一个平面的4个顶点的坐标,6个平面所以有6列,每个平面4个顶点所以有4行。但是这样表示就会出现一个问题——一些点重复出现,这是由于相邻两个平面之间有2个顶点是共用的。可以看到上图中实际上给出了24个顶点的坐标,但是真正“独立”的只有8个,也就是立方体实际的顶点数,很明显这样做不是一个明智的做法。于是,另一种方法出来了。

这种表示方法可以只提供8个顶点的坐标(Vertices 属性)而不是24个,给这些顶点编号1~8,然后指定每个平面上的4个顶点(Faces 属性,要按连接顺序),再指定颜色就可以了。顶点编号可以自己任意(最好有顺序,方便记忆理解),如下图,

然后平面和顶点的对应关系如下图所示,

Vertices每一行是一个顶点的坐标,共8行8个顶点;Faces每一行是一个平面的4个顶点,共6行6个平面。

下面举个例子。

假设我要通过其他方法已经得到平行六面体的8个顶点的分别是 (0,0,0),(1,0,0),(1,1,0),(0,1,0),(0,0.5,0),(1,0.5,1),(1,1.5,1),(0,1.5,1),代码如下,

clear
vertices_matrix = [ 0 0 0;1 0 0;1 1 0;0 1 0;0 0.5 1;1 0.5 1;1 1.5 1;0 1.5 1];
faces_matrix = [1 2 6 5;  2 3 7 6;3 4 8 7;4 1 5 8;1 2 3 4;5 6 7 8];
h = patch('Vertices',vertices_matrix,'Faces',faces_matrix,...
    'FaceVertexCData',hsv(8), ...
    'FaceColor','interp');
view(3)
xlabel('x', 'FontName', 'courier new', 'FontWeight', 'bold');
ylabel('y', 'FontName', 'courier new', 'FontWeight', 'bold');
zlabel('z', 'FontName', 'courier new', 'FontWeight', 'bold');
grid on;
box on; 

结果如下图所示,

在MATLAB中运行出来后可以通过旋转来从不同角度观察。

另外,假如flat_matrix我写成这样,

faces_matrix = [1 2 5 6;  2 3 7 6;3 4 8 7;4 1 5 8;1 2 3 4;5 6 7 8];

也就是说我把第一个面的顶点顺序稍微变一下,5和6交换,那么结果如下图,

明显第一个面不对,有部分“镂空”了,这是因为改变点的连接顺序,所以一定要注意顶点顺序。

每个面的颜色可以不一样,例如第一个面黄色,第二个面蓝色……,可以通过 hsv 函数或者指定 FaceVertexCData 属性值并设置 FaceColor 属性值为flat实现,具体可以参考本文开始给出的 Multifaceted Patches 链接和 Specifying Patch Coloring

OK,昨晚睡不着跟室友聊到2点,现在好困,不过还是不想睡。扯远了~这篇就这样了~

Good Luck,Alan!

本文参与腾讯云自媒体分享计划,欢迎正在阅读的你也加入,一起分享。

发表于

我来说两句

0 条评论
登录 后参与评论

相关文章

来自专栏智能计算时代

数据科学家用得最多的十种数据挖掘算法

Latest KDnuggets Poll asked Which methods/algorithms you used in the past 12 mon...

3567
来自专栏数据结构与算法

2039 骑马修栅栏

题目描述 Description Farmer John每年有很多栅栏要修理。他总是骑着马穿过每一个栅栏并修复它破损的地方。 John是一个与其他农民一样懒的人...

38111
来自专栏数据结构与算法

4768 跳石头

4768 跳石头 时间限制: 1 s 空间限制: 128000 KB 题目等级 : 黄金 Gold 题目描述 Description 一年一度的“...

27810
来自专栏趣学算法

数据结构 第17讲 沟通无限校园网——最小生成树(kruskal算法)

构造最小生成树还有一种算法,Kruskal算法:设G=(V,E)是无向连通带权图,V={1,2,…,n};设最小生成树T=(V,TE),该树的初始状态为只有n...

1732
来自专栏编程理解

数据结构(十二):最短路径(Dijkstra算法)

Dijkstra 算法使用贪心策略计算从起点到指定顶点的最短路径,通过不断选择距离起点最近的顶点,来逐渐扩大最短路径权值,直到覆盖图中所有顶点。

2932
来自专栏小樱的经验随笔

从零基础学三分查找

今晚是我们学长第二次讲课,讲了一个三分!认真听了一下,感觉不是很难,可能会比二分还简单些!我就把上课讲的内容归纳为一篇文章概述吧!以后也会重点讲解的! 简单点说...

44110
来自专栏机器之心

教程 | 如何使用TensorFlow中的高级API:Estimator、Experiment和Dataset

选自Medium 作者:Peter Roelants 机器之心编译 参与:李泽南、黄小天 近日,背景调查公司 Onfido 研究主管 Peter Roelant...

7187
来自专栏程序生活

Leetcode-Easy 70. Climbing Stairs

21. Merge Two Sorted Lists 描述: 有n阶楼梯,每步只能走1个或2个台阶,请问到达第n阶楼梯一共有多少走法? ? ...

3296
来自专栏云霄雨霁

子字符串查找----各种算法总结

2450
来自专栏数值分析与有限元编程

可视化 | Tecplot作三角形单元后处理工具

对于有限元分析的后处理,除了单的信息,还包括单元信息,比如一个单元由哪些结点组成。Tecplot可以处理的单元类型有三角形单元,四边形单元,四面体单元和六面体单...

4034

扫码关注云+社区

领取腾讯云代金券