Metody matematyczne fizyki/Funkcje sferyczne w matematyce

Metody matematyczne fizyki
Metody matematyczne fizyki
Funkcje sferyczne w matematyce

Licencja
Autor: Mirosław Makowiecki
Absolwent UMCS Fizyki Komputerowej Uniwersytetu Marii Curie-Skłodowskiej w Lublinie
Email: miroslaw(kropka)makowiecki(małpa)gmail(kropka)pl
Dotyczy: książki, do której należy ta strona, oraz w niej zawartych stron i w nich podstron, a także w nich kolumn, wraz z zawartościami.
Użytkownika książki, do której należy ta strona, oraz w niej zawartych stron i w nich podstron, a także w nich kolumn, wraz z zawartościami nie zwalnia z odpowiedzialności prawnoautorskiej nieprzeczytanie warunków licencjonowania.
Umowa prawna: Creative Commons: uznanie autorstwa, na tych samych warunkach, z możliwością obowiązywania dodatkowych ograniczeń.
Autor tej książki dołożył wszelką staranność, aby informacje zawarte w książce były poprawne i najwyższej jakości, jednakże nie udzielana jest żadna gwarancja, czy też rękojma. Autor nie jest odpowiedzialny za wykorzystanie informacji zawarte w książce, nawet jeśli wywołaby jakąś szkodę, straty w zyskach, zastoju w prowadzeniu firmy, przedsiębiorstwa lub spółki bądź utraty informacji, niezależnie czy autor (a nawet Wikibooks) został powiadomiony o możliwości wystąpienie szkód. Informacje zawarte w książce mogą być wykorzystane tylko na własną odpowiedzialność.


Następny rozdział: Funkcje Bessela. Poprzedni rozdział: Wprowadzenie do wielomianów ortogonalnych.

Podręcznik: Metody matematyczne fizyki.

W tym rozdziale zajmować się będziemy funkcjami sferycznymi, która ta funkcja zależy od dwóch współrzędnych kątowych, którego to zapis Y=Ylm. Funkcje Ylm są wzajemnie ortogonalne do siebie, tzn. to zachodzi związek mówiącej o tej ortogonalności.

(10.1)

Iloczynem skalarnym dwóch funkcji zależnych od funkcji kątowych określamy wedle warunku:

(10.2)

Równanie różniczkowe, którego rozwiązaniem są funkcje sferyczne

edytuj

Równaniem, które się tutaj zajmować się tutaj będziemy, jest to równanie zależne od zmiennych kulistego układu współrzędnych, czyli zależy od zmiennych (r,φ.ξ), określmy przez równanie:

(10.3)

Funkcję f piszemy jako iloczyn funkcji potęgowej rl i funkcji kulistej Y(θ,φ):

(10.4)

Tożsamość (10.4) podstawiamy do równości różniczkowej (10.3), wtedy dostajemy następnie równanie różniczkowe, to dowód tego przejścia pokazujemy:



(10.5)

Takie samo równanie, co równanie końcowe wynikające przeprowadzonej w punkcie (10.5) otrzymamy z (10.3), tzn. gdy weżniemy dla ujemnych potęg r-l-1, czyli dla rozwiązania:

(10.6)

Co pokazują obliczenia, które wynikają z dodawania czynników wynikające z dwóch pierwszych pochodnych (pierwsza i druga pochodna funkcji) przy czynniku rl-3:

Funkcja Laplace'a jest funkcją, która jest kombinacją dwóch wyrazów z potęgą ujemną oraz dodatnią, czyli jest kombinacją liniową funkcji (10.4) i funkcji (10.6). Funkcja Y(θ,φ) jest iloczynem dwóch wyrazów, tzn. T(ξ) i Φ(θ), co opis tego rozwiązania zapisujemy:

(10.7)

Możemy podstawić rozwiązanie równania kulistego (10.7) do równania różniczkowego (10.5), w takim przypadku możemy napisać od razu po dokonaniu tejże opisywanej operacji:

(10.8)

Równanie różniczkowe (10.8) dzielimy obustronnie przez wyrażenie zapisane wzorem (10.7), następnie po krótkich przekształceniach i po krótkich przenosinach pewnych wyrazów w taki sposób by po na prawej i lewej stronie rozważanego równania znajdowały się różne funkcje zależne od różnych zmiennych, wtedy dostajemy tożsamość:

(10.9)

Funkcje znajdujące się po obu stronach równania (10.9) są zależne od innych zmiennych, obie te funkcje są tożsamościowo równe sobie, zatem obie jego strony możemy przyrównać do stałej. Aby tak się stało musimy napisać:

(10.10)

Aby końcowa tożsamość wynikowa (10.10) była poprawnie sformułowana, to ona na pewno powinna mieć takie same wartości dla φ=0 i 2π:

(10.11)

Wedle obliczeń przeprowadzonych w punkcie (10.11) funkcja (10.10) przyjmuje ostateczną końcową postać:

(10.12)

Jeśli równanie na Θ(θ) (10.12) podstawimy do równania (10.9), otrzymujemy następne równanie różniczkowe na funkcję T, w takim przypadku:


(10.13)

W równaniu (10.13), która określa pewną funkcję T w zależności od parametru m i zmiennej ξ. dokonujemy podstawienia, która jest definicją wyrazu (1-ξ2)m/2 i pewnej funkcji Plm.

