cad生成最小包围盒lisp_点集的最小包围圆(用LISP求解)

;;;************************************

;;;求最小包围圆的lisp程序--------------

;;;其算法为参见了有关文献--------------

;;;这种算法在退化很严重的情况结果也正确

;;;其中程序主段是核心算法，其他的附加程

;;;序为取点，画点，画圆和半径，用来测试

;;;************************************

(defun C:test (/ )

;;取点，画点，并对函数用时计算-------

(setq ss (ssget '((0 . "POINT,LINE,POLYLINE,LWPOLYLINE"))))

(setq pts (ssgetpoint ss))

(setq t0 (getvar "TDUSRTIMER"))

;;(setq x (mincir pts))

(setq x (mdesc pts nil))

(princ "n用时")

(princ (* (- (getvar "TDUSRTIMER") t0) 86400)) ;结束计时

(princ "秒")

(if (null x)

(alert "点的有效数目太小，请重新输入!")

(progn

(setq cen (car x)

rad (cdr x)

;;rad (cadr x)

;;ptmax (caddr x)

)

;;画圆及半径，列出圆的圆心半径值

(make-circle cen rad)

;;(make-line cen ptmax)

)

)

(princ)

)

;;;************************************

;;;求最小包围圆的函数，空集返回空集，否

;;;则返回最小圆的圆心，半径和圆上的一点

;;;这是程序的主段----------------------

;;;************************************

(defun mincir (ptlist / CEN CEN_R P1 P2 P3 PTMAX R rad X i)

;;判断有效点个数-