Неповна гамма-функція

Матеріал з Вікіпедії — вільної енциклопедії.
Перейти до навігації Перейти до пошуку
Нижня неповна гамма-функція в 3D
Анімація нижньої неповної гамма-функції
Верхня неповна гамма-функція в 3D

У математиці верхня неповна гамма-функція і нижня неповна гамма-функція є типом спеціальних функцій, які виникають при розв'язанні різноманітних математичних задач, таких як деякі інтеграли.

Їх відповідні назви випливають з їх інтегральних визначень, які визначаються аналогічно гамма-функції, іншим типом спеціальної функції, але з різними або "неповними" інтегральними межами. Гамма-функція визначається як інтеграл від нуля до нескінченності. Це відрізняється від нижньої неповної гамма-функції, яка визначається як інтеграл від нуля до змінної верхньої межі. Відповідно, верхня неповна гамма-функція визначається як інтеграл від змінної нижньої межі до нескінченності.

Визначення[ред. | ред. код]

Верхня неповна гамма-функція визначається як:

в той час як нижня неповна гамма-функція визначається як:

Властивості[ред. | ред. код]

В обох випадках s є складним параметром, таким, що дійсна частина s є позитивною.

Інтегруванням частинами знаходимо рекурентне співвідношення

та

Оскільки звичайна гамма-функція визначається як

маємо

та

Продовження для комплексних значень[ред. | ред. код]

Нижня неповна гамма-функція та верхня неповна гамма-функція, що визначені вище для дійсних позитивних s і x, можуть бути розвинені в голоморфній функції s, по відношенню як до x, так і до s, визначені для майже всіх комбінацій складних x та s.[1] Комплексний аналіз показує, що властивості дійсних неповних гамма-функцій поширюються на їх голоморфні аналоги.

Голоморфне розширення[ред. | ред. код]

Повторне застосування рекурентного відношення до нижньої неповної гамма-функції призводить до розширення степеневого ряду: [2]

Враховуючи стрімке зростання абсолютної величини Γ(z + k) при k → ∞, а також той факт, що взаємна Γ( z ) є цілою функцією, коефіцієнти у правій частині суми є чітко визначеними, а локально сума сходиться рівномірно для всіх комплексних s та x. За теоремою Вейерштрасса, гранична функція іноді позначається як ,

[3]

є цілою щодо z (для фіксованих s) та s (для фіксованих z) [4] , і, таким чином, голоморфна на ℂ×ℂ за теоремою Хартогса[5]. Отже, наступне розвинення

[6],

розширює дійсну неповну гамма-функцію в голоморфну функцію, як спільно, так і окремо в z і s. З властивостей zs та гамма-функції випливає, що перші два чинники фіксують особливі точки γ (при z = 0 або s - не додатнє ціле число), тоді як останній чинник сприяє його нулям.

Багатозначність[ред. | ред. код]

Комплексний логарифм log z = log |z| + i arg z визначається лише кратним 2πi, що робить його багатозначним. Функції, що включають складний логарифм, як правило, успадковують цю властивість. Серед них - складна потужність, і оскільки z s з'являється в його розкладі, то γ-функція теж.

Невизначеність багатозначних функцій призводить до ускладнень, оскільки необхідно вказати, як вибрати значення. Стратегії для вирішення цього є:

  • (найбільш загальний спосіб) замінити домен ℂ багатозначних функцій відповідним різновидом у ℂ × ℂ, що називається поверхнею Рімана. У той час як це усуває багатозначність, потрібно знати теорію, яка лежить в основі [7];
  • обмежуйте такий домен, при якому багатозначна функція розкладається на окремі однозначні гілки, які можна обробляти індивідуально.
Поведінка поблизу точки розгалуження[ред. | ред. код]

Поряд з z=0 γ поводиться асимптотично:

