python 编程练习——使用蒙特卡洛方法计算阴影部分面积

Python YX ⋅ 于 2020-05-10 12:02:44 ⋅ 最后回复由 姜姜 2020-05-13 17:22:00 ⋅ 335 阅读
使用python代码,完成如图习题

file

file

参考文献

中学数学建模与赛题集锦
file

成为第一个点赞的人吧 :bowtie:
回复数量: 2
  • 热心市民小杨
    2020-05-13 17:00:25
    import numpy as np
    import matplotlib.pyplot as plt 
    %matplotlib inline
    n=2000
    c=0
    x_a=[]
    y_a=[]
    x=np.random.random(2000)
    y=np.random.random(2000)
    for i in range(2000):
        if x[i]**2 > y[i]:
            x_a.append(x[i])
            y_a.append(y[i])
            c+=1
    rate=c/n
    print(rate)# 点落在D区域的概率
    x1=np.arange(0,1,0.01)
    y1=x1**2
    plt.plot(x1,y1,'r')
    plt.scatter(x_a,y_a)
  • 姜姜
    2020-05-13 17:22:00
    import random
    
    def integral():
        n = 1000000
        count = 0
    
        x_min, x_max = 0.0, 1.0
        y_min, y_max = 0.0, 1.0
    
        for i in range(0, n):
            x = random.uniform(x_min, x_max)
            y = random.uniform(y_min, y_max)
    
            if x*x > y:
                count += 1
    
        integral_value = count / float(n)
        print(integral_value)
    
    integral()
暂无评论~~
  • 请注意单词拼写,以及中英文排版,参考此页
  • 支持 Markdown 格式, **粗体**、~~删除线~~、`单行代码`, 更多语法请见这里 Markdown 语法
  • 支持表情,使用方法请见 Emoji 自动补全来咯,可用的 Emoji 请见 :metal: :point_right: Emoji 列表 :star: :sparkles:
  • 上传图片, 支持拖拽和剪切板黏贴上传, 格式限制 - jpg, png, gif
  • 发布框支持本地存储功能,会在内容变更时保存,「提交」按钮点击时清空
Ctrl+Enter