Wykres funkcji G Meigera G(((a 1,...,an),(a n+1,...,ap)),((b 1,...,bm),(b m+ 1,...,bq)),z) wejście ((½),()),((⅓),()) na płaszczyźnie zespolonej od -2-2i do 2+2i
Dzięki nowoczesnej definicji większość ustalonych funkcji specjalnych można przedstawić za pomocą funkcji G Meijera. Godną uwagi właściwością jest zamknięcie zbioru wszystkich funkcji G nie tylko przy różniczkowaniu, ale także przy całkowaniu nieokreślonym. W połączeniu z równaniem funkcyjnym , które pozwala wyzwolić z funkcji G G ( z ) dowolny czynnik z ρ będący stałą potęgą jego argumentu z , domknięcie implikuje, że ilekroć funkcja jest wyrażalna jako funkcja G stałej wielokrotności pewnej stałej potęgi argumentu funkcji, f ( x ) = G ( cx γ ), pochodna i funkcja pierwotna tej funkcji są wyrażalne, więc zbyt.
Szeroki zakres funkcji specjalnych nadaje również moc zastosowaniom funkcji G Meijera innym niż reprezentacja i manipulowanie pochodnymi i funkcjami pierwotnymi. Na przykład całka oznaczona po dodatniej osi rzeczywistej dowolnej funkcji g ( x ), którą można zapisać jako iloczyn G 1 ( cx γ )· G 2 ( dx δ ) dwóch funkcji G o wartościach wymiernych γ / δ równa się po prostu innej funkcji G, a uogólnienia transformacji całkowych , takich jak transformata Hankla i transformata Laplace'a oraz ich odwrotności, są wynikiem zastosowania odpowiednich par funkcji G jako jąder transformacji.
Jeszcze bardziej ogólną funkcją, która wprowadza dodatkowe parametry do funkcji G Meijera, jest funkcja H Foxa i jest używana do transformacji Matrix przez Ram Kishore Saxena
Jednym z zastosowań funkcji G Meijera było widmo cząstek promieniowania z inercyjnego horyzontu w modelu ruchomego lustra dynamicznego efektu Casimira ( Good 2020 ).
0 ≤ m ≤ q i 0 ≤ n ≤ p , gdzie m , n , p i q są liczbami całkowitymi
a k − b j ≠ 1, 2, 3, ... dla k = 1, 2, ..., n i j = 1, 2, ..., m , co implikuje, że żaden biegun dowolnego Γ( b j − s ), j = 1, 2, ..., m , pokrywa się z dowolnym biegunem dowolnego Γ(1 − a k + s ), k = 1, 2, ..., n
z ≠ 0
Należy zauważyć, że ze względów historycznych pierwszy dolny i drugi górny indeks odnoszą się do górnego rzędu parametrów, podczas gdy drugi dolny i pierwszy górny indeks odnoszą się do dolnego rzędu parametrów. Często spotyka się następującą, bardziej syntetyczną notację wykorzystującą wektory :
Implementacje funkcji G w systemach algebry komputerowej zazwyczaj wykorzystują oddzielne argumenty wektorowe dla czterech (prawdopodobnie pustych) grup parametrów a 1 ... a n , a n +1 ... a p , b 1 ... b m , i b m +1 ... b q , a zatem można pominąć rzędy p , q , n i m jako zbędne.
L w całce reprezentuje ścieżkę, którą należy podążać podczas całkowania . Dla tej ścieżki możliwe są trzy opcje:
1. L biegnie od − i ∞ do + i ∞ tak, że wszystkie bieguny Γ( b j − s ), j = 1, 2, ..., m , są po prawej stronie ścieżki, podczas gdy wszystkie bieguny Γ (1 − za k + s ), k = 1, 2, ..., n , są po lewej stronie. Następnie całka jest zbieżna dla |arg z | < δ π , gdzie
oczywistym warunkiem wstępnym jest δ > 0. Całka dodatkowo zbiega się dla | arg z | = δ π ≥ 0 jeśli (q - p) ( σ + 1 ⁄ 2 ) > Re ( ν ) + 1, gdzie σ reprezentuje Re ( s ), gdy zmienna całkowa s zbliża się zarówno do + i ∞, jak i - i ∞, i gdzie
arg z | = δ π i p = q całka jest zbieżna niezależnie od σ zawsze, gdy Re( ν ) < −1.
2. L jest pętlą zaczynającą się i kończącą w +∞, otaczającą wszystkie bieguny Γ( b j − s ), j = 1, 2, ..., m , dokładnie raz w kierunku ujemnym, ale nie otaczającym żadnego bieguna Γ(1 − a k + s ), k = 1, 2, ..., n . Wtedy całka jest zbieżna dla wszystkich z , jeśli q > p ≥ 0; jest również zbieżny dla q = p > 0, o ile | z | < 1. W tym drugim przypadku całka dodatkowo jest zbieżna dla | z | = 1 jeśli Re( ν ) < −1, gdzie ν jest zdefiniowane jak dla pierwszej ścieżki.
3. L jest pętlą zaczynającą się i kończącą w −∞ i okrążającą wszystkie bieguny Γ(1 − a k + s ), k = 1, 2, ..., n , dokładnie raz w kierunku dodatnim, ale nie okrążającą żadnej biegun Γ( b jot - s ), jot = 1, 2, ..., m . Teraz całka jest zbieżna dla wszystkich z , jeśli p > q ≥ 0; jest również zbieżny dla p = q > 0 dopóki | z | > 1. Jak zauważono również dla drugiej ścieżki, w przypadku p = q całka jest również zbieżna dla | z | = 1 gdy Re( ν ) < −1.
Warunki zbieżności można łatwo ustalić, stosując asymptotyczne przybliżenie Stirlinga do funkcji gamma w całce. Kiedy całka jest zbieżna dla więcej niż jednej z tych ścieżek, można wykazać, że wyniki całkowania są zgodne; jeśli zbiega się tylko dla jednej ścieżki, to jest to jedyna, którą należy wziąć pod uwagę. W rzeczywistości numeryczne całkowanie ścieżki na płaszczyźnie zespolonej stanowi praktyczne i rozsądne podejście do obliczania funkcji G Meijera.
W konsekwencji tej definicji funkcja G Meijera jest funkcją analityczną z z możliwym wyjątkiem pochodzenia z = 0 i okręgu jednostkowego | z | = 1.
Za fundamentalny zbiór rozwiązań tego równania w przypadku p ≤ q można przyjąć:
i podobnie w przypadku p ≥ q :
Te konkretne rozwiązania są analityczne, z wyjątkiem możliwej osobliwości w z = 0 (jak również możliwej osobliwości w z = ∞), aw przypadku p = q także nieuniknionej osobliwości w z = (-1) p - m - rz . Jak zobaczymy obecnie, można je utożsamiać z uogólnionymi funkcjami hipergeometrycznymi p F q −1 argumentu (−1) p − m − n z które są mnożone przez potęgę z b h , oraz z uogólnionymi funkcjami hipergeometrycznymi q F p −1 argumentu (−1) q − m − n z −1 które są mnożone odpowiednio przez potęgę z a h −1 .
Związek między funkcją G a uogólnioną funkcją hipergeometryczną
Jeśli całka jest zbieżna, gdy jest oceniana wzdłuż drugiej ścieżki wprowadzonej powyżej, i jeśli nie ma zlewających się biegunów między Γ ( b j - s ), j = 1, 2, ..., m , to można wyrazić funkcję G Meijera jako suma reszt pod względem uogólnionych funkcji hipergeometrycznych p F q −1 (twierdzenie Slatera):
Gwiazdka wskazuje, że termin odpowiadający j = h został pominięty. Aby całka była zbieżna wzdłuż drugiej ścieżki, trzeba mieć albo p < q , albo p = q i | z | < 1 i aby bieguny były różne, żadna para między b j , j = 1, 2, ..., m , nie może różnić się o liczbę całkowitą lub zero. Gwiazdki w relacji przypominają nam o zignorowaniu wkładu o indeksie j = h następująco: W iloczynie sprowadza się to do zamiany Γ(0) na 1, aw argumencie funkcji hipergeometrycznej, jeśli przypomnimy sobie znaczenie zapisu wektorowego,
sprowadza się to do skrócenia długości wektora z q do q -1.
Zauważ, że gdy m = 0, druga ścieżka nie zawiera żadnego bieguna, więc całka musi zniknąć identycznie,
jeśli albo p < q , albo p = q i | z | < 1.
Podobnie, jeśli całka jest zbieżna, gdy jest oceniana wzdłuż trzeciej ścieżki powyżej, i jeśli między Γ (1 - a k + s ), k = 1, 2, ..., n nie pojawiają się zlewające się bieguny , k = 1, 2, ..., n , to funkcja G może wyrazić jako:
W tym celu albo p > q , albo p = q i | z | > 1 są wymagane i żadna para spośród a k , k = 1, 2, ..., n , nie może różnić się o liczbę całkowitą lub zero. Dla n = 0 mamy zatem:
jeśli albo p > q , albo p = q i | z | > 1.
Z drugiej strony każdą uogólnioną funkcję hipergeometryczną można łatwo wyrazić za pomocą funkcji G Meijera:
gdzie skorzystaliśmy z notacji wektorowej:
Dzieje się tak, chyba że niedodatnia liczba całkowita co najmniej jednego z jej parametrów a p redukuje funkcję hipergeometryczną do skończonego wielomianu, w którym to przypadku prefaktor gamma jednej z funkcji G znika, a zestawy parametrów funkcji G naruszają wymaganie a k − b j ≠ 1, 2, 3, ... dla k = 1, 2, ..., n i j = 1, 2, ..., m z definicji powyżej. Poza tym ograniczeniem, zależność obowiązuje zawsze wtedy, gdy uogólniony szereg hipergeometryczny p F q ( z ) jest zbieżny, tj. dla dowolnego skończonego z , gdy p ≤ q , oraz dla | z | < 1 gdy p = q + 1. W tym drugim przypadku związek z funkcją G automatycznie zapewnia analityczną kontynuację p F q ( z ) do | z | ≥ 1 z gałęzią przeciętą od 1 do ∞ wzdłuż osi rzeczywistej. Wreszcie, relacja dostarcza naturalnego rozszerzenia definicji funkcji hipergeometrycznej na rzędy p > q + 1. Za pomocą funkcji G możemy zatem rozwiązać uogólnione hipergeometryczne równanie różniczkowe również dla p > q + 1.
Przypadki wielomianowe
Aby wyrazić wielomianowe przypadki uogólnionych funkcji hipergeometrycznych za pomocą funkcji G Meijera, ogólnie potrzebna jest liniowa kombinacja dwóch funkcji G:
gdzie h = 0, 1, 2, ... jest równe stopniowi wielomianu p +1 F q ( z ). Rzędy m i n można wybierać dowolnie w przedziałach 0 ≤ m ≤ q i 0 ≤ n ≤ p , co pozwala uniknąć sytuacji, w której określone wartości całkowite lub różnice całkowite między parametrami a p i b q wielomianu powodują rozbieżności gamma w prefaktorze lub w konflikcie z definicja funkcji G. Zauważ, że pierwsza funkcja G znika dla n = 0, jeśli p > q , podczas gdy druga funkcja G znika dla m = 0, jeśli p < q . Ponownie wzór można zweryfikować, wyrażając dwie funkcje G jako sumy reszt ; nie trzeba tutaj wykluczać żadnych przypadków zlewających się biegunów dozwolonych przez definicję funkcji G.
Podstawowe własności funkcji G
Jak widać z definicji funkcji G , jeśli wśród a p i b q występują równe parametry określające czynniki w liczniku i mianowniku całki, ułamek można uprościć, a tym samym rząd funkcji być zredukowanym. To, czy rząd m czy n zmniejszy się, zależy od konkretnej pozycji danych parametrów. Zatem jeśli jedno z a k , k = 1, 2, ..., n , jest równe jednemu z b j , j = m + 1, ..., q , funkcja G obniża swoje rzędy p , q i n :
Z tego samego powodu, jeśli jedno z a k , k = n + 1, ..., p , jest równe jednemu z b j , j = 1, 2, ..., m , to funkcja G obniża swoją rzędy p , q i m :
Wychodząc z definicji, możliwe jest również wyprowadzenie następujących właściwości:
Jeśli chodzi o pochodne funkcji G, można znaleźć następujące zależności:
Z tych czterech równoważnych relacji można wywnioskować, po prostu oceniając pochodną po lewej stronie i nieco manipulując. Otrzymuje się np.:
Co więcej, dla pochodnych dowolnego rzędu h , mamy
które zachowują się również dla h <0, umożliwiając w ten sposób uzyskanie funkcji pierwotnej dowolnej funkcji G równie łatwo, jak pochodną. Wybierając jeden lub drugi z dwóch wyników podanych w dowolnym wzorze, zawsze można zapobiec naruszeniu przez zestaw parametrów wyniku warunku a k − b j ≠ 1, 2, 3, ... dla k = 1, 2, ..., n i j = 1, 2, ..., m narzucone przez definicję funkcji G . Zauważ, że każda para wyników staje się nierówna w przypadku h < 0.
Przyrównując różne wyrażenia dla pochodnych pierwszego rzędu, dochodzi się do następujących 3-członowych relacji rekurencyjnych między ciągłymi funkcjami G:
Podobne relacje dla diagonalnych par parametrów a 1 , b q i b 1 , a p wynikają z odpowiedniej kombinacji powyższych. Ponownie, z tych relacji powtarzalności można wyprowadzić odpowiednie właściwości funkcji hipergeometrycznych i innych funkcji specjalnych.
Twierdzenia o mnożeniu
Zakładając, że z ≠ 0, zachodzą następujące zależności:
Należy zauważyć, że ograniczenia, w ramach których istnieje ta całka, zostały tutaj pominięte. Oczywiście nie jest niespodzianką, że transformata Mellina funkcji G powinna prowadzić z powrotem do całki występującej w powyższej definicji .
Obszerne ograniczenia, w ramach których istnieją te całki, można znaleźć na str. 417 „Tabele przekształceń całkowych”, tom. II (1954), pod redakcją A. Erdelyi. Należy zauważyć, że ze względu na ich wpływ na funkcję G, całki te można wykorzystać do zdefiniowania operacji całkowania ułamkowego dla dość dużej klasy funkcji ( operatory Erdélyi – Kobera ).
Wynikiem o fundamentalnym znaczeniu jest to, że iloczyn dwóch dowolnych funkcji G zintegrowanych na dodatniej osi rzeczywistej może być reprezentowany przez inną funkcję G (twierdzenie o splocie):
Ograniczenia, w ramach których istnieje całka, można znaleźć w Meijer, CS, 1941: Nederl. Akad. Wetensch, Proc. 44, s. 82–92. Zwróć uwagę, jak transformata Mellina wyniku po prostu łączy czynniki gamma z transformacji Mellina dwóch funkcji w całce.
Formułę splotu można wyprowadzić, zastępując definiującą całkę Mellina-Barnesa jedną z funkcji G, odwracając kolejność całkowania i oceniając wewnętrzną całkę z transformacją Mellina. Powyższe całki typu Eulera następują analogicznie.
gdzie Re( ω ) > 0. To jest transformata Laplace'a funkcji G ( ηx ) pomnożona przez potęgę x − α ; jeśli umieścimy α = 0, otrzymamy transformatę Laplace'a funkcji G. Jak zwykle, odwrotna transformata jest następnie dana przez:
gdzie c jest rzeczywistą stałą dodatnią, która umieszcza ścieżkę całkowania na prawo od dowolnego bieguna w całce.
Inny wzór na transformatę Laplace'a funkcji G to:
gdzie ponownie Re( ω ) > 0. W obu przypadkach pominięto szczegóły ograniczeń, w jakich istnieją całki.
Transformacje całkowe oparte na funkcji G
Ogólnie rzecz biorąc, dwie funkcje k ( z , y ) i h ( z , y ) nazywane są parą jąder transformacji, jeśli dla dowolnej odpowiedniej funkcji f ( z ) lub dowolnej odpowiedniej funkcji g ( z ) zachodzą jednocześnie następujące dwie zależności :
Mówimy, że para jąder jest symetryczna, jeśli k ( z , y ) = h ( z , y ).
są asymetryczną parą jąder transformacji, gdzie γ > 0, n - p = m - q > 0, oraz:
wraz z dalszymi warunkami konwergencji. W szczególności, jeśli p = q , m = n , a j + b j = 0 dla j = 1, 2, ..., p i c j + re j = 0 dla j = 1, 2, ..., m , to para jąder staje się symetryczna. Dobrze znana transformata Hankla jest symetrycznym przypadkiem szczególnym transformaty Naraina ( γ = 1, s = q = 0, m = n = 1, do 1 = - re 1 = ν ⁄ 2 ).
Transformacja mięczaka
Jet Wimp ( 1964 ) pokazał, że te funkcje są asymetryczną parą jąder transformacji:
gdzie funkcja A (·) jest zdefiniowana jako:
Uogólniona transformata Laplace'a
Transformatę Laplace'a można uogólnić w ścisłej analogii z uogólnieniem transformaty Hankla przez Naraina:
0 gdzie γ > , p ≤ q , oraz:
i gdzie stała c > 0 umieszcza drugą ścieżkę całkowania na prawo od dowolnego bieguna w całce. Dla γ = 1 / 2 , ρ = 0 i p = q = 0 odpowiada to znanej transformacie Laplace'a.
Transformacja Meijera
Dwa szczególne przypadki tego uogólnienia podał CS Meijer w latach 1940 i 1941. Przypadek wynikający dla γ = 1, ρ = − ν , p = 0, q = 1 i b 1 = ν można zapisać (Meijer 1940 ):
a przypadek uzyskany dla γ = 1 ⁄ 2 , ρ = − m − k , p = q = 1, a 1 = m − k i b 1 = 2 m można zapisać (Meijer 1941a ):
Łukasz, Yudell L. (1969). Funkcje specjalne i ich przybliżenia, tom. ja . Nowy Jork: prasa akademicka. ISBN 978-0-12-459901-7 . (patrz Rozdział V, „Uogólniona funkcja hipergeometryczna i funkcja G”, s. 136)
Meijera, CS (1936). "Über Whittakersche bzw. Besselsche Funktionen und deren Produkte". Nieuw Archief voor Wiskunde (2) (w języku niemieckim). 18 (4): 10–39. JFM 62.0421.02 .
Meijera, CS (1940). „Über eine Erweiterung der Laplace-Transformation - I, II” . Proceedings of the Section of Sciences, Koninklijke Akademie van Wetenschappen (Amsterdam) (w języku niemieckim). 43 : 599-608 i 702-711. JFM 66.0523.01 .
Meijer, CS (1941a). „Eine neue Erweiterung der Laplace-Transformation - I, II” . Proceedings of the Section of Sciences, Koninklijke Akademie van Wetenschappen (Amsterdam) (w języku niemieckim). 44 : 727-737 i 831-839. JFM 67.0396.01 .
Meijer, CS (1941b). „Multiplikationstheorem für die Funktion ". Proceedings of the Section of Sciences, Koninklijke Akademie van Wetenschappen (Amsterdam) (w języku niemieckim). 44 : 1062-1070. JFM 67.1016.01 .
Prudnikow, AP; Mariczew, OI; Bryczkow, Yu. A. (1990). Całki i serie, tom. 3: więcej funkcji specjalnych . Newark, NJ: Gordon i naruszenie. ISBN 978-2-88124-682-1 . (patrz § 8.2, „Funkcja G Meijera”, s. 617)
Mathai, Saxena, AM i RK (1973). Uogólnione funkcje hipergeometryczne z zastosowaniami w statystyce i naukach fizycznych . Springera . ISBN 9780387064826 .