円周率の近似のソースを表示

{{Expand English|Approximations of π|date=2023-5}}
{{円周率}}
本記事では、[[数学定数]]のひとつである'''円周率の近似'''（えんしゅうりつのきんじ）について詳述する。

[[円周率]] ''π'' は[[無理数]]であるため、[[小数]]部分は循環せず無限に続く。さらに、円周率 ''π'' は[[超越数]]でもあるため、その[[連分数]]表示は循環しない。その[[近似値]]は何千年にも亘り世界中で計算されてきた。

== 紀元前5000年ごろ～紀元前3世紀ごろまで ==

=== 古代の円周率 ===
人類は少なくとも紀元前5000年ごろから、車輪のような円形の物体を使って重いものを運ぶなど、円を役立つ形として認識し、利用してきた。特に円に関して高い数学的知識を持っていたのが、現在のイラク南部に当たる地域に紀元前2500年ごろから住んでいたバビロニア人。彼らは、円の研究を進め、円周の長さが円の直径に比例すること、すなわち円周率の存在に気付いていた。

=== 古代エジプトでの円周率 ===
古代エジプト人は円周率を3.125とした。彼らの結論にたどり着くには、簡単な実験をするだけでよい。まず、ひもと木の棒2本を用意し、ひもの両端に木の棒をくくりつける。そして、コンパスのように円を描く。そのあと、そのままひもを定規の代わりにし、半径の何倍になる調べる。そうすると約6.25倍になることがわかる。あとは2で割ると直径の約3.125倍になることがわかる。

<math>\pi=3.125</math>

=== 古代バビロニアでの円周率 ===
古代バビロニア人は円周率を3.16049...とした。彼らの結論にたどり着く方法は、彼らが残した記録を見ることでわかる。

<math>\pi=3.16049\cdots</math>

== 紀元前3世紀～16世紀前半まで ==

=== アルキメデスの手法 ===
紀元前3世紀に入ると、円周率の真の値に限りなく近づくことができる画期的な方法を考え出した人が現れた。古代ギリシャの数学者で物理学者の[[アルキメデス]]。彼の考えだした方法は次の通り。
 円に内接する正六角形の外周(=3)<直径1の円の外周(=π)<円に外接する正六角形の外周(=3.4661...)　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　この円に内接する正多角形と円に外接する正多角形の辺の数を無限に増やせば円周率が求められる。
{| class="wikitable"
|+アルキメデスの手法で求めた円周率の記録
!人物
!時代
!記録
|-
|アルキメデス
|紀元前287年～212年
|小数点以下2桁
|-
|祖 沖之
|430年～501年
|小数点以下7桁
|-
|ルドルフ・ファン・ケーレン
|1540年～1610年
|小数点以下35桁
|-
|関 孝和
|1642年～1708年
|小数点以下10桁
|-
|建部 賢弘
|1664年～1739年
|小数点以下40桁
|}
アルキメデスの方法では正96角形を使って3.14までを決定しているがこれを見るとかなり精度が悪いことがわかる。また、[[東京大学]]の「<math>\pi>3.05</math>を示せ。」という問題も正8角形や正12角形を用いて示すことができる。また、周でなく面積でも求めることができるが、こちらの方が精度が悪く、面積を使って<math>\pi>3</math>を示すためには正12角形を用いる必要がある。<math>\pi>3.05</math>を示すためには正24角形などを用いる必要がある。

== 16世紀なかば～16世紀末まで ==

=== ヴィエトの公式 ===
<math>\frac{2}{\pi} = \prod_{n=1}^{\infty} \cos \frac{90^\circ}{2^n} = \sqrt{\frac{1}{2}}\sqrt{\frac{1}{2}+\sqrt{\frac{1}{2}}}\sqrt{\frac{1}{2}+\sqrt{\frac{1}{2}+\sqrt{\frac{1}{2}}}}\cdots</math>

