13
高校数学解説
文献あり

正弦関数の極限の循環論法,解決

3377
0

まえがき

こんにちは!的場 沙雪です.この前友達に「円周率が一定であることの証明ってあるんかな?」と質問されて,その時は何も調べず思いついた解法で答えたのですが,その後この方法が有名な循環論法(limh0sinhh=1を導くのに使われる円板の面積公式は三角関数の微積分公式から従い,それはlimh0sinhh=1を使う)を回避出来ることに気が付いたので記事にした次第です.

前提知識は数Ⅲの「曲線の長さ」まで,とあとは…… 全射 単射 くらいは知っておいてほしいですかね(逆関数のところで習ったっけな……?).なるべく高校生でもわかるように書きました(書いたつもり).まあわからない部分は読み飛ばしても大体何とかなります(「極座標」は読み飛ばしてもいいかも).「とりあえずそうなるんだな」で一旦置いといて,後でまた帰ってきましょう.数学って大体そういうもんです.

私は開区間を],[で書きます.特にこの記事では順序対(,)も開区間も多数出現し,同じ記号だと大変紛らわしいので書き分けます.なるべく「開区間],[」と書くので安心していただきたいです.

忙しい人向け

頭の中に時限爆弾がセットされててそれどころじゃない人向けの概略です.本編と多少異なる部分がありますが言ってることは大体同じなので見逃してください.

円の方程式より,原点中心,半径rの円の円周L(r)
L(r)=2rr1+(ddx(r2x2))2dx=2rrrdxr1(xr)2=2r11dt1t2(xr=t)
なので,円周率の定義(円周率は円周と直径の比)よりπ=L(r)2r=11dt1t2です.
よって原点中心,半径rの円板の面積S(r)
S(r)=2rrr2x2dx=rr(1r2x2+x2x2r2x2)dx+r2rrdxr1(xr)2=[xr2x2]x=rr+r211dt1t2(xr=t)=πr2
と求まり,いつもの面積比較でh>0のとき
121sinh<h2πS(1)<121tanhsinh<h<tanhsinh<hh<tanhsinhh<1cosh<sinhhcosh<sinhh<1
なのでlimh0sinhh=1を得て,limh0sinhh=limh0sin(h)h=1と合わせてlimh0sinhh=1を得ます.循環論法は回避出来ました!

円周率

まず,円周率とは何でしょうか?循環論法を避けるためにsinをテイラー展開で定義して,sinの正の零点のうち最小のもの……なんて定義もありますが,うーんわかりにくい!

そもそも「円周率」という言葉の意味は何でしょうか?答えは小学校算数にありますね.「円周と直径の比」でした.ということは「円周が直径の定数倍であることの証明」をすれば自ずと円周率が出てきますね.つまり円周を求めればよいことになります.

ではまず円の定義を考えなければなりません.中心P,半径rrは正実数です)の円とは,「点Pからの(通常はユークリッド)距離がrである点全体の集合」として定義されます.今は二次元平面,xy平面を考えましょう.点P(a,b)とします.点(x,y)と点Pとの距離は,定義から(xa)2+(yb)2です.ということは,円の定義は次のようになります(つまり中心P,半径rの円の方程式は(xa)2+(yb)2=rということです).

正実数r,実数abに対して,
S:={(x,y)R2:(xa)2+(yb)2=r}
により定義されるR2の部分集合Sを「中心(a,b),半径rの円」とよび,円Sの中心は(a,b)であるといい,円Sの半径はrであるといいます.

(xa)2+(yb)2=ryについて解いてみましょう!すると,
(xa)2+(yb)2=r(xa)2+(yb)2=r2(r>0)(yb)2=r2(xa)2yb=±r2(xa)2y=b±r2(xa)2
となります.円の方程式を陽関数で表示することが出来ました!fa,b,r(x):=b+r2(xa)2ga,b,r(x):=br2(xa)2としておきます.

fa,b,rga,b,rの定義域を考えてみましょう.これは根号を含むので,根号の中身が非負であるという条件から定義域がわかります.r2(xa)20xについて解くと,
r2(xa)20(xa)2r2rxararxa+rx[ar,a+r]
となります.定義域は閉区間[ar,a+r]ですね.

