不完全ガンマ関数

sの値のいくつかに対する上側不完全ガンマ関数:0(青)、1(赤)、2(緑)、3(オレンジ)、4(紫)。
Mathematica 13.1 の ComplexPlot3D 関数を使用して作成した、複素平面における正規化不完全ガンマ関数 Q(2,z) の -2-2i から 2+2i までの色付きプロット
Mathematica 13.1 の ComplexPlot3D 関数を使用して作成した、複素平面における正規化不完全ガンマ関数 Q(2,z) の -2-2i から 2+2i までの色付きプロット

数学において、上側不完全ガンマ関数と下側不完全ガンマ関数は、特定の積分などのさまざまな数学の問題の解として生じる特殊関数の一種です。

それぞれの名称は、積分の定義に由来しています。ガンマ関数と同様に定義されますが、積分極限はガンマ関数と異なります(つまり「不完全」です)。ガンマ関数は、ゼロから無限大までの積分として定義されます。これは、ゼロから可変の上限までの積分として定義される下側不完全ガンマ関数とは対照的です。同様に、上側不完全ガンマ関数は、可変の下限から無限大までの積分として定義されます。

定義

上側不完全ガンマ関数は次のように定義されます。 一方、下側不完全ガンマ関数は次のように定義されます。 どちらの場合も、 sは複素パラメータであり、 sの実部は正です

性質

部分積分により、漸化式 が求められます 。 通常のガンマ関数は 次 のように定義されるので 、

複素数値への継続

上記で定義した実数sxに対する下側不完全ガンマ関数と上側不完全ガンマ関数は、 xsの両方に関して、複素数xsのほぼすべての組み合わせに対して定義される正則関数に展開できます。[ 1 ]複素解析は、実数不完全ガンマ関数の特性が、その正則関数にどのように拡張されるかを示しています。

下側不完全ガンマ関数

正則拡大

下側不完全ガンマ関数に対する漸化式を繰り返し適用すると、次式で表されるべき級数展開が得られる。[ 2 ] k → ∞のときにΓ( z + k )絶対値 が急激に増加すること、およびΓ( z )の逆数が整関数であることから、右端の和の係数は明確に定義され、局所的にはすべての複素sおよびxについて和は一様に収束する。ワイエルシュトラスの定理により、[ 3 ]極限関数 (と表記されることもある) [ 4 ]は、z (固定のsの場合) およびs (固定のzの場合)の両方に関してで あり、 [ 1 ]ハートッグスの定理によりC × C上で正則である。[ 5 ]したがって、次の分解[ 1 ]は、 実下側不完全ガンマ関数を、 zおよびsについて、共同でも個別にも、正則関数として拡張する。およびΓ関数の特性から、最初の 2 つの因子が(z = 0またはsが正でない整数) の特異点を捉え、最後の因子がその零点に寄与することがわかります。

多価性

複素対数log z = log | z | + i arg zは2πiの倍数までしか決まらないため、多価となります。複素対数を含む関数は、通常この性質を継承します。その中には、複素べき乗や、その分解にz s現れるγ関数など があります

多価関数の不確定性は、値の選択方法を明示的に指定する必要があるため、複雑な問題を引き起こします。これに対処するための戦略は以下のとおりです。

  • (最も一般的な方法)多価関数の領域Cを、 C × C内の適切な多様体、すなわちリーマン面によって置き換える。これにより多価性は解消されるが、その背後にある理論を理解する必要がある。[ 6 ]
  • 多値関数が個別の単一値ブランチに分解され、個別に処理できるようにドメインを制限します。

このセクションの式を正しく解釈するには、以下の規則に従う必要があります。特に断りのない限り、以下の規則が適用されます。

セクター

Cにおけるz = 0を頂点とするセクターは、複素表現の適切な定義域であることがしばしば証明されています。セクターD は、 z ≠ 0かつαδ < arg z < α + δ (ただし、 α0 < δπ )を満たすすべての複素zから構成されます。多くの場合、α は任意に選択でき、その場合は指定されません。δ が指​​定されていない場合πであると仮定され、セクターは実際には平面C全体ですが、 z = 0を起点としてαの方向を指す半直線(通常は分岐切断線として機能します)は例外です。注:多くの応用例やテキストでは、αは暗黙的に 0 とみなされ、セクターは正の実軸を中心に配置されます

支店

特に、そのようなセクター D 上には、虚数部が範囲( αδα + δ )に制限された単値の正則対数が存在する。このような制限された対数に基づくと、z sおよび不完全ガンマ関数は、 D (またはC × D )上の単値の正則関数に縮小され、 D 上の多値関数の枝と呼ばれる。α2 πの倍数を追加すると、同じセットD上に相関する枝の異なるセットが得られる。ただし、ここでの特定のコンテキストでは、αは固定であると想定され、関連するすべての枝はそれに関連付けられている。| α | < δの場合、枝は、正の実軸上の実際の類似物に等しいため、枝と呼ばれる。注: 多くのアプリケーションおよびテキストでは、公式は主枝に対してのみ有効です。

