ヤコビの楕円関数

数学において、**ヤコビの楕円関数**（ヤコビのだえんかんすう、）とは、基本的な楕円関数の一群であり、追加でテータ関数を含むこともあり、歴史的に重要な関数からなる。これらの関数は重要な構造を持っていて、さらに直接関連した応用も存在する。三角関数との類似性も便利で、*sin* に対応する関数を *sn* と表記する。実用的な問題にはヴァイエルシュトラスの楕円函数よりもヤコビの楕円関数のほうがよく用いられる。これは複素解析の概念を使わずに定義し考察できるからである。これらの関数はにより導入された。

## 導入

ヤコビの楕円関数は全部で12種類ある。これら12種は長方形のある頂点から他の頂点へ引いた矢印に対応している。ここでは、この頂点を順に s、c、d、n と呼ぶことにする。この長方形を複素平面に配置して、s は原点に、c は実軸上の *K* の位置に、d は *K*&nbsp;+&nbsp;*iK* の位置に、n は虚軸上の *iK* の位置になるようにする。実数 *K* と *K' *はと呼ばれる。このとき、 s、c、d、n から異なる2文字を選んで「p」と「q」とすると、ヤコビの楕円関数は「pq」と書くことができる。

ヤコビの楕円関数は二重周期を持つ有理型関数で、次の性質を満たす唯一のものをいう。
* 頂点 p に一位の零点を持ち、頂点 q には一位の極を持つ
* p から q までが関数 pq&nbsp;*u* の半周期となる。つまり、関数 pq&nbsp;*u* は pq の向きの周期を持っており、その周期は p と q の距離の倍である。さらに、関数 pq&nbsp;*u* はほかの2方向についても周期的であり、p から残りの頂点への距離が 1/4 周期である。
* 関数 pq&nbsp;*u* を各頂点で *u* について展開すると、先頭の項の係数は 1 となる。言い換えると、関数 pq&nbsp;*u* を、頂点 p において展開した場合の先頭の項は *u* であり、頂点 q では 1/*u* であり、残りの頂点では 1 である。

より一般的には、長方形である必要はなく、平行四辺形でもよい。しかし、*K* と *iK' * をそれぞれ実軸と虚軸に置いておくと、ヤコビの楕円関数 pq&nbsp;*u* は、*u* が実数のとき実数値を取る。

## 記法
楕円関数には様々な記法があり、無用な混乱を引き起こしている。楕円関数は2変数の関数である。最初の変数は**振幅**φを使って表すこともあるが、一般には、以下のように *u* を使う。二番目の変数は**パラメタ** *m* を使ったり、***' ''k'' を使って表す。ここで、''k''<sup>2</sup> = ''m'' である。他にも、***' &alpha; を使うこともある。*m*&nbsp;= &nbsp;sin2&nbsp;&alpha; である。これら別記法の定義や発展した話題、*相補的*な対応概念については、楕円積分やの記事を見られたい。

## 楕円積分の逆関数による定義
 

上記のように、特定の性質を持つ唯一の有理型関数として定義するのは非常に抽象的である。より単純で、完全に同値な定義として、第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&nbsp;*u*&nbsp;=&nbsp;&Delta;(*u*) を**delta amplitude**と呼ぶこともある。上の値 *m* は自由なパラメタで、通常は実数を使い、0&nbsp;≤&nbsp;*m*&nbsp;≤&nbsp;1 とする。これにより、楕円関数は、ふたつの変数、振幅\phiとパラメタ*m*の関数だと考えることができる。

残りの9つの楕円関数は上記の3つから簡単に構成することができ、下の節で与えられる。

注意として、\phi=\pi/2 の場合は、*u* は四半周期 *K* となる。

## 三角法による定義

*r*&nbsp;=&nbsp;1 の単位円上で \cos \theta, \sin \thetaを定義できたのと同様に、ヤコビの楕円関数は、*a*&nbsp;=&nbsp;1 の単位楕円上で定義できる。

:
\begin{align}
& x^2 + \frac{y^2}{b^2} = 1, \quad  b >  1, \\
& m = 1 - \frac{1}{b^2}, \quad 0 

とすると、

: 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*(&tau;)の式で定義されたので、これを反転して τ を *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) &pi;) ~ 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 と解釈する。

（この記法は、とグレイシャーによるもので、ヤコビの元々の記法にはない）

## 加法定理
ヤコビの楕円関数が持つ代数的な関係式として

:\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*&nbsp;+&nbsp;*m*1&nbsp;=&nbsp;1 and *m*&nbsp;=&nbsp;*k*2である。

pq2 &middot; 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&nbsp;\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)である。これらの逆関数は楕円積分で表すことができる。また、冪級数でも表現できる。

## 地図投影法
パース・クインカンシャル図法はヤコビの楕円関数を用いた投影法である。また各国で標準的に使われているガウス・クリューゲル図法でも、投影式をヤコビの楕円関数で表すことができる。

## 脚注

### 注釈

### 出典

## 参考文献
* 
* 
* N. I. Akhiezer, *Elements of the Theory of Elliptic Functions*, (1970) Moscow, translated into English as *AMS Translations of Mathematical Monographs Volume 79* (1990) AMS, Rhode Island ISBN 0-8218-4532-2
* A. C. Dixon  The elementary properties of the elliptic functions, with examples (Macmillan, 1894)
* A. G. Greenhill The applications of elliptic functions (London, New York, Macmillan, 1892)
* H. Hancock Lectures on the theory of elliptic functions (New York, J. Wiley & sons, 1910)
*
*
* E. T. Whittaker and G. N. Watson **, (1940, 1996) Cambridge University Press. ISBN 0-521-58807-3
*  P. Appell and E. Lacour Principes de la théorie des fonctions elliptiques et applications (Paris, Gauthier Villars, 1897)
*  Traité des fonctions elliptiques et de leurs applications (vol. 1) (Paris, Gauthier-Villars, 1886–1891)
*  G. H. Halphen Traité des fonctions elliptiques et de leurs applications (vol. 2) (Paris, Gauthier-Villars, 1886–1891)
*  G. H. Halphen Traité des fonctions elliptiques et de leurs applications (vol. 3) (Paris, Gauthier-Villars, 1886–1891)
*   and J. Molk 	Eléments de la théorie des fonctions elliptiques. Tome I, Introduction. Calcul différentiel. Ire partie (Paris : Gauthier-Villars et fils, 1893)
*  J. Tannery and J. Molk  	Eléments de la théorie des fonctions elliptiques. Tome II, Calcul différentiel. IIe partie (Paris : Gauthier-Villars et fils, 1893)
*  J. Tannery and J. Molk  	Eléments de la théorie des fonctions elliptiques. Tome III, Calcul intégral. Ire partie, Théorèmes généraux. Inversion (Paris : Gauthier-Villars et fils, 1893)
*  J. Tannery and J. Molk  	Eléments de la théorie des fonctions elliptiques. Tome IV, Calcul intégral. IIe partie, Applications (Paris : Gauthier-Villars et fils, 1893)
*  C. Briot and J. C. Bouquet Théorie des fonctions elliptiques (  Paris : Gauthier-Villars, 1875)

## 関連項目
* 楕円積分
* 楕円曲線
* 
* 
* ヴァイエルシュトラスの楕円函数
* ヤコビのテータ関数
* 

## 外部リンク
* 
* 
* ヤコビ楕円関数 - サンプルソースコード C/C++
* ヤコビ楕円関数 - サンプルソースコード Fortran