数学建模--蒙特卡洛模型的Python实现

目录

1.算法思想简介

2.算法应用1:问题一阐述

3.算法应用1:问题一解决

4.算法应用2:问题二阐述

5.算法应用2:问题二解决


1.算法思想简介

复制代码
#蒙特卡洛算法思想
"""
蒙特卡洛方法的理论其实很类似于概率论中一个比较重要的定理---大数定律。
其基本原理简单描述是先大量模拟,然后计算一个事件发生的次数,再通过这个发生次数除以总模拟次数,得到想要的结果。
这个时候我们在样本足够大的情况下我们认为这个事件发生的概率等于实验所得到的结果
"""

2.算法应用1:问题一阐述

复制代码
"""
Question1:
计算圆周率pi(Π)值
实验原理:
在正方形内部有一个相切的圆,圆面积/正方形面积之比是(PixRxR)/(2Rx2R)= Pi/4。
在这个正方形内随机产生n个点,假设点落在圆内的概率为P,那么P=圆面积/正方形面积,P= Pi/4。
如何计算点落在圆内的概率P?
可以计算点与中心点的距离,判断是否落在圆的内部,若这些点均匀分布,用M表示落到圆内投点数,N表示总的投点数,则圆周率Pi=4P=4xM/N。
"""

3.算法应用1:问题一解决

python 复制代码
#第一题求解代码
#我们假设r=1,a=0,b=0
import numpy as np
def Qustion1(n):
    r=1.0
    a=0
    b=0
    x_min=a-r
    x_max=a+r
    y_min=b-r
    y_max=b+r
    m=0#m是计算落在圆内的点的个数
    for i in range(n+1):
        x=np.random.uniform(x_min,x_max)
        y=np.random.uniform(y_min,y_max)
        if x*x+y*y<=1:
            m=m+1
    print("第一题的答案是:",4*(m/float(n)))
ans=Qustion1(1000000)
#cout:第一问题的答案是: 3.141468

4.算法应用2:问题二阐述

python 复制代码
""""
Question2:
计算函数定积分值:
实验原理:
若要求函数f(x)从a到b的定积分,我们可以用一个比较容易算得面积的矩型包围在函数的积分区间上(假设其面积为Area),
定积分值其实就是求曲线下方的面积。随机地向这个矩形框里面投点,统计落在函数f(x)下方的点数量占所有点数量的比例为P,
那么就可以据此估算出函数f(x)从a到b的定积分为Area*P。
此处我们将a和b设为0和1,函数f(x)=x^2。
"""

5.算法应用2:问题二解决

python 复制代码
#第二题求解代码
#此处我们将a和b设为0和1,函数f(x)=x^2。
def Question2(n):
    x_min=0
    y_min=0
    x_max=1
    y_max=1
    m=0
    for i in range(n+1):
        x=np.random.uniform(x_min,x_max)
        y=np.random.uniform(y_min,y_max)
        if y-x*x<0:
            m=m+1
    print("第二题的答案是:",m/float(n))
ans2=Question2(1000000)
#第二题的答案是: 0.334016
相关推荐
一个天蝎座 白勺 程序猿43 分钟前
Python爬虫(8)Python数据存储实战:JSON文件读写与复杂结构化数据处理指南
爬虫·python·json
q_q王1 小时前
dify对接飞书云文档,并且将图片传入飞书文档
python·大模型·飞书·dify·智能体·图片展示
noravinsc1 小时前
django filter 排除字段
后端·python·django
zandy10112 小时前
嵌入式BI开发指南:如何通过衡石API将分析能力集成到业务系统?
开发语言·python·嵌入式
曲幽2 小时前
零基础快速搭建AI绘画网站!用Gradio玩转Stable Diffusion
python·ai作画·stable diffusion·gradio·diffusers·webui
2401_890665862 小时前
免费送源码:Java+ssm+HTML 三分糖——甜品店网站设计与实现 计算机毕业设计原创定制
java·python·微信小程序·html·php·课程设计·android-studio
noravinsc2 小时前
django filter 日期大于当前日期的
python·django
悲喜自渡7212 小时前
pytorch & python常用指令
人工智能·pytorch·python
Star abuse3 小时前
Python爬虫课程实验指导书
开发语言·爬虫·python
秋名RG3 小时前
简单了解Java的I/O流机制与文件读写操作
java·开发语言·python