一些向量运算的子程序

陨落

话不多说，直接上源码

2. (defun num*vec(num vec);;;;;向量数乘
   
3.   (mapcar '(lambda(x) (* x num)) vec))
   
4. (defun vec+vec(vec1 vec2);;;;;向量相加
   
5.   (if (= (length vec1) (length vec2))
   
6.     (mapcar '(lambda (x y)(+ x y))vec1 vec2)
   
7.     (princ "\n出错：向量长度不一致")))
   
8. (defun vec-vec(vec1 vec2);;;;;向量相减
   
9.   (if (= (length vec1) (length vec2))
   
10.     (mapcar '(lambda (x y)(- x y))vec1 vec2)
    
11.     (princ "\n出错：向量长度不一致")))
    
12. 
    
13. (defun vec*vec(vec1 vec2);;;;;向量数量积
    
14.   (if (= (length vec1) (length vec2))
    
15.     (apply '+(mapcar '(lambda (x y)(* x y))vec1 vec2))
    
16.     (princ "\n出错：向量长度不一致")))
    
17. (defun vecmod(vec);;;;;;向量的模
    
18.   (sqrt(apply '+(mapcar '(lambda(x)(* x x)) vec))))
    
19. (defun pervec(pt1 pt2);;;;两点之间的单位向量
    
20.   (setq vecs(vec-vec pt2 pt1))
    
21.   (num*vec(/ 1(vecmod vecs))vecs)
    
22.   )

(defun num*vec(num vec);;;;;向量数乘<br />   (mapcar '(lambda(x) (* x num)) vec))<br /> (defun vec+vec(vec1 vec2);;;;;向量相加<br />   (if (= (length vec1) (length vec2))<br />     (mapcar '(lambda (x y)(+ x y))vec1 vec2)<br />     (princ "\n出错：向量长度不一致")))<br /> (defun vec-vec(vec1 vec2);;;;;向量相减<br />   (if (= (length vec1) (length vec2))<br />     (mapcar '(lambda (x y)(- x y))vec1 vec2)<br />     (princ "\n出错：向量长度不一致")))<br /> <br /> (defun vec*vec(vec1 vec2);;;;;向量数量积<br />   (if (= (length vec1) (length vec2))<br />     (apply '+(mapcar '(lambda (x y)(* x y))vec1 vec2))<br />     (princ "\n出错：向量长度不一致")))<br /> (defun vecmod(vec);;;;;;向量的模<br />   (sqrt(apply '+(mapcar '(lambda(x)(* x x)) vec))))<br /> (defun pervec(pt1 pt2);;;;两点之间的单位向量<br />   (setq vecs(vec-vec pt2 pt1))<br />   (num*vec(/ 1(vecmod vecs))vecs)<br />   )

kkq0305

(defun vec+vec(vec1 vec2);;;;;向量相加
  (if (= (length vec1) (length vec2))
    (mapcar '+ vec1 vec2)
    (princ "\n出错：向量长度不一致")))
(defun vec-vec(vec1 vec2);;;;;向量相减
  (if (= (length vec1) (length vec2))
    (mapcar '- vec1 vec2)
    (princ "\n出错：向量长度不一致")))

(defun vec*vec(vec1 vec2);;;;;向量数量积
  (if (= (length vec1) (length vec2))
    (apply '+(mapcar '* vec1 vec2))
    (princ "\n出错：向量长度不一致")))
(defun vecmod(vec);;;;;;向量的模
  (sqrt(apply '+(mapcar '* vec vec))))  可以这样写 不用lambda建立临时函数

ymcui

学习一下,谢谢   

ynhh

感谢大师的分享
不是这些用在什么地方啊
请指点

陨落

ynhh 发表于 2015-12-9 15:32
感谢大师的分享
不是这些用在什么地方啊
请指点

我也不是大师
在运算一些比较复杂的位置的时候，我感觉还是向量用得比较顺手

feng83

学习一下,谢谢楼主的源码

20060510412

请问如何实现向量旋转90度呢

lijunfa12345

学习学习，谢谢分享