已知两圆的圆心半径，求交点坐标（几何代数两种解法）——CAD VBA 解决

已于 2024-09-11 11:54:07 修改
阅读量374 收藏
点赞数 16
分类专栏： CAD VBA 文章标签： java 前端 javascript
于 2024-09-10 18:05:56 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/yongshiqq/article/details/142105409
版权

本文用几何和代数两种算法，分别求出两圆相交点坐标。

一、几何法

如下图， dwg图中若干图形，运行代码后提示选择两个圆，然后判断两个圆位置关系和相交点坐标:

本例难点在于通过几何知识求出交点坐标。

几何背景

假设有两个圆：
- 圆1：圆心 ( O_1(x_1, y_1) )，半径 ( r_1 )
- 圆2：圆心 ( O_2(x_2, y_2) )，半径 ( r_2 )

圆心 ( O_1 ) 和 ( O_2 ) 之间的距离为 d ，交点位于两圆的公共弦上。我们的目标是通过代数推导找到公共弦与两圆圆心的几何关系，并证明 a 的代数式。

几何分析

两个圆的交点（如果有两个）在公共弦上，且公共弦的中垂线经过两个圆心 O_1  和O_2 的连线。我们定义  P_0 为公共弦的中点，且它在两个圆心连线 O_1O_2 上。定义：
 a 是圆心 O_1到点 P_0  的距离。
h 是P_0 到交点的垂直距离。

因此，我们可以将 a 定义为从 O_1到公共弦 即  P_0 的距离。

利用余弦定理推导 a 

利用两圆的交点与圆心的几何关系，首先计算 \( a \) 的代数表达式。

