ヤコビの楕円関数
数学において、ヤコビの楕円関数(ヤコビのだえんかんすう、Jacobi elliptic functions)とは、基本的な楕円関数の一群であり、追加でテータ関数を含むこともあり、歴史的に重要な関数からなる。これらの関数は重要な構造を持っていて、さらに直接関連した応用も存在する。三角関数との類似性も便利で、sin に対応する関数を sn と表記する。実用的な問題にはヴァイエルシュトラスの楕円函数よりもヤコビの楕円関数のほうがよく用いられる。これは複素解析の概念を使わずに定義し考察できるからである。これらの関数はにより導入された。
ヤコビの楕円関数は全部で12種類ある。これら12種は長方形のある頂点から他の頂点へ引いた矢印に対応している。ここでは、この頂点を順に s、c、d、n と呼ぶことにする。この長方形を複素平面に配置して、s は原点に、c は実軸上の K の位置に、d は K + iK の位置に、n は虚軸上の iK の位置になるようにする。実数 K と K' は四半周期と呼ばれる。このとき、 s、c、d、n から異なる2文字を選んで「p」と「q」とすると、ヤコビの楕円関数は「pq」と書くことができる。
ヤコビの楕円関数は二重周期を持つ有理型関数で、次の性質を満たす唯一のものをいう。
頂点 p に一位の零点を持ち、頂点 q には一位の極を持つ
p から q までが関数 pq u の半周期となる。つまり、関数 pq u は pq の向きの周期を持っており、その周期は p と q の距離の倍である。さらに、関数 pq u はほかの2方向についても周期的であり、p から残りの頂点への距離が 1/4 周期である。
関数 pq u を各頂点で u について展開すると、先頭の項の係数は 1 となる。言い換えると、関数 pq u を、頂点 p において展開した場合の先頭の項は u であり、頂点 q では 1/u であり、残りの頂点では 1 である。
より一般的には、長方形である必要はなく、平行四辺形でもよい。しかし、K と iK'  をそれぞれ実軸と虚軸に置いておくと、ヤコビの楕円関数 pq u は、u が実数のとき実数値を取る。
楕円関数には様々な記法があり、無用な混乱を引き起こしている。楕円関数は2変数の関数である。最初の変数は振幅φを使って表すこともあるが、一般には、以下のように u を使う。二番目の変数はパラメタ m を使ったり、母数 (楕円関数) k を使って表す。ここで、k2 = m である。他にも、modular angle α を使うこともある。m =  sin2 α である。これら別記法の定義や発展した話題、相補的な対応概念については、楕円積分や四半周期の記事を見られたい。
上記のように、特定の性質を持つ唯一の有理型関数として定義するのは非常に抽象的である。より単純で、完全に同値な定義として、第1種不完全楕円積分の逆関数として定義することができる。まず、
u=\int_0^\phi \frac{\mathrm d\theta} {\sqrt {1-m \sin^2 \theta}
とする。楕円関数 sn u は
\operatorname {sn}\; u = \sin \phi\,
で与えられ、cn u は
\operatorname {cn}\; u = \cos \phi
で与えられ、dn u は
\operatorname {dn}\; u = \sqrt {1-m\sin^2 \phi}
で与えられる。
ここで、角 \phi を振幅と呼ぶ。dn u = Δ(u) をdelta amplitudeと呼ぶこともある。上の値 m は自由なパラメタで、通常は実数を使い、0 ≤ m ≤ 1 とする。これにより、楕円関数は、ふたつの変数、振幅\phiとパラメタmの関数だと考えることができる。
残りの9つの楕円関数は上記の3つから簡単に構成することができ、下の節で与えられる。
注意として、\phi=\pi/2 の場合は、u は四半周期 K となる。
r = 1 の単位円上で \cos \theta, \sin \thetaを定義できたのと同様に、ヤコビの楕円関数は、a = 1 の単位楕円上で定義できる。
\begin{align}
& x^2 + \frac{y^2}{b^2} = 1, \quad  b >  1, \\
& m = 1 - \frac{1}{b^2}, \quad 0 < m < 1, \\
& x = r \cos \theta, \quad y = r \sin \theta
\end{align}
とすると、
r( \theta,m) =  \frac{1} {\sqrt {1-m \sin^2 \theta}
である。
楕円に沿った弧長 u を計算すると、
u = u(\varphi,m)=\int_0^\varphi r(\theta,m) \, d\theta
となる。ここで、 \varphi = arccos(x) である。単位円でおなじみの関係式
\cos \theta = x , \quad \sin \theta = y
を楕円の場合に一般化して、
\operatorname{cn}(u,m)=x,\quad \operatorname{sn}(u,m) = \frac{y}{b}
を定義とする。
さらに、
\operatorname{dn}(u,m) = \frac 1 {r(\varphi,m)}
と定める。
ヤコビによるテータ関数を用いて定義することもできる。\vartheta(0;\tau)を\varthetaと略記し、\vartheta_{01}(0;\tau), \vartheta_{10}(0;\tau), \vartheta_{11}(0;\tau)もそれぞれ、\vartheta_{01}, \vartheta_{10}, \vartheta_{11}(テータ定数と呼ばれる)と略記する。このとき、母数 k はk=\left({\vartheta_{10} \over \vartheta}\right)^2 となる。u = \pi \vartheta^2 zとおくと、
\mathrm{sn}(u; k) = -{\vartheta \vartheta_{11}(z;\tau) \over \vartheta_{10} \vartheta_{01}(z;\tau)}
\mathrm{cn}(u; k) = {\vartheta_{01} \vartheta_{10}(z;\tau) \over \vartheta_{10} \vartheta_{01}(z;\tau)}
\mathrm{dn}(u; k) = {\vartheta_{01} \vartheta(z;\tau) \over \vartheta \vartheta_{01}(z;\tau)}
となる。
これで、ヤコビの関数が母数k(τ)の式で定義されたので、これを反転して τ を k の式で表す必要がある。初めに、余母数k' = \sqrt{1-k^2}を考える。これは τ の関数として、
k'(\tau) = \left({\vartheta_{01} \over \vartheta}\right)^2
と書ける。
次に、
\ell = {1 \over 2} {1-\sqrt{k'} \over 1+\sqrt{k'} =
{1 \over 2} {\vartheta - \vartheta_{01} \over \vartheta + \vartheta_{01}
と定める。
そして、ノーム q を q = \exp (\pi i \tau) と定義し、ノーム q に関して、\ellを冪級数の商に展開すると、
\ell = {q + q^9 + q^{25} + \cdots \over 1 + 2q^4 + 2q^{16} + \cdots}
となる。
級数の反転を行うと、
q = \ell + 2\ell^5 + 15\ell^9 + 150\ell^{13} + 1707\ell^{17} + 20910\ell^{21} + 268616\ell^{25} + \cdots
を得る。
τの虚部が 1/2 sqrt(3) 以上の場合に帰着すればよいので、q の絶対値は exp(-1/2 sqrt(3) π) ~ 0.0658 以下である場合だけを考えればよい。この値は小さいため、上の級数は急速に収束し、簡単にqの値を計算することができる。
関数名の文字の順序を反転した場合は、元の関数の逆数となる。
\begin{align}
\operatorname{ns}(u) & = \frac{1}{\operatorname{sn}(u)} \\[8pt]
\operatorname{nc}(u) & = \frac{1}{\operatorname{cn}(u)} \\[8pt]
\operatorname{nd}(u) & = \frac{1}{\operatorname{dn}(u)}
\end{align}
同様に、上の3つの関数どうしの比は、分子の関数の最初の文字と分母の関数の最初の文字を繋げたものになる。
\begin{align}
\operatorname{sc}(u) & = \frac{\operatorname{sn}(u)}{\operatorname{cn}(u)} \\[8pt]
\operatorname{sd}(u) & = \frac{\operatorname{sn}(u)}{\operatorname{dn}(u)} \\[8pt]
\operatorname{dc}(u) & = \frac{\operatorname{dn}(u)}{\operatorname{cn}(u)} \\[8pt]
\operatorname{ds}(u) & = \frac{\operatorname{dn}(u)}{\operatorname{sn}(u)} \\[8pt]
\operatorname{cs}(u) & = \frac{\operatorname{cn}(u)}{\operatorname{sn}(u)} \\[8pt]
\operatorname{cd}(u) & = \frac{\operatorname{cn}(u)}{\operatorname{dn}(u)}
\end{align}
もっと簡単に、
\operatorname{pq}(u)=\frac{\operatorname{pr}(u)}{\operatorname{qr}(u)}
と書くことができる。ここで、p、q、r は、s、c、d、n の任意の文字で、ss = cc = dd = nn = 1 と解釈する。
(この記法は、クリストフ・グーデルマンとJames Whitbread Lee Glaisherによるもので、ヤコビの元々の記法にはない)
\operatorname{cn}^2(u,k) + \operatorname{sn}^2(u,k) = 1
\operatorname{dn}^2(u,k) + k^2 \ \operatorname{sn}^2(u,k) = 1
がある。
これらの方程式で定まる2つの二次曲面の共通部分は楕円曲線であり、(cn, sn, dn)は楕円曲線のパラメタ表示を与えることが分かる。ヤコビの楕円関数の加法定理により、この楕円曲線の点は群となる。
\begin{align}
\operatorname{cn}(x+y) & =
{\operatorname{cn}(x)\;\operatorname{cn}(y)
- \operatorname{sn}(x)\;\operatorname{sn}(y)\;\operatorname{dn}(x)\;\operatorname{dn}(y)
\over {1 - k^2 \;\operatorname{sn}^2 (x) \;\operatorname{sn}^2 (y)}, \\[8pt]
\operatorname{sn}(x+y) & =
{\operatorname{sn}(x)\;\operatorname{cn}(y)\;\operatorname{dn}(y) +
\operatorname{sn}(y)\;\operatorname{cn}(x)\;\operatorname{dn}(x)
\over {1 - k^2 \;\operatorname{sn}^2 (x)\; \operatorname{sn}^2 (y)}, \\[8pt]
\operatorname{dn}(x+y) & =
{\operatorname{dn}(x)\;\operatorname{dn}(y)
- k^2 \;\operatorname{sn}(x)\;\operatorname{sn}(y)\;\operatorname{cn}(x)\;\operatorname{cn}(y)
\over {1 - k^2 \;\operatorname{sn}^2 (x)\; \operatorname{sn}^2 (y)}.
\end{align}
-\operatorname{dn}^2(u)+m_1= -m\;\operatorname{cn}^2(u) = m\;\operatorname{sn}^2(u)-m
-m_1\;\operatorname{nd}^2(u)+m_1= -mm_1\;\operatorname{sd}^2(u) = m\;\operatorname{cd}^2(u)-m
m_1\;\operatorname{sc}^2(u)+m_1= m_1\;\operatorname{nc}^2(u) = \operatorname{dc}^2(u)-m
\operatorname{cs}^2(u)+m_1=\operatorname{ds}^2(u)=\operatorname{ns}^2(u)-m
ここで、m + m1 = 1 and m = k2である。
pq2 · qp2 = 1 と pq = pr / qr を使うことで、二乗についてのさらなる関係式を得られる。ここで、p、q、rは s、c、d、n の任意の文字であり、ss = cc = dd = nn = 1 とする。
ノームをq=\exp(-\pi K'/K)とし、引数をv=\pi u /(2K)と変換する。このときの関数の展開はランベルト級数になる。
\operatorname{sn}(u)=\frac{2\pi}{K\sqrt{m}
\sum_{n=0}^\infty \frac{q^{n+1/2}{1-q^{2n+1} \sin ((2n+1)v),
\operatorname{cn}(u)=\frac{2\pi}{K\sqrt{m}
\sum_{n=0}^\infty \frac{q^{n+1/2}{1+q^{2n+1} \cos ((2n+1)v),
\operatorname{dn}(u)=\frac{\pi}{2K} + \frac{2\pi}{K}
\sum_{n=1}^\infty \frac{q^{n}{1+q^{2n} \cos (2nv).
ヤコビの楕円関数を微分すると、
\frac{\mathrm{d}{\mathrm{d}z}\, \mathrm{sn}\,(z) = \mathrm{cn}\,(z)\, \mathrm{dn}\,(z)
\frac{\mathrm{d}{\mathrm{d}z}\, \mathrm{cn}\,(z) = -\mathrm{sn}\,(z)\, \mathrm{dn}\,(z)
\frac{\mathrm{d}{\mathrm{d}z}\, \mathrm{dn}\,(z) = - k^2 \mathrm{sn}\,(z)\, \mathrm{cn}\,(z)
となる。
上記の加法定理を使うと、k が 0 < k < 1 を満たす場合は、ヤコビの楕円関数は下記の非線形常微分方程式の解となる。
\mathrm{sn}\,(x) は次の微分方程式の解である
\frac{\mathrm{d}^2 y}{\mathrm{d}x^2} + (1+k^2) y - 2 k^2 y^3 = 0
\left(\frac{\mathrm{d} y}{\mathrm{d}x}\right)^2 = (1-y^2) (1-k^2 y^2)
\mathrm{cn}\,(x) は次の微分方程式の解である
\frac{\mathrm{d}^2 y}{\mathrm{d}x^2} + (1-2k^2) y + 2 k^2 y^3 = 0
\left(\frac{\mathrm{d} y}{\mathrm{d}x}\right)^2 = (1-y^2) (1-k^2 + k^2 y^2)
\mathrm{dn}\,(x) は次の微分方程式の解である
\frac{\mathrm{d}^2 y}{\mathrm{d}x^2} - (2 - k^2) y + 2 y^3 = 0
\left(\frac{\mathrm{d} y}{\mathrm{d}x}\right)^2 = (y^2 - 1) (1 - k^2 - y^2)
ヤコビの楕円関数の逆関数は逆三角関数と同様のやり方で定義される。x=\mathrm{sn}(\xi, k)に対して、\xi=\mathrm{arcsn}(x, k)である。これらの逆関数は楕円積分で表すことができる。また、冪級数でも表現できる。
パース・クインカンシャル図法はヤコビの楕円関数を用いた投影法である。また各国で標準的に使われているガウス・クリューゲル図法でも、投影式をヤコビの楕円関数で表すことができる。