cost 73 ms
在 C++ 的并行 PRNG 中取得领先 - Jumping ahead in parallelised PRNGs in C++

我正在实施蒙特卡洛模拟,我需要运行一些动力学的多个实现,然后对所有模拟的最后 state 取平均值。 由于实现的数量很大,我使用 OpenMP 并行运行它们。 每个实现都从相同的初始条件开始,然后在每个时间步以给定的概率发生一个过程,并确定我从均匀分布中抽取随机数的过程。 我想确保所有模拟在统计上都 ...

如何在matplotlib上具体分配Python中的X和Y轴? - How to specifically assign X and Y Axis on matplotlib in Python?

我正在尝试创建一个蒙特卡罗模拟来模拟股票的价格。 股票的价格每天都在变化。 变化由随机变量决定。 天数 (numDays) 内的股票价格被捕获在列表 stock_price_list 中。 我创建了一个数组 monte_list 来存储一堆不同的 stock_price_list。 我想在同一张图 ...

优化蒙特卡洛算法 | 减少 GPU 和特征值问题的操作 | 多体问题 - Optimising Monte-Carlo algorithm | Reduce operation on GPU & Eigenvalues problem | Many-body problem

这个问题提醒一些典型的多体问题,但有一些额外的计算。 我正在研究广义 Metropolis Monte-Carlo 算法,用于对大量任意量子系统(例如磁性离子)进行经典交互建模。 但这实际上与问题无关。 有超过 100000 个交互对象,每个对象都可以通过一个坐标和一组描述其当前 state r_i ...

如何改进 VBA 蒙特卡洛模拟代码? - How to improve VBA Code for Monte Carlo Simulation?

我已经编写了一个代码来将蒙特卡洛模拟的个别结果存储在一张纸中。 宏所做的基本上是在数组中的每次迭代后将 B6:DS6 的值存储在模拟输出 (1) 表中,并最终将它们写入第 6 行下方。 我不是 VBA 方面的专家,但我注意到计算仍然需要相当长的时间,例如 10.000 次迭代需要 15 分钟。 此 ...

计算正方形和三角形网格中所选站点的最近邻居 - Computing the nearest neighbours of a choosen site in both, a square and triangular lattice

我是编程的新手,我一直在尝试制作函数(在 C 中)来计算 NxN 站点的方格和三角格中站点i的最近邻居。 至此,以方形格子为例。 我们有 4 个最近的邻居很容易,站点i的前四个最近的邻居可以写为: i+1 ,对于正确的邻居, i-1 ,对于左邻居,对于上邻居和下邻居依此类推。 问题是当我考虑三角 ...

当 n 增加时用求和证明中心极限定理 - Prove the central limit theorem with a sum when n increases

我正在尝试使用蒙特卡罗模拟来展示当样本维数增加时均匀样本的总和如何呈正态分布。 更准确地说:定义 $X ~ U[2,3]$ 其中 $X_1,...,X_n$ 是来自 X 和 $S = \sum_{1}^{n}(X_i) 的独立同分布样本。 我想使用蒙特卡洛模拟来证明当 n 很大时 S 的分布近似正态 ...

如何使用Excel中的仿真工具解决以下概率相关的问题? - How can I use simulation tool in Excel for solving the following problem related to probability?

试用号 1 2 3 4 5........ 2000000(两百万) 第 n 次尝试成功 12 4 21 5 10 12 注意:想象一下掷一个骰子,每个结果的概率都是 1/10(不是骰子通常的 1/6)。 对我们来说,“成功”意味着投出“3”。 对于每次试验(见上文),我们不断掷骰子,直到得到“3 ...

如何在没有正态性假设的情况下模拟新值? - How to simulate new values without normality assumption?

我有以下列表: series的mean为 4.74,其np.std等于:3.101 我想从series中生成 1000 个观察结果,所以我使用了以下方法: 问题上面的方法看起来不错,但是它是在series normally的假设下工作的。 目标我的目标是找到一种在不对原始series进行任何假设 ...

Monte Carlo Python 问题:模拟三个骰子,总和>10,返回True,否则返回False - Monte Carlo Python Question: Simulation of three dice in which the sum > 10, returns True, otherwise return False

这是一个蒙特卡洛模拟问题。 这是我的代码。 接下来的步骤是使用 for 循环并调用模拟 function 来总结最终为 True 的结果数。 我这样做的想法: 但是,我如何确定什么是 True 或 False 结果? 这段代码有意义吗? ...

julia中的量子蒙特卡洛方法能否并行运算? - Can we do a parallel operation for Quantum Monte Carlo method in julia?

这是我并行操作的主要代码: 现在,我需要存储一个关于时间的变量: 请注意,“步骤”和“路径”分别表示时间序列和轨迹总数。 但是,如果我定义这个变量,我会遇到 out of memory 问题,因为 Dim=10000,steps=600,paths=1000,尽管我可以使用多个内核来实现并行操作。 ...

为 Python 中的期权定价绘制蒙特卡洛模拟 - Plotting Monte Carlo Simulations for option pricing in Python

我试图在 x 轴上显示不同模拟次数的蒙特卡洛障碍价格。 到目前为止,这是我尝试过的方法,但出现错误 -> ValueError:x 和 y 必须具有相同的第一维,但形状为 (10,) 和 (5,)。 我是 python 的新手,尽管我很努力,但我找不到错误import numpy as np ...

Monte-Carlo 方法的实现以找到 python 中的积分值 - Implementation of Monte-Carlo method to find integration value in python

我发现这段代码可以集成给定的 function 并使用蒙特卡洛方法给出答案。 但是,我想在 python 中实现它,但我不知道如何使“srand(time(Null))”和 Rand_Max 部分发生在 python 中。我还想在这样的情况下使用“func(x)”这样我就可以输入不同的 functi ...

在 Python 中使用 Monte Carlo 方法估算球体的体积时无法计算小于 1 的数字 - Cant count up numbers < 1 while estimating the volume of a sphere using the Monte Carlo method in Python

对于一项作业,我需要编写一个程序,该程序可以使用蒙特卡洛方法估算 n 维球体的体积。 我认为我非常接近完成它,但是有一个问题。 我觉得一切正常,除了我需要计算 &lt;1 的sum_x数量的部分。 它似乎并没有将它们全部加在一起。 我怎样才能让它计算所有 &lt;1 的sum_x ? 我该怎么做才能 ...

生成 for 循环以计算具有不同样本量的 montecarlo 积分的误差 - Producing for loop to compute errors for montecarlo integration with different sample sizes

基本上,我正在尝试创建一个 for 循环,该循环针对不同的样本量计算以下积分的估计值和误差。: 这是单次计算的代码: 但是,我想为 M = 2^i 产生一个错误,其中 i = 1、2、... 10 我的尝试如下所示: 当我运行代码时出现错误: TypeError: only size-1 arr ...


 
粤ICP备18138465号  © 2020-2024 STACKOOM.COM