最小二乘拟合空间平面

 题目的意思是求一个平面P，使得 ∑(Distance(pi,P)^2)最小化, pi是给定点集里的点。 Right？
对于a*x+b*y+c*z+m=0,m!=0的情况,写成a*x+b*y+c*z+1=0
(x,y,z)记做p,
(a,b,c)记做n,
则上述平面可记做n·p+1=0
点pi到平面的有向距离di=n·pi+1

那么∑(Distance(pi,P)^2) = ∑(n·pi + 1)^2  := M(n)
其实就是确定适当的n使得M(n)最小化

联列
∂M/∂a = 0
∂M/∂b = 0
∂M/∂c = 0

推到后面似乎就是解一个3元一次线性方程组，用一下克拉默法则就直接得解了

补充，对于m=0的情况，你不妨也列出来去算，算出来一个m=1时候的最优(a,b,c)，一个m=0的最优(a,b,c)然后检查一下哪个产生的∑Distance最小就取哪个。

zz：http://bbs.bccn.net/thread-287643-1-1.html的leeco