1
大学数学基礎解説
文献あり

楕円関数の変換の解析的な解

150
0

はじめに

 この記事では後の記事に向けて楕円関数の変換という問題の解き方について解説していきます。
 また参考文献の少ない中筆者が自力で考察した内容となっているので一部間違い等があるかもしれませんがご了承ください。

楕円関数の変換

 楕円関数論の祖の一人であるヤコビはその著書"Fundamenta Nova Theoriæ Functionum Ellipticarum"の冒頭で
dyA+By+Cy2+Dy3+Ey4=dxA+Bx+Cx2+Dx3+Ex4
を満たすようなxの有理関数
y=y(x)=U(x)V(x)(U(x),V(x))
を決定するという問題を考えました。そしてヤコビはこの問題を次のような形に帰着させます。

楕円関数の変換

 微分方程式
dy(1y2)(1l2y2)=dxM(1x2)(1k2x2)(M)
を満たすような有理関数yを考える問題を楕円関数の変換と呼ぶ。またこの定数Mのことを乗法子(Multiplier)と言う。

 以下でこの問題を、いくつかの簡単な仮定を設けながら、解析的に解いていきます。

問題の言い換え

 まずx=0の付近で
0ydt(1t2)(1l2t2)=1M0xdt(1t2)(1k2t2)
が成り立ってくれると嬉しいのでy(0)=0を仮定します。このとき以下の主張が成り立ちます。

 楕円関数の変換はsn(u/M,l)sn(u,k)の有理式で表す問題と言い換えられる。

 仮定よりx=0の付近で
u=0xdt(1t2)(1k2t2),v=0ydt(1t2)(1k2t2)
つまり
x=sn(u,k),y=sn(v,l)
という置き換えは可逆となり、
du=dx(1x2)(1k2x2)
に注意すると楕円関数の変換は
dv=duM
つまり
v=uM+C
という関係に置き換わる。
またy(0)=0よりu=0のときv=0となるのでC=0つまりv=u/Mを得る。
 この変換によって微分方程式は満たされるので残る問題はy=sn(u/M,l)x=sn(u,k)の有理式として表すこととなる。

 いまy(0)=0からsn(u,k)sn(u/M,l)の零点の位置を考えると、格子2KZ+2iKZは格子M(2LZ+2iLZ)に埋め込まれる、つまり以下の主張が成り立つことになります。

 ある整数行列AM2(Z)があって
1M(iKK)=A(iLL)
が成り立つ。またA=pA1とおくと
M(iLL)=Ap(iKK)
が成り立つ。この関係によって得られる変換をp次の変換と言う。

モジュラー方程式

 p次の変換
M(iLL)=Ap(iKK)
に対してこの比を取ることで得られるk,lについての関係式
iLL=aiK+bKciK+dK(adbc=p)
p次のモジュラー方程式と言う。
 k=k(τ),l=k(τ)とおくと この記事 の定理10より
τ=aτ+bcτ+d=Aτ
k,lの関係が決定する。特にτ,τHよりdetA=adbc>0でなければならないことがわかる。

 逆にモジュラー方程式が成り立っているとき
M=ciK+dKpL=aiK+bKpiL
とおくことでy=sn(u/M,l)x=sn(u,k)の関係は決定します(ただしそれが有理的な関係であるかはまだわかりません)。
 これによって楕円関数の変換は行列AM2(Z)に対して
sn[uM,Aτ]=fA(sn[u,τ],τ)(sn[u,τ]=θ3(τ)θ2(τ)θ1(u/πθ3(τ)2,τ)θ4(u/πθ3(τ)2,τ))
が成り立つような有理式fA(x,τ)を考える問題に置き換えることができます。

 ちなみにA=(n00n)とおくとAτ=nτ/n=τなので
M=nKn2L=1n
となります。つまり楕円関数の変換はsnの倍数公式
sn[nu,τ]=fA(sn[u,τ],τ)
を一般化した問題だということができます。
 cosの倍角を多項式fによって
cosnx=f(cosx)
と表す問題(cf. チェビシェフ多項式)の類似と考えるとわかりやすいと思います。

 以下(A)τ=AτおよびMA=MAよりfA=fAとなるのが面倒なので諸々の符号についてはあまりよく考えないこととします。

補題その1

 まず楕円関数の変換y=y(x)が存在したとき、yAが満たす性質について考えます。

 yはある互いに素な多項式P,Q(P(0)=Q(0)=1)を用いて
y=xMP(x2)Q(x2)
と表せる。

 x=sn(u,k),y=sn(u/M,l)という置き換えからyxについての奇関数となる。よって
