設(shè)有寬度為1的正方形，其中有一半徑為0.5的圓內(nèi)接這個正方形，隨機生成一定數(shù)量的點，落在圓內(nèi)的點的數(shù)量和落在圓外點的數(shù)量的比值k即為π/4，我們只需要求出k的值，即可根據(jù)π=k*4即可求得圓周率。
參考資料：http://www.ruanyifeng.com/blog/2015/07/monte-carlo-method.html

代碼如下：

% kkk=[]; % for i=1:10clc dota=1000; x=rand(1,dota); y=rand(1,dota); plot(x,y) dis=[]; for i=1:length(x); k=sqrt((x(i)-0.5)^2+(y(i)-0.5)^2); dis=[dis,k]; end count=0; for i=1:length(dis) if dis(i)<0.5count=count+1;plot(x(i),y(i),'color','red','marker','o','linestyle','none')hold on elseplot(x(i),y(i),'color','blue','marker','o','linestyle','none') end end qw=count/length(x); disp('隨機點個數(shù)：') disp(dota) disp('圓周率:') disp(qw*4) % kkk=[]; % kkk=[kkk,qw*4]; % end % disp(kkk) % jugg=((abs(pi-kkk))./pi)*100; % plot(jugg) % ylabel('相對誤差百分比（%）')

總結(jié)

以上是生活随笔為你收集整理的matlab用蒙特卡洛法计算圆周率的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。