Matlab判断空间点是否在三角形内部

        在几何相关的数学建模问题中，会遇到已知三角形顶点空间坐标，判断空间任意点是否在三角形内部的情况，本文给出了一种判断空间点是否在三角形内的简单算法，并通过Matlab进行实现，同时支持矩阵运算。

        算法基本思想很简单，如图所示，和平面内任意点P，若P落在

内部 ，则必有；反之，若P落在

外部， 则必有。本算法即基于此基本原理进行实现，并可以支持输入为n*3的矩阵坐标集合。

 

        总代码如下所示：

% function result=Is_in_triangle(p,a,b,c)
% 计算点是否位于三角形内的函数,p可为n*3矩阵，对应输出result为n维向量
function result=Is_in_triangle(p,a,b,c)
% 通过面积计算是否在三角形内
Sabc=0.5*sqrt(sum((cross(b-a,c-a,2).^2)')');
s1=0.5*sqrt(sum((cross(a-p,b-p,2).^2)')');
s2=0.5*sqrt(sum((cross(a-p,c-p,2).^2)')');
s3=0.5*sqrt(sum((cross(c-p,b-p,2).^2)')');
result = (abs(Sabc-s1-s2-s3)<=0.0001);
end