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

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

import random

import matplotlib.pyplot as plt

repeats, tosses = 60, 8

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

# 返回當(dāng)前平衡參數(shù)p的情況下，8次實(shí)驗(yàn)中正面的次數(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 = [] # 存放每個(gè)p下，60次重復(fù)試驗(yàn)得到的h的值

p = 0 # 初始化p,從0開(kāi)始

while p < 1.01:

hh = []

for t in range(0, repeats):

h = heads(tosses,p)

# 添加微小抖動(dòng)值，防止點(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()

'''

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

# 開(kāi)始沒(méi)做出理想的圖，以為一個(gè)x,對(duì)應(yīng)多個(gè)y沒(méi)法畫(huà)(其實(shí)就是忘記點(diǎn)的抖動(dòng))...

# 這里循環(huán)plot也是可以實(shí)現(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实现简单随机模拟——抛呀抛硬币的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

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