最小覆盖圆

求三角形外接圆圆心。判断点到圆心的距离,即该点是否在圆内或外。

求外接圆圆心

void  circle_center(Point  center,Point pt[3],double &radiu)    //参数是圆心和半径的引用
{
    double  x1,x2,x3,y1,y2,y3;
    double  x  =  0;
    double  y  =  0;

    x1  = pt[0].x;
    x2  = pt[1].x;
    x3  = pt[2].x;
    y1  = pt[0].y;
    y2  = pt[1].y;
    y3  = pt[2].y;

    x=((y2-y1)*(y3*y3-y1*y1+x3*x3-x1*x1)-(y3-y1)*(y2*y2-y1*y1+x2*x2-x1*x1))/(2*(x3-x1)*(y2-y1)-2*((x2-x1)*(y3-y1)));
    y=((x2-x1)*(x3*x3-x1*x1+y3*y3-y1*y1)-(x3-x1)*(x2*x2-x1*x1+y2*y2-y1*y1))/(2*(y3-y1)*(x2-x1)-2*((y2-y1)*(x3-x1)));

    center.x  =  x  ;
    center.y  =  y  ;
    radiu  =  sqrt((pt[0].x  -  x)*(pt[0].x  -  x)  +  (pt[0].y  -  y)*(pt[0].y  -  y));
}

求两点间距离

double getDis(Point a, Point b)
{
    return sqrt((a.x - b.x) * (a.x - b.x) + (a.y - b.y) * (a.y - b.y));
}

判断是否为钝角

bool judgeDunJiao(Point pt[3], Point center, double &r)
{
    double c = getDis(pt[0], pt[1]);
    double b = getDis(pt[0], pt[2]);
    double a = getDis(pt[1], pt[2]);
    if(a * a > b * b + c * c)
    {
        center.x = (pt[1].x + pt[2].x) / 2.0;
        center.y = (pt[1].y + pt[2].y) / 2.0;
        r = a / 2.0;
        return true;
    }
    else if(b * b > a * a + c * c)
    {
        center.x = (pt[0].x + pt[2].x) / 2.0;
        center.y = (pt[0].y + pt[2].y) / 2.0;
        r = a / 2.0;
        return true;
    }
    else if(c * c > a * a + b * b)
    {
        center.x = (pt[1].x + pt[0].x) / 2.0;
        center.y = (pt[1].y + pt[0].y) / 2.0;
        r = a / 2.0;
        return true;
    }
    else
        return false;

}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值