纯数学方法的两个圆相关的函数

marting

1. ;;;圆和圆的交点
   
2. (defun c_int_c (ce1 r1 ce2 r2 / dis a b)
   
3.   ;;by w_kai
   
4.   (setq dis (distance ce1 ce2))
   
5.   (cond ((equal dis (+ r1 r2))
   
6.          (list (polar ce1 (angle ce1 ce2) r1))
   
7.         )
   
8.         ((equal dis (abs (- r1 r2)))
   
9.          (if (minusp (- r1 r2))
   
10.            (list (polar ce2 (angle ce2 ce1) r2))
    
11.            (list (polar ce1 (angle ce1 ce2) r1))
    
12.          )
    
13.         )
    
14.         ((and (> dis (abs (- r1 r2))) (< dis (+ r1 r2)))
    
15.          (setq a (/ (- (+ (* dis dis) (* r1 r1)) (* r2 r2)) (* 2. dis)))
    
16.          (setq b (sqrt (- (* r1 r1) (* a a))))
    
17.          (list (polar (polar ce1 (angle ce1 ce2) a)
    
18.                       (+ (angle ce1 ce2) (/ pi 2))
    
19.                       b
    
20.                )
    
21.                (polar (polar ce1 (angle ce1 ce2) a)
    
22.                       (- (angle ce1 ce2) (/ pi 2))
    
23.                       b
    
24.                )
    
25.          )
    
26.         )
    
27.   )
    
28. )
    

;;;圆和圆的交点 <br /> (defun c_int_c (ce1 r1 ce2 r2 / dis a b)<br />   ;;by w_kai <br />   (setq dis (distance ce1 ce2))<br />   (cond ((equal dis (+ r1 r2))<br />          (list (polar ce1 (angle ce1 ce2) r1))<br />         )<br />         ((equal dis (abs (- r1 r2)))<br />          (if (minusp (- r1 r2))<br />            (list (polar ce2 (angle ce2 ce1) r2))<br />            (list (polar ce1 (angle ce1 ce2) r1))<br />          )<br />         )<br />         ((and (> dis (abs (- r1 r2))) (< dis (+ r1 r2)))<br />          (setq a (/ (- (+ (* dis dis) (* r1 r1)) (* r2 r2)) (* 2. dis)))<br />          (setq b (sqrt (- (* r1 r1) (* a a))))<br />          (list (polar (polar ce1 (angle ce1 ce2) a)<br />                       (+ (angle ce1 ce2) (/ pi 2))<br />                       b<br />                )<br />                (polar (polar ce1 (angle ce1 ce2) a)<br />                       (- (angle ce1 ce2) (/ pi 2))<br />                       b<br />                )<br />          )<br />         )<br />   )<br /> )<br />

2. 
   
3. ;;;
   
4. ;;; 三点圆函数
   
5. (defun 3pcircle (p0 p1 p2 / x0 y0 x1 y1 x2 y2 dx1 dy1 dx2 dy2 d 2d c1 c2 ce)
   
6.   ;;by highflybird
   
7.   (setq x0  (car p0)
   
8.         y0  (cadr p0)
   
9.         x1  (car p1)
   
10.         y1  (cadr p1)
    
11.         x2  (car p2)
    
12.         y2  (cadr p2)
    
13.         dx1 (- x1 x0)
    
14.         dy1 (- y1 y0)
    
15.         dx2 (- x2 x0)
    
16.         dy2 (- y2 y0)
    
17.   )
    
18.   (setq d (- (* dx1 dy2) (* dx2 dy1)))
    
19.   (if (/= d 0.0)
    
20.     (progn (setq 2d (+ d d)
    
21.                  c1 (+ (* dx1 (+ x0 x1)) (* dy1 (+ y0 y1)))
    
22.                  c2 (+ (* dx2 (+ x0 x2)) (* dy2 (+ y0 y2)))
    
23.                  ce (list (/ (- (* c1 dy2) (* c2 dy1)) 2d)
    
24.                           (/ (- (* c2 dx1) (* c1 dx2)) 2d)
    
25.                     )
    
26.            )
    
27.            (list ce (distance ce p0) (list p0 p1 p2))
    
28.     )
    
29.   )
    
30. )
    

