问题已解决

桑子

理论是知道的

假设点 p1(x1,y1) p2(x2,y2) p3(x3,y3)

计算 p1, p2 的中垂线方程 f1,

计算 p2, p3 的中垂线方程 f2,

计算 f1 ,f2 方程的交点,就是圆心

圆心和任意一点的距离就是半径

不过我要做批量计算

怎么着能输入三点坐标,直接圆半径就出来了?

[ 本帖最后由 桑子 于 2008-6-10 14:32 编辑 ]

桑子

平时都是用excel这种简单的东西，不知道该怎么对付这个问题了

马克

QQ,
MSN也行

schaefer

还是可以用Exel啊，编辑公式，然后批量计算啊

桑子

原帖由 schaefer 于 2008-6-10 14:00 发表
还是可以用Exel啊，编辑公式，然后批量计算啊

唉，耐心研究公式，一个个输入去了:noway:

桑子

设任意三点坐标（a1，b1），（a2，b2），（a３，b３）
其中任意二点坐标（a1，b1），（a2，b2）
连接
取其中点
不妨设a1>＝a2，b1>=b2
则中点m（a1-a2，b1-b2)
y=kx+b
代入（a1，b1），（a2，b2）
得解析式y1=kx+b
将直线旋转９０度
得解析式y2=-1/k*x + c
代入中点m（a1-a2，b1-b2),得解析式y2
再连接另二点
同理,得解析式y3
y2=y3
kx+b=-1/k*x+c
得圆心o（x，y2）
又因为半径到圆周等长
所以
运用两点距离公式得
r=根号(x-a1)^2 + (y2-b1)^2

公式:
x0=((y3-Y1)*(Y2*Y2-Y1*Y1+X2*X2-X1*X1)+(Y2-Y1)*(Y1* Y1 - y3*y3+X1*X1-x3*x3))/(2*(X2-X1)*(y3-Y1)-2*(x3-X1)*(Y2-Y1))

y0 = ((x3 - X1) * (X2 * X2 - X1 * X1 + Y2 * Y2 - Y1 * Y1) + (X2 - X1) * (X1 * X1 - x3 * x3 + Y1 * Y1 - y3 * y3)) / (2 * (Y2 - Y1) * (x3 - X1) - 2 * (y3 - Y1) * (X2 - X1))

r=sqr((x1-x0)*(x1-x0)+(y1-y0)*(y1-y0))

公式搜索出来了，哪个数学好的帮我看看，x0,y0是不是圆心的计算公式，是的话我就在excel里慢慢输入去了:blackman:

桑子

赞一下google，以后偶一定先认真搜索再上来提问

帖上来，做个备份好了

#include<iostream.h>
#include<math.h>

int main()
{
   int x1,y1,x3,y3;
   double a,b,c,d,e,f;
   double r,k1,k2,x,y,x2,y2;
   cout<<"请输入x1,y1,x2,y2,x3,y3"<<endl;
   cin>>x1>>y1>>x2>>y2>>x3>>y3;
   if((y1==y2)&&(y2==y3))
{
   cout<<"三点不构成圆!"<<endl;
   return 0;
}
   if((y1!=y2)&&(y2!=y3))
{
   k1=(x2-x1)/(y2-y1);
   k2=(x3-x2)/(y3-y2);
}
   if(k1==k2)
{
   cout<<"三点不构成圆!"<<endl;
   return 0;
}
   a=2*(x2-x1);
   b=2*(y2-y1);
   c=x2*x2+y2*y2-x1*x1-y1*y1;
   d=2*(x3-x2);
   e=2*(y3-y2);
   f=x3*x3+y3*y3-x2*x2-y2*y2;
   x=(b*f-e*c)/(b*d-e*a);
   y=(d*c-a*f)/(b*d-e*a);
   cout<<"圆心为("<<x<<","<<y<<")"<<endl;
   r=sqrt((x-x1)*(x-x1)+(y-y1)*(y-y1));
   cout<<"半径为"<<r<<endl;
   return 0;
}

驴子

则中点m（a1-a2，b1-b2)


中点不是（（a1+a2）/2, (b1+b2)/2 ）么？

桑子

你说a和b的值么？

还没有仔细看，打算两个对比一下，看看哪个对

驴子

其中任意二点坐标（a1，b1），（a2，b2）
连接
取其中点
不妨设a1>＝a2，b1>=b2
则中点m（a1-a2，b1-b2)

两点中点坐标是这两点坐标的平均值，不是差值