周易六十四卦蓍草卦模拟及概率研究

aoe1981

　　又做了一次1亿爻的测试，感觉家里的电脑有点慢，也感觉肯定哪儿错了，频率呈现出了某种稳定……
　　如图：
　　
　　

aoe1981

终于发现了错误，是我的概率计算方法出错了，现在予以修正：

aoe1981

　　概率测试结果如下：
　　
　　关于6、7、8、9（老阴、少阳、少阴、老阳）爻出现的基础概率问题的结帖，附件如下：
　　 六十四卦的概率计算（aoe1981重算）（更正了错误）.rar (24.29 KB, 下载次数: 64)

2014-9-28 19:17 上传

点击文件名下载附件

　　 周易六十四卦（蓍草卦模拟）概率测试版.rar (21.74 KB, 下载次数: 51)

2014-9-29 20:31 上传

点击文件名下载附件

　　后面将着重讨论各变爻出现的概率……

aoe1981

还有一个问题就是彻底搞清六十四卦各卦出现的概率情况究竟是什么样的……
http://club.excelhome.net/forum. ... 863&pid=7088301

fayehoo

非常感谢。太好了。

aoe1981

fayehoo 发表于 2014-9-28 20:28
非常感谢。太好了。

貌似您是用易的人？

aoe1981

　　各变爻概率问题已解决：
　　
　　附件如下：
　　 六十四卦的概率计算.rar (31.82 KB, 下载次数: 32)

2014-9-28 21:53 上传

点击文件名下载附件

aoe1981

代码如下：

1. Option Explicit
   
2. Public Sub bygl() '变爻概率
   
3. Dim i1%, i2%, i3%, i4%, i5%, i6%, gl(6 To 9) As Double, n%, by(0 To 6) As Double, t!
   
4. t = Timer
   
5. gl(6) = 5.17357222844345E-02
   
6. gl(7) = 0.288736860925781
   
7. gl(8) = 0.448361203383929
   
8. gl(9) = 0.211166213405855
   
9. For i1 = 6 To 9
   
10.     For i2 = 6 To 9
    
11.         For i3 = 6 To 9
    
12.             For i4 = 6 To 9
    
13.                 For i5 = 6 To 9
    
14.                     For i6 = 6 To 9
    
15.                         n = 0 '累计变爻数
    
16.                         If i1 = 6 Or i1 = 9 Then n = n + 1
    
17.                         If i2 = 6 Or i2 = 9 Then n = n + 1
    
18.                         If i3 = 6 Or i3 = 9 Then n = n + 1
    
19.                         If i4 = 6 Or i4 = 9 Then n = n + 1
    
20.                         If i5 = 6 Or i5 = 9 Then n = n + 1
    
21.                         If i6 = 6 Or i6 = 9 Then n = n + 1
    
22.                         by(n) = by(n) + gl(i1) * gl(i2) * gl(i3) * gl(i4) * gl(i5) * gl(i6)
    
23.                     Next i6
    
24.                 Next i5
    
25.             Next i4
    
26.         Next i3
    
27.     Next i2
    
28. Next i1
    
29. Sheets("概率研究").Range("b2").Resize(7) = Application.Transpose(by)
    
30. MsgBox "用时：" & Format(Timer - t, "0.0000") & "秒", , "友情提示"
    
31. End Sub
    

复制代码

aoe1981

　　下面是对27楼理论概率的实测支持：
　　
　　附件如下：
　　 周易六十四卦（蓍草卦模拟）概率测试版.rar (21.8 KB, 下载次数: 30)

2014-9-28 22:43 上传

点击文件名下载附件

aoe1981

20000000卦测试代码如下：

1. Public Sub byglcs() '变爻概率测试
   
2. Dim t%, d%, r%, i%, j%, k&, scs%, ys1%, ys2%, js(0 To 6) As Long, n%, ti! '天、地、人、蓍草数、余数、计数
   
3. ti = Timer
   
4. For k = 1 To 20000000 '两千万卦测试
   
5.     n = 0 '累计变爻数
   
6.     For i = 1 To 6 '六爻成一卦
   
7.         scs = 49: r = 0
   
8.         Randomize
   
9.         For j = 1 To 3 '三变成一爻
   
10.             r = r + 1
    
11.             t = Int((scs - 1 - 1) * Rnd()) + 1
    
12.             d = scs - 1 - t
    
13.             ys1 = t Mod 4: If ys1 = 0 Then ys1 = 4
    
14.             ys2 = d Mod 4: If ys2 = 0 Then ys2 = 4
    
15.             t = t - ys1
    
16.             d = d - ys2
    
17.             r = r + ys1 + ys2
    
18.             scs = t + d
    
19.         Next j
    
20.         If scs / 4 = 6 Or scs / 4 = 9 Then n = n + 1
    
21.     Next i
    
22.     js(n) = js(n) + 1
    
23. Next k
    
24. Sheets("概率测试").Range("c11").Resize(7) = Application.Transpose(js)
    
25. MsgBox "用时：" & Format(Timer - ti, "0.0000") & "秒", , "友情提示"
    
26. End Sub

复制代码