问题标签 [stochastic-process]

For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.

0 投票
1 回答
413 浏览

matlab - 使用 randsample 的随机样本?

我想为 P 和我做一个随机步骤模拟,randsample就像下面这个简单的一样。

%简单的方法

%使用randsample

这是不正确的。你将如何有效地做到这一点?

这就是我正在尝试做的,但我认为这不太对... 在此处输入图像描述

更新与竞争种群 I 和 P 代码基于这组方程。

在此处输入图像描述

0 投票
2 回答
640 浏览

c++ - 使用 QuantLib 的泊松随机变量

好吧,你好,谁能告诉我是否有在 QuantLib 中实现的泊松分布随机变量的随机数生成器?如果是,我在哪里可以找到代码?我正在尝试模拟 Jump-Diffusion 过程并且需要时间步之间的跳转次数(即每个时间间隔 [t_(i-1);t_i[。有没有办法直接在 QuantLib 中执行此操作,还是我需要使用 boost 库?提前致谢!

ps 或者你会建议通过生成指数分布的数字来使用实际的跳转到达时间吗?

0 投票
0 回答
100 浏览

simulation - 比较模拟生成的随机过程的速度?

我有一个基于代理的模拟(使用基于 Java 的 Repast),它在其输出中为我的不同治疗生成时间序列。我正在通过时间测量性能,并且在每次滴答声中,性能都是 30 次运行(30 个样本)的平均值。在所有处理中,性能从接近 0 开始,以 100% 结束,但速度不同。我想知道是否有任何随机模型或概率方法来比较这些时间序列的速度或增长。我想找出哪个“显着”增长得更快。

谢谢。

0 投票
2 回答
271 浏览

matlab - 我们如何使用 Matlab 在以下模型中生成 kappa 和 delta?

我有一个以下随机模型,描述了一个过程(Y)在空间和时间上的演化。Ds 和 Dt 是空间域(具有 x 和 y 轴的 2D)和时间域(具有 t 轴的 1D)。该模型通常称为混合效应模型或变异分量模型

在此处输入图像描述

我目前正在开发 Y 如下:

我的问题是:

  1. 如何在 Y 的表达式中产生 delta 以及 kappa 和 delta 的差异?
  2. 通过使用 Matlab 的一些说明帮助解释我是否正确生成 Y?

如果您需要更多信息/解释,请告诉我。谢谢。

0 投票
1 回答
317 浏览

matlab - 取出维纳过程的随机时间

首先,我是 matlab 和这个论坛的新手,所以请原谅我的无知。

我通过以下方式生成了一个标准的维纳过程(我自己编的,所以如果它是愚蠢的或错误的,我想知道)。

现在我想在一些随机时间找到它的值,比如:

我想(用谷歌和耐心)做类似的事情

然后简单地从 Y 中取出值,但我想找到更直接的方法 - 你有什么建议吗?

0 投票
2 回答
872 浏览

haskell - Haskell 中的随机模拟教程

我想使用 Haskell 进行随机模拟,但我不知道如何。我已经阅读了 Hutton 的“在 Haskell 中编程”,并且我很擅长编写确定性函数程序。但是,我不知道如何开始编写在 R 或 python 等命令式语言中很容易的那种随机模拟。是否有我可以阅读的教程或入门书,或者任何人都可以提供一些入门提示?

0 投票
2 回答
1170 浏览

matlab - 没有均值零的白噪声会做一些改变吗?

我的问题再次与白噪声有关,但含义不同。让我们比较以下两个代码。首先

使用generate(3,500,10)

此代码的图表如下

在此处输入图像描述

但是让我们改变我们的代码,让它在白噪声的情况下为零均值

图表如下 在此处输入图像描述

如果我们比较这两张图片,我们可以说它几乎相同,只是有一些变化,那么我们是否使均值为零是否重要?对于实际分析,例如寻找峰值等。非常感谢

更新:有更新的代码

在此处输入图像描述

这是图片

0 投票
1 回答
3119 浏览

r - 通过 MLE 进行 CIR 模型估计

我想通过 R 中的 ML 来估计 CIR 模型参数。它如下所示:

dr=(theta1-theta2*r) + theta3*sqrt(r)*dW。

该方法在 Iacus 的书“Option Pricing and Estimation of Financial Models with R”附带的 sde packege 中实现。

在那里,在示例(第 5 章)中,实现了速率估计并计算了系数 theta1-3。现在我想用我的数据集(X2)做同样的事情。

我将非常感谢任何帮助!

0 投票
1 回答
142 浏览

numpy - 用numpy按给定概率的随机位数组

我有一个确定性神经网络,我想让它成为随机的。

两个问题:

  1. 我不确定这是否意味着我需要使用 sigmoid 的结果来确定输出的概率,或者概率是否只是神经元输入,而 sigmoid 函数现在是多余的。
  2. 如何使用 numpy 有效地做到这一点?我知道如何生成随机位,但是你如何在一个大数组中使用给定的概率来做到这一点?(如果重要的话,我当前的 sigmoid 函数是 tanh)
0 投票
1 回答
640 浏览

r - Gillespie Stochastic Simulation in Discrete Time using R

I'm simulating a Stochastic Simulation for Epidemiology. How do I simulate it in a discrete time? I managed to obtain for continuous time using the coding below.

How should I alter the coding to get discrete time? Thanking in advance.