设为首页收藏本站CAD论坛CAD教程CAD下载明经币充值
切换到宽版

明经CAD社区

 找回密码
 注册

QQ登录

只需一步，快速开始

搜索
明经CAD社区»社区 › 开发篇 › AutoLISP/Visual LISP 编程技术 › 关于矩阵的探索之路
返回列表 发新帖
查看: 620|回复: 7

[源码] 关于矩阵的探索之路

[复制链接]
尘缘一生
发表于 2024-12-19 18:09:59 | 显示全部楼层 |阅读模式
本帖最后由 尘缘一生 于 2024-12-19 18:34 编辑

闲来无事，看到本坛的说，COMMAND  MOVE   ROTATE
过往啊，我也尝试，探讨许多次，
总觉得，还是COMMAND的快，并不慢，到底怎么快，拿不准主意
鉴于矩阵的话，比较晦涩难懂，弄一段时间，放弃了，
所以啊，高飞大师这种数学好的，是基础科学，好比华为任正非讲的：“俄罗斯科学家那个故事....”
闲也无聊呢，弄几个函数，测试还是不快呢！
难道移动，旋转最快的方式我们找不到吗？特别是和GRREAD整合做起来，频繁刷新时候，能不卡的。
非得C ARX VB去写吗？ 可惜我早已忘记C++20几年了，拾不起来！
  2. ;实体,选择集,实体表->矩阵移动---(一级)---
  3. (defun sl:move (ss p1 p2)
  4.   (sl:mov-ang ss p1 p2 nil)
  5. )
  6. ;选择集ss以基点p0 旋转ang(弧度)---(一级)-----
  7. ;(sl-rot (ssget) (getpoint) (* pi 0.25))
  8. (defun sl-ssrot (ss p0 ang / mat)
  9.   (setq mat (mat:rotation (trans p0 1 0) ang))
  10.   (sl:sstransformby mat ss)
  11.   (sl-wzgz ss)
  12.   (princ)
  13. )
  14. ;;实体,选择集,实体表->矩阵从p1到p2移动并旋转ang度-----(一级)------
  15. (defun sl:mov-ang (ss p1 p2 ang / mat)
  16.   (if (setq mat (mat:translateby2pang p1 p2 ang))
  17.     (cond
  18.       ((= (type ss) 'ENAME) (vla-transformby (vlax-ename->vla-object ss) (vlax-tmatrix mat))) ;图元
  19.       ((= (type ss) 'PICKSET) (sl:sstransformby mat ss)) ;集
  20.       ((= (type ss) 'LIST) (sl:sstransformby mat (sl:pickset-fromlist ss))) ;实体表
  21.     )
  22.   )
  23. )
  24. ;平移+旋转变换矩阵---(一级)---
  25. ;p1 基点 p2 目标点 并旋转ang(nil 仅平移)度
  26. ;(mat:translateby2pang p1 p2 nil) 平移
  27. ;(mat:translateby2pang p1 p1 ang) 旋转
  28. ;(mat:translateby2pang p1 p2 ang) 平移+旋转
  29. (defun mat:translateby2pang (p1 p2 ang / q mat)
  30.   (setq mat (mat:translateby2p p1 p2))
  31.   (if ang
  32.     (if mat
  33.       (setq q (mat:rotation (trans p1 1 0) ang) mat (mat:mxm mat q)) ;转动矩阵;矩阵相乘:移动+转动
  34.       (cond
  35.         ((and p1 (is_pt p1))
  36.           (setq mat (mat:rotation (trans p1 1 0) ang))
  37.         )
  38.         ((and p2 (is_pt p2))
  39.           (setq mat (mat:rotation (trans p2 1 0) ang))
  40.         )
  41.       )
  42.     )
  43.   )  
  44.   mat  
  45. )
  46. ;选择集矩阵变换---(一级)---
  47. (defun sl:sstransformby (mat ss / i e)
  48.   (setq mat (vlax-tmatrix mat))
  49.   (setq i -1)
  50.   (while (setq e (ssname ss (setq i (1+ i))))
  51.     (vla-transformby (vlax-ename->vla-object e) mat)
  52.   )
  53. )
  54. ;平移变换矩阵---(一级)---
  55. ;p1基点 p2目标点
  56. (defun mat:translateby2p (p1 p2 / mat)
  57.   (if (and p1 p2 (is_pt p1) (is_pt p2) (> (distance p1 p2) 0.01))
  58.     (setq mat (mat:translation (mapcar '- p2 p1)))
  59.   )
  60.   mat
  61. )
  62. ;平移变换矩阵v 位移矢量---(一级)---
  63. (defun mat:translation (v)
  64.   (list
  65.     (list 1. 0. 0. (car v))
  66.     (list 0. 1. 0. (cadr v))
  67.     (list 0. 0. 1. (caddr v))
  68.     (list 0. 0. 0. 1.)
  69.   )
  70. )

这一切，只不过是想弄掉几个command代码罢了，举例三领函数之一；然而，效率不行，矩阵再高级，也绝不可入三领集成的。
  2. ;SLdesign 三领  by 尘缘一生  QQ 15290049
  3. ;;移动-拷贝-旋转-镜像-递增-放缩------(一级)---------
  4. ;;ss:实体、实体表、选择集 p0 移动起点 nil | k: t(move) nil(copy)
  5. (defun sldomov (ss p0 k / e_lst p1 p2 p3 pt pt1 loop ang dis f3 f8 bb p00 a a0 b e1 ee kk px py xx dqxx mk str)
  6.   (setq e_lst (sysvar '("QAFLAGS" "CMDECHO" "NOMUTT" "OSMODE" "ORTHOMODE")) p1 (cadr (grread 5)))
  7.   (cond
  8.     ((= (type ss) 'ENAME) (setq ss (ssadd ss))) ;图元
  9.     ((= (type ss) 'LIST) (setq ss (sl:pickset-fromlist ss))) ;实体表
  10.   )
  11.   (command "_.undo" "be") ;;舍弃高级代码,加速
  12.   (setvar "CMDECHO" 0) ;;命令显示关闭
  13.   (setvar "NOMUTT" 1)
  14.   (setvar "QAFLAGS" 0)
  15.   (if (= p0 nil)
  16.     (setq p0 (ssmpt ss))
  17.     (setq p0 (polar p1 (angle p1 p0) (distance p1 p0)))
  18.   )
  19.   (princ
  20.     (setq str
  21.       (slmsg
  22.         "\n->[逆转90度(TAB)/取角(A)/LIN.(E)/<-mir->(D)/↑mir ↓(S)/mir-F(4)/放大(Q)/缩小(W)/递增(Z)/递减(J)/大一倍(+)/小一半(-)/复制当前(C)/基点(X)/正交(F8)/扑捉(F3)](空格..复制)(左键>当前)(右键>退出)"
  23.         "\n->[癴锣90(TAB)/à(A)/LIN.(E)/<-mir->(D)/◆mir □(S)/mir-F(4)/(Q)/罽(W)/患糤(Z)/患搭(J)/(+)/(-)/確讽玡(C)/膀翴(X)/タユ(F8)/汲(F3)](..確)(オ龄>讽玡)(/龄>癶)"
  24.         "\n->[Rever90(TAB)/Takang(A)/LIN.(E)/<-mir->(D)/↑mir ↓(S)/mir-F(4)/Zoom1.1(Q)/Zoom0.9(W)/Increase(Z)/Decrease(J)/Zoom2.0(+)/Zoom0.5(-)/Copy current(C)/Base point(X)/Orth(F8)/Osnap(F3)]Spaced..Copy/Right>Exit"
  25.       )
  26.     )
  27.   )
  28.   (setq loop t pt1 p0 p2 p0 f8 (getvar "ORTHOMODE") f3 (getvar "OSMODE") dqxx 1)
  29.   (while loop
  30.     (setq bb (grread t 15 2) p00 (cadr bb))
  31.     (cond
  32.       ((equal bb '(2 6));F3切换捕捉开关
  33.         (cond
  34.           ((and (< f3 16384) (/= f3 0))
  35.             (setq f3 (+ f3 16384))
  36.             (prompt (slmsg "\n <对象捕捉 关>" "\n <癸禜 闽>" "\n <OSnap Off>"))
  37.           )
  38.           ((or (= f3 0) (>= f3 16384))
  39.             (setq f3 16383)
  40.             (prompt (slmsg "\n <对象捕捉 开>" "\n <癸禜 秨>" "\n <OSnap On>"))
  41.           )
  42.         )
  43.         (setvar "OSMODE" f3) (redraw)
  44.       )   
  45.       ((equal bb '(2 15))    ;F8切换正交开关
  46.         (if (= f8 0)
  47.           (progn (setq f8 1) (prompt (slmsg "\n <正交 开>" "\n <タユ 秨>" "\n <Orth open>")))
  48.           (progn (setq f8 0) (prompt (slmsg "\n <正交 关>" "\n <タユ 闽>" "\n <Orth off>")))
  49.         )
  50.         (setvar "ORTHOMODE" f8) (redraw)
  51.       )
  52.       ((= (car bb) 5)
  53.         (redraw)
  54.         (if (= f8 1)
  55.           (progn
  56.             (setq ang (angle p2 p00))
  57.             (cond
  58.               ((or (and (> ang pi4) (< ang 3pi4)) (and (> ang 5pi4) (< ang 7pi4)))
  59.                 (setq pt (list (car p2) (cadr p00)))
  60.               )
  61.               (t
  62.                 (setq pt (list (car p00) (cadr p2)))
  63.               )
  64.             )
  65.           )
  66.           (setq pt p00)
  67.         )
  68.         (if (and (<= f3 16384) (> f3 0) (/= f8 1))
  69.           (setq pt (slosnappt ss pt) kk t)
  70.         )
  71.         (if (and (not mk) (> (distance p0 pt) 0.001))
  72.           (command "move" ss "" "_non" p0 "_non" pt)
  73.         )  
  74.         (setq p0 pt)
  75.       )
  76.       ((equal bb '(2 52))  ;; 4mir
  77.         (if mk (setq mk nil) (setq mk t p3 nil a0 nil))
  78.       )
  79.       ((member bb '((2 9)))      ;;table 键
  80.         (command "ROTATE" ss "" "_non" pt 90)
  81.         (setq ss (sl-wzgz ss))
  82.       )
  83.       ((member bb '((2 65) (2 97)))   ;;A 旋转定角
  84.         (ss-rotang ss 0) ;选择集先转0度水平
  85.         (setq a (ss9pt ss 1))
  86.         (setq b (polar a 0 500))
  87.         (command ".rotate" ss "" "_non" a "r" "_non" a "_non" (polar pt (angle pt b) (distance pt b)) pause)
  88.         (setq ss (sl-wzgz ss))
  89.       )
  90.       ((member bb '((2 115) (2 83)))  ;;S s 上下翻
  91.         (command "mirror" ss "" "_non" pt "_non" (mapcar '- pt '(1 0)) "Y")
  92.       )
  93.       ((member bb '((2 100) (2 68))) ;;D d 左右翻
  94.         (command "mirror" ss "" "_non" pt "_non" (mapcar '- pt '(0 1)) "Y")
  95.       )
  96.       ((or (member bb '((2 13))) (= (car bb) 3));;左键、回车
  97.         (setq mk nil)
  98.         (if (= k t) ;;移动
  99.           (setq loop nil)
  100.           (progn
  101.             (command "copy" ss "" "_non" '(0 0) "_non" '(0 0))
  102.             (setq ang (angle pt1 pt) dis (distance pt1 pt) pt1 pt)
  103.           )
  104.         )
  105.       )
  106.       ((member bb '((2 69) (2 101)))   ;;E e 齐线
  107.         (princ (slmsg "\n 请选择要对齐的(线或实体)-->" "\n 叫匡拒璶癸霍(絬┪龟砰)-->" "\n Please select the (line or entity) to align -->"))
  108.         (setq ss (ss-rotang ss 0))
  109.         (setq e1 (car (setq ee (entsel ))))
  110.         (setq ang (e-ang e1 (cadr ee)))
  111.         (setq ss (ss-rotang ss ang))
  112.         (setq ss (sl-wzgz ss))
  113.         (princ str)
  114.       )
  115.       ((member bb '((2 67) (2 99)))   ;;C c 复制在当前
  116.         (command "copy" ss "" "_non" '(0 0) "_non" '(0 0))
  117.         (setq ang (angle pt1 pt) dis (distance pt1 pt) pt1 pt)
  118.       )
  119.       ((and (member bb '((2 32))) ang dis) ;;空格键
  120.         (command "copy" ss "" "_non" '(0 0) "_non" '(0 0))
  121.         (command "MOVE" ss "" "_non" pt1 "_non" (setq pt1 (polar pt1 ang dis))) ;;移位
  122.       )
  123.       ((member bb '((2 87) (2 119)))         ;;缩小  W w
  124.         (command "scale" ss "" "_non" pt "0.9")
  125.       )
  126.       ((member bb '((2 81) (2 113)))         ;;放大  Q q
  127.         (command "scale" ss "" "_non" pt "1.1")
  128.       )
  129.       ((equal bb '(2 45))         ;;缩小一半 -
  130.         (command "scale" ss "" "_non" pt "0.5")
  131.       )
  132.       ((member bb '((2 43) (2 61))) ;;放大一倍 +
  133.         (command "scale" ss "" "_non" pt "2.0")
  134.       )
  135.       ((member bb '((2 90) (2 122)))  ;;递增 Z z
  136.         (setq ss (ss++ ss 1))
  137.       )
  138.       ((member bb '((2 74) (2 106)))  ;;递减 J j
  139.         (setq ss (ss++ ss -1))
  140.       )
  141.       ((member bb '((2 88) (2 120)))  ;;基点 X x
  142.         (command "MOVE" ss "" "_non" pt "_non" p0)  ;;移回去
  143.         (princ (slmsg "\n 指定基点-->" "\n ﹚膀翴-->" "\n Specify Base Point -->"))
  144.         (slslx p0 0)
  145.         (setvar "OSMODE" 16383)
  146.         (setq p1 (getpoint p0))
  147.         (command "MOVE" ss "" "_non" p1 "_non" p0)
  148.         (setvar "OSMODE" f3)
  149.         (princ str)
  150.       )
  151.       ((member (car bb) '(11 25));;右键
  152.         (if (= k t) ;;移动
  153.           (command "MOVE" ss "" "_non" pt "_non" p0)  ;;移回去
  154.           (sl:erase ss)
  155.         )
  156.         (setq loop nil mk nil)
  157.       )
  158.     )
  159.     (if mk
  160.       (progn
  161.         (if (or (not p3) (not a0)) (setq a (ss9pt ss nil) p3 (nth 5 a) a0 (angle (car a) (caddr a))))
  162.         (setq xx (sl-4p p3 p0) px (polar p3 a0 50) py (polar p3 (+ a0 pi2) 50))
  163.         (if (/= dqxx xx)
  164.           (cond
  165.             ((or (= dqxx 1) (= dqxx 4))
  166.               (if (or (= xx 2) (= xx 3)) (command "mirror" ss "" "_non" p3 "_non" py "y"))
  167.               (if (and (/= xx 2) (/= xx 3)) (command "mirror" ss "" "_non" p3 "_non" px "y"))
  168.             )
  169.             ((or (= dqxx 2) (= dqxx 3))
  170.               (if (or (= xx 1) (= xx 4)) (command "mirror" ss "" "_non" p3 "_non" py "y"))
  171.               (if (and (/= xx 1) (/= xx 4)) (command "mirror" ss "" "_non" p3 "_non" px "y"))
  172.             )
  173.           )
  174.         )
  175.         (setq dqxx xx)
  176.       )
  177.     )
  178.   );while
  179.   (redraw)
  180.   (if kk (sl:erase (ssget "X" (list (cons 8 "f-i-n-d")))))
  181.   (command "_.undo" "e")
  182.   (mapcar 'eval e_lst)
  183.   (princ)
  184. )
以上函数，我还是采用了COMMAND多，总觉得，还有未知之路，还有永不止境的完善空间......

本帖子中包含更多资源

您需要 登录 才可以下载或查看，没有账号？注册

x

评分

参与人数 2明经币 +2 收起 理由
tranque + 1 赞一个!
yanshengjiang + 1

查看全部评分

"觉得好，就打赏"
打赏
还没有人打赏，支持一下
回复

使用道具 举报

你有种再说一遍
发表于 2024-12-19 18:23:49 | 显示全部楼层
矩阵不难懂,难的是Lisp敲矩阵,
学之前你得知道数组吧,结果你不知道,
数组可以O(1)寻址,结果你不知道,
别人都是有C++经验来教你,
他是降为打击啊,说到降维,
如果连多维数组就是一维数组都不知道,那么就白学了,
一维数组只需要array[i*j+m]就是二维数组了,
这样矩阵求值只需要面对一维数组,而不需要考虑多维.
再加强一点,就是CPU预读和硬件加速了.
https://www.cnblogs.com/JJBox/p/14300098.html
回复 支持 反对

使用道具 举报

尘缘一生
 楼主| 发表于 2024-12-19 18:28:38 | 显示全部楼层
本帖最后由 尘缘一生 于 2024-12-19 18:36 编辑
你有种再说一遍 发表于 2024-12-19 18:23
矩阵不难懂,难的是Lisp敲矩阵,
学之前你得知道数组吧,结果你不知道,
数组可以O(1)寻址,结果你不知道,

能不能提速，移动，旋转“大选择集”问题，顾版的哪个函数也是快不起来呢。
回复 支持 反对

使用道具 举报

你有种再说一遍
发表于 2024-12-19 18:36:59 | 显示全部楼层
尘缘一生 发表于 2024-12-19 18:28
能不能提速，移动，旋转“大选择集”问题，顾版的函数也是快不起来呢。

你的这个问题并不用矩阵提速,
这个问题更多是Lisp访问图元之前会进行锁符号表,
这些权限工作会令你大规模操作时候发生瓶颈,
每个图元都锁一次表,和一次锁定再使用,速度就不一样了,
即使同样的代码,在C#敲一遍速度也完全不一样了,
所以想要提速还是要转变一下思维,
尤其是C#有并行,并发,等等许许多多操作.

点评

尘缘一生
“并行”二字，我明白了。  发表于 2024-12-19 18:39
回复 支持 反对

使用道具 举报

ddyer
发表于 2024-12-19 18:39:23 | 显示全部楼层
每天看会天书，嗯，总有一天会起飞的。
回复 支持 反对

使用道具 举报

liuhe
发表于 2024-12-20 09:39:41 | 显示全部楼层
lsp  不适用于追求速度这种，它的特点灵活简单好用，大家学习高级算法，只是便于理解各种底层逻辑。高飞鸟大佬，人家其实玩的是c++  ，学习矩阵是必须过程
回复 支持 反对

使用道具 举报

寒潮大冬瓜
发表于 2024-12-21 09:35:56 | 显示全部楼层
飞诗的提示

本帖子中包含更多资源

您需要 登录 才可以下载或查看，没有账号？注册

x
回复 支持 反对

使用道具 举报

返回列表 发新帖
高级模式
B Color Image Link Quote Code Smilies
您需要登录后才可以回帖 登录 | 注册

本版积分规则

小黑屋|手机版|CAD论坛|CAD教程|CAD下载|联系我们|关于明经|明经通道 ( 粤ICP备05003914号 )  
©2000-2023 明经通道 版权所有 本站代码，在未取得本站及作者授权的情况下，不得用于商业用途

GMT+8, 2025-1-22 21:48 , Processed in 0.192157 second(s), 27 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表