我正在实施蒙特卡洛模拟,我需要运行一些动力学的多个实现,然后对所有模拟的最后 state 取平均值。 由于实现的数量很大,我使用 OpenMP 并行运行它们。 每个实现都从相同的初始条件开始,然后在每个时间步以给定的概率发生一个过程,并确定我从均匀分布中抽取随机数的过程。 我想确保所有模拟在统计上都 ...
我正在实施蒙特卡洛模拟,我需要运行一些动力学的多个实现,然后对所有模拟的最后 state 取平均值。 由于实现的数量很大,我使用 OpenMP 并行运行它们。 每个实现都从相同的初始条件开始,然后在每个时间步以给定的概率发生一个过程,并确定我从均匀分布中抽取随机数的过程。 我想确保所有模拟在统计上都 ...
我正在尝试创建一个蒙特卡罗模拟来模拟股票的价格。 股票的价格每天都在变化。 变化由随机变量决定。 天数 (numDays) 内的股票价格被捕获在列表 stock_price_list 中。 我创建了一个数组 monte_list 来存储一堆不同的 stock_price_list。 我想在同一张图 ...
这个问题提醒一些典型的多体问题,但有一些额外的计算。 我正在研究广义 Metropolis Monte-Carlo 算法,用于对大量任意量子系统(例如磁性离子)进行经典交互建模。 但这实际上与问题无关。 有超过 100000 个交互对象,每个对象都可以通过一个坐标和一组描述其当前 state r_i ...
我已经编写了一个代码来将蒙特卡洛模拟的个别结果存储在一张纸中。 宏所做的基本上是在数组中的每次迭代后将 B6:DS6 的值存储在模拟输出 (1) 表中,并最终将它们写入第 6 行下方。 我不是 VBA 方面的专家,但我注意到计算仍然需要相当长的时间,例如 10.000 次迭代需要 15 分钟。 此 ...
正如问题所述,我必须使用蒙特卡洛(随机性)来解决给定的问题。 我正在运行模拟 1,000,000 次。import java.util.*; public class MonteCarlo { public static void main(String[] args) { ...
在 R 中使用模拟,我试图测试选择标准的一致性,但我不断收到以下错误: “ colnames<- ( *tmp* , value = make.names(np)) 中的错误:尝试在小于二维的对象上设置‘colnames’” 我该如何解决这个错误? ...
我是编程的新手,我一直在尝试制作函数(在 C 中)来计算 NxN 站点的方格和三角格中站点i的最近邻居。 至此,以方形格子为例。 我们有 4 个最近的邻居很容易,站点i的前四个最近的邻居可以写为: i+1 ,对于正确的邻居, i-1 ,对于左邻居,对于上邻居和下邻居依此类推。 问题是当我考虑三角 ...
我已经在 Python 中实现了 2D ISING model,使用 NumPy 和 Numba 的 JIT: 当然,哪个有效: 我有两个主要问题: 还有什么可以优化的吗? 我知道 ISING model 很难模拟,但是看看下表,我好像漏掉了什么... 我尝试基于不使用字典一遍又一遍地计算指数来 ...
我正在尝试使用蒙特卡罗模拟来展示当样本维数增加时均匀样本的总和如何呈正态分布。 更准确地说:定义 $X ~ U[2,3]$ 其中 $X_1,...,X_n$ 是来自 X 和 $S = \sum_{1}^{n}(X_i) 的独立同分布样本。 我想使用蒙特卡洛模拟来证明当 n 很大时 S 的分布近似正态 ...
试用号 1 2 3 4 5........ 2000000(两百万) 第 n 次尝试成功 12 4 21 5 10 12 注意:想象一下掷一个骰子,每个结果的概率都是 1/10(不是骰子通常的 1/6)。 对我们来说,“成功”意味着投出“3”。 对于每次试验(见上文),我们不断掷骰子,直到得到“3 ...
我有以下列表: series的mean为 4.74,其np.std等于:3.101 我想从series中生成 1000 个观察结果,所以我使用了以下方法: 问题上面的方法看起来不错,但是它是在series normally的假设下工作的。 目标我的目标是找到一种在不对原始series进行任何假设 ...
这是一个蒙特卡洛模拟问题。 这是我的代码。 接下来的步骤是使用 for 循环并调用模拟 function 来总结最终为 True 的结果数。 我这样做的想法: 但是,我如何确定什么是 True 或 False 结果? 这段代码有意义吗? ...
这就是问题所在:模拟掷两个骰子的平均值。 这是我到目前为止的代码: 接下来的步骤是: 现在,我已经调用了 simulate() function,但我对如何将它添加到我拥有的总变量中感到有点困惑。 ...
这是我并行操作的主要代码: 现在,我需要存储一个关于时间的变量: 请注意,“步骤”和“路径”分别表示时间序列和轨迹总数。 但是,如果我定义这个变量,我会遇到 out of memory 问题,因为 Dim=10000,steps=600,paths=1000,尽管我可以使用多个内核来实现并行操作。 ...
我试图在 x 轴上显示不同模拟次数的蒙特卡洛障碍价格。 到目前为止,这是我尝试过的方法,但出现错误 -> ValueError:x 和 y 必须具有相同的第一维,但形状为 (10,) 和 (5,)。 我是 python 的新手,尽管我很努力,但我找不到错误import numpy as np ...
我发现这段代码可以集成给定的 function 并使用蒙特卡洛方法给出答案。 但是,我想在 python 中实现它,但我不知道如何使“srand(time(Null))”和 Rand_Max 部分发生在 python 中。我还想在这样的情况下使用“func(x)”这样我就可以输入不同的 functi ...
对于一项作业,我需要编写一个程序,该程序可以使用蒙特卡洛方法估算 n 维球体的体积。 我认为我非常接近完成它,但是有一个问题。 我觉得一切正常,除了我需要计算 <1 的sum_x数量的部分。 它似乎并没有将它们全部加在一起。 我怎样才能让它计算所有 <1 的sum_x ? 我该怎么做才能 ...
基本上,我正在尝试创建一个 for 循环,该循环针对不同的样本量计算以下积分的估计值和误差。: 这是单次计算的代码: 但是,我想为 M = 2^i 产生一个错误,其中 i = 1、2、... 10 我的尝试如下所示: 当我运行代码时出现错误: TypeError: only size-1 arr ...
我正在尝试评估以下多维积分: 不可缺少的 答案是 155/6。 但是我使用的 Python 代码: 尝试的代码 不起作用。 我该如何更正或我应该改写什么代码? ...
我正在尝试实施蒙特卡洛方法来求解积分,但无法计算出我的 function 的近似值。 当我替换 0.25 时,近似结果必须是 0.83 - 0.83i(使用在线计算器)但是在我的 C++ 代码中,它的结果是 1。我做错了什么? 近似代码: ...