[求助]：求三维中两条直线的夹角用什么函数

ssggg

谢谢
比较急

雨箭风刀

2. 
   
3. ;; a2l 求空间两直线夹角－－－ｂｙ　雨箭风刀.2006.5
   
4. ;; 前提：两直线必须有交点．
   
5. (defun a2l (/ e1 e2 p1 p2 p3 p4 px pu)
   
6.   (setq        e1 (car (entsel))
   
7.         e2 (car (entsel))
   
8.         p1 (vlax-curve-getstartpoint e1)
   
9.         p2 (vlax-curve-getendpoint e1)
   
10.         p3 (vlax-curve-getstartpoint e2)
    
11.         p4 (vlax-curve-getendpoint e2)
    
12.         px (cal "ill(p1,p2,p3,p4)")
    
13.         pa (cal "plt(p1,p1,0.5)")
    
14.         pb (cal "plt(p3,p4,0.5)")
    
15.         pu (cal "nor(p1,p2,pb)")
    
16.   )
    
17.   (if (not cal)(arxload "geomcal"))
    
18.   (princ(cal "ang(px,pa,pb,pu)"))
    
19. )
    
20. (vlax-add-cmd "a2l")
    
21. (vlax-add-cmd "a2l" 'a2l "a2l" 1)
    

<br /> ;; a2l 求空间两直线夹角－－－ｂｙ　雨箭风刀.2006.5<br /> ;; 前提：两直线必须有交点．<br /> (defun a2l (/ e1 e2 p1 p2 p3 p4 px pu)<br />   (setq        e1 (car (entsel))<br />         e2 (car (entsel))<br />         p1 (vlax-curve-getstartpoint e1)<br />         p2 (vlax-curve-getendpoint e1)<br />         p3 (vlax-curve-getstartpoint e2)<br />         p4 (vlax-curve-getendpoint e2)<br />         px (cal "ill(p1,p2,p3,p4)")<br />         pa (cal "plt(p1,p1,0.5)")<br />         pb (cal "plt(p3,p4,0.5)")<br />         pu (cal "nor(p1,p2,pb)")<br />   )<br />   (if (not cal)(arxload "geomcal"))<br />   (princ(cal "ang(px,pa,pb,pu)"))<br /> )<br /> (vlax-add-cmd "a2l")<br /> (vlax-add-cmd "a2l" 'a2l "a2l" 1)<br />

可透明使用

xxming

用CAL命令求得夹角
command:'cal
表达式：ang(int,end,end)
捕捉交角、和两个角点求得夹角

fsxm

向量夹角的余弦公式算出余弦值再反函数算出角就行了！

knightone

三维的很少用啊。抱歉了

fsxm

先用mapcar算出两个向量 A和B
cos@=|a1*a2+b1*b2+c1*c2|/(dist A)*(dist B)

winabcd

空间两直线不一定共面，其夹角是如何定义的？
我想半天没想出来，请教各位。

fsxm

上楼的老大~算偶不才 只是高在毕业
也知道空间两直线夹角是指平移至共面后所夹锐角啦!