请问用lisp怎么求两个线段的交角

pros

谢谢

ZZXXQQ

程序：

2. ;求两直线的最小夹角 ZZXXQQ 2005.11.2
3. (DEFUN C:2LA (/ GETLINE S S2 EN1 EN2 P1 P2 P3 P4 P5 PT1 PT2 AN1 AN2 ANG)
4. (DEFUN GETLINE (MSG / A1)
5.   (INITGET 1)
6.   (SETQ A1 (CAR (ENTSEL MSG)))
7.   (WHILE (/= (CDR (ASSOC 0 (ENTGET A1))) "LINE")
8.    (PRINC "\n您选的不是线图元,请再选一次...")
9.    (INITGET 1)
10.    (SETQ A1 (CAR (ENTSEL MSG)))
11.   )
12.   A1
13. )
14. (SETVAR "CMDECHO" 0)
15. (SETQ S1 (GETLINE "\nSelect first Line :")
16.        S2 (GETLINE "\nSelect second Line :")
17.        EN1 (ENTGET S1) P1 (CDR (ASSOC 10 EN1)) P2 (CDR (ASSOC 11 EN1))
18.        EN2 (ENTGET S2) P3 (CDR (ASSOC 10 EN2)) P4 (CDR (ASSOC 11 EN2)))
19. (IF (SETQ P5 (INTERS P1 P2 P3 P4 nil)) (PROGN
20.   (SETQ PT1 (IF (< (DISTANCE P1 P5) (DISTANCE P2 P5)) P2 P1)
21.         PT2 (IF (< (DISTANCE P3 P5) (DISTANCE P4 P5)) P4 P3)
22.         AN1 (ANGLE P5 PT1) AN2 (ANGLE P5 PT2)
23.         AN1 (IF (> AN1 PI) (- AN1 PI PI) AN1)
24.         AN2 (IF (> AN2 PI) (- AN2 PI PI) AN2)
25.         ANG (- (MAX AN1 AN2) (MIN AN1 AN2)))
26.   (PRINC "\nLine angle is : ")
27.   (PRINC (* (/ ANG PI) 180))
28.   (PRINC " DEG")
29. ))
30. (SETVAR "CMDECHO" 1)
31. (PRINC)
32. )

复制代码

xyp1964

2. (load "xyp_lib")
3. ;|加载通用函数(可在签名栏直接下载)
4. 如果已经下载xyp_lib并放到搜索路径下可以不再下载!
5. 利用以下任何一种方式（首选第一种）即可加载和运行通用函数内的所有子程序：
6. ★1·在acad.lsp中增加(load"xyp_lib")
7. ■2·在每个程序内增加(load"xyp_lib")
8. ■3·在command下，输入(load"xyp_lib")
9. ■4·在菜单.mnl中增加(load"xyp_lib")
10. ■5·将xyp_lib.vlx文件直接拽到cad屏幕
11. [COLOR=red] ★通用函数下载地址：[/COLOR]
12. [url]http://www.xdcad.net/forum/attachment.php?s=&postid=1606661[/url]
13. |;
14. ;;;求两直线的最小夹角
15. (defun c:test ()
16.   (setq s1 (usel1 0 "LINE""LINE线1")
17. s2 (usel1 0 "LINE""LINE线2")
18. ang1(xyp-get-AngleAtPoint (car s1) (osnap (cadr s1) "near"))
19. ang2(xyp-get-AngleAtPoint (car s2) (osnap (cadr s2) "near"))
20.   )
21.   (princ"\n两直线最小夹角 = ")
22.   (princ (abs(-(rem (abs(rad2ang(- ang1 ang2)))90.0)90)))
23.   (princ)
24. )

(load "xyp_lib") ;|加载通用函数(可在签名栏直接下载) 如果已经下载xyp_lib并放到搜索路径下可以不再下载! 利用以下任何一种方式（首选第一种）即可加载和运行通用函数内的所有子程序： ★1·在acad.lsp中增加(load"xyp_lib") ■2·在每个程序内增加(load"xyp_lib") ■3·在command下，输入(load"xyp_lib") ■4·在菜单.mnl中增加(load"xyp_lib") ■5·将xyp_lib.vlx文件直接拽到cad屏幕 [COLOR=red] ★通用函数下载地址：[/COLOR] [url]http://www.xdcad.net/forum/attachment.php?s=&postid=1606661[/url] |; ;;;求两直线的最小夹角 (defun c:test ()   (setq s1 (usel1 0 "LINE""LINE线1") s2 (usel1 0 "LINE""LINE线2") ang1(xyp-get-AngleAtPoint (car s1) (osnap (cadr s1) "near")) ang2(xyp-get-AngleAtPoint (car s2) (osnap (cadr s2) "near"))   )   (princ"\n两直线最小夹角 = ")   (princ (abs(-(rem (abs(rad2ang(- ang1 ang2)))90.0)90)))   (princ) )

无痕

 (princ (abs(-(rem (abs(rad2ang(- ang1 ang2)))90.0)90)))

是不是有点问题？