Для додатних дійсних x, y та s, xy/y → 0, при (x, y) → (0, s). Цим обґрунтовується, що γ(s, 0) = 0 для дійсних s > 0. Проте, в комплексній області все трохи інакше. Тільки якщо (a) дійсна частина s є додатною, та (b) значення uv беруться з лише кінцевого набору гілок, то вони гарантовано сходяться до нуля, оскільки (u, v) → (0, s), так само і для γ(u, v). γ(b) виконується на одиничній гілці, і тому γ(s, 0) = 0 для s з додатною дійсною частиною є неперевною межею. Таке продовження в жодному разі не є аналітичним.

Алгебраїчні співвідношення[ред. | ред. код]

Усі алгебраїчні співвідношення та диференціальні рівняння, які спостерігаються за допомогою дійсної γ(s, z) зберігаються і для її голоморфного аналога. Це є наслідком теореми про тотожність[8], яка вказує, що рівняння між голоморфними функціями, чинними на дійсному інтервалі, зберігаються скрізь. Зокрема, у відповідних гілках, зберігаються рекурентні співвідношення [9] та ∂γ(s,z)/∂z = zs−1 ez [10] зберігаються на відповідних гілках.

Інтегральне зображення[ред. | ред. код]

Останнє співвідношення говорить нам, що для фіксованого значення s, γ є первісною голоморфної функції zs−1 ez. Отже, [11], для будь-яких комплексних u, v ≠ 0,

має місце до тих пір, поки шлях інтегрування, цілком лежить в області підінтегральної гілки. Якщо додатково, дійсна частина s є додатною, то застосовується межа γ(s, u) → 0 при u → 0, остаточно приходячи до комплексного інтегрального визначення γ

[12]

Будь-який шлях інтегрування, який містить 0 тільки в його початку, а в іншому випадку обмежується областю підінтегральної гілки, є справедливим, наприклад, пряма лінія, що з'єднує 0 і z.

Дійсні значення[ред. | ред. код]

Враховуючи інтегральне подання головної гілки γ, для будь-якого позитивного дійсного s, x виконується таке рівняння:[13]

Комплексна s[ред. | ред. код]

Цей результат поширюється на комлексну s. Припустимо спочатку 1 ≤ Re(s) ≤ 2 nf 1 < a < b. Потім

де

[14]

був використаний посередині. Оскільки кінцевий інтеграл стає довільно малим, тільки якщо a досить великий, γ(s, x) рівномірно сходиться для x → ∞ на проміжку 1 ≤ Re(s) ≤ 2 відношенню до голоморфної функції, яка повинна дорівнювати Γ(s) за теоремою про тотожність [15]. Беручи межу в рекурентному співвідношенні γ(s,x) = (s − 1)γ(s − 1,x) − xs −1 ex, і відзначаючи, що lim xn ex = 0 при x → ∞ і всіх n, бачимо, що γ(s,x) також сходиться за межею проміжку до функції, яка відповідає рекурентному відношенню Γ-функції. З цього слідує

для всіх комплексних s недодатніх цілих чисел, x дійсне та γ основна.

Порожня збіжність[ред. | ред. код]

Нехай u з проміжку |arg z| < δ < π/2 з деякими фіксованими δ (α = 0), γ є основною галуззю в цьому секторі, бачимо, що

Як видно вище, перша різниця може бути довільно мала, якщо |u| є достатньо великим. Друга різниця дозволяє зробити наступну оцінку

де ми використали інтегральне зображення γ та формулу про |zs|. Якщо ми інтегруємо вздовж дуги з радіусом R = |u| навколо 0, що зв'язує u і |u|, тоді останній інтеграл дорівнює

де M = δ(cos δ)−Re s eIm є постійною незалежною u або R. Знову посилаючись на поведінку xn ex для великого x, ми бачимо, що останній вираз наближається до 0, оскільки R збільшується до ∞. Отже,маємо:

якщо s не є невід'ємним цілим числом, 0 < ε < π/2 є довільно малим, але фіксованим, а γ позначає основну гілку в цій області.

Загальне зображення[ред. | ред. код]

is:

Верхня неповна гамма-функція[ред. | ред. код]

Що стосується Верхньої неповної гамма-функції, голоморфне розширення, по відношенню до z або s, задається

[16]

в точках (s, z), де права частина існує. Оскільки багатозначна, те саме для , але обмеження до основних значень дає лише однозначну основну гілку .

