[求助]：过两点并与另一条线相切的圆

xyp1964 · 发表于 2005-5-29 17:10:21

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？立即注册

×

过两点并与另一条线相切的圆”，对这个类似几何的题，征求解法！

狂刀 · 发表于 2005-5-29 18:04:18

(command "circle" "3p" "end" pause "end" pause "tan" pause)

xyp1964 · 发表于 2005-5-29 18:09:54

简单明了，不负“狂刀”的大名！多谢！
如果用几何方法怎么解决呢？有点钻牛角了，呵呵……

xiao_longxin · 发表于 2005-6-1 00:56:10

纯几何方法解决：

;;By Xiao_longxin
;;po_star 直线起点 po_end　直线终点
;;p1、p2　所要求的圆上两点
(defun errexit (s)
(command "._undo" "_E")
(princ "\nError: ")
(princ s)
)
(defun twopo_line_qieyuan (po_star po_end p1 p2
/ cr_list jd1　jd2 ang_zxcz
ang1 ang2 ang3 ang4
dist1 dist2 p1_zx p2_zx
p1p2 cr_pt cr_r ent
oldcmdecho　oldosmode bl_gd
)
;;(setq po_star (getpoint))
;;(setq po_end (getpoint))
;;(setq p1 (getpoint))
;;(setq p2 (getpoint))
;;(command "point" po_star)
;;(command "point" po_end)
;;(command "point" p1)
;;(command "point" p2)
;;(command "pline" po_star po_end "")
(setq oldcmdecho (getvar "cmdecho"))
(setvar "cmdecho" 0)
(setq oldosmode (getvar "osmode"))
(setvar "osmode" 0)
(setq ang1 (angle po_star po_end)) ;直线的方位角
(setq ang2 (angle p1 p2)) ;二点的方位角
(setq ang_zxcz (rem (+ ang1 (/ pi 2)) (* 2 pi))) ;与直线垂直的方位角
;;求p1点与直线的垂距，以此判断是否在直线上
(setq jd1 (polar p1 ang_zxcz 100))
(setq jd1 (inters po_star po_end p1 jd1 nil))
(setq dist1 (distance p1 jd1))
(setq ang3 (angle jd1 p1))
(if (< dist1 0.001)
(setq p1_zx t)
)
;;求p2点与直线的垂距，以此判断是否在直线上
(setq jd2 (polar p2 ang_zxcz 100))
(setq jd2 (inters po_star po_end p2 jd2 nil))
(setq dist2 (distance p2 jd2))
(setq ang4 (angle jd2 p2))
(if (< dist2 0.001)
(setq p2_zx t)
)
;判断直线是否在两点间
(if (and (> (abs (- ang3 ang4)) 0.00000001)
(and (not p2_zx) (not p1_zx))
)
(setq p1p2 t)
)
(if (or p1p2 (and p2_zx p1_zx))
(progn ;如两点都在线上则中止程序
(prompt "\n两点在直线上或者点在直线两边，没有切圆")
(setvar "cmdecho" oldcmdecho)
(setvar "osmode" oldosmode)
(setq pr_list nil)
)
(progn
(cond
((or (< (abs (- ang1 ang2)) 0.000005)
(and (< (abs (- ang1 ang2)) (+ pi 0.000005))
(> (abs (- ang1 ang2)) (- pi 0.000005))
)
) ;p1、p2点的连线与直线平行
(setq dist1 (/ (distance p1 p2) 2))
(setq jd1 (polar p1 ang2 dist1))
(setq jd2 (polar jd1 ang_zxcz 100))
(setq jd1 (inters jd1 jd2 po_star po_end nil))
(command "._circle" "3p" p1 p2 jd1)
(setq ent (entget (entlast)))
(setq cr_pt (cdr (assoc 10 ent)))
(setq cr_r (cdr (assoc 40 ent)))
(setvar "cmdecho" oldcmdecho)
(setvar "osmode" oldosmode)
(setq cr_list (list cr_pt cr_r))
) ;end p1、p2点的连线与直线平行
((or p2_zx p1_zx) ;p1、p2点有一点在直线上
(if p1_zx
(progn
(setq jd1 p1)
(setq jd2 p2)
)
(progn
(setq jd1 p2)
(setq jd2 p1)
)
)
(setq p1 (polar jd1 ang1 100))
(if (= 100 (setq dist2 (distance jd2 p1)))
(progn
(setq p1 (polar jd1 ang1 200))
(setq dist2 (distance p1 jd2))
)
)
(setq dist1 (distance p1 jd1))
(setq dist1 (/ (* dist1 dist1) dist2))
(setq ang3 (angle p1 jd2))
(setq p1 (polar p1 ang3 dist1))
(command "._circle" "3p" jd1 jd2 p1)
(setq ent (entget (entlast)))
(setq cr_pt (cdr (assoc 10 ent)))
(setq cr_r (cdr (assoc 40 ent)))
(setvar "cmdecho" oldcmdecho)
(setvar "osmode" oldosmode)
(setq cr_list (list cr_pt cr_r))
) ;end p1、p2点有一点在直线上
(t ;一般情况
(setq jd1 (inters po_star po_end p1 p2 nil))
(setq dist1 (distance p1 jd1))
(setq dist2 (distance p2 jd1))
(setq dist1 (sqrt (* dist1 dist2)))
(if (or (< (abs (- ang_zxcz ang2)) 0.000000005)
(and (< (abs (- ang_zxcz ang2)) (+ pi 0.000000005))
(> (abs (- ang_zxcz ang2)) (- pi 0.000000005))
)
)
(progn
(setq ang1 (angle jd1 po_star))
(setq jd2 (polar jd1 ang1 dist1))
(command "._circle" "3p" jd2 p2 p1)
(setq ent (entget (entlast)))
(setq cr_pt (cdr (assoc 10 ent)))
(setq cr_r (cdr (assoc 40 ent)))
(setq cr_list (list cr_pt cr_r))
(setq ang1 (+ ang1 pi))
(setq jd2 (polar jd1 ang1 dist1))
(command "._circle" "3p" jd2 p2 p1)
(setq ent (entget (entlast)))
(setq cr_pt (cdr (assoc 10 ent)))
(setq cr_r (cdr (assoc 40 ent)))
(setvar "cmdecho" oldcmdecho)
(setvar "osmode" oldosmode)
(setq cr_list (append cr_list (list cr_pt cr_r)))
)
(progn
(setq ang1 (angle jd1 jd2))
(setq jd2 (polar jd1 ang1 dist1))
(command "._circle" "3p" jd2 p2 p1)
(setq ent (entget (entlast)))
(setq cr_pt (cdr (assoc 10 ent)))
(setq cr_r (cdr (assoc 40 ent)))
(setvar "cmdecho" oldcmdecho)
(setvar "osmode" oldosmode)
(setq cr_list (list cr_pt cr_r))
)
) ;end if
) ;end 一般情况
)
)
)
)

