无标题文章

import random

import numpy as np

import matplotlib.pyplot as plt

import math

def randpath(tt):

b=[]

bb=[]

b.append(0)

bb.append(0)

b2=[]

for i in range(tt):

c=random.random()

bb.append(i+1)

b.append(b[-1]-1+2*c)

for j in range(len(b)):

b2.append(b[j])

return b,bb,b2

a,b,c=randpath(1000)

a1,b1,c1=randpath(1000)

for k in range(len(b)):

plt.scatter([b[k],],[c[k],], 7, color ='red')

plt.scatter([b1[k],],[c1[k],], 7, color ='blue')

plt.xlim(0,1000)

plt.xlabel('time/step number')

plt.ylabel('X')

plt.title('random walk in one dimesion')

plt.show()

"""

c22=[]

c11=[0]*101

for i in range(100000):

a1,b1,c1=randpath(100)

for j in range(len(a1)):

c11[j]=c11[j]+c1[j]/100000.0

for k in range(len(b1)):

plt.scatter([b1[k],],[c11[k],],7,color='blue')

plt.scatter([b1[k],],[a1[k],], 3, color ='orange')

print len(b1),len(c11),len(a1)

z1=np.polyfit(b1, c11,1)

for i in range(len(b1)):

c22.append(z1[0]*b1[i]+z1[1])

plt.plot(b1,c22,color='red',label='y=0.3354*x+0.0102')

plt.xlim(0,100)

plt.ylim(0,)

plt.xlabel('time/step number')

plt.ylabel('the average square distance')

plt.title('random walk in one dimesion')

plt.legend(loc='upper left', frameon=True)

print z1

plt.show()

"""

Contact GitHub
API
Training
Shop
Blog
About

[

](https://github.com/)