真性乱数による円周率実験ツール

正方形にランダムな点を散らし、四分円の内側に入る点の割合から円周率を推定します。

モンテカルロ円周率実験
エントロピー源とサンプル数を選び、実験の推定値を円周率と比較します。

RANDOM.ORG の大気雑音から取得した 128 ビットを実験のシードに使います。

円周率実験について

円周率は、円周の長さと直径の一定の比です。小数は無限に続き循環しませんが、簡単な確率実験でその値を推定できます。単位正方形の中に半径1の四分円があると考えてください。正方形の面積は1、四分円の面積は円周率を4で割った値です。点が正方形に一様に分布すれば、四分円に入る割合は円周率を4で割った値に近づきます。観測した割合を4倍すると、モンテカルロ法による円周率の推定値になります。 本ツールはゼロから1までの座標の組を生成し、各組について x の二乗と y の二乗の和が1以下かを判定します。この不等式で四分円の内側または境界上の点を特定できます。ランダムサンプリングには自然なばらつきがあるため、小規模な実験ではおなじみの 3.14159 から上下に大きくずれる場合があります。点を増やすと通常は精度が改善します。典型的な誤差はサンプル数の平方根の逆数に比例して減るため、小数点以下をもう1桁正確にするには、およそ100倍の観測が必要になることがあります。 乱数源は2種類です。真性乱数モードは実行ボタンを押すと RANDOM.ORG から大気雑音のエントロピーを 128 ビット直接取得し、高品質な疑似乱数列のシードとして全座標を生成します。大気中の無線雑音はアルゴリズムの数列ではなく、物理現象です。ネットワークサービス、利用枠、応答データが利用できない場合は、代替を使ったことを明示し、ブラウザーの暗号学的に安全な乱数生成器で実験を完了します。ローカルモードはこの CSPRNG を直接使い、ネットワークリクエストを一切行いません。 この実験は学習用であり、円周率を最も効率よく求める方法ではありません。数学的アルゴリズムなら、確定的な保証のもとではるかに速く何兆桁も計算できます。モンテカルロ法の価値は、幾何、面積、確率、シミュレーション、大数の法則を目に見える1つの計算で結び付けることです。同じ原理は、不規則な領域の面積、多次元積分、金融リスク、粒子の挙動など、厳密な式が難しい問題にも応用できます。同じサンプル数で繰り返してばらつきを観察し、その後で点を増やして推定が安定する様子を確認してください。

円周率実験の例

実験典型的な結果わかること
1,000 個のランダムな点約 3.0 から 3.3サンプリングのばらつきが見える短時間の実演です。
10,000 個のランダムな点多くの場合 3.14 付近速さと実用的な推定精度のバランスが取れています。
100,000 個のランダムな点通常は 3.1416 により近い値観測数を増やすと、一般にランダム誤差が減ります。

円周率実験の実行方法

  1. 乱数源に真性乱数の大気雑音またはブラウザーのローカル CSPRNG を選びます。
  2. 100 から 1,000,000 までの整数で点の数を入力します。
  3. 「円周率実験を実行」を選択します。大気エントロピーの取得時は少しお待ちください。
  4. 円周率の推定値と内側の点数を、期待される幾何学的な比率と比較します。

円周率実験のよくある質問

この実験で円周率を推定できるのはなぜですか?

単位四分円の面積は円周率を4で割った値で、それを囲む単位正方形の面積は1です。一様なランダム点が四分円に入る確率も円周率を4で割った値になります。

毎回結果が変わるのはなぜですか?

実行のたびに新しいランダムサンプルを使うため、観測割合は自然に変動します。大きなサンプルほど変動は小さくなる傾向がありますが、有限回の実行でより近くなる保証はありません。

大気雑音による乱数とは何ですか?

RANDOM.ORG が予測できない大気中の無線雑音を測定し、乱数に変換します。本ツールはその値をシード付き生成器のエントロピーとして使い、シミュレーションに必要な多数の座標を生成します。

RANDOM.ORG に接続できない場合はどうなりますか?

ブラウザーの暗号学的に安全な乱数生成器へ切り替え、正直に警告を表示します。実験はローカルで継続し、生成源のバッジには実際に使用したものが表示されます。

点はいくつ使えばよいですか?

1万点なら授業で素早く実演でき、10万点なら通常はより安定します。モンテカルロ誤差の減少は遅いため、非常に高い精度にはこの対話型ツールの上限をはるかに超える点が必要です。