从坐标集精确地对中心点进行三角测量

Question

如何准确地从坐标集计算中心点？

例如，从一组观察结果：

point   xcoord  ycoord  accuracy    time
1       x1      y1      2000        t1
2       x2      y2      2000        t2
3       x3      y3      2000        t3
n       xn      yn      2000        tn

从技术上讲，所有这些点都是具有精确半径的圆的中心点，我们知道我们在时间内至少在一个圆内。

我想计算圆的交叉区域的中心点和精度。

这是一个粗略的例子。绿线试图表示准确性

Answer 1

+1非常有趣的问题......这是我在C ++中的方法

用于存储所有需要的全局表和变量：

const int    N=5;       // number of measurement circles
struct _circle { double x,y,r; _circle() { x=0; y=0; r=0.0; } } circle[N],circle0,circle1;

• circle[] - 测量位置
• circle0 - 原始的真实位置（不为您所知）
• circle1 - 计算出更准确的位置

这是我的随机测量值（复制您的测量值）：

int i;
double r,a,x0,x1,y0,y1;
// set real position
circle0.x=50.0;
circle0.y=50.0;
circle0.r=25.0;
// set meassured positions
Randomize();
for (i=0;i<N;i++)
    {
    r=circle0.r*(0.2+(0.3*double(Random(101))*0.01));
    a= 2.0*M_PI*          double(Random(101))*0.01;
    circle[i]=circle0;
    circle[i].x+=r*cos(a);
    circle[i].y+=r*sin(a);
    }

这是计算平均样式位置的方法：

// compute more accurate position (average style)
circle1.x=0.0;
circle1.y=0.0;
circle1.r=circle0.r;
for (i=0;i<N;i++)
    {
    circle1.x+=circle[i].x;
    circle1.y+=circle[i].y;
    }
circle1.x/=N;
circle1.y/=N;
for (i=0;i<N;i++)
    {
    r=circle1.x-circle[i].x; r*=r;
    a=circle1.y-circle[i].y; a*=a;
    r=circle[i].r-sqrt(a+r);
//  if (circle1.r>r) circle1.r=r;   // best accuracy
    if (circle1.r<r) circle1.r=r;   // worse accuracy
    }

• 在最后两个ifs中选择了你想要的精度半径......
• 位置基于所有测量的平均值

这是如何计算几何样式位置：

// compute more accurate position (geometry style)
x0=circle[i].x; x1=x0;
y0=circle[i].y; y1=y0;
for (i=0;i<N;i++)
    {
    a=circle[i].x; if (x0>a) x0=a; if (x1<a) x1=a;
    a=circle[i].y; if (y0>a) y0=a; if (y1<a) y1=a;
    }
circle1.x=0.5*(x0+x1); x1-=x0; x1*=x1;
circle1.y=0.5*(y0+y1); y1-=y0; y1*=y1;
circle1.r=1.0*sqrt(x1+y1);

• 位置是占用区域的中心

以下是我的代码输出的一些预览：