在我看来，您只是随机选择$n$在每个时间步，给定一些分布，你应该抽样$n$从，并使用您的动力学函数从时间传播电压$t$到时间$t + \Delta t$.

然后，您将运行从初始时间到最终时间的多个蒙特卡洛模拟，您基本上选择一个随机值$n$在每个时间步。使用最终结果，您可以计算统计信息，例如某个时间的电压预期值$\hat{t}$, 你的电压在时间的变化$\hat{t}$， 等等。

如果为此编写代码并按照设置的方式使用常量的值，则可能会生成以下类型的图形：

如果$\Delta t$那么足够小$n$只能为 0 或 1。否则，您必须减少$\Delta t$或将您的公式更改为$v(t+\Delta t) = v(t)e^{-\gamma \Delta t} + n \times h_{*}$（顺便说一句，这也适用于$n=0$）。

而如果$n$只能是 0 或 1，没有分布的概念：只有概率。的概率$n$为 1 是$p_0 \times \Delta t$（一阶，所以$\Delta t$必须很小）。

然后只需在 0 和 1 之间（等可能）绘制一个（伪）随机数：如果它小于$p_0 \times \Delta t$然后$n:=1$， 否则$n:=0$. 0到1之间的随机数小于的概率$p_0 \times \Delta t$是...$p_0 \times \Delta t$，这正是您要设置的概率$n$为 1。