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

aeo

最初由 四季豆 发布
[B]

我前面已经说过了。说详细一点吧，如图... [/B]

好象不完全.比如角度为-30 度,这方法就不行了.
我说要通用.

-30度应该是算60度的,第一点到了第四点上了.
还有0度?90度?

四季豆

写程序我不会丫。
不过，我认为知道了怎么做，写程序只不过是程序化的事情了。

cy956

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

2. 
   
3. cbx老师看来最近比较忙，没来。我来
   
4. 说说bx老师的看法：   :D
   
5. 
   
6. (defun c:qcr()
   
7.   (setq e0(entlast))
   
8.   (prompt "\n选择共圆的三点，如程序意外出错，表示三点共线.\n")
   
9.   (command "circle" "3p" pause pause pause)
   
10.   (setq e1(entlast))
    
11.   (if (/= e0 e1)(progn
    
12.     (setq en1(entget e1)
    
13.           r(dxf 40 en1)
    
14.           p1(dxf 10 en1))
    
15.     (entdel e1)
    
16.     (princ (strcat "\n三点圆半径＝" (rtos r 2 2) ",圆心为"))
    
17.     (princ p1)(princ)
    
18.     )
    
19.   )
    
20. )

<br /> cbx老师看来最近比较忙，没来。我来<br /> 说说bx老师的看法：   :D<br /> <br /> (defun c:qcr()<br />   (setq e0(entlast))<br />   (prompt "\n选择共圆的三点，如程序意外出错，表示三点共线.\n")<br />   (command "circle" "3p" pause pause pause)<br />   (setq e1(entlast))<br />   (if (/= e0 e1)(progn<br />     (setq en1(entget e1)<br />           r(dxf 40 en1)<br />           p1(dxf 10 en1))<br />     (entdel e1)<br />     (princ (strcat "\n三点圆半径＝" (rtos r 2 2) ",圆心为"))<br />     (princ p1)(princ)<br />     )<br />   )<br /> )

再看前面的帖子，发现龙龙仔居然和cbx老师一个思路。
老菜倒也。:)

cy956

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

2. 
   
3. (defun left(p1 p2 p3 / p c B C P)    ;;;p3点位于p1p2的关系
   
4.   (setq p p3)
   
5.   (apply '+ (mapcar '(lambda (b)
   
6.     (setq c(- (* (car p) (cadr b)) (* (cadr p) (car b))) p b )c )
   
7.     (list p1 p2 p3)))
   
8. )
   
9. (defun cy_left(p1 p2 p3 / k p c B C P);;;p3点位于p1p2的关系
   
10.   (setq p p3)
    
11.   (setq k(apply '+ (mapcar '(lambda (b)
    
12.     (setq c(- (* (car p) (cadr b)) (* (cadr p) (car b))) p b )c )
    
13.     (list p1 p2 p3))))
    
14.   (cond
    
15.     ((equal k 0.0 0.000001) 0.0);;上,此精度太小会出错
    
16.     ((> k 0) 1);;左
    
17.     (t -1);;右
    
18.   )
    
19. )
    

<br /> (defun left(p1 p2 p3 / p c B C P)    ;;;p3点位于p1p2的关系<br />   (setq p p3)<br />   (apply '+ (mapcar '(lambda (b)<br />     (setq c(- (* (car p) (cadr b)) (* (cadr p) (car b))) p b )c )<br />     (list p1 p2 p3)))<br /> )<br /> (defun cy_left(p1 p2 p3 / k p c B C P);;;p3点位于p1p2的关系<br />   (setq p p3)<br />   (setq k(apply '+ (mapcar '(lambda (b)<br />     (setq c(- (* (car p) (cadr b)) (* (cadr p) (car b))) p b )c )<br />     (list p1 p2 p3))))<br />   (cond <br />     ((equal k 0.0 0.000001) 0.0);;上,此精度太小会出错<br />     ((> k 0) 1);;左<br />     (t -1);;右<br />   )<br /> )<br />

aeo

用几何求很简单.
(inters p3(polar(xdrx_midp p1 p2)(+ 1.757(angle p1 p2))1e-9)
(polar p1(angle p1 p2)1e99) (polar p2(angle p2 p1)1e99)  )

相交p1 p2 p3逆时.