//判断三角形顶点的方向。返回1表示顺时针,0表示三点共线,-1表示逆时针。
int ccw(double start_x,double start_y,double end_x,double end_y,double center_x,double center_y) {
float ax = start_x, ay = start_y;
float bx = end_x, by = end_y;
float cx = center_x, cy = center_y;
float v = (bx - ax) * (cy - ay) - (cx - ax) * (by - ay);
return v < 0 ? -1 : v > 0 ? 1 : 0;
}
//x1,y11是终点坐标,x2,y2是终点坐标,R是半径
void Circle_Center(double x1,double y11, double x2,double y2,double R)
{
double c1 = 0, c2 = 0, A = 0, B = 0, C = 0, y01 = 0, x01 = 0, x02 = 0, y02 = 0;
c1 = (pow(x2, 2) - pow(x1, 2) + pow(y2, 2) - pow(y11, 2)) / 2 /(x2 - x1);
c2 = (y2 - y11) / (x2 - x1);
A = 1.0 + pow(c2, 2);
B = 2 * (x1 - c1) * c2 - 2 * y11;
C = pow((x1 - c1), 2) + pow(y11, 2) - pow(R, 2);
center_y_1 = (-B + sqrt(B*B - 4 * A * C)) / 2 / A;
center_x_1 = c1 - c2 * center_y_1;
center_y_2 = (-B - sqrt(B*B - 4 * A * C)) / 2 / A;
center_x_2 = c1 - c2 * center_y_2;
cout << "圆心坐标1为: (" << center_x_1 << ", " << center_y_1 << ")" << endl;
cout << "圆心坐标2为: (" << center_x_2 << ", " << center_y_2 << ")" << endl;
}
参考:
1.//https://blog.csdn.net/jet_lee01/article/details/40403895?utm_source=copy