[[フランソワ・ビエト]]による。アルキメデスの方法と本質的に同じで[[三角関数の公式の一覧#倍角公式|三角関数の半角の公式]]を用いて導出することができる。

=== ウォリスの公式 ===
<math>\frac{\pi}{2}=\frac{2\times2\times4\times4\times6\times6\cdots}{1\times3\times3\times5\times5\times7\cdots}</math>

[[正弦関数]]の[[三角関数の無限乗積展開|無限乗積展開]]を用いて導出することができる。

=== ブラウンカーの公式 ===
=== <math>\pi={\cfrac  {4}{1+{\cfrac  {1^{2}}{3+{\cfrac  {2^2 }{5 +{\cfrac  {\cdots}{\cdots+{\cfrac  {n^{2} }{(2m+1)+\cdots}} }}}}}}}}</math> ===
ウォリスの公式を変形することで導出できる。

== 17世紀～現代 ==
17世紀からは無限級数(無限和)を使った円周率を求める方法が考え出された。

=== 1.arctanを使った公式 ===

==== マーダヴァ・グレゴリー・ライプニッツ(Madhava-Gregory-Leibniz)級数 ====
[[ライプニッツ]]は[[arctan]]の[[マクローリン展開]]

<math>\arctan{x}=\sum^\infty_{k=0}\frac{x^{2k+1}} {2k+1}=x-\frac{x^3} {3}+\frac{x^5} {5}-\cdots</math>

に<math>x=1</math>を代入し、

<math>\arctan{1}=\frac{\pi}{4}=\frac{1}{1}-\frac{1}{3}+\frac{1}{5}-\frac{1}{7}+\frac{1}{9}-\frac{1}{11}+\frac{1}{13}\cdots</math>

を得た([[ライプニッツの公式]])。

arctanのマクローリン展開はすでに[[マーダヴァ]]が、ライプニッツの公式は独立に[[ジェームス・グレゴリー|グレゴリー]]が発見していたため、マーダヴァ・グレゴリー・ライプニッツ級数と呼ばれるようになった。

==== マチン(Machin)の公式 ====
[[ジョン・マチン]]はarctanのマクローリン展開から

<math>\frac{\pi}{4}=4\arctan{\frac{1}{5}}-\arctan{\frac{1}{239}}</math>

を得た。この公式は収束が速く、マチン自身も円周率を100桁求めている。また、[[ウィリアム・シャンクス]]は707桁目まで求めたが、その後<math>\arctan{\frac{1}{5}}</math>が527桁目までしか一致しておらず、正しかったのは527桁目までということがわかった。

==== クリンゲンシュティルナ(Klingenstierna)の公式 ====
<math>\frac{\pi}{4}=8\arctan{\frac{1}{10}}-\arctan{\frac{1}{239}}-4\arctan{\frac{1}{515}}</math>

==== ガウス(Gauss)の公式 ====
<math>\frac{\pi}{4}=12\arctan{\frac{1}{18}}+8\arctan{\frac{1}{57}}-\arctan{\frac{1}{239}}</math>

==== 高野喜久雄の公式 ====
<math>\frac{\pi}{4}=12\arctan{\frac{1}{49}}+32\arctan{\frac{1}{57}}-5\arctan{\frac{1}{239}}+12\arctan{\frac{1}{110443}}</math>

==== オイラー(Euler)の公式 ====
<math>\frac{\pi}{4}=\arctan{\frac{1}{2}}+\arctan{\frac{1}{3}}</math>

<math>\frac{\pi}{4}=5\arctan{\frac{1}{7}}+2\arctan{\frac{3}{79}}</math>

2つ目の公式はオイラー変換と呼ばれる変換でarctanを無限乗積展開すると<math>\frac{2}{100}</math>や<math>\frac{3792}{100000}</math>が現れるため、[[十進数]]での計算がしやすい。実際、[[レオンハルト・オイラー|オイラー]]自身もこの公式を使って1時間程度で20桁求めたといわれている。

==== ストーマー(Störmer)の公式 ====
<math>\frac{\pi}{4}=6\arctan{\frac{1}{8}}+2\arctan{\frac{1}{57}}+\arctan{\frac{1}{239}}</math>

<math>\frac{\pi}{4}=44\arctan{\frac{1}{57}}+7\arctan{\frac{1}{239}}-12\arctan{\frac{1}{682}}+24\arctan{\frac{1}{12943}}</math>

==== ヴェガ(Vega)の公式 ====
<math>\frac{\pi}{4}=2\arctan{\frac{1}{2}}-\arctan{\frac{1}{7}}</math>

<math>\frac{\pi}{4}=4\arctan{\frac{1}{5}}-2\arctan{\frac{1}{408}}+\arctan{\frac{1}{1393}}</math>

==== クラウゼン(Clausen)の公式 ====
<math>\frac{\pi}{4}=2\arctan{\frac{1}{3}}+\arctan{\frac{1}{7}}</math>

==== ダース(Dahse)の公式 ====
<math>\frac{\pi}{4}=\arctan{\frac{1}{2}}+\arctan{\frac{1}{5}}+\arctan{\frac{1}{8}}</math>

==== ラザフォード(Rutherford)の公式 ====
<math>\frac{\pi}{4}=4\arctan{\frac{1}{5}}-\arctan{\frac{1}{70}}+\arctan{\frac{1}{99}}</math>

=== 2.ラマヌジャン型公式…複雑だが、収束が速い公式 ===

==== ラマヌジャン(Ramanujan)の公式 ====
<math>\frac{1}{\pi}=\frac{2\sqrt{2}}{99^2}\sum_{n=0}^\infty \frac{(4n)!(1103+26390n)}{(4^n99^nn!)^4}</math>

<math>\frac{4}{\pi}=\sum_{n=0}^\infty \frac{{(-1)}^n(21460n+1123)(4n!)}{{882}^{2n+1}{(4^{n}n!)}^4}</math>

この級数は1項ごとに約8桁ずつ正確な桁が増える。

==== チュドノフスキー(Chudnovsky)級数 ====
<math>\frac{1}{\pi}=12\sum_{n=0}^\infty \frac{{(-1)}^{n}(6n)!}{(3n)!{(n!)}^3} \frac{13591409+545140134n}{{640320}^{3n+\frac{3}{2}}}</math>

この級数は1項ごとに約14桁ずつ正確な桁が増える。

=== 3.反復公式 ===

==== ガウス=ルジャンドルのアルゴリズム ====
「[[ガウス＝ルジャンドルのアルゴリズム]]」も参照

は{{Math|2}}つの数値の[[算術幾何平均]]を求めるために、それぞれの数値を[[算術平均]](相加平均)と[[幾何平均]](相乗平均)で置き換えていくものである。

===== 初期値の設定 =====
: <math>a_0 =1\qquad b_0 =\frac{1}{\sqrt{2}} \qquad t_0 =\frac{1}{4} \qquad p_0 =1</math>

===== 反復式 =====
{{Math|''a'',&nbsp;''b''}} が希望する精度（桁数）になるまで以下の計算を繰り返す。小数第{{Math|''n''}}位まで求めるとき　{{Math|log{{sub|2}}&thinsp;''n''}} 回程度の反復でよい。

: <math>\begin{align}
a_{n+1} &=\frac{a_n +b_n}{2} \\
b_{n+1} &=\sqrt{a_n b_n} \\
t_{n+1} &=t_n -p_n(a_n - a_{n+1})^2 \\
p_{n+1} &= 2p_n
\end{align}</math>

===== ''{{π}}'' の算出 =====
円周率 ''{{π}}'' は、{{Math|''a'',&nbsp;''b'',&nbsp;''t''}} を用いて以下のように近似される。

: <math>\pi \approx \frac{(a+b)^2}{4t}</math>

最初の3回の反復で得られる数値（最後の桁は真値とは異なる）は以下の通りである。

: <math>3.140\dots</math>　（小数点以下{{Math|2}}桁目までが正しい）
: <math>3.14159 264\dots</math>　（小数点以下{{Math|7}}桁目までが正しい）
: <math>3.14159 26535 89793 2382\dots</math>（小数点以下{{Math|18}}桁目までが正しい）

この計算過程は二次収束する。ガウス自身もこの式を用いて反復を{{Math|4}}回まで行って{{Math|12}}桁まで正しいことを確認したことが知られている。

== その他の近似方法 ==
<math>\sin{x}<x<\tan{x} (0<x<\frac{\pi}{2})</math>

を用いる近似方法もある。例えば

<math>\sqrt{\frac{1}{2}}<\frac{\pi}{4}<1</math>

より

<math>2\sqrt{2}<\pi<4</math>

がわかる。これで四角形の周を利用したものと同じ評価が得られた。

== 記録 ==
Google社の[[岩尾エマはるか]]が[[2019年]][[3月14日]]に円周率を31兆4159億2653万5897桁計算したと発表した。所要時間は約120日。岩尾エマはるかはその後100兆桁計算したと発表。さらに2024年3月、StorageReviewの編集者らが105兆桁まで計算した。所要時間は約70日。いずれもChudnovskyの公式が使われている。

== 関連項目 ==
*[[円周率が22/7より小さいことの証明]]
*[[円周率の日]]

{{DEFAULTSORT:えんしゆうりつのきんし}}
[[Category:数学に関する記事]]
[[Category:近似アルゴリズム]]
[[Category:円周率]]
[[Category:超越数]]