還是在上次提到的數(shù)據(jù)之魅那本書，看到模擬這章，有個python模擬腳本，但書上不全，就自己簡單寫了下。

流程：在不同的平衡參數(shù)p（為0.5時為均勻的）下，模擬60次實驗，每次投硬幣8次，統(tǒng)計正面朝上的次數(shù)，并作圖。

import random

import matplotlib.pyplot as plt

repeats, tosses = 60, 8

# p為平衡參數(shù)，tosses為每次重復(fù)試驗中投擲硬幣的次數(shù)

# 返回當(dāng)前平衡參數(shù)p的情況下，8次實驗中正面的次數(shù)

def heads(tosses, p):

h = 0

for x in range(0, tosses):

if random.random() < p:

h += 1

return h

x = [] # 存放平衡參數(shù)p

y = [] # 存放每個p下，60次重復(fù)試驗得到的h的值

p = 0 # 初始化p,從0開始

while p < 1.01:

hh = []

for t in range(0, repeats):

h = heads(tosses,p)

# 添加微小抖動值，防止點(diǎn)的重疊

h += (random.random()/4)*random.choice([-1,1])

hh.append(h)

#print(p,'\t',heads(tosses,p))

y.append(hh)

x.append(p)

p += 0.05

# 箭頭標(biāo)注p=0.5，即硬幣平衡的點(diǎn)

plt.plot(x,y,'g^')

plt.annotate('balence', xy=(0.5, 5), xytext=(0.15, 8),

arrowprops=dict(facecolor='black', shrink=0.05),

)

plt.show()

'''

# 一次有意思的嘗試...

# 開始沒做出理想的圖，以為一個x,對應(yīng)多個y沒法畫(其實就是忘記點(diǎn)的抖動)...

# 這里循環(huán)plot也是可以實現(xiàn)的.

for i in range(60):

q = []

for n in range(len(y)):

yn = y[n][i]

q.append(yn)

plt.plot(x,q,'ro')

'''

輸出圖形如下：

可以看到，在p=0.5，即硬幣均勻的情況下，投擲后，硬幣正面出現(xiàn)次數(shù)大部分在4的附近。

總結(jié)

以上是生活随笔為你收集整理的python模拟抛硬币_python实现简单随机模拟——抛呀抛硬币的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。