乱数を用いて円周率算出

ゆとり教育

戦後の詰め込み教育を反省して、ゆとり教育が叫ばれた時期がありました。そのとき、小数演算は難しいので円周と直径の比率は3.14…のところ、3で計算してもよいという通達があったと聞きました。

これでは誤差が大きいのではと気をもんでいましたが、ゆとり教育そのものがとん挫したようです。

乱数を用いて円周率を算出

乱数を用いて円周率を算出する資料は、山のように溢れていますが、仕組みを解説している記事は多くありません。今回、図を描いて丁寧に説明をします。

乱数の種類

乱数には、自然乱数、一様乱数、正規乱数などがありますが、ここではある有限の区間を区切って、その区間内で全ての実数が同じ確率(濃度)で現れるような乱数である一様乱数を採用します。

今回、JavaScriptを使ってプログラミングしますが、一様乱数はMath.random関数が用意されています。

算出の仕組み

一様乱数は0~1(0 <= r < 1)の区間で、同じ確率で発生します。一辺が1の正方形の中に2つの乱数を発生させ、xとyに割り振るとP点が求まります。

そのP点が1/4円の中に入る確率は正方形と1/4円の面積比であるπ/4になります。算出した確率を4倍すれば円周率が求まる仕掛けです。

P点が1/4円の中にあるかどうかは、xとyから半径zを求め、その値が1以下ならば円内にあると判断します。

正方形の面積:1、1/4円の面積:π/4

正しくはx、yの二乗和の平方根の値を求めておく必要があるが、

a、bがともに正で、a<bならば√a<√b
の式が成り立つので平方根を求める必要はありません。

検証

正確な円周率に近づけるにはMAXの値を大きくしますが、この方法で精度の高い円周率を求めることはできません。精度の高い円周率を算出するにはマクロリーン級数による近似値計算などがあります。

サンプルコード

<html>
<script type="text/javascript" charset="Shift_JIS">
// 乱数を用いて円周率を算出する
const MAX = 1000000;
var count = 0;

for(var i=0; i<MAX; ++i){
	var x = Math.random();
	var y = Math.random();
	var z = x * x + y * y;
	if(z < 1) ++count;
}
  /* 円周率を計算・出力 */
  var pi = count / MAX * 4;
  alert("PI="+pi);
</script>
</html>

動作結果

さいごに

乱数を用いて円周率を算出してみました。紀元前2000年ころには古代バビロニアで22/7=3.142857が使われていたといいます。なんぼなんでもπ≒3は暴挙としか思えません。

 

コメントを残す

下記項目は必須ですがメールアドレスは公開されません。名前はニックネームをお使いください。

CAPTCHA


このサイトはスパムを低減するために Akismet を使っています。コメントデータの処理方法の詳細はこちらをご覧ください