两个相离圆上各取一点P、Q，求线段PQ中点M的轨迹。

本文链接：https://blog.csdn.net/liuxiang3/article/details/120375182

中点轨迹多种表现形式

来源：【南通】张海洋(1531861160) 9/17/2021 6:07:47 PM
两圆上两点的中点轨迹

题目：两个相离圆上各取一点P、Q，求线段PQ中点M的轨迹。

1. 构造坐标系

不妨取一个圆的圆心为原点，另一个圆的圆心为 $(a, 0),$ 并设圆1的半径为 $r_1,$ 圆2的半径为 $r_2$ 。
得到

圆1的方程为: $x^2+y^2= r_1^2\qquad\qquad\;\; (1)$

圆2的方程为: $(x-a)^2+y^2= r_2^2\qquad (2)$

2. 解题过程（笛卡尔坐标系）

设圆1上任意一点 $P(x_1,y_1),$ 圆2上任意一点 $Q(x_2,y_2), PQ$ 的中点 $M (x, y),$ 则有：
$\begin{cases} 2x = x_1+x_2 &(3)\\ 2y = y_1+y_2 &(4)\\ x_1^2+y_1^2=r_1^2 &(5)\\ (x_2-a)^2+y_2^2=r_2^2 &(6)\end{cases}$

消元法，先消去 $x_1,y_1, (3),(4)$ 代入 $(5)$ 得到：
$2x-x_2)^2+(2y-y_2)^2=r_1^2$
$(x-\frac{x_2}{2})^2+(y-\frac{y_2}{2})^2 = (\frac{r_1}{2})^2 \qquad (7)$

或者消去 $x_2,y_2, (3),(4)$ 代入 $(5)$ 得到：
$2x-x_1-a)^2+(2y-y_`)^2=r_2^2$
$(x-\frac{x_1+a}{2})^2+(y-\frac{y_1}{2})^2 = (\frac{r_2}{2})^2 \qquad (8)$

显然以上 $(6), (7)$ 组合与 $(5), (8)$ 组合后，应该等价（TODO：如何验证？）。

动态演示参见Geogebra官方网站上的示例：
速度控制轨迹生成演示
速度控制轨迹生成.PNG

3. 解题过程（极坐标系）

不妨设 $P(r_1;v), Q(a+r_2\cos u,r_2\sin u)$

由 $(7)$ 式得到：
$(x-\frac{a+r_2\cos{u}}{2})^2+(y-\frac{r_2\sin{u}}{2})^2 = (\frac{r_1}{2})^2 \qquad (7')$

$\frac{r_{1}^{2}}{4} + \left(- \frac{r_{2} \sin{\left(u \right)}}{2} + y\right)^{2} + \left(- \frac{a}{2} - \frac{r_{2} \cos{\left(u \right)}}{2} + x\right)^{2}=0$

由 $(8)$ 式得到：
$(x-\frac{a+r_1\cos{v}}{2})^2+(y-\frac{r_1\sin{v}}{2})^2 = (\frac{r_2}{2})^2 \qquad (8')$

$\frac{r_{2}^{2}}{4} + \left(- \frac{r_{1} \sin{\left(v \right)}}{2} + y\right)^{2} + \left(- \frac{a}{2} - \frac{r_{1} \cos{\left(v \right)}}{2} + x\right)^{2}=0$

图形演示参见GeoGebra上的案例：
角度为参数的中点轨迹方程
角度参数中点轨迹方程.PNG

4. Python 源码表示

>>> from sympy import *
>>> r_1,r_2,a = symbols("r_1 r_2 a",real=true,positive=true)
>>> a.is_real and a.is_positive
True
>>> u,v = symbols('u v', real = true)
>>> x,y = symbols('x y')

>>> f1=(x-(a+r_2*cos(u))/2)**2+(y-r_2*sin(u)/2)**2-(r_1/2)**2
>>> f2=(x-(a+r_1*cos(v))/2)**2+(y-r_1*sin(v)/2)**2-(r_2/2)**2

>>> print(latex(f1))
- \frac{r_{1}^{2}}{4} + \left(- \frac{r_{2} \sin{\left(u \right)}}{2} + y\right)^{2} + \left(- \frac{a}{2} - \frac{r_{2} \cos{\left(u \right)}}{2} + x\right)^{2}
>>> print(latex(f2))
- \frac{r_{2}^{2}}{4} + \left(- \frac{r_{1} \sin{\left(v \right)}}{2} + y\right)^{2} + \left(- \frac{a}{2} - \frac{r_{1} \cos{\left(v \right)}}{2} + x\right)^{2}