【转载】蒙特卡洛方法

蒙特卡罗方法是一种计算方法。原理是通过大量随机样本，去了解一个系统，进而得到所要计算的值。

它非常强大和灵活，又相当简单易懂，很容易实现。对于许多问题来说，它往往是最简单的计算方法，有时甚至是唯一可行的方法。

1 圆周率的计算

 

 

 

staus=10;
for i=1:5 %4次模拟
point=staus.^i; %模拟的随机点数
RandData=rand(2,point); %产生在x~[0,1],y~[0,1]上的随机数根据随机点数，产生随机的(x,y)散点,不明白可以试试 %scatter(RandData(1,:),RandData(2,:))
Below=find(RandData(1,:).^2 + RandData(2,:) .^2 <= 1);%寻找位于圆内部的点
Outcome(i)=length(Below)/length(RandData);%最终结果的表示
end
% BelowData=RandData(:,Below);
hold on
scatter(BelowData(1,:),BelowData(2,:))

 

 

 

 2 定积分的计算

staus=10;
for i=1:4 %4次模拟
point=staus.^i; %模拟的随机点数
RandData=rand(2,point); %根据随机点数，产生随机的(x,y)散点,不明白可以试试 %scatter(RandData(1,:),RandData(2,:))
Below=find(RandData(1,:).^2>RandData(2,:));%寻找位于曲线下的散点
Outcome(i)=length(Below)/length(RandData);%最终结果的表示
end
BelowData=RandData(:,Below);
hold on
scatter(BelowData(1,:),BelowData(2,:))