y=U(x)V(x)=y(x)y(x)2=U(x)V(x)U(x)V(x)2V(x)V(x)
となるのでこの分母分子の係数に注意するとある互いに素な多項式P,Qがあって
y=AxP(x2)Q(x2)
と書けることがわかる(y(0)=0よりQ(0)0に注意する)。
 また問題の微分方程式においてx=0とおくことで
dydx(0)=1M=AP(0)Q(0)
がわかるので主張を得る。

 いまA=(abcd)とおいたとき
1M(iKK)=(dbca)(iLL)
が成り立つことに注意します。

 a1,b0(mod2)が成り立つ。

 u=Kにおける挙動を考えると
KM=aLciL
よりy(1)の候補としては0,±1,±1/k,がある。
 y(1)=0,とするとyは分母か分子に1xを因数に持つことになるが、1sn(u,k)は重根を持つのでsn(u/M,l)は二位以上の零点・極を持たないことに矛盾する。
 よってa1(mod2)を得る。

 またu=iKにおける挙動を考えると
sn(u/M,l)=sn(u,k)degUdegVO(1)(uiK)
および偶奇の違いから
degUdegV0
が成り立つのでsn(u/M,l)u=iKにおいて零点か極を取ることになる。
 よって
iKM=bL+diL
からb0(mod2)を得る。

補題その2

 次に行列Aに対してfAを決定するための補題を示しておきます。

 行列A,BM2(Z)に対し
fAB(x,τ)=fA(fB(x,τ),Bτ)
が成り立つ。特に任意のγΛに対し
fγA(x,τ)=fA(x,τ)
fAγ(x,τ)=fA(x,γτ)
となる。

sn[uMA(Bτ)MB(τ),ABτ]=fA(sn[uMB,Bτ],Bτ)=fA(fB(sn[u,τ],τ),Bτ)
より
fAB(x,τ)=fA(fB(x,τ),Bτ)
がわかる。
 また 前の記事 の命題13として確認したように
sn[u,τ]=sn[u,γτ]
つまり
fγ(x,τ)=x
が成り立つことと合わせて主張を得る。

 任意の複素数cに対し方程式
sn(u,k)c=0
は基本領域
{2Ks+iKts,t(1,1]}
内に重複度込みで丁度2つの解をもち、その解をu1,u2とおくと
u1+u2=2Kまたはu1+u2=2K
が成り立つ。特にc±1のときは重解を持たない。

 リウヴィルの第三定理( この記事 の命題5)から楕円関数sn(u,k)cの基本領域における零点と極の個数は等しく、極の個数はu=iK,2K+iK2つのみであるので零点の個数も丁度2つとなる。
 また
sn(u,k)=sn(±2Ku,k)
が成り立つことから
u1+u2=±2K
を得る。

 このことからsn(u,k)=cC上で解u=u1,u2を持つとき、
u1u24KZ+2iKZまたはu1+u22KZ+2iKZ
が成り立つことにも注意します。

素数次の変換

 命題4での議論から
y(1)=±1,±1/k
が成り立つので、簡単のため条件y(1)=±1を追加して考えます。これはc0(mod2)によって満たされます。
 このとき
γA=(abcd)j(adbc=p,a1,b,c0(mod2))
であったので 前の記事 の命題7よりあるγΛによって
γA=(a2j0d)(ad=p,a1(mod2),0j<d)
と変形できます。さらに
γA=(100d)(12j01)(a001)
と分解することで命題5から素数次の変換を考えれば任意の楕円関数の変換が得られます。
 またAの符号を調整することで以下a1(mod4)つまりy(1)=1とします。

2次の変換

 a1(mod2)に注意すると2次の変換は
A=(1002)
の一通りとなる。このとき
M(iLL)=(12001)(iKK)
より基本領域
{2Ks+iKts,t(1,1]}
におけるsn[u/M,τ/2]
零点はu=0,2K,iK,2K+iK
極はu=±iK/2,2K±iK/2
4つずつとなる。よって補題6や 前の記事 の定理6系に注意して
f2(x,τ)=x11/sn2[iK/2,τ]1x2/sn2[iK/2,τ]=(1+k)x1+kx2
とおくと
sn[u/M,τ/2]f2(sn[u,τ],τ)
は零点も極も持たない楕円関数となる。よってリウヴィルの第一定理( この記事 の命題3)からこれは定数関数であり、f2(1,τ)=1よりf2が求める有理関数となる。
 また命題3や 2次モジュラー方程式 から以下の主張が得れらる。

 2次の変換は
f2(x,τ)=(1+k)x1+kx2
となる。このとき
M=11+k,l=2k1+k
が成り立つ。

