标签：范例 ... 漫步 import print walks np position NumPy

import random
import numpy as np
import matplotlib.pyplot as plt

position = 0
walk = [position]
steps = 100
x = []
y = []
for i in range(steps):
    step = 1 if random.randint(0,1) else -1
    position += step
    x.append(i)
    y.append(position)
    walk.append(position)

plt.plot(x,y)
plt.show()

一次模拟多个随机漫步

import numpy as np

nwalks = 5000
nsteps = 1000
draws = np.random.randint(0,2,size=(nwalks,nsteps))
steps = np.where(draws>0,1,-1)
walks = steps.cumsum(1)
print(walks)
'''
[[  1   0   1 ... -10 -11 -12]
 [  1   0  -1 ...  -4  -5  -6]
 [ -1   0   1 ... -14 -13 -14]
 ...
 [  1   2   3 ...  80  81  82]
 [  1   0   1 ...  38  37  36]
 [  1   0   1 ...  72  73  74]]
'''
print(walks.max()) # 140
print(walks.min()) #-108

hits30 = (np.abs(walks)>=30).any(1)
print(hits30) # [ True False False ...  True  True  True]
print(hits30.sum()) #3408

crossing_times = (np.abs(walks[hits30])>30).argmax(1)
print(crossing_times.mean()) # 算术平均数 486.78372093023256

标签：范例,...,漫步,import,print,walks,np,position,NumPy
来源： https://www.cnblogs.com/nicole-zhang/p/12931204.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。