ある程度円の概形を把握しておきましょう.fa,b,rga,b,rは関数の形から連続で,
fa,b,r(ar)=b+r2(ara)2=b+r2r2=b+0=b
ga,b,r(ar)=br2(ara)2=br2r2=b0=bより
fa,b,r(ar)=b=ga,b,r(ar)であり,
fa,b,r(a+r)=b+r2(a+ra)2=b+r2r2=b+0=b
ga,b,r(a+r)=br2(a+ra)2=br2r2=b0=bより
fa,b,r(a+r)=b=ga,b,r(a+r)であり,ar<x<a+rを満たすxに対して
fa,b,r(x)=b+r2(xa)2=b+(r(xa))(r+(xa))=b+((a+r)x)(x(ar))>b
ga,b,r(x)=br2(xa)2=b(r(xa))(r+(xa))=b((a+r)x)(x(ar))<bより
fa,b,r(x)>b>ga,b,r(x)
ですから,円は単純閉曲線であることがわかります.一度通った点を通らずに一筆で書いて出発した点に戻ってこれるという意味です.

また
fa,b,r(x)=2(xa)2r2(xa)2=xar2(xa)2
よりa<x<a+rのときfa,b,rxに関して単調減少であり,
fa,b,r(x)=r2(xa)2(xa)2(xa)2r2(xa)2r2(xa)2=r2(r2(xa)2)r2(xa)2<0
なので,上側の陽関数fa,b,r(x)によって描かれる曲線は上に凸です.この結果は後で使います.

円周と円周率

準備がそろいました.では円周の長さを求めてみましょう.中心(a,b),半径rの円の円周の長さはL(a,b,r)としておきます.
L(a,b,r)=ara+r1+(ddxfa,b,r(x))2dx+ara+r1+(ddxga,b,r(x))2dx=ara+r1+(ddx(b+r2(xa)2))2dx+ara+r1+(ddx(br2(xa)2))2dx=ara+r1+(2(xa)2r2(xa)2)2dx+ara+r1+(2(xa)2r2(xa)2)2dx=ara+r1+(xa)2r2(xa)2dx+ara+r1+(xa)2r2(xa)2dx=2ara+r1+(xa)2r2(xa)2dx=2rara+rdxr1(xar)2=2r11dt1t2(xar=t)
となります.なんかよくわからない積分11dt1t2が出てきましたね.被積分関数が端点で発散しています.こういう積分は広義積分といいます.果たして本当に収束するのでしょうか?それとも,円周の長さは無限大!!?
11dt1t2=111t21t2dt+11t21t2dt=111t2dt11t2t21t2dt=111t2dt[t1t2]t=11+111t2dt=2111t2dt
なので,ちゃんと収束しますね(きちんと説明すると,積分区間を開区間]1+ε,1ε[として後でε0の極限を取ることで,閉区間で連続な関数は可積という定理( 理工系の微分積分学 p.106定 理 3)と,可積な関数の不定積分は連続であるという定理( 理工系の微分積分学 p.109定 理 4)が保証してくれます).ということで,円周率L(a,b,r)2rは一定であることがわかりました!これで円周率πL(a,b,r)2r,つまり11dt1t2で定義出来ます!

なお,弧長積分の公式を既知として扱いましたが,その証明に三角関数の極限や微積分の公式は使いません.証明は こちらのサイト や, 理工系の微分積分学 p.123(Ⅱ)を参照してください.

円板の面積

こんどは円板の面積を求めてみましょう.円板の定義もしておきます.

円板

正実数r,実数abに対して,
S:={(x,y)R2:(xa)2+(yb)2r}
により定義されるR2の部分集合Sを「中心(a,b),半径rの円板」とよび,円板Sの中心は(a,b)であるといい,円板Sの半径はrであるといいます.

中心(a,b),半径rの円板の面積はS(a,b,r)としておきます.先程と同じように変形すると,
(xa)2+(yb)2r(xa)2+(yb)2r2(r>0)(yb)2r2(xa)2r2(xa)2ybr2(xa)2br2(xa)2yb+r2(xa)2ga,b,r(x)yfa,b,r(x)
となることに注意してやれば,
S(a,b,r)=ara+r(fa,b,r(x)ga,b,r(x))dx=ara+r((b+r2(xa)2)(br2(xa)2))dx=ara+r2r2(xa)2dx=2r2ara+r1(xar)2dxr(r>0)=2r2111t2dt(xar=t)=r211dt1t2=πr2
となり,円板の面積公式S(a,b,r)=πr2が導出されます.ここで驚いてほしいのです!三角関数の微積分どころか定義さえ無いまま円板の面積公式が導出出来ました!!!これは,循環論法に陥らずにlimh0sinhh=1が導出できることを示唆しています!