<br /> ;;;<br /> ;;; 三点圆函数<br /> (defun 3pcircle (p0 p1 p2 / x0 y0 x1 y1 x2 y2 dx1 dy1 dx2 dy2 d 2d c1 c2 ce)<br />   ;;by highflybird<br />   (setq x0  (car p0)<br />         y0  (cadr p0)<br />         x1  (car p1)<br />         y1  (cadr p1)<br />         x2  (car p2)<br />         y2  (cadr p2)<br />         dx1 (- x1 x0)<br />         dy1 (- y1 y0)<br />         dx2 (- x2 x0)<br />         dy2 (- y2 y0)<br />   )<br />   (setq d (- (* dx1 dy2) (* dx2 dy1)))<br />   (if (/= d 0.0)<br />     (progn (setq 2d (+ d d)<br />                  c1 (+ (* dx1 (+ x0 x1)) (* dy1 (+ y0 y1)))<br />                  c2 (+ (* dx2 (+ x0 x2)) (* dy2 (+ y0 y2)))<br />                  ce (list (/ (- (* c1 dy2) (* c2 dy1)) 2d)<br />                           (/ (- (* c2 dx1) (* c1 dx2)) 2d)<br />                     )<br />            )<br />            (list ce (distance ce p0) (list p0 p1 p2))<br />     )<br />   )<br /> )<br />

newer

使用ARX的几何库，就简洁的多

1、两圆交点

1. (defun c_int_c (cen1 r1 cen2 r2 / g1 g2 mint ints)
   
2.   (and (setq g1 (xdge::constructor "kcircarc3d" cen1 '(0 0 1) r1)) ;构建圆1几何实体
   
3.        (setq g2 (xdge::constructor "kcircarc3d" cen2 '(0 0 1) r2)) ;构建圆2几何实体
   
4.        (setq mint (xdge::constructor "kcurvecurveint3d" g1 g2)) ;构建求交几何实体
   
5.        (setq ints (xdge::getpropertyvalue mint "intpoints"))
   
6.   )
   
7.   (xdrx_object_release g1 g2 mint)
   
8.   ints
   
9. )
   
10. 
    

(defun c_int_c (cen1 r1 cen2 r2 / g1 g2 mint ints)<br />   (and (setq g1 (xdge::constructor "kcircarc3d" cen1 '(0 0 1) r1)) ;构建圆1几何实体<br />        (setq g2 (xdge::constructor "kcircarc3d" cen2 '(0 0 1) r2)) ;构建圆2几何实体<br />        (setq mint (xdge::constructor "kcurvecurveint3d" g1 g2)) ;构建求交几何实体<br />        (setq ints (xdge::getpropertyvalue mint "intpoints"))<br />   )<br />   (xdrx_object_release g1 g2 mint)<br />   ints<br /> )<br /> <br />

2、三点圆

1. (defun 3pcircle (p0 p1 p2 / data)
   
2.   (and (setq g1 (xdge::constructor "kcircarc3d" p0 p1 p2)) ;构建圆1几何实体
   
3.        (setq data (append (xdge::getpropertyvalue g1 "center" "radius") (list (list p0 p1 p2))))
   
4.   )
   
5.   (xdrx_object_release g1)
   
6.   data
   
7. )
   
8. 
   

(defun 3pcircle (p0 p1 p2 / data)<br />   (and (setq g1 (xdge::constructor "kcircarc3d" p0 p1 p2)) ;构建圆1几何实体<br />        (setq data (append (xdge::getpropertyvalue g1 "center" "radius") (list (list p0 p1 p2))))<br />   )<br />   (xdrx_object_release g1)<br />   data<br /> )<br /> <br />