在固定区域内生成均匀随机点

aimisiyou · 发表于 2020-2-12 01:54:44

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

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

×

本帖最后由 aimisiyou 于 2020-2-12 08:22 编辑

如图，区域内随机生成的点均匀性不是很好，这跟随机函数有很大关系。
(defun fun_point (p1 p2 p3)
  (defun rnd ()
   (*(rem (getvar "cputicks") 1e4) 1e-4)
  )
(command "pline" p1 p2 p3 "c")
  (setq ptlst nil)
  (repeat 200
   (setq r1 (rnd) r2 (rnd))
      (setq ptlst (cons (list
                              (+ (* (- 1 r2) (- 1 r1) (car p1)) (* (- 1 r2) r1 (car p2)) (* r2 (car p3)))
                                 (+ (* (- 1 r2) (- 1 r1) (cadr p1)) (* (- 1 r2) r1 (cadr p2)) (* r2 (cadr p3)))
                           )
                                       ptlst
                     )
      )
)
(mapcar '(lambda (x)
         (entmake
                           (list
                              '(0 . "point")
                              '(100 . "AcDbEntity")
                              '(100 . "AcDbPoint")
                              '(67 . 0)
                              (cons 10 (list (car x) (cadr x)))
                              (cons 210 '(0 0 1))
                           )
                  )
                     )
                     ptlst
)
)
(fun_point '(0 0) '(100 100) '(200 0))

aimisiyou · 发表于 2020-2-12 10:10:43

本帖最后由 aimisiyou 于 2020-2-12 10:21 编辑

既然改变不了随机函数的特性，那么我们就依次调整下三个点的排列位置，尽量得到"较均匀分布"的点。
(defun fun_point (p1 p2 p3)
  (defun rnd ()
   (*(rem (getvar "cputicks") 1e4) 1e-4)
  )
(command "pline" p1 p2 p3 "c")
  (setq ptlst nil)
  (repeat 100
   (setq r1 (abs (cos (* pi (/ (atan (log (rnd))) (/ pi 2.0))))))
   (setq r2 (abs (cos (* pi (/ (atan (log (rnd))) (/ pi 2.0))))))
      (setq ptlst (cons (list
                              (+ (* (- 1 r2) (- 1 r1) (car p1)) (* (- 1 r2) r1 (car p2)) (* r2 (car p3)))
                                 (+ (* (- 1 r2) (- 1 r1) (cadr p1)) (* (- 1 r2) r1 (cadr p2)) (* r2 (cadr p3)))
                           )
                                       ptlst
                     )
      )
   (setq r1 (abs (cos (* pi (/ (atan (log (rnd))) (/ pi 2.0))))))
   (setq r2 (abs (cos (* pi (/ (atan (log (rnd))) (/ pi 2.0))))))
      (setq ptlst (cons (list
                              (+ (* (- 1 r2) (- 1 r1) (car p2)) (* (- 1 r2) r1 (car p3)) (* r2 (car p1)))
                                 (+ (* (- 1 r2) (- 1 r1) (cadr p2)) (* (- 1 r2) r1 (cadr p3)) (* r2 (cadr p1)))
                           )
                                       ptlst
                     )
      )
   (setq r1 (abs (cos (* pi (/ (atan (log (rnd))) (/ pi 2.0))))))
   (setq r2 (abs (cos (* pi (/ (atan (log (rnd))) (/ pi 2.0))))))
      (setq ptlst (cons (list
                              (+ (* (- 1 r2) (- 1 r1) (car p3)) (* (- 1 r2) r1 (car p1)) (* r2 (car p2)))
                                 (+ (* (- 1 r2) (- 1 r1) (cadr p3)) (* (- 1 r2) r1 (cadr p1)) (* r2 (cadr p2)))
                           )
                                       ptlst
                     )
      )
)
(mapcar '(lambda (x)
         (entmake
                           (list
                              '(0 . "point")
                              '(100 . "AcDbEntity")
                              '(100 . "AcDbPoint")
                              '(67 . 0)
                              (cons 10 (list (car x) (cadr x)))
                              (cons 210 '(0 0 1))
                           )
                  )
                     )
                     ptlst
)
)
(fun_point '(0 0) '(100 100) '(200 0))

aimisiyou · 发表于 2020-2-12 10:05:52

为了拉开两个连续随机数的距离，可采用log函数。经改进后，效果图有些进步，但仍不理想，靠近P3处密集，P2处稀疏，P1处介于两者之间。

aisong220 · 发表于 2020-2-14 00:04:52

aimisiyou · 发表于 2020-2-16 10:41:53

可应用到平面找点的问题。随机生成点后，可获得各点对应的函数值，整体有一个最优值，定义局部半径r，局部内的点向局部最优值移动，局部最优值向整体最优值移动，整体最优值随机小范围扰动。

账号		自动登录	找回密码
密码			立即注册

[研讨] 在固定区域内生成均匀随机点

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。