1.基本原理
已知三点p1,p2,p3,求其平面方程,p1=[x1,y1,z1];p2=[x2,y2,z2];p3=[x3,y3,z3];那么求解下面矩阵q行列式就是了:q=[ones(4,1),[x,y,z];p1;p2;p3];detb=det(q);最后令q=0;
这里的求解方法是四阶矩阵的行列式求解,可以从三阶矩阵的行列式推导出来,那就是常规的求解方法:方程组写成[p1;p2;p3]=D,其中D=[d,d,d]'表示平面方程的常数项,先求q=[p1;p2;p3]的行列式,再分别求D与三个坐标两个分量合成的三阶矩阵的行列式,设为q1,q2,q3;那么所求的平面方程的四个参数就是A=q1/q;B=q2/q;c=q3/q;D=d;同时乘以q/d所得结果本质上就是上面四阶矩阵行列式解。
2.Matlab代码
<span sty