matlab绘制凸多面体,MATLAB计算多边形面积，比polyarea更加健壮

matlab自带了一个用于计算多边形面积的函数polyarea，但是这个函数使用起来比较让人纠结，它要求多边形的顶点必须指定为连续的，但在很多时候，这个是很难做到的，尤其是顶点数比较多的情况下。所以这个函数的健壮性就比较差。

比如，使用polyarea计算下面这些顶点定义的多边形的面积，得到50，这是正确的。-5 0

0 -5

5 0

0 5

但如果将其中间两个点交换后，再用polyarea计算，得到的面积是0。这显然是错误的。因为我只是交换了顶点的顺序，多边形并没有改变。

针对以上情况，我从另外的思路写了更加健壮的计算多边形面积的函数，实际上与计算多面体体积的函数非常类似。

代码 mypolyarea.mfunction area = mypolyarea(vertex)

% calculate the area of a polygon defined by the vertexes.

% Input:

%vertex: vertex of the polygon.

% Output:

%area: area of the polygon.

% 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

[~, area] = convexHull(dt);

用这个函数计算多边形面积时，对于点的顺序没有任何要求。需要注意，输入参数的格式略有不同，我个人不喜欢polyarea那种输入方式。

对于这个算法，虽然很简单，但是我还是想解释两点：1. DelaunayTri是用于将点集进行分解，得到很多三角形，对于点的顺序是没有任何要求的。

2. convexhull是一个凸包处理的函数，它可以通过对这些小三角形进行分析，得到最外围的凸包，并按着顺序对凸包边界上的点进行排列。当然，它同时也直接给出了凸包的体积或者面积。

如果你有更好的方法，请回复告诉我。谢谢。