本文介绍了如何使用scipy.optimize.leastsq进行球面、圆柱面和平面的最小二乘拟合，强调了该方法在处理柱面拟合时对初值的敏感性，并提供了误差函数的定义。对于圆柱面拟合，建议根据已知信息限制初值以提高拟合效果。
scipy.optimize.leastsq官方文档
之前自己用最小二乘进行过球面、柱面和平面拟合，球面和平面效果都很好，但是柱面效果一直不好（对初值十分敏感，迭代过程中经常出现奇异矩阵，也有可能是自己写的算法没有考虑数值误差）。所以采用scipy自带的最小二乘拟合进行曲面拟合，当对初值进行一定限制时，拟合效果很好。
球面拟合
球面方程：
( x − x 0 ) 2 + ( y − y 0 ) 2 + ( z − z 0 ) 2 = r 0 2 (x-x_0)^2+(y-y_0)^2+(z-z_0)^2=r_0^2 (x−x0)2+(y−y0)2+(z−z0)2=r02
误差函数：
F ( x 0 , y 0 , z 0 , r 0 ) = ( x − x 0 ) 2 + ( y − y 0 ) 2 + ( z − z 0 ) 2 − r 0 2 F(x_0,y_0,z_0,r_0)=(x-x_0)^2+(y-y_0)^2+(z-z_0)^2-r_0^2 F(x0,y0,z0,r0)=(x−x
最低0.47元/天 解锁文章
扫一扫
1万+
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
