Loading [MathJax]/jax/output/CommonHTML/jax.js

ヤコビの楕円関数3(二重周期性・零点・極)

前回はこちら:

ヤコビの楕円関数2(定義域の拡張・半角公式・倍角公式・展開)

前回はヤコビの楕円関数3つを複素数へ拡張しました。複素関数となったからには、零点や極を知っておきたいところ。また、一般の楕円関数を定義づける二重周期性についても解説します。

復習:ヤコビの楕円関数の複素バージョン

定理35(再掲)

sn(u+iv)=snudnv+icnudnusnvcnvcn2v+k2sn2usn2vcn(u+iv)=cnucnvisnudnusnvdnvcn2v+k2sn2usn2vdn(u+iv)=dnucnvdnvik2snucnusnvcn2v+k2sn2usn2v

ただしダッシュは母数が k という意味です(微分ではない)。左辺の母数はすべて k ですが省略しています。

二重周期性

定理35で v=K を考えます。命題29を参考にすると、例えばdnv=dn(v,k)=dn(K,k)=kであり、同様に snv=1 , cnv=1 となることに注意してsn(u+iK)=ksnuk2sn2ucn(u+iK)=iksnudnuk2sn2udn(u+iK)=ik2snucnuk2sn2uを経て

系36

sn(u+iK)=1ksn(u,k)cn(u+iK)=idn(u,k)ksn(u,k)dn(u+iK)=icn(u,k)sn(u,k)

左辺の母数 k は省略しています。また同様の計算をすることで

系37

sn(u+2iK)=sn(u,k)cn(u+2iK)=cn(u,k)dn(u+2iK)=dn(u,k)

cn(u+4iK)=cn(u,k)dn(u+4iK)=dn(u,k)

よって snu は周期 2iK(すなわち 2K(k), すなわち 2K(1k2))を、cnu は周期 4iK を、dnu は周期 4iK を持つことが分かります。ところでやはり命題29から snu は周期 4K , cnu は周期 4K , dnu は周期 2K をもつのでした。このことから、ヤコビの楕円関数は実数方向および虚数方向に周期をもつことになります。例えばsn(u+8K6iK)=snuといった具合です。

一般にf(z)=f(z+mω1+nω2)m,nZなる ω1,ω2C が存在するとき、関数 f二重周期性をもつといいます。ただし I(ω1/ω2)>0 という条件付きです。複素平面で ω1ω2 という2点を取ったとき、0, ω1, ω2, ω1+ω2 でできる平行四辺形が基本単位で、あとはそれが繰り返されます。

cn に関してだけは cn(u+2K+2iK)=cnu であることに注意すると、

系38

snu は周期 4K , 2iK をもつ。

cnu は周期 4K , 2K+2iK をもつ。

dnu は周期 2K , 4iK をもつ。

つまり snu0, 4K, 2iK, 4K+2iK でできる長方形の中のみを論じれば、周期性によって複素平面全体について分かることになります。

cnu は先ほど見たように 4K, 4iK を周期にもちますが、これを基本単位としてしまうと周期 2K+2iK を含めることができません。逆に系38のようにしておけば、(4K)+2(2K+2iK)=4iKとなって周期 4iK を表すことができます。

零点と極

実数の範囲においては mZsn2mK=0cn(2m+1)K=0が零点で、dnu0 なのでした。

このことから snu の零点と特異点を求めてみましょう。定理35よりsn(u+iv)cnv=0 かつ snu=0 で発散します(cnv=0 なら snv0 なので!)。よって u=2mK かつ v=(2n+1)K となって 2mK+(2n+1)iKsnu の特異点となります。

