mc1.m

来自「这是一段利用matlab仿真抛骰子问题的蒙特卡洛算法」· M 代码 · 共 36 行

% 蒲丰氏问题   蒙特卡洛仿真
close all;
clear all;
clc;

a = 4 ;         % 参数设置
l = 3 ;
sum = 0 ;
 
N2 = 5000;
    
for j = 1:N2 

epsl1 = rand;   % 产生两个随机变量
epsl2 = rand;

x=a*epsl1;
ct=pi*epsl2;
if (x<l*sin(ct))
    s(j)=1;
else
    s(j)=0;
end 

sum = s(j) + sum ;

end;

result = 2*l/(a*sum/ (N2) )    %  计算所得的 pi
error = result -3.1415926