\def\bm#1{{\boldsymbol{#1}}} \def\coloneqq{{:=}}
楕円関数
前々回に議論したように第一種楕円積分\displaystyle u=\int^z\frac{dz}{\sqrt{\varphi\left(z\right)}}の逆関数をz=f\left(u\right)で表せば、
f\left(u+mP\left(A\right)+nP\left(B\right)\right)=f\left(u\right)\hspace{1cm}\left(m,nは整数\right)
の関係がある事が知られています。しかしこれを用いてただちにf\left(u\right)が二重周期関数であるとは断定出来ません。なぜならP\left(A\right)とP\left(B\right)がある同一の数の倍数であってf\left(u\right)は実は単一周期関数であるのかも知れないからです。
今、仮にf\left(u\right)が単一周期関数であるとして、その基本周期を\omegaとすれば、R面上の一点zに対するuのすべての値は
\int_c^z\frac{dz}{\sqrt{\varphi\left(z\right)}}=u+n\omega\hspace{1cm}\left(nは整数\right)
で表されます。よって今
\varPhi\left(z\right)=e^{2\pi i\frac{u}{\omega}}
というzの関数を考えれば、これはR面上において一価です。なぜならば今zがR面上において一定ならば積分路が変わってもuがu+n\omegaに変わるに過ぎず、\displaystyle2\pi i\frac{u}{\omega}が2n\pi iだけ変化するに過ぎないから、\varPhi\left(z\right)は不変となります。そうしてuは第一種積分であるからzのいかなる値においても正則です。それならば\varPhi\left(z\right)はR面上で一価でいたる所正則であるから、定数でなければなりません(第5回の最後の定理)。したがってuも定数ということになります。これは今我々の考えている主張と反する結果です。
もしまたP\left(A\right),P\left(B\right)がともに0となってf\left(u\right)が無周期となれば、積分u自身がR面上で一価正則となります。したがって定数となってやはり不都合です。
これによってf\left(u\right)は確かに二重周期関数であると分かります。
注意: P\left(A\right)、P\left(B\right)のいずれか一方が0となったり、または比P\left(A\right)\colon P\left(B\right)が実数になったりすることはありません。それは上の証明から明らかです。
今数平面上に任意の一点uをとり、これとu_0+P\left(A\right)及びu_0+P\left(B\right)をそれぞれ結ぶ線分を二辺とする平行四辺形を作り、これと合同な平行四辺形をその上下左右に隙間なく並べて全平面を覆ったと考えれば、f\left(u\right)は各平行四辺形内で全く同様な値の分布を示すはずです。このような平行四辺形の一つずつを周期平行四辺形といいます。特にu_0=0とし、0,P\left(A\right),P\left(B\right),P\left(A\right)+P\left(B\right)を頂点とする周期平行四辺形を基本周期平行四辺形といいます。
一般に、楕円積分の逆という考えを離れて、一価解析関数で二重周期をもち、かつ無限遠点以外において有理型のもの(極以外の特異点をもたないもの)を楕円関数といいます。
我々の考えたz=f\left(u\right)は一つの楕円関数です。
sn 関数
u=\int_0^z\frac{dz}{\sqrt{\left(1-z^2\right)\left(1-k^2z^2\right)}}\tag{1}
をとり、その逆関数を
z=\mathrm{sn}\ u
と名付けます。ただし被積分関数の根号の値はz=0のときに+1となる方の分枝をとるものと定めます。
三角関数の記号\sinは\mathrm{sine}の略字ですが、\mathrm{sn}は何の略字でもなくただの関数記号です。したがって読むときには文字通り「エスエヌ」と読むことにします。
この関数の性質は被積分関数に含まれるkの数に関係し、kを\mathrm{sn}関数の母数(Modulus)といいます。母数kを明示する必要があるときには
z=\mathrm{sn}\left(u,k\right)
と書きます。
1-k^2=k^{‘2}とおくと、このk’というものも後にしばしば使用される重要な数となります。これを補母数(complementary modulus)といいます。
特にk=0のときは、(1)は
u=\int_0^z\frac{dz}{\sqrt{1-z^2}}=\sin^{-1}z
となります。ゆえに
\mathrm{sn}\left(u,0\right)=\sin u
またk=1のときは、(1)は
u=\int_0^z\frac{dz}{1-z^2}=\frac{1}{2}\log\frac{1+z}{1-z}
となります。ゆえに
\mathrm{sn}\left(u,1\right)=\frac{e^{2u}-1}{e^{2u}+1}=\tanh u
さて\mathrm{sn}\ uの性質としてまず第一に挙げるべきものは二重周期性です。\mathrm{sn}\ uが楕円関数であって、P\left(A\right)及びP\left(B\right)の整数倍の和が周期となることは今までの議論で明らかになりましたが、このP\left(A\right),P\left(B\right)の値についてさらに詳しく考えてみましょう。
第6回で証明したように、
P\left(A\right)=\int_{\left(B\right)}\ ,\ \ P\left(B\right)=\int_{\left(-A\right)}
ですが、P\left(B\right)が周期ならば-P\left(B\right)もまた周期だから、今ここでは二つの積分\displaystyle\int_{\left(B\right)}及び\displaystyle\int_{\left(A\right)}を考えることにします。A及びBは前に述べた通り、閉曲線です。しかし積分路は被積分関数の正則な範囲内では任意に変形しても積分の値に影響がないため、今、Aの代わりに二点1,-1を結ぶ線分を二重にとることとし(-1から1に行き、また-1に帰ると考える)Bの代わりに二点-1、\displaystyle-\frac{1}{k}を結ぶ線分(すなわちR面の交線)を一度はR面の一葉において一度は他葉において前と反対の方向に進むと考えた経路をとることとします。そうすれば
\int_{\left(A\right)}=2\int_{-1}^1\frac{dz}{\sqrt{\left(1-z^2\right)\left(1-k^2z^2\right)}}=4\int_0^1\frac{dz}{\sqrt{\left(1-z^2\right)\left(1-k^2z^2\right)}}
また
\int_{\left(B\right)}=2\int_{-\frac{1}{k}}^{-1}\frac{dz}{\sqrt{\left(1-z^2\right)\left(1-k^2z^2\right)}}=2\int_1^\frac{1}{k}\frac{dz}{\sqrt{\left(1-z^2\right)\left(1-k^2z^2\right)}}
ここで
z=\frac{1}{\sqrt{1-k^{‘2}t^2}}\ ,\ \ k^{‘2}=1-k^2
とおけば、
\int_{\left(B\right)}=2i\int_0^1\frac{dt}{\sqrt{\left(1-t^2\right)\left(1-k^{‘2}t^2\right)}}
となります。よって、二つの積分
\int_0^1\frac{dz}{\sqrt{\left(1-z^2\right)\left(1-k^2z^2\right)}}\ ,\ \ \int_0^1\frac{dz}{\sqrt{\left(1-z^2\right)\left(1-k^{‘2}z^2\right)}}
の値をそれぞれK、K’で表せば、
4K\ ,\ \ 2iK’
の整数倍の和はすべて\mathrm{sn}\ uの周期です。
今、原点と4K及び2iK’の点をそれぞれ結ぶ二つの線分を二辺として平行四辺形を作り、その中で\mathrm{sn}\ uのとる値の分布状態を調べてみましょう。
(1)によって明らかなように、zを-zに変えればuは-uとなります。ゆえに
\mathrm{sn}\left(-u\right)=-z=-\mathrm{sn}\ u\tag{2}
すなわち\mathrm{sn}は奇関数です。今この平行四辺形内の任意の一点をuとすれば、-uは図形内にはないが、これと同じu’=4K+2iK’-uは図形内にあって、uとu’はちょうど平行四辺形の中心2K+iK’に関して対称の位置にあります。そうしてこの二点において\mathrm{sn}\ u’=-\mathrm{sn}\ uの関係が成立するので、平行四辺形内の値の分布を調べるには全図形の半分(例えば2Kと中心を通る直線で切断したもの)だけについて調べればよいので、これを中心の周りに180^\circ回転して、かつその値にすべて-1を掛ければ残りの半分の分布を得ることになります。
次は図形内のいかなる二点おいて\mathrm{sn}\ uの値が等しくなるかを考えてみましょう。それには(1)においてzの同一の値にいかなるuの値が対応するかを調べればよいことになります。zの一つの値(\pm1,\displaystyle\pm\frac{1}{k}を除く)にはR面上の二点が対応します。確かにR面は二枚の数平面から成るものでその各葉の上にそれぞれ与えられたzの値を代表する点があるからです。(分岐点\pm1,\displaystyle\pm\frac{1}{k}だけにおいては両葉が融合して各唯一の点となっている。)今(1)の積分において、ある一葉の上の点0から同葉の点zにいたる種々の積分路に対する値を
u+4mK+2niK’\hspace{1cm}\left(m,nは整数\right)
とすれば、同じ0から他葉の点zにいたる積分は
2K-u+4mK+2niK’
で表されます(第9回参照)。なぜならば、0から他葉のzに行くにはまず0から1に行き、1の周りを小さく回ってR面の他葉に入り、したがって被積分関数の根号の符号を変えて1から0(前の0とは違う、他葉の0)に行き、この0からzに行けばよいです。これを式で書けば
\begin{eqnarray*}
\int_0^1\frac{dz}{\sqrt{\hspace{5mm}}}+\int_1^0\frac{dz}{-\sqrt{\hspace{5mm}}}+\int_0^z\frac{dz}{-\sqrt{\hspace{5mm}}}&=&2\int_0^1\frac{dz}{\sqrt{\hspace{5mm}}}-\int_0^z\frac{dz}{\sqrt{\hspace{5mm}}}\\
&=&2K-\int_0^z\frac{dz}{\sqrt{\hspace{5mm}}}
\end{eqnarray*}
よって前述のような値を得ることになります。
同一のzに対するuの値は以上の他にありませんが、この中で今考える平行四辺形内に入るものは二つしかありません。その一つをuとすれば他の一つはu^{”}=2K+2iK’-uです。この二点はK+iK’に関して互いに対称の位置にあります。ゆえに結局\left(0,2K\right)及び\left(0,iK’\right)を二辺とする平行四辺形(全図形の四分の一)内の値の分布が分かれば、対称性より全平行四辺形内の模様が分かることになります。
\mathrm{sn}\ u=\mathrm{sn}\ u^{”}=\mathrm{sn}\left(2K+2iK’-u\right)=\mathrm{sn}\left(2K-u\right)\ ,\hspace{1cm}すなわち
\mathrm{sn}\left(2K-u\right)=\mathrm{sn}\ u\tag{3}
これは重要な公式です。ここでuに-uを代入し、(2)を利用すれば
\mathrm{sn}\left(2K+u\right)=-\mathrm{sn}\ u\tag{4}
またKの定義により
\mathrm{sn}\ K=1
これと\mathrm{sn}\ 0=0から、(4)により
\mathrm{sn}\ 2K=0\ ,\ \ \mathrm{sn}\ 3K=-1\ ,\ \ \mathrm{sn}\ 4K=0
上の議論によって分かったように、\mathrm{sn}\ uは一つの平行四辺形内で同じ値を二度しか取りません。ゆえに\mathrm{sn}\ u零点は0,2Kの二つだけです。(4Kは0と合同であるから別に数えない。)
では\mathrm{sn}\ uの極はどこにあるのでしょうか、これを求めるには積分
U=\int_0^\infty\frac{dz}{\sqrt{\left(1-z^2\right)\left(1-k^2z^2\right)}}
を計算しなければなりません。ここに
z=\frac{it}{\sqrt{1-t^2}}
の置換を行えばただちに
U=i\int_0^1\frac{dt}{\sqrt{\left(1-t^2\right)\left(1-k^{‘2}t^2\right)}}=iK’
を得ます。ゆえに
\mathrm{sn}\ iK’=\infty
したがってまた(4)により
\mathrm{sn}\left(2K+iK’\right)=\infty
ゆえに極はiK’及び2K+iK’の二点です。
ここでついでに前の\displaystyle\int_{\left(B\right)}の計算の結果を思い出してみると
\int_1^\frac{1}{k}\frac{dz}{\sqrt{\left(1-z^2\right)\left(1-k^2z^2\right)}}=iK’
でした。この左辺を書き直して
\int_0^\frac{1}{k}-\int_0^1=\int_0^\frac{1}{k}-K
としてみると
\int_0^\frac{1}{k}\frac{dz}{\sqrt{\left(1-z^2\right)\left(1-k^2z^2\right)}}=K+iK’
ゆえに
\mathrm{sn}\left(K+iK’\right)=\frac{1}{k}
したがって(4)により
\mathrm{sn}\left(3K+iK’\right)=-\frac{1}{k}
cn, dn 関数
再び前回の(1)を書くと、
u=\int_0^z\frac{dz}{\sqrt{\left(1-z^2\right)\left(1-k^2z^2\right)}}\tag{1}
したがって
\frac{du}{dz}=\frac{1}{\sqrt{\left(1-z^2\right)\left(1-k^2z^2\right)}}
ゆえに
\frac{d}{du}\mathrm{sn}\ u=\frac{dz}{du}=\sqrt{\left(1-z^2\right)\left(1-k^2z^2\right)}=\sqrt{\left(1-\mathrm{sn}^2u\right)\left(1-k^2\mathrm{sn}^2u\right)}\tag{2}
の結果を得ます。ここで今
\sqrt{1-\mathrm{sn}^2u}\tag{3}
の関数を考えてみましょう。これは根号を帯びているから二価関数のように見えます。分岐点の位置は、
1-\mathrm{sn}^2u=0
すなわち
u=K,3K
でなければなりません。そこで試しにu=Kを中心として(3)をTaylor級数に展開してみることにします。まず
\varphi\left(u\right)=1-\mathrm{sn}^2u
とおき、(2)を利用して計算すれば
\begin{eqnarray*}
\varphi’\left(u\right)&=&-2\mathrm{sn}\ u\sqrt{\left(1-\mathrm{sn}^2u\right)\left(1-k^2\mathrm{sn}^2u\right)}\\
\varphi^{”}\left(u\right)&=&-2\left\{1-2\left(1+k^2\right)\mathrm{sn}^2u+3k^2\mathrm{sn}^4u\right\}等
\end{eqnarray*}
したがって
\varphi\left(K\right)=0,\ \varphi’\left(K\right)=0,\ \varphi^{”}\left(K\right)=2k^{‘2}等
ゆえに
\varphi\left(u\right)=k^{‘2}\left(u-K\right)^2+c_4\left(u-K\right)^4+\cdots\left(c_4,\cdots は定数\right)
の形の展開式を得ます。よって
\begin{eqnarray*}
\sqrt{1-\mathrm{sn}^2u}&=&k’\left(u-K\right)\sqrt{1+\frac{c_4}{k^{‘2}}\left(u-K\right)^2+\cdots}\\
&=&\pm k’\left(u-K\right)\left\{1+\frac{c_4}{2k^{‘2}}\left(u-K\right)^2+\cdots\right\}
\end{eqnarray*}
すなわちこれによるとu=Kは(3)の関数の分岐点ではありません。二つの離れた分枝が偶然に共に0の値をとるに過ません。他のu=3Kの点についても全く同様の結果を得ます。ゆえに(3)は二価関数ではない、実は二つの一価関数を表すものであることが分かります。(例えば\sqrt{1-\sin^2u}が\pm\cos uを表す。)
ここでこの関数を
\sqrt{1-\mathrm{sn}^2u}=\mathrm{cn}\ u
とします。ただし根式の二つの値の中でu=0のときに+1となる方をとるものです。全く同様の議論により
\sqrt{1-k^2\mathrm{sn}^2u}=\mathrm{dn}\ u
もまた一価関数であることが証明できます(やはり\mathrm{dn}\ 0=1とする)。\mathrm{cn},\ \mathrm{dn}はそれぞれ「シーエヌ」、「ディーエヌ」と読みます。\mathrm{sn}と同様にこれも略字ではなく、ただの名前です。
これらの新しい関数を用いて(2)を書き直せば次の式を得ます。
\frac{d}{du}\mathrm{sn}\ u=\mathrm{cn}\ u\ \mathrm{dn}\ u\tag{4}
\mathrm{cn}\ u、\mathrm{dn}\ uを微分すれば次の結果を得ます。
\left.\begin{array}{l}
\frac{\displaystyle d}{\displaystyle du}\mathrm{cn}\ u=\frac{\displaystyle d}{\displaystyle du}\sqrt{1-\mathrm{sn}^2u}=\frac{\displaystyle -\mathrm{sn}\ u\ \mathrm{cn}\ u\ \mathrm{dn}\ u}{\displaystyle \sqrt{1-\mathrm{sn}^2 u}}=-\mathrm{sn}\ u\ \mathrm{dn}\ u\\
\frac{\displaystyle d}{\displaystyle du}\mathrm{dn}\ u=\frac{\displaystyle d}{\displaystyle du}\sqrt{1-k^2\mathrm{sn}^2u}=\frac{\displaystyle -k^2\mathrm{sn}\ u\ \mathrm{cn}\ u\ \mathrm{dn}\ u}{\displaystyle \sqrt{1-k^2\mathrm{sn}^2 u}}=-k^2\mathrm{sn}\ u\ \mathrm{cn}\ u
\end{array}\right\}\tag{5}
(4)、(5)を反復使用すれば高次の導関数も求められます。したがってまた次のような展開式を作ることも出来ます。
\begin{eqnarray*}
\mathrm{sn}\ u&=&u-\left(1+k^2\right)\frac{u^3}{3!}+\left(1+14k^2+k^4\right)\frac{u^5}{5!}-\cdots\\
\mathrm{cn}\ u&=&1-\frac{u^2}{2!}+\left(1+4k^2\right)\frac{u^4}{4!}-\cdots\\
\mathrm{dn}\ u&=&1-k^2\frac{u^2}{2!}+\left(4k^2+k^4\right)\frac{u^4}{4!}-\cdots
\end{eqnarray*}
参考文献
参考文献は以下の通り。
[1]竹内端三,『楕円関数論』,岩波書店,1936
出版社在庫無し、著作権消失済み。
[2]E.T. Whittaker, et al., A Course of Modern Analysis (AMS PRESS, 1927)
著作権消失済み。
[3]戸田盛和,『楕円関数入門』,日本評論社,2001
[4]戸田盛和,『臨時別冊・数理科学SGC ライブラリ49 ソリトンと物理学』,サイエンス社,2006
同出版社より電子書籍の形で復刊済み。
[5]Landau・Lifshitz,『力学』,東京図書,2017