嘿，这个问题在几何计算里算是高频需求了，我来给你拆解清楚——从传统的坐标公式推导，到更简洁的向量代数解法，再附上可直接复用的代码。

一、传统几何坐标公式推导 #

首先明确已知条件：圆心坐标 C(Cx, Cy)，圆半径 a，圆外点 P(Px, Py)。咱们先算几个关键的中间值：

• 点P到圆心C的距离：d = sqrt((Px - Cx)² + (Py - Cy)²)，必须满足 d > a，否则点在圆内/圆上，无实切线（圆上的话切点就是自身）。
• 切线长（P到切点的距离）：h = sqrt(d² - a²)，这是勾股定理来的——CP是斜边，半径a和切线h是直角边。

基于这些，两个切点的坐标公式如下：
先计算 dx = Px - Cx，dy = Py - Cy，d_squared = dx² + dy²（避免重复开根号，提升精度）：

t1x = Cx + (a²*dx - a*dy*h) / d_squared
t1y = Cy + (a²*dy + a*dx*h) / d_squared

t2x = Cx + (a²*dx + a*dy*h) / d_squared
t2y = Cy + (a²*dy - a*dx*h) / d_squared

这个公式是通过旋转向量推导来的，本质和向量解法一致，但用坐标形式写出来更直观。

二、更简便的向量代数解法（无需联立切线方程） #

其实刚才的公式已经隐含了向量思想，但直接用向量运算会更简洁，全程不用解二元方程组，步骤清晰：

1. 计算向量 CP = (Px - Cx, Py - Cy)（从圆心指向外点的向量）
2. 计算CP的模长 d，验证 d > a（边界情况后面说）
3. 把CP单位化：u = CP / d（长度为1的向量，方向和CP一致）
4. 生成垂直于u的单位向量 v = (-u.y, u.x)（另一个垂直向量是 -v）
5. 计算比例系数：k = sqrt(d² - a²)/d（其实就是CP与切点连线的夹角余弦值）
6. 两个切点的向量表达式：

   T1 = C + a*(k*u + sqrt(1 - k²)*v)
   T2 = C + a*(k*u - sqrt(1 - k²)*v)
   

这里 sqrt(1 - k²) 等于 a/d，代入后和传统公式完全等价，但向量写法更符合编程思维，代码实现时逻辑更清晰。

三、代码实现（Python示例） #

下面是封装好的函数，包含边界情况处理（点在圆内、圆上），用浮点数精度容错避免计算误差：

import math

def calculate_tangent_points(center_x, center_y, radius, point_x, point_y):
    # 计算圆心到外点的向量差
    dx = point_x - center_x
    dy = point_y - center_y
    distance_squared = dx ** 2 + dy ** 2
    distance = math.sqrt(distance_squared)
    
    # 处理边界情况
    epsilon = 1e-8  # 浮点数精度容错
    if distance < radius - epsilon:
        # 点在圆内，无实切点
        return None
    elif abs(distance - radius) <= epsilon:
        # 点在圆上，切点就是自身
        return ((point_x, point_y), (point_x, point_y))
    
    # 计算切线长的平方和开方
    tangent_length_squared = distance_squared - radius ** 2
    tangent_length = math.sqrt(tangent_length_squared)
    
    # 计算两个切点坐标
    t1_x = center_x + (radius**2 * dx - radius * dy * tangent_length) / distance_squared
    t1_y = center_y + (radius**2 * dy + radius * dx * tangent_length) / distance_squared
    t2_x = center_x + (radius**2 * dx + radius * dy * tangent_length) / distance_squared
    t2_y = center_y + (radius**2 * dy - radius * dx * tangent_length) / distance_squared
    
    return ((round(t1_x, 6), round(t1_y, 6)), (round(t2_x, 6), round(t2_y, 6)))

# 测试用例：圆心(0,0)，半径1，外点(2,0)，预期切点(0.5, √3/2)和(0.5, -√3/2)
if __name__ == "__main__":
    tangents = calculate_tangent_points(0, 0, 1, 2, 0)
    if tangents:
        print(f"第一个切点：{tangents[0]}")
        print(f"第二个切点：{tangents[1]}")
    else:
        print("该点在圆内，不存在实切线切点")

代码里对结果做了四舍五入保留6位小数，方便查看，你可以根据需求调整精度。

内容的提问来源于stack exchange，提问作者cegprakash