;;By Xiao_longxin 
;;po_star 直线起点           po_end　直线终点 
;;p1、p2　所要求的圆上两点 
 
 
(defun errexit (s) 
  (command "._undo" "_E") 
  (princ "\nError:  ") 
  (princ s) 
) 
 
(defun twopo_line_qieyuan (po_star    po_end         p1            p2 
                           /              cr_list         jd1　jd2   ang_zxcz 
                           ang1              ang2         ang3            ang4 
                           dist1      dist2         p1_zx            p2_zx 
                           p1p2              cr_pt         cr_r            ent 
                           oldcmdecho　oldosmode bl_gd 
                          ) 
 
  ;;(setq po_star (getpoint)) 
  ;;(setq po_end (getpoint)) 
  ;;(setq p1 (getpoint)) 
  ;;(setq p2 (getpoint)) 
  ;;(command "point" po_star) 
  ;;(command "point" po_end) 
  ;;(command "point" p1) 
  ;;(command "point" p2) 
  ;;(command "pline" po_star po_end "") 
 
  (setq oldcmdecho (getvar "cmdecho")) 
  (setvar "cmdecho" 0) 
  (setq oldosmode (getvar "osmode")) 
  (setvar "osmode" 0) 
 
  (setq ang1 (angle po_star po_end))        ;直线的方位角 
  (setq ang2 (angle p1 p2))                ;二点的方位角 
  (setq ang_zxcz (rem (+ ang1 (/ pi 2)) (* 2 pi))) ;与直线垂直的方位角 
 
  ;;求p1点与直线的垂距，以此判断是否在直线上 
  (setq jd1 (polar p1 ang_zxcz 100)) 
  (setq jd1 (inters po_star po_end p1 jd1 nil)) 
  (setq dist1 (distance p1 jd1)) 
  (setq ang3 (angle jd1 p1)) 
  (if (< dist1 0.001) 
    (setq p1_zx t) 
  ) 
  ;;求p2点与直线的垂距，以此判断是否在直线上 
  (setq jd2 (polar p2 ang_zxcz 100)) 
  (setq jd2 (inters po_star po_end p2 jd2 nil)) 
  (setq dist2 (distance p2 jd2)) 
  (setq ang4 (angle jd2 p2)) 
  (if (< dist2 0.001) 
    (setq p2_zx t) 
  ) 
                                        ;判断直线是否在两点间 
  (if (and (> (abs (- ang3 ang4)) 0.00000001) 
           (and (not p2_zx) (not p1_zx)) 
      ) 
    (setq p1p2 t) 
  ) 
 
  (if (or p1p2 (and p2_zx p1_zx)) 
    (progn                                ;如两点都在线上则中止程序 
      (prompt "\n两点在直线上或者点在直线两边，没有切圆") 
      (setvar "cmdecho" oldcmdecho) 
      (setvar "osmode" oldosmode) 
      (setq pr_list nil) 
    ) 
    (progn 
      (cond 
        ((or (< (abs (- ang1 ang2)) 0.000005) 
             (and (< (abs (- ang1 ang2)) (+ pi 0.000005)) 
                  (> (abs (- ang1 ang2)) (- pi 0.000005)) 
             ) 
         )                                ;p1、p2点的连线与直线平行 
         (setq dist1 (/ (distance p1 p2) 2)) 
         (setq jd1 (polar p1 ang2 dist1)) 
         (setq jd2 (polar jd1 ang_zxcz 100)) 
         (setq jd1 (inters jd1 jd2 po_star po_end nil)) 
         (command "._circle" "3p" p1 p2 jd1) 
         (setq ent (entget (entlast))) 
         (setq cr_pt (cdr (assoc 10 ent))) 
         (setq cr_r (cdr (assoc 40 ent))) 
         (setvar "cmdecho" oldcmdecho) 
         (setvar "osmode" oldosmode) 
         (setq cr_list (list cr_pt cr_r)) 
        )                                ;end   p1、p2点的连线与直线平行 
        ((or p2_zx p1_zx)                ;p1、p2点有一点在直线上 
         (if p1_zx 
           (progn 
             (setq jd1 p1) 
             (setq jd2 p2) 
           ) 
           (progn 
             (setq jd1 p2) 
             (setq jd2 p1) 
           ) 
         ) 
         (setq p1 (polar jd1 ang1 100)) 
         (if (= 100 (setq dist2 (distance jd2 p1))) 
           (progn 
             (setq p1 (polar jd1 ang1 200)) 
             (setq dist2 (distance p1 jd2)) 
           ) 
         ) 
         (setq dist1 (distance p1 jd1)) 
         (setq dist1 (/ (* dist1 dist1) dist2)) 
         (setq ang3 (angle p1 jd2)) 
         (setq p1 (polar p1 ang3 dist1)) 
         (command "._circle" "3p" jd1 jd2 p1) 
         (setq ent (entget (entlast))) 
         (setq cr_pt (cdr (assoc 10 ent))) 
         (setq cr_r (cdr (assoc 40 ent))) 
         (setvar "cmdecho" oldcmdecho) 
         (setvar "osmode" oldosmode) 
         (setq cr_list (list cr_pt cr_r)) 
        )                                ;end  p1、p2点有一点在直线上 
        (t                                ;一般情况 
         (setq jd1 (inters po_star po_end p1 p2 nil)) 
         (setq dist1 (distance p1 jd1)) 
         (setq dist2 (distance p2 jd1)) 
         (setq dist1 (sqrt (* dist1 dist2))) 
         (if (or (< (abs (- ang_zxcz ang2)) 0.000000005) 
                 (and (< (abs (- ang_zxcz ang2)) (+ pi 0.000000005)) 
                      (> (abs (- ang_zxcz ang2)) (- pi 0.000000005)) 
                 ) 
             ) 
           (progn 
             (setq ang1 (angle jd1 po_star)) 
             (setq jd2 (polar jd1 ang1 dist1)) 
             (command "._circle" "3p" jd2 p2 p1) 
             (setq ent (entget (entlast))) 
             (setq cr_pt (cdr (assoc 10 ent))) 
             (setq cr_r (cdr (assoc 40 ent))) 
             (setq cr_list (list cr_pt cr_r)) 
             (setq ang1 (+ ang1 pi)) 
             (setq jd2 (polar jd1 ang1 dist1)) 
             (command "._circle" "3p" jd2 p2 p1) 
             (setq ent (entget (entlast))) 
             (setq cr_pt (cdr (assoc 10 ent))) 
             (setq cr_r (cdr (assoc 40 ent))) 
             (setvar "cmdecho" oldcmdecho) 
             (setvar "osmode" oldosmode) 
             (setq cr_list (append cr_list (list cr_pt cr_r))) 
           ) 
           (progn 
             (setq ang1 (angle jd1 jd2)) 
             (setq jd2 (polar jd1 ang1 dist1)) 
             (command "._circle" "3p" jd2 p2 p1) 
             (setq ent (entget (entlast))) 
             (setq cr_pt (cdr (assoc 10 ent))) 
             (setq cr_r (cdr (assoc 40 ent))) 
             (setvar "cmdecho" oldcmdecho) 
             (setvar "osmode" oldosmode) 
             (setq cr_list (list cr_pt cr_r)) 
           ) 
         )                                ;end if 
        )                                ;end 一般情况 
      ) 
    ) 
  ) 
)

xyp1964 · 发表于 2005-6-1 01:32:38

利用以上程序，可行：
[php]
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
(defun c:test ()
  (setq pt1 (getpoint "\n线端点1 : ")
pt2 (getpoint "\n线端点2 : ")
pt3 (getpoint "\n点1 : ")
pt4 (getpoint "\n点2 : ")
  )
  (twopo_line_qieyuan pt1 pt2 pt3 pt4)
  (princ)
)[/php]

xiao_longxin · 发表于 2005-6-1 09:11:39

请问版主，我如何把代码帖成你你这样子的php代码？？

账号		自动登录	找回密码
密码			立即注册

[求助] [求助]：过两点并与另一条线相切的圆

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

相关帖子