求解两圆相交的交点坐标

今天遇到了一个问题,要求两个相交圆的交点坐标,用了三种不同方法求解,其中纯代数方法求解比较麻烦,这里没有贴出来。暂时只推导了公式,明天再写程序,程序写好后会更新到博客,如果大家有其它的方法,欢迎一起探讨。因为博客不好编辑公式,只能截图。

这里写图片描述
这里写图片描述
这里写图片描述
这里写图片描述
这里写图片描述
这里写图片描述

本程序是两圆相交交点,输入信息为两圆坐标和半径值。如有更好的算法,欢迎交流!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
在Python中求解两圆交点可以使用numpy库中的intersect函数或者自己编写计算交点的函数。下面是一个示例代码: ```python import numpy as np def circle_intersection(x0, y0, r0, x1, y1, r1): d = np.sqrt((x1 - x0)**2 + (y1 - y0)**2) if d > r0 + r1: return None # 两圆相交 if d < abs(r0 - r1): return None # 一个包含另一个 if d == 0 and r0 == r1: return None # 两圆重合 a = (r0**2 - r1**2 + d**2) / (2 * d) h = np.sqrt(r0**2 - a**2) x2 = x0 + a * (x1 - x0) / d y2 = y0 + a * (y1 - y0) / d x3_1 = x2 + h * (y1 - y0) / d y3_1 = y2 - h * (x1 - x0) / d x3_2 = x2 - h * (y1 - y0) / d y3_2 = y2在Python中求解两圆交点可以使用numpy库中的intersect函数或者自己编写计算交点的函数。下面是一个示例代码: ```python import numpy as np def circle_intersection(x0, y0, r0, x1, y1, r1): d = np.sqrt((x1 - x0)**2 + (y1 - y0)**2) if d > r0 + r1: return None # 两圆相交 if d < abs(r0 - r1): return None # 一个包含另一个 if d == 0 and r0 == r1: return None # 两圆重合 a = (r0**2 - r1**2 + d**2) / (2 * d) h = np.sqrt(r0**2 - a**2) x2 = x0 + a * (x1 - x0) / d y2 = y0 + a * (y1 - y0) / d x3_1 = x2 + h * (y1 - y0) / d y3_1 = y2 - h * (x1 - x0) / d x3_2 = x2 - h * (y1 - y0) / d y3_2 = y2 + h * (x1 - x0) / d return (x3_1, y3_1), (x3_2, y3_2) ``` 其中,x0、y0、r0分别表示第一个坐标和半径,x1、y1、r1分别表示第二个坐标和半径。函数返回两个元组,分别表示两个交点坐标。如果两个相交,则返回None。
评论 17
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值