加载中…
个人资料
  • 博客等级:
  • 博客积分:
  • 博客访问:
  • 关注人气:
  • 获赠金笔:0支
  • 赠出金笔:0支
  • 荣誉徽章:
正文 字体大小:

求两圆的相交面积

(2011-07-07 22:49:39)
标签:

杂谈

分类: 计算几何
  • 假设半径小的圆为c1,半径大的圆为c2
  • c1的半径r1,圆心坐标(x1,y1)c2的半径r2,圆心坐标(x2,y2)
  • d为两圆圆心连线的长度。
  • 相交面积为S
  • d=sqrt((x1-x2)^2+(y1-y2)^2)
  • 1)如果r1+r2<=d
  • 那么两圆相离,相交面积S=0
  • 2)如果r2-r1>=d
  • 那么半径小的圆内含半径大的圆,那么相交面积为小圆的面积S=pi*r1*r1
  • 3)既非(1)也非(2
  • 在图上画两个相交圆,结合图像看。
  • 那么两圆相交,连接小圆的圆心与两个圆的交点,连接大圆的圆心和两个圆的交点。
  • 可以发现形成的图形被两个圆心的连线平分成2个全等三角形。
  • 由小圆圆心和交点所连两条线(长度为半径)以及在大圆之内的弧所形成的扇形为S1
  • 由大圆圆心和交点所连两条线(长度为半径)以及在小圆之内的弧所形成的扇形为S2
  • 由小圆圆心和交点所连两条线以及由大圆圆心和交点所连两条线所形成的四边形的面积为S3
  • 可见相交面积S=S1+S2-S3
  • 要求出扇形的面积,要知道扇形的圆心角。
  • 小圆包含的扇形的圆心角为2*a1(考虑一个三角形)
  • a1=acos((r1^2+d^2-r2^2)/(2.0*r1*d)) 余弦定理
  • a2=acos((r2^2+d^2-r1^2)/(2.0*r2*d)) 
  • S1=pi*r1*r1*2*a1/(2*pi)=a1*r1*r1
  • 同理
  • S2=a2*r2*r2
  • S3为一个三角形面积的2
  • S3=2*r1*d*sin(a1)/2=r1*d*sin(a1)
  • S=a1*r1*r1+a2*r2*r2-r1*d*sin(a1)

代码:

 

  • #define pi acos(-1.0)
  • struct node{
  •  double x
  •  double y;
  •  double r;
  • } c[maxn];
  • double area(int i,double r1,int jdouble r2){
  •  double d=
  • sqrt((c[i].x-c[j].x)*(c[i].x-c[j].x)+(c[i].y-  c[j].y)*(c[i].y-c[j].y));//圆心距
  •  if(r1>r2){
  •  double temp=r1;
  •  r1=r2;
  •  r2=temp;
  •  }//r1取小
  •  if(r1+r2<=d) 
  • return 0;//相离
  •  else if(r2-r1>=d) 
  • return pi*r1*r1;//内含
  •  else {
  •  
  •  double a1=acos((r1*r1+d*d-r2*r2)/(2.0*r1*d));
  •  double a2=acos((r2*r2+d*d-r1*r1)/(2.0*r2*d));
  •  return (a1*r1*r1+a2*r2*r2-r1*d*sin(a1));
  •  }//相交
  • }

 

0

阅读 收藏 喜欢 打印举报/Report
  

新浪BLOG意见反馈留言板 欢迎批评指正

新浪简介 | About Sina | 广告服务 | 联系我们 | 招聘信息 | 网站律师 | SINA English | 产品答疑

新浪公司 版权所有