枝間の関係

複素べき乗関数と下側不完全ガンマ関数の両方の異なる枝の値は、 kを適切な整数 として、[ 1 ]を乗算することによって互いに導くことができます

分岐点付近の挙動

上記の分解はさらに、γ がz = 0付近で漸近的に次のように振舞うことを示しています。

正の実数xysに対して、( xy ) → (0、s )のとき、 x y /y → 0となります。これは、実数s > 0に対してγ ( s、 0) = 0と設定することを正当化するように見えます。しかし、複素数領域では状況が多少異なります。 (a) sの実部が正であり、 (b) 値u v が有限の枝の集合からのみ取得される場合に限り、それらは( uv ) → (0、s )としてゼロに収束することが保証され、 γ ( uv )も同様です。 γ ( b )の単一のでは自然に満たされるため、正の実部を持つsに対してγ ( s、 0) = 0連続極限となります。また、このような連続は決して解析的なものではないことにも注意してください。

代数関係

実数γ ( s , z )に見られるすべての代数関係と微分方程式は、その正則関数にも成り立ちます。これは、実区間上で有効な正則関数間の方程式はどこでも成り立つという恒等定理の帰結です。特に、再帰関係[ 2 ]∂γ ( s , z )/ ∂z = z s −1 e z [ 2 ]は、対応する枝において保存されます

積分表現

最後の関係式は、s を固定した場合、γ は正則関数z s −1 e zの原始微分または反微分であることを示しています。したがって、任意の複素uv ≠ 0に対して、 積分経路が被積分関数の枝の領域に完全に含まれる 限り、が成立します。さらに、 sの実部が正の場合、 u → 0に対して極限γ ( s , u ) → 0が適用され、最終的にγの複素積分定義に到達します[ 1 ]

ここでは、先頭にのみ 0 を含む積分パス(それ以外は積分対象のブランチのドメインに制限)が有効です。たとえば、 0zを結ぶ直線などです。

z → +∞の限界
実数値

γの主枝の積分表現が与えられた場合、すべての正の実数sxに対して次の式が成り立つ:[ 7 ]

複雑

この結果は複素数sにも拡張される。まず1 ≤ Re( s ) ≤ 2かつ1 < a < bと仮定する。すると と なる。 ここで[ 8 ] は途中で用いられている。a が十分に大きい場合、最終積分は任意に小さくなるため、 γ ( s , x ) はストリップ 1 ≤ Re(s) ≤ 2 上で x → ∞ に対して一様収束し、正則関数となる [ 3 ]恒等定理よりこの関数Γ ( s )なければならない。再帰関係γ ( s , x ) = ( s − 1) γ ( s − 1, x ) − x s − 1 e xの極限をとり、 x → ∞およびすべてのnに対してlim x n e x = 0となることから、γ ( s , x )はストリップの外側でも、Γ関数の再帰関係に従う関数へと収束することが示される。これは、 s が非正整数でなく、x が実数で、γ が主関数で あるすべての複素数に対して成り立つ。

セクターごとの収束

ここで、uをセクター| arg z | < δ < π /2から、 δを固定値( α = 0 ) として、γ をこのセクターの主枝とし、

上で示したように、 | u |が十分に大きい場合、最初の差は任意に小さくすることができます。 2 番目の差により、次の推定が可能になります。 ここで、 γの積分表現と上記の| z s |に関する式 を使用しました。 u| u |を結ぶ半径R = | u |の円弧に沿って 0 の周りを積分すると、最後の積分は次のようになります。 ここで、M = δ (cos δ ) −Re s e Im はuRに依存しない定数です。 再び、大きなxに対するx n e xの挙動を参照すると、 R がに向かって増加するにつれて、最後の式が 0 に近づくことがわかります。 合計すると、次のようになります。 sが非負の整数でない 場合、 0 < ε < π /2は任意に小さいですが固定であり、γ はこのドメインの主枝を示します。

概要

は:

  • 固定された正の整数sに対してz整数
  • 整数でない固定のsに対してzの多値正則性を持ち、分岐点z = 0です。
  • 各枝は固定のz ≠ 0に対してsに関して有理型であり、非正の整数 s に単純な極を持つ。

上側不完全ガンマ関数

上側不完全ガンマ関数に関しては、右辺が存在する 点( s , z )において、 zまたはsに関する正則拡大が[ 1 ]によって与えられる。は多価関数であるため、 についても同様であるが、主値に制限すると の主枝は一価のみとなる。

