三点求平面方程式的算法

// 点

//

//     P1(X1,Y1,Z1), P2(X2,Y2,Z2), P3(X3,Y3,Z3)

//

// 设法向量为n

//

//     P2-P1  [（X2-X1),(Y2-Y1),(Z2-Z1)]

//     P3-P1  [（X3-X1),(Y3-Y1),(Z3-Z1)]

//

// 叉乘

//          |                                     |

//           |        i         j         k          |

//    n =  |  X2-X1   Y2-Y1   Z2-Z1   |

//           |  X3-X1   Y3-Y1   Z3-Z1   |

//          |                                      |

//

// 得

//     [(Y2-Y1)*(Z3-Z1) - (Z2-Z1)*(Y3-Y1)]*i

//     [(Z2-Z1)*(X3-X1) - (X2-X1)*(Z3-Z1)]*j

//     [(X2-X1)*(Y3-Y1) - (X3-X1)*(Y2-Y1)]*k

//

//     n = [(Y2-Y1)*(Z3-Z1) - (Z2-Z1)*(Y3-Y1)]*i + [(Z2-Z1)*(X3-X1) - (X2-X1)*(Z3-Z1)]*j + [(X2-X1)*(Y3-Y1) - (X3-X1)*(Y2-Y1)]*k

//

// 代入点法式

// 

//     [(Y2-Y1)*(Z3-Z1) - (Z2-Z1)*(Y3-Y1)]*(X-X1) +

//     [(Z2-Z1)*(X3-X1) - (X2-X1)*(Z3-Z1)]*(Y-Y1) +

//     [(X2-X1)*(Y3-Y1) - (X3-X1)*(Y2-Y1)]*(Z-Z1) = 0

//

//  化简

//

//     (Y2*Z3 - Y2*Z1 - Y1*Z3 - Y3*Z2 + Y1*Z2 + Y3*Z1)*(X-X1) +

//

//     (X3*Z2 - X1*Z2 - X3*Z1 - X2*Z3 + X2*Z1 + X1*Z3)*(Y-Y1) +

//

//     (X2*Y3 - X2*Y1 - X1*Y3 - X3*Y2 + X3*Y1 + X1*Y2)*(Z-Z1) = 0
// 

// 化简

//

//     X*Y2*Z3 - X*Y2*Z1 - X*Y1*Z3 - X*Y3*Z2 + X*Y1*Z2 + X*Y3*Z1 - X1*Y2*Z3 + X1*Y3*Z2 +

//

//     Y*X3*Z2 - Y*X1*Z2 - Y*X3*Z1 - Y*X2*Z3 + Y*X2*Z1 + Y*X1*Z3 - X3*Y1*Z2 + X2*Y1*Z3 +

//

//     Z*X2*Y3 - Z*X2*Y1 - Z*X1*Y3 - Z*X3*Y2 + Z*X3*Y1 + Z*X1*Y2 - X2*Y3*Z1 + X3*Y2*Z1  = 0;

//

// 一般式

//

//     (Y2*Z3 - Y2*Z1 - Y1*Z3 - Y3*Z2 + Y1*Z2 + Y3*Z1)*X +

//

//     (X3*Z2 - X1*Z2 - X3*Z1 - X2*Z3 + X2*Z1 + X1*Z3)*Y +

//

//     (X2*Y3 - X2*Y1 - X1*Y3 - X3*Y2 + X3*Y1 + X1*Y2)*Z + 
//
//    X1*Y3*Z2 + X2*Y1*Z3 + X3*Y2*Z1 - X1*Y2*Z3 - X3*Y1*Z2 - X2*Y3*Z1 = 0;