Коли s є не додатнім цілим числом у вищезгаданому рівнянні, то ні одна частина різниці не визначена, а обмеження процесу, розвинене для s → 0, заповнює відсутні значення. Комплексний аналіз гарантує голоморфність, оскільки виявляється обмежена в сусіди цього обмеження для фіксованого z[17].

Щоб визначити межу, потужний ряд при z = 0 виявляється корисним. Замінивши своїм послідовним рядком у інтегральному визначенні , одержуємо (припускаємо x,s додатні дійсні числа в даний момент):

або

[18]

який як серійне представлення всієї функції , збігається для всього комплексу x (а весь комплекс s не є не додатнім цілим числом).

З його обмеженням до справжніх значень, серія дозволяє робити розширення:

Де s → 0:

,

( Стала Ейлера—Маскероні ), отже,

є граничною функцією верхньої неповної гамма-функції як s → 0, також відомою як .[2]

За допомогою рекурентного співвідношення, з цього результату можна вивести значення для натуральних чисел n[3]

тому верхня неповна гамма-функція виявляється існуючою і голоморфною, як для z , так і для s для всіх s і z ≠ 0.

is:

  • Ціла функція в z для фіксованого додатного інтеграла s;
  • = для s з додатної дійсною частиною z = 0 (границя, коли ),

але це безперервне продовження, а не аналітичне (не тримається для дійсного s <0!);

Спеціальні значення[ред. | ред. код]

  • якщо s додатне ціле,
  • якщо s додатне ціле,[4]
  • ,
  • ,
  • ,
  • for ,
  • ,
  • ,
  • .

Тут, це інтегральна показникова функція, це узагальнена інтегрально-показникова функція, це функція помилок, і це доповнююча функція помилок, .

Асимптотична поведінка[ред. | ред. код]

  • якщо ,
  • якщо and (для дійсних s, помилка Γ(s, x) ~ −xs / s знаходиться на прядку O(xmin{s + 1, 0}) if s ≠ −1 та O(ln(x)) якщо s = −1),
  • коли ,
  • коли ,
  • як асимптотичний розклад при та .[5]

Оціночні формули[ред. | ред. код]

Нижню неповну гамма-функцію можна оцінити за допомогою розширення потужності рядка:[19]

Альтернативне розширення - це

де M це вироджена гіпергеометрична функція Куммера.

Зв'язок з виродженою гіпергеометричною функцією Куммера[ред. | ред. код]

Коли дійсна частина z додатна,

де

має нескінченний радіус збіжності.

Знову ж таки, з виродженими гіпергеометричними функціями і використанням тотожності Куммера,,

Для фактичного обчислення числових значень ланцюговий дріб Гауса забезпечує корисне розширення:

Цей неперервний дріб збігається для всіх комплексних z, тоді і тільки тоді, коли s є додатнім цілим числом.

Верхня неповна гамма-функція має неперервний дріб

[6]

та

[джерело?]

Теорема множення[ред. | ред. код]

Існує наступна теорема множення:

Регуляризовані гамма-функції та випадкові величини Пуассона[ред. | ред. код]

Дві пов'язані функції — це регуляризовані гамма-функції:

— це функції розподілу ймовірностей для випадкових змінних гамма-функції з параметром форми та коефіцієнтом масштабу 1.

Коли ціле, є сукупною функцією розподілу для змінних Пуассона: Якщо є випадковою величиною, тоді

Ця формула може бути отримана шляхом повторного інтегрування частинами.

Похідні[ред. | ред. код]

Похідна верхньої неповної гамма-функції по відношенню до x добре відома. Вона задається від'ємним підсумком від його інтегрального визначення (від оцінюваного на нижній межі):

Похідна по відношенню до її першого аргументу задана як[7]

і друга похідна як

де функція є особливим випадком G-функції Мейєра

Цей конкретний окремий випадок має власні властивості внутрішнього закриття, тому він може використовуватися для виразу всіх послідовних похідних. В загальному,

де це перестановка, визначена символом Похгамера:

Усі такі похідні можуть бути сформовані послідовно з:

