蒙特卡罗法

##本文来自网课笔记整理

---

方法简介

无法求得精确解时，进行随机抽样，根据统计试验求近似解

简单例题

圆周率π没有精确值，如何求出其近似值？

分析：

1.π常用于面积求解，通过面积的近似值入手来求其近似解

2.我们取半径为1的圆，作一个边长为4的外切正方形

3.通过随机撒大量点来模拟面积分布

4.统计意义上，圆内点数/总点数=圆面积/正方形面积

5.我们知道，圆面积=π*（r^2）

代码求解

%求解圆周率近似值
clc,clear;
%参数初始化，投放足够多的点（这里为10000）,圆半径为1，圆心坐标为（1,1）
%初始化时圆内有0个点
p=10000;r=1;x0=1;y0=1;n=0;
hold on%保持绘图窗口，多次绘图
for i=1:p%对于要投放的总共p个点
    %rand函数产生（0,1）之间的随机数
    px=rand*2;%随机产生该点的横坐标
    py=rand*2;%随机产生该点的纵坐标
%如果点在园内，颜色设为蓝色并且n加1，否则设为红色
if (px-1)^2+(py-1)^2<1
    plot(px,py,'.b');
    n=n+1;
else
    plot(px,py,'.r');

end
end
axis equal%绘图时横纵坐标的单位长度一致，便于观察圆
s=(n/p)*4;
pi0=s

 运算结果如下

以上为简单示例，不同问题具有不同的概率分布，本题是均匀分布，常见的还有泊松分布，正态分布，指数分布等等

注意：

无法求得精确解时，根据统计试验求近似解

随机性（rand（n）来生成随机数）

统计性（n不能太少）

近似解（如果可以求精确，就不要求近似）