Autolisp函数求救

SQ215

己知一个圆弧的两个端点及圆弧的半径,怎么设计一个Autolisp函数返回这个圆弧圆心？求求大家给我想想办法，谢了！

winabcd

解一个二元二次方程组就可以了，一般情况下有两个解，当然也可能只有一个解或无解。

 

Andyhon

(defun acos (c b)
    (atan (sqrt (- (* c c) (* b b))) b)
 )

 ;若有二解:
 (setq vv (angle p1 p2)
       aa (acos (* rr 2) (distance p1 p2))
      po1 (polar p1 (+ vv aa) rr)
      po2 (polar p1 (- vv aa) rr)
 )

ZZXXQQ

2. ;;已知圆弧的两个端点及圆弧的半径，求圆弧的圆心。
3. (DEFUN TEST (P1 P2 R / D ANG DF PC1 PC2)
4. (SETQ D (/ (DISTANCE P1 P2) 2))
5. (IF (< D R) (PROGN
6.   (SETQ ANG (ANGLE P1 P2)
7.         DF (SQRT (- (* R R) (* D D)))
8.         PC1 (POLAR (POLAR P1 ANG D) (+ (/ PI 2) ANG) DF)
9.         PC2 (POLAR PC1 (+ (/ PI 2) ANG) R))
10.   (LIST PC1 PC2)
11. )
12.   (PRINC "\nData Error!")
13. )
14. )

复制代码

随便逛一逛

2点: p1 p2

半径:R

(COMMAND "ARc" p1 "e" p2 "R" r)

(entlast)

-->求圆心?

SQ215

半圆形圆弧的起点160.4, 1.0终点 162.401, 1.0半径1.0004

可得到的圆心却是:161.4, 1.0    这个结果是圆心到两端点的距离不一样?

这样的问题应怎样解决?

ZZXXQQ

有点错，改了一下。

2. (DEFUN TEST (P1 P2 R / D ANG DF PC1 PC2)
3.   (SETQ D (/ (DISTANCE P1 P2) 2))
4.   (IF (< D R) (PROGN
5.    (SETQ ANG (ANGLE P1 P2)
6.          DF (SQRT (- (* R R) (* D D)))
7.          PC1 (POLAR (POLAR P1 ANG D) (+ (/ PI 2) ANG) DF)
8.          PC2 (POLAR PC1 (- ANG (/ PI 2)) (+ DF DF)))
9.    (LIST PC1 PC2)
10.   )
11.    (PRINC "\nData Error!")
12.   )
13. )

复制代码

Andyhon

半径1.0004 ==> 直径 2.0008
两点距 2.0010

直径 < 弦长 ==> 无解

 

SQ215

是半圆,应该两端点的距离=直径,而现在却出现了半圆两端点的距离NO=直径这种现象,在实际应用当中出现了这种误差应如何处理?而输出某些其它软件的文件时,却出现因为这种原因而造成其它软件无法接受.这个小问题害苦我了~~

Andyhon

 (if (equal 两点距 直径 可接受的公差)
  (setq po (mapcar
             (function (lambda (a b) (/ (+ a b) 2)))
             p1
             p2
  )        )
)