三角関数

三角関数を図形的に定義したいので,角度を定義しておかなければなりません.といっても図形の角度まで定義する必要は無く,偏角だけで十分です(もし図形の角度も定義したいのであれば,偏角同士の差の絶対値と,それを2πから引いた値のうち小さい方かなんかで定義すればいいでしょう).

偏角

角度の定義と聞いて,まずラジアンのあの定義が思い浮かびます.「半径1の扇形の円弧の長さが1であるとき,その中心角を1[rad]と定める」です.この定義,実は不十分です.なぜなら1[rad]しか定義されていないし,ラジアンがデシベルやマグニチュードやペーハーや等級のように対数を表した単位である可能性が残っているためです.それに,そもそも角度とは何かを今定義したいのに,これでは「中心角」がなんなのかわかりません.どうしましょう……

そういえば,こんな性質がありました.扇形の弧長をl,半径をr,中心角をθとするとl=rθが成り立ちます.それなら逆にθ:=lrで角度を定義してやるのはどうでしょうか?弧長は先程のように積分を使えば表現出来そうですね.物理でも[rad]はSI組立単位で[m/m]と書くようですし,これでよさそうです.

ということで,偏角を定義していきます.まずxy平面上の点Q(c,d)を考え,これは点Pでないとします(なので(ca)2+(db)2>0です).中心が点Pで点Qを通る円はただ一つ存在します(半径を(ca)2+(db)2で定めるしかないので).この円上の点(a+(ca)2+(db)2,b)から点Qまで反時計回りに沿った弧長を半径(ca)2+(db)2で割って,点Pに対する点Qの偏角とします.

偏角

(a,b)(c,d)を満たすような,実数abcdに対して,
dbθ:=1(ca)2+(db)2ca+(ca)2+(db)21+(ddxfa,b,(ca)2+(db)2(x))2dx
及び
d<bθ:=1(ca)2+(db)2(a(ca)2+(db)2a+(ca)2+(db)21+(ddxfa,b,(ca)2+(db)2(x))2dx+a(ca)2+(db)2c1+(ddxga,b,(ca)2+(db)2(x))2dx)
により定義される実数θを「(a,b)に対する点(c,d)の偏角」とよびます.

極座標

三角関数を図形で定義するためには偏角を変数として扱えなければなりません.つまるところ,点Pを除く直交座標を点Pを極とする極座標に移す関数が逆関数を持たなければなりません.それが証明できれば,この関数の逆関数を用いて三角関数を定義することが出来ます!

まず偏角がどのような値をとり得るのか調べておきましょう.点Pに対する点Qの偏角を計算すると,次のようになります(計算は先程と同じ手順です).

  • dbのときca(ca)2+(db)21dt1t2
  • d<bのときπ+1ca(ca)2+(db)2dt1t2

