也就是说,我想尽量减少其中C是球的中心,r是它的半径,每个P是我的n个点集合中的一个点。这些变量显然是Cx、Cy、Cz和r.,在我的例子中,我可以事先获得一个已知的 r ,只留下 C 的组件作为变量。我真的不想使用任何形式的迭代极小化(例如牛顿法、Levenberg-Marquardt法等)--我更喜欢一组线性方程组或一个显式使用SVD的解。
对于线性矩阵方程的求解,可以使用numpy.linalg.solve来实现LAPACK例程。DGESV computes the solution to a real system of linear equations where A is an N-by-N matrix and X and B are N-by-NRHS matrices.
used to factor A as
A = P * L