三点求平面法线向量等计算几何代码

marting

1. ;向量叉积
   
2. (defun vcv (u v)
   
3.   (list
   
4.     (- (* (cadr u) (caddr v)) (* (cadr v) (caddr u)))
   
5.     (- (* (car v) (caddr u)) (* (car u) (caddr v)))
   
6.     (- (* (car u) (cadr v)) (* (car v) (cadr u)))
   
7.   )
   
8. )
   
9. ;向量乘标量
   
10. (defun vxs (v s)
    
11.   (mapcar '(lambda (n) (* n s)) v)
    
12. )
    
13. ;向量长度
    
14. (defun nrm (v)
    
15.   (sqrt (apply '+ (mapcar '(lambda (n) (* n n)) v)))
    
16. )
    
17. ;单位向量
    
18. (defun one (v)
    
19.   (vxs v (/ 1.0 (nrm v)))
    
20. )
    
21. ;法线向量
    
22. (defun normal (p1 p2 p3)
    
23.   (one (vcv (mapcar '- p3 p2) (mapcar '- p3 p1)))
    
24. )
    

;向量叉积<br /> (defun vcv (u v)<br />   (list<br />     (- (* (cadr u) (caddr v)) (* (cadr v) (caddr u)))<br />     (- (* (car v) (caddr u)) (* (car u) (caddr v)))<br />     (- (* (car u) (cadr v)) (* (car v) (cadr u)))<br />   )<br /> )<br /> ;向量乘标量<br /> (defun vxs (v s)<br />   (mapcar '(lambda (n) (* n s)) v)<br /> )<br /> ;向量长度<br /> (defun nrm (v)<br />   (sqrt (apply '+ (mapcar '(lambda (n) (* n n)) v)))<br /> )<br /> ;单位向量<br /> (defun one (v)<br />   (vxs v (/ 1.0 (nrm v)))<br /> )<br /> ;法线向量<br /> (defun normal (p1 p2 p3)<br />   (one (vcv (mapcar '- p3 p2) (mapcar '- p3 p1)))<br /> )<br />

ynhh

感谢大师
不知能不能
有个实例学习一下具体应用？

aeo

小心点重合

newer

ynhh 发表于 2017-4-1 10:47

上面那个normal函数，就是上面所有函数的应用，你给三个点，就求出三点组成平面的法线向量了。

kqqt6236

啥也不说了，感谢楼主分享哇！

sicky111

啥也不说了，感谢楼主分享哇！