ガウスの連分数

提供: testwiki
ナビゲーションに移動 検索に移動

複素解析におけるガウスの連分数(ガウスのれんぶんすう、テンプレート:Lang-en-short)は、超幾何関数から導出される特別なクラスのテンプレート:仮リンクである。これは数学史上最も早く見出された解析的な連分数の一つであり、いくつかの重要な初等関数およびより複雑な超越関数の表現に用いることができる。

歴史

ヨハン・ハインリヒ・ランベルトはこの形式の連分数のいくつかの例を1768年に発表し、またレオンハルト・オイラージョゼフ=ルイ・ラグランジュは同様の構造についての研究を行った[1]。しかし、以下の節に記すような算法を基にして、これらの連分数の一般的な形式を導いたのはカール・フリードリヒ・ガウスであった(1813年)[2]

ガウスは連分数の形式を与えはしたが、その収束性についての証明は与えなかった。ベルンハルト・リーマン[3]と L.W.トーメ[4]は部分的な結果を得ていたが、これらの連分数の収束性に関して最終的な結論がまとめられたのは1901年、テンプレート:仮リンク[5]によってであった。

導出

f0,f1,f2,解析関数の列で、任意の i>0 に対し

fi1fi=kizfi+1

を満たすものとする。ここで各 ki は定数である。

このとき

fi1fi=1+kizfi+1fi, and so fifi1=11+kizfi+1fi

gi=fi/fi1 とおくと

gi=11+kizgi+1

より、

g1=f1f0=11+k1zg2=11+k1z1+k2zg3=11+k1z1+k2z1+k3zg4=. 

これを無限に続けると、連分数展開

f1f0=11+k1z1+k2z1+k3z1+

が得られる。

ガウスの連分数では、関数列 fi0F1, 1F1, 2F1 の形の超幾何関数とする。等式 fi1fi=kizfi+1 は整数部分だけずらした関数間の恒等式から生じる。これらの恒等式は、級数展開をして係数を比較したり、何度か微分して消えることを確かめるといった方法で示すことができる。

0F1 による展開

最も簡単なものは、関数

0F1(a;z)=1+1a1!z+1a(a+1)2!z2+1a(a+1)(a+2)3!z3+

を使った恒等式

0F1(a1;z)0F1(a;z)=za(a1)0F1(a+1;z)

から始めるもので、

fi=0F1(a+i;z),ki=1(a+i)(a+i1)

と置けば次の展開が得られる。

0F1(a+1;z)0F1(a;z)=11+1a(a+1)z1+1(a+1)(a+2)z1+1(a+2)(a+3)z1+

または

0F1(a+1;z)a0F1(a;z)=1a+z(a+1)+z(a+2)+z(a+3)+

これらの展開は、2つの収束級数の比による有理型関数に収束する(ただし当然、a は0や負の整数であってはいけない)。

1F1 による展開

次は、関数

1F1(a;b;z)=1+ab1!z+a(a+1)b(b+1)2!z2+a(a+1)(a+2)b(b+1)(b+2)3!z3+

から生じる2通りの恒等式

1F1(a;b1;z)1F1(a+1;b;z)=(ab+1)zb(b1)1F1(a+1;b+1;z)
1F1(a;b1;z)1F1(a;b;z)=azb(b1)1F1(a+1;b+1;z)

を交互に用いて、

f0(z)=1F1(a;b;z),
f1(z)=1F1(a+1;b+1;z),
f2(z)=1F1(a+1;b+2;z),
f3(z)=1F1(a+2;b+3;z),
f4(z)=1F1(a+2;b+4;z)

…と置く。これより恒等式群 fi1fi=kizfi+1 (ここで k1=abb(b+1),k2=a+1(b+1)(b+2),k3=ab1(b+2)(b+3),k4=a+2(b+3)(b+4), …)が定まって次の展開が得られる。

1F1(a+1;b+1;z)1F1(a;b;z)=11+abb(b+1)z1+a+1(b+1)(b+2)z1+ab1(b+2)(b+3)z1+a+2(b+3)(b+4)z1+

または

