6

1-1/3+1/5-1/7+…=π/4の視覚化

259
0

はじめに

 この記事では、先日私が X(旧Twitter) にポストした次の画像について解説します。


1-1/3+1/5-1/7+…=π/4の視覚化 1-1/3+1/5-1/7+…=π/4の視覚化

この図を作ったモチベーション

 みなさんはグレゴリー・ライプニッツ級数とかライプニッツの公式と呼ばれる次の式をご存じでしょうか?

グレゴリー・ライプニッツ級数(ライプニッツの公式)

113+1517±=π4

 規則的に分数を足したり引いたりするだけの簡単な式で円周率を計算できるという、初見では信じがたい形をした級数です。
 私は中学生のときに初めてこの式を見て衝撃を受けました。「こんな簡単な式で円周率が計算できるわけないやん!」と・・・
 覚えたばかりのBASICで計算してみて、確かに 3.14 に収束しそうなのを見て「ほんまや・・・」となってから、私はこの級数が大好きです。

 ただし、世間での扱いはそれほどではないように思います。その理由は私の想像では、収束の遅さがハンパないからではないかと思います。
 例えば第 100 万項あたりまで計算しても、誤差約 200 万分の 1 の幅で上下にフラフラしてるわけですから、桁数としては小数点以下 6 ケタがやっと。
 そのことは、式の形をみればなんとなくわかると思います。
 この収束の遅さ、他のいろいろな円周率公式と比べても圧倒的に遅い。

 しかし、収束の遅さを差し引いても、この見た目の美しさは最高だと思います。「出来の悪い子ほどカワイイ」みたいなものかもしれませんが。

 ともかく、私はこの式の「見た目」が好きなのです。しかし、世に広まっているこの式の「証明」は、この「見た目」とあまりにかけ離れているものばかりだと感じていました。「もっと視覚的に納得したい!」と。「1313 のままで、1515 のままで直接この式を感じたい!」と。

 というわけで作ってみたのが先の図というわけでした。
 先の図では、「1tan(π/4)=1」から始めて、項をどんどん 2 つずつに「分割」していき、分割を繰り返した極限から 1315 が現れ、ライプニッツの公式に至る様子を視覚化しています。

ざっくり導出

 それではここから、どのようにして「1tan(π/4)=1」からライプニッツの公式を導出したのか説明します。

 まず、tan についての次の補題を見てください。

tan(θ)=1tan(π2θ)

tanの関係 tanの関係

tan(θ)=sin(θ)cos(θ)=cos(π/2θ)sin(π/2θ)=1tan(π2θ)

 補題 1 を使って次の補題 2 を示します。

1tan(2θ)=12tan(θ)12tan(π/2θ)

1tan(2θ)=cos(2θ)sin(2θ)=cos2(θ)sin2(θ)2sin(θ)cos(θ)=cos(θ)2sin(θ)sin(θ)2cos(θ)=12tan(θ)tan(θ)2=12tan(θ)12tan(π/2θ)

 補題 2 を使って、1tan(π/4) を「分割」します。

1=1tan(π/4)=12tan(π/8)12tan(3π/8)

 もう一回「分割」すると

12tan(π/8)12tan(3π/8)(=(122tan(π/16)122tan(7π/16))(122tan(3π/16)122tan(5π/16)))=122tan(π/16)122tan(3π/16)+122tan(5π/16)122tan(7π/16)

 こうなります。
 分割するたびに項数が 2 倍になり、分母の数が規則的に大きくなっていることがわかります。
 
 このように次々と「分割」することで、次々と新しい式が得られます。
 どの式も値は 1 であることに注意してください。

1tan(π/4)=1

12tan(π/8)12tan(3π/8)=1

122tan(π/16)122tan(3π/16)+122tan(5π/16)122tan(7π/16)=1

12ktan(π/2k+2)12ktan(3π/2k+2)±12ktan((2k+11)π/2k+2)=1

 このような「分割」が常に可能であることを式で確認しておきましょう。
 k を非負整数とし、

 Ik=n=12k(1)n+12ktan((2n1)π/2k+2)

とします。Ik=1 の左辺を「分割」して Ik+1=1 の式が得られることを確認します。

