平面方程

點法式：
若平面E法向量n=(a,b,c)且過點A(x₀,y₀,z₀)，則平面E的方程式為a(x-x₀)+b(y-y₀)+c(z-z₀)=0。

一般式：
將方程式a(x-x₀)+b(y-y₀)+c(z-z₀)=0 化簡可得ax+by+cz+d=0的方程式。
我們將ax+by+cz+d=0稱為一般式。

一般式ax+by+cz+d=0的法向量為n=(a,b,c)

 (c)求平面方程：已知三点

// ------------- 由三点确定平面方程的个系数---------------- //

//参数: pDots ------------> 体数据中的三点

//               pdCoef -----------> 三点确定的平面方程系数指针

//返回值: 成功返回, 失败返回.

long  CMathFunCls::GetPanelCoef(double pPanelDots1[3],double pPanelDots2[3],double pPanelDots3[3],double pdCoef[4])

{

         long lRes =1;

         double      x21,x31,y21,y31,z21,z31;

 

         //计算中间量..---法向量

         x21=pPanelDots2[0]-pPanelDots1[0];//x21表示x2-x1,以下以此类推..

         x31=pPanelDots3[0]-pPanelDots1[0];

         y21=pPanelDots2[1]-pPanelDots1[1];

         y31=pPanelDots3[1]-pPanelDots1[1];

         z21=pPanelDots2[2]-pPanelDots1[2];

         z31=pPanelDots3[2]-pPanelDots1[2];

         //这里用点法式求解方程，最初自己理解的使用矩阵，是错误的。

         //已知三点坐标计算平面方程的公式(三点式)...

         //平面方程Ax+By+Cz+D=0

         pdCoef[0]=y21*z31-z21*y31;

         pdCoef[1]=z21*x31-x21*z31;

         pdCoef[2]=x21*y31-y21*x31;

         Normalise3D(pdCoef);//法向量单位化

 

         // 常数项

         pdCoef[3] = - pdCoef[0] * pPanelDots1[0] - pdCoef[1] * pPanelDots1[1] - pdCoef[2] * pPanelDots1[2];

 

         for(long i=0; i<4; i++) {

                   if(fabs(pdCoef[i])<1.0E-10)

                            pdCoef[i]=0.0;

         }

 

         return lRes;

}