數(shù)學求解原理請看:求兩圓相交的交點的方法
實現(xiàn)代碼:
package main
import (
"math"
)
//代表一個點鲤孵,包含橫縱坐標
type Point struct {
X, Y float64
}
//代表一個圓,包含橫縱坐標及半徑
type Circle struct {
Point
R float64
}
//創(chuàng)建圓對象
func NewCircle(x, y, r float64) *Circle {
return &Circle{Point{x, y}, r}
}
//求兩圓相交的交點包颁,交點個數(shù)可能有0,1,2
func Intersect(a *Circle, b *Circle) (p []Point) {
dx, dy := b.X - a.X, b.Y - a.Y
lr := a.R + b.R //半徑和
dr := math.Abs(a.R - b.R) //半徑差
ab := math.Sqrt(dx * dx + dy * dy) //圓心距
if ab <= lr && ab > dr {
theta1 := math.Atan(dy / dx)
ef := lr - ab
ao := a.R - ef / 2
theta2 := math.Acos(ao / a.R)
theta := theta1 + theta2
xc := a.X + a.R * math.Cos(theta)
yc := a.Y + a.R * math.Sin(theta)
p = append(p, Point{xc, yc})
if ab < lr { //兩個交點
theta3 := math.Acos(ao / a.R)
theta = theta3 - theta1
xd := a.X + a.R * math.Cos(theta)
yd := a.Y - a.R * math.Sin(theta)
p = append(p, Point{xd, yd})
}
}
return
}