一个迭代的问题

[zihen520]

设实际圆心坐标为（A，B）。为了求某一截面圆的半径，我们可按下列步骤进行迭代计算。
首先，假定该截面圆的圆心坐标初始值为：A0=0,B0=0;


代码如下：
void CFitCircleDlg::Iteration(double &x,double &y, double &R)
{
double A0=0.0,B0=0.0,R0=0.0,R1,A1,B1,Ri,r;
double X=0.0,Y=0.0,X1=0.0,Y1=0.0;
int count,i,j,k;
Point node,point;
count=m_PointDataArray.GetSize();
R1=0.0;
for (i=0;i<count;i++)
{
node=m_PointDataArray.GetAt(i);
R1+=sqrt((node.East-A0)*(node.East-A0)+(node.North-B0)*(node.North-B0));
}
r=R0;
while(fabs(R1-r)>0.0001)
{
R1=0.0;
for (k=0;k<count;k++)
{
node=m_PointDataArray.GetAt(k);
R1+=sqrt((node.East-A0)*(node.East-A0)+(node.North-B0)*(node.North-B0));
}

X=0.0;
X1=0.0;
Y=0.0;
Y1=0.0;
for(j=0;j<count;j++)
{
point=m_PointDataArray.GetAt(j);
Ri=sqrt((point.East-A0)*(point.East-A0)+(point.North-B0)*(point.North-B0));
X+=point.East;
X1+=(point.East-A0)/Ri;
Y+=point.North;
Y1+=(point.North-B0)/Ri;
}
A1=X/count-R1/count*X1;
B1=Y/count-R1/count*Y1;

if (fabs(R1-R0)>0.0001)
{
r=R0;
R0=R1;
A0=A1;
B0=B1;
}
else
break;
}
R=R1;
x=A1;
y=B1;
}
不知道在哪里弄错了，往各位指点。

[zihen520]

我看得到啊，怎么回事

[zihen520]

本来图片上是一些相关的公式，弄错了，一会再贴出来

[zihen520]

再贴

[zihen520]

谢谢各位了！

[ouyh12345]

此图片来自qq空间
未经允许不可引用

[arong1234]

你代码也太乱了，难道不知道CSDN有贴代码这种方法，可以保持排版的?
估计你是根据一堆采样点估计圆在哪（包括原心和半径），不过看不出你递归的规则

一个建议，愿心的开始迭代点应该是所有采样点的均值，这样至少能加速迭代

[dch4890164]

找两条切线，切线中垂线的焦点既是圆心
代码看不懂