-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathrandomWalks.py
65 lines (47 loc) · 1.91 KB
/
randomWalks.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
from random import choice
import matplotlib.pyplot as plt
class RandomWalk():
"""A class that generates random walks"""
def __init__(self, num_points=5000):
"""Initialise attributes for walk"""
self.num_points = num_points
# All walks start at (0,0).
self.x_values = [0]
self.y_values = [0]
def fill_walk(self):
"""Calculate all the points in the walk."""
# Keep taking steps until the walk reaches the desired length
while len(self.x_values) < self.num_points:
# decide which direction to go and how far to go in that direction
x_direction = choice([1, -1])
x_distance = choice([0, 1, 2, 3, 4])
x_step = x_direction * x_distance
y_direction = choice([1, -1])
y_distance = choice([0, 1, 2, 3, 4])
y_step = y_direction * y_distance
# Reject moves that go nowhere
if x_step == 0 and y_step == 0:
continue
# calculate the next x and y values
next_x = self.x_values[-1] + x_step
next_y = self.y_values[-1] + y_step
self.x_values.append(next_x)
self.y_values.append(next_y)
while True:
# Make a random_walk,and plot the points
rw = RandomWalk(50000)
rw.fill_walk()
#set size of the plotting window
plt.figure(figsize=(10,6))
point_numbers=list(range(rw.num_points))
plt.scatter(rw.x_values, rw.y_values, s=1, c=point_numbers, cmap=plt.cm.PuBuGn ,edgecolors='none')
#emphasiz the first and last points
plt.scatter(0,0,c='green',edgecolors='none',s=100)
plt.scatter(rw.x_values[-1],rw.y_values[-1],c='red',edgecolors='none',s=100)
#remove axes
plt.axes().get_xaxis().set_visible(False)
plt.axes().get_yaxis().set_visible(False)
plt.show()
reply=input("Wanna make another random walk?(y/n)")
if reply=='n':
break