技術メモ

役に立てる技術的な何か、時々自分用の覚書。幅広く色々なことに興味があります。

2次元のランダムウォーク

import numpy as np
import matplotlib.pyplot as plt

time = 2000

state = np.zeros(2)
orbit = np.zeros(2)
for t in range(time):
    speed =  1
    randomvec = np.random.random(2) - 0.5
    randomvec *= speed/np.linalg.norm(randomvec)
    state += randomvec
    orbit = np.vstack([orbit, state])

orbit = orbit.T

plt.scatter(0, 0, color="r")
plt.plot(*orbit, lw=0.5)
plt.scatter(*state, color="r")
plt.axis("equal")
plt.show()

f:id:swdrsker:20170421210847p:plain