问题标签 [stochastic]

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 投票
2 回答
643 浏览

random - SAT 求解:DPLL 与?

现在我正在写关于解决 SAT 的文章,但我被困在了一个点上。我希望你能帮助我。

我想描述一些解决 SAT 问题的方法。现在我有三种不同的方法:

  1. 蛮力
  2. 随机(天真)
  3. DPLL(具有不同的启发式)
  4. ? 失踪 ?
  5. ...

我的问题是唯一有效的算法是 DPLL(以及其他一些与 DPLL 略有不同的算法)。因此我没有什么可以比较 DPLL 的。

我的问题:如果您能告诉我一些不基于 DPLL (DP) 的算法,我可以将其与之进行比较,那就太好了。

以下是我发现的一些,但无法确定它们是否是一个不错的选择,或者是否有更好的选择:

  • 莫尼恩-斯佩肯迈尔
  • Dantsin、Goerdt、Hirsch 和 Schöning
  • Paturi-Pudlák-Zane-算法
  • 霍夫迈斯特、舍宁、舒勒和渡边

谢谢你的帮助。

0 投票
1 回答
107 浏览

events - How to code a arrival generator with a varying intensity rate

This is for a simulation model:

Most questions I've come about deal with how to code an generator with exponential arrival times.

But I'm currently stuck on how to program a generator where the arrival rate can change within a discrete event simulation.

In particular I'm stuck with the following case: my generator has an input port which accepts an arrival rate (double). If this rate change arrives exactly when an entity is generated, I can simply create the entity, update the rate parameter for the distribution and sample a new arrival time.

But what should I do when the generator at time t1 receives a new rate input event and is already scheduled to create an entity in the future t2 -

Should I a) Abort the creation at t2 and schedule a new creation time using the new rate parameter or b) Just update the rate parameter, let the generator create the entity at t2, and then sample a new arrival time

0 投票
0 回答
33 浏览

time-series - 两种模型中的趋势解读

我估计了以下两个模型:

(注意,yt 是月交易量的对数。)

我如何解释每个人如何为趋势建模?

0 投票
1 回答
294 浏览

python - 用矩阵条目制作直方图?

今天我的任务是制作一个直方图来表示 A^n 的运算,其中 A 是一个矩阵,但仅适用于矩阵中的特定条目。

例如,假设我有一个矩阵,其中行总和为 1。第一个条目是一些特定的十进制数。但是,如果我将该矩阵提高到 2 次方,那么第一个条目就会变成别的东西,如果我将该矩阵提高到 3 次方,它会再次发生变化,等等 - 令人作呕,这就是我需要绘制的。

现在我的尝试是创建一个空列表,然后使用 for 循环将矩阵乘法产生的条目添加到列表中。但是,它所做的只是将最终矩阵乘法的结果打印到列表中,而不是在每次迭代时打印其值。

这是我正在谈论的特定代码:

这是完整的代码,以及我得到的输出。

http://pastebin.com/EkfQX2Hu

0 投票
0 回答
244 浏览

machine-learning - 如何检查或验证 RBM(受限玻尔兹曼机)模型?

我正在尝试实现 RBM,然后我用打网球盒来测试 rbm。

我之前尝试过自动编码器,结果很好。实际上,我对 RBM 它本身的功能感到困惑,我认为它就像自动编码器一样,对输入(每个实例)进行编码以进行特征提取,然后我们可以测试或验证模型(网络)以尝试编码和解码一些实例。

但我面临的问题是 RBM 中某些功能的结果,这似乎很奇怪。

例如 gibbs 采样的结果,采样数据的结果与实际数据非常接近。效果是来自采样数据的 h(x) 的结果,并且实际数据足够接近。

因此,当我尝试将隐藏层中的所有单元解码结果与实际值进行比较时,结果很糟糕,每个特征(单元)的结果几乎相同,大约为 0.4 到 0.5。

然后 f(x) = 1/m*sigma(log(p(x))) 丢失的功能它本身只是大约 0.07142857142857142,它永远不会改变(变化大约是 0.00000000000000001 或 0.00000000000000002)。

我对每个特征都使用了 continue 值,使用标准归一化,因此输入范围值为 0 到 1。

有人有建议吗?

*对不起我的英语不好:D

0 投票
0 回答
212 浏览

matlab - Matlab中使用欧拉法的随机微分方程

我正在尝试使用“sde_euler”求解两个随机微分方程组。这是我的代码:

指定 t0,tf,x0,F0 的值后,我运行代码并遇到以下错误:

我不明白我的错误。我的 TSPAN 长度大于 2。以下是工具箱的来源:

https://github.com/horchler/SDETools/blob/master/SDETools/sde_euler.m

请问有什么帮助吗?谢谢你!

0 投票
2 回答
11412 浏览

r - R:如何生成嘈杂的正弦函数

我对整个 R 事物还是很陌生。

我有以下目标;我有一个正弦函数,可以描述随时间变化的钙粒子数:类似于 y = a * sin (b*t) + c

由于实际上钙的产生和去除是在随机事件中描述的,我想在我的函数中添加一个随机噪声项(最好在平均噪声幅度上进行扩展)。

像 z = y + random*Amplitude

你能帮我吗?

最好的

0 投票
0 回答
209 浏览

python-3.x - 几何分布中的随机数,使得它们的总和等于 SUM

我想从指数形状的值分布中绘制 k 个随机数 i_1,...,i_k,min <= i <= max,m,std 是人口值的中值和标准。sum(i1,..,ik) 应该等于给定参数 SUM。

示例:
给定

k = 9
总和 = 175
最小值 = 8
最大值 = 40
m = 14

所需:[9、10、11、12、14、17、23、30、39]

我不知道这是否真的可能,而不取决于运气来绘制满足 SUM 规则的组合。我将不胜感激任何帮助或评论。谢谢你。

编辑:在以前的版本中,我写了关于不可能精确解的指数分布,而是指具有离散值的指数分布,例如几何分布。

EDIT2:更正了示例中的数字 k。

0 投票
2 回答
171 浏览

algorithm - 实现随机 ACO 算法

我正在尝试实现一个随机蚁群优化算法,但我无法确定如何根据概率实现运动选择。

到目前为止,我已经实现的标准(贪婪)版本是,在图的顶点处的蚂蚁是m边集,将根据以下标准选择下一个顶点:iG = (V,E)E(i, j)j

我遇到的问题是尝试实现它的随机版本,因此现在选择新顶点的标准j是:

我了解它背后的数学原理,我只是无法弄清楚我应该如何实际实现它。

如果说,我有 3 个顶点连接到i,每个顶点的概率为 0.2、0.3、0.5 - 进行选择的最佳方法是什么?我应该只是随机选择一个顶点,然后在 (0,1) 范围内j生成一个随机数,如果,选择顶点?或者,还有更好的方法?rr >= P(j)j

0 投票
1 回答
499 浏览

python - 随机 python 模拟的校准

我在 python 中有一个(基于代理的)模拟。结果是随机的。我想校准模拟以最小化平方距离的总和。有哪些 Python 算法可用于最小化随机模拟?