1F1(a+1;b+1;z)b1F1(a;b;z)=1b+(ab)z(b+1)+(a+1)z(b+2)+(ab1)z(b+3)+(a+2)z(b+4)+

同様に

1F1(a;b+1;z)1F1(a;b;z)=11+ab(b+1)z1+ab1(b+1)(b+2)z1+a+1(b+2)(b+3)z1+ab2(b+3)(b+4)z1+

または

1F1(a;b+1;z)b1F1(a;b;z)=1b+az(b+1)+(ab1)z(b+2)+(a+1)z(b+3)+(ab2)z(b+4)+

1F1(0;b;z)=1 であるから、最初の連分数で a を0と置き、b + 1 を b で置き換えると次の単純化された特別な連分数が得られる。

1F1(1;b;z)=11+zb+z(b+1)+bz(b+2)+2z(b+3)+(b+1)z(b+4)+

2F1 による展開

最後は

2F1(a,b;c;z)=1+abc1!z+a(a+1)b(b+1)c(c+1)2!z2+a(a+1)(a+2)b(b+1)(b+2)c(c+1)(c+2)3!z3+

から、再び2つの恒等式を交互に用いる。

2F1(a,b;c1;z)2F1(a+1,b;c;z)=(ac+1)bzc(c1)2F1(a+1,b+1;c+1;z)
2F1(a,b;c1;z)2F1(a,b+1;c;z)=(bc+1)azc(c1)2F1(a+1,b+1;c+1;z)

これら2式は ab の入れ替えに対し本質的に不変である。

f0(z)=2F1(a,b;c;z),
f1(z)=2F1(a+1,b;c+1;z),
f2(z)=2F1(a+1,b+1;c+2;z),
f3(z)=2F1(a+2,b+1;c+3;z),
f4(z)=2F1(a+2,b+2;c+4;z),

…とおく。これより恒等式群 fi1fi=kizfi+1 (ここで k1=(ac)bc(c+1),k2=(bc1)(a+1)(c+1)(c+2),k3=(ac1)(b+1)(c+2)(c+3),k4=(bc2)(a+2)(c+3)(c+4), …)が定まって次の展開が得られる。

2F1(a+1,b;c+1;z)2F1(a,b;c;z)=11+(ac)bc(c+1)z1+(bc1)(a+1)(c+1)(c+2)z1+(ac1)(b+1)(c+2)(c+3)z1+(bc2)(a+2)(c+3)(c+4)z1+

または

2F1(a+1,b;c+1;z)c2F1(a,b;c;z)=1c+(ac)bz(c+1)+(bc1)(a+1)z(c+2)+(ac1)(b+1)z(c+3)+(bc2)(a+2)z(c+4)+

2F1(0,b;c;z)=1 であるから、a を 0 と置き、c + 1 を c で置き換えると次の単純化された特別な連分数が得られる。

2F1(1,b;c;z)=11+bzc+(bc)z(c+1)+c(b+1)z(c+2)+2(bc1)z(c+3)+(c+1)(b+2)z(c+4)+

収束性

この節では、1個以上のパラメータが負の整数である場合は考えないことにする。そのような場合、超幾何関数は定義されないか、または多項式に退化するために連分数展開が有限回で止まるからである。その他の自明な状況も排除するものとする。

0F11F1 の場合、級数は任意の点で収束し、左辺の分数関数は有理型関数になる。右辺の連分数はを含まない任意の有界な集合上一様に収束する[6]

2F1 の場合、展開の収束半径は 1 で、左辺の関数はこの円板の内部で有理型関数を表す。右辺の連分数はこの円板内部の任意の点で収束する。

円板の外部では、連分数は テンプレート:Math から無限遠点までを除いた実軸に沿って解析接続された関数を表す。テンプレート:Math分岐点テンプレート:Mathから無限遠点への実軸上の半直線が分岐截線とされることが多い。右辺の連分数はこの領域で有理型関数へ収束し、また極を含まない任意の有界閉集合上で収束は一様である[7]

応用

0F1 による展開

cosh(z)=0F1(12;z24)
sinh(z)=z0F1(32;z24)

より、

