已知3点坐标,求平面ax+by+cz+d=0;
得出的abcd均在等式的一侧,如上表达式
void get_panel(TPoint p1,TPoint p2,TPoint p3,double &a,double &b,double &c,double &d)
{
a = ( (p2.y-p1.y)*(p3.z-p1.z)-(p2.z-p1.z)*(p3.y-p1.y) );
b = ( (p2.z-p1.z)*(p3.x-p1.x)-(p2.x-p1.x)*(p3.z-p1.z) );
c = ( (p2.x-p1.x)*(p3.y-p1.y)-(p2.y-p1.y)*(p3.x-p1.x) );
d = ( 0-(a*p1.x+b*p1.y+c*p1.z) );
}
这样的:
设 已知三点A(x1,y1,z1),B(x2,y2,z2),C(x3,y3,z3)
任意找在这个面的两个不平行的向量,
BA=(x1-x2,y1-y2,z1-z2)=(v1[0],v1[1],v1[2])
CB=(x2-x3,y2-y3,z2-z3)=(v2[0],v2[1],v2[2])
法向量为同时垂直于这两个向量的一个向量。
利用叉乘可以直接得到
n = BA×CB
=[ i , j , k ]
[v1[0],v1[1],v1[2]]
[v2[0],v2[1],v2[2]]
= (v1[1]*v2[2]-v1[2]*v2[1],v1[2]*v2[0]-v1[0]*v2[2],v1[0]*v2[1]-v1[1]*v2[0])
本文代码段为CSDN博主「S小超」,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接及本声明。
原文链接:https://blog.csdn.net/qq_48240095/article/details/123224513