求两圆的交点,Solve[(x - a1)*(b2 - b1) == (y - b1)*(a2 - a1) && (x - c1)*(d2 -
d1) == (y - d1)*(c2 - c1), {x, y}],当然这个很恐怖
即使化简也很恐怖{{x -> (1/(
2 ((a1 - a2)^2 + (b1 - b2)^2)))(a1^3 - a1^2 a2 + a2^3 +
a1 (-a2^2 + (b1 - b2)^2 - r1^2 + r2^2) -
Sqrt[-(b1 - b2)^2 (((a1 - a2)^2 + (b1 - b2)^2 - r1^2)^2 -
2 ((a1 - a2)^2 + (b1 - b2)^2 + r1^2) r2^2 + r2^4)] +
a2 ((b1 - b2)^2 + (r1 - r2) (r1 + r2))),
y -> (1/(2 ((a1 - a2)^2 + (b1 - b2)^2) (b1 -
b2)))(a1^2 (b1 - b2) (b1 + b2) +
a2^2 (b1 - b2) (b1 + b2) + (b1 - b2)^2 (b1^2 - b2^2 - r1^2 +
r2^2) - a2 Sqrt[-(b1 -
b2)^2 (((a1 - a2)^2 + (b1 - b2)^2 - r1^2)^2 -
2 ((a1 - a2)^2 + (b1 - b2)^2 + r1^2) r2^2 + r2^4)] +
a1 (2 a2 (-b1^2 + b2^2) +
Sqrt[-(b1 - b2)^2 (((a1 - a2)^2 + (b1 - b2)^2 - r1^2)^2 -
2 ((a1 - a2)^2 + (b1 - b2)^2 + r1^2) r2^2 +
r2^4)]))}, {x -> (1/(
2 ((a1 - a2)^2 + (b1 - b2)^2)))(a1^3 - a1^2 a2 + a2^3 +
a1 (-a2^2 + (b1 - b2)^2 - r1^2 + r2^2) +
Sqrt[-(b1 - b2)^2 (((a1 - a2)^2 + (b1 - b2)^2 - r1^2)^2 -
2 ((a1 - a2)^2 + (b1 - b2)^2 + r1^2) r2^2 + r2^4)] +
a2 ((b1 - b2)^2 + (r1 - r2) (r1 + r2))),
y -> (1/(
2 ((a1 - a2)^2 + (b1 - b2)^2) (b1 -
b2)))(a1^2 (b1 - b2) (b1 + b2) +
a2^2 (b1 - b2) (b1 + b2) + (b1 - b2)^2 (b1^2 - b2^2 - r1^2 +
r2^2) + a2 Sqrt[-(b1 -
b2)^2 (((a1 - a2)^2 + (b1 - b2)^2 - r1^2)^2 -
2 ((a1 - a2)^2 + (b1 - b2)^2 + r1^2) r2^2 + r2^4)] -
a1 (2 a2 (b1 - b2) (b1 + b2) +
Sqrt[-(b1 - b2)^2 (((a1 - a2)^2 + (b1 - b2)^2 - r1^2)^2 -
2 ((a1 - a2)^2 + (b1 - b2)^2 + r1^2) r2^2 + r2^4)]))}},
束手无策,
但是人类智慧是无穷的,把圆圆心相对于{a1,b1}偏移就简单多了
Solve[(x)^2 + (y)^2 == r1^2 && (x - a)^2 + (y - b)^2 == r2^2, {x,
y}] a=a2-a1,b=b2-b1
{{x -> (a^3 + a (b^2 + r1^2 - r2^2) - Sqrt[
b^2 (-(a^2 + b^2 - r1^2)^2 + 2 (a^2 + b^2 + r1^2) r2^2 - r2^4)])/(
2 (a^2 + b^2)),
y -> (b^4 + b^2 (a^2 + r1^2 - r2^2) +
a Sqrt[b^2 (-(a^2 + b^2 - r1^2)^2 + 2 (a^2 + b^2 + r1^2) r2^2 -
r2^4)])/(2 b (a^2 + b^2))}, {x -> (
a^3 + a (b^2 + r1^2 - r2^2) + Sqrt[
b^2 (-(a^2 + b^2 - r1^2)^2 + 2 (a^2 + b^2 + r1^2) r2^2 - r2^4)])/(
2 (a^2 + b^2)),
y -> (b^4 + b^2 (a^2 + r1^2 - r2^2) -
a Sqrt[b^2 (-(a^2 + b^2 - r1^2)^2 + 2 (a^2 + b^2 + r1^2) r2^2 -
r2^4)])/(2 b (a^2 + b^2))}}
化简 c=a^2+b^2
x (a(c+ r1^2 - r2^2)+-b((r1 + r2)^2-c))/(2c)
y (b(c+ r1^2 - r2^2)-+a((r1 + r2)^2-c))/(2c)
x=x+a1
y=y+b1