求教各位老师，怎么用LISP写一个求两个圆的交点的代码

cd342978119

最近刚学lisp，想用程序画一个已知三边长的任意三角形，我的思路是先画L1直线，然后在L1直线的起点跟终点分别画半径为L2和L3的圆，然后找到两个圆的交点，然后把3个点连起来，再把圆删掉就完成了，现在就是求两个圆的交点的程序不知道怎么写，希望各位老师能提点我一下，谢谢了

王鹏_pBZlo

不是应该发求助板块?

819534890

1. ;声明：（原创）调用或修改请保留信息行
   
2. ;求两个曲线的交点表. ----无痕.2004.2
   
3. ;调用方法:(x@_int (car (entsel)) (car(entsel)))
   
4. (defun x@_int (e1 e2 / ptlst ptints)
   
5.     (setq obj1 (vlax-ename->vla-object e1)
   
6.         obj2 (vlax-ename->vla-object e2)
   
7.               ptlst (xl-div (vlax-invoke obj1 'IntersectWith obj2 0) 3))
   
8.             ;;此版本未提供2dpolyline正确的自身交点过滤方法，请自行加入.
   
9. )
   
10. ;按项数分割列表.
    
11. (defun xl-div (lst x / lst2)
    
12.     (foreach n lst
    
13.         (if (and    lst2 (/= x (length (car lst2))))      
    
14.         (setq lst2 (cons (append (car lst2) (list n)) (cdr lst2)))
    
15.         (setq lst2 (cons (list n) lst2))
    
16.         )
    
17.     )(reverse lst2)
    
18. )
    
19. 
    
20. 
    
21. 
    
22. 
    
23. (defun x_inters (ename1 ename2)
    
24.   (setq ptlst (vlax-safearray->list   (vlax-variant-value     (vla-intersectwith (vlax-ename->vla-object ename1) (vlax-ename->vla-object ename2) acextendnone)   )
    
25.   )
    
26.   )
    
27.   (if (> (length ptlst) 3)   (progn (setq #ptlst '())   (repeat (/ (length ptlst) 3)     (setq x   (car ptlst)    ptlst  (cdr ptlst)    y   (car ptlst)    ptlst  (cdr ptlst)    z   (car ptlst)    ptlst  (cdr ptlst)    pt0   (list x y z)    #ptlst (cons pt0 #ptlst)     )   )   (setq #ptlst (reverse #ptlst))   )
    
28.   )
    
29.   #ptlst
    
30. )

;声明：（原创）调用或修改请保留信息行<br /> ;求两个曲线的交点表. ----无痕.2004.2<br /> ;调用方法:(x@_int (car (entsel)) (car(entsel))) <br /> (defun x@_int (e1 e2 / ptlst ptints) <br />     (setq obj1 (vlax-ename->vla-object e1) <br />         obj2 (vlax-ename->vla-object e2) <br />               ptlst (xl-div (vlax-invoke obj1 'IntersectWith obj2 0) 3)) <br />             ;;此版本未提供2dpolyline正确的自身交点过滤方法，请自行加入. <br /> ) <br /> ;按项数分割列表. <br /> (defun xl-div (lst x / lst2) <br />     (foreach n lst <br />         (if (and    lst2 (/= x (length (car lst2))))       <br />         (setq lst2 (cons (append (car lst2) (list n)) (cdr lst2))) <br />         (setq lst2 (cons (list n) lst2)) <br />         ) <br />     )(reverse lst2) <br /> ) <br /> <br /> <br /> <br /> <br /> (defun x_inters (ename1 ename2)<br />   (setq ptlst (vlax-safearray->list   (vlax-variant-value     (vla-intersectwith (vlax-ename->vla-object ename1) (vlax-ename->vla-object ename2) acextendnone)   )<br />   )<br />   )<br />   (if (> (length ptlst) 3)   (progn (setq #ptlst '())   (repeat (/ (length ptlst) 3)     (setq x   (car ptlst)    ptlst  (cdr ptlst)    y   (car ptlst)    ptlst  (cdr ptlst)    z   (car ptlst)    ptlst  (cdr ptlst)    pt0   (list x y z)    #ptlst (cons pt0 #ptlst)     )   )   (setq #ptlst (reverse #ptlst))   )<br />   )<br />   #ptlst<br /> )

Humanway

1. (defun c:triangle (/ a aa b bb c cc code pt1 pt2 pt3)
   
2.   (setq        pt1  (getpoint "输入三角形的插入点")
   
3.         code T
   
4.         )
   
5.   (while code
   
6.     (initget 7)
   
7.     (setq a (getreal "输入三角形第1条边的长度"))
   
8.     (initget 7)
   
9.     (setq b (getreal "输入三角形第2条边的长度"))
   
10.     (initget 7)
    
11.     (setq c (getreal "输入三角形第3条边的长度"))
    
12.     (if        (or (<= (+ a b) c) (<= (+ a c) b) (<= (+ c b) a))
    
13.       (alert "输入的3条边长无法满足三角形的几何条件！")
    
14.       (setq code nil)
    
15.       )
    
16.     )
    
17.   (setq        pt2 (polar pt1 0.0 a)
    
18.         aa  (* a a)
    
19.         bb  (* b b)
    
20.         cc  (* c c)
    
21.         )
    
22.   (cond
    
23.     ((= (+ aa bb) cc)
    
24.      (setq pt3 (polar pt2 1.5708 b))
    
25.      )
    
26.     ((> (+ aa bb) cc)
    
27.      (setq code        (/ (- (+ aa cc) cc) (* 2 a))
    
28.            pt3        (polar pt1 0.0 code)
    
29.            code        (expt (- cc (* code code)) 0.5)
    
30.            pt3        (polar pt3 1.5708 code)
    
31.            )
    
32.      )
    
33.     (T
    
34.      (setq code        (/ (- cc (+ aa bb)) (* 2 a))
    
35.            pt3        (polar pt2 0.0 code)
    
36.            code        (expt (- bb (* code code)) 0.5)
    
37.            pt3        (polar pt3 1.5708 code)
    
38.            )
    
39.      )
    
40.     )
    
41.   (setq a (list pt1 pt2 pt3))
    
42.   (entmake (append (list '(0 . "LWPOLYLINE")
    
43.                          '(100 . "AcDbEntity")
    
44.                          '(100 . "AcDbPolyline")
    
45.                          (cons 90 3)
    
46.                          '(70 . 1)
    
47.                          )
    
48.                    (mapcar '(lambda (pt) (cons 10 pt)) a)
    
49.                    )
    
50.            )
    
51.   (princ)
    
52.   )

(defun c:triangle (/ a aa b bb c cc code pt1 pt2 pt3)<br />   (setq        pt1  (getpoint "输入三角形的插入点")<br />         code T<br />         )<br />   (while code<br />     (initget 7)<br />     (setq a (getreal "输入三角形第1条边的长度"))<br />     (initget 7)<br />     (setq b (getreal "输入三角形第2条边的长度"))<br />     (initget 7)<br />     (setq c (getreal "输入三角形第3条边的长度"))<br />     (if        (or (<= (+ a b) c) (<= (+ a c) b) (<= (+ c b) a))<br />       (alert "输入的3条边长无法满足三角形的几何条件！")<br />       (setq code nil)<br />       )<br />     )<br />   (setq        pt2 (polar pt1 0.0 a)<br />         aa  (* a a)<br />         bb  (* b b)<br />         cc  (* c c)<br />         )<br />   (cond<br />     ((= (+ aa bb) cc)<br />      (setq pt3 (polar pt2 1.5708 b))<br />      )<br />     ((> (+ aa bb) cc)<br />      (setq code        (/ (- (+ aa cc) cc) (* 2 a))<br />            pt3        (polar pt1 0.0 code)<br />            code        (expt (- cc (* code code)) 0.5)<br />            pt3        (polar pt3 1.5708 code)<br />            )<br />      )<br />     (T<br />      (setq code        (/ (- cc (+ aa bb)) (* 2 a))<br />            pt3        (polar pt2 0.0 code)<br />            code        (expt (- bb (* code code)) 0.5)<br />            pt3        (polar pt3 1.5708 code)<br />            )<br />      )<br />     )<br />   (setq a (list pt1 pt2 pt3))<br />   (entmake (append (list '(0 . "LWPOLYLINE")<br />                          '(100 . "AcDbEntity")<br />                          '(100 . "AcDbPolyline")<br />                          (cons 90 3)<br />                          '(70 . 1)<br />                          )<br />                    (mapcar '(lambda (pt) (cons 10 pt)) a)<br />                    )<br />            )<br />   (princ)<br />   )

Humanway

Humanway 发表于 2021-10-11 14:52

2. ((> (+ aa bb) cc)
   
3.      (setq code        (/ (- (+ aa cc) cc) (* 2 a))
   
4.            pt3        (polar pt1 0.0 code)
   
5.            code        (expt (- cc (* code code)) 0.5)
   
6.            pt3        (polar pt3 1.5708 code)
   
7.            )
   
8.      )

((> (+ aa bb) cc)<br />      (setq code        (/ (- (+ aa cc) cc) (* 2 a))<br />            pt3        (polar pt1 0.0 code)<br />            code        (expt (- cc (* code code)) 0.5)<br />            pt3        (polar pt3 1.5708 code)<br />            )<br />      )

这一段代码写错了，本来打算直接修改，系统却说有不合法参数，直接无语了，应该是

2. ((> (+ aa bb) cc)
   
3.      (setq code        (/ (- (+ aa cc) bb) (* 2 a))
   
4.            pt3        (polar pt1 0.0 code)
   
5.            code        (expt (- cc (* code code)) 0.5)
   
6.            pt3        (polar pt3 1.5708 code)
   
7.            )
   
8.      )

((> (+ aa bb) cc)<br />      (setq code        (/ (- (+ aa cc) bb) (* 2 a))<br />            pt3        (polar pt1 0.0 code)<br />            code        (expt (- cc (* code code)) 0.5)<br />            pt3        (polar pt3 1.5708 code)<br />            )<br />      )