та

Ця функція може бути обчислена з представлення її серії, дійсної для ,

з умовою, що s - це не від'ємне ціле число чи нуль. У такому випадку треба використовувати обмеження. Результати для можна отримати за допомогою аналітичного продовження. Деякі особливі випадки цієї функції можуть бути спрощені. Наприклад, , , де це Інтегральна показникова функція. Ці похідні та функція дають точні розв'язання для ряду інтегралів шляхом повторного диференціювання інтегрального визначення верхньої неповної гамма-функції. До прикладу,

Ця формула може бути далі наповнена або узагальнена на величезний клас перетворення Лапласа та перетворення Мелліна. У поєднанні з системою комп'ютерної алгебри експлуатація спеціальних функцій забезпечує потужний метод вирішення певних інтегралів, зокрема тих, що зустрічаються у практичних інженерних застосуваннях.

Невизначені і визначені інтеграли[ред. | ред. код]

Наступні невизначені інтеграли легко отримати при використанні інтегрування частинами:

Нижня і верхня неповна гамма-функція поєднуються через перетворення Фур'є:

Це випливає, наприклад, з відповідною спеціалізацією від (Gradshteyn та Ryzhik, 2015, §7.642)

Примітки[ред. | ред. код]

  1. DLMF, Incomplete Gamma functions, analytic continuation
  2. DLMF, Incomplete Gamma functions, Special Values, 8.4.4
  3. DLMF, Incomplete Gamma functions, Special Values, 8.4.15
  4. Weisstein, Eric W. Incomplete Gamma Function(англ.) на сайті Wolfram MathWorld. (equation 2)
  5. DLMF, Incomplete Gamma functions, 8.11(i)
  6. Abramowitz and Stegun p. 263, 6.5.31
  7. [K.O. Geddes, M.L. Glasser, R.A. Moore and T.C. Scott, Evaluation of Classes of Definite Integrals Involving Elementary Functions via Differentiation of Special Functions, AAECC (Applicable Algebra in Engineering, Communication and Computing), vol. 1, (1990), pp. 149–165, [1][недоступне посилання]

Див. також[ред. | ред. код]

Джерела[ред. | ред. код]

  • people.math.sfu.ca Incomplete Gamma function. §6.5.
  • Allasia, Giampietro; Besenghi, Renata (1986). Numerical calculation of incomplete gamma functions by the trapezoidal rule. Numer. Math. 50 (4): 419—428. doi:10.1007/BF0139666.
  • Amore, Paolo (2005). Asymptotic and exact series representations for the incomplete Gamma function. Europhys. Lett. 71 (1): 1—7. doi:10.1209/epl/i2005-10066-6. MR 2170316.
  • G. Arfken and H. Weber. Mathematical Methods for Physicists. Harcourt/Academic Press, 2000. (See Chapter 10.)
  • DiDonato, Armido R.; Morris, Jr., Alfred H. (Dec 1986). Computation of the incomplete gamma function ratios and their inverse. ACM Transactions on Mathematical Software. 12 (4): 377—393. doi:10.1145/22721.23109.
  • Barakat, Richard (1961). Evaluation of the Incomplete Gamma Function of Imaginary Argument by Chebyshev Polynomials. Math. Comp. 15 (73): 7—11. doi:10.1090/s0025-5718-1961-0128058-1. MR 0128058.
  • Carsky, Petr; Polasek, Martin (1998). Incomplete Gamma F_m(x) functions for real and complex arguments. J. Comput. Phys. 143 (1): 259—265. doi:10.1006/jcph.1998.5975. MR 1624704.
  • Chaudhry, M. Aslam; Zubair, S. M. (1995). On the decomposition of generalized incomplete Gamma functions with applications to Fourier transforms. J. Comput. Appl. Math. 59 (101): 253—284. doi:10.1016/0377-0427(94)00026-w. MR 1346414.
  • DiDonato, Armido R.; Morris, Jr., Alfred H. (Sep 1987). ALGORITHM 654: FORTRAN subroutines for computing the incomplete gamma function ratios and their inverse. ACM Transactions on Mathematical Software. 13 (3): 318—319. doi:10.1145/29380.214348. (See also www.netlib.org/toms/654).
  • Früchtl, H.; Otto, P. (1994). A new algorithm for the evaluation of the incomplete Gamma Function on vector computers. ACM Trans. Math. Softw. 20 (4): 436—446. doi:10.1145/198429.198432.
  • Gautschi, Walter (1998). The incomplete gamma function since Tricomi. Atti Convegni Lincei. 147: 203—237. MR 1737497.
  • Gautschi, Walter (1999). A Note on the recursive calculation of Incomplete Gamma Functions. ACM Trans. Math. Softw. 25 (1): 101—107. doi:10.1145/305658.305717. MR 1697463.
  • Gradshteyn, Izrail Solomonovich; Ryzhik, Iosif Moiseevich; Geronimus, Yuri Veniaminovich; Tseytlin author-first5=Alan, Michail Yulyevich; Jeffrey (2015). 8.35.. У Zwillinger, Daniel; Moll, Victor Hugo (ред.). Table of Integrals, Series, and Products (English) (вид. 8). [Academic Press, Inc.] с. 908—911. ISBN 978-0-12-384933-5. LCCN 2014010276.
  • Jones, William B.; Thron, W. J. (1985). On the computation of incomplete gamma functions in the complex domain. J. Comp. Appl. Math. 12—13: 401—417. doi:10.1016/0377-0427(85)90034-2. MR 0793971.
  • Mathar, Richard J. (2004). Numerical representation of the incomplete gamma function of complex-valued argument. Numerical Algorithms. 36 (3): 247—264. doi:10.1023/B:NUMA.0000040063.91709.5. MR 2091195.
  • Miller, Allen R.; Moskowitz, Ira S. (1998). On certain Generalized incomplete Gamma functions. J. Comput. Appl. Math. 91 (2): 179—190. doi:10.1016/s0377-0427(98)00031-4.
  • Paris, R. B. (2002). A uniform asymptotic expansion for the incomplete gamma function. J. Comput. Appl. Math. 148 (2): 323—339. doi:10.1016/S0377-0427(02)00553-8. MR 1936142.
  • Press, WH; Teukolsky, SA; Vetterling, WT; Flannery, BP (2007). Section 6.2. Incomplete Gamma Function and Error Function. Numerical Recipes: The Art of Scientific Computing (вид. 3). New York: Cambridge University Press. ISBN 978-0-521-88068-8.
  • Takenaga, Roy (1966). On the Evaluation of the Incomplete Gamma Function. Math. Comp. 20 (96): 606—610. doi:10.1090/S0025-5718-1966-0203911-3. MR 0203911.
  • Temme, Nico (1975). Uniform Asymptotic Expansions of the Incomplete Gamma Functions and the Incomplete Beta Function. Math. Comp. 29 (132): 1109—1114. doi:10.1090/S0025-5718-1975-0387674-2. MR 0387674.
  • Terras, Riho (1979). The determination of incomplete Gamma Functions through analytic integration. J. Comp. Phys. 31: 146—151. doi:10.1016/0021-9991(79)90066-4. MR 0531128.
  • Tricomi, Francesco G. (1950). Sulla funzione gamma incompleta. Ann. Mat. Pura Appl. 31: 263—279. doi:10.1007/BF02428264. MR 0047834.
  • Tricomi, F. G. (1950). Asymptotische Eigenschaften der unvollst. Gammafunktion. Math. Zeitsch. 53 (2): 136—148. doi:10.1007/bf01162409. MR 0045253.
  • van Deun, Joris; Cools, Ronald (2006). A stable recurrence for the incomplete gamma function with imaginary second argument. Numer. Math. 104: 445—456. doi:10.1007/s00211-006-0026-1. MR 2249673.
  • Winitzki, Serge (2003). Computing the incomplete gamma function to arbitrary precision. Lect. Not. Comp. Sci. 2667: 790—798. doi:10.1007/3-540-44839-x_83. MR 2110953.
  • Weisstein, Eric W. Incomplete Gamma Function(англ.) на сайті Wolfram MathWorld.