(10.14)

Następnym naszym krokiem jest policzenie składnika pierwszego, która znajduje się w równaniu końcowym wynikowym (10.13), zatem policzmy ten wyraz korzystając z definicji funkcji T poprzez iloczyn dwóch funkcji zależnych od ξ (10.14):




(10.15)

Za pierwszy składnik sumy występującej w sumie po lewej stronie równania różniczkowego (10.13) podstawiamy wyrażenie wynikłe z obliczeń (10.15), to:


(10.16)

W celu dalszego przekształcenia równania (10.16) należy dokonać obliczeń dla wyrazów stojących przy czynniku Plm, ale będziemy liczyć bez tego czynnika, co jest dla nas oczywiste dlaczego tak robimy.



(10.17)

Do równania (10.16) możemy podstawić wynik obliczeń przeprowadzonych w punkcie (10.17) dla wyrazów stojących przy Plm, wtedy mamy wniosek:

(10.18)

Równanie różniczkowe (10.18) dzielimy obustronnie przez wyrażenie , co można tak zrobić przy założeniu, że mamy , w takim przypadku równanie (10.18) przyjmuje postać:

(10.19)

W równaniu (10.19), gdy m zastąpimy przez zero, wtedy otrzymamy równanie Legendre'a wprowadzone w punkcie (9.16), zatem możemy powiedzieć, że Pl0 są to wielomiany Legendre'a. Równanie różniczkowe (10.19) możemy zróżniczkować względem zmiennej ξ, w takim przypadku mamy równość wynikową w porównaniu do poprzedniej równości:



(10.20)

Ze wzoru końcowego możemy powiedzieć, że spełniona jest tożsamość:

(10.21)

Wykorzystując wzór na wielomiany Legendre'a, które są zdefiniowane wzorem (9.33) oraz wzór na funkcje sferyczne (10.21) dochodzimy wtedy do wniosku, że wielomiany, które są rozwiązaniami równania (10.19), które to ogólnie przedstawienie przez Plm możemy napisać:

(10.22)

Z otrzymanego wzoru (10.22) dochodzimy do wniosku, że parametr m występujący we wzorze na równanie różniczkowe (10.19) przebiega wartości:

(10.23)

Końcowy wzór na rozwiązanie równania (10.13) jest w postaci:

(10.24)

Funkcje sferyczne, a jego ortogonalność

edytuj

Mając funkcje sferyczne (10.24), można udowodnić, że te funkcje wraz z odpowiednio dobraną stałą N jest funkcją ortogonalną, zatem z definicji iloczynu skalarnego i wspomnianej funkcji Ylm możemy napisać tożsamość:

(10.25)

Ostatnia całka oznaczona występująca we wzorze na iloczyn skalarny dwóch funkcji kulistych jest równa zero, gdy , a jest różna od zera, gdy m=m' i jest ona równa 2π, w takim razie (10.25) możemy przepisać w postaci:

(10.26)

Wyrażenie (10.25) całkujemy przez części i wyrazy stojące pod całką są równe zero dla ξ=-1,1, zatem takim bądź razem dostajemy wniosek:

(10.27)

We wzorze (10.27), gdy zachodzi l'>l, toteż całka oznaczona występująca w tym samym wzorze jest równa zero, co wynika że każdy wielomian ortogonalny ma iloczyn skalarny równy zero, dla stopnia niższego niż posiadany stopień przez wielomian Pl, w takim przypadku całka występująca we wzorze (10.27) jest różna od zera, zatem w takim razie możemy napisać tożsamość, czemu jest równe (10.25):

(10.28)

Normowanie funkcji sferycznych

edytuj

Gdy we wzorze (10.26) podstawimy l=l' i m=m', wtedy ten wspomniany wzór przedstawia się:

(10.29)

Pod całką wystepuje wielomian , który jest rzędu l, ale również wielomian Pl, też jest takiego rzędu, w takim razie całka występująca we wzorze (10.27) nie jest równa zero. Rozważmy wielomian W(ξ) występujący we wzorze (10.29), zatem jest taki, że potęgi mniejsze niż l, są wycinane przez wielomian Legendre'a , bo to jest wielomian ortogonalny, w takim przypadku w całce należy posługiwać się tylko największą potęgą tegoż wielomianu, bo pozostałe potęgi są wycinane.

Możemy teraz określić największy wyraz stojący przy wielomianie W(ξ), w takim przypadku możemy napisać czynnik stojący przy największej potędze w wspominanym wielomianie, zatem ten współczynnik:

(10.30)

Liczby całkowite występujące w pierwszym wyrazie pochodzą z różniczkowania wielomianu Pl, a liczby całkowite w drugim czynniku pochodzą z różniczkowania całego wyrażenia w dużym nawiasie. Wzór (10.30) możemy napisać w postaci zwartej:

(10.31)

Zatem norma funkcji kulistej Ylm (10.29), wiedząc że w iloczynie skalarnym iloczyn alξl można uzupełnić do pełnego wielomianu Pl i zarazem znając normę wielomianu Legendre'a (9.57), określamy według:

(10.32)

Ze wzoru wynikającego z (10.32) możemy napisać stałą Nlm, bo norma funkcji kulistej jest równa jeden, jako:

(10.33)

Pełny wzór (10.24) uzyskujemy podstawiając do niego tożsamość za stałą Nlm (10.33):

(10.34)