蒙特卡罗解题归结为三个主要步骤:
构造或描述概率过程;
实现从已知概率分布抽样;
建立各种估计量。
接下来我们介绍3个简单的案例,看一下如何在实际问题中应用这3个步骤进行求解。...正方形内部有一个相切的圆,假设圆形的半径为 ,那么正方形的边长则为 ,因此,圆形和正方形的面积之比为:
image.png
现在,在正方形内部,随机产生一系列的点,计算他们与中心点的距离,从而判断是否落在圆的内部...把所有点画出来:
# 绘制正方形
fig = plt.figure(figsize=(10, 10))
axes = fig.add_subplot(1, 1, 1)
plt.plot(x, y, 'ro...', markersize=1)
plt.axis('equal')
# 绘制圆形
circle = Circle(xy=(a, b), radius=r, alpha=0.5, color='gray...把所有点画出来:
# 绘制正方形
fig = plt.figure(figsize=(10, 10))
axes = fig.add_subplot(1, 1, 1)
plt.plot(x, y, 'ro