123
返回列表 上一主題 發帖

自己寫 自己的方式產生亂數

sunnyso

帖子
605 
主題
92 
精華
0 
積分
648 
點名
0  
作業系統
windows 
軟體版本
7 
閱讀權限
50 
性別
男 
來自
macau 
註冊時間
2013-4-5 
最後登錄
2019-2-10 

熱烈參與勳章好學求知勳章
11# 跳轉到 »
發表於 2013-10-30 12:40 | 只看該作者
回復 10# sjgau

thanks in advance
ss

TOP

sjgau

帖子
145 
主題
45 
精華
0 
積分
190 
點名
0  
作業系統
windows 
軟體版本
7 
閱讀權限
100 
性別
男 
來自
台北市 
註冊時間
2013-10-18 
最後登錄
2014-3-29 
12#
發表於 2013-10-30 13:47 | 只看該作者
回復 11# sunnyso


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

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

所以，必須使用 double,

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

當然，必須 反覆測試 每一個 細節
天下為公

TOP

stillfish00

帖子
1018 
主題
15 
精華
0 
積分
1058 
點名
0  
作業系統
win7 32bit 
軟體版本
Office 2016 64-bit 
閱讀權限
50 
性別
男 
來自
桃園 
註冊時間
2012-5-9 
最後登錄
2022-9-28 
13#
發表於 2013-10-31 14:14 | 只看該作者
回復 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
複製代碼

TOP

stillfish00

帖子
1018 
主題
15 
精華
0 
積分
1058 
點名
0  
作業系統
win7 32bit 
軟體版本
Office 2016 64-bit 
閱讀權限
50 
性別
男 
來自
桃園 
註冊時間
2012-5-9 
最後登錄
2022-9-28 
14#
發表於 2013-11-1 08:16 | 只看該作者
本帖最後由 stillfish00 於 2013-11-1 08:20 編輯

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

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

TOP

c_c_lai

帖子
2035 
主題
24 
精華
0 
積分
2031 
點名
0  
作業系統
Win7 
軟體版本
Office2010 
閱讀權限
100 
性別
男 
註冊時間
2012-3-22 
最後登錄
2024-2-1 
15#
發表於 2013-11-1 09:11 | 只看該作者
回復 14# stillfish00
這是我測試的結果：
A.  rng = myRnd(sngVal)

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

略加增修部分：
  1. Option Explicit

  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

  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有時有需要符合的條件。
複製代碼

TOP

stillfish00

帖子
1018 
主題
15 
精華
0 
積分
1058 
點名
0  
作業系統
win7 32bit 
軟體版本
Office 2016 64-bit 
閱讀權限
50 
性別
男 
來自
桃園 
註冊時間
2012-5-9 
最後登錄
2022-9-28 
16#
發表於 2013-11-1 09:51 | 只看該作者
回復 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裡面提到的生命週期

TOP

sjgau

帖子
145 
主題
45 
精華
0 
積分
190 
點名
0  
作業系統
windows 
軟體版本
7 
閱讀權限
100 
性別
男 
來自
台北市 
註冊時間
2013-10-18 
最後登錄
2014-3-29 
17#
發表於 2013-11-1 13:18 | 只看該作者
比較好的 亂數產生器的 條件
一 產生的計算速度 越快越好

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


三 盡可能的 接近均勻分配

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


五 循環的週期要 固定

六 循環的 週期要 越大 越好
天下為公

TOP

sunnyso

帖子
605 
主題
92 
精華
0 
積分
648 
點名
0  
作業系統
windows 
軟體版本
7 
閱讀權限
50 
性別
男 
來自
macau 
註冊時間
2013-4-5 
最後登錄
2019-2-10 

熱烈參與勳章好學求知勳章
18#
發表於 2013-11-10 09:50 | 只看該作者
回復 10# sjgau

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

有空 麻煩指導這個亂數

謝謝
ss

TOP

sjgau

帖子
145 
主題
45 
精華
0 
積分
190 
點名
0  
作業系統
windows 
軟體版本
7 
閱讀權限
100 
性別
男 
來自
台北市 
註冊時間
2013-10-18 
最後登錄
2014-3-29 
19#
發表於 2013-11-10 12:17 | 只看該作者
回復 18# sunnyso


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

TOP

sjgau

帖子
145 
主題
45 
精華
0 
積分
190 
點名
0  
作業系統
windows 
軟體版本
7 
閱讀權限
100 
性別
男 
來自
台北市 
註冊時間
2013-10-18 
最後登錄
2014-3-29 
20#
發表於 2013-11-10 12:18 | 只看該作者
回復 19# sjgau


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

TOP

        靜思自在 : 要用心，不要操心、煩心。
123
返回列表 上一主題