円周率

出典: フリー百科事典『ウィキペディア(Wikipedia)』
移動先: 案内検索

円周率(えんしゅうりつ)は、周長の直径に対する比率として定義される数学定数である。通常、ギリシア文字 π(パイ、ピー、ラテン文字表記: pi 英語発音: [pai])で表される。数学をはじめ、物理学工学といった様々な科学分野に出現し、最も重要な数学定数とも言われる。

円周率は無理数であり、その小数展開は循環しない。小数点以下35桁までの値は次の通りである。

π = 3.14159 26535 89793 23846 26433 83279 50288 …

円周率は、無理数であるのみならず、超越数でもある。

基礎[編集]

表記と呼び方[編集]

円周率を表すギリシア文字 π は、ギリシア語 περίμετρος[1][2][3](ペリメトロス)あるいは περιφέρεια[4](ペリペレイア)の頭文字から取られた[注 1]。いずれも周辺・円周・周などを意味する。文字 πウィリアム・オートレッド1631年に著した著書において半円の円弧部分の長さを表す文字として用い、アイザック・バローは論文において半径 R の円周の長さとして用いた[5]ウィリアム・ジョーンズ英語版 (1706) やレオンハルト・オイラーらにより(現代と同じく)円周の直径に対する比率を表す記号として用いられ、それが広まった[1][2][5]日本では「パイ」と発音する。

π を指す言葉には、日本・中国・韓国における「円周率(圓周率)」、ドイツの「Kreiszahl」(Kreis は円(周)、Zahl は数の意)の他、それを計算した人物の名前を取った「アルキメデス数」(: Archimedes' constant)、「ルドルフ数」(: Ludolph's constant: Ludolphsche Zahl)などがある。一般にドイツ語を除いたヨーロッパの諸言語には「円周率」に対応する単語はない[3][6]

なお、「π」の字体は、表示環境によってはキリル文字п に近い π などと表示されることがある。

また、文字「π」は、数学では他に素数計数関数基本群ホモトピー群にも用いられる。またある種の写像を表すときにも慣習的に用いられる。

定義[編集]

直径 1 の円の周長は π

平面幾何学において、円周率 π は、周長の直径に対する比率として定義される。すなわち、円の周長を C, 直径を d としたとき、

である。全ての円は互いに相似なので、この比率は円の大きさに依らず一定である。

ところが、この定義は円の周長を用いているため、曲線の長さを最初に定義していない解析学などの分野では、π が現れる際に問題となることがある。この場合、円の周長に言及せず、解析学などにおける性質の一つを π の定義とすることが多い[7]。この際の π の定義の一般なものとして、三角関数 cos x0 を取るような x > 0 の最小値の2倍とするもの、級数による定義、定積分による定義などがある。後述の#円周率に関する式も参照。

歴史[編集]

円に内接する正多角形による π の近似
円に内接・外接する正多角形による π の近似。アルキメデスによる計算。

古代[編集]

円周の直径に対する比率が円の大きさに依らず一定であり、それが 3 より少し大きい程度だということは古代エジプトバビロニアインドギリシアの幾何学者たちにはすでに知られていた。また、古代インドやギリシアの数学者たちの間では半径 r の円板の面積が πr2 であることも知られていた。さらに、アルキメデスは半径 rの体積が (4/3)πr3 であることや、この球の表面積4πr2(その球の大円による切り口の面積の4倍)であることを示した。

2千年紀[編集]

14世紀インド数学者天文学者であるサンガマグラマのマーダヴァは次のような π級数表示を見いだしている(ライプニッツの公式):

これは逆正接関数 Arctan xテイラー展開x = 1 での実現になっている。マーダヴァはまた、

を用いて π の値を小数点以下11桁まで求めている。

18世紀フランスの数学者アブラーム・ド・モアブルは、コインを 2n 回投げたときに表が x 回出る確率は、n が十分大きいとき、ある定数 C を取ると、

近似できることを、n = 900 における数値計算により見いだした。この正規分布の概念は1738年に出版されたド・モアブルの『巡り合わせの理論』に現れている。ド・モアブルの友人のジェイムズ・スターリングは後に、C = 1/2π であることを示した。

1751年ヨハン・ハインリッヒ・ランベルトは、x0 でない有理数ならば正接関数 tan x の値は無理数であることを示し、その系として π は無理数であることを導いた。さらに1882年フェルディナント・フォン・リンデマンπ が超越数であることを示し、円積問題(与えられた長さを半径とする円と等積の正方形作図する問題)は解くことができないことを導いた。

