数模方法论-蒙特卡洛法

一、基本概念

蒙特卡洛法是一种基于随机抽样的数值计算方法,主要用于估计复杂系统的数值解。其基本原理是通过生成大量随机样本,来模拟系统的行为或求解特定的数学问题,比如积分、概率和优化等。

在应用上,蒙特卡洛法可以用于金融风险评估、物理模拟、运筹学等领域。例如,在金融中,它可以模拟不同市场条件下的资产价格变化,从而帮助投资者评估风险和回报。通过随机生成样本并进行统计分析,最终得出问题的近似解,从而有效处理不确定性。

假设我们想估计一个单位圆的面积。可以通过蒙特卡洛法来实现:

  1. 在一个边长为2的正方形内随机生成大量点(例如,10000个)。
  2. 计算这些点中落在单位圆内的点的比例。
  3. 单位圆的面积为π,正方形的面积为4,因此圆的面积可以通过比例来估算:面积 ≈ (落在圆内的点数 / 总点数) × 4。

通过增加点的数量,可以逐步提高估算的准确性。

二、实例求解

例题一

Matlab求解

Matlab 复制代码
clc, clear
x=unifrnd(0,12,[1,10000000]);
y=unifrnd(0,9,[1,10000000]);
pinshu=sum(y<x.^2 & x<=3)+sum(y<12-x & x>=3);
area_appr=12*9*pinshu/10^7

Python求解

python 复制代码
import numpy as np

# 生成随机数
x = np.random.uniform(0, 12, 10000000)
y = np.random.uniform(0, 9, 10000000)

# 计算 pinshu
pinshu = np.sum((y < x**2) & (x <= 3)) + np.sum((y < 12 - x) & (x >= 3))

# 计算近似面积
area_appr = 12 * 9 * pinshu / 10**7

print(area_appr)

例题二

Matlab求解

Matlab 复制代码
rand('state',sum(clock));  %初始化随机数发生器
p0=0;
tic    %计时开始
for i=1:10^6
   x=randi([0,99],1,5); %产生一行五列的区间[0,99]上的随机整数
   [f,g]=mengte(x);
   if all(g<=0)
       if p0<f
           x0=x; p0=f; %记录下当前较好的解
       end
   end
end
x0,p0
toc    %计时结束

Python求解

python 复制代码
import numpy as np
import time

# 初始化随机数发生器
np.random.seed(int(time.time()))

p0 = 0
start_time = time.time()  # 计时开始

def mengte(x):
    # 这里你需要根据实际的 mengte 函数实现来定义它
    # 作为示例,这里假设返回随机的 f 和 g 值
    f = np.random.random() * 100  # 示例 f
    g = np.random.random(5) - 0.5  # 示例 g,可能有负值
    return f, g

for i in range(10**6):
    x = np.random.randint(0, 100, 5)  # 产生一行五列的区间[0,99]上的随机整数
    f, g = mengte(x)
    if np.all(g <= 0):
        if p0 < f:
            x0 = x.copy()  # 记录下当前较好的解
            p0 = f

print("最佳解:", x0)
print("最佳值:", p0)
print("耗时:", time.time() - start_time)  # 计时结束
相关推荐
搞不懂语言的程序员4 分钟前
装饰器模式详解
开发语言·python·装饰器模式
mosquito_lover19 分钟前
矿山边坡监测预警系统设计
人工智能·python·深度学习·神经网络·视觉检测
王禄DUT9 分钟前
化学方程式配平 第33次CCF-CSP计算机软件能力认证
开发语言·c++·算法
船长@Quant9 分钟前
PyTorch量化进阶教程:第二章 Transformer 理论详解
pytorch·python·深度学习·transformer·量化交易·sklearn·ta-lib
wuqingshun31415911 分钟前
蓝桥杯 XYZ
数据结构·c++·算法·职场和发展·蓝桥杯
DreamByte24 分钟前
C++菜鸟教程 - 从入门到精通 第五节
开发语言·c++·算法
南玖yy31 分钟前
数据结构C语言练习(两个队列实现栈)
c语言·数据结构·算法
攻城狮7号41 分钟前
【第一节】Python爬虫基础-HTTP基本原理
爬虫·python·python爬虫
IT乐手41 分钟前
adb logcat 写文件乱码的解决方案
android·python
Python测试之道1 小时前
Deepseek API+Python 测试用例一键生成与导出 V1.0.6(加入分块策略,返回更完整可靠)
开发语言·python·测试用例