一种产生随机数的方法

Highflybird

对Alisp来说，没有生成随机数的函数，一般来说，有如下方法可以产生它：
1、用vbscript 或者Jscirpt等内部的随机函数产生，然后在lisp用调用脚本语言。
2、利用CAD的系统变量，如cputicks,TDUSRTIMER等产生随机数。
3、利用种子数产生随机数。可以是纯LISP的函数。

这里介绍随机数产生的另外一种方法。

(setq rndObj (vlax-create-object "System.Random"))

(vlax-invoke rndobj 'nextdouble)==>此处产生一个0到1之间的随机小数
(vlax-invoke rndobj 'next) ==>此处产生一个32位的正整数。

如何利用这个做成一个一定范围内的随机函数，请大家自己试验。

Highflybird

关于其用法，可以参见如下链接：
https://docs.microsoft.com/zh-cn ... =netframework-4.7.2
对LISP来说，有一定限制， 貌似只能使用 next 和nextdouble 方法。
另外  vlax-invoke 和 vlax-get 都可以产生结果。
譬如也可以使用 (vlax-get rndobj 'nextdouble)  

marting

Highflybird 发表于 2019-3-13 23:18
关于其用法，可以参见如下链接：
https://docs.microsoft.com/zh-cn/dotnet/api/system.random?redirected ...

这里有个随机数生成的函数，用REPEAT循环，极短时间内生成的随机数看着分布挺好的

[Actionscript3] 纯文本查看 复制代码

(defun rand ()
  (setq seed (if seed
           (rem (+ (* seed 15625.7) 0.21137152) 1)
           0.3171943
         )
  )
)

(repeat 10 (princ "\n")(princ(rand)))

0.0126488
0.858366
0.778777
0.152224
0.814014
0.756342
0.587726
0.849062
0.392444
0.421463

Highflybird

marting 发表于 2019-3-13 23:20
这里有个随机数生成的函数，用REPEAT循环，极短时间内生成的随机数看着分布挺好的

[mw_shl_code=actio ...

随机数的函数关键要看它的均匀性，满足正态分布。
你上面的这个有一定局限性，我以前测试过。譬如用你说的这个函数，产生随机点，如果数量多了，很可能产生重复点。

marting

Highflybird 发表于 2019-3-13 23:37
随机数的函数关键要看它的均匀性，满足正态分布。
你上面的这个有一定局限性，我以前测试过。譬如用你说 ...

随机数其实都是“伪随机”，上面的函数简单，乘以100000，出现重复的概率应该很低，可以通过别的方法规避掉。

Highflybird

marting 发表于 2019-3-13 23:46
随机数其实都是“伪随机”，上面的函数简单，乘以100000，出现重复的概率应该很低，可以通过别的方法规避 ...

简单的不一定好，虽然可能速度快，但也许会产生问题。
我以前测试的是10万个点，用这种随机方法，居然产生了5个重复的点，而且多次试验，都多多少少有几个重复的。
用system.random其实速度很快，
留意我引用的链接里面的一段话：
当前实现Random类基于 Donald E.Knuth 删减随机数生成器算法的修改版本。 有关详细信息，请参阅 D.e。 Knuth。 计算机编程，卷 2 的技术：Seminumerical 算法。 Addison-wesley，Reading，MA，第三个版本，1997年。

sicky111

又研究起随机数了。

Highflybird

1. ;;;
   
2. ;;; https://docs.microsoft.com/zh-cn ... =netframework-4.7.2
   
3. (defun C:RND()
   
4.   (if
   
5.     (or
   
6.       (setq wcs (vlax-create-object "Aec32BitAppServer.AecScriptControl.1"))
   
7.       (setq wcs (vlax-create-object "ScriptControl"))
   
8.     )
   
9.     (vlax-put-property wcs "language" "VBScript")
   
10.   )
    
11.   (setq str "Randomize
    
12.              Function Rand(x,y)\n
    
13.                  Rand=x+Rnd*(y-x)\n
    
14.              End Function"
    
15.   )       
    
16.   (vlax-invoke wcs 'addcode str)
    
17. 
    
18.   (defun wcsRnd(a b)
    
19.     (vlax-invoke wcs 'run "rand" a b)
    
20.   )
    
21. 
    
22.   (setq rndObj (vlax-create-object "System.Random"))
    
23.   (setq x (vlax-invoke rndobj 'nextdouble))
    
24.   (setq y (vlax-invoke rndobj 'next))
    
25. 
    
26.   (defun sysRnd (x y)
    
27.     (+ x (* (vlax-invoke rndobj 'nextdouble) (- y x)))
    
28.   )
    
29. 
    
30.   (defun SEEDrnd (x y)
    
31.     (setq seed (if seed
    
32.            (rem (+ (* seed 15625.7) 0.21137152) 1)
    
33.            0.3171943
    
34.          )
    
35.     )
    
36.     (+ x (* seed (- y x)))
    
37.   )
    
38. 
    
39.   (uti:bench 100000
    
40.     (list
    
41.       (list 'wcsRnd 1 100)
    
42.       (list 'sysRnd 1 100)
    
43.       (list 'seedRnd 1 100)
    
44.     )
    
45.   )
    
46.   (princ)
    
47. )
    

;;; <br /> ;;; <a href=\"https://docs.microsoft.com/zh-cn/dotnet/api/system.random?redirectedfrom=MSDN&view=netframework-4.7.2\" target=\"_blank\">https://docs.microsoft.com/zh-cn ... =netframework-4.7.2</a><br /> (defun C:RND()<br />   (if<br />     (or<br />       (setq wcs (vlax-create-object "Aec32BitAppServer.AecScriptControl.1"))<br />       (setq wcs (vlax-create-object "ScriptControl"))<br />     )<br />     (vlax-put-property wcs "language" "VBScript")<br />   )<br />   (setq str "Randomize<br />              Function Rand(x,y)\n<br />                  Rand=x+Rnd*(y-x)\n<br />              End Function"<br />   )        <br />   (vlax-invoke wcs 'addcode str)<br /> <br />   (defun wcsRnd(a b)<br />     (vlax-invoke wcs 'run "rand" a b)<br />   )<br /> <br />   (setq rndObj (vlax-create-object "System.Random"))<br />   (setq x (vlax-invoke rndobj 'nextdouble))<br />   (setq y (vlax-invoke rndobj 'next))<br /> <br />   (defun sysRnd (x y)<br />     (+ x (* (vlax-invoke rndobj 'nextdouble) (- y x)))<br />   )<br /> <br />   (defun SEEDrnd (x y)<br />     (setq seed (if seed<br />            (rem (+ (* seed 15625.7) 0.21137152) 1)<br />            0.3171943<br />          )<br />     )<br />     (+ x (* seed (- y x)))<br />   )<br /> <br />   (uti:bench 100000<br />     (list<br />       (list 'wcsRnd 1 100)<br />       (list 'sysRnd 1 100)<br />       (list 'seedRnd 1 100)<br />     )<br />   )<br />   (princ)<br /> )<br />

上面的代码是为测试三种代码的效率用的，经过验证，用同余的速度最快，vbs最慢，system.random中间，但同余法的效率并未超过system.random很多，仅仅是几倍的差距。而system.random虽也是伪随机的，但是正确性比同余法高。

zhystar

学习一下，
回来亲自验证

arker

弄个随机买彩票的函数把。

yangjian001

做一个热心并受欢迎的人

qyming2000

学习                                      

evayleung

谢谢分享。。。。。。。。。。。。。。。。。。。。。。。。。。。。

kqqt6236

谢谢分享。