sunnyso

回復 10# sjgau

thanks in advance

sjgau

回復 11# sunnyso


    integer 的範圍，只有到 +32767, 不行使用

long 的範圍，雖然有到 +2147483647
但是，下一個 亂數，等於 前一個亂數，乘上 a, a= 7^5= 16807
就爆掉了

所以，必須使用 double,

single 的精確度，只有 2^24, 大約是 10^7, 不夠使用
所以，必須使用 double

當然，必須 反覆測試 每一個 細節

stillfish00

回復 8# sjgau
自己寫真的可以發現很多問題，學到很多，
希望能拋磚引玉，我是採用Microsoft Visual Basic 的那組參數，仿造rnd函數一樣傳回single。

1. Function myRnd(Optional Number) As Single
   
2.   Dim tmp As Double
   
3.   'LCG Algorithm：
   
4.   '  X(n+1) = (a * X(n) + c) mod m
   
5.   Const a As Long = 1140671485
   
6.   Const c As Long = 12820163
   
7.   Const m As Long = 16777216  '2^24
   
8.   Const x0 As Long = 1     '根據你使用哪組參數，x0有時有需要符合的條件。
   
9.   
   
10.   Static seed As Variant
    
11.   '利用IsEmpty判斷，來初始化一 static變數為非零值。
    
12.   If IsEmpty(seed) Then seed = x0
    
13.   
    
14.   '利用IsMissing判斷 optional 參數是否傳入，若函數有傳入參數時，設為種子。
    
15.   If Not IsMissing(Number) Then seed = Number
    
16.    
    
17.   'myRnd = (a * seed + c) Mod m
    
18.   '(a * seed + c)過大時，mod 會發生溢位，只能用其他方法求。
    
19.   
    
20.   tmp = a * seed + c  'double
    
21.   seed = tmp - m * Fix(tmp / m)   '即計算 (a * seed + c) Mod m，此為產生的亂數值，並設為下次種子
    
22.   myRnd = seed / m  '使傳回亂數值介於0到1之間
    
23. End Function

複製代碼

stillfish00

回復 13# stillfish00
實測後仍不滿意，呼叫m次後，不會回到第一組產生的隨機數，
主因是 tmp = a * seed + c 即便是 double 數仍存有精確度誤差...   

LongLong又只有在64位元平台上可以使用

c_c_lai

回復 14# stillfish00
這是我測試的結果：
A.  rng = myRnd(sngVal)

下載 (18.04 KB)

2013-11-1 09:09

B.  rng = myRnd(sngVal, 2 ^ 32)

下載 (17.97 KB)

2013-11-1 09:09

略加增修部分：

1. Option Explicit
   
2. 
   
3. Sub Ex()
   
4.     Dim sngVal As Single, rng As Range, xi As Long
   
5.    
   
6.     With Sheets(1)
   
7.         .Cells.Clear
   
8.         
   
9.         For xi = 1 To 5
   
10.             Set rng = .Range(Chr(64 + xi) & "1")
    
11.             For sngVal = 1 To 20
    
12.                 '  rng = myRnd(sngVal)
    
13.                 rng = myRnd(sngVal, 2 ^ 32)
    
14.                 Set rng = rng.Offset(1)
    
15.             Next sngVal
    
16.         Next xi
    
17.         .Columns.AutoFit
    
18.     End With
    
19. End Sub
    
20. 
    
21. Function myRnd(Optional Number, Optional m As Double = 2 ^ 24) As Single
    
22.     Dim tmp As Double
    
23.     '  LCG Algorithm：
    
24.     '      X(n+1) = (a * X(n) + c) mod m
    
25.     Const a As Long = 1140671485
    
26.     Const c As Long = 12820163
    
27.     '  Const m As Long = 16777216    '  2^24
    
28.     Const x0 As Long = 1         '  根據你使用哪組參數，x0有時有需要符合的條件。

複製代碼

stillfish00

回復 15# c_c_lai
其實我的意思是這樣：

1. Sub EX()
   
2.   Dim a As Single, b As Single
   
3.   Dim i
   
4.   
   
5.   a = Rnd
   
6.   For i = 2 To 16777216
   
7.     Rnd
   
8.   Next
   
9.   b = Rnd
   
10.   Debug.Print a, b
    
11.   
    
12.   a = myRnd
    
13.   For i = 2 To 16777216
    
14.     myRnd
    
15.   Next
    
16.   b = myRnd
    
17.   Debug.Print a, b
    
18.   
    
19. End Sub

複製代碼

像rnd一樣，這m次取的亂數都會不一樣，
但rnd取完m次後，再取下一次時，會回到第一次取的值，再一次循環。
就像LCG裡面提到的生命週期

sjgau

比較好的 亂數產生器的 條件
一 產生的計算速度 越快越好

二 亂數的範圍和 個數，越大越好


三 盡可能的 接近均勻分配

四 在產生下一個 相同的亂數之前，要把所有的亂數都跑過一遍


五 循環的週期要 固定

六 循環的 週期要 越大 越好

sunnyso

回復 10# sjgau

不知大大的日本之旅如何？

有空 麻煩指導這個亂數

謝謝

sjgau

回復 18# sunnyso


    11/8 回國
這兩天，在台大有授課的 行程
需要 明天 才能 開始

sjgau

回復 19# sjgau


    http://www.csie.ntu.edu.tw/train/teacher_display.php?num=18
我在 台大ˇ呃 課程介紹網頁