用蒙特卡罗法计算圆周率（附matlab程序）(详细)

用蒙特卡罗法计算圆周率（附matlab程序）(详细)

蒙特卡罗法是一种以概率统计理论为指导的一类非常重要的数值计算方法。是指使用随机数（或更常见的伪随机数）来解决很多计算问题的方法。

现在正方形里有一个与之相切的圆形，建立直角坐标系如图所示：

在正方形内随机投点，设投了n次点，若投的点足够得多，多到投的点都布满了整个正方形（点与点不可重叠），在圆内的点数与整个正方形点数之比就等于圆形与正方形的面积之比，随后即可得到圆形的面积，利用圆的面积公式，可以得到pi的值。

实现的matlab算法：

%蒙特卡洛法求圆周率

%任意定一个正方形的的面积，建直角坐标系，其原点定在正方形的中点（0，0）
a_length=2;%定正方形的边长为2
n=10^7;%进行随机过程的次数
r=a_length/2;%与正方形相切圆的半径
s_tangle=a_length^2;%正方形面积

%随机打点求圆形面积
j=0;
for i=1:1:n
x=rand(1,1)a_length-a_length/2;%产生[a,b]的随机数：x=rand(1,1)(b-a)+a
y=rand(1,1)*a_length-a_length/2;
if x^2+y2<=r^2
j=j+1;
end
end

%正方形与圆形的面积之比k
k=j/n;
s_circle=s_tangle*k;
pi_estimate=s_circle/r/r

---

结果：
pi_estimate =

3.1426

若有不足之处还恳请各位网友批评指正！