0.50
1
18

ビュフォンの針についての詳細

ビュフォンの針

長さ l の針を、間隔 d で平行線が引かれた床に落とします。針が間隔より長くなければ、線と交わる確率は 2l/(πd) です。PictorX が投下を代わりに行い、1本ずつ描き、数え方を逆にして同じ実験から π を推定します。 再生を押すと針がランダムに落ちます。交差から π を推定します。針の長さが間隔の半分なら、π はおよそ投下数 ÷ 交差数です。床の上をドラッグすると針を1本手で置けます。その針は数えられません。

使い方

  1. このページを開いて再生を押します。針が線の上にランダムに落ちます。
  2. π が落ち着いていく様子を見ます。針が間隔の半分なら、π は投下数 ÷ 交差数です。
  3. サイン図を開く、針を1本手で置く、または1,000本を一度に落としてから、PNG を保存します。

2つの数が交差を決める

1回の投下は、針の中心から最も近い線までの距離 x と、針と線のなす角 θ で決まります。x は 0 から d/2 まで、θ は 0 から π までです。針には頭がないので、半回転ですべての向きをカバーできます。x ≤ (l/2) sin θ のとき、針は線と交わります。

短い針の公式

l ≤ d のとき、交差の確率は P = 2l/(πd) です。組 (θ, x) は幅 π、高さ d/2 の長方形を埋めます。交差にあたるのはサイン曲線の下の面積で、その面積は l です。2つの面積を割ると、分母に π が残ります。l = d/2 では確率が 1/π になるので、π はおよそ投下数 ÷ 交差数です。

π を求める遅い方法

式を変形すると π ≈ 2l·N / (d·N_hit) になります。これはモンテカルロ推定で、しかも効率の悪いものです。数千回投げれば、たいてい最初の数桁が見えてきます。次の1桁にはずっと多くの回数が必要です。表示の横の ± は比の標準誤差1つ分です。床の上をドラッグすると針を1本置いて、その x と θ を確認できます。その針は描かれますが、数には含まれません。

よくある質問

ビュフォンの針とは?

間隔 d の平行線の上に落とす長さ l の針です。短い針(l ≤ d)なら、線と交わる確率は 2l/(πd) です。

なぜ答えに π が入るのですか?

針の向きは 0 から π までの角度です。起こりうるすべての落ち方はその幅の長方形を埋め、交差はサイン曲線の下の面積になります。これらの面積を割ると、分母に π が残ります。

この実験はどうやって π を推定するのですか?

投下数 N と交差数 N_hit を数えると、π ≈ 2l·N / (d·N_hit) です。針が間隔の半分なら、これは投下数 ÷ 交差数になります。

手で置いた針はなぜ数えないのですか?

推定にはランダムな位置とランダムな角度が必要です。ドラッグして置いた針は x と θ が見えるように描かれますが、N にも N_hit にも加えられません。

試してみますか?

ブラウザで PictorX を開きます。インストール不要、始めるのにクレジットカードも不要です。

エディターへ