p次の変換(p:奇素数)

 p次の変換は
Ap=(p001),Aj=(12j0p)(0j<p)
p+1通りある。
 いまjが奇数のとき
(1201)Aj=(12(j+p)01)
によってb0(mod4)とみなすとp=adbcd1(mod2)および
sn(u+iK,k)=1ksn(u,k)iKM=bL+diL
に注意すると
fAj(1kx,τ)=1lfAj(x,τ)
が成り立つ。このことから
fAj(x,τ)=xMm(1x2/am2)n(1x2/bn2)
とおくとk2an2bn2=1つまり
fAj(x,τ)=xMn1x2/an21k2an2x2
となることがわかる。
 そして4KZ+2iKZ=4KZ+(2iK+4jK)Zを法としたsn[u/M,Ajτ]の零点は

  • j=pのとき
    u=2nKp(p<np)
    つまり|n|(p+1)/2のときを少しずらすことで
    u=2nKp,2K2nKp(|n|p12)
    と表せ、
  • jpのとき
    u=2mK+4jK+2iKpn(m=0,1,|n|p12)
    つまり
    u=4jK+2iKpn,2K4jK+2iKpn(|n|p12)
    と表せるので

ωp=2Kp,ωj=4jK+2iKp(jp)
とおくことで
fAj(x,τ)=xMn=1p121x2/an21k2an2x2(an=sn[ωjn,τ])
を得る。
 またfAj(1,τ)=1fAj(1/kx,τ)=1/lfAj(x,τ)としていたことから以下の主張が得られる。

ωp=2Kp,ωj=4jK+2iKp(jp)
とおくと、p次の変換(p:奇素数)は
fj(x,τ)=xMn=1p121x2/an21k2an2x2(an=sn[ωjn,τ])
p+1通りとなる。このとき
M=n=1p12(1an2)an2(1k2an2)=(1)p12n=1p12cn2[ωjn,τ]sn2[ωjn,τ]dn2[ωjn,τ]l=M2kn=1p12k2an4=kpn=1p12cn4[ωjn,τ]dn4[ωjn,τ]
が成り立つ。

ヤコビの解

 ヤコビは少なくとも一方はpと互いに素な整数m,mに対し
ω=4mK+4miKp
とおき
f(x,τ)=xMn=1p121x2/an21k2an2x2(an=sn[ωn,τ])
p次の変換の解としました。
 これらはfの取れ方によって
ω=4Kp,4mK+4iKp(m=0,1,,p1)
p+1通りに帰着できます。
 これと上で定めたp次の変換は一見異なるように見えますが、sn[u/M,Aτ]の零点の整理の仕方を変えることで以下のような対応があることがわかります。

  • j=pのとき零点は
    u=2nKp,2K2nKp(|n|p12)
    であったのでnが偶数のときn=2nmが奇数のときn=±p2nとすることで
    u=4nKp,2K4nKp(|n|p12)
    と整理できる。
  • jpのとき
    m={2j2jp12jp2jp+1
    とすると零点は
    u=2mK+2iKpn,2K2mK+2iKpn(|n|p12)
    であったのでnが偶数のときn=2nnが奇数のときn=±p2nとおくことで
    u=4mK+4iKpn,2K4mK+4iKpn(|n|p12)
    と整理できる。

おわりに

 はい。
 というわけで楕円関数の変換という問題について解析的に考えることで「楕円関数の変換の具体形」がわかりました。次の記事ではこれを利用し楕円関数の変換を代数的に考えることで「楕円関数の変換の具体値」を求めます。それによってp次の変換におけるklの関係やMの値が具体的に求まっていくのでお楽しみに。
 とりあえず今回はこんなところで。
 では。

参考文献

[1]
A. G. Greenhill, The applications of elliptic functions, Dover, New York, 1959, pp. 305-328
[2]
Jacobi, 高瀬正仁 訳, ヤコビ楕円関数原論, 講談社, 2012
投稿日:20221220
OptHub AI Competition

この記事を高評価した人

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

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

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

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

投稿者

子葉
子葉
1070
263431
主に複素解析、代数学、数論を学んでおります。 私の経験上、その証明が簡単に探しても見つからない、英語の文献を漁らないと載ってない、なんて定理の解説を主にやっていきます。 同じ経験をしている人の助けになれば。最近は自分用のノートになっている節があります。

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中
  1. はじめに
  2. 楕円関数の変換
  3. 問題の言い換え
  4. 補題その1
  5. 補題その2
  6. 素数次の変換
  7. 2次の変換
  8. p次の変換(p:奇素数)
  9. おわりに
  10. 参考文献