【越飞越高讲堂15】用LISP论矩阵

cheng5276

高人！明经的灵魂之一

004

谢谢高飞老师.

jishiman

留个位置。。线性代数都忘的差不多了。。

tm20038175

太帅了，帅呆了！…

waterchen

先收藏，再学习，感叹LZ真高人

highflybir

zdqwy19 发表于 2013-1-5 20:00
矩阵求逆的网页没有帐号打不开，又无法注册。能否把代码贴出来。

这里是那几个链接的代码。

zdqwy19

highflybir 发表于 2013-1-9 12:49
这里是那几个链接的代码。

谢谢了，高兄！

springwillow

纯支持一下吧，要想玩这个还得补一下线性代数的知识，忘光光了!

LLXXZZ

必须下载学习...........

highflybir

再来举一个例子：
3D 的环形阵列。
选择物体，然后对物体集进行3d环形阵列。
很类似命令:3darray.
经过测试，比3darray快了很多，而且不出错。

1. ;;;3d环形阵列测试
   
2. (defun c:PAT(/ ss N P1 P2 ActDoc ObjLst)                   ;c:PolarArrayTest
   
3.   (setq ss (ssget))
   
4.   (initget 7)
   
5.   (setq N (getInt "\n数量:"))
   
6.   (initget 9)
   
7.   (setq P1 (getpoint "\n中心点:"))
   
8.   (initget 9)
   
9.   (setq p2 (getpoint P1 "\n另一点:"))
   
10.   (if ss
    
11.     (progn
    
12.       (setq ActDoc (vla-get-ActiveDocument (vlax-get-acad-object)))
    
13.       (vla-StartUndoMark ActDoc)
    
14.       (setq Ojblst (Ent:Ents->Objs ss))
    
15.       (Mat:3dPolarArray ojblst N (+ pi pi) T (trans p1 1 0) (trans p2 1 0))
    
16.       (vla-EndUndoMark ActDoc)
    
17.       (vlax-release-object ActDoc)
    
18.       (princ)
    
19.     )
    
20.   )
    
21. )
    
22. 
    
23. 
    
24. ;;;-----------------------------------------------------------;;
    
25. ;;; 三维环形阵列                                              ;;
    
26. ;;; 输入：Objlst    -- 物体集                                 ;;
    
27. ;;;       Number    -- 要阵列的个数（包含自身在内）           ;;
    
28. ;;;       FillAngle -- 旋转角度                               ;;
    
29. ;;;       IsCCW     -- 是否逆时针                             ;;
    
30. ;;;       P1        -- 阵列中心点                             ;;
    
31. ;;;       P2        -- 阵列轴线的另一点                       ;;
    
32. ;;; 输出：阵列的物体列表                                      ;;
    
33. ;;;-----------------------------------------------------------;;
    
34. (defun Mat:3dPolarArray (Objlst Number FillAngle IsCCW P1 P2 / lst1 lst2 ANG MAT NEW)
    
35.   (if IsCCW
    
36.     (setq FillAngle (float FillAngle))
    
37.     (setq FillAngle (- FillAngle (+ pi pi)))
    
38.   )
    
39.   (setq ang (/ FillAngle (1- Number)))
    
40.   (setq mat (MAT:RotateBy2P P1 P2 ang))
    
41.   (setq mat (vlax-tmatrix mat))
    
42.   (repeat (1- Number)
    
43.     (setq lst1 nil)
    
44.     (foreach obj ObjLst
    
45.       (setq new (vla-copy obj))
    
46.       (vla-transformby new mat)
    
47.       (setq lst1 (cons new lst1))
    
48.     )
    
49.     (setq objLst (reverse lst1))
    
50.     (setq lst2 (cons objLst lst2))
    
51.   )
    
52.   (reverse lst2)
    
53. )
    
54. 
    
55. ;;;-----------------------------------------------------------;;
    
56. ;;;获取物体Objects                                            ;;
    
57. ;;;-----------------------------------------------------------;;
    
58. (defun Ent:Ents->Objs (sel / i e o l)
    
59.   (setq i (sslength sel))
    
60.   (repeat i
    
61.     (setq e (ssname sel (setq i (1- i))))
    
62.     (setq o (vlax-ename->vla-object e))
    
63.     (setq l (cons o l))
    
64.   )
    
65. )
    
66. 
    

;;;3d环形阵列测试<br /> (defun c:PAT(/ ss N P1 P2 ActDoc ObjLst)                   ;c:PolarArrayTest<br />   (setq ss (ssget))<br />   (initget 7)<br />   (setq N (getInt "\n数量:"))<br />   (initget 9)<br />   (setq P1 (getpoint "\n中心点:"))<br />   (initget 9)<br />   (setq p2 (getpoint P1 "\n另一点:"))<br />   (if ss<br />     (progn<br />       (setq ActDoc (vla-get-ActiveDocument (vlax-get-acad-object)))<br />       (vla-StartUndoMark ActDoc)<br />       (setq Ojblst (Ent:Ents->Objs ss))<br />       (Mat:3dPolarArray ojblst N (+ pi pi) T (trans p1 1 0) (trans p2 1 0))<br />       (vla-EndUndoMark ActDoc)<br />       (vlax-release-object ActDoc)<br />       (princ)<br />     )<br />   )<br /> )<br /> <br /> <br /> ;;;-----------------------------------------------------------;;<br /> ;;; 三维环形阵列                                              ;;<br /> ;;; 输入：Objlst    -- 物体集                                 ;;<br /> ;;;       Number    -- 要阵列的个数（包含自身在内）           ;;<br /> ;;;       FillAngle -- 旋转角度                               ;;<br /> ;;;       IsCCW     -- 是否逆时针                             ;;<br /> ;;;       P1        -- 阵列中心点                             ;;<br /> ;;;       P2        -- 阵列轴线的另一点                       ;;<br /> ;;; 输出：阵列的物体列表                                      ;;<br /> ;;;-----------------------------------------------------------;;<br /> (defun Mat:3dPolarArray (Objlst Number FillAngle IsCCW P1 P2 / lst1 lst2 ANG MAT NEW)<br />   (if IsCCW<br />     (setq FillAngle (float FillAngle))<br />     (setq FillAngle (- FillAngle (+ pi pi)))<br />   )<br />   (setq ang (/ FillAngle (1- Number)))<br />   (setq mat (MAT:RotateBy2P P1 P2 ang))<br />   (setq mat (vlax-tmatrix mat))<br />   (repeat (1- Number)<br />     (setq lst1 nil)<br />     (foreach obj ObjLst <br />       (setq new (vla-copy obj))<br />       (vla-transformby new mat)<br />       (setq lst1 (cons new lst1))<br />     )<br />     (setq objLst (reverse lst1))<br />     (setq lst2 (cons objLst lst2))<br />   )<br />   (reverse lst2)<br /> )<br /> <br /> ;;;-----------------------------------------------------------;;<br /> ;;;获取物体Objects                                            ;;<br /> ;;;-----------------------------------------------------------;;<br /> (defun Ent:Ents->Objs (sel / i e o l)<br />   (setq i (sslength sel))<br />   (repeat i<br />     (setq e (ssname sel (setq i (1- i))))<br />     (setq o (vlax-ename->vla-object e))<br />     (setq l (cons o l))<br />   )<br /> )<br /> <br />