Нижня неповна гамма-функція в 3DАнімація нижньої неповної гамма-функціїВерхня неповна гамма-функція в 3D
У математиціверхня неповна гамма-функція і нижня неповна гамма-функція є типом спеціальних функцій, які виникають при розв'язанні різноманітних математичних задач, таких як деякі інтеграли.
Їх відповідні назви випливають з їх інтегральних визначень, які визначаються аналогічно гамма-функції, іншим типом спеціальної функції, але з різними або "неповними" інтегральними межами. Гамма-функція визначається як інтеграл від нуля до нескінченності. Це відрізняється від нижньої неповної гамма-функції, яка визначається як інтеграл від нуля до змінної верхньої межі. Відповідно, верхня неповна гамма-функція визначається як інтеграл від змінної нижньої межі до нескінченності.
Визначення
Верхня неповна гамма-функція визначається як:
в той час як нижня неповна гамма-функція визначається як:
Властивості
В обох випадках s є складним параметром, таким, що дійсна частина s є позитивною.
Нижня неповна гамма-функція та верхня неповна гамма-функція, що визначені вище для дійсних позитивних s і x, можуть бути розвинені в голоморфній функції s, по відношенню як до x, так і до s, визначені для майже всіх комбінацій складних x та s.[1] Комплексний аналіз показує, що властивості дійсних неповних гамма-функцій поширюються на їх голоморфні аналоги.
Голоморфне розширення
Повторне застосування рекурентного відношення до нижньої неповної гамма-функції призводить до розширення степеневого ряду:
[2]
Враховуючи стрімке зростання абсолютної величини Γ(z + k) при k → ∞, а також той факт, що взаємна Γ( z ) є цілою функцією, коефіцієнти у правій частині суми є чітко визначеними, а локально сума сходиться рівномірно для всіх комплексних s та x. За теоремою Вейерштрасса, гранична функція іноді позначається як ,
розширює дійсну неповну гамма-функцію в голоморфну функцію, як спільно, так і окремо в z і s. З властивостей zs та гамма-функції випливає, що перші два чинники фіксують особливі точки γ (при z = 0 або s - не додатнє ціле число), тоді як останній чинник сприяє його нулям.
Багатозначність
Комплексний логарифм log z = log |z| + i arg z визначається лише кратним 2πi, що робить його багатозначним. Функції, що включають складний логарифм, як правило, успадковують цю властивість. Серед них - складна потужність, і оскільки zs з'являється в його розкладі, то γ-функція теж.
Невизначеність багатозначних функцій призводить до ускладнень, оскільки необхідно вказати, як вибрати значення. Стратегії для вирішення цього є:
(найбільш загальний спосіб) замінити домен ℂ багатозначних функцій відповідним різновидом у ℂ × ℂ, що називається поверхнею Рімана. У той час як це усуває багатозначність, потрібно знати теорію, яка лежить в основі [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−1e−z[10] зберігаються на відповідних гілках.
Інтегральне зображення
Останнє співвідношення говорить нам, що для фіксованого значення s, γ є первісною голоморфної функції zs−1e−z. Отже, [11], для будь-яких комплексних u, v ≠ 0,
має місце до тих пір, поки шлях інтегрування, цілком лежить в області підінтегральної гілки. Якщо додатково, дійсна частина s є додатною, то застосовується межа γ(s, u) → 0 при u → 0, остаточно приходячи до комплексного інтегрального визначення γ
Будь-який шлях інтегрування, який містить 0 тільки в його початку, а в іншому випадку обмежується областю підінтегральної гілки, є справедливим, наприклад, пряма лінія, що з'єднує 0 і z.
Дійсні значення
Враховуючи інтегральне подання головної гілки γ, для будь-якого позитивного дійсного s, x виконується таке рівняння:[13]
Комплексна s
Цей результат поширюється на комлексну s. Припустимо спочатку 1 ≤ Re(s) ≤ 2 nf 1 < a < b. Потім
був використаний посередині. Оскільки кінцевий інтеграл стає довільно малим, тільки якщо a досить великий, γ(s, x) рівномірно сходиться для x → ∞ на проміжку 1 ≤ Re(s) ≤ 2 відношенню до голоморфної функції, яка повинна дорівнювати Γ(s) за теоремою про тотожність [15]. Беручи межу в рекурентному співвідношенні γ(s,x) = (s − 1)γ(s − 1,x) − xs
−1e−x, і відзначаючи, що lim xne−x = 0 при x → ∞ і всіх n, бачимо, що γ(s,x) також сходиться за межею проміжку до функції, яка відповідає рекурентному відношенню Γ-функції. З цього слідує
для всіх комплексних s недодатніх цілих чисел, x дійсне та γ основна.
Порожня збіжність
Нехай u з проміжку |arg z| < δ < π/2 з деякими фіксованими δ (α = 0), γ є основною галуззю в цьому секторі, бачимо, що
Як видно вище, перша різниця може бути довільно мала, якщо |u| є достатньо великим. Друга різниця дозволяє зробити наступну оцінку
де ми використали інтегральне зображення γ та формулу про |zs|. Якщо ми інтегруємо вздовж дуги з радіусом R = |u| навколо 0, що зв'язує u і |u|, тоді останній інтеграл дорівнює
де M = δ(cos δ)−Re seIm sδ є постійною незалежною u або R. Знову посилаючись на поведінку xne−x для великого x, ми бачимо, що останній вираз наближається до 0, оскільки R збільшується до ∞.
Отже,маємо:
якщо s не є невід'ємним цілим числом, 0 < ε < π/2 є довільно малим, але фіксованим, а γ позначає основну гілку в цій області.
Загальне зображення
is:
Ціла функція в z для фіксованого додатного інтеграла s;
в точках (s, z), де права частина існує. Оскільки багатозначна, те саме для , але обмеження до основних значень дає лише однозначну основну гілку .
Коли s є не додатнім цілим числом у вищезгаданому рівнянні, то ні одна частина різниці не визначена, а обмеження процесу, розвинене для s → 0, заповнює відсутні значення. Комплексний аналіз гарантує голоморфність, оскільки виявляється обмежена в сусіди цього обмеження для фіксованого z[17].
Щоб визначити межу, потужний ряд при z = 0 виявляється корисним. Замінивши своїм послідовним рядком у інтегральному визначенні , одержуємо (припускаємо x,s додатні дійсні числа в даний момент):
Коли ціле, є сукупною функцією розподілу для змінних Пуассона: Якщо є випадковою величиною, тоді
Ця формула може бути отримана шляхом повторного інтегрування частинами.
Похідні
Похідна верхньої неповної гамма-функції по відношенню до x добре відома. Вона задається від'ємним підсумком від його інтегрального визначення (від оцінюваного на нижній межі):
Похідна по відношенню до її першого аргументу задана як[7]
Цей конкретний окремий випадок має власні властивості внутрішнього закриття, тому він може використовуватися для виразу всіх послідовних похідних. В загальному,
Усі такі похідні можуть бути сформовані послідовно з:
та
Ця функція може бути обчислена з представлення її серії, дійсної для ,
з умовою, що s - це не від'ємне ціле число чи нуль. У такому випадку треба використовувати обмеження. Результати для можна отримати за допомогою аналітичного продовження. Деякі особливі випадки цієї функції можуть бути спрощені. Наприклад, , , де це Інтегральна показникова функція. Ці похідні та функція дають точні розв'язання для ряду інтегралів шляхом повторного диференціювання інтегрального визначення верхньої неповної гамма-функції.
До прикладу,
Ця формула може бути далі наповнена або узагальнена на величезний клас перетворення Лапласа та перетворення Мелліна. У поєднанні з системою комп'ютерної алгебри експлуатація спеціальних функцій забезпечує потужний метод вирішення певних інтегралів, зокрема тих, що зустрічаються у практичних інженерних застосуваннях.
↑[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][недоступне посилання]
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. MR2170316.
G. Arfken and H. Weber. Mathematical Methods for Physicists. Harcourt/Academic Press, 2000. (See Chapter 10.)
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. MR1624704.
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. MR1346414.
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. MR1737497.
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. MR1697463.
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. ISBN978-0-12-384933-5. LCCN2014010276. {{cite book}}: Пропущено вертикальну риску в: |author-last4= (довідка)
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. MR0793971.
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. MR2091195.
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. MR1936142.
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. MR0531128.
Tricomi, Francesco G. (1950). Sulla funzione gamma incompleta. Ann. Mat. Pura Appl. 31: 263—279. doi:10.1007/BF02428264. MR0047834.
Tricomi, F. G. (1950). Asymptotische Eigenschaften der unvollst. Gammafunktion. Math. Zeitsch. 53 (2): 136—148. doi:10.1007/bf01162409. MR0045253.
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. MR2249673.
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. MR2110953.