v matematika, Hankelova transformace vyjadřuje jakoukoli danou funkci F(r) jako vážený součet nekonečného počtu Besselovy funkce prvního druhuJν(kr). Besselovy funkce v součtu jsou všechny stejného řádu ν, ale liší se v měřítku k podél r osa. Potřebný koeficient Fν každé Besselovy funkce v součtu, jako funkce měřítka k představuje transformovanou funkci. Hankelova transformace je integrální transformace a byl poprvé vyvinut matematikem Hermann Hankel. Je také známá jako Fourier-Besselova transformace. Stejně jako Fourierova transformace pro nekonečný interval souvisí s Fourierova řada přes konečný interval, takže Hankelova transformace přes nekonečný interval souvisí s Série Fourier – Bessel v konečném intervalu.
The Hankelova transformace řádu funkce F(r) darováno
kde je Besselova funkce prvního druhu objednávky s . Inverzní Hankelova transformace Fν(k) je definován jako
které lze snadno ověřit pomocí níže popsaného vztahu ortogonality.
Doména definice
Invertování Hankelovy transformace funkce F(r) je platný v každém okamžiku, kdy F(r) je spojitý, za předpokladu, že funkce je definována v (0, ∞), je po částech spojitý a má omezenou variabilitu v každém konečném podintervalu v (0, ∞) a
Stejně jako Fourierova transformace však lze doménu rozšířit argumentem hustoty, aby zahrnoval některé funkce, jejichž výše uvedený integrál není konečný, například .
Alternativní definice
Alternativní definice říká, že Hankelova transformace G(r) je[1]
Tyto dvě definice spolu souvisejí:
Li , pak
To znamená, že stejně jako u předchozí definice je takto definovaná Hankelova transformace také svou vlastní inverzí:
Zjevná doména má nyní podmínku
ale toto může být prodlouženo. Podle výše uvedeného odkazu můžeme vzít integrál jako limit, protože horní limit jde do nekonečna (an nesprávný integrál spíše než a Lebesgueův integrál ), a tímto způsobem Hankelova transformace a její inverzní práce pro všechny funkce v L2 (0, ∞).
Transformace Laplaceovy rovnice
K transformaci a řešení lze použít Hankelovou transformaci Laplaceova rovnice vyjádřeno ve válcových souřadnicích. V rámci Hankelovy transformace se Besselův operátor stává násobením .[2] V osově symetrickém případě je parciální diferenciální rovnice transformována jako
což je obyčejná diferenciální rovnice v transformované proměnné .
je zvláštní případ Plancherelovy věty. Tyto věty lze dokázat pomocí vlastnosti ortogonality.
Vztah k vícerozměrné Fourierově transformaci
Hankelova transformace se objeví, když se zapíše vícerozměrná Fourierova transformace hypersférické souřadnice, což je důvod, proč se Hankelova transformace často objevuje ve fyzických problémech s válcovou nebo sférickou symetrií.
Zvažte funkci a -dimenzionální vektor r. Své -dimenzionální Fourierova transformace je definována jako
Abychom jej přepsali na hypersférické souřadnice, můžeme použít rozklad rovinné vlny na -dimenzionální hypersférické harmonické :[4]
kde a jsou množiny všech hypersférických úhlů v -prostor a -prostor. To dává následující výraz pro -dimenzionální Fourierova transformace v hypersférických souřadnicích:
Pokud se rozšíříme a v hypersférických harmonických:
Fourierova transformace v hypersférických souřadnicích se zjednodušuje
To znamená, že funkce s úhlovou závislostí ve formě hypersférické harmonické ji zachovávají na vícerozměrné Fourierově transformaci, zatímco radiální část prochází Hankelovou transformací (až na některé další faktory, jako je ).
Speciální případy
Fourierova transformace ve dvou dimenzích
Pokud je dvourozměrná funkce F(r) je rozšířen v a vícepólová řada,
pak je jeho dvourozměrná Fourierova transformace dána
kde
je -Hankelova transformace řádu (v tomto případě hraje roli momentu hybnosti, který byl označen v předchozí části).
pak je jeho trojrozměrná Fourierova transformace dána
kde
je Hankelova transformace řádu .
Tento druh Hankelovy transformace poločíselného řádu je také známý jako sférická Besselova transformace.
Fourierova transformace d rozměry (radiálně symetrický případ)
Pokud d-dimenzionální funkce F(r) nezávisí na úhlových souřadnicích, pak na d-dimenzionální Fourierova transformace F(k) také nezávisí na úhlových souřadnicích a je dán vztahem[5]
což je Hankelova transformace řádu až na faktor .
2D funkce uvnitř omezeného poloměru
Pokud je dvourozměrná funkce F(r) je rozšířen v a vícepólová řada a koeficienty roztažnosti Fm jsou dostatečně hladké v blízkosti počátku a nula vně poloměru R, radiální část F(r)/rm lze rozšířit na výkonovou řadu 1- (r / R) ^ 2:
tak, že dvojrozměrná Fourierova transformace F(r) se stává
kde poslední rovnost vyplývá z §6,567.1 ze dne.[6] Koeficienty roztažnosti Fm, t jsou přístupné pomocí diskrétní Fourierova transformace techniky:[7] pokud je radiální vzdálenost zmenšena
koeficienty Fourier-Čebyševovy řady G se objeví jako
Pomocí opětovné expanze
výnosy Fm, t vyjádřeno jako součet Gm, j.
To je jedna příchuť technik rychlé Hankelovy transformace.
Vztah k Fourierovým a Ábelovým transformacím
Hankelova transformace je jedním z členů FHA cyklus integrálních operátorů. Ve dvou dimenzích, pokud definujeme A jako Ábelova transformace operátor, F jako Fourierova transformace operátor a H jako operátor Hankelovy transformace s nulovým řádem, pak speciální případ věta o projekci a řezu pro kruhově symetrické funkce uvádí, že
Jinými slovy, použití Abelovy transformace na jednorozměrnou funkci a následné použití Fourierovy transformace na tento výsledek je stejné jako použití Hankelovy transformace na tuto funkci. Tuto koncepci lze rozšířit na vyšší dimenze.
Numerické hodnocení
Jednoduchý a efektivní přístup k numerickému vyhodnocení Hankelovy transformace je založen na pozorování, že může být sesláno ve formě konvoluce logaritmickou změnou proměnných[8]
V těchto nových proměnných čte Hankelova transformace
kde
Nyní lze integrál vypočítat numericky pomocí složitost použitím rychlá Fourierova transformace. Algoritmus lze dále zjednodušit použitím známého analytického výrazu pro Fourierovu transformaci :[9]
Optimální volba parametrů záleží na vlastnostech , zejména jeho asymptotické chování při a .
Tento algoritmus je známý jako „kvazi-rychlá Hankelova transformace“ nebo jednoduše „rychlá Hankelova transformace“.
^Avery, James Emil, autor. Hypersférické harmonické a jejich fyzikální aplikace. ISBN978-981-322-930-3. OCLC1013827621.CS1 maint: více jmen: seznam autorů (odkaz)
^Gradshteyn, I. S .; Ryzhik, I. M. (2015). Zwillinger, Daniel (ed.). Tabulka integrálů, sérií a produktů (Osmá ed.). Akademický tisk. p. 687. ISBN978-0-12-384933-5.
^Cree, M. J.; Bones, P.J. (červenec 1993). Msgstr "Algoritmy pro numerické vyhodnocení Hankelovy transformace". Počítače a matematika s aplikacemi. 26 (1): 1–12. doi:10.1016/0898-1221(93)90081-6. ISSN0898-1221.
^Papoulis, Athanasios (1981). Systémy a transformace s aplikacemi na optiku. Florida USA: Krieger Publishing Company. str. 140–175. ISBN978-0898743586.
Magni, Vittorio; Cerullo, Giulio; De Silverstri, Sandro (1992). "Vysoce přesná rychlá Hankelova transformace pro šíření optického paprsku". J. Opt. Soc. Dopoledne. A. 9 (11): 2031–2033. Bibcode:1992JOSAA ... 9.2031M. doi:10.1364 / JOSAA.9.002031.
Agnesi, A .; Reali, Giancarlo C .; Patrini, G .; Tomaselli, A. (1993). "Numerické vyhodnocení Hankelovy transformace: poznámky". Journal of the Optical Society of America A. 10 (9): 1872. Bibcode:1993JOSAA..10.1872A. doi:10.1364 / JOSAA.10.001872.
Barakat, Richard (1996). „Numerické vyhodnocení Hankelovy transformace nultého řádu pomocí Filonovy kvadraturní filozofie“. Aplikovaná matematická písmena. 9 (5): 21–26. doi:10.1016/0893-9659(96)00067-5. PAN1415467.
Zhang, D. W .; Yuan, X.-C .; Ngo, N. Q .; Shum, P. (2002). "Rychlá Hankelova transformace a její aplikace pro studium šíření válcových elektromagnetických polí". Opt. Vyjádřit. 10 (12): 521–525. Bibcode:2002Oexpr..10..521Z. doi:10,1364 / oe. 10.000521. PMID19436390.
Gizar-Sicairos, Manuel; Guitierrez-Vega, Julio C. (2004). "Výpočet kvazi-diskrétní Hankelovy transformace celočíselného řádu pro šíření polí optických vln". J. Opt. Soc. Dopoledne. A. 21 (1): 53–58. Bibcode:2004JOSAA..21 ... 53G. doi:10.1364 / JOSAA.21.000053. PMID14725397.