15-5

最新推荐文章于 2023-03-14 09:29:12 发布

weixin_42167016

最新推荐文章于 2023-03-14 09:29:12 发布

阅读量142

点赞数

本文链接：https://blog.csdn.net/weixin_42167016/article/details/80759276

版权

random_walk.py

from random import choice

class RandomWalk():
    def __init__(self,num_points = 5000):
        self.num_points = num_points

        self.x_values = [0]
        self.y_values = [0]

    def fill_walk(self):
        while len(self.x_values) < self.num_points:

            x_step = self.get_step()
            y_step = self.get_step()

            if x_step == 0 and y_step == 0:
                continue

            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)

    def get_step(self):
        direction = choice([1, -1])
        distance = choice([0, 1, 2, 3, 4, 5, 6, 7, 8])
        step = direction * distance

        return step

rw_visual.py

import matplotlib.pyplot as plt
from random_walk import RandomWalk

while True:
    fw = RandomWalk(50000)
    rw = RandomWalk()
    rw.fill_walk()
    plt.figure(dpi=128,figsize=(10,6))
    point_numbers = list(range(rw.num_points))
    plt.scatter(rw.x_values,rw.y_values,c=point_numbers,cmap=plt.cm.Blues,edgecolors='none',s=15)
    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)

    plt.axes().get_xaxis().set_visible(False)
    plt.axes().get_yaxis().set_visible(False)
    plt.show()

    keep_running = input("Make another walk?(y/n):")
    if keep_running == 'n':
        break