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 条评论
登录 后参与评论

相关文章

来自专栏小二的折腾日记

LeetCode-55-Jump-Game

由题可知,数组的位置表示从该位置可以像前跳的步数,看最终能否跳到结尾。乍一看,这像是一个动态规划的问题,dp数组内存储每一个位置能够走的最远的位置,但是仔细一想...

1003
来自专栏C语言及其他语言

[每日一题]1453: [蓝桥杯][历届试题]翻硬币

题目描述 小明正在玩一个“翻硬币”的游戏。 桌上放着排成一排的若干硬币。我们用 * 表示正面,用 o 表示反面(是小写字母,不是零)。 比如,可...

26911
来自专栏云霄雨霁

加权无向图----Prim算法实现最小生成树

2430
来自专栏简书专栏

Python数据科学库-小测验

答:np.arange、np.array、np.ones、np.zeros、np.full

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

一个例子教你如何与出题人斗智斗勇

我以前出过一道题,卡了10种贪心,但还是被第11种贪心A了,  一道题不会做?贪嘛,能怎么贪怎么贪,想怎么贪怎么贪! 现在NOIP题目的数据给的不...

2646
来自专栏算法channel

除自身累乘算法题,又有创意解法了

一个数组,求除了某元素自身位置之外的其他元素累积相乘,返回一个同长度的数组。有两个要求比较苛刻: 1) 不能用除法 2) 时间复杂度O(n),空间复杂度O(1)...

940
来自专栏null的专栏

挑战数据结构和算法——栈的push、pop序列

题目来源“数据结构与算法面试题80道”。在此给出我的解法,如你有更好的解法,欢迎留言。 ? 问题分析:本题考查栈的基本操作,栈是一种“先进后出”的数据结构。判...

2966
来自专栏HansBug's Lab

算法模板——KMP字符串匹配

功能:输入一个原串,再输入N个待匹配串,在待匹配串中找出全部原串的起始位置 原理:KMP算法,其实这个东西已经包含了AC自动机的思想(fail指针/数组),只不...

3315
来自专栏数据魔术师

运筹学教学 | 十分钟教你求解分配问题(assignment problem)

biu~ biu~ biu~ 我们的运筹学教学推文又出新文拉 还是熟悉的配方,熟悉的味道 今天向大家推出的是 运筹学教学--第六弹 分配问题(Assignmen...

4948
来自专栏计算机视觉与深度学习基础

Leetcode 162 Find Peak Element

A peak element is an element that is greater than its neighbors. Given an inpu...

1927

扫码关注云+社区