1. 定义圆心距离 d：

   d = sqrt{(x_2 - x_1)^2 + (y_2 - y_1)^2

   
2. **两圆相交**：假设两个圆有两个交点，公共弦 \( AB \) 将连线 \( O_1O_2 \) 分成两部分：从 \( O_1 \) 到公共弦的距离 \( a \)，和从 \( O_2 \) 到公共弦的另一段距离。

3. **两圆的关系**：根据几何原理，有：

   r_1^2 = a^2 + h^2

   r_2^2 = (d - a)^2 + h^2 

   其中，h 是从 P_0 到交点的垂直距离。

4. 消去  h^2：从公式 (1) 和 (2) 可以消去 \( h^2 \)，得到：

   r_1^2 - a^2 = r_2^2 - (d - a)^2

   
5. **展开并整理**：

   r_1^2 - a^2 = r_2^2 - (d^2 - 2ad + a^2)

   r_1^2 - a^2 = r_2^2 - d^2 + 2ad - a^2

   r_1^2 - r_2^2 + d^2 = 2ad

   
6. 解出  a ：

         a = (r1 ^ 2 - r2 ^ 2 + d ^ 2) / (2 * d)
         h = Sqr(r1 ^ 2 - a ^ 2)
       

通过这个公式，我们可以进一步计算出交点的坐标，根据三角函数，详见代码。

附部分计代码如下：

#If VBA7 Then
  ' 64位系统声明
  Declare PtrSafe Sub Sleep Lib "kernel32" (ByVal dwMilliseconds As Long)
#Else
  ' 32位系统声明
  Declare Sub Sleep Lib "kernel32" (ByVal dwMilliseconds As Long)
#End If
#If VBA7 Then
    ' 64位系统声明
    Declare PtrSafe Function GetAsyncKeyState Lib "user32" (ByVal vKey As Long) As Integer
#Else
    ' 32位系统声明
    Declare Function GetAsyncKeyState Lib "user32" (ByVal vKey As Long) As Integer
#End If
Sub SelectTwoCircles()
'yngqq@2024年9月10日21:40:46
    Dim escapePressed As Boolean
    escapePressed = False
    Dim ent As AcadEntity
    Dim circle1 As AcadCircle
    Dim circle2 As AcadCircle
    Dim selectionCount As Integer
    Dim center1 As Variant
    Dim center2 As Variant
    Dim x1 As Double, y1 As Double, r1 As Double
    Dim x2 As Double, y2 As Double, r2 As Double
    selectionCount = 0

2000:
    Do While selectionCount < 2
        ' 如果按下ESC键，退出循环
        If GetAsyncKeyState(vbKeyEscape) <> 0 Then
        If GetAsyncKeyState(vbKeyEscape) <> 0 Then
            ThisDrawing.Utility.Prompt "检测到ESC键，退出循环 " & vbCrLf
            MsgBox "已按下Esc键,退出程序", , "CopyRight@yngqq"
                GoTo errocontrol
        End If
        End If
        DoEvents
        ThisDrawing.Utility.Prompt "请选择第" & (selectionCount + 1) & "个圆: "
'        If Err Then
'            Err.Clear
'            GoTo 2000
'        End If
        On Error Resume Next
        ThisDrawing.Utility.GetEntity ent, basePnt, " "
        If Err Then
            Err.Clear
            GoTo 2000
        End If
        ' 判断用户是否选择了一个圆
        If TypeOf ent Is AcadCircle Then
            selectionCount = selectionCount + 1
            If selectionCount = 1 Then
                ' 第一个圆
                Set circle1 = ent
            ElseIf selectionCount = 2 Then
                ' 第二个圆
                Set circle2 = ent
            End If
        Else
            ThisDrawing.Utility.Prompt "选择的不是圆，请重新选择。" & vbCrLf
        End If
        
    Loop
    ' 获取圆心坐标和半径
   ' On Error GoTo 0
    center1 = circle1.Center
    center2 = circle2.Center
    x1 = center1(0): y1 = center1(1): r1 = circle1.Radius
    x2 = center2(0): y2 = center2(1): r2 = circle2.Radius
    Call FindCircleIntersection(x1, y1, r1, x2, y2, r2)
errocontrol:
End Sub

Public Function FindCircleIntersection(x1 As Double, y1 As Double, r1 As Double, x2 As Double, y2 As Double, r2 As Double) As Variant
   
    Dim d As Double
    d = Sqr((x2 - x1) ^ 2 + (y2 - y1) ^ 2)
    
    ' 判断圆的关系
    If d > r1 + r2 Then
        MsgBox "两个圆不相交"
    ElseIf d < Abs(r1 - r2) Then
        MsgBox "一个圆在另一个圆内，且不相交"
    ElseIf d = 0 And r1 = r2 Then
        MsgBox "两个圆重合"
    Else
        ' 圆相交，计算交点
        
        ' 计算 a 和 h
        Dim a As Double, h As Double
        a = (r1 ^ 2 - r2 ^ 2 + d ^ 2) / (2 * d)
        h = Sqr(r1 ^ 2 - a ^ 2)
        
        ' 计算中间点 P0
        Dim P0x As Double, P0y As Double
        P0x = x1 + a * (x2 - x1) / d
        P0y = y1 + a * (y2 - y1) / d
        
        ' 计算两个交点
        Dim x3_1 As Double, y3_1 As Double
        Dim x3_2 As Double, y3_2 As Double
        
        x3_1 = P0x + h * (y2 - y1) / d
        y3_1 = P0y - h * (x2 - x1) / d
        
        x3_2 = P0x - h * (y2 - y1) / d
        y3_2 = P0y + h * (x2 - x1) / d
        
        ' 输出交点坐标
        If d = r1 + r2 Or d = Abs(r1 - r2) Then
            MsgBox "两个圆相切，交点坐标为：" & vbCrLf & "(" & x3_1 & "   ,   " & y3_1 & ")"
        Else
            MsgBox "两个圆相交，交点坐标为：" & vbCrLf & "(" & x3_1 & "   ,   " & y3_1 & ")" & vbCrLf & "和" & vbCrLf & "(" & x3_2 & "   ,   " & y3_2 & ")"
        End If
    End If
End Function

二、代数法：

已知：圆1坐标（a1,b1），半径a，圆1坐标（a2,b2），半径b:
那么两个圆的标准方程式 (即所有符合条件的xy在圆上):
方程1：(x-a1）^2  + (y-b1）^2  = r1^2
方程2：(x-a2）^2  + (y-b2）^2  = r2^2
两个圆交点的xy坐标即满足方程1，又满足方程2,所以交点坐标可以把联立方程1、2求得：
（
    ①：(x-a1）^2  + (y-b1）^2  = r1^2
    ②：(x-a2）^2  + (y-b2）^2  = r2^2
）
①②联系得③: y = AX + B
简化下公式:
③代入①，即一元二次方程，根据求根公式，即可得到相交点坐标x1，x2
y可带入公式 y = A * x + B求解
 y1 = A * x1 + B
 y2 = B * x2 + B
注意：这里圆心坐标a1,a2,程序里圆心坐标是x1,x2，注意转换问题。
程序运行如下：

 附部分代码：

#If VBA7 Then
'64   位系统声明
  Declare PtrSafe Sub Sleep Lib "kernel32" (ByVal dwMilliseconds As Long)
#Else
'32   位系统声明
  Declare Sub Sleep Lib "kernel32" (ByVal dwMilliseconds As Long)
#End If
#If VBA7 Then
'64     位系统声明
    Declare PtrSafe Function GetAsyncKeyState Lib "user32" (ByVal vKey As Long) As Integer
#Else
'32     位系统声明
    Declare Function GetAsyncKeyState Lib "user32" (ByVal vKey As Long) As Integer
#End If

Sub SelectTwoCircles()
'yngqq@2024年9月10日21:40:46
    Dim escapePressed As Boolean
    escapePressed = False
    Dim ent As AcadEntity
    Dim circle1 As AcadCircle
    Dim circle2 As AcadCircle
    Dim selectionCount As Integer
    Dim center1 As Variant
    Dim center2 As Variant
    Dim x1 As Double, y1 As Double, r1 As Double
    Dim x2 As Double, y2 As Double, r2 As Double
    selectionCount = 0

2000:
    Do While selectionCount < 2
         '如果按下ESC键，退出循环
        If GetAsyncKeyState(vbKeyEscape) <> 0 Then
        If GetAsyncKeyState(vbKeyEscape) <> 0 Then
            ThisDrawing.Utility.Prompt "检测到ESC键，退出循环 " & vbCrLf
            MsgBox "已按下Esc键,退出程序", , "CopyRight@yngqq"
              
        End If
        End If
        DoEvents
        ThisDrawing.Utility.Prompt "请选择第" & (selectionCount + 1) & "个圆: "
        On Error Resume Next
        ThisDrawing.Utility.GetEntity ent, basePnt, " "
        If Err Then
            Err.Clear
            GoTo 2000
        End If
        ' 判断用户是否选择了一个圆
        If TypeOf ent Is AcadCircle Then
            selectionCount = selectionCount + 1
            If selectionCount = 1 Then
                ' 第一个圆
                Set circle1 = ent
            ElseIf selectionCount = 2 Then
                 '第二个圆
                Set circle2 = ent
            End If
        Else
            ThisDrawing.Utility.Prompt "选择的不是圆，请重新选择。" & vbCrLf
        End If

    Loop
     '获取圆心坐标和半径
    On Error GoTo 0
    center1 = circle1.Center
    center2 = circle2.Center
    x1 = center1(0): y1 = center1(1): r1 = circle1.Radius
    x2 = center2(0): y2 = center2(1): r2 = circle2.Radius
    Call FindCircleIntersection(x1, y1, r1, x2, y2, r2)
End Sub
Public Function FindCircleIntersection(x1 As Double, y1 As Double, r1 As Double, x2 As Double, y2 As Double, r2 As Double) As Variant

    Dim d As Double
    d = Sqr((x2 - x1) ^ 2 + (y2 - y1) ^ 2)

    ' 判断圆的关系
    If d > r1 + r2 Then
        MsgBox "两个圆不相交"
    ElseIf d < Abs(r1 - r2) Then
        MsgBox "一个圆在另一个圆内，且不相交"
    ElseIf d = 0 And r1 = r2 Then
        MsgBox "两个圆重合"
    Else
        ' 圆相交，计算交点
        ' 计算两个交点
        Dim x3_1 As Double, y3_1 As Double
        Dim x3_2 As Double, y3_2 As Double
   Dim A As Double, B As Double, AA As Double, BB As Double, CC As Double '定义四个中间参数
  A = (x2 - x1) / (y1 - y2)
  B = (r1 ^ 2 - r2 ^ 2 + x2 ^ 2 + y2 ^ 2 - x1 ^ 2 - y1 ^ 2) / (2 * (y2 - y1))
  AA = 1 + A ^ 2
  BB = (2 * A * B - 2 * x1 - 2 * A * y1)
  CC = x1 ^ 2 + B ^ 2 - 2 * y1 * B + y1 ^ 2 - r1 ^ 2
  MsgBox (BB ^ 2 - 4 * AA * CC)
  x3_1 = (-1 * BB + Sqr(BB ^ 2 - 4 * AA * CC)) / (2 * AA)
  x3_2 = (-1 * BB - Sqr(BB ^ 2 - 4 * AA * CC)) / (2 * AA)
  'y可带入公式 y = A * x + B求解
  y3_1 = A * x3_1 + B
  y3_2 = A * x3_2 + B

        ' 输出交点坐标
        If d = r1 + r2 Or d = Abs(r1 - r2) Then
            MsgBox "两个圆相切，交点坐标为：" & vbCrLf & "(" & x3_1 & "   ,   " & y3_1 & ")"
        Else
            MsgBox "两个圆相交，交点坐标为：" & vbCrLf & "(" & x3_1 & "   ,   " & y3_1 & ")" & vbCrLf & "和" & vbCrLf & "(" & x3_2 & "   ,   " & y3_2 & ")"
        End If
    End If
End Function
'
'已知：圆1坐标（a1,b1），半径a，圆1坐标（a2,b2），半径b:
'那么两个圆的标准方程式 (即所有符合条件的xy在圆上):
'方程1：(x-a1）^2  + (y-b1）^2  = r1^2
'方程2：(x-a2）^2  + (y-b2）^2  = r2^2
'两个圆交点的xy坐标即满足方程1，又满足方程2,所以交点坐标可以把联立方程1、2求得：
'（
'    ①：(x-a1）^2  + (y-b1）^2  = r1^2
'    ②：(x-a2）^2  + (y-b2）^2  = r2^2
'）
'①②联系得③: y = AX + B
'简化下公式:
'③代入①，即一元二次方程，根据求根公式，即可得到相交点坐标x1，x2
'y可带入公式 y = a * x + B求解
' y1 = A * x1 + B
' y2 = B * x2 + B
'注意：这里圆心坐标a1,a2,程序里圆心坐标是x1,x2，注意转换问题。

草稿纸：

CAD二次开发、插件、代码代写，详情见下方↓

yngsqq
关注
  • 16
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
专栏目录
已知圆上任意两点求圆心和半径_vb中，已知圆上两点和半径值，求圆心坐标的程序。...
weixin_28912017的博客
01-17 1723
匿名用户1级2016-09-02 回答问题：已知圆上的两点坐标和半径，求圆心。数学分析：这个题目，涉及到简单的数学问题，但是计算比较繁琐。假设已知圆上的两点坐标分别为N(X1，Y1)和M(X2，Y2)，半径为R，圆心坐标为o(a,b),根据数学知识可得到：(x1-a)^2 + (y1-b)^2 = R^2----(1)式(x2-a)^2 + (y2-b)^2 = R^2----(2)式分别展开上述...
已知圆上任意两点求圆心和半径_圆曲线(已知两点和半径)座标计算
weixin_28722991的博客
12-25 3576
for4850P圆曲线(已知两点和半径)座标计算程序for4850P一、编制原理：因景观工程、房建和市政工程的曲线大多是圆曲线并且设计仅给两端点和半径而无曲线要素的前提下，利用圆曲线两端点和半径，基于圆心点座标数学精确解析圆弧上任意点座标。二、程序清单：Lbl0:CLS:{ABCGR}:A"X1":B"Y1":C"X2":G"Y2":RLbl1:J=0:I=POL(C-A,G-B):J<0=...
已知两圆圆心坐标及半径求两圆交点 (C语言|参数方程求解)
科技浪人
09-10 7597
已知两圆圆心坐标及半径求两圆交点 (C语言|参数方程求解)在一个二维平面上给定两个圆的圆心横纵坐标、半径共6个参数, 求交点. 这个问题无非是解二元二次方程组.普通二元二次方程联立消元求解的困难在于, 中间过程里的系数会变得非常复杂, 从而导致容易出错—因为公式毕竟还是要人来推导,人的出错率比计算机要高得多得多—改用圆的参数方程求解, 可以在显著地减轻这个负担. 现在谈谈这问题的求解过程. 选
已知两圆圆心坐标及半径求两圆交点
What's more
05-27 6073
在一个二维平面上给定两个圆的横纵坐标、半径共6个参数，如果两圆不相交，那么输出“NO INTERSECTION”，如果两圆重合，则输出“THE CIRCLES ARE THE SAME”，否则输出交点，相切的话输出1个交点坐标，相交则输出两个，而且要保证横坐标大的在前，如果横坐标一样，则纵坐标大的在前。前两种情况很判断，也很好输出，求交点，无非是解二元二次方程组。然而，这个高中生就能够熟练解决的问
C#——已知两点及半径，求圆心坐标
Ro969668074的博客
08-20 4300
已知A点和B点及半径，其中B点为所在圆的切点，求该圆心坐标 float radius;//半径 float angleO;//圆心角 Vector3 m_yuanxin;//圆心坐标 /// <summary> /// 计算圆心坐标 /// </summary> /// &...
matlab知三点求圆心坐标,已知圆心坐标及半径求圆上各点坐标
weixin_39876739的博客
03-24 2133
functionresult=myfitcicledemo1();y=[223.955000000000,219.390000000000,215.211000000000,211.554000000000,209.053000000000,207.667000000000,207.669000000000,208.964000000000,211.374000000000,215.7730000...
已知圆上两点坐标和半径求圆心
MH3709的博客
07-24 3385
牛客题目链接 已知两点坐标可求出直线方程，当然也包括斜率了，已知斜率tan a，最近刚学的高数上，可根据三角函数的反函数（C++ atan2(y,x)）求出角度（小心为double类型），用到圆心到两点中点这条边可求出圆心。 /**遇到此类型的题目就是先确定圆心再去判断其他店是否在圆上或者圆内**/ #include <map> #include <queue> #include <string> #include<iostream> #include&
给定两个点的坐标和一个半径，求能组成的圆的圆心坐标
weixin_66431784的博客
12-20 2276
给定两个点的坐标和一个半径，求能组成的圆的圆心。 给定两个点坐标分别为(x1, y1), (x2, y2)和半径r，求出圆心的坐标(a, b)。
Python - 两圆相交求交点坐标
Max's Note
05-16 6952
Python - 两圆相交求交点坐标 Max.Bai 2016-05-16 Python - 两圆相交求交点坐标 三轴机械臂求坐标问题，其实转化为平面问题就是两圆相交求交点问题，交点算出来就可以用反三角函数算出各个机械臂的夹角。 已知圆1， 半径R， 坐标（x, y) 圆2， 半径S， 坐标（a, b) 求两圆交点x3, y3   x
C++ 已知两点坐标和半径求圆心坐标程序
积累点滴，保持自我
10-10 5724
已知圆上的两点坐标和半径，求圆心。 数学分析：这个题目，涉及到简单的数学问题，但是计算比较繁琐。 假设已知圆上的两点坐标分别为N(X1，Y1)和M(X2，Y2)，半径为R，圆心坐标为o(a,b),根据数学知识可得到： (x1-a)^2 + (y1-b)^2 = R^2----(1)式 (x2-a)^2 + (y2-b)^2 = R^2----(2)式 分别展开上述两个式子得到 (x1)^2 -...
已知两点坐标和半径求圆心坐标程序
09-18
这里给出一个已知两点坐标和半径，求元新坐标的程序源码，供大家参考，代码涉及大量数学推导，如有疑问，来信请寄taiyangshen80293@sohu.com.
spring事务详细讲解-深入浅出
小电玩
09-08 443
Spring 事务是本身就是对数据库的事务的支持，没有数据库的事务 Spring 是本身无法实现事务的。Spring只提供了统一事务管理接口，具体的实现还是由各数据库自己实现。在使用 Spring 事务时，可以通过注解或编程方式将需要进行事务管理的方法和代码块标记为事务性操作，当这些操作被执行时，Spring 会负责开启时根据当前环境中设置的隔离级别，调整数据库隔离级别。
Java多线程：深入探索与详细解析
m0_63550220的博客
09-08 1437
作为Java中的基本执行单元，线程是轻量级的进程，由线程ID、程序计数器、Java虚拟机栈、本地方法栈、和线程私有内存等部分组成。每个线程都有独立的执行路径，但共享进程的资源（如内存）。：进程是系统资源分配的基本单位，它包含了一个或多个线程以及这些线程运行所需的资源。在Java中，JVM（Java虚拟机）就是一个进程，而运行在JVM上的多个线程则共享JVM的内存空间。：并发指的是多个任务在同一时间段内交替执行，而并行则指的是多个任务在同一时刻真正同时执行。
Java集合：Stack详解
最新发布
yang_brother的博客
09-10 465
Java 中的Stack类是一个后进先出（LIFO, Last In First Out）的数据结构，它继承自Vector类。尽管Stack是一个可用的类，但它被认为是遗留的，并且在新代码中不推荐使用。通常建议使用Deque或ArrayDeque作为替代。
自定义Stater
m0_63624484的博客
09-07 701
我们平时在导入依赖的时候，大部分导入的都是xxx-spring-boot-stater。那是因为它们都基于spring的规则将写好的框架定义成一个stater，这样方便springboot引用它们写好的功能。那我们为什么也要自定义stater呢？
币安/OK现货合约量化系统APP开发
I592O929783的博客
09-08 516
后端：考虑系统需要处理大量实时数据和高频交易的特点，选择高性能的编程语言和技术框架，如Python的Django或Flask框架，Java的Spring Boot等。用户交互界面：设计直观易用的用户交互界面，包括交易界面、策略配置界面、风险管理界面等，方便用户进行交易操作、策略配置和风险管理。市场数据接入：开发市场数据接入模块，从币安等交易所或数据服务商获取实时的市场数据，包括行情数据、订单簿数据、成交数据等。前端：选择适合开发复杂界面的前端框架，如React、Vue.js等，以提供用户友好的交互界面。
Linux服务器Java启动脚本
Jack魏
09-08 542
本文章介绍了如何在Linux服务器上执行Java并启动jar包，通常我们会使用nohup直接启动，但是还是需要手动停止然后再次启动，那如何更优雅的在服务器上启动jar包呢，让我们一起探讨一下吧。
MyBatis面试
Tony_yitao的博客
09-09 644
MyBatis面试
springboot数据库连接由localhost改成IP以后访问报错500（2024/9/7
不起眼小菜菜的博客
09-07 313
步骤很详细，直接上教程……
已知圆上两个点和半径求圆心坐标
05-28
设圆心坐标为 $(x_0,y_0)$，已知圆上两点 $(x_1,y_1)$ 和 $(x_2,y_2)$，以及半径 $r$，则有以下方程组： $$ \begin{cases} (x_1-x_0)^2+(y_1-y_0)^2=r^2\\ (x_2-x_0)^2+(y_2-y_0)^2=r^2 \end{cases} $$ 将方程组展开，化简可得： $$ \begin{cases} x_1^2-2x_1x_0+x_0^2+y_1^2-2y_1y_0+y_0^2=r^2\\ x_2^2-2x_2x_0+x_0^2+y_2^2-2y_2y_0+y_0^2=r^2 \end{cases} $$ 将两个方程相减，消去 $x_0^2$ 和 $y_0^2$，整理可得： $$ \begin{aligned} x_1^2-2x_1x_0+y_1^2-2y_1y_0 - x_2^2+2x_2x_0-y_2^2+2y_2y_0 &=0\\ (x_1-x_2)x_0+(y_1-y_2)y_0 &=\frac{x_1^2+y_1^2-x_2^2-y_2^2}{2} \end{aligned} $$ 这是一个二元一次方程组，可以使用高斯消元等方法求解圆心坐标 $(x_0,y_0)$。
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值