これをよく見ると,どちらもca(ca)2+(db)2という部分があり,そこ以外にcdは無いです.ということはこれを調べればよさそうですね.Xa,b(c,d):=ca(ca)2+(db)2としておきます.
(i)d=bのとき
簡単な方から考えます.(c,d)(a,b)d=bよりcaであることに注意して,Xa,b(c,d)=Xa,b(c,b)=ca(ca)2+(bb)2=ca|ca|=sgn(ca)となります(sgn 符号関数 とよばれます).ということはXa,b(c,d)11のどちらかのみをとりますね.
(ii)d>bのとき
関数のとり得る値を調べたいので,増減がわかればよいですね.ということでcで微分してみます.
cXa,b(c,d)=(ca)2+(db)2(ca)2(ca)2(ca)2+(db)2(ca)2+(db)2=(db)2((ca)2+(db)2)(ca)2+(db)2>0
より,Xa,b(c,d)cに関して狭義単調増加です.
これとlimcXa,b(c,d)=limcca(ca)2+(db)2=1limcXa,b(c,d)=limcca(ca)2+(db)2=1Xa,b(c,d)cに関して連続であることより,Xa,b(c,d)のとり得る値の範囲は開区間]1,1[です.
(iii)d<bのとき
(ii)と同様にXa,b(c,d)のとり得る値の範囲は開区間]1,1[です.

Θ+(x):=x1dt1t2Θ(x):=1xdt1t2とし(xは閉区間[1,1]の元です)て,Θ+(x)Θ(x)のとり得る値も調べておきます.どちらも被積分関数が端点を除いて正なので,Θ+は狭義単調減少,つまり単射で,Θは狭義単調増加,つまり単射です.また不定積分は連続であるという定理(大学の学部一年で習います)があるので,これよりΘ+Θは連続となります.

Θ+(1)=0Θ+(1)=πであり,Θ+について,Θ+は狭義単調減少かつ連続なので,1<x<1のとき,Θ+(x)のとり得る値の範囲は開区間]Θ+(1),Θ+(1)[,つまり開区間]0,π[です.Θについても同じように,1<x<1のとき,Θ(x)のとり得る値の範囲は開区間]Θ(1),Θ(1)[,つまり開区間]0,π[です.

準備がそろいました.それでは点Pを除く直交座標を点Pを極とする極座標に移す関数φa,bは全単射であることを証明しましょう.関数φa,bを,

  • ybのときφa,b(x,y):=((xa)2+(yb)2,Θ+(Xa,b(x,y)))
  • y<bのときφa,b(x,y):=((xa)2+(yb)2,π+Θ(Xa,b(x,y)))

として定義します.ただし(x,y)(a,b)で定義されていて,終域は第1変数が正実数]0,[,第2変数が半開区間[0,2π[とします.

単射性は簡単なのでまずはこちらから.φa,b(x1,y1)=φa,b(x2,y2)とします.もしy1b>y2なら,第2変数を見比べてΘ+(Xa,b(x1,y1))=π+Θ(Xa,b(x2,y2))となりますが,Θ+(Xa,b(x1,y1))のとり得る値の範囲は閉区間[0,π]なのにπ+Θ(Xa,b(x2,y2))のとり得る値の範囲は開区間]π,2π[です.矛盾ですね.y1<by2も同様に矛盾です(同様というか文字を入れ替えただけなのでそりゃそう).ということはy1y2b以上か,y1y2b未満です.
(I)y1bかつy2bのとき
2変数を見比べてΘ+(Xa,b(x1,y1))=Θ+(Xa,b(x2,y2))となります.Θ+は単射なので,Xa,b(x1,y1)=Xa,b(x2,y2),つまりx1a(x1a)2+(y1b)2=x2a(x2a)2+(y2b)2となります.ここで第1変数を見比べて(x1a)2+(y1b)2=(x2a)2+(y2b)2となります.従って,x1a=x2a,つまりx1=x2を得ます.あとは,これと(x1a)2+(y1b)2=(x2a)2+(y2b)2から(y1b)2=(y2b)2となりますが,y1bかつy2bですから,y1b=y2b,つまりy1=y2を得ます.よって(x1,y1)=(x2,y2)を得ます.
(II)y1<bかつy2<bのとき
(I)と同様にして(x1,y1)=(x2,y2)を得ます.
以上(I)(II)より(x1,y1)=(x2,y2)なので,φa,bは単射です.

次に全射性です.正実数r0以上2π未満の実数θをとります.
(a)θ=0のとき
φa,b(a+r,b)を考えると
φa,b(a+r,b)=((a+ra)2+(bb)2,Θ+(Xa,b(a+r,b)))=(r2,Θ+(Xa,b(a+r,b)))=(r,Θ+(Xa,b(a+r,b)))=(r,Θ+(a+rar))=(r,Θ+(1))=(r,0)=(r,θ)
となり,(x0,y0)(a,b)及びφa,b(x0,y0)=(r,θ)となる実数x0y0は存在することがわかります.
(b)0<θ<πのとき
(x,y)(a,b)y>bのとき,Xa,b(x,y)のとり得る値の範囲は開区間]1,1[です.1<x<1のとき,Θ+(x)のとり得る値の範囲は開区間]0,π[なので,Θ+(Xa,b(x,y))のとり得る値の範囲は開区間]0,π[です.よってΘ+(Xa,b(x,y))=θとなる実数xyは存在し,(x,y)(a,b)y>bを満たします.そのような実数xyに対してφa,b(a+r(xa)(xa)2+(yb)2,b+r(yb)(xa)2+(yb)2)を考えると((a+r(xa)(xa)2+(yb)2,b+r(yb)(xa)2+(yb)2)=(a,b)+r(xa)2+(yb)2((x,y)(a,b))(a,b)+r(xa)2+(yb)2((a,b)(a,b))=(a,b)に注意してください),
φa,b(a+r(xa)(xa)2+(yb)2,b+r(yb)(xa)2+(yb)2)=((a+r(xa)(xa)2+(yb)2a)2+(b+r(yb)(xa)2+(yb)2b)2,Θ+(Xa,b(a+r(xa)(xa)2+(yb)2,b+r(yb)(xa)2+(yb)2)))=(r2(xa)2(xa)2+(yb)2+r2(yb)2(xa)2+(yb)2,Θ+(Xa,b(a+r(xa)(xa)2+(yb)2,b+r(yb)(xa)2+(yb)2)))=(r,Θ+(Xa,b(a+r(xa)(xa)2+(yb)2,b+r(yb)(xa)2+(yb)2)))=(r,Θ+(a+r(xa)(xa)2+(yb)2ar))=(r,Θ+(xa(xa)2+(yb)2))=(r,Θ+(Xa,b(x,y)))=(r,θ)
となり,(x0,y0)(a,b)及びφa,b(x0,y0)=(r,θ)となる実数x0y0は存在することがわかります.
(c)θ=πのとき
φa,b(ar,b)を考えると
φa,b(ar,b)=((ara)2+(bb)2,Θ+(Xa,b(ar,b)))=(r2,Θ+(Xa,b(ar,b)))=(r,Θ+(Xa,b(ar,b)))=(r,Θ+(arar))=(r,Θ+(1))=(r,π)=(r,θ)
となり,(x0,y0)(a,b)及びφa,b(x0,y0)=(r,θ)となる実数x0y0は存在することがわかります.
(d)π<θ<2πのとき
(b)と同様にして(x0,y0)(a,b)及びφa,b(x0,y0)=(r,θ)となる実数x0y0は存在することがわかります.
以上(a)(b)(c)(d)より(x0,y0)(a,b)及びφa,b(x0,y0)=(r,θ)となる実数x0y0は存在するので,φa,bは全射です.

ということで,φa,bは全単射です.従って逆関数φa,b1が存在します.

三角関数とその極限

さあ三角関数を定義しましょう!三角関数はよく単位円で定義するので,そうしておきましょう.中心(0,0),半径1の円x2+y2=1上の点のうち,点(0,0)に対する偏角がθ2πθ2πである点のx座標,y座標をそれぞれcosθsinθとします( 床関数 です.要はガウス記号と同じです).

三角関数

θR, (cosθ,sinθ):=φ0,01(1,θ2πθ2π)
により関数cos:RRsin:RRを定義します.

θ{xR:cosx0}, tanθ:=sinθcosθ
により関数tan:{xR:cosx0}Rを定義します.

φa,bの全射性の証明で,(b)内でφa,b(x0,y0)=(r,θ)なるx0y0としてy0>bになるものを選んでいました.ということは0<θ<πのときsinθ>0であることがわかります.同様にπ<θ<2πのときsinθ<0です.π<θ<0のときsinθ<0もいえます.

三角関数の極限でよく使う公式,0<h<πのときcos(h)=coshかつsin(h)=sinhを証明しておきます(一般の実数hに対しても言えますが,今回は必要ありません.定義からほぼ自明な三角関数の周期性を使えば証明できます).まず定義よりφ0,0(cosθ,sinθ)=(1,θ2πθ2π)であり,第1成分を見比べてcos2θ+sin2θ=1つまりcos2θ+sin2θ=1を得ます.0<h<πのとき,φ0,0(cosh,sinh)=(1,h)の第2成分を見比べてΘ+(cosh)=hφ0,0(cos(h),sin(h))=(1,h+2π)の第2成分を見比べてπ+Θ(cos(h))=h+2πより
Θ(cos(h))=Θ+(cosh)+π=cosh1dt1t2+11dt1t2=1coshdt1t2=Θ(cosh)
であり,Θは単射なのでcos(h)=coshを得ます.よってsin(h)=sin2(h)=1cos2(h)=1cos2h=sin2h=sinhよりsin(h)=sinhを得ます.

limh0sinhh=1を示すのにlimh0cosh=1を使っていましたね.これも証明しておきます.0<h<πのとき,φ0,0(cosh,sinh)=(1,h)の第2成分を見比べてΘ+(X0,0(cosh,sinh))=Θ+(cosh)=hを得ます.よってlimh0Θ+(cosh)=limh0h=0=Θ+(1)であり,Θ+は連続なのでΘ+(limh0cosh)=Θ+(1)となり,さらにΘ+は単射なのでlimh0cosh=1を得ます.limh0cosh=limh0cos(h)=limh0cosh=1と合わせてlimh0cosh=1を得ます.なお今回はめんどくさいので証明しませんが,導関数を求める際に必要になる三角関数の連続性も同様の方法で証明できます.

さあお待たせしました,limh0sinhh=1を証明しましょう!0<h<π2とし,点(0,0)P,点(cosh,sinh)Q,点(1,0)R,点(1,tanh)Sとします(0<h<π2,つまりΘ+(1)<Θ+(cosh)<Θ+(0)より0<cosh<1なので,cosh0であることに注意してください).中心P,半径1の円を考えると,点QRはこの円上にあります.この円上の点Rから点Qまで反時計回りに沿った弧をQRとします.

円の上側の陽関数f0,0,1によって描かれる曲線は上に凸でしたので,cosh<t<1を満たす実数tについて,線分QR上のx座標がtである点のy座標より弧QR上のx座標がtである点のy座標の方が大きいことがわかります.またa<x<a+rのときf0,0,1xに関して単調減少でした.0<cosh<1よりsinhtanhなので,直線QSの傾きは正です.ということは弧QR上のx座標がtである点のy座標より線分QS上のx座標がtである点のy座標の方が大きいことがわかります.これで三角形PQRの面積,扇形PQRの面積,三角形PSRの面積の順に大きくなることがわかり,
121sinh<12coshsinh+cosh1f0,0,1(t)dt<121tanhを得ます.これを変形すると
121sinh<12coshsinh+cosh1f0,0,1(t)dt<121tanhsinh<coshsinh+2cosh11t2dt<tanhsinh<coshsinh+cosh11t2dt+cosh11t2dt<tanhsinh<coshsinh+cosh11t21t2dt+[t1t2]t=cosh1cosh1t2t21t2dt<tanhsinh<coshsinh+cosh111t2dt+[t1t2]t=cosh1<tanhsinh<coshsinh+hcosh1cos2h<tanhsinh<h<tanhsinh<hh<tanhsinhh<1cosh<sinhhcosh<sinhh<1
なので,limh0cosh=1とはさみうちの原理よりlimh0sinhh=1です!limh0sinhh=limh0sin(h)h=limh0sinhh=1と合わせてlimh0sinhh=1が得られました!!!!!!

あとがき

いや~長かったですね(ほとんど極座標のせいですが).これにて循環論法問題解決です!これで例の大阪大学の問題が解けますね()

一度諦めても,ふとしたことで解決策が出ることもあるものです.過去に悩んでいたことを思い返してみると今なら解決できるかもしれないですね!ここまで読んでくださりありがとうございました.

参考文献

[1]
吹田 信之,新保 経彦, 理工系の微分積分学, 学術図書出版社, 1987, p.106,p.109,p.123
投稿日:2021327
OptHub AI Competition

この記事を高評価した人

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

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

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

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

投稿者

微分積分学,数理論理学,順序数解析が好きです.ここでは主に微積や級数の話題をすると思います.記事まとめは下のリンクからどうぞ.

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中
  1. まえがき
  2. 忙しい人向け
  3. 円周率
  4. 円周と円周率
  5. 円板の面積
  6. 三角関数
  7. 偏角
  8. 極座標
  9. 三角関数とその極限
  10. あとがき
  11. 参考文献