在 Julia 中可视化一维随机游走

Question

我想知道如何绘制 1000 个“步行者”，每个“步行者”走 100 步。
我不确定如何绘制包含不止一个步行者的图表，更不用说每次步行 100 步 1000 人了。

下面是一个“步行者”走 10 步的代码。

function walk(N, init::Int=0)
   trace = Int[init]
   for t in 1:N
       if randn() > 0
          push!(trace, trace[end] + 1)
       else
          push!(trace, trace[end] - 1)
       end
   end

   trace
end
walk(10)

[0, 1, 0, 1, 0, -1, 0, 1, 2, 3, 2]

Answer 1

有很多方法可以做到。我会这样做（这不是最有效的方法，但我发现它很容易写）：

using Plots
init = 0
walkers = 1000
walk_length = 100
walks = cumsum(vcat(fill(init, 1, walkers), # initial state
                    rand([-1, 1], walk_length, walkers)), # vertically append move direction
               dims=1) # cumulative sum over the first dimension to get one walk per column
plot(walks, legend=nothing)

在代码中 walks 是一个 Matrix 将每次行走存储在一列中（因此在您的情况下它有 101 行和 1000 列）

在 Julia 中可视化一维随机游走

Visualizing the 1D random walk in Julia

julia

random-walk