[求助]：已知三点求外接圆的圆心和半径

aeo · 发表于 2003-4-15 23:05:25

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

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

×

有没有很简洁的写法.
帮我写一下吧.

XDSoft · 发表于 2003-4-15 23:45:19

最初由 aeo 发布
[B]有没有很简洁的写法.
帮我写一下吧. [/B]

任意两条线中垂线的交点不就是圆心吗? 圆心和任意一个点的距离就是半径了.

aeo · 发表于 2003-4-16 00:07:13

我用角度求,求死我了.
我刚也想到了.

有了圆心,起始角和终止角也有了.
当然,用entmake生成这arc,还要判断时针.

第二个问题,已知一长方体的和X轴角度和它的boundingbox,求这长方体的四个点.(唯一吗?)

cy956 · 发表于 2003-4-16 02:13:15

最初由 aeo 发布
[B]我用角度求,求死我了.
我刚也想到了.

有了圆心,起始角和终止角也有了.
当然,用entmake生成这arc,还要判断时针.

第二个问题,已知一长方体的和X轴角度和它的boundingbox,求这长方体的四个点.(唯一吗?) [/B]

你说是不是唯一？

lxfbj · 发表于 2003-4-16 08:50:20

这种方法好

LUCAS · 发表于 2003-4-16 11:59:46

3點畫弧再求圓心,半徑

四季豆 · 发表于 2003-4-19 00:23:53

应该是唯一解。
CY画的那个是特例，外包是正方形，里面的矩形倾角为45度，这种情况是有无数解的。（如果倾角是45度，外包肯定是个正方形）

求法也很简单，一个一次方程而已。

eachy · 发表于 2003-4-19 01:12:58

最初由 XDSoft 发布
[B]

任意两条线中垂线的交点不就是圆心吗? 圆心和任意一个点的距离就是半径了. [/B]

(defun 3ptocircle (p1 p2 p3 / midp midp1 midp11 midp2 midp21 pcen)
;;点是否在LINE上
(Defun ponln (p p1 p2 / l)
(foreach x (list p1 p2 p)
(setq l (cons (distance p x) l)
p x
)
)
(if (equal (+ (car l) (caddr l)) (cadr l))
'T
nil
)
)
;;求中点
(defun midp (p1 p2)
(mapcar '(lambda (x) (/ x 2)) (mapcar '+ p1 p2))
)
(if (not (ponln p1 p2 p3));不共线
(progn
(setq midp1 (midp p1 p2)
midp2 (midp p2 p3)
angle1 (angle p1 p2)
angle2 (angle p2 p3)
midp11 (polar midp1 (+ angle1 (/ pi 2)) 1.)
midp21 (polar midp2 (+ angle2 (/ pi 2)) 1.)
pcen (inters midp1 midp11 midp2 midp21 nil)
)
(list pcen (distance pcen p1))
)
)
)

lijiao · 发表于 2003-4-19 09:30:01

我认为LUCAS 的方法最好，既简单，又实用。

eachy · 发表于 2003-4-19 12:08:33

最初由 lijiao 发布
[B]我认为LUCAS 的方法最好，既简单，又实用。 [/B]

画园也要判断是否共线，如果共线就不会生成园。

aeo · 发表于 2003-4-19 17:46:07

我要的是第二个问题!
能不能写个通用的.已知 box '(p1 p2 p3 p4 ) 和ang.

求圆心只要想起来就简单了,我开始用代数求.

四季豆 · 发表于 2003-4-26 12:38:54

最初由 aeo 发布
[B]我要的是第二个问题!
能不能写个通用的.已知 box '(p1 p2 p3 p4 ) 和ang.

求圆心只要想起来就简单了,我开始用代数求. [/B]

我前面已经说过了。说详细一点吧，如图。解一个二元一次方程组即可。

LUCAS · 发表于 2003-4-26 16:23:38

;;;判断是否共线
;;;Are Points CoLinear
;;;ARG -> 3 Points and Fuzz
;;;RET -> T nil
;;; BY David Bethel
(defun IS_PT_COLINEAR (P1 P2 P3 FZ)
(or
(equal
(distance P1 P3)
(+
(distance P2 P1)
(distance P2 P3)
)
FZ
)
(equal
(distance P1 P2)
(+
(distance P3 P1)
(distance P3 P2)
)
FZ
)
(equal
(distance P2 P3)
(+
(distance P1 P2)
(distance P1 P3)
)
FZ
)
)
)

梦断江南 · 发表于 2003-4-26 19:00:47

最初由 eachy 发布
[B]
画园也要判断是否共线，如果共线就不会生成园。 [/B]

另外:
;;点是否在LINE上
  (Defun ponln (p p1 p2 / l)
(foreach x (list p1 p2 p)
   (setq l (cons (distance p x) l)
      p x
   )
)
(if (equal (+ (car l) (caddr l)) (cadr l))
   'T
   nil
)
  )
这个判断方法似有问题.

[FONT=courier new]
;;已知3点求外接圆圆心半径.--------lxx2003.4
(defun 3pcr ( p1 p2 p3 / oe ne ent ptcen rad)
(setq oerr *error*
oe (entlast))
(defun *error* (msg) (princ "\n3点无法画圆,请检查是否共线或有重复点")(setq *error* oerr)(princ))
(command ".circle" "3p" p1 p2 p3)
(if (entnext oe)
(progn
(setq entl (entget (entlast))
ptcen (cdr (assoc 10 entl))
rad (cdr (assoc 40 entl))
)
(princ "圆心=")(princ ptcen)
(princ " 半径=")(princ rad)
(entdel (entlast))
)
)(princ)
)
[/FONT]

再探讨一下判断共线的方法
1.共线则总有其中某两点距离等于两两点位距离和的一半.

[FONT=courier new]
;;点是否在LINE上
(Defun ponln (p p1 p2 / l all key)
(setq l (list (distance p p1) (distance p1 p2) (distance p2 p))
all (+(car l) (cadr l)(caddr l)
)
(foreach x l
(if (equal x (/ all 2)) (setq key 'T))
)
)
[/FONT]

2.共线某点与其它两点的角度之差为0,或者180度*n
3.共线则某两点距离等于其它两点两两间距离之和(参考上贴)

[FONT=courier new]
(setq d1 (distance p p1)
d2 (distance p1 p2)
d3 (distance p2 p)
)
(if (or (= d1 (+ d2 d3)) (= d2 (+ d1 d3))(= d3 (+ d1 d2)))
(setq key 'T)
)
[/FONT]

4.共线则用inters求交点为nil

[FONT=courier new]
(if (and (not (equal p1 p2)) (not (equal p1 p3)));重点时inters不为nil
(if (inters p1 p3 p1 p2) (setq key 'T))
)
[/FONT]

aeo · 发表于 2003-4-27 10:44:01

最初由梦断江南发布
[B]
另外:
;;点是否在LINE上
(Defun ponln (p p1 p2 / l)
(foreach x (list p1 p2 p)
(setq l (c... [/B]

还有种方法
p3到p1 p2连线距离为0.

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

[求助] [求助]：已知三点求外接圆的圆心和半径

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

相关帖子

Re: [求助]：已知三点求外接圆的圆心和半径

Re: Re: [求助]：已知三点求外接圆的圆心和半径