一方、零点については定理35よりsn(u+iv)=0{snudnv=0cnudnusnvcnv=0dn はゼロになりませんので snu=0 が必要です。つまり u=2mK であり、sn(2mK+iv)=i(1)msnvcnvこのとき snv=0 すなわち v=2nK なら零点を与えます。よって u=2mK+2niKsnu の零点と分かります。定理33よりこの零点は1位です(面倒なので証明は省きますが snu は後述の極を除いて解析的なので、べき展開可能)。

話を戻して先ほどの特異点 u=2mK+(2n+1)iK について。この特異点の性質を調べるには、周期性から m=n=0 のみ論じてよいです。要は sn(u+iK)u0 としたときのようすを探ります。 系36からsn(u+iK)=1ksn(u,k)です。snuu=0 で1位の零点をもつので、sn(u+iK) としては1位の極となります。よって u=2mK+(2n+1)iKsnu の1位の極です。

cnudnu についても同様に考えると

定理39
零点
snu2mK+2niK2mK+(2n+1)iK
cnu(2m+1)K+2niK2mK+(2n+1)iK
dnu(2m+1)K+(2n+1)iK2mK+(2n+1)iK
零点・極はすべて1位

この結果はHancock [4]を参考にしました。

結局、特異点がすべて極なのでヤコビの楕円関数は有理型です。そもそも一般の楕円関数の定義は「2方向に周期をもつ有理型関数」であり、ヤコビの楕円関数はその例ということになります。

零点・極まわりでの展開と留数について概略

零点や極のまわりの展開は複素解析等で必要になります。いまは snz に絞って考えましょう。

周期性より snz の零点は z=0,2K を考えれば十分です。z=0 まわりの展開はマクローリン展開で、すでに導出を終えています。snz=z+O(z3)また前回の定理28を使えばsnz=sn(z2K)=(z2K)+O((z2K)3)です。

極は z=iK,2K+iK の2つを考えれば十分です。iK のローラン展開を考えると留数はa1=limz0zsn(z+iK)=1kここで系36を使いました。2K+iK も系36と定理28を使って同様にすると、留数は 1/k となります。

定理39-b

snz の零点と極まわりの展開はsnz=z+O(z3)snz=(z2K)+O((z2K)3)snz=1k(ziK)+O(1)snz=1k(z2KiK)+O(1)

雑多な公式集

以上で書きたいことは終わりですが、余興としていくつか公式を導出しましたので書き留めておきます。

系32と定理34を再掲すると

系32(再掲)

snK2=11+k,cnK2=k1+k,dnK2=k

定理34(再掲)

sn(iu,k)=isc(u,k)cn(iu,k)=nc(u,k)dn(iu,k)=dc(u,k)

これらより

命題40

sniK2=ikcniK2=1+kkdniK2=1+k

命題40と定理28より

命題41

sn(K+iK2)=1kcn(K+iK2)=i1kkdn(K+iK2)=1k

snu の定義からu=sn(u,k)0dx(1x2)(1k2x2)kx=y と置換してku=ksn(u,k)0dy(1y2)(1y2k2)一方で snu の定義からku=sn(ku,1/k)0dx(1x2)(1x2k2)このようにして次の公式ができあがります。

命題42

sn(ku,1k)=ksn(u,k)cn(ku,1k)=dn(u,k)dn(ku,1k)=cn(u,k)

参考文献

[1] 武部尚志. (2019). 楕円積分と楕円関数 おとぎの国の歩き方. 日本評論社.楽天はココ

楕円積分・楕円関数のことが平易に書かれています。

[2] Whittaker, E. T., & Watson, G. N. (2021). A course of modern analysis. Cambridge University Press.

第5版です。いわずと知れた名著。ヤコビの楕円関数に1章を割いています。

[3] Borwein,J.M., Borwein,P.B. (1987) "Pi and the AGM : a study in analytic number theory and computational complexity"
楕円積分に関する定理がいろいろあります。興味深い式が多く導出されており、一見の価値あり。

[4] Hancock, H. (1917). Elliptic Integrals

次回はこちら:

楕円関数の定義と基本的性質

応援のおねがい

Please support me!

まめしば
まめしば

記事を気に入って下さった方、「応援してあげてもいいよ」という方がいらっしゃったら15円から可能なので支援していただければ幸いです。情報発信を継続していくため、サーバー維持費などに充てさせていただきます。

ご支援いただいた方は、こちらで確認できます。

Amazonギフトの場合、
Amazonギフト券- Eメールタイプ – Amazonベーシック
より、金額は空白欄に適当に(15円から)書きこんで下さい。受取人は「mamekebiamazonあっとgmail.com」です(あっとは@に置き換えてください)。贈り主は「匿名」等でOKです。全額がクリエイターに届きます。

OFUSEは登録不要で、100円から寄付できます。金額の90%がクリエイターに届きます。

codocは登録不要で、100円から寄付できます。金額の85%がクリエイターに届きます。

コメントを残す

メールアドレスが公開されることはありません。 が付いている欄は必須項目です

CAPTCHA