from __future__ import division
import numpy as np
import matplotlib.pyplot as plt
steps = np.linspace(0, 100, 101)
x_ave = np.zeros(101)
x_y0 = np.zeros(101)
x_now = np.zeros(500)
for i in range(100):
for j in range(500):
ruler = np.random.rand()
if ruler<=0.5:
x_now[j] = x_now[j] + 1
else:
x_now[j] = x_now[j] - 1
average = sum(x_now)/500
x_ave[i+1] = average
plt.scatter(steps, x_ave)
plt.plot(steps, x_y0)
plt.xlim(0,100)
plt.ylim(-1,1)
plt.grid(True)
plt.xlabel('step number(= time)')
plt.ylabel('<x>')
plt.title('<x> of 500 walkers')
plt.show()
代码1
最后编辑于 ：
©著作权归作者所有,转载或内容合作请联系作者
- 文/潘晓璐 我一进店门，熙熙楼的掌柜王于贵愁眉苦脸地迎上来，“玉大人，你说我怎么就摊上这事。” “怎么了？”我有些...
- 文/花漫 我一把揭开白布。 她就那样静静地躺着，像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上，一...
- 文/苍兰香墨 我猛地睁开眼，长吁一口气：“原来是场噩梦啊……” “哼！你这毒妇竟也来了？” 一声冷哼从身侧响起，我...
推荐阅读更多精彩内容
- Pythonic的定义 遵循Pythonic的代码，看起来就像是伪代码，所有的伪代码都可以轻易转换为可执行的Pyt...
- 新版本1.4添加了全选功能添加了不同的变量命名方式可供选择（匈牙利 tv_name、驼峰 tvName、m驼峰 m...
- 今天看见一个博客分享一个很有用的小工具 在线JS转C#代码http://www.m2h.nl/files/js_t...
