高分子链体系的蒙特卡罗模拟
是一种计算方法，用于模拟高分子链的行为和性质。它基于蒙特卡罗方法，通过随机抽样和统计分析来模拟高分子链的构象、动力学和热力学性质。
高分子链体系的蒙特卡罗模拟可以帮助研究人员深入了解高分子材料的结构和性能，以及它们在不同条件下的行为。通过模拟，可以预测高分子链的构象、链的运动方式、相变行为等重要信息，为高分子材料的设计和应用提供指导。
在实际应用中，高分子链体系的蒙特卡罗模拟可以用于研究聚合物的力学性质、热学性质、输运性质等。例如，可以模拟聚合物的拉伸行为，预测其力学性能；可以模拟聚合物的熔融过程，研究其热学性质；可以模拟聚合物在溶液中的扩散行为，研究其输运性质。
腾讯云提供了一系列与高性能计算相关的产品和服务，可以支持高分子链体系的蒙特卡罗模拟。其中，推荐的产品包括：
- 弹性计算服务（Elastic Compute Service，ECS）：提供高性能的计算实例，可用于运行蒙特卡罗模拟的计算任务。 产品介绍链接：https://cloud.tencent.com/product/ecs
- 弹性高性能计算（Elastic High-Performance Computing，EHPC）：提供高性能计算集群，可用于并行计算和大规模蒙特卡罗模拟。 产品介绍链接：https://cloud.tencent.com/product/ehpc
- 云硬盘（Cloud Block Storage，CBS）：提供高性能的块存储服务，可用于存储蒙特卡罗模拟的输入数据和输出结果。 产品介绍链接：https://cloud.tencent.com/product/cbs
- 对象存储（Cloud Object Storage，COS）：提供高可靠性、高可扩展性的对象存储服务，可用于存储大规模的蒙特卡罗模拟数据。 产品介绍链接：https://cloud.tencent.com/product/cos
通过腾讯云的产品和服务，研究人员可以方便地进行高分子链体系的蒙特卡罗模拟，并获得准确、高效的模拟结果。
相关·内容
2回答
我想开发风险板游戏，这将包括一个电脑玩家的人工智能。Moreovor，我读了两篇关于它的文章，和，我意识到我必须学习蒙特卡罗模拟和马尔可夫链技术。我认为我必须把这些技术结合起来，但我想它们是不同的技术，可以用来计算关于过渡态的概率。
那么，有谁能解释一下它们之间的重要区别和优缺点呢？
最后，如果你要为风险游戏实现人工智能，你会更喜欢哪种方式？
你可以找到简单的确定概率的战斗的结果在风险棋盘游戏，和蛮力算法使用。有一个树形图，它指定所有可能的状态。我应该在这棵树上使用蒙特卡罗还是马尔可夫链？
浏览 4提问于2013-05-07得票数 8
回答已采纳
NSTableView声明了这个可选的委托方法：
- (CGFloat)tableView:(NSTableView *)tableView sizeToFitWidthOfColumn:(NSInteger)column;
实现它很容易，但是当您的表中填充了数千个条目时，您就会看到性能问题，调用这个方法会使应用程序失去响应性。
以下是文档中关于实现该方法的说明：
如果要控制当用户双击调整大小分配器时列的宽度，请实现此方法。默认情况下，NSTableView迭代表中的每一行，通过preparedCellAtRow:column:访问一个单元格，并请求"cellSize“找到要使用的适
浏览 1提问于2014-10-14得票数 1
回答已采纳
7回答
使用带有GPU的随机数
、、、
我正在研究使用nvidia GPU进行蒙特卡罗模拟。但是，我想使用gsl随机数生成器和并行随机数生成器(如SPRNG )。有人知道这是否可能吗？
更新
我用GPU和RNG玩过。目前还没有一个很好的解决办法。随SDK而来的Mersenne并不适合(我的)蒙特卡罗模拟，因为它需要相当长的时间才能产生种子。
NAG图书馆更有希望。您可以在批处理或单个线程中生成in。然而，目前只支持少数几种分布--均匀、指数和正态分布。
浏览 11提问于2010-06-26得票数 9
回答已采纳
2回答
我很难用python 3创建我自己的蒙特卡罗模拟器..。我正试图预测未来一年的价格。我仍然在学习如何在python中编程。
我的数据：
Price = np.array([29429,30426,32513,40605,52806,57581])
Year = np.array([1970,1971,1972,1973,1974,1975])
尝试蒙特卡罗模拟器..。
N_pts_per_fit = 30
indices= random.randint(0,Price,N_pts_per_fit)
for i in range[N_pts_per_fit]:
newdata = in
浏览 10提问于2016-11-30得票数 0
我试着测量函数使用的时间，比如：
start = time.time()
function("does_something")
end = time.time()
function_time = end - start
然而，每次我运行这个程序时，我都会得到一个(稍微)不同的解决方案。原因是什么？我知道我可以做一个蒙特卡罗模拟来使这个测量更可靠，但是我很想知道为什么结果会有这种差异。
浏览 2提问于2022-06-04得票数 2
回答已采纳
2回答
在for循环中启动多进程应用程序时，Python内存不足。每增加一个循环，分配的内存就会变得更大。如何解决此问题？ 我正在执行大型蒙特卡洛模拟，每个模拟都有数千个模拟。为了获得性能提升，我使用了multiprocessing模块，并在10个内核上并行运行各个模拟。每个蒙特卡洛模拟使用相同的模型，但不同的模型输入。我基本上是在循环一个输入列表，并在上一个蒙特卡罗模拟完成后，用下一个模型的输入开始一个新的蒙特卡罗模拟。 奇怪的是，之前的蒙特卡洛模拟中分配的内存并没有在它完成后释放。每增加一次蒙特卡洛模拟，分配的内存就会越来越大，直到python耗尽内存。上一次蒙特卡洛模拟中的每个python对象
浏览 44提问于2019-02-06得票数 0
1回答
阶小于E-1的随机数
、、、
我目前正在用Python进行动态蒙特卡罗模拟，所以我需要生成在0到1之间的伪随机数，它们遵循均匀分布，但是我正在搜索一些小于E-1的数字，例如，0.001。我知道这并不常见，因为它们毕竟是随机数，但我的问题开始于蒙特卡罗的实现。
简单地说，蒙特卡罗模拟的问题是，我有一个constraint_probability=c和一个伪随机的number=r，但是由于问题的动态性，我的代码只能在r< c的情况下运行。如果我做了c=0.001，我的代码就永远不会运行。
我的第一个想法是使用random.uniform，但我发现得到的数字至少是E-2的。
后来，我认为可以通过设置以下函数来解决我的问题：
浏览 2提问于2021-04-23得票数 0
在这种情况下，有人能帮我吗?我想用蒙特卡罗模拟生成一系列服从MATLAB中给定分布的随机数。
例如:指数分布p(x)=exp(-x) x>0
谢谢。
浏览 4提问于2011-12-04得票数 1
2回答
我正在尝试将一些python代码翻译成C++。代码所做的是运行蒙特卡罗模拟。我认为Python和C++的结果可能非常接近，但似乎发生了一些有趣的事情。
下面是我在Python中所做的工作：
self.__length = 100
self.__monte_carlo_array=np.random.uniform(0.0, 1.0, self.__length)
下面是我在C++中所做的工作：
int length = 100;
std::random_device rd;
std::mt19937_64 mt(rd());
std::uniform_real_distribution<
浏览 11提问于2013-08-12得票数 9
回答已采纳
我有一项任务是使用蒙特卡罗方法来评估一个不公平的硬币翻转，并确定在n个模拟中从n个翻转中得到n个头部的概率。我的教授提供的指南(底部)对于如何根据不公平的硬币调整它是没有意义的。我已经包括了我开始编码的内容，但它还远未完成。我怎样才能调整蒙特卡罗来为一枚不公平的硬币工作呢？
问题：
假设你有一枚硬币是不公平的，有一个头的概率是p。我们感兴趣的是，在总共的nflips翻转中，能看到确切的人头的概率。编写一个函数，将这些参数用于在nsim模拟中计算这个概率(总共有四个参数)，并在几个不同的值上测试它。提示:在示例函数中使用prob参数。
nsim <- 10000
nheads <-
浏览 4提问于2022-04-21得票数 1
回答已采纳
3回答
我对n个样本做了蒙特卡罗模拟。对于每一个样本，我需要计算值Xi，因此，我将得到的结果如下：
X= X1，X2，.，Xn
(在这里，Xi可以是矩阵或数字)。
现在我要计算这些样本的平均值，我称之为Xmean。所以我需要得到这样的东西：
X1，(X1+X2)/2，(X1+X2+X3)/3 .，(X1+X2+...+Xn)/n
在Python中，我编写了一个代码：
for i in range(N):
for j in range(i+1):
Xmean(i) = Xmean(i) + X(j)
Xmean(i) = Xmean(i) / (i+1)
浏览 5提问于2014-01-18得票数 0
我正在使用Matlab脚本对Simulink模型运行蒙特卡洛模拟，该脚本大致如下所示：
model = 'modelName';
load_system(model)
for ii = 1 : numberOfMC
% Some set_param...
% Some values are set
sim(model);
results{ii, 1} = numberOfMC;
% ect...
end
close_system(model,0);
随着蒙特卡罗试验次数的增加，一次模拟的时间也会增加，如n^2。
有没有一个简单的解释
浏览 0提问于2013-12-04得票数 6
1回答
我试图在卡方test.In R中绘制零分布，用以下代码进行蒙特卡罗模拟得到经验p值是可行的：
chisq.test(d,simulate.p.value=TRUE,B=10000)
但是它不返回分布图。有没有办法使R返回测试的模拟值？
浏览 5提问于2016-04-21得票数 3
回答已采纳
如果我只想计算一个标准的正常变量，那么使用：
samples=norminv(rand(N),0,1);
和
samples=normrnd(0,1,N,1);
无论是在处理时间上还是在蒙特卡罗模拟中使用时的收敛性？
我问这个问题的原因是，我想使用类似于norminv()的Halton数这样的准蒙特卡罗技术来复制normrnd()，但首先我要确保我理解它们之间的关系。
我想中心问题之一是:在rand()和normrnd中，随机数是如何产生的？是同一种方法吗？
如果它们完全等同，为什么重复呢？
浏览 0提问于2014-01-23得票数 2
回答已采纳
1回答
对于BPSK的误码率和信噪比分析，我使用MATLAB中的Simulink和BERtool。下面显示的是Simulink模型：
在AWGN组件中，参数设置如下：
然后我打开蒙特卡洛的BERtool ->。此处，Eb/No的取值范围为0:1:10
它不会为每次迭代更改AWGN组件中的Eb/No值。Eb/No值在蒙特卡罗模拟中的确切作用是什么？
浏览 43提问于2017-02-18得票数 0
1回答
我正在我的4核PC上使用GNU Octave 4.0.0计算蒙特卡罗模拟.模拟大约需要4个小时来计算脚本50,000次(特定于我的问题)，这是一个很长的计算时间。我想知道是否有办法同时在多个核上运行Octave，以减少计算时间。
提前谢谢。
浏览 0提问于2016-05-18得票数 3
回答已采纳
我不确定我是否应该在交叉验证或这里打开这个问题。但是由于这个问题属于与数据科学相关的一个更大的项目，所以我选择了这个。
我将介绍我的工作项目的一个简化版本，因为原始版本过于复杂，而且是特定的领域。
假设我们有1小时(60分钟)的时间线。在此期间，作业开始运行，并以随机点创建用户通知。我已经写了一个蒙特卡罗模拟来研究这个过程。
我的主要问题是：
这些通知在60分钟内传播得有多好？
我们是否有这个时间线中不包含通知的部分，其余的是否聚集在特定的时间内？
在执行中改变随机函数，会对上述答案产生怎样的影响？
蒙特卡罗模拟实际代码的伪代码是：
Repeat one million times
n
浏览 0提问于2017-01-31得票数 4
回答已采纳
1回答
(Monte粒子)是一种用蒙特卡罗方法模拟辐射效应的古老而著名的Fortran程序。可用于辐射安全计算。
要获得MCNP套件是非常困难和昂贵的，而且代码很难阅读。
什么是MCNP的开源替代方案？
浏览 2提问于2013-10-08得票数 3
回答已采纳
我计划做一个蒙特卡罗模拟，但是AnyLogic页面上的介绍并没有真正帮助我。我希望通过统一(0，1)更改一个参数，并运行模型1000次，并将结果保存在excel表中。我的想法是有一个名为"durchlaufzeit“的变量，它需要模型在水槽上花费时间。所以在水槽里我说durchlaufzeit = Math.round(time())。这是获得总延迟时间的正确方法吗？或者有什么可以替代的呢？
接收器中的命令如下所示：
因此，在离散事件仿真结束时，我想将仿真的整个延迟时间写入我的ExcelSheet ergebnisse中。我用以下方式创建了蒙特卡罗实验：
在我看来，
浏览 7提问于2022-03-14得票数 -2
回答已采纳
1回答
我正在寻找一种方法，用蒙特卡罗模拟从正态分布中生成两个相关变量。特别是，我想定义这两个变量之间的不同相关性(即R = .30，.60，.90)。非常感谢!
浏览 2提问于2021-02-22得票数 0
回答已采纳
我试着做一个时间序列蒙特卡罗模拟的可靠性分析。我的系统的故障率/年.The是一个常数，所以我从指数分布的exprnd(1/λ)生成随机数，计算出未来可能的故障次数(年份)，-k.The蒙特卡罗模拟周期将是30年，试验的次数将由用户定义，每次迭代的结果将存储在矩阵中，每列都将是我的系统在其30年生命中的故障次数，cycle.The问题是，我不知道loop.The代码在一个30年内的确切迭代次数将与以下类似。
MTTF=0.6; % enter mean time to failure 1/λ
A=zeros(300,1); %preallocate a large enough
浏览 2提问于2014-01-08得票数 0
回答已采纳
2回答
我有一个关系式:y=a+b+c
我有a，b和c的平均值和标准差，我想通过蒙特卡罗模拟得到y的概率分布。
有没有什么函数、包或者简单的方法可以让我这样做呢？
浏览 0提问于2012-06-29得票数 1
回答已采纳
3回答
我最近正在研究一些关于优化的算法。我发现了一种非常有趣的算法，叫做模拟退火。根据我的理解，我们的初始点变成了一个分布，而不是一个实际的点。它也与初始“温度”有关，但其主要思想是“降温”。这是当你到达局部最优时，尝试某些情况下的“偏移”。从而克服局部最优的陷阱。我想知道的是模拟退火法也是蒙特卡罗方法的一种。不同之处在于蒙特卡洛算法是全局的，而模拟退火算法是局部的。我的想法正确吗？
浏览 1提问于2019-05-15得票数 1
1回答
从二项分布生成随机数
、、、、
在我的蒙特卡罗模拟中，关键部分是从参数n= size，p= 0.5的二项分布中生成一个随机数。以下是我当前的实现
#include <stdio.h>
#include <stdint.h>
#include <x86intrin.h>
int64_t rbinom(int64_t size) {
if (!size) {
return 0;
}
int64_t result = 0;
while (size >= 64) {
uint64_t random64;
wh
浏览 0提问于2019-03-23得票数 3
我正试图对系统gmm估计的模型进行蒙特卡罗模拟。因此，我需要从python ()中的pydynpd包中提取模型的系数。我正在搜索一个命令/函数，它与使用fit().params返回状态模型一样，也就是数组中的系数。
浏览 19提问于2022-05-04得票数 0
1回答
我试图改变我的Y轴上的直方图，我已经创建了一个蒙特卡罗模拟。Y轴目前正在输出与百分比相反的事件数.下面是我的代码，我假设需要有一个函数"For Y in Range"，但我不确定实现它的确切方法。
浏览 15提问于2022-01-12得票数 0
2回答
我创建了一个，它使用蒙特卡罗模拟来计算π，目前有3202959700000‘点’和2515023487209在圆圈中。将其插入公式4 * 2515023487209 / 3202959700000返回3.1408743446993728956377440527897，这是非常不正确的。这是错的，因为这里没有足够的“点”，还是我只是把它们放错了地方？
浏览 0提问于2018-06-15得票数 0
回答已采纳
1回答
我有三个变量，假设a，b和c，我用正态分布产生了10000个随机数。现在，我想运行蒙特卡罗模拟，以得到1000次的结果。
我的模型是，Y=0.5*a+0.4*b+0.6*c。
有谁能帮忙处理matlab的代码吗？
提前感谢
浏览 4提问于2012-05-14得票数 0
1回答
我想用蒙特卡罗方法模拟光脉冲在云层中的传播和它的后向散射激光脉冲波形。我能够模拟传播，但我的问题是如何模拟激光脉冲的检测波形。激光和探测器位于云环境中的$z=0$，我的问题是如何对探测器探测到的光子的轨迹进行编程。值得一提的是，在模拟中使用了Henyey-Greenstein散射函数。在此之前，非常感谢您。
浏览 2提问于2020-01-15得票数 0
1回答
在游戏中有5个玩家，每个人手上有8张牌。没有其他卡片存在。
使用蒙特卡洛(MCTS)，我运行了游戏从当前状态到比赛结束的许多模拟。在这些模拟中，对手使用非常简单的愚蠢规则。当生成一个新的模拟时，我给它们分配了一个均匀随机分布的未知卡片。我说如果一张卡片不在我手里，它不在桌子上，它以前也没有被丢弃，那么它一定是在别人的手里。
然而，均匀的随机分布并不是分发卡片的最佳方式。事实上，有时候对人类玩家来说，这样的推理是很简单的：
如果他有3♥，他可能会在前一轮放弃它，因为这是一个非常好的移动。那么，他不太可能在下面丢弃3♥。
那么，在放弃之前，在生成蒙特卡罗模拟之前，我会给每个对手分配一组未知的卡片及
浏览 0提问于2014-01-09得票数 3
我试图实现蒙特卡罗模拟方法来估计R中的积分，但我还是得到了错误的答案。我的代码如下：
f <- function(x){
((cos(x))/x)*exp(log(x)-3)^3
}
t <- integrate(f,0,1)
n <- 10000 #Assume we conduct 10000 simulations
int_gral <- Monte_Car(n)
int_gral
浏览 6提问于2020-04-21得票数 0
回答已采纳
我需要为一些金融交易创建一个蒙特卡罗模拟器。这些投入将是：
最终盈利的交易的平均百分比。
每笔交易的平均利润
每个时间段的交易数量
我已经看过 Perl模块，但不确定如何制定模拟器的输入。
考虑到我正在工作的输入，有人能给我一些关于开始的建议吗？
浏览 1提问于2010-09-29得票数 4
1回答
我正在使用蒙特卡罗模拟来计算失败的概率，我想在散射直方图中画出代表故障的点(红色)。我能.我可以吗？代码如下
% Histograms
hist(S,20)
hold on
hist(R,40)
set(findobj('Type','patch'),'Facecolor','none','Edgecolor','black')
set(gca,'Fontsize',18,'Fontname','euclid')
xlabel('R
浏览 5提问于2014-10-14得票数 0
回答已采纳
1回答
我正在寻找一本用R编写参数化多变量蒙特卡罗模拟程序的指南。
我花了很多时间寻找、阅读和研究R中的多变量模拟指南。然而，它们都是从一个数据样本开始的，并通过从样本中随机选择来执行模拟。
我还看到了用于编程参数单变量模拟的方法。
这两个都不能回答我的用例:我有一个没有实际样本的多变量问题。相反，我有几个用它们的平均值和标准差描述的正态变量，并且我有一个变量的相关矩阵。我想做的是直接根据参数编写一个蒙特卡洛模拟程序。有谁能给我指个方向吗？
谢谢。
浏览 4提问于2019-09-29得票数 0
Dask分布式支持，它可以加快计算速度，提高计算的鲁棒性，但是每个任务可以运行不止一次。
在这里，我要求一个方法来“整理”工人的结果，他们并没有对最终的结果作出贡献。为了说明我所要求的：
让我们假设，每个工作人员都在进行蒙特卡罗模拟，并将~10 in的模拟结果保存在一个结果文件夹中。当工作被窃取时，仿真结果会被存储几次，因此只保留其中的一个就好了。实现这一目标的最佳途径是什么？dask.distributed能自动调用一些“清理”过程来处理那些最终不会对最终结果做出贡献的任务吗？
编辑:我目前使用以下代码启动模拟：
c = distributed.Client(myserver)
mytask
浏览 1提问于2017-02-01得票数 0
回答已采纳
2回答
我在VBA中做了蒙特卡罗模拟。客户端希望(不问为什么)修正随机数序列，即每次运行模型时，序列将保持不变。我设法修复了描述的随机种子。但在不同的个人电脑上却不一样。知道为什么以及如何在不同的机器上修复它吗？
浏览 8提问于2014-04-07得票数 2
回答已采纳
在中，根据蒙特卡罗政策评估，他在第111页提到了note that the computational expense of estimating the value of a single state is independent of the number of states。然而，自蒙特卡罗以来：
状态的平均回报是从第一次遇到状态到插曲结束时计算的。
州越多，事件结束的时间越长。
关于那份声明我遗漏了什么？
浏览 4提问于2017-05-09得票数 0
回答已采纳
1回答
我想知道是否有一种更有效的方法来加速下面的代码。这个函数的意思是使用Mersenne进行蒙特卡罗模拟，用剩馀的牌填充一组扑克牌。
void HandEvaluator::RandomFill(std::vector<std::shared_ptr<Card>>& _Set, std::vector<std::shared_ptr<Card>>& _Dead, unsigned int _Target){
//Add the cards that are currently in Set as dead cards
浏览 0提问于2018-12-01得票数 3
都市光交通是大都会黑斯廷斯算法的一个应用.其变异只在使用的“突变策略”中有所不同。
是否正在研究其他的马尔可夫链蒙特卡罗算法(具有更快的收敛速度)用于射线追踪？
如果是的话，我对任何报纸都感兴趣。
浏览 0提问于2018-08-26得票数 3
我试图计算我的估计量的蒙特卡罗标准误差。我首先根据蒙特卡罗模拟的重复得到了我的估计量：
library(VGAM)
# Log-likelihood function of Gumbel distribution
ll_fn <- function(par, z, m){
mu <- par[1]
sigma <- par[2]
-m * log(sigma) - sum((z - mu)/sigma) - sum(exp(-((z - mu)/sigma)))
}
set.seed(101)
N <- 500
#sample size 20
MLE &
浏览 2提问于2022-11-17得票数 0
回答已采纳
1回答
B <- 10000
results <- replicate(B, {
hand <- sample(hands1, 2)
(hand[1] %in% aces & hand[2] %in% facecard) | (hand[2] %in% aces & hand[1] %in% facecard)
})
mean(results)
这段代码工作得很完美，并完成了所需的操作。
这是蒙特卡罗模拟。我不明白他们把花括号{}放在复制函数中的方式。我可以理解这些代码的功能，但我不能理解他们将代码放入的方式。
浏览 3提问于2021-08-01得票数 2
回答已采纳
1回答
使用loess(y~x)和predict，我在R中找到了一个随机变量的黄土估计。现在，我想进行100次蒙特卡罗模拟来找出标准误差。我使用了下面的R代码，但它在所有100个模拟中都返回了相同的值。请告诉我哪里出错了。
set.seed(6)
pn = replicate(100,predict(fit,data=data.frame(X=78)),simplify=FALSE)
plot(pn)
sd(pn)
浏览 0提问于2020-08-24得票数 0
2回答
有谁知道如何最好地实现一个随机算法(蒙特卡罗模拟来逼近一个值)在一个固定的时间(5分钟)运行？也就是说，算法将继续采集新的样本，直到执行时间在5分钟左右？
该算法可以在一定的时间内比较不同的随机方法。
谢谢。
浏览 0提问于2015-06-30得票数 0
回答已采纳
1回答
我想从MATLAB的代码定义如下的均匀分布中随机抽取:-
pd1 = makedist('Uniform','lower',-0.0319,'upper',0.0319); % X1
在MATLAB中，通常的命令是random()，但是帮助文件告诉我它只用于Guassian混合分布。因此，它也可以推广到均匀分布中，或者在蒙特卡罗模拟中是否有其他函数可以显式地随机绘制。
浏览 6提问于2020-11-27得票数 0
回答已采纳
我正在为蒙特卡罗模拟编写一个c++代码。因此，我需要生成许多均匀分布在[0，1]之间的数字。我包括了从获取的生成我的数字的以下代码：
// uniform_real_distribution
#include <iostream>
#include <random>
std::default_random_engine generator;
std::uniform_real_distribution<double> distribution(0.0,1.0);
int main()
{
double number = distributio
浏览 4提问于2015-05-28得票数 2
回答已采纳
1回答
如何在多进程中共享对象
、、、、
感谢你阅读这篇文章。我已经开发了国际象棋AI，就像李世石或AlphaGo AlphaGo。我用过Python和tensorflow。国际象棋人工智能由蒙特卡罗树搜索、策略网络和价值网络组成。
我为蒙特卡罗树搜索的策略和价值网络做了学习。没有问题。但是，Montecarlo树搜索中的每个模拟都太慢了。所以我想提高每个模拟的速度。
我知道因为GIL，python不能共享对象。我真的需要帮上忙。如果你们有在python多进程中共享对象的经验，请分享你们的经验。
我在这个页面下面发布了摘要代码。
附言:我的英语不好。所以，如果你在阅读这一页时感到不舒服，那是我的错。请理解这一点。
class mont
浏览 6提问于2018-03-31得票数 1
1回答
我想添加完整的蒙特卡罗体积散射到我的路径追踪器，但我有一个困难的时间来研究如何做。让我解释一下我想做什么：
📷
一条射线进入一种物质，我们应用BTDF，然后经过一段距离后，发生体积散射事件，之后(在各向同性情况下)，射线在球体中的任意方向散射。这会重复，直到射线用另一个BTDF离开物质为止。
我的问题如下：
如何选择散射事件之间的距离？直觉告诉我，应该有某种散射pdf，它给出了一定距离后散射的概率？。
这是对的吗？
pdf是否是各向同性材料的线性函数？
这个函数有名字或者我可以搜索的东西吗？
Beer还会在散射事件之间应用吗？
我可不这么想。因为Beer是对实际散射计算的简化。
再说一遍，也
浏览 0提问于2016-06-04得票数 9
回答已采纳
我在我的R代码中使用future.apply软件包对一个优化问题进行100蒙特卡罗模拟的并行处理。
现在，让我们假设每个模拟的优化输出都是True (找到了最优解)或False (没有找到最优解)。我是否可以设置条件，以便如果前10次模拟返回False，则并行过程将被终止？
编辑:，我在下面包含了一个小例子--我可以在运行10次之后终止它吗？
library(future.apply)
plan(multiprocess)
outputs <- future_lapply(1:100, function(i){
temp <- FALSE
return(temp)
})
浏览 5提问于2020-08-12得票数 3
我正在运行一个蒙特卡罗实验，使用“随机种子”选项进行几次迭代，以进行独特的模拟运行。看起来，在其中一些运行中会发生奇怪的事情，这是由输出数据所指示的。
我想收集每一次蒙特卡洛运行的种子值，以便使用“固定种子”选项来重现实验后的奇怪运行。因为我没有找到一个‘getSeedValue’-函数-我怎么能这样做？
我已经在：中尝试了建议的解决方案，但是得到了一个空指针异常。
浏览 4提问于2021-09-17得票数 2
回答已采纳
从概念上讲，我理解像蒙特卡罗这样的数值方法是如何用来求解一个确定的积分的。由于函数的积分是以曲线为界的面积，曲线内随机点与点总数的比值是积分的值。
从概念上讲，有人能解释一个非数学的人，我们如何能够解决PDE/ODE使用蒙特卡罗？
浏览 0提问于2019-04-02得票数 2
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
腾讯会议
实时音视频活动推荐
腾讯云开发者