コンピュータによる計算の時代[編集]

20世紀以降、コンピュータの発達により、計算された円周率の桁数は飛躍的に増大した。1949年に、ジョン・フォン・ノイマンはコンピュータ ENIAC を使い70時間かけて、円周率を2037桁まで計算した[8]。その後の数十年間、さまざまな計算機科学者によって計算は進められ、1973年には100万桁を超えた。この進歩は高速なハードウェアの開発だけによるものではなく、効率のよいアルゴリズムが考案されたためである。そのうちの最も重要な発見の一つとして、1960年代高速フーリエ変換がある。これにより、多倍精度の演算が高速に実行できるようになった。

2014年の時点では、円周率は小数点以下13.3桁まで計算されている[9]

性質[編集]

無理性[編集]

π無理数である。つまり、2つの整数の商で表すことはできず、小数展開は循環しない。このことは1761年ヨハン・ハインリッヒ・ランベルトが証明したが、厳密性に欠けた部分があった。その部分は1806年ルジャンドルによって補われた。

したがって、円周率のコンピュータによる計算や暗唱10進法における各数字 (0, 1, ..., 9) の出現頻度は、興味の対象となる。

超越性[編集]

さらに、π超越数である。つまり、有理数係数の代数方程式の根とはならない。これは1882年フェルディナント・フォン・リンデマンによって証明された(リンデマンの定理)。特に、整数から四則演算冪根をとる操作だけを有限回組み合わせて π の正確な値を求めることはできないことが分かる。

π が超越数であることより、古代ギリシアの三大作図問題の内の一つ「円積問題」(与えられた長さを半径とする円と等積の正方形を作図すること)が不可能であることが従う。

ランダム性[編集]

π は現在小数点以下10兆桁を超える桁まで計算されている。0 から 9 までの数字がランダムに現れているようには見えるが、実際は、π正規数であるかどうかは分かっていない。例えば π10進表示において、各桁を順に取り出して得られる数列

3, 1, 4, 1, 5, 9, 2, 6, 5, 3, 5, …(オンライン整数列大辞典の数列 A796

には、0 から 9 が均等に現れるのかどうか、すなわち、この数列が乱数列になっているかどうかは分かっていない。それどころか、0 から 9 がそれぞれ無数に現れるのかどうかすら分かっていない。

したがって、10兆桁以降の桁についてもランダムであるかどうかは、現在分かっていないのである。

ベイリーとクランドールの2000年の発表によると、ベイリー=ボールウィン=プラウフの式を用いて2進表示で様々な桁の計算をした結果では、各数字の出現率はカオス理論に基づいていると推測できるようである。

5兆桁までの数字の出現回数は以下の通りである。全てほぼ等しく(約0.0005%の違いに収まる)、最も多いのは 8 で、最も少ないのは 6 である。

0:4999億9897万6328回
1:4999億9996万6055回
2:5000億0070万5108回
3:5000億0015万1332回
4:5000億0026万8680回
5:4999億9949万4448回
6:4999億9893万6471回
7:5000億0000万4756回
8:5000億0121万8003回
9:5000億0027万8819回


未解決問題[編集]

円周率に関する式[編集]

π についての式は非常に多い。ここではその一部を紹介する。数式によってはそれ自体が π の定義になり得るし、π近似値の計算などにも使われてきた。

幾何[編集]

解析[編集]

  • ライプニッツの公式#2千年紀も参照)
  • #2千年紀も参照)
  • ウォリス
  • 1735年オイラーバーゼル問題ゼータ関数
  • オイラー
  • Bnベルヌーイ数
  • ガウス積分
  • ガンマ関数
  • [10]
  • 逆三角関数主値を取るものとすると
  • 逆三角関数(逆正弦関数)の公式より
  • 逆三角関数(逆正接関数)の公式より
初期値の設定:
反復式:an, bn が希望する精度(桁数)になるまで以下の計算を繰り返す。小数第 n 位まで求めるとき log2 n 回程度の反復でよい。
π の算出:円周率 π は、an, bn, tn を用いて以下のように近似される。
非常に収束が早く[注 2]金田康正が1995年に42億桁、2002年に1.24兆桁を計算したスーパーπに使われていた。
  • スターリングの近似f(n) ∼ g(n) を表す)
  • φ(k)オイラーのφ関数
  • eπi + 1 = 0オイラーの等式
