给定三个点求其可以确定的圆的半径（C语言）

本是古典何须时尚 2023-01-20 08:01 27阅读 0赞

#include<stdio.h>
    #include<math.h>
    
    int main()
    
    {
        
    	double x1 = 5., y1 = 0., x2 = 0., y2 = 5., x3 = 10., y3 = 5.; //定义double型变量
    	double x12, y12, x13, y13, x12q, y12q, x13q, y13q;
    	double a, b, xm, ym, R;
    
    
    	printf("GEBEN SIE DIE KOORDINATEN x1, y1 DES ERSTEN PUNKTES EIN \n");//控制台显示输入坐标
    
    	printf("x1 = ");//控制台显示输入x1
    	scanf_s("%lf", &x1);//读取输入的坐标x1
    
    	printf("y1 = ");
    	scanf_s("%lf", &y1);
    
    	printf("GEBEN SIE DIE KOORDINATEN x2, y2 DES ZWEITEN PUNKTES EIN \n");
    
    	printf("x2 = ");
    	scanf_s("%lf", &x2);
    
    	printf("y2 = ");
    	scanf_s("%lf", &y2);
    
    	printf("GEBEN SIE DIE KOORDINATEN x3, y3 DES DRITTEN PUNKTES EIN \n");
    
    	printf("x3 = ");
    	scanf_s("%lf", &x3);
    
    	printf("y3 = ");
    	scanf_s("%lf", &y3);
    
    	x12 = x1 - x2; //圆心坐标公式很繁琐，避免出错，引入一些变量帮助计算。
    	x13 = x1 - x3;
    
    	y12 = y1 - y2;
    	y13 = y1 - y3;
    
    	x12q = x1 * x1 - x2 * x2;
    	x13q = x1 * x1 - x3 * x3;
    
    	y12q = y1 * y1 - y2 * y2;
    	y13q = y1 * y1 - y3 * y3;
    
    	a = x12 / y12 - x13 / y13;
    	b = (x12q + y12q) / (2. * y12) - (x13q + y13q) / (2. * y13);
    
    	xm = b / a; ym = (-x12q + 2. * x12 * xm - y12q) / (-2. * y12);
    
    	R = (x1 - xm) * (x1 - xm) + (y1 - ym) * (y1 - ym);
    
    	R = sqrt(R);
    
    	printf("Der Mittelpunkt (xm,ym) = (%lf,%lf) \n", xm, ym);//控制台显示xm和ym
    	printf("Der Radius R = %lf \n", R);//控制台显示半径
    }

本文章为转载的，保存便于自己学习，原链接已丢失