三点确定一个圆的计算方法

三点确定一个圆的计算方法

最近在写的一个软件需要根据三个坐标点来计算一个圆。因此花了点时间推导了相关的公式。这个推导不算太难，放在这里主要是做个备忘。

我们设一个圆的圆心坐标为 $(x_0, y)$，半径为 $r$。那么这个圆的方程可以写为：

$$(x - x_0)^2 + (y - y_0)^2 = r^2$$

在这个圆上随便取三个点，设这三个点的坐标分别是 $(x_1, y_1)$，$(x_2, y_2)$，$(x_3, y_3)$。那么有：

$$\begin{cases} (x_1 - x_0)^2 + (y_1-y_0)^2 = r^2 & (1)\\ (x_2 - x_0)^2 + (y_2-y_0)^2 = r^2 & (2)\\ (x_3 - x_0)^2 + (y_3-y_0)^2 = r^2 & (3) \end{cases}$$

公式（1）（2）相减，（1）（3）相减之后经过化简可以得到：

$$(x_1 - x_2) x_0 + (y_1- y_2) y_0 = \frac{(x_1^2-x_2^2)-(y_2^2-y_1^2)}{2}\\ (x_1 - x_3) x_0 + (y_1- y_3) y_0 = \frac{(x_1^2-x_3^2)-(y_3^2-y_1^2)}{2}\\ $$

$x_0,y_0$有唯一解的条件是系数行列式不为 $0$：

$$\begin{vmatrix} (x_1 - x_2) & (y_1- y_2) \\ (x_1 - x_3) & (y_1- y_3) \end{vmatrix} \neq 0$$

简单变变型也就是：

$$\frac{x_1 - x_2}{y_1-y_2} \neq \frac{x_1 - x_3}{y_1-y_3}$$

这样写几何含义就很明显了，三点不能共线。

设:

$$a = x_1-x_2\\ b = y_1-y_2\\ c = x_1-x_3\\ d = y_1-y_3\\ e = \frac{(x_1^2-x_2^2)-(y_2^2-y_1^2)}{2}\\ f = \frac{(x_1^2-x_3^2)-(y_3^2-y_1^2)}{2}$$

那么 :

$$x_0 = -\frac{d e-b f}{b c-a d}\\ y_0 = -\frac{a f-c e}{b c-a d}$$

有了 $x_0$ 和 $y_0$ 的值后，带入(1) 式就可以得到 r <script type="math/tex" id="MathJax-Element-17">r</script>的值。至此，三点确定圆的问题就解决了。

下面是个 C++ 代码（用到了Qt 的 QPointF 类型）：

    #include <math.h>
    #include <limits>
    #include <QPoint>
    #include <QDebug>
    QPointF tcircle(QPointF pt1, QPointF pt2, QPointF pt3, double &radius)
    {
        double x1 = pt1.x(), x2 = pt2.x(), x3 = pt3.x();
        double y1 = pt1.y(), y2 = pt2.y(), y3 = pt3.y();
        double a = x1 - x2;
        double b = y1 - y2;
        double c = x1 - x3;
        double d = y1 - y3;
        double e = ((x1 * x1 - x2 * x2) + (y1 * y1 - y2 * y2)) / 2.0;
        double f = ((x1 * x1 - x3 * x3) + (y1 * y1 - y3 * y3)) / 2.0;
        double det = b * c - a * d;
        if( fabs(det) < 1e-5)
        {
            radius = -1;
            return QPointF(0,0);
        }

        double x0 = -(d * e - b * f) / det;
        double y0 = -(a * f - c * e) / det;
        radius = hypot(x1 - x0, y1 - y0);
        return QPointF(x0, y0);
    }

用 {1, 1}, {1, 2}, {3, 4} 这三个点测试了一下。结果没有问题。