Biotechn3rnew, Technologia chemiczna pw, 2 rok, stata


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

(*) 0x01 graphic
,

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

0x01 graphic
)

ma osobliwo regularn w punkcie 0. Przykadem takiego równania jest np. tzw. równanie Bessela:

0x01 graphic
;

tutaj k jest tzw. wskanikiem, k"0. Tutaj f1(x)=1, f2(x)=-k2+0x+1x2.

Niech wspomnianymi szeregami potgowymi bd

0x01 graphic
.

Istnieje twierdzenie, e przy tych zaoeniach, rozwaane równanie ma przynaj-mniej jedno rozwizanie w postaci “uogólnionego szeregu potgowego”, tzw. szeregu Frobeniusa -

0x01 graphic
;

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

0x01 graphic

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

0x01 graphic

Po podstawieniu do równania otrzymujemy

0x01 graphic

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: 0x01 graphic
, czyli

(**) 0x01 graphic
,

bo poszukujemy rozwizania, w którym a0"0;

dla n=1: 0x01 graphic
, czyli

0x01 graphic
;

wreszcie dla n"2

0x01 graphic
.

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

0x01 graphic
,

gdzie z kolei u jest dane w postaci szeregu Frobeniusa 0x01 graphic
, gdzie r - mniejszy pierwiastek, wzgldnie - pierwiastek podwójny.

Dla równania Bessela

0x01 graphic
.

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 0x01 graphic
, tzn

0x01 graphic
.

Poniewa 2k jest tu >=0, wic

0x01 graphic
.

Poniewa a1=0, wic mamy a3=a5=...=a2n+1=0; natomiast dla parzystych n mamy

0x01 graphic

itd., ogólnie

0x01 graphic
,

czyli

0x01 graphic
.

Przyjmujc w szczególnoci 0x01 graphic
otrzymujemy funkcj, któr oznaczamy przez Jk i nazywamy funkcj Bessela pierwszego rodzaju rzdu k:

0x01 graphic
.

(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

0x01 graphic
.

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

0x01 graphic
.

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

0x01 graphic

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

0x01 graphic
,

gdzie 0x01 graphic
. Mona wykaza, e z dokadnoci do staej, drugie rozwizanie jest postaci

0x01 graphic

gdzie 0x01 graphic

Tutaj k=1,2,3,... . Rozwizanie ogólne jest wic wtedy postaci

0x01 graphic
.

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

0x01 graphic

dla x>0, gdzie z kolei

0x01 graphic

(zbieny dla wszystkich x).

Uwaga.

0x01 graphic
;

W praktyce najczciej wykorzystuje si J0 i J1, dlatego wypiszmy jeszcze

0x01 graphic

.

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.

0x01 graphic

D.:

0x01 graphic

Podobnie wykazuje si drug z równoci.

Z pierwszej z powyszych równoci otrzymujemy:

0x01 graphic
;

dzielc przez xk mamy

0x01 graphic

i podobnie

(*) 0x01 graphic
.

Std w miejscach zerowych

0x01 graphic
.

Mamy take wzór rekurencyjny

0x01 graphic

oraz wzór na róniczkowanie

0x01 graphic

Wemy pod uwag dwa równania Bessela

0x01 graphic

Cakami tych równa, ograniczonymi w otoczeniu punktu x=0 przy zaoeniu, e k>-1/2, s funkcje

0x01 graphic

Mnoc pierwsze z tych równa przez -z/x a drugie przez y/x i dodajc stronami, otrzymujemy

0x01 graphic
,

czyli

0x01 graphic

Cakujemy od 0 do 1; wobec ograniczonoci y i z w otoczeniu punktu 0, dostajemy

0x01 graphic
,

czyli

(**) 0x01 graphic

tzn.

0x01 graphic
.

Niech teraz alfa i beta bd dwoma rónymi pierwiastkami funkcji Jk(x), czyli

0x01 graphic
, gdzie n"m. Róniczkujc wzór (**) wzgldem  otrzymujemy

0x01 graphic
;

przyjmujc teraz =, dzielc przez 2 i przejmujc jako = miejsce zerowe funkcji Jk(x), otrzymujemy

0x01 graphic
,

gdzie ostatnia równo wynika z faktu, e

0x01 graphic
,

co z kolei wynika ze wzoru (*).



Wyszukiwarka

Podobne podstrony:
Biotechn3rnew, Technologia chemiczna pw, 2 rok, stata
Biotechn3rnew, Technologia chemiczna pw, 2 rok, stata
Biotechn3rnew, Technologia chemiczna pw, 2 rok, stata
Biotechn3rnew, Technologia chemiczna pw, 2 rok, stata
Biotechn3rnew, Technologia chemiczna pw, 2 rok, stata
Biotechn3rnew, Technologia chemiczna pw, 2 rok, stata
Biotechn3rnew, Technologia chemiczna pw, 2 rok, stata
Biotechn3rnew, Technologia chemiczna pw, 2 rok, stata
Biotechn3rnew, Technologia chemiczna pw, 2 rok, stata
30 Egzamin ECW 2006-01-30, Technologia chemiczna pw, 2 rok, stata
31 Egzamin ECW 2006-02-06, Technologia chemiczna pw, 2 rok, stata
Redoksometria, Technologia chemiczna pw, 2 rok, anality
Analiza straceniowa, Technologia chemiczna pw, 2 rok, anality
co gdzie jest, Technologia chemiczna pw, 2 rok, infa
kol2, Technologia chemiczna pw, 2 rok, infa
kolos1, Technologia chemiczna pw, 2 rok, infa
Opracowanko zestawuf, Technologia chemiczna pw, 2 rok, anality

więcej podobnych podstron