【飞鸟集】最小包围圆的最佳算法

chlh_jd

hifly兄这个源码确实有更高的效率，没有仔细推敲，不知能否利用交叉引用递归来提高效率
暂时发现in2这个函数改为递归效率会更高些

3. ;;判断点集是否在圆内----------------------
4. ;;;改为递归写法
5. (defun in2 (ptl cen r)
6.   (cond ((null (cadr ptl))
7.   (in1 (car ptl) cen r)
8. )
9. ((in1 (car ptl) cen r)
10.   (in2 (cdr ptl) cen r)
11. )
12. (t nil)
13.   )
14. )

;;判断点集是否在圆内---------------------- ;;;改为递归写法 (defun in2 (ptl cen r)   (cond ((null (cadr ptl))   (in1 (car ptl) cen r) ) ((in1 (car ptl) cen r)   (in2 (cdr ptl) cen r) ) (t nil)   ) )

chlh_jd

另外在4pc函数里面，有个疑问foreach里面是否需要重新定义minr为r

3. ;;定义四点的最小圆圆心半径，并返回三点坐标
4. (defun 4pc (p1 p2 p3 ptmax / pts mind minr r 4ps)
5.   (setq pts (list (3pc p1 p2 ptmax)
6.     (3pc p1 p3 ptmax)
7.     (3pc p2 p3 ptmax)
8.      )
9.   )
10.   (setq 4ps (list p1 p2 p3 ptmax))
11.   (setq minr 1e308)
12.   (foreach n pts
13.     (setq r (cadr n))
14.     (if (and (< r minr)
15.       (in2 4ps (car n) r)
16. )
17.       (setq mind n
18.      minr r);_这里是否需要加上重定义minr为r
19.     )
20.   )
21.   mind
22. )

;;定义四点的最小圆圆心半径，并返回三点坐标 (defun 4pc (p1 p2 p3 ptmax / pts mind minr r 4ps)   (setq pts (list (3pc p1 p2 ptmax)     (3pc p1 p3 ptmax)     (3pc p2 p3 ptmax)      )   )   (setq 4ps (list p1 p2 p3 ptmax))   (setq minr 1e308)   (foreach n pts     (setq r (cadr n))     (if (and (< r minr)       (in2 4ps (car n) r) )       (setq mind n      minr r);_这里是否需要加上重定义minr为r     )   )   mind )

highflybird

chlh_jd发表于2010-9-20 0:26:00另外在4pc函数里面，有个疑问foreach里面是否需要重新定义minr为r 以下内容为程序代码: ;;定义四点的最小圆圆心半径，并返回三点坐标 (defun 4pc (p1 p2 p3 ptmax / pts m

chlh_jd

也祝你中秋快乐！

顺祝大家中秋快乐！~

仲文玉

楼主，如果是最小正三角形有没有试过，或者正方形，正N边形

highflybird

仲文玉发表于2010-10-18 8:17:00楼主，如果是最小正三角形有没有试过，或者正方形，正N边形

没问题，就算是正多边形，算法也成立。

highflybird

chlh_jd发表于2010-9-20 0:26:00另外在4pc函数里面，有个疑问foreach里面是否需要重新定义minr为r 以下内容为程序代码: ;;定义四点的最小圆圆心半径，并返回三点坐标 (defun 4pc (p1 p2 p3 ptmax / pts m

461045462

  来明经这么久，今天才闯到了这里。

似乎又找到了一个芝麻开门的地方

来的太迟了

wwwcncomcn

highflybir 发表于 2009-2-27 02:41
原以为这样的文章可能没多少实际用途，被冷藏了这么久。想不到还有点用处啊。下面的程序可以给你解决这个烦 ...

楼主啊，我新来的，没币了，能不能施舍一下你的附件。。。。

monkeylzx

多谢分享，来学习一下