楕円積分と楕円関数
楕円積分や楕円関数を使うときは,特殊関数の場合以上に引数の指定には注意を要する.これらの関数の記述や使い方に関しては,相違した規約がしばしば使われる.その多くは,引数の呼び方において違いがあったり,引数と引数の間に入れる区切り記号としてコンマ以外の記号が使われたりする.
• 振幅 (Mathematica で用いられる,ラジアンで指定) |
• 引数 u(Mathematica で用いられる): =am (u) で振幅に関連する |
• デルタ振幅 ( ): |
• 座標 x:x=sin ( ) |
| • 特性 n(第三種楕円積分において Mathematica で用いられる) |
• パラメータ m(Mathematica で用いられる):I ( m)のように垂直バー ( )に続く |
| • 補パラメータ m1:m1=1-m |
• 母数 (Modulus) k:I ( , k); m=k2 のようにコンマに続く |
• モジュラ角 :I ( \ ); m=sin2 ( )のようにバックスラッシュ (\ )に続く |
• ノーム q: 関数でコンマに続く; q=exp[- K (1-m)/K{m)]=exp (i  / ) |
| • 不変数 g2,g3(Mathematica で用いられる) |
• 半周期 , :g2=60 w-4,g3=140 w-6,ここで w=2r +2s  |
• 周期の比 : = / |
• 判別式 : |
| • 曲線のパラメータ a,b(Mathematica で用いられる) |
| • 座標 y(Mathematica で用いられる):y2=x3+ax2+b x と関連している |
楕円積分と楕円関数の共通的な引数に関する規約
異なった引数規約の間における変換
楕円積分
楕円積分
R を有理関数,
y2を
x の三次または四次の多項式とするとき,
R (x, y)
x で表される積分を楕円積分と呼ぶ.あらゆる楕円積分は,3種類あるルジャンドル・ヤコビの楕円積分で表すことができる.
第一種楕円積分
EllipticF[
, m]は,
-
/2<
<
/2において


で与えられる.この楕円積分は,単振子の運動状態を表した方程式の解法に使われる.また,第一種の不完全楕円積分としても知られる.
文献によっては,
Mathematica で使われる引数の順序とは逆の順序で引数が与えられることがあるので注意する.
第一種完全楕円積分
EllipticK[m]は,

で与えられる.ここで,
K は第一種完全楕円積分を表し,
F は不完全楕円積分を表すことに注目されたい.実際の用途では,
m が省略れ,
K (m) は単に
K と記述されることがよくある.補第一種完全楕円積分
K
(m)は,
K (1-m)で与えられる.この積分はしばしば
K
とも表記される.
K と
iK
は,対応するヤコビの楕円積分(
「楕円関数」を参照)の「実部」と「虚部」の各々からなる4分の1周期を与える.
第二種楕円積分
EllipticE[
, m] は,
-
/2<
<
/2の範囲において,


で与えられる.
第二種完全楕円積分
EllipticE[m]は,

で与えられ,
E とも表記される.この補形は
E
(m)=E (1-m)となる.
ヒューマン(Heuman)のラムダ関数は,

で与えられる.
第三種楕円積分
EllipticPi[n,
, m]は,

で与えられる.
第三種完全楕円積分
EllipticPi[n, m]は,

