由搓麻將想到的貝葉斯概率計算【matlab】
本期專欄是由畢導(dǎo)近期相關(guān)視頻所聯(lián)想到的,本UP沒玩過麻將,因此從牌堆剩余牌數(shù)反推畢導(dǎo)視頻中單人擁有一萬的最大牌數(shù)的先驗概率分布的計算推導(dǎo)如下:
視頻鏈接如下:如何成為麻將高手,包你一學(xué)就會,一打就廢!
假設(shè)剩余65張牌,一人一局可擁有13張牌,因此一人有0張(即一張都沒有)一萬的概率為超幾何概率分布。
給出如下定義:

matlab中計算組合數(shù)的函數(shù)是nchoosek(n,k),其中n>k,即從n個數(shù)抽取k個數(shù)的組合數(shù)。
即

至于pi(theta1,2,3,4) 以此類推。
我們按照畢導(dǎo)給出的前提假設(shè),如下:

由此可得邊緣概率分布:

由貝葉斯計算公式,可計算:


以此類推。
下面給出matlab代碼:

主函數(shù)為:

輸出的excel表格示例為:

標(biāo)簽: