weixin_52014585的博客

原创 【MATLAB】蒙特卡洛估计pi值

x =-1+ 2*rand(1e4,1);y =-1+ 2*rand(1e4,1);%搜寻落到半径为1的圆内的个数z = x((x.^2 + y.^2 )<=1);%其实是先判断，后索引符合条件的x值pitest = 4*length(z)/length(x); %原理: 所有随机数落到半径为1的圆内的概率等err = abs(pitest - pi)/pi;fprintf( '2.梦塔卡罗模拟π的值为%3.4f,相对误差为%3.4f\n\n ' ,pitest,err)figur.