如何計(jì)算π的值

1、蒙特卡羅(Monte Carlo )法

思想：

取一正方形A ，以A 的一個(gè)頂點(diǎn)為圓心，A 的邊長為半徑畫圓，取四分之一圓(正方形內(nèi)的四分之一圓)為扇形B 。已知A 的面積，只要求出B 的面積與A 的面積之比B A

S k S =，就能得出B S ，再由B 的面積為圓面積的四分之一，利用公式2=S R π圓即可求出π的值。因此，我

們的目的就是要找出k 的值。

可以把A 和B 看成是由無限多個(gè)點(diǎn)組成，而B 內(nèi)的所有點(diǎn)都在A 內(nèi)。隨機(jī)產(chǎn)生n 個(gè)點(diǎn)，若落在B 內(nèi)的有m 個(gè)點(diǎn)(假定A 的邊長為1，以扇形圓心為坐標(biāo)系原點(diǎn)。則只要使隨機(jī)產(chǎn)生橫縱坐標(biāo)x 、y 滿足

221x y +≤的點(diǎn)，就是落在B 內(nèi)的點(diǎn))

，則可近似得出k 的值，即m k n =，由此就可以求出π的值。

程序(1)：

i=1;m=0;n=1000;

for i=1:n

a=rand(1,2);

if a(1)^2+a(2)^2<=1

m=m+1;

end

end

p=vpa(4*m/n,30)

總結(jié)

以上是生活随笔為你收集整理的matlab计算π,如何计算π的值(MATLAB)的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

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