(オイラーの等式の一般式)
後者は 2 以上の任意の整数 n に対して成り立ち、1n 乗根全ての和は 0 であることを意味している。n = 2 とするとオイラーの等式を得る。
  • ラマヌジャン
  • (ラマヌジャン)
  • ビエト
  • マチン、1709年)
ただし arctan x主値
を取るものとする。
と書かれることもある。収束が早いため、コンピュータ初期の高精度整数演算の練習問題として、よく使われた。
  • 連分数表示:
  • (チュドノフスキー)
(各項の素因数分解:
13591409 = 13 × 1045493,
545140134 = 2 × 32 × 7 × 11 × 19 × 127 × 163,
640320 = 26 × 3 × 5 × 23 × 29.
  • (チュドノフスキー)
  • (David Bailey, Peter Borwein and Simon Plouffe、俗称“BBP”)
  • (Adamchik and Wagon)
  • (Fabrice Bellard)[11][12]

数論[編集]

力学系・エルゴード理論[編集]

ロジスティック写像 xi+1 = 4xi(1 − xi) により帰納的に定まる数列 {xi} を考える。初期値 x00 以上 1 以下に取るとき、そのほとんど全てで、次が成り立つ。

統計[編集]

  • 正規分布確率密度関数
  • 1 の無数の平行線の上から長さ 1/2 の針を落とすとき、その針が直線と共有点を持つ確率は 1/π である(ビュフォンの針)。

その他[編集]

  • 河川の長さの水源河口間の直線距離に対する比率は、平均すると円周率に近い[13]

暗唱[編集]

語呂合わせ[編集]

π の桁を記憶術に頼らずに暗記する方法が各種存在している。

日本語では、語呂合わせにより、長い桁を暗記するのも比較的簡単である。有名なものとして、以下がある。

産医師異国ニ向コー、産後厄無ク産婦御社ニ虫サンザン闇ニ鳴ク

マーティン・ガードナー著、金沢養訳、『現代の娯楽数学 新しいパズル・マジック・ゲーム』(白揚社、1960年)144頁
かう さん ざん
3. 1 4 1 59 2 6 5 3 5 89 7 9 3 2 3 8 4 6 2 64 3 3 83 2 7 9 (30桁))

英語圏では語呂合わせがうまくいかないため、単語の文字数で覚える方法がある。

Yes, I have a number.
3. 1 4 1 6 (小数点以下4桁までで四捨五入)
Can I find a trick recalling Pi easily?
3. 1 4 1 5 9 2 6 (7桁、また「π を簡単に思い出せるトリックってある?」という文章自体がその質問の答えにもなっている)
How I want a drink, alcoholic of course, after the heavy lectures involving quantum mechanics!
3. 1 4 1 5 9 2 6 5 3 5 8 9 7 9 (14桁)
And if the lectures were boring or tiring, then any odd thinking was on quartic equations again.
3 2 3 8 4 6 2 6 4 3 3 8 3 2 7 9 5 (上に続けて、31桁)S. ボトムリー

これらのような覚え方は多くあり、日本語では上記のものの改編で90桁までのものや、歌に合わせたもの、数値を文字に置き換えて1,000桁近く覚える方法などがある。

暗唱記録[編集]

ギネス世界記録』によれば、円周率暗唱の世界記録は2005年11月20日に6万7890桁を暗唱した中国人、呂超(西北農林科技大学大学院生)が記録したものである[14][15]

2004年9月25日原口證が8時間45分かけて円周率5万4000桁の暗唱に成功し、従来の世界記録を更新した。しかしながら、実際はより多くの桁を覚えていたため、2005年7月1日 - 7月2日に再挑戦し、8万3431桁までの暗唱に成功した。2006年10月3日午前9時 - 10月4日午前1時30分(16時間30分)の挑戦で円周率10万桁の暗唱に成功した。ギネス世界記録に申請中である。

文化的影響[編集]

ベルリン工科大学数学科の近くにあるタイル

