0
点赞
收藏
分享

微信扫一扫

MATLAB 蒙特卡罗法求π

舍予兄 2022-03-30 阅读 82

n = [1:10000];
PI=[];
for i = n
    m=0;
    for j = 1:i
	x = rand;
	y = rand;
    
    if (x^2 + y^2) < 1   
		m = m + 1;
    end
    
    end
    p=4 * m/i;
    PI=[PI,p];
end
plot(PI)
fprintf('圆周率值是:%d',4 * m/n);  

思路:

1 蒙特卡罗法:使用类似抛针法,统计投入圆内的数量

2 建立一个抛掷次数的list,储存次数(1~10000),用于用图像表明随着实验次数的增加,π值趋近于真实值

3 嵌套循环:对于实验次数n,我们建立一个子循环,通过if语句判断每次产生的实验结果,最后在母循环里进行统计

举报

相关推荐

0 条评论