Równania róniczkowe drugiego rzdu z osobliwoci regularn w punkcie 0.
Funkcje Bessela (walcowe) - definicja i wasnoci.
Rozwamy równanie róniczkowe drugiego rzdu postaci
(*)
,
gdzie f1 i f2 rozkadaj si w szeregi potgowe w pewnym otoczeniu punktu x0=0. Mówimy wtedy, e to równanie (lub, równowane mu równanie
)
ma osobliwo regularn w punkcie 0. Przykadem takiego równania jest np. tzw. równanie Bessela:
;
tutaj k jest tzw. wskanikiem, k"0. Tutaj f1(x)=1, f2(x)=-k2+0x+1x2.
Niech wspomnianymi szeregami potgowymi bd
.
Istnieje twierdzenie, e przy tych zaoeniach, rozwaane równanie ma przynaj-mniej jedno rozwizanie w postaci “uogólnionego szeregu potgowego”, tzw. szeregu Frobeniusa -
;
szereg ten jest zbieny przynajmniej w obszarze wspólnej zbienoci szeregów na f1 i f2, z wyjtkiem co najwyej punktu x0 = 0.
Poszukajmy rozwizania równania Bessela w postaci szeregu Frobeniusa. Niech
W ogólnym przypadku podstawienie do równania (*) wymaga mnoenia szeregów, a wic grupowania wyrazów danego stopnia w iloczynie szeregów. W przypadku równania Bessela, wszystko si znacznie upraszcza, poniewa zarówno f1, jak i f2 s wielomianami. Mamy wic
Po podstawieniu do równania otrzymujemy
Dzielc przez xr (x"0) i korzystajc z faktu, e na to, aby szereg potgowy zerowa si dla x z pewnego ssiedztwa punktu x0=0 potrzeba (i wystarcza), aby wszystkie jego wspóczynniki byy równe zeru, otrzymujemy:
dla n=0:
, czyli
(**)
,
bo poszukujemy rozwizania, w którym a0"0;
dla n=1:
, czyli
;
wreszcie dla n"2
.
Rozwamy przede wszystkim pierwsze z tych równa, tzn. (**). Odgrywa ono rol równania charakterystycznego. Jeeli równanie to ma pierwiastki rzeczywiste, to ilo liniowo niezalenych rozwiza w postaci szeregu Frobeniusa zaley od pierwiastków tego równania, a dokadniej, od tego, czy s
a) dwa pierwiastki nie rónice si o liczb cakowit (wtedy otrzymujemy dwa liniowo niezalene rozwizania w postaci szeregów Frobeniusa);
b) dwa pierwiastki rónice si o liczb cakowit, powiedzmy r (mniejszy) i r+d (wikszy); wtedy,
ba) albo przy wziciu mniejszego pierwiastka równanie na ad jest sprzeczne, i otrzymujemy tylko rozwizanie otrzymane dla wikszego pierwiastka,
bb) albo te, przy wziciu mniejszego pierwiastka, równanie na ad ma posta 0=0, wtedy otrzymujemy dwa rozwizania - zalene od staych dowolnych a0 i ad (rozwizanie otrzymane dla r+d jest wtedy liniowo zalene od nich i nie bierzemy go ju wicej pod uwag);
c) jeden pierwiastek podwójny - wtedy otrzymujemy jedno rozwizanie w postaci szeregu Frobeniusa.
We wszystkich przypadkach, kiedy powysze procedury daj nam tylko jedno rozwizanie (powiedzmy y1), drugiego rozwizania szuka si w postaci
,
gdzie z kolei u jest dane w postaci szeregu Frobeniusa
, gdzie r - mniejszy pierwiastek, wzgldnie - pierwiastek podwójny.
Dla równania Bessela
.
Rozwamy najpierw wikszy pierwiastek r=k; drugie równanie daje
(2k+1)a1=0, czyli a1=0 (2k+1>0, jako e k jest >=0).
Nastpnie
, tzn
.
Poniewa 2k jest tu >=0, wic
.
Poniewa a1=0, wic mamy a3=a5=...=a2n+1=0; natomiast dla parzystych n mamy
itd., ogólnie
,
czyli
.
Przyjmujc w szczególnoci
otrzymujemy funkcj, któr oznaczamy przez Jk i nazywamy funkcj Bessela pierwszego rodzaju rzdu k:
.
(Szereg ten jest zbieny dla wszystkich x, co wynika z kryterium d'Alemberta.)
Analogicznie dla k " 0, 1, 2, ... otrzymujemy drugie rozwizanie i jest to dokadnie
.
Tak wic dla k"0,1,2 otrzymujemy dwa liniowo niezalene rozwizania w postaci szeregu Frobeniusa Jk i J-k, a rozwizaniem ogólnym równania jest
.
Chocia dla k=0,1,2,... szereg na J-k(x) formalnie traci sens, mona mu nada sens przyjmujc, e wyrazy, dla których funkcja gamma w mianowniku dy do nieskoczonoci s równe zeru, i wtedy otrzymujemy
Tak wic dla k=0,1,2,... funkcje Jk i J-k s liniowo zalene i rozwizanie ogólne równania Bessela nie moe mie takiej postaci, jak dla k niecakowitych.
Dla k=1,2,... drugie rozwizanie jest postaci
,
gdzie
. Mona wykaza, e z dokadnoci do staej, drugie rozwizanie jest postaci
gdzie
Tutaj k=1,2,3,... . Rozwizanie ogólne jest wic wtedy postaci
.
Funkcj Nk nazywamy funkcj Neumana rzdu k lub te funkcj Bessela drugiego rodzaju rzdu k. Powysze rozwizanie jest wane take dla k=0, jeeli zdefiniowa
dla x>0, gdzie z kolei
(zbieny dla wszystkich x).
Uwaga.
;
W praktyce najczciej wykorzystuje si J0 i J1, dlatego wypiszmy jeszcze
.
Wasnoci funkcji Bessela
A. Miejsca zerowe funkcji Bessela.
1. Jeeli x1 i x2 s miejscami zerowymi funkcji Jk(x), to w przedziale (x1,x2) istnieje miejsce zerowe funkcji Jk-1(x) i Jk+1(x).
2. Funkcja Bessela J0(x) ma miejsce zerowe w kadym przedziale o dugoci .
3. Kada funkcja Bessela Jk(x), k=1,2, ... ma nieskoczenie wiele miejsc zerowych w przedziale (0,").
4. Jeeli k>1/2, to rónica midzy dwoma kolejnymi miejscami zerowymi Jk(x) jest wiksza ni .
5. Jeeli k>1/2, to rónica midzy dwoma kolejnymi miejscami zerowymi Jk(x) dy do , gdy x dy do nieskoczonoci.
6. Pierwsze dodatnie miejsce zerowe Jk(x) jest wiksze ni k.
7. Funkcja Bessela Jk(x) ma tylko rzeczywiste miejsca zerowe.
B. Róne zwizki pomidzy Jk, Jk-1 i Jk+1.
D.:
Podobnie wykazuje si drug z równoci.
Z pierwszej z powyszych równoci otrzymujemy:
;
dzielc przez xk mamy
i podobnie
(*)
.
Std w miejscach zerowych
.
Mamy take wzór rekurencyjny
oraz wzór na róniczkowanie
Wemy pod uwag dwa równania Bessela
Cakami tych równa, ograniczonymi w otoczeniu punktu x=0 przy zaoeniu, e k>-1/2, s funkcje
Mnoc pierwsze z tych równa przez -z/x a drugie przez y/x i dodajc stronami, otrzymujemy
,
czyli
Cakujemy od 0 do 1; wobec ograniczonoci y i z w otoczeniu punktu 0, dostajemy
,
czyli
(**)
tzn.
.
Niech teraz alfa i beta bd dwoma rónymi pierwiastkami funkcji Jk(x), czyli
, gdzie n"m. Róniczkujc wzór (**) wzgldem otrzymujemy
;
przyjmujc teraz =, dzielc przez 2 i przejmujc jako = miejsce zerowe funkcji Jk(x), otrzymujemy
,
gdzie ostatnia równo wynika z faktu, e
,
co z kolei wynika ze wzoru (*).