[风之影][Lisp大挑战第二季]光路计算

VBALISPER

哇,又有看的了.

prolee2008

cabinsummer 发表于 2011-11-7 22:50
先对prolee2008的程序表示赞赏，不过程序有点小小的瑕疵。在上面的图片中，曲线为折射面，入射点位于左 ...

真的很不好意思呀，折射面现在还只能是样条线.圆弧忘了写进去了。

jicqj

坐等高手 我来看热闹

prolee2008

尽快出第三季吧

cabinsummer

prolee2008 发表于 2011-11-9 11:40
尽快出第三季吧

第二季还没完全搞定，不过第三季已经在酝酿中了。你可以先去挑战一下第一季，感觉你很厉害啊

cabinsummer

提供思路
prolee2008给了我们一个很好的开端，但是其有不足之处。比如，他的程序中采用了(command "line")和(command "rotate")这样的函数。凡是用command直接在屏幕上画线的，其效率肯定提不起来，因为要不断的刷显存，确保显示后才能进行下一步命令。所以对于大数据量的计算，最好不要用command，把所有的数据准备好后，再用entmake可以显著提高效率。

cabinsummer

这是asen在2008年4月27日发表的帖子中的一个函数，见
http://bbs.mjtd.com/forum.php?mod=viewthread&tid=67300

1. (defun c:n (/ obj p1 p2 p3 li)
   
2.   (setq obj (vlax-ename->vla-object (car (entsel "\ncurve: ")))
   
3.          p1  (trans (getpoint "\npoint 1: ") 1 0)
   
4.          p2  (trans (getpoint p1 "\npoint 2: ") 1 0)
   
5.         li (vla-addline (vla-get-ModelSpace (vla-get-ActiveDocument (vlax-get-acad-object)))(vlax-3d-point p1)(vlax-3d-point p2))
   
6.         p3  (vlax-invoke obj 'IntersectWith li acExtendNone)
   
7.   )
   
8.   (if p3
   
9.     (vla-put-EndPoint li (vlax-3d-point (list (car p3) (cadr p3) (caddr p3))))
   
10.     (vla-delete li)
    
11.   )
    
12.   (princ)
    
13. )

(defun c:n (/ obj p1 p2 p3 li)<br />   (setq obj (vlax-ename->vla-object (car (entsel "\ncurve: ")))<br />          p1  (trans (getpoint "\npoint 1: ") 1 0)<br />          p2  (trans (getpoint p1 "\npoint 2: ") 1 0)<br />         li (vla-addline (vla-get-ModelSpace (vla-get-ActiveDocument (vlax-get-acad-object)))(vlax-3d-point p1)(vlax-3d-point p2))<br />         p3  (vlax-invoke obj 'IntersectWith li acExtendNone)<br />   )<br />   (if p3<br />     (vla-put-EndPoint li (vlax-3d-point (list (car p3) (cadr p3) (caddr p3))))<br />     (vla-delete li)<br />   )<br />   (princ)<br /> )

函数可以直接求出直线与曲线的交点，但缺点和prolee2008一样，都是求最远的交点，而不是最近的交点。

cabinsummer

我发现挑战还是有一定难度的，如何求出直线与曲线最近的交点？

jxcaixiaomeng

很好的资料啊！呵呵      

vormittag

cabinsummer 发表于 2011-11-10 20:09
这是asen在2008年4月27日发表的帖子中的一个函数，见
http://bbs.mjtd.com/forum.php?mod=viewthread&tid= ...

vla-get-intersectwith 获得的是交点列表，每三个数是一个交点的三个坐标，第一组三个数未必是最远交点的坐标，这是和曲线方向有关系的，如果把光源点和折射面沿光路传播方向mirror一下，可能第一个交点就是最近交点了。