で与えられる.
| Out[1]= |  |
|
| Out[2]= |  |
|
プロットすると,楕円積分が複素平面上で複雑な構造を持つことが分かる.
| Out[1]= |  |
|
楕円関数
| JacobiAmplitude[u,m] | 振幅関数 am (u m) |
| JacobiSN[u,m], JacobiCN[u,m], etc. |
| ヤコビの楕円関数 sn (u m)等 |
| InverseJacobiSN[v,m], InverseJacobiCN[v,m], etc. |
| 逆ヤコビの楕円関数 sn-1 (v m)等 |
| EllipticTheta[a,u,q] | 楕円シータ関数 a (u, q) (a=1, ..., 4) |
| EllipticThetaPrime[a,u,q] | 楕円シータ導関数 (a=1, ..., 4) |
SiegelTheta[ ,s] | シーゲルのシータ関数 ( , s) |
SiegelTheta[v, ,s] | シーゲルのシータ関数 [v] ( , s) |
| WeierstrassP[u,{g2,g3}] | ワイエルシュトラスの楕円関数 (u;g2, g3) |
| WeierstrassPPrime[u,{g2,g3}] |
| ワイエルシュトラスの楕円導関数  (u;g2, g3) |
| InverseWeierstrassP[p,{g2,g3}] |
| 逆ワイエルシュトラスの楕円関数 |
| WeierstrassSigma[u,{g2,g3}] | ワイエルシュトラスのシグマ関数 (u;g2, g3) |
| WeierstrassZeta[u,{g2,g3}] | ワイエルシュトラスのゼータ関数 (u;g2, g3) |
楕円関数とこれに関連した関数
二次形式の平方根を含む有理関数は,逆三角関数を使って積分することができる.このため,三角関数をこれらの積分で得られる関数の逆として定義することもできる.
類推から,楕円関数は楕円積分で得られる関数の逆として定義することができる.
ヤコビの楕円関数の振幅
JacobiAmplitude[u, m]は,第一種楕円積分の逆関数である.
u=F (
m)ならば,
=am (u
m)と書ける.ヤコビの楕円関数を用いるとき,しばしば
m が省略され,
am (u
m)が
am (u)と書かれる.
ヤコビの楕円関数
JacobiSN[u, m]および
JacobiCN[u, m]は,それぞれ
sn (u)=sin (
)および
cn (u)=cos (
)で与えられる.ただし,
=am (u
m)とする.また,
JacobiDN[u, m]は,

で与えられる.
ヤコビの楕円関数
JacobiPQ[u, m]は全部で12種類存在する.それらは,
P と
Q を
S,
C,
D,
Nから組み合せることで作れる.それぞれのヤコビの楕円関数
JacobiPQ[u, m]は,
pq (u)=pn (u)/qn (u)を満足する.ただし,この式では
nn (u)=1とする.
三角関数の間にはさまざまな関係式があるが,それに似た関係がヤコビの楕円関数の間にも数多くある.ヤコビの楕円関数は,極限を取ったとき,三角関数に還元する.このため,
sn (u
0)=sin (u),
sn (u
1)=tanh (u),
cn (u
0)=cos (u),
cn (u
1)=sech (u),
dn (u
0)=1,
dn (u
1)=sech (u)というような関係が成り立つ.
Pq (u)の記述が積分

を表すのに使われることがある.
「楕円積分」で定義したヤコビのゼータ関数を使って表せる.
楕円関数の重要な特性のひとつとして,引数の複素数値について二重の周期性を持つことが挙げられる.通常の三角関数は単一な周期性を持つ.つまり,任意の整数
s について
f (z+s
)=f (z) の周期性を持つ.これに対して,楕円関数は二重周期性を持つ.つまり,
f (z+r
+s
)=f (z)の周期性が
r と
s の任意の整数の対について成り立つ.
sn (u
m)のようなヤコビの楕円関数は,複素
u 平面上で2重の周期性を持つ.周期には
=4K (m)と

=4iK (1-m)が含まれる.ここで,
K は第一種完全楕円積分を表す.
ヤコビの楕円関数の記述
pq (u
m)においてpとqのどちらを選択するかは,4分の1周期の点で得られる
Kと
iK
の関数値から判明する.
ヤコビの楕円関数  の絶対値が2つの方向に沿って基本周期を持つのがグラフから分かる.
| Out[3]= |  |
|
Mathematica には,逆ヤコビの楕円関数である
InverseJacobiSN[v, m]や
InverseJacobiCN[v, m]等も組み込まれている.例えば,逆関数の1つ,
sn-1 (v
m)は,
v=sn (u
m)が成り立つための
u の値を与える.逆ヤコビの楕円関数は楕円積分に関連している.
4種類の楕円シータ関数
a (u, q)は,
aを
1,
2,
3または
4として
EllipticTheta[a, u, q]から作ることができる.これらの関数は,
1 (u, q)=2q1/4
(-1)nqn (n+1)sin[ (2n+1)u],
2 (u, q)=2q1/4
qn (n+1)cos[ (2n+1)u],
3 (u, q)=1+2
qn2cos (2nu),
4 (u, q)=1+2
(-1)nqn2cos (2nu)で定義される.シータ関数はパラメータ
q を明示しないで
a (u)と書かれることもある.シータ関数は,
(u
m)と書かれることもあるが,その場合に
m と
q は
q=exp[-
K (1-m)/K (m)]という式で関係付けられている.また,
q は
q=ei

