怎么产生随机数

小妹丁

meflying发表于2004-10-24 18:27:00我用如下程序测试: (defun c:test() (repeat 10 (princ \"\n\") (princ (rnd 1 10)) )) 结果如下,试了几次,发现还是不行的,连续...

这是因为计算机运行太快，要让每次取种子数的时间差＞0.1s以上，当然还可以改变范围等办法 。 dispbbs.asp?BoardID=3&ID=15667

xyp1964

随机数程序：加载sjsh.fas，运行sjsh。 随机数个数为30的结果： 12
15
25
42
60
78
13
28
47
81
88
50
63
24
40
2
18
80
19
74
8
60
92
11
52
77
91
96
97
22 不知是否有规律，再修改后公布源程序。

meflying

请10楼和11楼再看清楚前面的帖子,现在讨论的内容已经转到"连续"产生随机数了,楼上的这个虽然连续产生了,不过速度太慢,其实是程序已经将连续内容断开成不连续的了...而且提供FAS在这里好象没有什么意义... 结果也不尽如人意: Command: sjsh
随机数个数: 100 92
6
17
31
53
8
67
99
8
36
45
69
92
16
15
15
15
16
24
24
23
23
25
25
25
25
25
25
25
25
25
25
25
25
25
27
27
27
26
27
26
26
27
27
26
28
28
28
28
28
28
28
28
28
28
28
28
3
29
31
31
31
31
32
32
31
31
33
33
33
32
33
33
33
33
33
33
33
33
33
33
34
34
34
34
34
34
34
34
34
36
36
36
36
36
36
36
36
36
36

wdb

也可以用vc编一个程序，将连续产生的随机数写在注册表里或记事本里，再用VL读出来。


用start命令调用VC程序。


直接用VL连续产生随机数确实有点困难。

雪山飞狐_lzh

(defun c:t1( / a)
(setq a (tls-vbs))
(repeat 10
(print (vbs-eval a 'rnd '()))
)
) (defun tls-join(chrs fchr / pstr)
(setq pstr "")
(foreach i Chrs (setq pstr (strcat pstr i fchr)))
(substr pstr 1 (- (strlen pstr) (strlen fchr)))
) (defun tls-vbs()
(setq Tls-VBScript (vlax-create-Object "ScriptControl"))
(vlax-put-property Tls-VBScript "language" "vbs")
Tls-VBScript
) (defun vbs-eval(vbs funname args)
(vla-eval
vbs
(strcat
(vl-prin1-to-string funname)
"("
(tls-join
(mapcar 'vl-prin1-to-string args)
","
)
")"
)
)
)

小妹丁

这个结果出数较均匀，接近4／5。 命令: test
20
28
32
40
44
2
9
11
21
25
33
40
44
15
6
14
17
25
29
37
19
49
6
10
18
22
23
37
41
49
3
11
27
22
26
34
42
46
3
30
15
19
27
34
38
34
50
8
15
1919
命令: (defun c:test()
(repeat 50
(princ "\n")
(princ (rnd 1 50))
)
)
(defun rnd (n1 n2)
(setq ra 66791
rb 89；此值113也可得到较好的结果。
nn 66791
)
(setq imin (expt 2 31)
imax (1- imin)
)
(if (> n1 n2)
(setq n3 n1
n1 n2
n2 n3
)
)
(setq rn (atoi (substr (rtos (getvar "cdate") 2 7) 14)))
(repeat nn
(setq rn (+ (* rn ra) rb))
(if (minusp rn)
(setq rn (- rn imin))
)
)
(setq rn (rem rn imax))
(setq rn (/ rn imax 1.0))
(setq rn (+ (fix (* rn (1+ (- n2 n1)))) n1))
)

meflying

看看我13楼的话,有部分是针对你这样的情况的...


15楼的方法应该是比较好的了

hyinghao

这个办法不好

4. ;;---------------------------------------------------------------------------;;
5. ;; initialize the global
6. ;;---------------------------------------------------------------------------;;
7. (setq *SeedRand* nil)
8. ;;---------------------------------------------------------------------------;;
9. ;; the base function for random
10. ;;---------------------------------------------------------------------------;;
11. (defun rand16 (/ s)
12.   ;; when this is used for the first time, initialize the seed
13.   ;; from the system clock.
14.   (if (null *SeedRand*)
15.     (progn
16.       (setq s (getvar "date"))      
17.       (setq *SeedRand* (fix (* 86400 (- s (fix s)))))
18.     )
19.     ;; progn
20.   )
21.   ;; if
22.   ;; To generate a psudo-sandom number sequence
23.   ;; I use the routine described in Kernighan and Ritchie's
24.   ;; "C Programming Language" second edition, p46
25.   (setq *SeedRand* (+ (* *SeedRand* 1103515245) 12345))
26.   ;; trim off the bits left of the 16th bits      
27.   (logand (/ *SeedRand* 65536) 32767)
28. )
29. ;;---------------------------------------------------------------------------;;
30. ;; generates a random number between min and max.                            ;;
31. ;; min and max must be a non-negative integer smaller than 32678.            ;;
32. ;;---------------------------------------------------------------------------;;
33. (defun rand (umin umax / r16 range quotient remainder result)
34.   (setq r16 (rand16))
35.   ;; random number smaller than 32678
36.   (setq range (+ 1 (- umax umin)))
37.   ;; number of integers to be produced
38.   (setq quotient (/ r16 range))
39.   ;; result in non-neg. integer
40.   (setq remainder (- r16 (* quotient range)))
41.   (setq result (+ umin remainder))
42.   result
43. )

;;---------------------------------------------------------------------------;; ;; initialize the global ;;---------------------------------------------------------------------------;; (setq *SeedRand* nil) ;;---------------------------------------------------------------------------;; ;; the base function for random ;;---------------------------------------------------------------------------;; (defun rand16 (/ s)   ;; when this is used for the first time, initialize the seed   ;; from the system clock.   (if (null *SeedRand*)     (progn       (setq s (getvar "date"))             (setq *SeedRand* (fix (* 86400 (- s (fix s)))))     )     ;; progn   )   ;; if   ;; To generate a psudo-sandom number sequence   ;; I use the routine described in Kernighan and Ritchie's   ;; "C Programming Language" second edition, p46   (setq *SeedRand* (+ (* *SeedRand* 1103515245) 12345))   ;; trim off the bits left of the 16th bits         (logand (/ *SeedRand* 65536) 32767) ) ;;---------------------------------------------------------------------------;; ;; generates a random number between min and max.                            ;; ;; min and max must be a non-negative integer smaller than 32678.            ;; ;;---------------------------------------------------------------------------;; (defun rand (umin umax / r16 range quotient remainder result)   (setq r16 (rand16))   ;; random number smaller than 32678   (setq range (+ 1 (- umax umin)))   ;; number of integers to be produced   (setq quotient (/ r16 range))   ;; result in non-neg. integer   (setq remainder (- r16 (* quotient range)))   (setq result (+ umin remainder))   result )

我用了上面的代码,比刚才稍微好一点,但也不行,郁闷啊
Autolisp做随机都不行

TANER

3. ;;;晓东空间上讨论的结果。
4. ;;; ==================================================================
5. ;;; 求随机数(0~1)－－ｂｙ　狂刀
6. (defun th-rnd ()                       ; 随机数种子
7.   (* (rem (getvar "cputicks") 1e3) 1e-3)
8. )
9. ;;; ==================================================================
10. ;;; 求ｎ个ａ～ｂ随机数列表----ｂｙ　狂刀
11. (defun th-rndx1-list (a b n  / c lst)   ; 随机数表
12.   (setq c (- b a))
13.   (repeat n
14.     (setq lst (cons (+ a (* (th-rnd) c)) lst))
15.   )
16. )
17. ;;; ==================================================================
18. ;;; 返回在a~b內的n個雜檆整數的列表
19. (defun th-rndx2-list (a b n  / lst c)   ; 随机整数表
20.   (setq c (- b a))
21.   (repeat n
22.     (setq lst (cons (fix (+ a (* (th-rnd) c))) lst))
23.   )
24. )
25. ;;; ==================================================================
26. (defun th-rnd-single (a b)             ; 单个随机整数
27.   (car (th-rndx2-list a b 1))
28. )
29. ;;; ==================================================================

;;;晓东空间上讨论的结果。 ;;; ================================================================== ;;; 求随机数(0~1)－－ｂｙ　狂刀 (defun th-rnd ()                       ; 随机数种子   (* (rem (getvar "cputicks") 1e3) 1e-3) ) ;;; ================================================================== ;;; 求ｎ个ａ～ｂ随机数列表----ｂｙ　狂刀 (defun th-rndx1-list (a b n  / c lst)   ; 随机数表   (setq c (- b a))   (repeat n     (setq lst (cons (+ a (* (th-rnd) c)) lst))   ) ) ;;; ================================================================== ;;; 返回在a~b內的n個雜檆整數的列表 (defun th-rndx2-list (a b n  / lst c)   ; 随机整数表   (setq c (- b a))   (repeat n     (setq lst (cons (fix (+ a (* (th-rnd) c))) lst))   ) ) ;;; ================================================================== (defun th-rnd-single (a b)             ; 单个随机整数   (car (th-rndx2-list a b 1)) ) ;;; ==================================================================

TANER

3. ;;;晓东上讨论的结果
4. ;;; ==================================================================
5. ;;; 求随机数(0~1)－－ｂｙ　狂刀
6. (defun th-rnd ()                       ; 随机数种子
7.   (* (rem (getvar "cputicks") 1e3) 1e-3)
8. )
9. ;;; ==================================================================
10. ;;; 求ｎ个ａ～ｂ随机数列表----ｂｙ　狂刀
11. (defun th-rndx1-list (a b n / c lst)   ; 随机数表
12.   (setq c (- b a))
13.   (repeat n
14.     (setq lst (cons (+ a (* (th-rnd) c)) lst))
15.   )
16. )
17. ;;; ==================================================================
18. ;;; 返回在a~b內的n個雜檆整數的列表
19. (defun th-rndx2-list (a b n / lst c)   ; 随机整数表
20.   (setq c (- b a))
21.   (repeat n
22.     (setq lst (cons (fix (+ a (* (th-rnd) c))) lst))
23.   )
24. )
25. ;;; ==================================================================
26. (defun th-rnd-single (a b)             ; 单个随机整数
27.   (car (th-rndx2-list a b 1))
28. )
29. ;;; ==================================================================

;;;晓东上讨论的结果 ;;; ================================================================== ;;; 求随机数(0~1)－－ｂｙ　狂刀 (defun th-rnd ()                       ; 随机数种子   (* (rem (getvar "cputicks") 1e3) 1e-3) ) ;;; ================================================================== ;;; 求ｎ个ａ～ｂ随机数列表----ｂｙ　狂刀 (defun th-rndx1-list (a b n / c lst)   ; 随机数表   (setq c (- b a))   (repeat n     (setq lst (cons (+ a (* (th-rnd) c)) lst))   ) ) ;;; ================================================================== ;;; 返回在a~b內的n個雜檆整數的列表 (defun th-rndx2-list (a b n / lst c)   ; 随机整数表   (setq c (- b a))   (repeat n     (setq lst (cons (fix (+ a (* (th-rnd) c))) lst))   ) ) ;;; ================================================================== (defun th-rnd-single (a b)             ; 单个随机整数   (car (th-rndx2-list a b 1)) ) ;;; ==================================================================