\u5bf9CAD\u4e2d\u5706\u548c\u5f27\u7684\u7814\u7a76\u3002\u5305\u62ec\u5982\u4e0b\u5185\u5bb9\uff1a
\n\u5706\u548c\u5f27\u7684\u57fa\u672c\u8981\u7d20\uff0c\u5706\u7684\u548c\u5f27\u7684\u5468\u957f \uff0c\u9762\u79ef\uff0c\u8d28\u5fc3\uff1b
\n\u5b83\u4eec\u7684\u79cd\u79cd\u4f5c\u56fe\u6cd5\uff0c\u8b6c\u5982\u4e09\u70b9\u753b\u5706\uff0c\u4e09\u5207\u7ebf\u753b\u5706\u7b49\uff1b
\n\u5b83\u4eec\u4e0e\u70b9\u3001\u5706\u548c\u76f4\u7ebf\u7684\u5173\u7cfb\u53ca\u5176\u5224\u65ad\u3002
\n\u7b49\u7b49 \u3002
\n\u4e0b\u9762\u662f\u76f8\u5173\u4ee3\u7801\uff1a
\n
\n[codesyntax lang=”lisp”];;;****************************************************;
\n;;;\u5706\u548c\u5f27\u90e8\u5206 ;
\n;;;****************************************************;<\/p>\n
;;;—————————————————-;
\n;;;\u5224\u65ad\u70b9\u8ddf\u5706\u7684\u4f4d\u7f6e\u5173\u7cfb ;
\n;;;\u8f93\u5165\uff1aPt,\u8981\u5224\u65ad\u7684\u70b9\uff0cCen,\u5706\u5fc3\uff0cRad\u534a\u5f84 ;
\n;;;\u8f93\u51fa\uff1a0 \u4e3a\u5728\u5706\u4e0a\uff0c-1\u4e3a\u5728\u5706\u5185\uff0c1\u4e3a\u5728\u5706\u5916 ;
\n;;;—————————————————-;
\n(defun CIR:GetPos (Pt Cen Rad \/ d)
\n (setq d (distance cen pt))
\n (cond
\n ( (equal d Rad 1e-8) 0)
\n ( (< d Rad) -1)
\n ( t 1)
\n )
\n)<\/p>\n
;;;—————————————————-;
\n;;;\u5224\u65ad\u4e24\u70b9\u5bf9\u67d0\u4e2a\u5706\u662f\u5426\u5f02\u4fa7 ;
\n;;;—————————————————-;
\n(defun CIR:Is_Opposite_Side (rad d1 d2)
\n (or (and (> d1 rad) (< d2 rad)) (and (< d1 rad) (> d2 rad)))
\n)<\/p>\n
;;;—————————————————-;
\n;;;\u5224\u65ad\u4e24\u70b9\u5bf9\u67d0\u4e2a\u5706\u662f\u5426\u540c\u4fa7 ;
\n;;;—————————————————-;
\n(defun CIR:Is_Same_Side (rad d1 d2)
\n (or (and (> d1 rad) (> d2 rad)) (and (< d1 rad) (< d2 rad)))
\n)<\/p>\n
;;;—————————————————-;
\n;;;This is known as the “secant-tangent theorem”, ;
\n;;;”intersecting chords theorem”, “Secant Theorem” ;
\n;;;or the “power-of-a-point theorem” ;
\n;;;\u5272\u5f26\u5b9a\u7406 ;
\n;;;\u8f93\u5165: \u5728\u5706\u4e0a\u7684\u4e00\u70b9Pc\uff0c\u4e00\u7ed9\u5b9a\u70b9P\uff0c\u7ed9\u5b9a\u5706\u5706\u5fc3C\u548c\u534a\u5f84R ;
\n;;;\u8f93\u51fa: \u7ecf\u8fc7\u5b9a\u70b9Pt\u548cPC\u7684\u5f26\u7684\u53e6\u4e00\u70b9 ;
\n;;;—————————————————-;
\n(defun CIR:Intersecting_Chords (Pc P C R \/ d L eps)
\n (setq eps 1e-6)
\n (setq d (distance P C))
\n (if (equal d R eps)
\n P
\n (polar P (angle P Pc) (\/ (* (+ d R) (- d R)) (distance P Pc)))
\n )
\n)<\/p>\n
;;;—————————————————-;
\n;;;\u6c42\u76f4\u7ebf\u548c\u5706\u7684\u4ea4\u70b9 ;
\n;;;\u8f93\u5165\uff1a\u5706\u5fc3\uff0c\u534a\u5f84\u548c\u76f4\u7ebf\u7684\u4e24\u4e2a\u7aef\u70b9 ;
\n;;;\u8f93\u51fa\uff1a\u76f4\u7ebf\u4e0e\u5706\u7684\u4ea4\u70b9\u96c6(nil,\u4e00\u4e2a\u70b9\u6216\u8005\u4e24\u4e2a\u70b9) ;
\n;;;—————————————————-;
\n(defun CIR:Inters_Circle_Line (Cen rad PtA PtB \/ p d L A eps)
\n (setq eps 1e-6)
\n (setq p (Line:Perpendicular_Foot Cen PtA PtB)) ;\u5706\u5fc3\u5230\u76f4\u7ebf\u7684\u5782\u8db3
\n (setq d (abs (car p))) ;\u5782\u8ddd
\n (setq p (cadr p)) ;\u5782\u8db3
\n (if (equal d Rad eps) ;\u5982\u679c\u5782\u8ddd\u7b49\u4e8e\u534a\u5f84
\n (list p) ;\u76f8\u5207
\n (if (< d Rad) ;\u5982\u679c\u5782\u8ddd\u5c0f\u4e8e\u534a\u5f84
\n (progn
\n (setq L (sqrt (* (+ Rad d)(- Rad d)))) ;\u534a\u5f26\u957f
\n (setq A (angle PtA PtB)) ;\u76f4\u7ebf\u7684\u4e0eX\u8f74\u7ebf\u7684\u4ea4\u89d2
\n (list (polar p a (- L)) (polar p a L)) ;\u6709\u4e24\u4e2a\u4ea4\u70b9
\n )
\n )
\n )
\n)<\/p>\n
;;;—————————————————-;
\n;;;\u76f4\u7ebf\u548c\u5706\u7684\u4ea4\u70b9\u5417\uff0c\u4e3aCIR:LCR\u76f8\u5207\u5706\u4fee\u6539 ;
\n;;;\u8f93\u5165\uff1a\u5706\u5fc3\uff0c\u534a\u5f84\uff0c\u89d2\u5ea6\u65b9\u5411\uff0c\u4e0e\u5706\u5fc3\u7684\u8ddd\u79bb ;
\n;;;\u8f93\u51fa\uff1a\u76f4\u7ebf\u4e0e\u5706\u7684\u4ea4\u70b9\u96c6(nil,\u4e00\u4e2a\u70b9\u6216\u8005\u4e24\u4e2a\u70b9) ;
\n;;;—————————————————-;
\n(defun CIR:Inters_Circle_Line_1 (cen rad An Dist \/ d p l s)
\n (setq d (abs dist))
\n (if (equal d Rad 1e-6)
\n (list (polar cen (- an (* pi 0.5)) Dist)) ;\u76f8\u5207
\n (if (< d Rad)
\n (setq p (polar cen (- an (* pi 0.5)) Dist) ;\u5782\u8db3
\n L (sqrt (* (+ rad d)(- rad d))) ;\u534a\u5f26\u957f
\n S (list (polar p an (- L))
\n (polar p an L)
\n )
\n )
\n )
\n )
\n)<\/p>\n
;;;—————————————————-;
\n;;;\u5706\u548c\u5706\u7684\u4ea4\u70b9 ;
\n;;;\u8f93\u5165\uff1a\u4e24\u4e2a\u5706\uff08\u7528\u5706\u5fc3\u548c\u534a\u5f84\u8868\u793a\uff09 ;
\n;;;\u8f93\u51fa\uff1a\u4e24\u4e2a\u5706\u7684\u4ea4\u96c6\uff08nil\uff0c\u4e00\u4e2a\u70b9\u6216\u8005\u4e24\u4e2a\u70b9\uff09 ;
\n;;;—————————————————-;
\n(defun CIR:Inters_Circle_Circle (c1 r1 c2 r2 \/ eps m n L d H A B)
\n (setq eps 1e-6)
\n (setq m (+ r2 r1))
\n (setq n (- r2 r1))
\n (setq L (distance c1 c2))
\n (cond
\n ( (equal L 0 eps) nil) ;\u5706\u5fc3\u91cd\u5408\uff0c\u65e0\u89e3\u6216\u65e0\u7a77\u89e3
\n ( (equal m l eps) ;\u4e24\u4e2a\u5706\u5916\u5207
\n (list (polar c1 (angle c1 c2) r1))
\n )
\n ( (equal (abs n) l eps) ;\u5185\u5207
\n (if (< n 0)
\n (list (polar C2 (angle c2 c1) (- r2)))
\n (list (polar C2 (angle c2 c1) r2))
\n )
\n )
\n ( (and (> m L) (< (abs n) L)) ;\u76f8\u4ea4
\n (setq d (+ (* r1 r1) (* (+ L r2) (- L r2)))) ;\u6839\u636e\u4f59\u5f26\u5b9a\u7406
\n (setq d (\/ d (+ L L)))
\n (setq H (sqrt (* (+ r1 d) (- r1 d)))) ;\u4ea4\u70b9\u4e0e\u5706\u5fc3\u8fde\u7ebf\u5782\u76f4\u8ddd\u79bb
\n (setq A (angle c1 c2)) ;\u5706\u5fc3\u8fde\u7ebf\u89d2
\n (setq B (atan h d))
\n (list (polar c1 (- A B) r1) ;\u6709\u4e24\u4e2a\u4ea4\u70b9
\n (polar c1 (+ A B) r1)
\n )
\n )
\n )
\n)<\/p>\n
;;;—————————————————-;
\n;;;\u5706\u5916\u7684\u4e00\u70b9\u5bf9\u5706\u5f62\u7684\u5207\u70b9 ;
\n;;;\u8f93\u5165: \u7ed9\u5b9a\u5706\u7684\u5706\u5fc3\u548c\u534a\u5f84\u548c\u5706\u5916\u7684\u4e00\u70b9Pt ;
\n;;;\u8f93\u51fa: \u8fd9\u70b9Pt\u5bf9\u5706\u7684\u5207\u70b9\u96c6(nil\uff0c\u4e00\u4e2a\u6216\u8005\u4e24\u4e2a) ;
\n;;;—————————————————-;
\n(defun CIR:Point_Tangent (cen rad pt \/ d l a b p q r eps)
\n (setq eps 1e-8)
\n (setq d (distance cen pt))
\n (if (equal d rad eps)
\n (list pt)
\n (if (> d rad)
\n (setq l (sqrt (* (+ d rad) (- d rad)))
\n a (atan l rad)
\n b (angle cen pt)
\n p (polar cen (+ b a) rad)
\n q (polar cen (- b a) rad)
\n r (list p q)
\n )
\n )
\n )
\n)<\/p>\n
;;;—————————————————-;
\n;;;\u4e24\u4e2a\u5706\u7684\u516c\u5207\u7ebf ;
\n;;;\u8f93\u5165: \u7ed9\u5b9a\u4e24\u4e2a\u5706\u7684\u5706\u5fc3\u548c\u534a\u5f84 ;
\n;;;\u8f93\u51fa: \u8fd9\u4e24\u4e2a\u5706\u7684\u516c\u5207\u7ebf(nil\uff0c\u4e00\u4e2a\u70b9\u6216\u8005\u70b9\u5bf9\u96c6\u5408) ;
\n;;;—————————————————-;
\n(defun CIR:Common_Tangent (c1 r1 c2 r2 \/ a b d eps p q m n u v ret L)
\n (setq d (distance c1 c2))
\n (setq a (angle c1 c2))
\n (setq eps 1e-8)
\n (cond
\n ( (equal d 0 eps) nil) ;\u5706\u5fc3\u91cd\u5408\uff0c\u65e0\u89e3\u6216\u8005\u65e0\u7a77\u89e3
\n ( (equal d (abs (setq b (- r1 r2))) eps) ;\u5185\u5207
\n (setq p (polar c1 a (Math:Sign_reversal r1 b))) ;\u5207\u70b9
\n (list (list p (polar p (+ a (* pi 0.5)) 1000))) ;\u5185\u5207\u7ebf
\n )
\n ( (> d b)
\n (foreach r (list b (+ r1 r2)) ;\u8003\u8651\u5185\u5916\u516c\u5207\u4e24\u79cd\u60c5\u51b5
\n (setq l (* (+ d r) (- d r)))
\n (cond
\n ( (equal l 0 eps) ;\u4e24\u4e2a\u5706\u5916\u5207
\n (setq p (polar c1 a r1))
\n (setq q (polar p (+ a (* pi 0.5)) 1000))
\n (setq ret (cons (list p q) ret))
\n )
\n ( (> l 0)
\n (setq u (atan (sqrt l) r))
\n (setq v (- a u))
\n (setq u (+ a u))
\n (setq p (polar C1 u r1))
\n (setq q (polar C1 v r1))
\n (if (eq r b)
\n (setq m (polar c2 u r2) ;\u5916\u516c\u5207\u7ebf
\n n (polar c2 v r2)
\n )
\n (setq m (polar c2 u (- r2)) ;\u5185\u516c\u5207\u7ebf
\n n (polar c2 v (- r2))
\n )
\n )
\n (setq ret (cons (list p m) ret))
\n (setq ret (cons (list q n) ret))
\n )
\n )
\n )
\n Ret
\n )
\n )
\n)<\/p>\n
;;;—————————————————-;
\n;;;\u4f5c\u7ed9\u5b9a\u534a\u5f84\u5e76\u7ecf\u8fc7\u4e24\u70b9\u7684\u5706 ;
\n;;;\u8f93\u5165: \u4e0d\u91cd\u5408\u7684\u4e24\u70b9\u548c\u534a\u5f84 ;
\n;;;\u8f93\u51fa: \u65e0\u89e3\u6216\u65e0\u7a77\u89e3\u4e3anil,\u5426\u5219\u7528\u4e2d\u5fc3\u534a\u5f84\u8868\u793a\u7684\u5706\u7684\u5217\u8868;
\n;;;—————————————————-;
\n(defun CIR:PPR (P1 P2 R \/ D M H A eps)
\n (setq eps 1e-6)
\n (setq M (GEO:Midpoint P1 P2))
\n (setq D (distance P1 M))
\n (cond
\n ( (equal D 0 eps) nil) ;\u4e24\u70b9\u91cd\u5408\uff0c\u65e0\u7a77\u89e3
\n ( (equal D R eps) ;\u4e24\u70b9\u4e3a\u76f4\u5f84\u7684\u4e24\u4e2a\u7aef\u70b9
\n (list (list M R))
\n )
\n ( (> R D) ;\u6ee1\u8db3\u4e24\u70b9\u8ddd\u79bb\u5c0f\u4e8e\u7ed9\u5b9a\u76f4\u5f84
\n (setq H (sqrt (* (+ R D) (- R D))))
\n (setq A (+ (angle P1 P2) (* pi 0.5)))
\n (list (list (polar M A H) R) ;\u6709\u4e24\u89e3
\n (list (polar M A (- H)) R)
\n )
\n )
\n )
\n)<\/p>\n
;;;—————————————————-;
\n;;;\u7ecf\u8fc7\u4e00\u70b9\u548c\u4e00\u76f4\u7ebf\u76f8\u5207\u7684\u5177\u6709\u7ed9\u5b9a\u534a\u5f84\u7684\u5706 ;
\n;;;\u8f93\u5165: \u70b9\u548c\u76f4\u7ebf\u7684\u4e24\u4e2a\u7aef\u70b9\u4ee5\u53ca\u7ed9\u5b9a\u7684\u534a\u5f84 ;
\n;;;\u8f93\u51fa: \u65e0\u89e3\u6216\u65e0\u7a77\u89e3\u4e3anil,\u5426\u5219\u7528\u4e2d\u5fc3\u534a\u5f84\u8868\u793a\u7684\u5706\u7684\u5217\u8868;
\n;;;—————————————————-;
\n(defun CIR:PLR (Pt0 PtA PtB Rad \/ H P D A B X Y eps)
\n (setq eps 1e-6)
\n (setq P (Line:Perpendicular_Foot Pt0 PtA PtB))
\n (setq H (abs (car P)))
\n (setq P (cadr P))
\n (setq D (+ Rad Rad))
\n (cond
\n ( (equal H D eps)
\n (list (list (GEO:Midpoint Pt0 P) Rad))
\n )
\n ( (equal H 0 eps)
\n (setq A (+ (angle PtA PtB) (\/ pi 2)))
\n (list
\n (list (polar Pt0 A Rad) Rad)
\n (list (polar Pt0 A (- Rad)) Rad)
\n )
\n )
\n ( (< H D)
\n (setq Y (- H Rad))
\n (setq X (sqrt (* (+ Rad Y) (- Rad Y))))
\n (setq B (atan X Y))
\n (setq A (angle Pt0 P))
\n (list
\n (list (polar Pt0 (+ A B) Rad) Rad)
\n (list (polar Pt0 (- A B) Rad) Rad)
\n )
\n )
\n )
\n)<\/p>\n
;;;—————————————————-;
\n;;;\u548c\u4e24\u76f4\u7ebf\u76f8\u5207\u7684\u5177\u6709\u7ed9\u5b9a\u534a\u5f84\u7684\u5706(\u4e00\u822c\u6765\u8bf4\u662f\u4e24\u76f8\u4ea4\u76f4\u7ebf);
\n;;;\u8f93\u5165: \u4e24\u76f4\u7ebf\u7684\u56db\u4e2a\u7aef\u70b9\u4ee5\u53ca\u7ed9\u5b9a\u7684\u534a\u5f84 ;
\n;;;\u8f93\u51fa: \u65e0\u89e3\u6216\u65e0\u7a77\u89e3\u4e3anil,\u5426\u5219\u7528\u4e2d\u5fc3\u534a\u5f84\u8868\u793a\u7684\u5706\u7684\u5217\u8868;
\n;;;—————————————————-;
\n(defun CIR:LLR (PtA PtB PtC PtD Rad \/ a1 a2 a3 a4 a5 L1 L2 Int)
\n (setq int (inters PtA PtB PtC PtD nil))
\n (if int
\n (progn
\n (setq a1 (angle pta ptb))
\n (setq a2 (angle ptc ptd))
\n (setq a3 (* (+ a1 a2) 0.5))
\n (setq a4 (* (- a1 a2) 0.5))
\n (setq a5 (+ a3 (* pi 0.5)))
\n (setq L1 (\/ Rad (sin a4)))
\n (setq L2 (\/ Rad (cos a4)))
\n (list
\n (list (polar int a3 L1) Rad)
\n (list (polar int a3 (- L1)) Rad)
\n (list (polar int a5 L2) Rad)
\n (list (polar int a5 (- L2)) Rad)
\n )
\n )
\n )
\n)<\/p>\n
;;;—————————————————-;
\n;;;\u548c\u4e00\u76f4\u7ebf\u53ca\u4e00\u5b9a\u5706\u76f8\u5207\u7684\u5177\u6709\u7ed9\u5b9a\u534a\u5f84\u7684\u5706 ;
\n;;;\u8f93\u5165: \u76f4\u7ebf\u7684\u4e24\u4e2a\u7aef\u70b9\uff0c\u5706\u7684\u5706\u5fc3\u548c\u534a\u5f84\u4ee5\u53ca\u7ed9\u5b9a\u7684\u534a\u5f84 ;
\n;;;\u8f93\u51fa: \u65e0\u89e3\u6216\u65e0\u7a77\u89e3\u4e3anil,\u5426\u5219\u7528\u4e2d\u5fc3\u534a\u5f84\u8868\u793a\u7684\u5706\u7684\u5217\u8868;
\n;;;—————————————————-;
\n(defun CIR:LCR (PtA PtB C0 R0 R1 \/ d a ret)
\n (setq d (Line:Perpendicular_Distance C0 PtA PtB))
\n (setq A (angle PtA PtB))
\n (foreach h (list (+ d r1) (- d r1))
\n (foreach r (list (+ r0 r1) (abs (- r0 r1)))
\n (foreach p (CIR:Inters_Circle_Line_1 C0 R A h)
\n (setq ret (cons (list p r1) ret))
\n )
\n )
\n )
\n (reverse ret)
\n)<\/p>\n
;;;—————————————————-;
\n;;;\u4e0e\u4e24\u4e2a\u5706\u76f8\u5207\u7684\u5177\u6709\u7ed9\u5b9a\u534a\u5f84\u7684\u5706 ;
\n;;;\u8f93\u5165: \u4e24\u4e2a\u5706\u7684\u5706\u5fc3\u548c\u534a\u5f84\u4ee5\u53ca\u7ed9\u5b9a\u7684\u534a\u5f84 ;
\n;;;\u8f93\u51fa: \u65e0\u89e3\u6216\u65e0\u7a77\u89e3\u4e3anil,\u5426\u5219\u7528\u4e2d\u5fc3\u534a\u5f84\u8868\u793a\u7684\u5706\u7684\u5217\u8868;
\n;;;—————————————————-;
\n(defun CIR:CCR (C1 R1 C2 R2 R0 \/ ret)
\n (foreach R (list (cons (+ R1 R0) (+ R2 R0))
\n (cons (+ R1 R0) (abs (- R2 R0)))
\n (cons (abs (- R1 R0)) (+ R2 R0))
\n (cons (abs (- R1 R0)) (abs (- R2 R0)))
\n )
\n (foreach p (CIR:Inters_Circle_Circle C1 (car R) C2 (cdr R))
\n (setq ret (cons (list p R0) Ret))
\n )
\n )
\n ret
\n)<\/p>\n
;;;—————————————————-;
\n;;;\u901a\u8fc7\u4e00\u5b9a\u70b9\u548c\u4e00\u5b9a\u5706\u76f8\u5207\u7684\u5177\u6709\u7ed9\u5b9a\u534a\u5f84\u7684\u5706 ;
\n;;;\u8f93\u5165: \u70b9\uff0c\u5706\u7684\u5706\u5fc3\u548c\u534a\u5f84\u4ee5\u53ca\u7ed9\u5b9a\u7684\u534a\u5f84 ;
\n;;;\u8f93\u51fa: \u65e0\u89e3\u6216\u65e0\u7a77\u89e3\u4e3anil,\u5426\u5219\u7528\u4e2d\u5fc3\u534a\u5f84\u8868\u793a\u7684\u5706\u7684\u5217\u8868;
\n;;;—————————————————-;
\n(defun CIR:PCR (P0 C0 R0 R1 \/ ret)
\n (foreach r (list (+ r0 r1) (abs (- r0 r1)))
\n (foreach p (CIR:Inters_Circle_Circle p0 r1 c0 r)
\n (setq ret (cons (list p R1) ret))
\n )
\n )
\n ret
\n)<\/p>\n
;;;—————————————————-;
\n;;;\u6c42\u4e09\u89d2\u5f62\u5916\u5fc3 TRI:CircumCenter\uff0cOr Called ExCenter ;
\n;;;\u5c3d\u7ba1\u8fd9\u6837\u5199\u5f88\u9ebb\u70e6\uff0c\u663e\u5f97\u4ee3\u7801\u5f88\u591a\uff0c\u4f46\u8fd0\u884c\u5374\u5f88\u5feb ;
\n;;;\u8f93\u5165: \u7ed9\u5b9a\u4e0d\u5171\u7ebf\u7684\u4e09\u4e2a\u70b9 ;
\n;;;\u8f93\u51fa: \u4e09\u70b9\u7684\u5916\u63a5\u5706(\u5706\u5fc3\u548c\u534a\u5f84\u8868\u793a),nil\u8868\u793a\u4e09\u70b9\u5171\u7ebf ;
\n;;;—————————————————-;
\n(defun CIR:PPP (P0 P1 P2 \/ X0 Y0 X1 Y1 X2 Y2 DX1 DY1 DX2 DY2 D 2D C1 C2 CE)
\n (setq X0 (car P0)
\n Y0 (cadr P0)
\n X1 (car P1)
\n Y1 (cadr P1)
\n X2 (car P2)
\n Y2 (cadr P2)
\n DX1 (- X1 X0)
\n DY1 (- Y1 Y0)
\n DX2 (- X2 X0)
\n DY2 (- Y2 Y0)
\n )
\n (setq D (- (* DX1 DY2) (* DX2 DY1)))
\n (if (equal D 0 1e-14)
\n nil
\n (progn
\n (setq 2D (+ D D)
\n C1 (+ (* DX1 (+ X0 X1)) (* DY1 (+ Y0 Y1)))
\n C2 (+ (* DX2 (+ X0 X2)) (* DY2 (+ Y0 Y2)))
\n CE (List (\/ (- (* C1 DY2) (* C2 DY1)) 2D)
\n (\/ (- (* C2 DX1) (* C1 DX2)) 2D)
\n )
\n )
\n (list CE (distance CE P0))
\n )
\n )
\n)<\/p>\n
;;;—————————————————-;
\n;;;\u4e09\u89d2\u5f62\u5185\u5fc3 Tri:InCenter ;
\n;;;(aX1+bx2+cx3)\/(a+b+c) (aY2+bY2+CY3)\/(a+b+c) ;
\n;;;\u8f93\u5165: \u7ed9\u5b9a\u4e09\u4e2a\u70b9 ;
\n;;;\u8f93\u51fa: \u4e09\u70b9\u7684\u5185\u5207\u5706(\u5706\u5fc3\u548c\u534a\u5f84\u8868\u793a) ;
\n;;;—————————————————-;
\n(defun Tri:InCenter (pa pb pc \/ a b c L I r)
\n (setq a (distance pb pc))
\n (setq b (distance pc pa))
\n (setq c (distance pa pb))
\n (setq L (+ a b c))
\n (if (\/= L 0.0)
\n (setq I (MAT:SxVs (list pa pb pc) (list (\/ a L) (\/ b L) (\/ c L)))
\n R (list I (abs (Line:Perpendicular_Distance I pa pb)))
\n )
\n (list pa 0)
\n )
\n)<\/p>\n
;;;—————————————————-;
\n;;;\u8ddf\u4e09\u6761\u76f4\u7ebf\u76f8\u5207\u7684\u5706\u5f62 ;
\n;;;\u8f93\u5165: \u7ed9\u5b9a\u4e09\u6761\u76f4\u7ebf ;
\n;;;\u8f93\u51fa: \u4e0e\u8fd9\u4e09\u6761\u76f4\u7ebf\u90fd\u76f8\u5207\u7684\u5706\uff08\u5185\u5207\u5706\u548c\u65c1\u5207\u5706\uff09 ;
\n;;;—————————————————-;
\n(defun CIR:LLL (P1a P1b P2a P2b P3a P3b \/ bs1 bs2 ret int rad)
\n (setq bs1 (LINE:Angular_Bisector P1a p1b p2a p2b)) ;\u89d2\u5e73\u5206\u7ebf1
\n (setq bs2 (LINE:Angular_Bisector P2a p2b p3a p3b)) ;\u89d2\u5e73\u5206\u7ebf2
\n (foreach p bs1
\n (foreach q bs2
\n (if (setq int (inters (car p) (cadr p) (car q) (cadr q) nil)) ;\u5982\u679c\u89d2\u5e73\u5206\u7ebf\u76f8\u4ea4
\n (setq rad (LINE:Perpendicular_Distance int p1a p1b) ;\u6c42\u51fa\u4ea4\u70b9\u5230\u4efb\u610f\u7ebf\u6bb5\u8ddd\u79bb\u5373\u534a\u5f84
\n ret (cons (list int (abs rad)) ret) ;\u52a0\u5165\u5230\u89e3\u96c6\u4e2d
\n )
\n )
\n )
\n )
\n ret
\n)<\/p>\n
;;;—————————————————-;
\n;;;\u70b9\u70b9\u5207\u7ebf\u5706,\u6c42\u901a\u8fc7\u4e24\u5b9a\u70b9\uff0c\u548c\u4e00\u76f4\u7ebf\u76f8\u5207\u7684\u5706 ;
\n;;;\u8f93\u5165: P1,P2,\u901a\u8fc7\u7684\u70b9,PA,PB\u76f4\u7ebf\u4e0a\u7684\u4e24\u70b9\u3002 ;
\n;;;\u8f93\u51fa: \u7b26\u5408\u6761\u4ef6\u7684\u5706\u7684\u96c6\u5408(nil\uff0c\u4e00\u4e2a\u6216\u591a\u4e2a\u5706\u7684\u96c6\u5408) ;
\n;;;—————————————————-;
\n(defun CIR:PPL (P1 P2 PA PB \/ AN D1 D2 P3 DD INT C1 C2 L1 L2 Mid NRM T1 T2 RT eps)
\n (setq eps 1e-6)
\n (if (and (not (equal p1 p2 eps)) (Line:IsSameSide P1 P2 PA PB)) ;\u540c\u4fa7\u4e14\u4e0d\u91cd\u5408
\n (progn
\n (setq int (inters p1 p2 pa pb nil))
\n (setq nrm (MAT:Rot90 (mapcar ‘- pa pb))) ;\u5782\u76f4PAPB\u7684\u77e2\u91cf
\n (setq Mid (GEO:Midpoint p1 p2)) ;P1P2\u7684\u4e2d\u70b9
\n (if int
\n (progn
\n (setq l1 (distance int p1))
\n (setq l2 (distance int p2))
\n (setq P3 (GEO:Rot90 Mid p1 p2))
\n (if (or (equal int p1 eps) (equal int p2 eps)) ;\u5982\u679cP1\u6216P2\u8ddf\u4ea4\u70b9\u91cd\u5408\u5219\u4ea4\u70b9\u4e3a\u5207\u70b9
\n (setq C1 (inters Mid P3 int (mapcar ‘+ int nrm) nil) ;\u4e24\u70b9\u7684\u5782\u76f4\u5e73\u5206\u7ebf\u8ddf\u4ea4\u70b9\u5904\u7684\u5782\u7ebf\u7684\u4ea4\u70b9
\n RT (list (list C1 (distance C1 P1))) ;\u4e3a\u5706\u5fc3\uff0c\u534a\u5f84\u81ea\u7136\u786e\u5b9a\u4e86\u3002
\n )
\n (setq an (angle pA pB)
\n dd (sqrt (* l1 l2)) ;\u5229\u7528\u5272\u5f26\u5b9a\u7406\u5f97\u5230\u4ea4\u70b9\u5230\u5207\u70b9\u8ddd\u79bb
\n t1 (polar int an dd) ;\u5207\u70b91
\n t2 (polar int an (- dd)) ;\u5207\u70b92
\n C1 (inters Mid P3 t1 (mapcar ‘+ t1 nrm) nil) ;\u4e24\u70b9\u7684\u5782\u76f4\u5e73\u5206\u7ebf\u8ddf\u8fc7\u5207\u70b91\u7684\u5782\u7ebf\u7684\u4ea4\u70b9
\n C2 (inters Mid P3 t2 (mapcar ‘+ t2 nrm) nil) ;\u4e24\u70b9\u7684\u5782\u76f4\u5e73\u5206\u7ebf\u8ddf\u8fc7\u5207\u70b92\u7684\u5782\u7ebf\u7684\u4ea4\u70b9
\n RT (list (list C1 (distance C1 P1)) ;\u753b\u57061
\n (list C2 (distance C2 P2)) ;\u753b\u57062
\n )
\n )
\n )
\n )
\n (if (setq t1 (inters Mid (mapcar ‘+ Mid nrm) PA PB nil)) ;\u4e24\u70b9\u8ddf\u76f4\u7ebf\u5e73\u884c\u7684\u60c5\u51b5
\n (list (CIR:PPP p1 p2 t1))
\n )
\n )
\n )
\n )
\n)<\/p>\n
;;;—————————————————-;
\n;;;\u70b9\u7ebf\u7ebf\u5207\u5706\uff0c\u6c42\u901a\u8fc7\u67d0\u5b9a\u70b9\u5e76\u4e0e\u4e24\u6761\u76f4\u7ebf\u76f8\u5207\u7684\u5706 ;
\n;;;\u8f93\u5165: P,\u901a\u8fc7\u7684\u70b9,PA,PB,PC,PD\u4e24\u76f4\u7ebf\u4e0a\u7684\u56db\u4e2a\u7aef\u70b9\u3002 ;
\n;;;\u8f93\u51fa: \u7b26\u5408\u6761\u4ef6\u7684\u5706\u7684\u96c6\u5408(nil\uff0c\u4e00\u4e2a\u6216\u591a\u4e2a\u5706\u7684\u96c6\u5408) ;
\n;;;—————————————————-;
\n(defun CIR:PLL-1 (p pA pB Int \/ H L X r1 r2 an) ;\u7279\u6b8a\u60c5\u51b5\u4e0b\u7684\u5904\u7406P\u70b9\u5728\u89d2\u5e73\u5206\u7ebf\u4e0a
\n (setq H (Line:Perpendicular_Distance p PA PB))
\n (setq H (abs H))
\n (setq L (distance p Int))
\n (setq X (* H L))
\n (setq r1 (\/ X (+ L H)))
\n (setq r2 (\/ X (- L H)))
\n (setq an (angle P int))
\n (list (list (polar P an r1) r1)
\n (list (polar p an (- r2)) r2)
\n )
\n)
\n(defun CIR:PLL (P pA pB pC pD \/ A A1 A2 D D1 D2 EPS H INT L P1 P2 PM R)
\n (setq eps 1e-6)
\n (if (setq int (inters pA pB pC pD nil)) ;\u5224\u65ad\u4e24\u76f4\u7ebf\u662f\u5426\u76f8\u4ea4
\n (if (equal int P eps) ;\u5982\u679c\u4ea4\u70b9\u8ddf\u5b9a\u70b9\u91cd\u5408
\n (list (list int 0)) ;\u89c6\u4f5c\u534a\u5f84\u4e3a0
\n (progn
\n (if (equal int PA eps) (setq p1 PA PA PB PB P1))
\n (if (equal int PC eps) (setq p2 PC PC PD PD P2))
\n (setq D1 (Line:Perpendicular_Distance p int PA));\u5b9a\u70b9\u5230\u76f4\u7ebf1\u8ddd\u79bb
\n (setq D2 (Line:Perpendicular_Distance p int PC));\u5b9a\u70b9\u5230\u76f4\u7ebf2\u8ddd\u79bb
\n (setq a1 (angle int PA))
\n (setq a2 (angle int PC))
\n (setq a1 (* (+ a1 a2) 0.5)) ;\u5185\u89d2\u5e73\u5206\u7ebf
\n (setq a2 (+ a1 (* pi 0.5))) ;\u5916\u89d2\u5e73\u5206\u7ebf
\n (setq p1 (GEO:Mirror2D P int a1)) ;\u5bf9\u5185\u89d2\u5e73\u5206\u7ebf\u7684\u955c\u50cf\u70b9
\n (setq p2 (GEO:Mirror2D P int a2)) ;\u5bf9\u5916\u89d2\u5e73\u5206\u7ebf\u7684\u955c\u50cf\u70b9
\n (if (or (equal p1 P eps) (equal p2 P eps)) ;\u5982\u679c\u8ddf\u89d2\u5e73\u5206\u7ebf\u91cd\u5408
\n (CIR:PLL-1 p pA pB Int) ;\u6309\u7279\u6b8a\u60c5\u51b5\u5904\u7406
\n (if (or (equal d1 0 eps) (equal d2 0 eps)) ;\u5982\u679c\u5b9a\u70b9\u5728\u4e00\u6761\u76f4\u7ebf\u4e0a
\n (append ;\u5219\u5185\u5916\u89d2\u5e73\u5206\u7ebf\u90fd\u8981\u8003\u8651
\n (CIR:PPL P p1 PA PB)
\n (CIR:PPL P p2 PA PB)
\n )
\n (if (MATH:Same_Sign d1 d2) ;\u5269\u4e0b\u7684\u8f6c\u5316\u4e3aCIR:PPL\u95ee\u9898
\n (CIR:PPL P p2 PA PB) ;\u540c\u53f7\u662f\u5916\u89d2\u5e73\u5206\u7ebf\u7684\u955c\u50cf\u70b9
\n (CIR:PPL P p1 PA PB) ;\u5f02\u53f7\u662f\u5185\u89d2\u5e73\u5206\u7ebf\u7684\u955c\u50cf\u70b9
\n )
\n )
\n )
\n )
\n )
\n (progn ;\u4ee5\u4e0b\u662f\u4e3a\u5e73\u884c\u7684\u4e24\u76f4\u7ebf\u8003\u8651
\n (setq D1 (Line:Perpendicular_Foot p PA PB)) ;\u5b9a\u70b9\u5230\u76f4\u7ebf\u7684\u8ddd\u79bb
\n (setq D2 (Line:Perpendicular_Foot p PC PD)) ;\u5b9a\u70b9\u5230\u76f4\u7ebf\u7684\u8ddd\u79bb
\n (setq D (Line:Perpendicular_Distance PA PC PD)) ;\u4e24\u5e73\u884c\u76f4\u7ebf\u7684\u8ddd\u79bb
\n (setq D (abs D))
\n (setq P1 (cadr D1))
\n (setq P2 (cadr D2))
\n (setq D1 (abs (car D1)))
\n (setq D2 (abs (car D2)))
\n (setq pM (GEO:Midpoint P1 P2)) ;\u4e24\u5782\u8db3\u7684\u4e2d\u70b9
\n (cond
\n ( (or (equal D1 0 eps) (equal D2 0 eps)) ;\u5982\u679c\u5b9a\u70b9\u5728\u67d0\u4e00\u76f4\u7ebf\u4e0a
\n (list (list PM (* D 0.5))) ;\u4e24\u5782\u8db3\u7684\u4e2d\u70b9\u662f\u5706\u5fc3
\n )
\n ( (equal D (+ d1 d2) eps) ;\u5b9a\u70b9\u53ea\u6709\u5728\u4e24\u6761\u76f4\u7ebf\u4e4b\u95f4\u624d\u6709\u89e3
\n (setq R (* D 0.5)) ;\u534a\u5f84\u662f\u4e24\u5e73\u884c\u76f4\u7ebf\u8ddd\u79bb\u7684\u4e00\u534a
\n (setq H (- d1 R))
\n (setq L (sqrt (- (* R R) (* H H))))
\n (setq A (angle PA PB))
\n (list (list (polar PM A L) R) ;\u6b64\u65f6\u6709\u4e24\u4e2a\u89e3
\n (list (polar pM A (- L)) R)
\n )
\n )
\n )
\n )
\n )
\n)<\/p>\n
;;;—————————————————-;
\n;;;\u70b9\u70b9\u5207\u7684\u5706(\u6c42\u901a\u8fc7\u4e24\u4e0d\u76f8\u540c\u7684\u70b9\u5e76\u4e0e\u4e00\u4e2a\u5706\u76f8\u5207\u7684\u5706) ;
\n;;;\u8f93\u5165: \u4e24\u5b9a\u70b9\u548c\u7ed9\u5b9a\u7684\u4e00\u4e2a\u5706\u7684\u5706\u5fc3\u53ca\u534a\u5f84 ;
\n;;;\u8f93\u51fa: \u7b26\u5408\u6761\u4ef6\u7684\u5706\u7684\u96c6\u5408(nil\uff0c\u4e00\u4e2a\u6216\u591a\u4e2a\u5706\u7684\u96c6\u5408) ;
\n;;;—————————————————-;
\n(defun CIR:PPC (p1 p2 cen rad \/ PA PB pC PD PM RT an d1 d2 dd eps)
\n (setq eps 1e-6)
\n (setq dd (distance p1 p2))
\n (setq d1 (distance cen p1))
\n (setq d2 (distance cen p2))
\n (setq pM (GEO:Midpoint p1 p2))
\n (setq an (+ (angle p1 p2) (\/ pi 2)))
\n (cond
\n ( (equal dd 0 eps) nil)
\n ( (or (equal d1 rad eps) (equal d2 rad eps))
\n (if (equal d1 rad eps)
\n (setq PA p1 PB p2)
\n (setq PA p2 PB p1)
\n )
\n (setq pC (polar PM an rad))
\n (setq pD (inters pM PC PA CEN nil))
\n (if PD
\n (list (list pD (distance pD p1)))
\n )
\n )
\n ( (equal d1 d2 eps)
\n (setq pA (polar cen an rad))
\n (setq pB (polar cen an (- rad)))
\n (foreach p (list PA PB)
\n (setq rt (cons (CIR:PPP p1 p2 p) rt))
\n )
\n )
\n ( (Cir:Is_Same_Side rad d1 d2)
\n (setq pA (CIR:Radical_Axis PM (\/ dd 2) Cen rad))
\n (setq PB (polar pA (+ (angle pA cen) (\/ pi 2)) rad))
\n (setq pC (inters pA pB p1 p2 nil))
\n (foreach p (CIR:Point_Tangent cen rad pC)
\n (setq rt (cons (CIR:PPP p1 p2 p) rt))
\n )
\n )
\n )
\n)
\n;;;sqrt(((dee\/dee)^2-1)*(dee^2))
\n;;;sqrt((x^2\/a^2-1)*b^2)<\/p>\n
;;;—————————————————-;
\n;;;\u628a\u70b9\u5706\u5706\u5207\u95ee\u9898\u8f6c\u5316\u4e3a\u70b9\u70b9\u5706\u5207\u95ee\u9898 ;
\n;;;Antihomologous Points ;
\n;;;—————————————————-;
\n(defun CIR:PCC->PPC (pt c1 r1 c2 r2 \/ eps k p q d s)
\n (setq eps 1e-8)
\n (foreach a (list (+ r1 r2) (- r1 r2))
\n (if (equal a 0 eps) ;\u534a\u5f84\u76f8\u7b49\u60c5\u51b5\u4e0b\u4f5c\u955c\u50cf
\n (setq p (GEO:Midpoint c1 c2)
\n q (GEO:Mirror2D pt p (+ (angle c1 c2) (\/ pi 2)))
\n s (cons q s)
\n )
\n (progn
\n (setq p (GEO:Scale c2 c1 (\/ r1 a))) ;\u5bf9\u4e8e\u5185\u5916\u5207\u7ebf\u4ea4\u70b9\u7684\u6620\u5c04
\n (setq d (distance pt p))
\n (and (= a (+ r1 r2)) (setq d (- d)))
\n (if (equal d 0 eps)
\n s
\n (setq k (\/ (distance c1 c2) a)
\n k (* r1 r2 (1- (* k k)))
\n q (polar p (angle p pt) (\/ k d))
\n s (cons q s)
\n )
\n )
\n )
\n )
\n )
\n)<\/p>\n
;;;—————————————————-;
\n;;;\u70b9\u5207\u5207\u7684\u5706(\u6c42\u901a\u8fc7\u4e00\u5b9a\u70b9\u5e76\u4e0e\u4e24\u7ed9\u5b9a\u5706\u76f8\u5207\u7684\u5706) ;
\n;;;\u8f93\u5165: \u4e00\u5b9a\u70b9\u548c\u7ed9\u5b9a\u7684\u4e24\u4e2a\u5706\u7684\u5706\u5fc3\u53ca\u534a\u5f84 ;
\n;;;\u8f93\u51fa: \u7b26\u5408\u6761\u4ef6\u7684\u5706\u7684\u96c6\u5408(nil\uff0c\u4e00\u4e2a\u6216\u591a\u4e2a\u5706\u7684\u96c6\u5408) ;
\n;;;—————————————————-;
\n(defun CIR:PCC (pt Cen1 rad1 Cen2 rad2 \/ HC an c EPS H PM v x y)
\n (setq eps 1e-8)
\n (setq HC (CIR:PCC->PPC pt cen1 rad1 cen2 rad2))
\n (setq pM (GEO:Midpoint Cen1 Cen2))
\n (foreach p HC
\n (if (equal p pt eps)
\n (progn
\n (setq H (Line:Perpendicular_Distance pt cen1 cen2))
\n (setq x (distance cen1 pM))
\n (setq x (* x x))
\n (setq a (+ (angle cen1 cen2) (\/ pi 2)))
\n (foreach k (list (+ H rad1) (- H rad1))
\n (if (not (equal k 0 eps))
\n (setq y (\/ (- (* K K) x) K 2)
\n c (polar pm a y)
\n v (cons (list c (distance pt c)) v)
\n )
\n )
\n )
\n )
\n (setq v (append (CIR:PPC p pt cen1 rad1) v))
\n )
\n )
\n v
\n)<\/p>\n
;;;—————————————————-;
\n;;;\u70b9\u7ebf\u5706\u5207\u7684\u5706(\u6c42\u901a\u8fc7\u4e00\u5b9a\u70b9\u5e76\u4e0e\u4e00\u76f4\u7ebf\u548c\u4e00\u5706\u90fd\u76f8\u5207\u7684\u5706);
\n;;;\u8f93\u5165: \u4e00\u5b9a\u70b9\u548c\u76f4\u7ebf\u7684\u4e24\u4e2a\u7aef\u70b9\u7ed9\u5b9a\u7684\u5706\u7684\u5706\u5fc3\u53ca\u534a\u5f84 ;
\n;;;\u8f93\u51fa: \u7b26\u5408\u6761\u4ef6\u7684\u5706\u7684\u96c6\u5408(nil\uff0c\u4e00\u4e2a\u6216\u591a\u4e2a\u5706\u7684\u96c6\u5408) ;
\n;;;—————————————————-;
\n(defun CIR:PLC_1 (P0 Pa Pb R0 C0 \/ eps f a b c d x y y0 y1 dy k k1 a1 a2 s p)
\n (setq eps 1e-8)
\n (setq f (Line:Equation Pa Pb))
\n (setq A (car f))
\n (setq B (cadr f))
\n (setq C (caddr f))
\n (setq D (sqrt (+ (* A A) (* B B))))
\n (setq y0 (\/ (+ (* A (car C0)) (* B (cadr C0)) C) D))
\n (setq y1 (\/ (+ (* A (car P0)) (* B (cadr P0)) C) D))
\n (setq dy (- y0 y1))
\n (setq k1 (- (* dy (+ y0 y1)) (* R0 R0)))
\n (setq a1 (angle Pa Pb))
\n (setq a2 (+ a1 (* pi 0.5)))
\n (setq s nil)
\n (foreach R (list R0 (- R0))
\n (setq k (+ dy R))
\n (setq k (+ k k))
\n (if (not (equal k 0 eps))
\n (progn
\n (setq y (\/ k1 k))
\n (setq p (polar C0 a2 (- y y0)))
\n (setq x (* y1 (- (+ y y) y1)))
\n (cond
\n ( (equal x 0 eps)
\n (setq s (cons (list p (abs y)) s))
\n )
\n ( (> x 0)
\n (setq x (sqrt x))
\n (setq s (cons (list (polar p a1 x) (abs y)) s))
\n (setq s (cons (list (polar p a1 (- x)) (abs y)) s))
\n )
\n )
\n )
\n )
\n )
\n s
\n)<\/p>\n
;;;—————————————————-;
\n;;;\u70b9\u7ebf\u5706\u5207\u7684\u5706(\u6c42\u901a\u8fc7\u4e00\u5b9a\u70b9\u5e76\u4e0e\u4e00\u76f4\u7ebf\u548c\u4e00\u5706\u90fd\u76f8\u5207\u7684\u5706);
\n;;;\u8f93\u5165: \u4e00\u5b9a\u70b9\u548c\u76f4\u7ebf\u7684\u4e24\u4e2a\u7aef\u70b9\u7ed9\u5b9a\u7684\u5706\u7684\u5706\u5fc3\u53ca\u534a\u5f84 ;
\n;;;\u8f93\u51fa: \u7b26\u5408\u6761\u4ef6\u7684\u5706\u7684\u96c6\u5408(nil\uff0c\u4e00\u4e2a\u6216\u591a\u4e2a\u5706\u7684\u96c6\u5408) ;
\n;;;—————————————————-;
\n(defun CIR:PLC (Pt0 PtA PtB Cen Rad \/ v1 v2 P1 P2 PO an q m i c s)
\n (setq v1 (mapcar ‘- cen Pt0))
\n (setq v2 (mapcar ‘- PtA ptb))
\n (if (equal (MAT:Dot v1 v2) 0 1e-6)
\n (CIR:PLC_1 Pt0 PtA PtB Rad Cen) ;\u5b9a\u70b9\u5728\u8fc7\u5706\u5fc3\u7684\u5782\u76f4\u76f4\u7ebf\u7684\u7ebf\u4e0a
\n (progn
\n (setq p1 (CIR:Radical_Axis_Point Cen Rad Pt0)) ;\u7b49\u5e42\u8f74\u7ebf
\n (setq p2 (GEO:Rot90 P1 cen Pt0)) ;\u7b49\u5e42\u8f74\u7ebf\u53e6\u5916\u4e00\u70b9
\n (setq pO (inters P1 P2 PtA PtB nil)) ;\u7b49\u5e42\u4e2d\u5fc3
\n (setq an (+ (angle PtA PtB) (\/ pi 2))) ;\u5782\u76f4\u76f4\u7ebf\u7684\u89d2\u5ea6
\n (setq p1 (polar cen an rad)) ;\u5706\u7684\u76f4\u5f84\u7aef\u70b91
\n (setq p2 (polar cen an (- rad))) ;\u5706\u7684\u76f4\u5f84\u7aef\u70b92
\n (foreach p (list p1 p2) ;\u5bf9\u6bcf\u4e2a\u76f4\u5f84\u7aef\u70b9
\n (setq q (CIR:Polar_Point_1 Cen Rad Pt0 p)) ;\u6c42\u51fa\u5706\u5bf9\u5230\u7aef\u70b9\u548cPt0\u7684\u76f4\u7ebf\u7684\u6781\u70b9
\n (foreach n (CIR:Inters_Circle_Line Cen rad pO q);\u5bf9\u6bcf\u4e2a\u6781\u70b9\u548c\u7b49\u5e42\u4e2d\u5fc3\u7684\u76f4\u7ebf\u4e0e\u5706\u7684\u4ea4\u70b9
\n (setq m (Geo:MidPoint Pt0 n)) ;\u4e2d\u70b9
\n (setq i (Geo:Rot90 m Pt0 n)) ;\u5782\u76f4\u5e73\u5206\u7ebf
\n (if (setq c (inters m i Cen n nil)) ;\u5982\u679c\u76f8\u4ea4\u5219\u4ea4\u70b9\u662f\u5706\u5fc3
\n (setq s (cons (list c (distance c n)) s)) ;\u52a0\u5165\u5230\u89e3\u96c6
\n )
\n )
\n )
\n s
\n )
\n )
\n)<\/p>\n
;;;—————————————————-;
\n;;;\u5207\u7ebf\u5207\u7ebf\u5207\u5706\uff08\u5207\u4e24\u5b9a\u76f4\u7ebf\u548c\u4e00\u5b9a\u5706\u7684\u5706\uff09 ;
\n;;;\u8f93\u5165: \u4e24\u5b9a\u76f4\u7ebf\u7684\u56db\u4e2a\u7aef\u70b9\u7ed9\u5b9a\u7684\u5706\u7684\u5706\u5fc3\u53ca\u534a\u5f84 ;
\n;;;\u8f93\u51fa: \u7b26\u5408\u6761\u4ef6\u7684\u5706\u7684\u96c6\u5408(nil\uff0c\u4e00\u4e2a\u6216\u591a\u4e2a\u5706\u7684\u96c6\u5408) ;
\n;;;—————————————————-;
\n(defun CIR:LLC (PA PB PC PD C0 R0 \/ Ret Pts eps Int a1 a2 a3 a4 Va Vc Vp Vx Vy L G S P1 P2 Pt c d r)
\n (setq eps 1e-6)
\n (defun CIR:Equal (P1 P2)
\n (equal p1 p2 1e-6)
\n )
\n (if (setq int (inters PA PB PC PD nil))
\n (progn
\n (if (equal pa int eps)
\n (mapcar ‘set ‘(pb pa) (list pa pb)) ;\u786e\u4fddPa\u4e0eint\u4e0d\u91cd\u5408
\n )
\n (if (equal pc int eps)
\n (mapcar ‘set ‘(pd pc) (list pc pd)) ;\u786e\u4fddPc\u4e0eint\u4e0d\u91cd\u5408
\n )
\n (setq a1 (angle int pa))
\n (setq a2 (angle int pc))
\n (setq Va (mapcar ‘- Pa Int)) ;\u4ea4\u70b9\u5230Pa\u7684\u77e2\u91cf
\n (setq Vc (mapcar ‘- Pc Int)) ;\u4ea4\u70b9\u5230Pc\u7684\u77e2\u91cf
\n (setq a3 (* (- a1 a2) 0.5)) ;\u89d2\u5ea6\u5dee\u7684\u4e00\u534a
\n (setq a4 (- a1 a3)) ;\u89d2\u5ea6\u548c\u7684\u4e00\u534a
\n (setq P1 (polar int a4 1000)) ;\u89d2\u5e73\u5206\u7ebf1
\n (setq P2 (polar int (+ a4 (* pi 0.5)) 1000)) ;\u89d2\u5e73\u5206\u7ebf2
\n (foreach a (list 0 (* pi 0.5) Pi (* pi -0.5))
\n (setq L (abs (\/ r0 (sin (+ a3 a))))) ;\u5706\u5fc3\u5230\u5706\u7684\u6781\u70b9\u7684\u8ddd\u79bb
\n (setq G (polar C0 (+ a4 a) L))
\n (setq S (CIR:Inters_Circle_Line C0 R0 Int G)) ;\u5982\u679c\u8ddf\u5706\u76f8\u4ea4
\n (if S
\n (progn
\n (setq Vp (mapcar ‘- G int))
\n (setq Vx (car (trans Vp 0 Va)))
\n (setq Vy (car (trans Vp 0 Vc)))
\n (if (or (equal vx 0 eps) (equal vy 0 eps)) ;\u5982\u679c\u8ddf\u5176\u4e2d\u4e00\u8fb9\u76f8\u5207
\n (setq Pt (list p1 p2)) ;\u5219\u5185\u89d2\u548c\u5916\u89d2\u5e73\u5206\u7ebf\u90fd\u8981\u8003\u8651
\n (if (MATH:Opposite_Sign vx vy) ;\u8ddd\u79bb\u65b9\u5411\u662f\u5426\u76f8\u53cd
\n (setq Pt (list p1)) ;\u5982\u679c\u8ddd\u79bb\u65b9\u5411\u76f8\u53cd\u5219\u662f\u5185\u89d2\u5e73\u5206\u7ebf
\n (setq Pt (list p2)) ;\u5982\u679c\u8ddd\u79bb\u65b9\u5411\u76f8\u540c\u5219\u662f\u5916\u89d2\u5e73\u5206\u7ebf
\n )
\n )
\n (foreach p S
\n (if (not (MISC:IsExist p Pts ‘CIR:Equal)) ;\u907f\u514d\u91cd\u590d
\n (progn
\n (setq pts (cons p pts))
\n (foreach q Pt
\n (setq c (inters p C0 Int q nil))
\n (if c
\n (setq ret (cons (list c (distance c p)) ret))
\n (setq ret (append (CIR:PLL-1 p PA PB Int) ret))
\n )
\n )
\n )
\n )
\n )
\n )
\n )
\n )
\n )
\n (progn
\n (setq d (LINE:Perpendicular_Distance PA PC PD))
\n (setq r (* d 0.5))
\n (setq l (LINE:offset PC PD r))
\n (grdraw (car l) (cadr l) 6)
\n (setq r (abs r))
\n (foreach x (list (+ r0 r) (abs (- r0 r)))
\n (foreach c (CIR:Inters_Circle_Line C0 x (car l) (cadr l))
\n (setq ret (cons (list c r) ret))
\n )
\n )
\n )
\n )
\n ret
\n)<\/p>\n
;;;—————————————————-;
\n;;;\u5207\u7ebf\u5706\u5706\u753b\u5706\uff08\u4e0e\u4e00\u76f4\u7ebf\u548c\u4e24\u4e2a\u5706\u90fd\u76f8\u5207\u7684\u5706\uff09 ;
\n;;;\u8f93\u5165: \u7ed9\u5b9a\u4e00\u6761\u76f4\u7ebf\u548c\u4e24\u4e2a\u5706 ;
\n;;;\u8f93\u51fa: \u7b26\u5408\u6761\u4ef6\u7684\u5706\u7684\u96c6\u5408(nil\uff0c\u4e00\u4e2a\u6216\u591a\u4e2a\u5706\u7684\u96c6\u5408) ;
\n;;;—————————————————-;
\n(defun CIR:LCC (Pa Pb C1 R1 C2 R2)
\n (if (LINE:Colinearity C1 Pa Pb)
\n (mapcar ‘set ‘(C1 R1 C2 R2) (list C2 R2 C1 R1))
\n )
\n (setq IC2 (CIR:Circle_Inversion C1 R1 C2 R2))
\n (setq IC3 (CIR:Line_Inversion C1 R1 Pa Pb))
\n (setq lst (append (list C1 R1) IC2 IC3))
\n (setq ret nil)
\n (foreach c (apply ‘CIR:CCC lst)
\n (setq ret (cons (CIR:Circle_Inversion C1 R1 (car c) (cadr c)) Ret))
\n )
\n)<\/p>\n
;;;—————————————————-;
\n;;;\u5207\u5207\u5207\u7684\u5706(\u4e0e\u4e09\u4e2a\u5706\u90fd\u76f8\u5207\u7684\u5706) Apollonius’ problem ;
\n;;;\u8f93\u5165: \u4e09\u4e2a\u7ed9\u5b9a\u7684\u5706\u7684\u5706\u5fc3\u53ca\u534a\u5f84 ;
\n;;;\u8f93\u51fa: \u7b26\u5408\u6761\u4ef6\u7684\u5706\u7684\u96c6\u5408(nil\uff0c\u4e00\u4e2a\u6216\u591a\u4e2a\u5706\u7684\u96c6\u5408) ;
\n;;;—————————————————-;
\n(defun CIR:CCC (c1 r1 c2 r2 c3 r3 \/ CC1 CC2 CC3 CCC CEN DAT DMC EX1 EX2 EX3 HC1 HC2 HC3 HCS EPS NT1 NT2 NT3
\n PC1 PC2 PC3 PS1 PS2 PS3 Pt1 Pt2 PT3 RR1 RR2 RR3 IN1 IN2 IN3 KA1 KB1 KC1
\n LN1 MT1 MT2 MT3 RAD RET ISALL I AN ARG C IDX INC LST)
\n (setq RET nil)
\n (setq eps 1e-6)
\n (cond
\n ( (and (equal c1 c2 eps) (equal c2 c3 eps)) nil)
\n ( (LINE:Colinearity c1 c2 c3)
\n (setq lst (list (list C1 R1) (list C2 R2) (list C3 R3)))
\n (setq idx (vl-sort-i (mapcar ‘cadr lst) ‘>))
\n (setq C1 (nth (car idx) lst))
\n (setq C2 (nth (cadr idx) lst))
\n (setq C3 (nth (caddr idx) lst))
\n (setq R1 (cadr C1))
\n (setq R2 (cadr C2))
\n (setq R3 (cadr C3))
\n (setq C1 (car C1))
\n (setq C2 (car C2))
\n (setq C3 (car C3))
\n (if (equal C1 C2 eps)
\n (setq an (angle c1 c3))
\n (setq an (angle c1 c2))
\n )
\n (setq InC (polar C1 (+ an (* pi 0.5)) (+ R1 R1)))
\n (setq arg nil)
\n (setq ret nil)
\n (foreach n (list (list C1 R1) (list C2 R2) (list C3 R3))
\n (setq C (CIR:Circle_Inversion InC R1 (car n) (cadr n)))
\n (setq arg (append c arg))
\n )
\n (foreach n (apply ‘CIR:CCC arg)
\n (setq C (CIR:Circle_Inversion InC R1 (car n) (cadr n)))
\n (if (numberp (cadr c))
\n (setq ret (cons c ret))
\n )
\n )
\n )
\n (t
\n (setq DMC (CIR:Power_Center c1 r1 c2 r2 c3 r3))
\n (setq HC1 (CIR:Homothetic_Center c2 r2 c3 r3))
\n (setq HC2 (CIR:Homothetic_Center c3 r3 c1 r1))
\n (setq HC3 (CIR:Homothetic_Center c1 r1 c2 r2))
\n (setq in1 (car HC1))
\n (setq in2 (car HC2))
\n (setq in3 (car HC3))
\n (setq ex1 (cadr HC1))
\n (setq ex2 (cadr HC2))
\n (setq ex3 (cadr HC3))<\/p>\n
(setq CCC (car DMC))
\n (setq dat (cadr DMC))
\n (setq Cen (list c1 c2 c3))
\n (setq Rad (list r1 r2 r3))<\/p>\n
(setq HCs (list (list ex1 ex2 ex3 T)
\n (list in3 in1 ex2)
\n (list in1 in2 ex3)
\n (list in2 in3 ex1)
\n )
\n )<\/p>\n
(foreach x HCs
\n (setq isAll (cadddr x))
\n (setq Pt1 (car x))
\n (setq Pt2 (cadr x))
\n (setq pt3 (caddr x))<\/p>\n
(setq cc1 (car cen))
\n (setq cc2 (cadr cen))
\n (setq cc3 (caddr cen))<\/p>\n
(setq rr1 (car rad))
\n (setq rr2 (cadr rad))
\n (setq rr3 (caddr rad))<\/p>\n
(setq cen (append (cdr cen) (list (car cen))))
\n (setq rad (append (cdr rad) (list (car rad))))<\/p>\n
(if (and (= (car Pt1) 1e400) (= (car Pt2) 1e400))
\n (setq pc1 cc1 pc2 cc2 pc3 cc3)
\n (progn
\n (if (= (car Pt1) 1e400)
\n (setq Ln1 (Line:Equation Pt2 pt3))
\n (if (= (car Pt2) 1e400)
\n (setq Ln1 (Line:Equation pt3 Pt1))
\n (setq Ln1 (Line:Equation Pt1 Pt2))
\n )
\n )
\n (setq kA1 (car Ln1)
\n kB1 (cadr ln1)
\n kC1 (caddr ln1)
\n pc1 (CIR:Polar_Point cc1 rr1 kA1 kB1 kC1)
\n pc2 (CIR:Polar_Point cc2 rr2 kA1 kB1 kC1)
\n pc3 (CIR:Polar_Point cc3 rr3 kA1 kB1 kC1)
\n )
\n )
\n )
\n (setq Ps1 (CIR:Inters_Circle_Line cc1 rr1 CCC PC1))
\n (setq Ps2 (CIR:Inters_Circle_Line cc2 rr2 CCC pc2))
\n (setq ps3 (CIR:Inters_Circle_Line cc3 rr3 CCC pc3))<\/p>\n
(setq mt1 (car ps1))
\n (setq mt2 (car ps2))
\n (setq mt3 (car ps3))<\/p>\n
(setq nt1 (cadr ps1))
\n (setq nt2 (cadr ps2))
\n (setq nt3 (cadr ps3))<\/p>\n
(if (and nt1 nt2 nt3)
\n (if isAll
\n (setq RET (cons (CIR:PPP mt1 mt2 mt3) RET)
\n RET (cons (CIR:PPP nt1 nt2 nt3) RET)
\n )
\n (setq RET (cons (CIR:PPP mt1 nt2 nt3) RET)
\n RET (cons (CIR:PPP nt1 mt2 mt3) RET)
\n )
\n )
\n (if (and mt1 mt2 mt3)
\n (progn
\n (setq nt1 mt3)
\n (and (equal mt1 mt3 eps) (setq nt1 mt2))
\n (setq RET (append (CIR:PPC mt1 nt1 cc1 rr1) RET))
\n )
\n )
\n )
\n )
\n )
\n )
\n RET
\n)<\/p>\n
;;;—————————————————-;
\n;;;\u6309\u7167\u67d0\u4e2a\u51fd\u6570\u5347\u5e8f\u6392\u5e8f ;
\n;;;—————————————————-;
\n(defun MISC:ASort (l f)
\n (vl-sort l (function (lambda (e1 e2) (< (f e1) (f e2)))))
\n)<\/p>\n
;;;—————————————————-;
\n;;;\u6309\u7167\u67d0\u4e2a\u51fd\u6570\u964d\u5e8f\u6392\u5e8f ;
\n;;;—————————————————-;
\n(defun MISC:DSort (l f)
\n (vl-sort l (function (lambda (e1 e2) (> (f e1) (f e2)))))
\n)<\/p>\n
;;;—————————————————-;
\n;;;\u53cd\u53f7 ;
\n;;;—————————————————-;
\n(defun MISC:RevSign (f)
\n (if (eq f +) – +)
\n)<\/p>\n
;;;—————————————————-;
\n;;;\u5207\u5207\u5207\u7684\u5706(\u4e0e\u4e09\u4e2a\u5706\u90fd\u76f8\u5207\u7684\u5706) \u53e6\u4e00\u79cd\u65b9\u6cd5 ;
\n;;;\u8f93\u5165: \u4e09\u4e2a\u7ed9\u5b9a\u7684\u5706\u7684\u5706\u5fc3\u53ca\u534a\u5f84 ;
\n;;;\u8f93\u51fa: \u7b26\u5408\u6761\u4ef6\u7684\u5706\u7684\u96c6\u5408(nil\uff0c\u4e00\u4e2a\u6216\u591a\u4e2a\u5706\u7684\u96c6\u5408) ;
\n;;;\u6ce8\u660e: \u6b64\u79cd\u4f5c\u6cd5\u5728\u67d0\u4e9b\u60c5\u51b5\u4e0b\u4e0d\u5b8c\u5168\uff0c\u8fd8\u662f\u4ee5\u524d\u9762\u65b9\u6cd5\u4e3a\u51c6;
\n;;;—————————————————-;
\n(defun CIR:CCC_1 (c1 r1 c2 r2 c3 r3 \/ Ret C D1 D2 EPS LST R S1 S2)
\n (setq eps 1e-6)
\n (setq lst (mapcar ‘cons (list c1 c2 c3) (list r1 r2 r3)))
\n (setq lst (MISC:ASort lst cdr))
\n (setq c1 (caar lst))
\n (setq c2 (caadr lst))
\n (setq c3 (caaddr lst))
\n (setq r1 (cdar lst))
\n (setq r2 (cdadr lst))
\n (setq r3 (cdaddr lst))
\n (foreach f (mapcar ‘cons (list + + – -) (list + – + -))
\n (setq s1 (car f))
\n (setq s2 (cdr f))
\n (foreach p (CIR:PCC c1 c2 (s1 r2 r1) c3 (s2 r3 r1))
\n (setq c (car p))
\n (setq r (cadr p))
\n (setq d1 (distance c c2))
\n (setq d2 (distance c c3))
\n (cond
\n ( (equal (abs (s1 d1 r2)) (abs (s2 d2 r3)) eps)
\n (setq ret (cons (list c (- r r1)) ret))
\n )
\n ( (equal (abs ((MISC:RevSign s1) d1 r2))
\n (abs ((MISC:RevSign s2) d2 r3))
\n eps
\n )
\n (setq ret (cons (list c (+ r r1)) ret))
\n )
\n )
\n )
\n )
\n ret
\n)<\/p>\n
;;;—————————————————-;
\n;;;\u70b9\u70b9Per\u753b\u5706 ;
\n;;;—————————————————-;
\n(defun CIR:PPT (Pt1 Pt2 PtM PtN \/ eps mid pt3 cen)
\n (setq eps 1e-8)
\n (setq Mid (Geo:MidPoint Pt1 Pt2))
\n (setq Pt3 (GEO:Rot90 Mid Pt1 Pt2))
\n (if (setq cen (inters Mid Pt3 PtM PtN nil))
\n (list (list cen (distance Cen Pt1)))
\n (if (LINE:Colinearity Mid PtM PtN)
\n (list (list Mid (distance Mid Pt1)))
\n )
\n )
\n)<\/p>\n
;;;—————————————————-;
\n;;;\u70b9\u5207\u7ebfPer\u753b\u5706 ;
\n;;;—————————————————-;
\n(defun CIR:PLT (Pt0 PtA PtB PtM PtN \/ eps ang Pt1)
\n (setq eps 1e-8)
\n (setq ang (angle PtM PtN))
\n (setq Pt1 (GEO:Mirror2D Pt0 PtM ang))
\n (if (equal pt1 pt0 eps)
\n (CIR:PLL Pt0 PtA PtB (GEO:Mirror2D PtA PtM ang) (GEO:Mirror2D PtB PtM ang))
\n (CIR:PPL Pt0 Pt1 PtA PtB)
\n )
\n)<\/p>\n
;;;—————————————————-;
\n;;;\u70b9\u5207\u5706Per\u753b\u5706 ;
\n;;;—————————————————-;
\n(defun CIR:PCT (Pt0 Cen Rad PM PN \/ eps AN C0 C1 C2 D1 D2 D3 D4 IsA IsB L PT PT1 S C R R1 R2)
\n (setq eps 1e-8)
\n (setq Pt1 (GEO:Mirror3D Pt0 PM PN))
\n (setq C0 (GEO:Mirror3D Cen PM PN))
\n (setq IsA (equal Pt0 Pt1 eps))
\n (setq IsB (equal Cen C0 eps))
\n (if (and IsA IsB)
\n (progn
\n (setq d1 (distance Pt0 Cen))
\n (if (equal d1 0 eps)
\n (setq an (angle PM PN))
\n (setq an (angle Pt0 cen))
\n )
\n (setq R1 (* (+ d1 rad) 0.5))
\n (setq R2 (* (- d1 rad) 0.5))
\n (setq C1 (polar Pt0 an R1))
\n (setq C2 (polar Pt0 an R2))
\n (list (list C1 R1) (list C2 (abs R2)))
\n )
\n (if IsA
\n (progn
\n (and (equal Pt0 PM eps) (setq Pt PM PM PN PN Pt))
\n (setq an (angle PM Pt0))
\n (setq Pt (trans (mapcar ‘- Cen Pt0) 0 (mapcar ‘- Pt0 pM)))
\n (setq d1 (car pt))
\n (setq d2 (caddr Pt))
\n (setq d3 (- d2 rad))
\n (setq d4 (+ d2 rad))
\n (setq L (+ (* d1 d1) (* d3 d4)))
\n (foreach d (list (+ d3 d3) (+ d4 d4))
\n (if (not (equal d 0 eps))
\n (setq r (\/ L d)
\n c (polar Pt0 an (\/ L d))
\n S (cons (list c (abs r)) S)
\n )
\n )
\n )
\n S
\n )
\n (CIR:PPC Pt0 Pt1 Cen Rad)
\n )
\n )
\n)<\/p>\n
;;;—————————————————-;
\n;;;\u5207\u7ebf\u5207\u7ebfPer\u753b\u5706 ;
\n;;;—————————————————-;
\n(defun CIR:LLT (PA PB PC PD PM PN \/ eps int p1 p2 an cen per ret)
\n (setq eps 1e-8)
\n (if (setq int (inters PA PB PC PD nil))
\n (progn
\n (setq an (* (+ (angle pa pb) (angle pc pd)) 0.5))
\n (foreach a (list an (+ an (* pi 0.5)))
\n (setq p1 (polar int a 100))
\n (if (setq cen (inters int p1 PM PN nil))
\n (setq Per (LINE:Perpendicular_Foot Cen PA PB)
\n ret (cons (list cen (abs (car per))) ret)
\n )
\n )
\n )
\n ret
\n )
\n (progn
\n (setq P1 (GEO:Midpoint PA PC))
\n (setq P2 (polar p1 (angle Pa Pb) 100))
\n (setq Cen (inters P1 P2 PM PN nil))
\n (setq Per (LINE:Perpendicular_Foot Cen PA PB))
\n (setq ret (list (list cen (abs (car per)))))
\n )
\n )
\n)<\/p>\n
;;;—————————————————-;
\n;;;\u5207\u7ebf\u5207\u5706Per\u753b\u5706 ;
\n;;;—————————————————-;
\n(defun CIR:LCT (PtA PtB Cen Rad PtM PtN \/ Pt1 Pt2)
\n (setq Pt1 (GEO:Mirror3d PtA PtM PtN))
\n (setq Pt2 (GEO:Mirror3d PtB PtM PtN))
\n (vl-remove-if-not
\n (function (lambda (x) (LINE:Colinearity (car x) PtM PtN)))
\n (CIR:LLC PtA PtB Pt1 Pt2 Cen Rad)
\n )
\n)<\/p>\n
;;;—————————————————-;
\n;;;\u4e24\u4e2a\u5706\u7684\u4f4d\u4f3c\u4e2d\u5fc3 Internal,External Homothetic Center;
\n;;;\u4e5f\u5c31\u662f\u516c\u5207\u7ebf\u7684\u4ea4\u70b9(\u5f53\u7136\u4e5f\u5305\u62ec\u6ca1\u6709\u516c\u5207\u7ebf\u65f6\u5019\u7684\u60c5\u51b5) ;
\n;;;\u8f93\u5165: \u4e24\u4e2a\u5706\u7684\u5706\u5fc3\u548c\u534a\u5f84 ;
\n;;;\u8f93\u51fa: \u7b2c\u4e00\u4e2a\u4e3a\u5185\u4f4d\u4f3c\u4e2d\u5fc3\uff0c\u7b2c\u4e8c\u4e2a\u4e3a\u5916\u4f4d\u4f3c\u4e2d\u5fc3 ;
\n;;;—————————————————-;
\n(defun CIR:Homothetic_Center (c1 r1 c2 r2 \/)
\n (if (equal r1 r2 1e-14)
\n (list (GEO:Midpoint c1 c2) ‘(1e400 1e400 0)) ;\u4e2d\u70b9\u548c\u65e0\u7a77\u8fdc\u70b9
\n (list (GEO:Scale c2 c1 (\/ r1 (+ r1 r2))) ;\u5185\u4f4d\u4f3c\u4e2d\u5fc3
\n (GEO:Scale c2 c1 (\/ r1 (- r1 r2))) ;\u5916\u4f4d\u4f3c\u4e2d\u5fc3
\n )
\n )
\n)<\/p>\n
;;;—————————————————-;
\n;;;\u70b9\u7684\u53cd\u6f14 ;
\n;;;\u8f93\u5165: \u5706\u5fc3\uff0c\u534a\u5f84\u548c\u4e00\u70b9 ;
\n;;;\u8f93\u51fa: nil\u6216\u8005\u53cd\u6f14\u70b9 ;
\n;;;—————————————————-;
\n(defun CIR:Inversion (c r p \/ d)
\n (setq d (distance c p))
\n (if (equal d 0 1e-8)
\n nil ;\u5706\u5fc3\u5904\u6ca1\u6709\u53cd\u6f14\u70b9
\n (polar c (angle c p) (\/ (* r r) d)) ;\u6839\u636e\u53cd\u6f14\u516c\u5f0f\u8ba1\u7b97
\n )
\n)<\/p>\n
;;;—————————————————-;
\n;;;\u76f4\u7ebf\u5bf9\u5706\u7684\u53cd\u6f14 ;
\n;;;\u8f93\u5165: \u5706\u5fc3\uff0c\u534a\u5f84\u548c\u4e00\u70b9 ;
\n;;;\u8f93\u51fa: \u4e00\u4e2a\u7528\u5706\u5fc3\u548c\u534a\u5f84\u8868\u793a\u7684\u5706\u6216\u8005\u76f4\u7ebf\u672c\u8eab ;
\n;;;—————————————————-;
\n(defun CIR:Line_Inversion (Cen Rad P1 P2 \/ d p q c r)
\n (setq p (Line:Perpendicular_Foot Cen P1 P2))
\n (setq d (car p))
\n (setq p (cadr p))
\n (if (equal d 0.0 1e-8)
\n (list P1 P2) ;\u7ecf\u8fc7\u5706\u5fc3\u7684\u76f4\u7ebf\u662f\u672c\u8eab
\n (setq q (CIR:Inversion Cen Rad p)
\n c (GEO:Midpoint q Cen)
\n r (list c (distance c Cen)) ;\u628a\u76f4\u7ebf\u53cd\u6f14\u6210\u5706
\n )
\n )
\n)<\/p>\n
;;;—————————————————-;
\n;;;\u5706\u5bf9\u5706\u7684\u53cd\u6f14 ;
\n;;;\u8f93\u5165: \u53cd\u6f14\u5706\u7684\u5706\u5fc3\u534a\u5f84\u548c\u88ab\u53cd\u6f14\u7684\u5706\u7684\u5706\u5fc3\uff0c\u534a\u5f84 ;
\n;;;\u8f93\u51fa: \u4e00\u4e2a\u7528\u5706\u5fc3\u548c\u534a\u5f84\u8868\u793a\u7684\u5706\u6216\u8005\u4e00\u6761\u76f4\u7ebf\u7684\u4e24\u4e2a\u7aef\u70b9;
\n;;;—————————————————-;
\n(defun CIR:Circle_Inversion (C0 R0 C R \/ an P1 P2 PM)
\n (if (equal C0 C 1e-8)
\n (list C0 (* (\/ R0 R) R0))
\n (progn
\n (setq an (angle C C0))
\n (setq p1 (CIR:Inversion C0 R0 (polar C an R)))
\n (setq p2 (CIR:Inversion C0 R0 (polar C an (- R))))
\n (if (and p1 p2)
\n (list (GEO:Midpoint p1 p2) (* (distance p1 p2) 0.5))
\n (progn
\n (setq PM (GEO:Midpoint C0 (CIR:Inversion C0 R0 C)))
\n (list PM (mapcar ‘+ PM (MAT:Rot90 (mapcar ‘- PM C0))))
\n )
\n )
\n )
\n )
\n)<\/p>\n
;;;—————————————————-;
\n;;;\u7ebf\u6bb5(\u7528\u7cfb\u6570\u8868\u793a)\u5bf9\u5706\u7684\u6781\u70b9 ;
\n;;;—————————————————-;
\n(defun CIR:Polar_Point (cen rad A B C \/ p d)
\n (setq p (Line:Perpendicular_Distance_1 Cen A B C))
\n (setq d (car p))
\n (setq p (cadr p))
\n (if (\/= d 0)
\n (polar cen (angle cen p) (\/ (* rad rad) (abs d)))
\n )
\n)<\/p>\n
;;;—————————————————-;
\n;;;\u7ebf\u6bb5(\u7528\u4e24\u70b9\u8868\u793a)\u5bf9\u5706\u7684\u6781\u70b9 ;
\n;;;—————————————————-;
\n(defun CIR:Polar_Point_1 (cen rad Pa Pb \/ p d)
\n (setq p (Line:Perpendicular_Foot Cen Pa Pb))
\n (setq d (car p))
\n (setq p (cadr p))
\n (if (\/= d 0)
\n (polar cen (angle cen p) (\/ (* rad rad) (abs d)))
\n )
\n)<\/p>\n
;;;—————————————————-;
\n;;;\u7b49\u5e42\u8f74 ;
\n;;;\u8f93\u5165: \u7ed9\u5b9a\u7684\u4e24\u4e2a\u5706\u7684\u5706\u5fc3\u548c\u534a\u5f84 ;
\n;;;\u8f93\u51fa: \u8fd9\u4e24\u4e2a\u5706\u7684\u7b49\u5e42\u8f74\u4e0e\u5b83\u4eec\u7684\u5706\u5fc3\u8fde\u7ebf\u7684\u4ea4\u70b9 ;
\n;;;—————————————————-;
\n(defun CIR:Radical_Axis (c1 r1 c2 r2 \/ l d)
\n (setq l (distance c1 c2))
\n (if (equal l 0 1e-8)
\n (setq d (if (> r1 r2) 1e400 -1e400))
\n (setq d (* 0.5 (+ L (\/ (* (+ r1 r2) (- r1 r2)) L))))
\n )
\n (polar c1 (angle c1 c2) d)
\n)<\/p>\n
;;;—————————————————-;
\n;;;\u7b49\u5e42\u4e2d\u5fc3(\u4e09\u4e2a\u5706\u5f62\u7684\u7b49\u5e42\u8f74\u7684\u4ea4\u70b9) Radical center, ;
\n;;;also called the power center of three circles ;
\n;;;\u8f93\u5165: \u4e09\u4e2a\u5706\uff08\u7528\u4e2d\u5fc3\u534a\u5f84\u8868\u793a\uff09 ;
\n;;;\u8f93\u51fa: \u8fd9\u4e09\u4e2a\u5706\u7684\u7b49\u5e42\u4e2d\u5fc3 ;
\n;;;—————————————————-;
\n(defun CIR:Power_Center (c1 r1 c2 r2 c3 r3 \/ p1 p2 q1 q2 CC d k)
\n (setq P1 (CIR:Radical_Axis c1 r1 c2 r2))
\n (setq P2 (CIR:Radical_Axis c2 r2 c3 r3))
\n (if (and p1 p2)
\n (progn
\n (setq q1 (GEO:Rot90 p1 c1 c2))
\n (setq q2 (GEO:Rot90 p2 c2 c3))
\n (setq CC (inters p1 q1 p2 q2 nil))
\n (if CC
\n (progn
\n (setq d (distance CC C1))
\n (setq k (* (+ d r1) (- d r1)))
\n (list cc K)
\n )
\n )
\n )
\n )
\n)<\/p>\n
;;;—————————————————-;
\n;;;\u963f\u6ce2\u7f57\u5c3c\u5965\u65af\u5706Apollonian circles ;
\n;;;—————————————————-;
\n(defun CIR:Apollonian_Circle (P1 P2 k \/ k k1 k2 d rt)
\n (if (and (> k 0) (\/= k 1))
\n (setq k1 (* k k)
\n k2 (\/ 1.0 (1- k1))
\n d (distance p1 p2)
\n Rt (list (mapcar (function (lambda (i j) (* k2 (- (* k1 j) i)))) p1 p2)
\n (abs (* d k k2))
\n )
\n )
\n )
\n)<\/p>\n
;;;—————————————————-;
\n;;;\u5171\u8f74\u5706 ;
\n;;;—————————————————-;
\n(defun CIR:Coaxal_Circle (c1 r1 c2 r2 \/ P1 P2 d1 d2 RT)
\n (foreach HC (CIR:Homothetic_Center c1 r1 c2 r2)
\n (if (= (car HC) 1e400)
\n (setq P1 (Geo:MidPoint C1 C2)
\n P2 (polar P1 (+ (angle C1 C2) (* pi 0.5)) 1000)
\n RT (cons (list P1 P2) RT)
\n )
\n (setq p1 (polar C1 (angle HC C1) r1)
\n p2 (polar c2 (angle HC C2) (- r2))
\n d1 (distance HC p1)
\n d2 (distance HC p2)
\n RT (cons (list HC (sqrt (* d1 d2))) RT)
\n )
\n )
\n )
\n)<\/p>\n
;;;—————————————————-;
\n;;;\u5171\u8f74\u5706\u7684\u7279\u6b8a\u60c5\u51b5\uff08\u5176\u4e2d\u4e00\u4e2a\u5706\u9000\u5316\u4e3a\u4e00\u70b9\uff09 ;
\n;;;—————————————————-;
\n(defun CIR:Radical_Axis_Point (C R P \/ l d)
\n (if (equal c p 1e-8)
\n nil
\n (progn
\n (setq l (distance c p))
\n (setq d (\/ (+ (* r r) (* l l)) 2 l))
\n (polar c (angle c p) d)
\n )
\n )
\n)<\/p>\n
;;;—————————————————-;
\n;;;\u5361\u65af\u8482\u90ce\u6c42\u89e3 ;
\n;;;Castillon’s Problem ;
\n;;;\u8f93\u5165: \u4e00\u7ed9\u5b9a\u5706\uff08\u5706\u5fc3\u534a\u5f84\u8868\u793a\uff09\u548c\u4e0d\u91cd\u5408\u7684\u4e09\u5b9a\u70b9\u3002 ;
\n;;;\u8f93\u51fa\uff1a\u5706\u4e0a\u4e09\u70b9\uff0c\u4f7f\u5f97\u4e09\u5df2\u77e5\u70b9\u5206\u522b\u901a\u8fc7\u8fd9\u4e09\u70b9\u5f62\u6210\u7684\u8fb9 ;
\n;;;—————————————————-;
\n(defun CIR:Castillon (Cen Rad pa pb pc \/ EPS INT1 INT2 L d1 d2 d3 LST P0 P1 P2 x0 y0 RET)
\n (defun Check (Func lst d eps)
\n (apply Func (mapcar (function (lambda (x) (equal x d eps))) lst))
\n )
\n (setq eps 1e-6)
\n (setq d1 (abs (- rad (distance pa cen))))
\n (setq d2 (abs (- rad (distance pb cen))))
\n (setq d3 (abs (- rad (distance pc cen))))
\n (if (Check ‘or (list d1 d2 d3) 0 eps)
\n (if (Check ‘and (list d1 d2 d3) 0 eps)
\n (list (list pa pb pc))
\n (setq l (vl-sort-i (list d1 d2 d3) ‘<)
\n lst (list pa pb pc)
\n pa (nth (car l) lst)
\n pb (nth (cadr l) lst)
\n pc (nth (caddr l) lst)
\n p0 Pa
\n p1 (CIR:Intersecting_Chords P0 pb Cen Rad)
\n p2 (CIR:Intersecting_Chords P1 PC Cen Rad)
\n ret (cons (list p0 p1 p2) ret)
\n p1 (CIR:Intersecting_Chords P0 pc Cen Rad)
\n p2 (CIR:Intersecting_Chords P1 pb Cen Rad)
\n ret (cons (list P0 p1 p2) ret)
\n )
\n )
\n (progn
\n (foreach an (list 1 2 3)
\n (setq x0 (+ (* rad (cos an)) (car Cen)))
\n (setq y0 (+ (* rad (sin an)) (cadr Cen)))
\n (setq p0 (list x0 y0 0))
\n (setq p1 p0)
\n (foreach p (list pa pb pc)
\n (setq P2 (CIR:Intersecting_Chords P1 P cen rad))
\n (setq p1 p2)
\n )
\n (setq lst (cons (cons p0 p2) lst))
\n )
\n ;;\u5230\u8fd9\u6b65\u5f97\u5230\u4e09\u4e2a\u4e8c\u91cd\u6620\u5c04\u70b9\u3002
\n (setq int1 (inters (caar lst) (cdadr lst) (cdar lst) (caadr lst) nil))
\n (setq int2 (inters (caar lst) (cdaddr lst) (cdar lst) (caaddr lst) nil))
\n (if (and int1 int2)
\n (foreach q (CIR:Inters_Circle_Line cen rad int1 int2)
\n (setq l (list q))
\n (setq p1 q)
\n (foreach p (list pa pb)
\n (setq p2 (CIR:Intersecting_Chords p1 p cen rad))
\n (setq l (cons p2 l))
\n (setq p1 p2)
\n )
\n (setq ret (cons (reverse l) ret))
\n )
\n )
\n )
\n )
\n)<\/p>\n
;;|—————————————————-;
\n;;;\u4ee5\u4e0b\u6837\u4f8b\u4ec5\u4f9b\u6d4b\u8bd5\u3002 ;
\n;;;—————————————————-;<\/p>\n
;;;—————————————————-;
\n;;;A sample for Apollonian circles ;
\n;;;\u963f\u6ce2\u7f57\u5c3c\u65af\u5706\u7684\u6d4b\u8bd5\u6837\u4f8b ;
\n;;;—————————————————-;
\n(defun c:cac (\/ p1 p2 k rt)
\n (initget 1)
\n (setq p1 (getpoint “n\u70b91:”))
\n (initget 2)
\n (setq p2 (getpoint p1 “n\u70b92:”))
\n (initget 7)
\n (setq k (getreal “n\u6bd4\u4f8b:”))
\n (Ent:Make_Line p1 P2)
\n (setq rt (CIR:Apollonian_Circle P1 P2 k))
\n (and rt (apply ‘Ent:Make_Circle rt))
\n)<\/p>\n
;;;—————————————————-;
\n;;;A sample for CIR:coaxl_Circle ;
\n;;;\u5171\u8f74\u5706\u7684\u6d4b\u8bd5\u6837\u4f8b ;
\n;;;—————————————————-;
\n(defun C:CXC (\/ sel d1 d2 c1 c2 r1 r1 rt)
\n (setq sel (ssget ‘((0 . “CIRCLE,ARC”))))
\n (if (and sel (>= (sslength sel) 2))
\n (progn
\n (setq d1 (entget (ssname sel 0)))
\n (setq d2 (entget (ssname sel 1)))
\n (setq c1 (cdr (assoc 10 d1)))
\n (setq c2 (cdr (assoc 10 d2)))
\n (setq r1 (cdr (assoc 40 d1)))
\n (setq r2 (cdr (assoc 40 d2)))
\n (setq rt (CIR:Coaxal_Circle c1 r1 c2 r2))
\n (foreach c rt
\n (if (numberp (cadr c))
\n (apply ‘Ent:Make_Circle c)
\n (apply ‘Ent:Make_Line c)
\n )
\n )
\n )
\n )
\n)<\/p>\n
;;;—————————————————-;
\n;;;A sample for CIR:Radical_Axis_Point ;
\n;;;\u5171\u8f74\u5706\u7279\u4f8b\u7684\u6d4b\u8bd5\u6837 ;
\n;;;—————————————————-;
\n(defun c:rax1(\/ e p d r c ret)
\n (setq e (car (entsel “n\u57061:”)))
\n (setq p (getpoint “n\u70b9:”))
\n (if (and e p)
\n (progn
\n (setq d (entget e))
\n (setq r (cdr (assoc 40 d))
\n c (cdr (assoc 10 d))
\n )
\n (Ent:Make_Point p)
\n (setq ret (CIR:Radical_Axis_Point C R p))
\n (and ret (Ent:Make_Point ret))
\n )
\n )
\n)<\/p>\n
;;;A sample for Castillon’s problem
\n(defun C:Castillon(\/ ss pa pb pc ent dxf rad cen ret)
\n (setq ss (ssget “:S” ‘((0 . “CIRCLE,ARC”))))
\n (initget 1)
\n (setq pa (getpoint “n\u70b91:”))
\n (initget 2)
\n (setq pb (getpoint “n\u70b92:”))
\n (initget 3)
\n (setq pc (getpoint “n\u70b93:”))
\n (if ss
\n (progn
\n (setq ent (ssname ss 0))
\n (setq dxf (entget ent))
\n (setq rad (cdr (assoc 40 dxf)))
\n (setq cen (cdr (assoc 10 dxf)))
\n (mapcar ‘Ent:Make_Point (list pa pb pc))
\n (setq ret (CIR:Castillon Cen Rad pa pb pc))
\n (foreach p ret
\n (apply ‘Ent:Make_Triangle p)
\n )
\n )
\n )
\n)<\/p>\n
;;;\u6d4b\u8bd5\u70b9\u70b9Per\u753b\u5706
\n(defun C:PPT (\/ s1 p1 p2 dxf pta ptb ret)
\n (prompt “n\u9009\u53d6\u7ebf\u6bb5: “)
\n (setq s1 (ssget “:S” ‘((0 . “LINE”))))
\n (initget 1)
\n (setq p1 (getpoint “n\u70b91:”))
\n (initget 1)
\n (setq p2 (getpoint “n\u70b91:”))
\n (if (and s1 p1 P2)
\n (progn
\n (setq dxf (entget (ssname s1 0)))
\n (setq ptA (cdr (assoc 10 dxf)))
\n (setq ptB (cdr (assoc 11 dxf)))
\n (grdraw p1 p2 1)
\n (setq ret (CIR:PPT P1 P2 PtA PtB))
\n (foreach p Ret
\n (apply ‘ENt:Make_circle p)
\n )
\n )
\n )
\n)<\/p>\n
;;;\u6d4b\u8bd5\u70b9\u70b9Per\u753b\u5706
\n(defun C:PLT (\/ sel Pt0 d1 d2 PtA PtB PtC PtD Ret)
\n (prompt “n\u9009\u53d6\u7ebf\u6bb5: “)
\n (setq sel (ssget ‘((0 . “LINE”))))
\n (initget 1)
\n (setq Pt0 (getpoint “n\u70b91:”))
\n (if (and sel Pt0 (>= (sslength sel) 2))
\n (progn
\n (setq d1 (entget (ssname sel 0)))
\n (setq d2 (entget (ssname sel 1)))
\n (setq ptA (cdr (assoc 10 d1)))
\n (setq ptB (cdr (assoc 11 d1)))
\n (setq ptC (cdr (assoc 10 d2)))
\n (setq ptD (cdr (assoc 11 d2)))
\n (setq ret (CIR:PLT Pt0 PtA PtB PtC PtD))
\n (foreach p Ret
\n (apply ‘ENt:Make_circle p)
\n )
\n )
\n )
\n)
\n(defun c:LCC(\/ s1 s2 d1 d2 d3 pa pb c1 r1 c2 r2 ret)
\n (prompt “n\u9009\u7ebf: “)
\n (setq s1 (ssget “_+.:E:S:L” ‘((0 . “LINE”))))
\n (prompt “n\u9009\u5706: “)
\n (setq s2 (ssget ‘((0 . “ARC,CIRCLE”))))
\n (if (and s1 s2 (>= (sslength s2) 2))
\n (progn
\n (setq d1 (entget (ssname s1 0)))
\n (setq d2 (entget (ssname s2 0)))
\n (setq d3 (entget (ssname s2 1)))
\n (setq pA (cdr (assoc 10 d1)))
\n (setq pB (cdr (assoc 11 d1)))
\n (setq C1 (cdr (assoc 10 d2)))
\n (setq R1 (cdr (assoc 40 d2)))
\n (setq C2 (cdr (assoc 10 d3)))
\n (setq R2 (cdr (assoc 40 d3)))
\n (command “undo” “be”)
\n (setq ReT (CIR:LCC Pa Pb C1 R1 C2 R2))
\n (foreach p ReT
\n (apply ‘ENt:Make_circle p)
\n )
\n (command “undo” “e”)
\n )
\n )
\n (princ)
\n)<\/p>\n
;;;\u6d4b\u8bd5\u70b9\u70b9Per\u753b\u5706
\n(defun C:LLT (\/ sel Pt0 d1 d2 PtA PtB PtC PtD Ret)
\n (prompt “n\u9009\u53d6\u7ebf\u6bb5: “)
\n (setq sel (ssget ‘((0 . “LINE”))))
\n (if (and sel (>= (sslength sel) 3))
\n (progn
\n (setq d1 (entget (ssname sel 0)))
\n (setq d2 (entget (ssname sel 1)))
\n (setq d3 (entget (ssname sel 2)))
\n (setq pA (cdr (assoc 10 d1)))
\n (setq pB (cdr (assoc 11 d1)))
\n (setq pC (cdr (assoc 10 d2)))
\n (setq pD (cdr (assoc 11 d2)))
\n (setq pM (cdr (assoc 10 d3)))
\n (setq pN (cdr (assoc 11 d3)))
\n (setq RT (CIR:LLT PA PB PC PD PM PN))
\n (foreach p RT
\n (apply ‘ENt:Make_circle p)
\n )
\n )
\n )
\n)<\/p>\n
;;;\u6d4b\u8bd5\u5207\u7ebf\u5207\u5706Per\u753b\u5706
\n(defun C:PCT (\/ ss1 ss2 pt0 dx1 dx2 ptM ptN Cen Rad Ret)
\n (prompt “n\u9009\u7ebf: “)
\n (setq ss1 (ssget “_+.:E:S:L” ‘((0 . “LINE”))))
\n (prompt “n\u9009\u5706: “)
\n (setq ss2 (ssget “_+.:E:S:L” ‘((0 . “ARC,CIRCLE”))))
\n (initget 1)
\n (setq Pt0 (getpoint “n\u9009\u53d6\u70b9:”))
\n (if (and ss1 ss2)
\n (progn
\n (setq dx1 (entget (ssname ss1 0)))
\n (setq dx2 (entget (ssname ss2 0)))
\n (setq ptM (cdr (assoc 10 dx1)))
\n (setq ptN (cdr (assoc 11 dx1)))
\n (setq cen (cdr (assoc 10 dx2)))
\n (setq rad (cdr (assoc 40 dx2)))
\n (Ent:Make_Point Pt0)
\n (setq Ret (CIR:PCT Pt0 Cen Rad PtM PtN))
\n (foreach p Ret
\n (apply ‘ENt:Make_circle p)
\n )
\n )
\n )
\n)<\/p>\n
;;;\u6d4b\u8bd5\u5207\u7ebf\u5207\u5706Per\u753b\u5706
\n(defun C:LCT (\/ s1 s2 d1 d2 d3 pa pb pc pd c0 r0 ret)
\n (prompt “n\u9009\u7ebf: “)
\n (setq s1 (ssget ‘((0 . “LINE”))))
\n (prompt “n\u9009\u5706: “)
\n (setq s2 (ssget “_+.:E:S:L” ‘((0 . “ARC,CIRCLE”))))
\n (if (and s1 s2 (>= (sslength s1) 2))
\n (progn
\n (setq d1 (entget (ssname s1 0)))
\n (setq d2 (entget (ssname s1 1)))
\n (setq d3 (entget (ssname s2 0)))
\n (setq pA (cdr (assoc 10 d1)))
\n (setq pB (cdr (assoc 11 d1)))
\n (setq pC (cdr (assoc 10 d2)))
\n (setq pD (cdr (assoc 11 d2)))
\n (setq c0 (cdr (assoc 10 d3)))
\n (setq r0 (cdr (assoc 40 d3)))
\n (command “undo” “be”)
\n (setq ret (CIR:LCT Pa Pb C0 R0 Pc Pd))
\n (foreach p ret
\n (apply ‘ENt:Make_circle p)
\n )
\n (command “undo” “e”)
\n )
\n )
\n (princ)
\n)
\n;;;\u6d4b\u8bd5\u5207\u7ebf\u5207\u5706Per\u753b\u5706
\n(defun C:LLC (\/ s1 s2 d1 d2 d3 pa pb pc pd c0 r0 ret)
\n (prompt “n\u9009\u7ebf: “)
\n (setq s1 (ssget ‘((0 . “LINE”))))
\n (prompt “n\u9009\u5706: “)
\n (setq s2 (ssget “_+.:E:S:L” ‘((0 . “ARC,CIRCLE”))))
\n (if (and s1 s2 (>= (sslength s1) 2))
\n (progn
\n (setq d1 (entget (ssname s1 0)))
\n (setq d2 (entget (ssname s1 1)))
\n (setq d3 (entget (ssname s2 0)))
\n (setq pA (cdr (assoc 10 d1)))
\n (setq pB (cdr (assoc 11 d1)))
\n (setq pC (cdr (assoc 10 d2)))
\n (setq pD (cdr (assoc 11 d2)))
\n (setq c0 (cdr (assoc 10 d3)))
\n (setq r0 (cdr (assoc 40 d3)))
\n (command “undo” “be”)
\n (setq ret (CIR:LLC PA PB PC PD C0 R0))
\n (foreach p ret
\n (apply ‘ENt:Make_circle p)
\n )
\n (command “undo” “e”)
\n )
\n )
\n (princ)
\n)<\/p>\n
;;; A sample for CIR:Inters_Circle_Line
\n;;; \u7ebf\u5706\u6c42\u4ea4\u793a\u4f8b
\n(defun C:CIL (\/ s1 s2 dx1 dx2 pta ptb cen rad ret)
\n (prompt “n\u9009\u53d6\u7ebf\u6bb5: “)
\n (setq s1 (ssget “:S” ‘((0 . “LINE”))))
\n (prompt “n\u9009\u53d6\u5706: “)
\n (setq s2 (ssget “:S” ‘((0 . “CIRCLE,ARC”))))
\n (if (and s1 s2)
\n (progn
\n (setq dx1 (entget (ssname s1 0)))
\n (setq dx2 (entget (ssname s2 0)))
\n (setq ptA (cdr (assoc 10 dx1)))
\n (setq ptB (cdr (assoc 11 dx1)))
\n (setq cen (cdr (assoc 10 dx2)))
\n (setq rad (cdr (assoc 40 dx2)))
\n (setq ret (CIR:Inters_Circle_Line cen rad PtA PtB))
\n (foreach p Ret
\n (Ent:Make_Point p)
\n )
\n )
\n )
\n)<\/p>\n
;;; A sample for CIR:LLR
\n;;; CIR:LLR\u6837\u4f8b
\n(defun c:LLR (\/ sel rad en1 en2 dx1 dx2 pta ptb ptc ptd ret)
\n (setq sel (ssget ‘((0 . “LINE”))))
\n (initget 1)
\n (setq rad (getdist “n\u534a\u5f84:”))
\n (if (and sel (>= (sslength sel) 2))
\n (progn
\n (setq en1 (ssname sel 0))
\n (setq en2 (ssname sel 1))
\n (setq dx1 (entget en1))
\n (setq dx2 (entget en2))
\n (setq PtA (cdr (assoc 10 d1)))
\n (setq PtB (cdr (assoc 11 d1)))
\n (setq PtC (cdr (assoc 10 d2)))
\n (setq PtD (cdr (assoc 11 d2)))
\n (if (setq ret (CIR:LLR PtA PtB PtC PtD Rad))
\n (foreach p ret
\n (apply ‘Ent:Make_Circle p)
\n )
\n )
\n )
\n )
\n)<\/p>\n
;;; A sample for CIR:PLR
\n;;; CIR:PLR\u6837\u4f8b
\n(defun c:PLR (\/ pt0 sel rad ent dxf pta ptb ret)
\n (initget 1)
\n (setq pt0 (getpoint “n1:”))
\n (initget 1)
\n (setq sel (ssget “:S” ‘((0 . “LINE”))))
\n (initget 7)
\n (setq rad (getdist “n\u534a\u5f84:”))
\n (if (and sel pt0 rad)
\n (progn
\n (setq ent (ssname sel 0))
\n (setq dxf (entget ent))
\n (setq PtA (cdr (assoc 10 dxf)))
\n (setq PtB (cdr (assoc 11 dxf)))
\n (if (setq ret (CIR:PLR Pt0 PtA PtB Rad))
\n (foreach p ret
\n (apply ‘Ent:Make_Circle p)
\n )
\n )
\n )
\n )
\n)<\/p>\n
;;; A sample for CIR:PPR
\n;;; CIR:PPR\u6837\u4f8b
\n(defun c:PPR (\/ pta ptb rad ret)
\n (initget 1)
\n (setq pta (getpoint “n1:”))
\n (initget 1)
\n (setq ptb (getpoint “n2:”))
\n (initget 7)
\n (setq rad (getdist “n\u534a\u5f84:”))
\n (if (setq ret (CIR:PPR pta ptb rad))
\n (foreach p ret
\n (apply ‘Ent:Make_Circle p)
\n )
\n )
\n)<\/p>\n
;;; A sample for CIR:LCR
\n;;; CIR:LCR\u6837\u4f8b
\n(defun c:LCR (\/ s1 s2 r1 d1 d2 pa pb c0 r0)
\n (prompt “n\u9009\u62e9\u76f4\u7ebf”)
\n (setq s1 (ssget “:S” ‘((0 . “LINE”))))
\n (prompt “n\u9009\u62e9\u5706”)
\n (setq s2 (ssget “:S” ‘((0 . “CIRCLE,ARC”))))
\n (initget 7)
\n (setq r1 (getdist “n\u534a\u5f84:”))
\n (if (and s1 s2 r1)
\n (progn
\n (setq d1 (entget (ssname s1 0)))
\n (setq d2 (entget (ssname s2 0)))
\n (setq pa (cdr (assoc 10 d1)))
\n (setq pb (cdr (assoc 11 d1)))
\n (setq c0 (cdr (assoc 10 d2)))
\n (setq r0 (cdr (assoc 40 d2)))
\n (foreach c (CIR:LCR PA PB C0 R0 R1)
\n (apply ‘Ent:Make_Circle c)
\n )
\n )
\n )
\n)<\/p>\n
;;; A sample for CIR:CCR
\n;;; CIR:CCR\u6837\u4f8b
\n(defun C:CCR (\/ s1 r0 d1 d2 c1 c2 r1 r2)
\n (prompt “n\u9009\u62e9\u4e24\u4e2a\u5706”)
\n (setq s1 (ssget ‘((0 . “CIRCLE,ARC”))))
\n (initget 7)
\n (setq r0 (getdist “n\u534a\u5f84:”))
\n (if (and s1 (>= (sslength s1) 2))
\n (progn
\n (setq d1 (entget (ssname s1 0)))
\n (setq d2 (entget (ssname s1 1)))
\n (setq c1 (cdr (assoc 10 d1)))
\n (setq r1 (cdr (assoc 40 d1)))
\n (setq c2 (cdr (assoc 10 d2)))
\n (setq r2 (cdr (assoc 40 d2)))
\n (foreach c (CIR:CCR C1 R1 C2 R2 R0)
\n (apply ‘Ent:Make_Circle c)
\n )
\n )
\n )
\n)<\/p>\n
;;; A sample for CIR:PCR
\n;;; CIR:PCR\u6837\u4f8b
\n(defun c:PCR (\/ p0 ss r1 dxf c0 r0)
\n (initget 1)
\n (setq p0 (getpoint “n\u70b9:”))
\n (setq p0 (trans p0 1 0))
\n (prompt “n\u9009\u62e9\u5706”)
\n (setq ss (ssget “:S” ‘((0 . “CIRCLE,ARC”))))
\n (initget 7)
\n (setq r1 (getdist “n\u534a\u5f84:”))
\n (if (and ss p0 r1)
\n (progn
\n (setq dxf (entget (ssname ss 0)))
\n (setq c0 (cdr (assoc 10 dxf)))
\n (setq r0 (cdr (assoc 40 dxf)))
\n (Ent:Make_Point P0)
\n (foreach c (CIR:PCR P0 C0 R0 R1)
\n (apply ‘Ent:Make_Circle c)
\n )
\n )
\n )
\n)<\/p>\n
;;;Sample
\n(defun C:CIV (\/ pnt cir dxf rad cen)
\n (initget 1)
\n (setq pnt (getpoint “n\u8bf7\u9009\u53d6\u4e00\u70b9:”))
\n (if (setq cir (car (entsel “n\u9009\u53d6\u5706\u6216\u8005\u5f27:”)))
\n (progn
\n (setq dxf (entget cir))
\n (setq rad (cdr (assoc 40 dxf)))
\n (setq cen (cdr (assoc 10 dxf)))
\n (Ent:Make_Point pnt)
\n (if (setq pt (CIR:Inversion cen rad pnt))
\n (Ent:Make_Point pt)
\n )
\n (foreach p (CIR:Point_Tangent cen rad pnt)
\n (Ent:Make_Point p)
\n )
\n )
\n )
\n)
\n;;;Sample
\n(defun C:CLI(\/ ent p1 p2 dxf rad cen ret)
\n (setq ent (car (entsel)))
\n (setq p1 (getpoint “nP1:”))
\n (setq p2 (getpoint “nP2:”))
\n (if (and ent p1 p2)
\n (progn
\n (setq dxf (entget ent))
\n (setq rad (cdr (assoc 40 dxf)))
\n (setq cen (cdr (assoc 10 dxf)))
\n (setq ret (CIR:Line_Inversion cen rad p1 p2))
\n (Ent:Make_Line p1 p2)
\n (apply ‘Ent:Make_Circle ret)
\n )
\n )
\n)<\/p>\n
;;;Sample\u5706\u5bf9\u5706\u7684\u53cd\u6f14
\n(defun C:CCI (\/ e0 e1 d0 d1 r0 r1 c0 c1 ret)
\n (setq e0 (car (entsel “n\u9009\u62e9\u57061:”)))
\n (setq e1 (car (entsel “n\u9009\u62e9\u57062:”)))
\n (if (and e0 e1)
\n (progn
\n (setq d0 (entget e0))
\n (setq d1 (entget e1))
\n (setq r0 (cdr (assoc 40 d0)))
\n (setq c0 (cdr (assoc 10 d0)))
\n (setq r1 (cdr (assoc 40 d1)))
\n (setq c1 (cdr (assoc 10 d1)))
\n (setq ret (CIR:Circle_Inversion C0 R0 C1 R1))
\n (if (numberp (cadr ret))
\n (apply ‘Ent:Make_Circle ret)
\n (apply ‘Ent:Make_Line ret)
\n )
\n )
\n )
\n)<\/p>\n
;;;A sample for CIR:Inters_Circle_Circle
\n;;;\u5706\u5706\u6c42\u4ea4\u70b9
\n;;;\u6d4b\u8bd5\u4f4d\u4f3c\u4e2d\u5fc3
\n(defun C:HC (\/ ent1 ent2 dxf1 dxf2 rad1 rad2 cen1 cen2)
\n (setq ent1 (car (entsel “n\u9009\u53d6\u57061:”)))
\n (setq ent2 (car (entsel “n\u9009\u53d6\u57062:”)))
\n ;(setq pt (getpoint “n\u70b9\u53d6:”))
\n (if (and ent1 ent2)
\n (progn
\n (setq dxf1 (entget ent1))
\n (setq dxf2 (entget ent2))
\n (setq rad1 (cdr (assoc 40 dxf1)))
\n (setq rad2 (cdr (assoc 40 dxf2)))
\n (setq cen1 (cdr (assoc 10 dxf1)))
\n (setq cen2 (cdr (assoc 10 dxf2)))
\n (setq i 1)
\n ;(Ent:Make_Point pt)
\n (setq ret1 (CIR:Radical_Axis cen1 rad1 cen2 rad2))
\n (setq ret2 (CIR:Common_Tangent cen1 rad1 cen2 rad2))
\n (setq ret3 (CIR:Homothetic_Center cen1 rad1 cen2 rad2))
\n (setq ret4 (CIR:Inters_Circle_Circle cen1 rad1 cen2 rad2))
\n (setq ret5 (CIR:Circle_Inversion cen1 rad1 cen2 rad2))
\n (if (numberp (cadr ret5))
\n (apply ‘Ent:Make_Circle ret5)
\n (apply ‘Ent:Make_Line ret5)
\n )
\n (and ret1 (Ent:Make_point ret1))
\n (and ret3 (mapcar ‘Ent:make_Point ret3))
\n (and ret4 (mapcar ‘Ent:make_Point ret4))
\n (foreach p ret2
\n (apply ‘Ent:make_line p)
\n )
\n )
\n )
\n)
\n;;;\u6d4b\u8bd5\u7a0b\u5e8f
\n(defun C:PPP (\/ pa pb pc pd ret)
\n (initget 1)
\n (setq pa (getpoint “n\u8f93\u5165\u7b2c\u4e00\u70b9:”))
\n (initget 1)
\n (setq pb (getpoint “n\u8f93\u5165\u7b2c\u4e8c\u70b9:”))
\n (initget 1)
\n (setq pc (getpoint “n\u8f93\u5165\u7b2c\u4e09\u70b9:”))
\n (initget 1)
\n (setq pd (getpoint “n\u8f93\u5165\u6620\u5c04\u70b9:”))<\/p>\n
(Ent:Make_Poly (list pa pb pc))
\n (setq InC (Tri:InCenter pa pb pc))
\n (setq CCC (CIR:PPP pa pb pc))
\n (setq Pd1 (Tri:Isogonal-Conjugate-Point Pd Pa Pb Pc))
\n (setq Otc (tri:OrthoCenter pa pb pc))
\n (setq 9pc (Tri:9P_Circle pa pb pc))
\n (mapcar ‘Ent:Make_Point (list pa pb pc pd Pd1 Otc))
\n (foreach n (list Inc CCC 9pc)
\n (apply ‘Ent:Make_Circle n)
\n )
\n (princ)
\n)<\/p>\n
;;;\u6d4b\u8bd5\u70b9\u70b9\u7ebf
\n(defun c:PPL(\/ p1 p2 ln dxf pa pb ret)
\n (setq p1 (getpoint “n\u7b2c\u4e00\u70b9:”))
\n (setq p2 (getpoint “n\u7b2c\u4e8c\u70b9:”))
\n (setq ln (car (entsel “n\u8bf7\u9009\u62e9\u7ebf\u6bb5”)))
\n (if (and p1 p2 ln)
\n (progn
\n (setq dxf (entget ln))
\n (setq ptA (cdr (assoc 10 dxf)))
\n (setq ptb (cdr (assoc 11 dxf)))
\n (setq ret (CIR:PPL p1 p2 PtA PtB))
\n (Ent:Make_Line p1 p2)
\n (if (VL-CONSP ret)
\n (foreach n ret
\n (and n (apply ‘Ent:Make_Circle n))
\n )
\n )
\n )
\n )
\n)<\/p>\n
;;;\u6d4b\u8bd5\u70b9\u7ebf\u7ebf\u5706
\n(defun C:PLL(\/ pt l1 l2 d1 d2 pa pb pc pd)
\n (setq pt (getpoint “n\u9009\u62e9\u70b9:”))
\n (setq l1 (car (entsel “n\u8bf7\u9009\u62e9\u7ebf\u6bb5”)))
\n (setq l2 (car (entsel “n\u8bf7\u9009\u62e9\u7ebf\u6bb5”)))
\n (if (and pt l1 l2)
\n (progn
\n (setq d1 (entget l1))
\n (setq d2 (entget l2))
\n (setq pa (cdr (assoc 10 d1)))
\n (setq pb (cdr (assoc 11 d1)))
\n (setq pc (cdr (assoc 10 d2)))
\n (setq pd (cdr (assoc 11 d2)))
\n (Ent:Make_Point pt)
\n (setq ret (CIR:PLL pt pa pb Pc Pd))
\n (if (VL-CONSP ret)
\n (foreach n ret
\n (and n (Ent:Make_Circle (car n) (cadr n)))
\n )
\n )
\n )
\n )
\n)<\/p>\n
;;;\u6d4b\u8bd5\u70b9\u5207\u5207\u5706
\n(defun C:PCC(\/ pt e1 e2 d1 d2 r1 r2 rt)
\n (initget 1)
\n (setq pt (getPoint “n\u70b91\uff1a”))
\n (setq e1 (ssget “:S” ‘((0 . “CIRCLE,ARC”))))
\n (setq e2 (ssget “:S” ‘((0 . “CIRCLE,ARC”))))
\n (if (and pt e1 e2)
\n (progn
\n (setq d1 (entget (ssname e1 0)))
\n (setq d2 (entget (ssname e2 0)))
\n (setq r1 (cdr (assoc 40 d1)))
\n (setq r2 (cdr (assoc 40 d2)))
\n (setq c1 (cdr (assoc 10 d1)))
\n (setq c2 (cdr (assoc 10 d2)))
\n (setq rt (CIR:PCC pt c1 r1 c2 r2))
\n (Ent:Make_Point pt)
\n (if rt
\n (foreach n rt
\n (apply ‘Ent:Make_Circle n)
\n )
\n )
\n )
\n )
\n)<\/p>\n
;;;\u6d4b\u8bd5\u70b9\u70b9\u7ebf\u5706
\n(defun C:PPC(\/ p1 p2 en dxf rad cen ret)
\n (initget 1)
\n (setq p1 (getPoint “n\u70b91\uff1a”))
\n (initget 1)
\n (setq p2 (getPoint “n\u70b92\uff1a”))
\n (setq en (car (entsel “n\u9009\u4e00\u4e2a\u5706\uff1a”)))
\n (if (and p1 p2 en)
\n (progn
\n (setq dxf (entget en))
\n (setq rad (cdr (assoc 40 dxf)))
\n (setq cen (cdr (assoc 10 dxf)))
\n (setq ret (CIR:PPC p1 p2 cen rad))
\n (Ent:Make_Point p1)
\n (Ent:Make_Point p2)
\n (if ret
\n (foreach n ret
\n (apply ‘Ent:Make_Circle n)
\n )
\n )
\n )
\n )
\n)<\/p>\n
;;;\u6d4b\u8bd5\u4e0e\u4e09\u5706\u90fd\u76f8\u5207\u7684\u5706(\u5373\u963f\u6ce2\u7f57\u5c3c\u65af\u5706)\u51fd\u6570
\n(defun c:rax (\/ e1 e2 e3 d1 d2 d3 c1 c2 c3 r1 r2 r3)
\n (setq e1 (car (entsel “n\u57061:”)))
\n (setq e2 (car (entsel “n\u57062:”)))
\n (setq e3 (car (entsel “n\u57063:”)))
\n (if (and e1 e2 e3)
\n (progn
\n (setq d1 (entget e1))
\n (setq d2 (entget e2))
\n (setq d3 (entget e3))
\n (setq r1 (cdr (assoc 40 d1))
\n r2 (cdr (assoc 40 d2))
\n r3 (cdr (assoc 40 d3))
\n c1 (cdr (assoc 10 d1))
\n c2 (cdr (assoc 10 d2))
\n c3 (cdr (assoc 10 d3))
\n )
\n (command “undo” “be”)
\n (vla-put-color (vlax-ename->vla-object e1) 1)
\n (vla-put-color (vlax-ename->vla-object e2) 2)
\n (vla-put-color (vlax-ename->vla-object e3) 3)
\n (MiSC:test
\n 100
\n ‘((CIR:CCC c1 r1 c2 r2 c3 r3)
\n (CIR:CCC_1 c1 r1 c2 r2 c3 r3)
\n )
\n )
\n (foreach p (CIR:CCC c1 r1 c2 r2 c3 r3)
\n (and p (apply ‘Ent:Make_Circle p))
\n )
\n (command “undo” “e”)
\n (princ)
\n )
\n )
\n)<\/p>\n
;;;\u6d4b\u8bd5\u70b9\u5207\u5207\u5706
\n(defun C:PLC(\/ pt e1 e2 d1 d2 r1 r2 ret)
\n (initget 1)
\n (setq pt0 (getPoint “nPoint:”))
\n (princ “nSelect a line:”)
\n (setq en1 (ssget “:S” ‘((0 . “LINE”))))
\n (princ “nSelect a circle:”)
\n (setq en2 (ssget “:S” ‘((0 . “CIRCLE”))))
\n (if (and pt0 en1 en2)
\n (progn
\n (setq d1 (entget (ssname en1 0)))
\n (setq d2 (entget (ssname en2 0)))
\n (setq ptA (cdr (assoc 10 d1)))
\n (setq ptB (cdr (assoc 11 d1)))
\n (setq rad (cdr (assoc 40 d2)))
\n (setq cen (cdr (assoc 10 d2)))
\n (Ent:Make_Point pt0)
\n (setq ret (CIR:PLC pt0 ptA ptB Cen Rad))
\n (foreach n ret
\n (apply ‘Ent:Make_Circle n)
\n )
\n )
\n )
\n)<\/p>\n
;;;\u6d4b\u8bd5\u6781\u70b9\u548c\u6781\u7ebf\u7684\u6837\u4f8b
\n(defun C:RAL(\/ e1 e2 d1 d1 c r f p1 p2 rt)
\n (setq e1 (car (entsel “n\u7ebf:”)))
\n (setq e2 (car (entsel “n\u5706:”)))
\n (if (and e1 e2)
\n (progn
\n (setq d1 (entget e1))
\n (setq d2 (entget e2))
\n (setq C (cdr (assoc 10 d2)))
\n (setq r (cdr (assoc 40 d2)))
\n (setq p1 (cdr (assoc 10 d1)))
\n (setq p2 (cdr (assoc 11 d1)))
\n (setq f (Line:Equation p1 p2))
\n (setq rt (CIR:Polar_Point c r (car f) (cadr f) (caddr f)))
\n (and rt (Ent:Make_Point rt))
\n )
\n )
\n)<\/p>\n
;;;\u6d4b\u8bd5\u4e09\u7ebf\u5207\u5706\u51fd\u6570
\n(defun c:LLL(\/ e1 e2 e3 d1 d2 d3 p1 p2 p3 p4 p5 p6 ret)
\n (prompt “n\u8bf7\u9009\u62e9\u4e09\u6761\u76f4\u7ebf:”)
\n (if (and (setq ss (ssget ‘((0 . “LINE”))))
\n (>= (sslength ss) 3)
\n )
\n (progn
\n (setq e1 (ssname ss 0))
\n (setq e2 (ssname ss 1))
\n (setq e3 (ssname ss 2))
\n (setq d1 (entget e1))
\n (setq d2 (entget e2))
\n (setq d3 (entget e3))
\n (setq p1 (cdr (assoc 10 d1)))
\n (setq p2 (cdr (assoc 11 d1)))
\n (setq p3 (cdr (assoc 10 d2)))
\n (setq p4 (cdr (assoc 11 d2)))
\n (setq p5 (cdr (assoc 10 d3)))
\n (setq p6 (cdr (assoc 11 d3)))
\n (setq ret (CIR:LLL p1 p2 p3 p4 p5 p6))
\n (foreach n ret
\n (apply ‘Ent:Make_Circle n)
\n )
\n )
\n )
\n)
\n;;;—————————————————-;
\n;;;\u7a0b\u5e8f\u7ed3\u5c3e ;
\n;;;————————————————–;|;
\n[\/codesyntax]<\/p>\n","protected":false},"excerpt":{"rendered":"
\u5bf9CAD\u4e2d\u5706\u548c\u5f27\u7684\u7814\u7a76\u3002\u5305\u62ec\u5982\u4e0b\u5185\u5bb9\uff1a \u5706\u548c\u5f27\u7684\u57fa\u672c\u8981\u7d20\uff0c\u5706\u7684\u548c\u5f27\u7684\u5468\u957f \uff0c\u9762\u79ef\uff0c\u8d28\u5fc3\uff1b \u5b83\u4eec\u7684\u79cd\u79cd\u4f5c\u56fe\u6cd5\uff0c\u8b6c<\/p>\n