上式においてsが非正の整数である場合、差のどちらの部分も定義されておらず、ここでs → 0に対して展開される極限過程によって欠損値が補われる。複素解析は が固定されたzに対してその極限の近傍有界であることが証明されるため、正則性 を保証する。

極限を決定するには、 z = 0におけるのべき級数が有用である。の積分定義において のべき級数を で置き換えると、次式が得られる(ここではxs は正の実数と仮定する)。 または[ 4 ] これは、関数全体の級数表現として、すべての複素数x(および非正整数でない すべての複素数s )に対して収束する。

実数値への制限がなくなると、この級数は次のように展開できるようになります。

s → 0 のとき: [ 9 ] (ここではオイラー・マスケロニ定数)したがって、 はs → 0 の上側不完全ガンマ関数の極限関数であり、指数積分としても知られています[ 10 ]

再帰関係により、正の整数nに対するの値は、この結果から導くことができる。[ 11 ] したがって、上側不完全ガンマ関数は、すべてのsおよびz ≠0に対して、 zsの両方に関して存在し、正則であることが証明される。

は:

  • 固定された正の整数sに対してz軸方向に整数
  • 固定されたsに対してzの多値正則性があり、分岐点z = 0です。
  • は、実部が正でz = 0 ( のときの極限) のsに対してと等しくなりますが、これは連続的な拡張であり、解析的な拡張ではありません(は実数s < 0に対しては成立しません)。
  • 各枝上で、z ≠0を固定してs全体になります。

特別な値

  • sが正の整数の場合、
  • sが正の整数の場合、[ 12 ]
  • について、

ここで、は指数積分、は一般指数積分、は誤差関数、は 相補誤差関数です

漸近的挙動

  • として、
  • として、(実数sに対して、Γ( s , x )~ −xs / sの誤差は、 s ≠−1の場合はO ( xmin { s +1,0} )s =−1の場合はO (ln( x ))のオーダーである)、
  • 漸近級数として表される。[ 13 ]
  • 漸近級数として、および、、、ここでオイラー・マスケローニ定数である。[ 13 ]
  • として、
  • として、
  • 漸近級数として表される。[ 14 ]

評価式

下ガンマ関数は、べき級数展開を用いて評価することができます。[ 15 ]ここ で、はポッホハマー記号です

別の展開は、 Mがクンマーの合流型超幾何関数 であるときです。

クンマーの合流型超幾何関数との関連

zの実部が正のとき、 は無限大の収束半径を持ちます。

再び合流超幾何関数とクンマーの恒等式を用いると、

実際の数値の計算には、ガウスの連分数が便利な展開を提供します。

この連分数は、sが負の整数でない 限り、すべての複素数zに対して収束します。

上側のガンマ関数は連分数[ 16 ] を持ち 、

乗法定理

次の乗法定理が成り立ちます

ソフトウェア実装

不完全ガンマ関数は、さまざまなコンピュータ 代数システムで利用できます

直接的に利用できない場合でも、スプレッドシート(​​およびコンピュータ代数パッケージ)に一般的に含まれる関数を使用することで、不完全な関数値を計算することは可能です。例えばExcelでは、ガンマ関数とガンマ分布関数を組み合わせることで計算できます。

  • 下側の不完全関数: 。 = EXP(GAMMALN(s))*GAMMA.DIST(x,s,1,TRUE)
  • 上の不完全関数: 。 = EXP(GAMMALN(s))*(1-GAMMA.DIST(x,s,1,TRUE))

これらはガンマ分布の累積分布関数の定義から導かれます。

Pythonでは、Scipy ライブラリが の不完全ガンマ関数の実装を提供していますが、最初の引数に負の値を指定することはできません。mpmathライブラリの scipy.special関数は、すべての複素引数をサポートします。gammainc

正規化ガンマ関数とポアソン確率変数

関連する 2 つの関数は、正規化されたガンマ関数です。 これは、形状パラメータースケール パラメーターが 1 であるガンマ ランダム変数累積分布関数です。

が整数のとき、 はポアソン確率変数の累積分布関数である。が確率変数のとき、

この式は部分積分を繰り返すことで導き出すことができます。

これらはscipyでは[ 17 ][ 18 ]のように実装されている。 gammaincgammaincc

デリバティブ

