小编Soo*_*kie的帖子

初学者Python蒙特卡洛模拟

我是Python的初学者,正在通过我们的讲师进行练习。我正在为这个问题而苦苦挣扎。

在Python编辑器中,编写一个Monte Carlo模拟以估计数字的值。具体来说,请按照以下步骤操作:A.产生两个数组,一个称为x,一个称为y,每个数组包含100个元素,它们在-1和1之间随机且均匀地分布于实数中。情节。相应地标记轴。C.写下一个数学表达式,该数学表达式定义哪些(x,y)对数据点位于以图的(0,0)为原点,半径为1的圆中。D.使用布尔蒙版来标识圆内的点,并以不同的颜色和标记大小在已在B中绘制的数据点之上对其进行过度绘制。

这就是我目前所拥有的。

import numpy as np
import math
import matplotlib.pyplot as plt
np.random.seed(12345)
x = np.random.uniform(-1,1,100) 
y = np.random.uniform(-1,1,100) 
plt.plot(x,y) //this works


for i in x:
    newarray = (1>math.sqrt(y[i]*y[i] + x[i]*x[i]))
plt.plot(newarray)
Run Code Online (Sandbox Code Playgroud)

有什么建议么?

python montecarlo

3
推荐指数
1
解决办法
246
查看次数

简单的Matlab优化

鉴于此问题,我被要求优化代码。我似乎找不到找到使它更快运行的方法。

nmax= 2000;
mmax= 2000;
    for m=1:mmax
        for n=1:nmax
            A(n,m)= n+m;
        end 
    end
Run Code Online (Sandbox Code Playgroud)

有任何想法吗?

matlab vectorization

-2
推荐指数
1
解决办法
48
查看次数

标签 统计

matlab ×1

montecarlo ×1

python ×1

vectorization ×1