で与えられる

に置き換えられることもある.これらすべてのシータ関数は拡散方程式
2
(u,
)/
u2=4
i 
(u,
)/

を満たす.
次元
p,ベクトル
s のリーマンの正方モジュラ行列

を持つシーゲルのシータ関数
SiegelTheta[
, s]は,複素次元
p までの楕円シータ関数を生成する.これは

で定義され,
n はすべての
p 次元整数ベクトル上を動く.特性関数
SiegelTheta[
,
, s]を持つシーゲルのシータ関数は,

で定義され,特性

は1対の
p 次元ベクトル
{
,
}を表す.
ヤコビの楕円関数はシータ関数の比としても表すことができる.
ワイエルシュトラス(Weierstrass)の楕円関数
WeierstrassP[u, {g2, g3}] は,楕円積分の逆関数と見ることができる.ワイエルシュトラスの関数
(u;g2, g3)は,

を満たす
x の値を与える.関数
WeierstrassPPrime[u, {g2, g3}]は,

で与えられる.
ワイエルシュトラスの関数は,基本半周期

と


を用いて書かれることもあり,この

と


は
WeierstrassHalfPeriods[{u, {g2, g3}]を使い不変数
g2 と
g3 から求まる.
逆ワイエルシュトラスの楕円関数
InverseWeierstrassP[p, {g2, g3}]は,
p=
(u;g2, g3)とする
u の2つの値の1つを返す.この値は複素半周期

と


によって定義される平行四辺形内にある.
逆ワイエルシュトラスの楕円関数
InverseWeierstrassP[{p, q}, {g2, g3}]は,
p=
(u;g2, g3)および
q=
(u;g2, g3)となる
u の一意的な値を返す.
u にそのような値が存在するには,
p と
q は
q2=4p3-g2p-g3の関係式を満たす必要がある.
ワイエルシュトラスのゼータ関数
WeierstrassZeta[u, {g2, g3}] とワイエルシュトラスのシグマ関数
WeierstrassSigma[u, {g2, g3}]は,ワイエルシュトラスの楕円関数に各々,

(z;g2, g3)=-
(z;g2, g3)と

(z;g2, g3)/
(z;g2, g3)=
(z;g2, g3)の関係式で関連している.
ワイエルシュトラスのゼータ関数とワイエルシュトラスのシグマ関数は周期性を持たないため,厳密には楕円関数ではない.
楕円モジュラ関数
楕円モジュラ関数
クライン(Klein)の不変モジュラ関数
KleinInvariantJ[
]とデデキント(Dedekind)のイータ関数
DedekindEta[
]は,

の関係式を満たす.
一般化された楕円積分と楕円関数
一般化された楕円積分と楕円関数
ここまで述べた楕円積分や関数の定義は伝統的な使い方に基づいている.現代の代数幾何学では,もう少し一般化された定義の方が都合がよい.
関数
EllipticLog[{x, y}, {a, b}]は,積分

の値として定義される.ただし,平方根の符号は

とする
y を与えることで指定する.

の形の積分は,常用対数(および,逆三角関数)を用いて表せる.平方根でくくられた多項式を三次にして,これを一般化したのが
EllipticLogと考えるとよい.
EllipticExp[u, {a, b}]は
EllipticLogの逆関数である.この関数は
EllipticLogに引数としてある
{x, y}のリストを返す.
EllipticExpは楕円関数で,複素
u 平面上で2重周期を持つ.
ArithmeticGeometricMean[a, b]は,2個の数
a と
b の算術幾何平均(AGM)を返す.この量は,楕円積分や他の関数を計算する数値アルゴリズムにおいて中心的な役割を果たす.
a0=a と
b0=b を初期値として,必要な精度の下で
an=bnとなるまで変換

,

を繰り返すことで,AGMが求まる.