steps = np.linspace(0, 100, 101)
x2_ave = np.zeros(101)
x_y0 = np.zeros(101)
x_now = np.zeros(500)
x2_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
x2_now[j] = x_now[j]**2
average2 = sum(x2_now)/500
x2_ave[i+1] = average2
para = np.polyfit(steps, x2_ave,1)
poly = np.poly1d(para)
y_fit = poly(steps)
plt.scatter(steps, x2_ave,s=2)
plt.plot(steps, y_fit, 'r', label = 'fit line')
plt.legend(loc='upper left')
plt.xlim(0,100)
plt.ylim(0,100)
plt.grid(True)
plt.xlabel('step number(= time)')
plt.ylabel('<x^2>')
plt.title('<x^2> of 500 walkers')
plt.show()
代码3
最后编辑于 ：
©著作权归作者所有,转载或内容合作请联系作者
- 文/潘晓璐 我一进店门，熙熙楼的掌柜王于贵愁眉苦脸地迎上来，“玉大人，你说我怎么就摊上这事。” “怎么了？”我有些...
- 文/花漫 我一把揭开白布。 她就那样静静地躺着，像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上，一...
- 文/苍兰香墨 我猛地睁开眼，长吁一口气：“原来是场噩梦啊……” “哼！你这毒妇竟也来了？” 一声冷哼从身侧响起，我...
推荐阅读更多精彩内容
- 本博文原文地址摸我 本篇博文是关于使用注解处理器生成java代码系列的第三篇也是最后一篇文章。在第一篇（在这里）...
- 审查日期：2017-08-09 审查人员：鸣人 代码文件：WebNewsImageShowActivity（点击显...
- （1）浮点型（float）保留小数 （2）在Inspector面板显示的变量，可以使用它来帮你给变量加个描述 未完...
