蒙特卡洛算法的MATLAB实现

一 解决问题

蒙特卡洛算法，即随机取样法，也称为计算机随机模拟方法，源于世界著名赌城——Monte Carlo，算法由此得名。

它是基于概率论的一种计算方法：即通过随机数或者伪随机数来计算比较难以解决的问题，因为设计概率问题，那么结果就不是一个准确的值，而是一个近似值，并且所使用的样本（统计点数）越多，所得到的结果越准确。例如：不规则图形的面积问题等。

二 基本原理

我们以面积计算的问题讲解蒙特卡洛算法的实现步骤如下：

1.画出图像

求要计算图像的基本形状。此处以y1=3x；y2=8-x为例，画出图像如下：

2.确定边界

通常为矩形，要求矩形必须完整包括所求图像。

如上图，因为图像交点为（2,6），因此，确定边界如红框所示：（0,0），（8,0），（8,6），（0,6）

3.随机统计

在（矩形）边界范围内随机产生点，并统计落在所求图像中的点。所用到的函数为unifrnd函数。本例子以10^7（7次方在保证运行速度的情况下，基本可以满足准确度ÿ