という日常でもよく知られた図形についての単純な定義でありながら、小数部分が無限に続くという不思議さから、数学における概念の中で最もよく知られたものの一つである。

  • 3月14日円周率の日および数学の日である。小数点以下が「永遠に続く」という意味にあやかり、3月14日に結婚するカップルもいる[16]。また7月22日は円周率近似値の日とされている(22/7 は円周率の近似値)。
  • 2012年8月14日、米国勢調査局が、米国の人口が円周率と同じ並びの3億1415万9265人に達したと発表した。アメリカには円周率の曲を作る人もいる[17]
  • 組版処理ソフトウェア TeX のバージョン番号は、3.14, 3.141, 3.1415, … というように、更新のたびに円周率に近づいていくように一桁ずつ増やされる。

[編集]

小数点以下1000桁までの値

π = 3.

14159 26535 89793 23846 26433 83279 50288 41971 69399 37510 58209 74944 59230 78164 06286 20899 86280 34825 34211 70679 82148 08651 32823 06647 09384 46095 50582 23172 53594 08128 48111 74502 84102 70193 85211 05559 64462 29489 54930 38196 44288 10975 66593 34461 28475 64823 37867 83165 27120 19091 45648 56692 34603 48610 45432 66482 13393 60726 02491 41273 72458 70066 06315 58817 48815 20920 96282 92540 91715 36436 78925 90360 01133 05305 48820 46652 13841 46951 94151 16094 33057 27036 57595 91953 09218 61173 81932 61179 31051 18548 07446 23799 62749 56735 18857 52724 89122 79381 83011 94912 98336 73362 44065 66430 86021 39494 63952 24737 19070 21798 60943 70277 05392 17176 29317 67523 84674 81846 76694 05132 00056 81271 45263 56082 77857 71342 75778 96091 73637 17872 14684 40901 22495 34301 46549 58537 10507 92279 68925 89235 42019 95611 21290 21960 86403 44181 59813 62977 47713 09960 51870 72113 49999 99837 29780 49951 05973 17328 16096 31859 50244 59455 34690 83026 42522 30825 33446 85035 26193 11881 71010 00313 78387 52886 58753 32083 81420 61717 76691 47303 59825 34904 28755 46873 11595 62863 88235 37875 93751 95778 18577 80532 17122 68066 13001 92787 66111 95909 21642 01989 …

注釈[編集]

[ヘルプ]
  1. ^ ただし、これは明らかな根拠がない話であり、適切に表現すれば定まらないというのが正しい、という主張も見られる[5]
  2. ^ 3回の反復で小数18位まで求めることができる

出典[編集]

  1. ^ a b 日本数学会 2007, pp. 94–95.
  2. ^ a b Boeing, Niels (14 March 2016). “Die Welt ist Pi” (ドイツ語). Zeit Online. http://www.zeit.de/zeit-wissen/2016/02/pi-tag-mathematik-pi-kreiszahl. 
  3. ^ a b 「円周率」 - 世界大百科事典 第2版(コトバンク)
  4. ^ Simon Singh The Simpsons and Their Mathematical Secrets
  5. ^ a b c πの文字使用について
  6. ^ 杉浦, 光夫 『解析入門I』 東京大学出版会、1980年、185頁。ISBN 4130620053
  7. ^ Rudin 1976, p. 183.
  8. ^ "An {ENIAC} Determination of pi and e to more than 2000 Decimal Places", Mathematical Tables and Other Aids to Computation, 4 (29), pp. 11–15. (January, 1950)
  9. ^ “y-cruncher - A Multi-Threaded Pi-Program”. http://www.numberworld.org/y-cruncher/ 2015年2月1日閲覧。 
  10. ^ 黒田 2002, p. 176.
  11. ^ A new formula to compute the n'th binary digit of pi - Fabrice Bellard
  12. ^ 円周率の公式集 暫定版 V er:3.141 - 松元隆二
  13. ^ サイモン・シン著、青木薫訳、『フェルマーの最終定理』、新潮社、2000年、ISBN 4-10-539301-4、42ページ
  14. ^ 陝西省の大学院生、円周率暗唱のギネス新記録樹立[リンク切れ]
  15. ^ レコードチャイナ:円周率6万桁以上を暗唱、世界記録に輝いた「記憶の達人」-陝西省楊凌市
  16. ^ 安田美沙子3・14結婚は『円周率=永遠』の意味だった[リンク切れ]スポニチアネックス 2014年3月16日(日)12時17分配信
  17. ^ 米国の人口が円周率と「同じ」に 3億1415万9265人 CNN 2012.08.15 Wed posted at 12:42 JST

参考文献[編集]

関連項目[編集]

数学に関するもの
教育に関するもの
社会に関するもの

外部リンク[編集]