■ドローネー集合(その89)

 商がaになる確率は

  log2(1+1/a)−log2(1+1/(a+1))

=log2((a+1)^2/((a+1)^2−1))

 商が1になる確率はlog2(4/3)=41.504%

 商が2になる確率はlog2(9/8)=16.993%

 商が3になる確率はlog2(16/15)=9.311%

 商が4になる確率はlog2(25/24)=5.889%

===================================

 ところで,n→∞のとき,

  Σlog2(n^2/(n^2−1))→1

は保証されているのだろうか?

  Σlog2(n^2/(n^2−1))

=Σ{log2(1+1/k)−log2(1+1/(k+1))}

=log22−log2(1+1/2)+・・・−log2(1+1/n)

→log22=1

===================================

【補】ベンフォードの法則

 1938年,GEの物理学者ベンフォードは対数表の対数表の最初が残りの部分よりもひどく汚れていることに気づき,「1ではじまる数が多いのはなぜか」という問題に説明を与えました.

 先頭の数字がどのような確率で出現するかを考えましょう.単純に各数字(0〜9)の出現確率が同じと考えれば,同じ確率1/9で現れるはずですが,実際には1から始まる数値が圧倒的に多く30%くらいもあります.

 たとえば,簡単な例として,2のベキ乗2^nを順に並べてそれぞれの最大桁の数を取り出すと

  2,4,8,16,32,64,128,256,512,1024,2048,・・・

  →2,4,8,1,3,6,1,2,5,1,2,・・・

となっているのですが,倍にした数が9で始まるためには,その前の数字が45−49で始まっていなければなりません.それに対して,5−9で始まる数はどれも倍にすると1で始まる数になります.そして,最大桁がk(1≦k≦9)である確率はn→∞のとき,

  log10((k+1)/k)

に収束することが知られています.

 したがって,最大桁の頻度は1が一番高く

  1→log102=0.3010,

以下,

  2→log103/2=0.1761,

  3→log104/3,

  ・・・・・・・・・,

  9→log1010/9=.0458

の順になるというわけです.

===================================