Ik=n=12k(1)n+12ktan((2n1)π/2k+2)=n=12k((1)n+12k+1tan((2n1)π/2k+3)(1)n+12k+1tan((2k+2(2n1))π/2k+3))=n=12k(1)n+12k+1tan((2n1)π/2k+3)n=2k+12k+1(1)n+12k+1tan((2n1)π/2k+3)=n=12k+1(1)n+12k+1tan((2n1)π/2k+3)=Ik+1

と変形できますので、この「分割」は何回でも繰り返すことができることが確認できました。

 さて、この分割を繰り返していくと、各項が特定の値へと収束していきます。収束先を式で確認してみましょう。Ik の第 n 項を ak,n とすると

 ak,n=(1)n+12ktan((2n1)π/2k+2)

で表すことができます。
※ n>2k のときは項が存在しませんが、ここでは k の極限を考えているので気にしないことにします。

limx0xtanx=1 を使って変形すると

 limkak,n=limk(1)n+12ktan((2n1)π/2k+2)=limk2k+2(2n1)π(1)n+1(2n1)π/2k+22ktan((2n1)π/2k+2)=(1)n+122(2n1)π1=(1)n+14(2n1)π

となりました。したがって、Ik の収束先を具体的に書き下すと

4π43π+45π47π±

となります。Ik=1 と合わせて次の式が得られました。

4π43π+45π47π±=1

 両辺に π4 を乗じることで所望の式が得られました!

グレゴリー・ライプニッツ級数(ライプニッツの公式)(再掲)

113+1517±=π4

もう少し厳密に

 ここまでの議論がなんだか厳密でなさそうでだまされてる気がする、と感じる人のためにもう少し厳密にこの式をみていくことにしましょう。

 Ik=n=12k(1)n+12ktan((2n1)π/2k+2)

の両辺を π4 倍したものを Ik とします。

 Ik=n=12k(1)n+1π2k+2tan((2n1)π/2k+2)

このとき、k にかかわらず必ず Ik=1 ですから、必ず Ik=π4 となることに注意します。
つぎに、ライプニッツの公式の第 2k 項までの部分和を Jk とします。

 Jk=n=12k(1)n+12n1

このとき、limk(IkJk)=0 となることを式で確認しましょう。

 まず、積分を使うことで、誤差項を評価することを考えます。k>0 のとき

π4=0π/4tan0xdx=0π/4tan0xdxcos2x0π/4(1cos2x1)tan0xdx=01u0du0π/4tan2xdx=10π/4tan2xdx=1(0π/4tan2xdxcos2x0π/4(1cos2x1)tan2xdx)=1(01u2du0π/4tan4xdx)=113+0π/4tan4xdx=113+150π/4tan6xdx=113+1517+0π/4tan8xdx=113+1517+12k+11+0π/4tan2k+1xdx

したがって

 Jk=Ik+0π/4tan2k+1xdx

積分の項を評価すると
  0<0π/4tan2k+1xdx<0π/4(4πx)2k+1dx=01u2k+1π4du=π4(2k+1+1)

なのではさみうちの原理により
 limk0π/4tan2k+1xdx=0

したがって

 limk(IkJk)=0

となることが確認できました!

おわりに

 というわけで、私の好きなライプニッツの公式の視覚化の解説でした。
 もう一度、最初の図を見てください。

1-1/3+1/5-1/7+…=π/4の視覚化 1-1/3+1/5-1/7+…=π/4の視覚化

 ……。いかがでしょうか。うまく言語化できないのであえて語りませんが、私がこの図にこめた思いが少しでも伝わっていればいいなあ、と思います。

 しかし、これがベストな視覚化だとは思っていません。
 もっと直接的に「1313 のままで、1515 のままで直接この式を感じられる」ような視覚化がどこかにありそうな気がしてなりません。
 何かアイデア等ありましたら教えていただければ幸いです!

投稿日:20231217
OptHub AI Competition

この記事を高評価した人

高評価したユーザはいません

この記事に送られたバッジ

バッジはありません。
バッチを贈って投稿者を応援しよう

バッチを贈ると投稿者に現金やAmazonのギフトカードが還元されます。

投稿者

apu_yokai
apu_yokai
489
67715

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中
  1. はじめに
  2. この図を作ったモチベーション
  3. ざっくり導出
  4. もう少し厳密に
  5. おわりに