python實(shí)現(xiàn)：鏈接如下：個(gè)人覺(jué)得拋硬幣并不是真正的隨機(jī)事件，和拋硬幣時(shí)候的各種狀態(tài)參量有關(guān)系，那么到底什么是真正的隨機(jī)？?www.zhihu.com

# -*- coding: utf-8 -*-

import numpy as np

import matplotlib.pyplot as plt

# 拋硬幣次數(shù)

n_flip = 100

# 硬幣彈跳次數(shù)

n_jump = 4

# 是否均勻初值

is_uniform = True

# 模擬硬幣彈跳

def sim_coin_jump(x):

if is_uniform:

# 均勻初值使用sin(1000x+1)

return np.sin(1000 * x + 1)

else:

# 非均勻初值使用sin(1000πx+1)

return np.sin(1000 * np.pi * x + 1)

# 模擬拋硬幣

def sim_flip_coin(x, n):

for i in range(n):

x = sim_coin_jump(x)

return x

# 初值均勻的隨機(jī)過(guò)程：使用等差序列模擬

def uniform(f):

# 均勻間隔

iterable = (1 + (10 ** -10) * n for n in range(f))

return np.fromiter(iterable, float)

# 生成非均勻間隔序列

def uneven(f):

# 非均勻間隔: 隨著n變化間隔會(huì)變小

iterable = (1 + (10 ** -10) * np.log(n + 1) for n in range(f))

return np.fromiter(iterable, float)

# 統(tǒng)計(jì)并畫(huà)圖

def stat_and_draw(R):

# 統(tǒng)計(jì)正負(fù)數(shù)個(gè)數(shù)

pos = np.where(R > 0)

npos = len(R[pos])

neg = np.where(R < 0)

nneg = len(R[neg])

# 計(jì)算相鄰數(shù)的乘積期望

R1 = R[:-1]

R2 = R[1:]

mean = np.mean(np.multiply(R1, R2))

stat = "Statistics: pos[%d] neg[%d] mean[%s]" % (npos, nneg, mean)

plt.xlabel(stat)

# 標(biāo)題

title = "Flip coin"

para = "(%s and coin_dump[%d])" % ('uniform' if is_uniform else 'uneven', n_jump)

title += para

plt.title(title)

print(title)

print(stat)

# X軸是拋硬幣次數(shù)編號(hào)

X = np.arange(n_flip)

# plt.plot(X, R)

# 紅圓圈標(biāo)記點(diǎn)

plt.plot(X, R, marker='.', markeredgecolor='red', markerfacecolor='red')

plt.show()

def flap_coin():

# 生成初值序列

S = uniform(n_flip) if is_uniform else uneven(n_flip)

# 開(kāi)始模擬拋硬幣

R = sim_flip_coin(S, n_jump)

# 統(tǒng)計(jì)并畫(huà)圖

stat_and_draw(R)

if __name__ == '__main__':

flap_coin()

總結(jié)

以上是生活随笔為你收集整理的python模拟抛硬币_python模拟抛硬币的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

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