tanh(z)=z0F1(32;z24)0F1(12;z24)=z/212+z2432+z2452+z2472+=z1+z23+z25+z27+

この展開は特にランベルトの連分数として知られており、初出は1768年にまで遡る[8]

これより次式がすぐに従う。

tan(z)=z1z23z25z27

tanh の展開は、ネイピア数の任意の整数乗 en無理数であることを証明するのに使うことができる(しかし残念なことに、e超越数であることを証明するのには力不足である)。tan の展開はランベルトとアドリアン=マリ・ルジャンドルによって円周率の無理性の証明に用いられた。

ベッセル関数 Jν は次のように表示できる。

Jν(z)=(12z)νΓ(ν+1)0F1(ν+1;z24)

これより、

Jν(z)Jν1(z)=z2νz22(ν+1)z22(ν+2)z22(ν+3)

これらの等式は全ての複素数 z に対し成り立つ。

1F1 による展開

ez=1F1(1;1;z), 1/ez=ez より

ez=11+z1+z2+z3+2z4+2z5+
ez=1+z1+z2+z3+2z4+2z5+

これを使って少し計算すると、e の簡単な連分数表示が得られる。

e=2+11+12+11+11+14+

次式で与えられる誤差関数 erf (z)

erf(z)=2π0zet2dt

もまた、クンマーの超幾何関数として

erf(z)=2zπez21F1(1;32;z2)

と計算でき、ガウスの連分数表示を用いると全ての z に対し成り立つ有用な公式が得られる[9]

π2ez2erf(z)=z1z232+z25232z272+2z29252z2112+3z213272z2152+

同様の議論によってフレネル積分テンプレート:仮リンク不完全ガンマ関数に対する連分数展開表示が得られる。簡易化版の議論からは指数関数の2通りの有用な連分数展開表示が得られる[10]

2F1 による展開

(1z)b=1F0(b;;z)=2F1(1,b;1;z),
(1z)b=11+bz1+(b1)z2+(b+1)z3+2(b2)z4+

から、arctan z の0を中心としたテイラー展開が容易に求まる。

arctanz=zF(12,1;32;z2)

これにガウスの連分数表示を用いると、展開

arctanz=z1+(1z)23+(2z)25+(3z)27+(4z)29+

が得られる。この展開は、虚軸上 i から無限遠点まで、−i から無限遠点までが切り取られた複素平面において逆正接関数のテンプレート:仮リンクに収束する[11]

この連分数は z = 1 で非常に高速に収束し、第9番目までの展開で π/4 の値を小数第7位まで与える。対応する級数

π4=11+122+322+522+=113+1517±

の収束は非常に遅く、小数第7位まで正確に計算するには百万項以上を必要とする[12]

これと同様の論法で、自然対数逆正弦関数一般化二項級数といった関数の連分数表示が得られる。

脚注

テンプレート:脚注ヘルプ テンプレート:Reflist

参考文献

  1. Jones & Thron (1980) p. 5
  2. C. F. Gauss (1813), Werke, vol. 3 pp. 134–38.
  3. B. Riemann (1863), "Sullo svolgimento del quoziente di due serie ipergeometriche in frazione continua infinita"(『2つの超幾何級数の比の、無限連分数の形での展開について』) in Werke. pp. 400–406. (Posthumous fragment).
  4. L. W. Thomé (1867), "Über die Kettenbruchentwicklung des Gauß'schen Quotienten ...,"(『ガウスの比の連分数展開については…』) Jour. für Math. vol. 67 pp. 299–309.
  5. E. B. Van Vleck (1901), "On the convergence of the continued fraction of Gauss and other continued fractions." Annals of Mathematics, vol. 3 pp. 1–18.
  6. Jones & Thron (1980) p. 206
  7. Wall, 1973 (p. 339)
  8. Wall (1973) p. 349.
  9. Jones & Thron (1980) p. 208.
  10. テンプレート:仮リンクの記事にある ez の連分数展開例を参照。
  11. Wall (1973) p. 343. i と −i が逆正接の分岐点であることに注意。
  12. Jones & Thron (1980) p. 202.