两个圆交点

标签: 算法
7人阅读 评论(0) 收藏 举报
分类:
问题描述:
给两个相交的圆,第一个圆的圆心为(x1,y1),半径为r1,第二个圆的圆心为(x2,y2),半径为r2,求两个圆的交点。


问题分析:
《训练指南》上求两圆交点的模板用了atan2,acos等库函数,精度损失比较严重。
下面介绍一种精度损失较小的做法:
原文地址


首先回顾一下圆的两种表示方法:


圆的标准方程:(x−x0)2+(y−y0)2=r2
圆的参数方程:{x=x0+r⋅cosθy=y0+r⋅sinθ
将第一个圆的参数方程带入第二个圆的标准方程:
(x1+r1cosθ−x2)2+(y1+r1sinθ−y2)2=r22
展开后得到:
2r1(x1−x2)cosθ+2r1(y1−y2)sinθ=r22−r21−(x1−x2)2−(y1−y2)2
令:
a=2r1(x1−x2)
b=2r1(y1−y2)
c=r22−r21−(x1−x2)2−(y1−y2)2
原式变为:
acosθ+bsinθ=c
令cosθ=x,sinθ=1−x2−−−−−√,关于sinθ的正负后面再判断。
代入方程,得到,ax+b1−x2−−−−−√=c
移项再两边平方,(ax−c)2=b2(1−x2)
整理得,(a2+b2)x2−2acx+c2−b2=0
下面就是解一元二次方程了。


将sinθ和cosθ代回到第一个圆的参数方程,能得到交点的坐标。
如果该点不在第二个圆上,说明sinθ是个负数,还需要对这个交点稍作调整。
还有一种特殊情况就是,如果已经确定有两个不同的交点,但解出来的cosθ值只有一个。
说明对应的sinθ值必然一正一负。


#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <vector>


typedef long double LD;


const LD eps = 1e-10;


int dcmp(LD x) {
    if(fabs(x) < eps) return 0;
    return x < 0 ? -1 : 1;
}


LD sqr(LD x) { return x * x; }


struct Point
{
    LD x, y;
    Point(LD x = 0, LD y = 0):x(x), y(y) {}
    void read() { cin >> x >> y; }
};


Point operator - (const Point& A, const Point& B) {
    return Point(A.x - B.x, A.y - B.y);
}


bool operator == (const Point& A, const Point& B) {
    return dcmp(A.x - B.x) == 0 && dcmp(A.y - B.x) == 0;
}


LD Dot(const Point& A, const Point& B) {
    return A.x * B.x + A.y * B.y;
}


LD Length(const Point& A) { return sqrt(Dot(A, A)); }


struct Circle
{
    Point c;
    LD r;
    Circle() {}
    Circle(Point c, LD r):c(c), r(r) {}
};


int getCircleIntersection(Circle C1, Circle C2) {
    LD &r1 = C1.r, &r2 = C2.r;
    LD &x1 = C1.c.x, &x2 = C2.c.x, &y1 = C1.c.y, &y2 = C2.c.y;
    LD d = Length(C1.c - C2.c);
    if(dcmp(fabs(r1-r2) - d) > 0) return -1;
    if(dcmp(r1 + r2 - d) < 0) return 0;
    LD d2 = Dot(C1.c - C2.c, C1.c - C2.c);
    LD a = r1*(x1-x2)*2, b = r1*(y1-y2)*2, c = r2*r2-r1*r1-d*d;
    LD p = a*a+b*b, q = -a*c*2, r = c*c-b*b;


    LD cosa, sina, cosb, sinb;
    //One Intersection
    if(dcmp(d - (r1 + r2)) == 0 || dcmp(d - fabs(r1 - r2)) == 0) {
        cosa = -q / p / 2;
        sina = sqrt(1 - sqr(cosa));
        Point p(x1 + C1.r * cosa, y1 + C1.r * sina);
        if(!OnCircle(p, C2)) p.y = y1 - C1.r * sina;
        inter.push_back(p);
        return 1;
    }
    //Two Intersections
    LD delta = sqrt(q * q - p * r * 4);
    cosa = (delta - q) / p / 2;
    cosb = (-delta - q) / p / 2;
    sina = sqrt(1 - sqr(cosa));
    sinb = sqrt(1 - sqr(cosb));
    Point p1(x1 + C1.r * cosa, y1 + C1.r * sina);
    Point p2(x1 + C1.r * cosb, y1 + C1.r * sinb);
    if(!OnCircle(p1, C2)) p1.y = y1 - C1.r * sina;
    if(!OnCircle(p2, C2)) p2.y = y1 - C1.r * sinb;
    if(p1 == p2)  p1.y = y1 - C1.r * sina;
    inter.push_back(p1);
    inter.push_back(p2);
    return 2;
}
查看评论

计算两个圆相交的交点坐标值

  • 2017年12月15日 11:00
  • 3.51MB
  • 下载

求解两圆相交的交点坐标

今天遇到了一个问题,要求两个相交圆的交点坐标,用了三种不同方法求解,其中纯代数方法求解比较麻烦,这里没有贴出来。暂时只推导了公式,明天再写程序,程序写好后会更新到博客,如果大家有其它的方法,欢迎一起探...
  • zx3517288
  • zx3517288
  • 2016-11-24 21:39:44
  • 1432

求两个相交圆的交点的公式

半径为R的圆心为A,坐标为(x,y)  半径为S的圆心为B,坐标为(a,b)  两圆交点为C,D   AB与CD的交点为E,坐标为(X0,Y0)  过C点垂线与过E点水平线交点为F   令L为AB长度...
  • qq_24378389
  • qq_24378389
  • 2016-12-19 11:24:07
  • 3913

java求两个圆相交坐标

最近由于项目需要,根据两个圆函数求出相交的坐标。实现代码如下,另感谢两圆求交点算法实现Java代码,虽然他所贡献的代码中存在问题,但仍有借鉴意义。1.两个圆相交的数学求法在中学数学中我们知道,一个圆可...
  • u013780605
  • u013780605
  • 2016-09-28 10:46:42
  • 1800

计算几何--两圆的位置关系(求交点个数及交点坐标)

两圆相交。假定圆心分别为C1和C2,半径为r1和r2,圆心距为d,根据余弦定理可以算出C1C2到C1P1的角da,根据向量C1C2的极角a,加减da就可以得到C1P1和C1P2的极角。有了极角,就可以...
  • u013021513
  • u013021513
  • 2015-02-28 13:16:24
  • 1555

求过圆心直线与圆的两个交点(JAVA & OC)

转自:http://blog.csdn.net/yang3wei/article/details/7521298
  • gf771115
  • gf771115
  • 2014-05-14 17:22:46
  • 3881

求两圆的交点(Java)

package com.hj.graphics;import java.util.ArrayList;import com.hj.tools.NumberTools;public class Circ...
  • u013056779
  • u013056779
  • 2015-04-02 15:27:32
  • 1182

zoj 1199 Point of Intersection 求两个圆公切线的交点

题目来源: http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=199 分析: 两个圆,不存在交点的情况。 1)在 半径相同的...
  • zn505119020
  • zn505119020
  • 2014-06-26 15:08:48
  • 394

一般方法编程求直线与圆交点

用直线方程与圆方程联立求方程组解即可,使用时直接套用下面函数中公式即可,(x1,y1)(x2,y2)即为所求! //下面这个函数得到交点,如果有2个就返回第一个求出的。 CPoint GetPoi...
  • zheyejs
  • zheyejs
  • 2013-04-10 09:51:03
  • 1712

C#关于圆与圆的交点问题

关于圆与圆的交点问题,因为求交点的算法比较多,最近在做关于三点定位的项目,搜了好多圆与圆交的程序,好多都是错误的,输入坐标输出的结果不对。要么就是考虑的不全面因此决定自己写一个关于圆与圆求交点的程序,...
  • laoyaozaici1023
  • laoyaozaici1023
  • 2016-06-25 23:22:52
  • 953
    个人资料
    等级:
    访问量: 184
    积分: 33
    排名: 196万+
    文章分类
    文章存档
    最新评论