上記の積分表現を使用すると、上不完全ガンマ関数のxに関する導関数は、 で表されます。 最初の引数に関する導関数は[ 19 ]で与えられ 、2 番目の導関数は で 表されます。 ここで、関数はMeijer G 関数の特殊なケースです。 この特殊なケースは、すべての連続する導関数を表現できるため、独自の内部閉包特性を持っています。 一般に、 はポッホハマー記号で定義される置換です 。この ようなすべての導関数は、次から連続して生成できます。 および この関数は、 が負の整数またはゼロではない ことを理解すれば、 に対して有効な級数表現から計算できます。 このような場合、極限を使用する必要があります。 の結果は、解析接続によって取得できます。 この関数のいくつかの特殊なケースは簡略化できます。たとえば、、で、は指数積分です。 これらの導関数と関数は、上不完全ガンマ関数の積分定義を繰り返し微分することにより、多数の積分の正確な解を提供します。[ 20 ] [ 21 ] 例えば、 この式はさらに拡張または一般化され、ラプラス変換メリン変換の膨大なクラスに拡張することができます。コンピュータ代数システムと組み合わせることで、特殊関数の利用は、特に実用的な工学アプリケーションで遭遇する定積分を解くための強力な手段となります(詳細については 記号積分を参照してください)。

不定積分と定積分

次の不定積分は、部分積分(どちらの場合も 積分定数は省略)を使用して簡単に得られます。 下側不完全ガンマ関数と上側不完全ガンマ関数は、フーリエ変換を介して接続されます。 これは、たとえば、(Gradshteyn et al. 2015、§7.642)の適切な特殊化によって得られます。

注釈

  1. ^ a b c d e f 「DLMF:§8.2 定義と基本的性質 ‣ 不完全ガンマ関数 ‣ 第8章 不完全ガンマ関数と関連関数」dlmf.nist.gov
  2. ^ a b c「DLMF: §8.8 再帰関係と導関数 ‣ 不完全ガンマ関数 ‣ 第8章 不完全ガンマ関数と関連関数」 . dlmf.nist.gov .
  3. ^ a b Donald E. Marshall (2009年秋). 「複素解析」(PDF) .数学534 (学生用配布資料). ワシントン大学. 定理3.9、p.56.オリジナル(PDF)から2011年5月16日時点のアーカイブ。 2011年4月23日閲覧
  4. ^ a b「DLMF: §8.7 級数展開 ‣ 不完全ガンマ関数 ‣ 第8章 不完全ガンマ関数と関連関数」 . dlmf.nist.gov .
  5. ^ Paul Garrett. 「ハートッグスの定理:個別解析性は結合解析性を意味する」(PDF) . cse.umn.edu . 2023年12月21日閲覧
  6. ^ C. Teleman. 「リーマン面」(PDF) . berkeley.edu . 2023年12月21日閲覧
  7. ^ 「DLMF: §5.2 定義 ‣ プロパティ ‣ 第5章 ガンマ関数」 . dlmf.nist.gov .
  8. ^ 「DLMF: §4.4 特殊な値と限界 ‣ 対数、指数、累乗 ‣ 第4章 基本関数」 . dlmf.nist.gov .
  9. ^最後の式を参照してください。
  10. ^ 「DLMF: §8.4 特殊値 ‣ 不完全ガンマ関数 ‣ 第8章 不完全ガンマ関数と関連関数」 . dlmf.nist.gov .
  11. ^ 「DLMF: 8.4 特殊な値」
  12. ^ Weisstein, Eric W.不完全ガンマ関数」MathWorld .(式2)
  13. ^ a b Bender & Orszag (1978). Advanced Mathematical Methods for Scientists and Engineers . Springer. Bibcode : 1978amms.book..... B
  14. ^ 「DLMF: §8.11 漸近近似と展開 ‣ 不完全ガンマ関数 ‣ 第8章 不完全ガンマ関数と関連関数」 . dlmf.nist.gov .
  15. ^ 「DLMF: §8.11 漸近近似と展開 ‣ 不完全ガンマ関数 ‣ 第8章 不完全ガンマ関数と関連関数」 . dlmf.nist.gov .
  16. ^アブラモヴィッツとステガンp. 263、6.5.31
  17. ^ "scipy.special.gammainc — SciPy v1.11.4 マニュアル" . docs.scipy.org .
  18. ^ "scipy.special.gammaincc — SciPy v1.11.4 マニュアル" . docs.scipy.org .
  19. ^ KO Geddes、ML Glasser、RA Moore、TC Scott、「特殊関数の微分による基本関数を含む定積分のクラスの評価」、AAECC(工学、通信、コンピューティングにおける適用可能代数)、第1巻、(1990年)、pp. 149–165、 [1]
  20. ^ Milgram, MS (1985). 「一般化積分指数関数」 . Math. Comp . 44 (170): 443– 458. doi : 10.1090/S0025-5718-1985-0777276-4 . MR 0777276 . 
  21. ^ Mathar (2009). 「1から無限大までのexp(i*pi*x)*x^(1/x)上の振動積分の数値評価」arXiv : 0912.3844 [ math.CA ].、App B

参考文献