javascript - 查找由 2 个点及其切线交点定义的圆心

Question

我有笛卡尔坐标 P、P0 和 P1，如何找到 C？

请注意，C 也是以 P 为顶点，P0 和 P1 为高的三 Angular 形的重心。

此代码将在图形渲染器中每秒运行数百次，在一个已经在处理复杂物理的程序中，因此它需要快速。我正在用它来计算二次曲线的控制点。如果可能的话，我想避免三 Angular 函数和计算 Angular ，并坚持使用纯矢量数学。

纯数学答案很好，即使我是用 Javascript 编写的。请记住，答案越简单越好，因为任何矢量运算都需要时间。

Answer 1

方法一:
向量CP0和PP0是垂直的，CP1和PP1是垂直的，所以我们可以写出这些向量对点积的两个方程组，并找到两个未知的CX和CY(中心坐标)

CP0 = (CX - P0.X, CY - P0.Y)
CP1 = (CX - P1.X, CY - P1.Y)
PP0 = (P.X - P0.X, P.Y - P0.Y)
PP1 = (P.X - P1.X, P.Y - P1.Y)

Dot(CP0, PP0 ) = 0
Dot(CP1, PP1 ) = 0

(CX - P0.X) * (P.X - P0.X) + (CY - P0.Y) * (P.Y - P0.Y) = 0
(CX - P1.X) * (P.X - P1.X) + (CY - P1.Y) * (P.Y - P1.Y) = 0

方法二:
P0和P1之间的中点M有坐标

  M.X = (P0.X + P1.X) / 2          ///1
  M.Y = (P0.Y + P1.X) / 2          ///2
  Distance P-M is DM
  DM_Squared = (P.X-M.X)^2 + (P.Y-M.Y)^2     ///3
  L, D distances are from your notation
  L_Squared = (P.X-P0.X)^2 + (P.Y-P0.Y)^2     ///4

Theta 是 Angular P0PC

 Cos(Theta) = Length(P - M) / Length(P - P0) = DM / L 
 D = Length(C - P) = L / Cos(Theta) = L * L / DM
 and finally
 C = P - L^2 * (P0 - M) / DM^2

 Coeff = L_Squared / DM_Squared           ///5
 C.X = P.X - Coeff * (P.X - M.X)         ///6
 C.Y = P.Y - Coeff * (P.Y - M.Y)         ///7

公式 1..7 可以立即使用