UTF8gbsn
本文介绍一个圆的拟合方法. 参照论文为 Coope, I. D. (1993). Circle fitting
by linear and nonlinear least squares. Journal of Optimization Theory
and Applications, 76(2), 381–388. https://doi.org/10.1007/BF00939613
问题描述
首先还是先来描述一下问题,
假如我们有一个采样数据A={a1,a2,⋯ ,aN}A=\{a_1,a_2, \cdots, a_N\}A={a1,a2,⋯,aN}, 其中aia_iai,
是n维空间中点的坐标向量. 我们希望去拟合一个圆心xxx, 和一个半径rrr,
注意xxx是一个n维向量而rrr是一个标量. 注意重申一下我们的已知集合AAA,
未知x,rx, rx,r.
求解过程
首先, 我们需要写出我们的最优化化方程
argminx,r=∑i=1N[∥x−ai∥22−r2]2\arg\min_{x,r}=\sum_{i=1}^{N}[\|x-a_i\|_2^2-r^2]^2argx,rmin=i=1∑N[∥x−ai∥22−r2]2
接下来我们来分析下中括号内部的这个东西怎么来写.
∥x−ai∥22−r2=(x−ai)T(x−ai)−r2=(xT−aiT)(x−ai)−r2=xTx−2aiTx+aiTai−r2\left.
\begin{aligned}
\|x-a_i\|_2^2-r^2&=(x-a_i)^T(x-a_i)-r^2\\
&=(x^T-a_i^T)(x-a_i)-r^2\\
&=x^