「ヤコビの楕円関数」cd(x, k) を用いて、
以下のように定義される関数をチェビシェフ有理関数(elliptic rational function)という。
R
n, k1(x
)
=
cd(cd-1(
x,
k
)
,
k
1)
k1' = √1 - k12
K1 = K(k1),
K1' = K(k1')
K = K(k),
K' = K(k')
k は、
= n を満たすように選ぶ(
k' = √1 - k2)。
ただし、
K(k) は第1種完全楕円積分である。
数値計算上、k は
k
=
q-1(
exp
(
-
))
で求める。
(q(k) は Jacobi のテータ関数のノーム。)
ヤコビの楕円関数の詳細は
「ヤコビの楕円関数」
参照。
関数のグラフ。
↑の導出の過程。
零点と極 → 有理式で表す。
複素平面状に以下のような経路
C = C1 + C2 + C3
を考える。
経路 C 上での cd(u, k) の値の変化の仕方を表1に、値を3次元的にグラフ化したものを図1に示す。
図1中では、
経路 C をオレンジ色、
C1 上での値を赤色、
C2 上での値を青色、
C3 上での値を緑色の線で表す。
cn(u, k)の値の変化の仕方
経路 |
式変形 |
値の変化の仕方 |
C1(K→0)
|
cd(K - u, k)
=
sn(u, k)
|
0→1の間で単調増加。
|
C2(0→i K')
|
cd(i u, k)
=
|
1→の間で単調増加。
|
C3(i K'→K + i K')
|
cd(u + i K', k)
=
|
→∞の間で単調増加。
|
cn(u, k) の C 上での値
cd(u, k) は C 上で単調増加 →
cd の逆関数
cd-1(x, k) の実軸上での値は経路 C に。
次に、経路 C に定数
を掛けた経路
D = D1 + D2 + D3
を考える。
-
D1 =
C1… 実軸上をn K1→0と進む。
-
D2 =
C2… 実軸上を0→i K1'と進む。 ( = K1')
-
D3 =
C3…Im(u) = i K1'上をi K1'→n K1 + i K1'と進む。
この経路 D 上での cd(u, k1) の値を3次元的にグラフ化したものを図2に示す。
図2中では、
経路 C を紫色、
C1 上での値を赤色、
C2 上での値を青色、
C3 上での値を緑色の線で表す。
cn(u, k1) の D 上での値
「
実軸 [0, ∞] →
(cd-1(x, k))→
経路 C →
(
を掛ける)→
経路
D →
(
cd(u, k1))→
Rn, k1(x)」
となるので、チェビシェフ有理関数
Rn, k1(x) のグラフは図3のようになる。
チェビシェフ有理関数のグラフ
Rn, k1(x) は n 個の零点と極で表される有理式となる。
まず、Rn, k1(x) の零点/極を求める。
cd(u, k) の零点は
u = (2m + 1)K(k)
なので、
Rn, k1(x) の零点は、
cd-1(
x,
k
)
=
(2m + 1
)K
1
したがって、
x
=
cd(K,
k
)
=
sn(K,
k
)
同様に、
cd(u, k) の極は
u = (2m + 1)K(k) + i K(k')
なので、
Rn, k1(x) の極は、
cd-1(
x,
k
)
=
(2m + 1
)K
1 + i K
1'
したがって、
と置くと、
sn の性質から、
ωm, n = -ωn - m - 1, n となる。
Rn, k1(x)は、
Rn, k1(x)
=
x
x2 - ωm, n2 |
k2 ωm, n2 x2 - 1 |
… n が奇数のとき Rn, k1(x)
=
x2 - ωm, n2 |
k2 ωm, n2 x2 - 1 |
… n が偶数のとき