判断三维空间中的一个点是否在三角形内，边上的一种算法

   假设三角形的三点用A（Xa,Ya,Za）,B(Xb,Yb,Zb),C(Xc,Yc,Zc)表示，判断P(Xp,Yp,Zp)与三角形ABC的关系

网上判断点是否处于三角内内部，边上的方法有很多，但是比较简洁的方法，比如面积法，即判断三角形PAB,PAC,PBC之和与三角形ABC的面积。如果S(PAB+PAC+PBC)=S(ABC),则点位于三角形内部，或边上。但是上述方法由于涉及到平方和开根号，计算量比较大。

   最近自己整理了一个计算量稍小的算法。具体如下：

令向量（B-A）=（X1,Y1,Z1）；向量(C-A)=(X2,Y2,Z2)；向量（P-A）=（X3,Y3,Z3）；

假设点P位于三角形ABC所在的平面，则向量（P-A）可以用向量（B-A）和（C-A）表示，即（P-A）=m*(B-A)+n(C-A);只要求出m,n的值即可判断点与三角形的关系。如果没有解，那么p不在三角形所在平面，如果有解，那么：                            若：          若：         若：

   m>0                         m=0;          n=0;         m>0;

   n>0;                        n>0;          m>0;         n>0;

   m+n<1                       m+n<1;         m+n<1;       m+n=1;

则P在三角形内部。               P在AC边上   P在AB边上      P在BC边上

 

所以，问题的关键是要计算m,n的值，从而可以得到点和三角形的空间关系。

具体步骤如下：

由（P-A）=m*(B-A)+n(C-A)得到（X3,Y3,Z3）=m*（X1,Y1,Z1）+n*（X2,Y2,Z2），写成矩阵形式，得到：

X1 X2   m    X3

Y1 Y2*    =  Y3

Z1 Z2   n    Z3

          X1X2                   X3

令矩阵K为  Y1Y2               b为 Y3

          Z1Z2                   Z3

      因为ABC为三角形，所以向量（B-A）=（X1,Y1,Z1）；向量(C-A)=(X2,Y2,Z2)线性不相关，即矩阵K为列满秩矩阵，所以

    m

        =K的加号逆*b;

    n

其中K的加号逆为