matlab如何计算多面体体积

一个多面体的体积应该如何计算?

这本来是一个比较复杂的问题,应该要用到积分,但在matlab中确实很简单的事情,因为有现成的函数,只要稍微组合一下就可以。

下面我把函数写成了被调用的格式,输入参数为多面体的顶点坐标。

代码 vpolyhedron.m

function volume = vpolyhedron(vertex)
% calculate the volume of a polyhedron defined by the vertexes.
% Input:
%	vertex: vertex of the polyhedron.
% Output:
%	volume: volume of the polyhedron.
% Poweed by Xianbao Duan
% Email: xianbao.d@gmail.com
% Website: http://www.52souji.net/

% generate the Delaunay Triangle
dt = DelaunayTri(vertex);
% calculate the volumn of the body
[~, volume] = convexHull(dt);

下面举个例子具体讲解如何调用。

首先准备好多面体的顶点坐标文件 configure.in

48.2134 49.1482 48.9853
48.2134 49.1482 51.8421
48.0908 51.1528 50.4137
52.4593 49.7390 48.9853
50.3190 47.7426 48.9853
50.3926 49.8086 47.5569
52.4594 49.7390 51.8421
50.3927 49.8087 53.2706
50.3190 47.7426 51.8421
52.3392 47.7448 50.4137
50.6354 51.8668 48.9853
50.6354 51.8668 51.8421
52.7074 51.7431 50.4137
50.0000 50.0000 50.4000

然后,可以使用如下代码调用。

% load the vertexes of the polyhedron
p_coord = load('configure.in');

% calculate the volume
V = vpolyhedron(p_coord);

% display the result
disp(['The volume is: ',num2str(V)]);

程序的健壮性我不能保证,但我想一般的凸包多面体应该都可以吧。

如果有什么疑问或建议,给我留言。

标签: matlab

相关文章推荐

添加新评论 (无需注册,可直接评论)

已有 6 条评论

  1. qwe

    一看就不靠谱,都没有Face信息或者Edge信息

    1. 本人不是专业搞拓扑的,不过这个代码至少解决了我当时的问题。:)

  2. IceWorld

    不要在这里误导人了,你只输入点的坐标有什么用?没有多面体的拓扑关系,同样的点,连接形式不同,会得到完全不同的多面体。你居然还用DelaunayTri。。。无语,你这是在计算多面体的体积吗?你只是给了一些点,让matlab自己建立多面体然后计算。

  3. 学生路人

    请问 存入数据的.in是在哪建立的,我在matlab中找了半天也找不到

    1. 我爱搜集网博主

      这个需要你自己创建,就是你多面体各个顶点的坐标呀。

  4. 我爱搜集网博主

    如果是计算多边形的面积,matlab提供了直接的函数 polyarea